Functional dissection of the Sox9-Kcnj2 locus identifies nonessential and instructive roles of TAD architecture

The genome is organized in three-dimensional units called topologically associating domains (TADs), through a process dependent on the cooperative action of cohesin and the DNA-binding factor CTCF. Genomic rearrangements of TADs have been shown to cause gene misexpression and disease, but genome-wide depletion of CTCF has no drastic effects on transcription. Here, we investigate TAD function in vivo in mouse limb buds at the Sox9-Kcnj2 locus. We show that the removal of all major CTCF sites at the boundary and within the TAD resulted in a fusion of neighboring TADs, without major effects on gene expression. Gene misexpression and disease phenotypes, however, were achieved by redirecting regulatory activity through inversions and/or the repositioning of boundaries. Thus, TAD structures provide robustness and precision but are not essential for developmental gene regulation. Aberrant disease-related gene activation is not induced by a mere loss of insulation but requires CTCF-dependent redirection of enhancer-promoter contacts

Promoter repression and 3D-restructuring resolves divergent developmental gene expression in TADs

Cohesin loop extrusion facilitates precise gene expression by continuously driving promoters to sample all enhancers located within the same topologically-associated domain (TAD). However, many TADs contain multiple genes with divergent expression patterns, thereby indicating additional forces further refine how enhancer activities are utilised. Here, we unravel the mechanisms enabling a new gene, Rex1, to emerge with divergent expression within the ancient Fat1 TAD in placental mammals. We show that such divergent expression is not determined by a strict enhancer-promoter compatibility code, intra-TAD position or nuclear envelope-attachment. Instead, TAD-restructuring in embryonic stem cells (ESCs) separates Rex1 and Fat1 with distinct proximal enhancers that independently drive their expression. By contrast, in later embryonic tissues, DNA methylation renders the inactive Rex1 promoter profoundly unresponsive to Fat1 enhancers within the intact TAD. Combined, these features adapted an ancient regulatory landscape during evolution to support two entirely independent Rex1 and Fat1 expression programs. Thus, rather than operating only as rigid blocks of co-regulated genes, TAD-regulatory landscapes can orchestrate complex divergent expression patterns in evolution

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

The little skate genome and the evolutionary emergence of wing-like fin appendages

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

Characterization of hundreds of regulatory landscapes in developing limbs reveals two regimes of chromatin folding

Complex regulatory landscapes control the pleiotropic transcriptional activities of developmental genes. For most genes, the number, location and dynamic of their associated regulatory elements are unknown. In this work we characterized the three-dimensional chromatin micro-architecture and regulatory landscape of 446 limb-associated gene loci in mouse using Capture-C, ChIP-seq and RNA-seq in forelimb, hindlimb at three developmental stages and midbrain. The fine mapping of chromatin interactions revealed a strong preference for functional genomic regions such as repressed or active domains. By combining chromatin marks and interaction peaks, we annotated more than one thousand putative limb enhancers and their associated genes. Moreover, the analysis of chromatin interactions revealed two regimes of chromatin folding, one producing interactions stable across tissues and stages and another one associated with tissue and/or stage-specific interactions. Whereas stable interactions associate strongly with CTCF/RAD21 binding, the intensity of variable interactions correlates with changes in underlying chromatin modifications, specifically at the viewpoint and at the interaction site. In conclusion, this comprehensive dataset provides a resource for the characterization of hundreds of limb-associated regulatory landscapes and a framework to interpret the chromatin folding dynamics observed during embryogenesis

Deletions, Inversions, Duplications: Engineering of Structural Variants using CRISPR/Cas in Mice

Structural variations (SVs) contribute to the variability of our genome and are often associated with disease. Their study in model systems was hampered until now by labor-intensive genetic targeting procedures and multiple mouse crossing steps. Here we present the use of CRISPR/Cas for the fast (10 weeks) and efficient generation of SVs in mice. We specifically produced deletions, inversions, and also duplications at six different genomic loci ranging from 1.1 kb to 1.6 Mb with efficiencies up to 42%. After PCR-based selection, clones were successfully used to create mice via aggregation. To test the practicability of the method, we reproduced a human 500 kb disease-associated deletion and were able to recapitulate the human phenotype in mice. Furthermore, we evaluated the regulatory potential of a large genomic interval by deleting a 1.5 Mb fragment. The method presented permits rapid in vivo modeling of genomic rearrangements