Search

SUMMARYCohesin 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.HIGHLIGHTSNew genes can emerge in evolution without taking on the expression pattern of their surrounding pre-existing TAD.Compartmentalisation can restructure seemingly evolutionarily stable TADs to control a promoter’s access to enhancers.Lamina-associated domains neither prevent transcriptional activation nor enhancer-promoter communication.Repression rather than promoter-specificity refines when genes respond to promiscuous enhancer activities in specific tissues.

The genome is organized in megabase-sized three-dimensional units, called Topologically Associated Domains (TADs), that are separated by boundaries. TADs bring distant cis-regulatory elements into proximity, a process dependent on the cooperative action of cohesin and the DNA binding factor CTCF. Surprisingly, genome-wide depletion of CTCF has little effect on transcription, yet structural variations affecting TADs have been shown to cause gene misexpression and congenital disease. Here, we investigate TAD functionin vivoin mice by systematically editing components of TAD organization at theSox9/Kcnjlocus. We find that TADs are formed by a redundant system of CTCF sites requiring the removal of all major sites within the TAD and at the boundary for two neighboring TADs to fuse. TAD fusion resulted in leakage of regulatory activity from theSox9to theKcnjTAD, but no major changes in gene expression. This indicates that TAD structures provide robustness and precision, but are not essential for developmental gene regulation. Gene misexpression and resulting disease phenotypes, however, were attained by re-directing regulatory activity through inversions and/or the re-positioning of boundaries. Thus, efficient re-wiring of enhancer promoter interaction and aberrant disease causing gene activation is not induced by a mere loss of insulation but requires the re-direction of contacts.

Rare diseases and their underlying molecular causes are often poorly studied, posing challenges for patient diagnosis and prognosis. The development of next-generation sequencing and its decreasing costs promises to alleviate such issues by supplying personal genomic information at a moderate price. Here, we used crowdfunding as an alternative funding source to sequence the genome of Lil BUB, a celebrity cat affected by rare disease phenotypes characterized by supernumerary digits, osteopetrosis and dwarfism, all phenotypic traits that also occur in human patients. We discovered that Lil BUB is affected by two distinct mutations: a heterozygous mutation in the limb enhancer of the Sonic hedgehog gene, previously associated with polydactyly in Hemingway cats; and a novel homozygous frameshift deletion affecting the TNFRSF11A (RANK) gene, which has been linked to osteopetrosis in humans. We communicated the progress of this project to a large online audience, detailing the ‘inner workings’ of personalized whole genome sequencing with the aim of improving genetic literacy. Our results highlight the importance of genomic analysis in the identification of disease-causing mutations and support crowdfunding as a means to fund low-budget projects and as a platform for scientific communication.

Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq) is a powerful technology to identify the genome-wide locations of transcription factors and other DNA binding proteins. Computational ChIP-seq peak calling infers the location of protein–DNA interactions based on various measures of enrichment of sequence reads. In this work, we introduce an algorithm, Q, that uses an assessment of the quadratic enrichment of reads to center candidate peaks followed by statistical analysis of saturation of candidate peaks by 5′ ends of reads. We show that our method not only is substantially faster than several competing methods but also demonstrates statistically significant advantages with respect to reproducibility of results and in its ability to identify peaks with reproducible binding site motifs. We show that Q has superior performance in the delineation of double RNAPII and H3K4me3 peaks surrounding transcription start sites related to a better ability to resolve individual peaks. The method is implemented in C+l+ and is freely available under an open source license.

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.