Marion Leleu

The Dynamic Architecture of Hox Gene Clusters

Sequential activation of Hox genes correlates with a transition of negative to positive three-dimensional chromosome structure. The spatial and temporal control of Hox gene transcription is essential for patterning the vertebrate body axis. Although this process involves changes in histone posttranslational modifications, the existence of particular three-dimensional (3D) architectures remained to be assessed in vivo. Using high-resolution chromatin conformation capture methodology, we examined the spatial configuration of Hox clusters in embryonic mouse tissues where different Hox genes are active. When the cluster is transcriptionally inactive, Hox genes associate into a single 3D structure delimited from flanking regions. Once transcription starts, Hox clusters switch to a bimodal 3D organization where newly activated genes progressively cluster into a transcriptionally active compartment. This transition in spatial configurations coincides with the dynamics of chromatin marks, which label the progression of the gene clusters from a negative to a positive transcription status. This spatial compartmentalization may be key to process the colinear activation of these compact gene clusters.

Conservation and Divergence of Regulatory Strategies at Hox Loci and the Origin of Tetrapod Digits

During development, expression of the Hoxa and Hoxd genes in zebrafish fins and mouse limbs are regulated via a conserved chromatin structure. However, zebrafish lack certain regulatory elements required to produce digits, revealing that radials—the fins bony elements—are likely not homologous to tetrapod digits.

Functional Analysis of CTCF During Mammalian Limb Development

CCCTC-binding factor (CTCF) is a nuclear zinc-finger protein that displays insulating activity in a variety of biological assays. For example, CTCF-binding sites have been suggested to isolate Hox gene clusters from neighboring transcriptional interference. We investigated this issue during limb development, where Hoxd genes must remain isolated from long-range effects to allow essential regulation within independent sub-groups. We used conditional Ctcf inactivation in incipient forelimbs and show that the overall pattern of Hoxd gene expression remains unchanged. Transcriptome analysis using tiling arrays covering chromosomes 2 and X confirmed the weak effect of CTCF depletion on global gene regulation. However, Ctcf deletion caused massive apoptosis, leading to a nearly complete loss of limb structure at a later stage. We conclude that, at least in this physiological context, rather than being an insulator, CTCF is required for cell survival via the direct transcriptional regulation of target genes critical for cellular homeostasis.

Temporal dynamics and developmental memory of 3D chromatin architecture at Hox gene loci

Hox genes are essential regulators of embryonic development. Their step-wise transcriptional activation follows their genomic topology and the various states of activation are subsequently memorized into domains of progressively overlapping gene products. We have analyzed the 3D chromatin organization of Hox clusters during their early activation in vivo, using high-resolution circular chromosome conformation capture. Initially, Hox clusters are organized as single chromatin compartments containing all genes and bivalent chromatin marks. Transcriptional activation is associated with a dynamic bi-modal 3D organization, whereby the genes switch autonomously from an inactive to an active compartment. These local 3D dynamics occur within a framework of constitutive interactions within the surrounding Topological Associated Domains, indicating that this regulation process is mostly cluster intrinsic. The step-wise progression in time is fixed at various body levels and thus can account for the chromatin architectures previously described at a later stage for different anterior to posterior levels. DOI: http://dx.doi.org/10.7554/eLife.02557.001

Clustering of mammalian Hox genes with other H3K27me3 targets within an active nuclear domain

Significance The development of an embryo from a single fertilized cell is orchestrated by a large set of key regulatory genes whose activities need to be precisely controlled. Proteins from the Polycomb group (PcG) family maintain the inactive state of these genes by modifying the surrounding histone H3 tails. Here we report that these inactive genes contact other PcG-rich regions within otherwise active environments in the cell nucleus, suggesting the presence of repressive microenvironments. PcG positive genes interact independently from neighboring genes, depending on their linear distance and more global chromosome folding characteristics, with only moderate changes during embryonic development. Embryogenesis requires the precise activation and repression of many transcriptional regulators. The Polycomb group proteins and the associated H3K27me3 histone mark are essential to maintain the inactive state of many of these genes. Mammalian Hox genes are targets of Polycomb proteins and form local 3D clusters centered on the H3K27me3 mark. More distal contacts have also been described, yet their selectivity, dynamics, and relation to other layers of chromatin organization remained elusive. We report that repressed Hox genes form mutual intra- and interchromosomal interactions with other genes located in strong domains labeled by H3K27me3. These interactions occur in a central and active nuclear environment that consists of the HiC compartment A, away from peripheral lamina-associated domains. Interactions are independent of nearby H3K27me3-marked loci and determined by chromosomal distance and cell-type–specific scaling factors, thus inducing a moderate reorganization during embryogenesis. These results provide a simplified view of nuclear organization whereby Polycomb proteins may have evolved to repress genes located in gene-dense regions whose position is restricted to central, active, nuclear environments.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

HTSstation: A Web Application and Open-Access Libraries for High-Throughput Sequencing Data Analysis

The HTSstation analysis portal is a suite of simple web forms coupled to modular analysis pipelines for various applications of High-Throughput Sequencing including ChIP-seq, RNA-seq, 4C-seq and re-sequencing. HTSstation offers biologists the possibility to rapidly investigate their HTS data using an intuitive web application with heuristically pre-defined parameters. A number of open-source software components have been implemented and can be used to build, configure and run HTS analysis pipelines reactively. Besides, our programming framework empowers developers with the possibility to design their own workflows and integrate additional third-party software. The HTSstation web application is accessible at http://htsstation.epfl.ch.

Structural Variation-Associated Expression Changes Are Paralleled by Chromatin Architecture Modifications

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this “neighboring effect”, we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples. We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype. GEO Series accession number: GSE33784, GSE33867.

Processing and analyzing ChIP-seq data: from short reads to regulatory interactions

Chromatin-immunoprecipitation and sequencing (ChIP-seq) is a rapidly maturing technology that draws on the power of high-throughput short-read sequencing to decipher chromatin states with unprecedented precision and breadth. Although some aspects of the experimental protocol require careful tuning, the bottleneck currently firmly lies with the downstream data analysis. We give an overview of the better-established aspects of genome mapping and data normalization and we describe the more recent progress in peak calling and their statistical analysis and provide a brief overview of popular follow-up analyses such as genomic feature categorization and motif search.

Large scale genomic reorganization of topological domains at the HoxD locus

BackgroundThe transcriptional activation of HoxD genes during mammalian limb development involves dynamic interactions with two topologically associating domains (TADs) flanking the HoxD cluster. In particular, the activation of the most posterior HoxD genes in developing digits is controlled by regulatory elements located in the centromeric TAD (C-DOM) through long-range contacts.ResultsTo assess the structure–function relationships underlying such interactions, we measured compaction levels and TAD discreteness using a combination of chromosome conformation capture (4C-seq) and DNA FISH. We assessed the robustness of the TAD architecture by using a series of genomic deletions and inversions that impact the integrity of this chromatin domain and that remodel long-range contacts. We report multi-partite associations between HoxD genes and up to three enhancers. We find that the loss of native chromatin topology leads to the remodeling of TAD structure following distinct parameters.ConclusionsOur results reveal that the recomposition of TAD architectures after large genomic re-arrangements is dependent on a boundary-selection mechanism in which CTCF mediates the gating of long-range contacts in combination with genomic distance and sequence specificity. Accordingly, the building of a recomposed TAD at this locus depends on distinct functional and constitutive parameters.