Delineation of metabolic gene clusters in plant genomes by chromatin signatures

Plants are a tremendous source of diverse chemicals, including many natural product-derived drugs. It has recently become apparent that the genes for the biosynthesis of numerous different types of plant natural products are organized as metabolic gene clusters, thereby unveiling a highly unusual form of plant genome architecture and offering novel avenues for discovery and exploitation of plant specialized metabolism. Here we show that these clustered pathways are characterized by distinct chromatin signatures of histone 3 lysine trimethylation (H3K27me3) and histone 2 variant H2A.Z, associated with cluster repression and activation, respectively, and represent discrete windows of co-regulation in the genome. We further demonstrate that knowledge of these chromatin signatures along with chromatin mutants can be used to mine genomes for cluster discovery. The roles of H3K27me3 and H2A.Z in repression and activation of single genes in plants are well known. However, our discovery of highly localized operon-like co-regulated regions of chromatin modification is unprecedented in plants. Our findings raise intriguing parallels with groups of physically linked multi-gene complexes in animals and with clustered pathways for specialized metabolism in filamentous fungi.

4C-ker: A Method to Reproducibly Identify Genome-Wide Interactions Captured by 4C-Seq Experiments

4C-Seq has proven to be a powerful technique to identify genome-wide interactions with a single locus of interest (or "bait") that can be important for gene regulation. However, analysis of 4C-Seq data is complicated by the many biases inherent to the technique. An important consideration when dealing with 4C-Seq data is the differences in resolution of signal across the genome that result from differences in 3D distance separation from the bait. This leads to the highest signal in the region immediately surrounding the bait and increasingly lower signals in far-cis and trans. Another important aspect of 4C-Seq experiments is the resolution, which is greatly influenced by the choice of restriction enzyme and the frequency at which it can cut the genome. Thus, it is important that a 4C-Seq analysis method is flexible enough to analyze data generated using different enzymes and to identify interactions across the entire genome. Current methods for 4C-Seq analysis only identify interactions in regions near the bait or in regions located in far-cis and trans, but no method comprehensively analyzes 4C signals of different length scales. In addition, some methods also fail in experiments where chromatin fragments are generated using frequent cutter restriction enzymes. Here, we describe 4C-ker, a Hidden-Markov Model based pipeline that identifies regions throughout the genome that interact with the 4C bait locus. In addition, we incorporate methods for the identification of differential interactions in multiple 4C-seq datasets collected from different genotypes or experimental conditions. Adaptive window sizes are used to correct for differences in signal coverage in near-bait regions, far-cis and trans chromosomes. Using several datasets, we demonstrate that 4C-ker outperforms all existing 4C-Seq pipelines in its ability to reproducibly identify interaction domains at all genomic ranges with different resolution enzymes.

How to rule the nucleus: divide et impera

Genome-wide molecular studies have provided new insights into the organization of nuclear chromatin by revealing the presence of chromatin domains of differing transcriptional activity, frequency of cis-interactions, proximity to scaffolding structures and replication timing. These studies have not only brought our understanding of genome function to a new level, but also offered functional insight for many phenomena observed in microscopic studies. In this review, we discuss the major principles of nuclear organization based on the spatial segregation of euchromatin and heterochromatin, as well as the dynamic genome rearrangements occurring during cell differentiation and development. We hope to unite the existing molecular and microscopic data on genome organization to get a holistic view of the nucleus, and propose a model, in which repeat repertoire together with scaffolding structures blueprint the functional nuclear architecture.

A single three-dimensional chromatin compartment in amphioxus indicates a stepwise evolution of vertebrate Hox bimodal regulation

The HoxA and HoxD gene clusters of jawed vertebrates are organized into bipartite three-dimensional chromatin structures that separate long-range regulatory inputs coming from the anterior and posterior Hox-neighboring regions. This architecture is instrumental in allowing vertebrate Hox genes to pattern disparate parts of the body, including limbs. Almost nothing is known about how these three-dimensional topologies originated. Here we perform extensive 4C-seq profiling of the Hox cluster in embryos of amphioxus, an invertebrate chordate. We find that, in contrast to the architecture in vertebrates, the amphioxus Hox cluster is organized into a single chromatin interaction domain that includes long-range contacts mostly from the anterior side, bringing distant cis-regulatory elements into contact with Hox genes. We infer that the vertebrate Hox bipartite regulatory system is an evolutionary novelty generated by combining ancient long-range regulatory contacts from DNA in the anterior Hox neighborhood with new regulatory inputs from the posterior side.

Visualizing the HoxD Gene Cluster at the Nanoscale Level

Transcription of HoxD cluster genes in limbs is coordinated by two topologically associating domains (TADs), neighboring the cluster and containing various enhancers. Here, we use a combination of microscopy approaches and chromosome conformation capture to assess the structural changes occurring in this global architecture in various functional states. We observed that despite their spatial juxtaposition, the TADs are consistently kept as distinct three-dimensional units. Hox genes located at their boundary can show significant spatial segregation over long distances, suggesting that physical elongation of the HoxD cluster occurs. The use of superresolution imaging (STORM [stochastic optical reconstruction microscopy]) revealed that the gene cluster can be in an either compact or elongated shape. The latter configuration is observed in transcriptionally active tissue and in embryonic stem cells, consistent with chromosome conformation capture results. Such morphological changes at HoxD in developing digits seem to be associated with its position at the boundary between two TADs and support the idea that chromatin dynamics is important in the establishment of transcriptional activity.

Determination of High-Resolution 3D Chromatin Organization Using Circular Chromosome Conformation Capture (4C-seq)

3D chromatin organization is essential for many aspects of transcriptional regulation. Circular Chromosome Conformation Capture followed by Illumina sequencing (4C-seq) is among the most powerful techniques to determine 3D chromatin organization. 4C-seq, like other modifications of the original 3C technique, uses the principle of "proximity ligation" to identify and quantify ten thousands of genomic interactions at a kilobase scale in a single experiment for predefined loci in the genome.In this chapter we focus on the experimental steps in the 4C-seq protocol, providing detailed descriptions on the preparation of cells, the construction of the circularized 3C library and the generation of the Illumina high throughput sequencing library. This protocol is particularly suited for the use of mammalian tissue samples, but can be used with minimal changes on circulating cells and cell lines from other sources as well. In the final section of this chapter, we provide a brief overview of data analysis approaches, accompanied by links to publicly available analysis tools.

Multiplexed analysis of chromosome conformation at vastly improved sensitivity

Methods for analyzing chromosome conformation in mammalian cells are either low resolution or low throughput and are technically challenging. In next-generation (NG) Capture-C, we have redesigned the Capture-C method to achieve unprecedented levels of sensitivity and reproducibility. NG Capture-C can be used to analyze many genetic loci and samples simultaneously. High-resolution data can be produced with as few as 100,000 cells, and single-nucleotide polymorphisms can be used to generate allele-specific tracks. The method is straightforward to perform and should greatly facilitate the investigation of many questions related to gene regulation as well as the functional dissection of traits examined in genome-wide association studies.

C-ing the Genome: A Compendium of Chromosome Conformation Capture Methods to Study Higher-Order Chromatin Organization

Three-dimensional organization of the chromatin has important roles in transcription, replication, DNA repair, and pathologic events such as translocations. There are two fundamental ways to study higher-order chromatin organization: microscopic and molecular approaches. In this review, we briefly introduce the molecular approaches, focusing on chromosome conformation capture or "3C" technology and its derivatives, which can be used to probe chromatin folding at resolutions beyond that provided by microscopy techniques. We further discuss the different types of data generated by the 3C-based methods and how they can be used to answer distinct biological questions.

Profiling of conserved non-coding elements upstream of SHOX and functional characterisation of the SHOX cis-regulatory landscape

Genetic defects such as copy number variations (CNVs) in non-coding regions containing conserved non-coding elements (CNEs) outside the transcription unit of their target gene, can underlie genetic disease. An example of this is the short stature homeobox (SHOX) gene, regulated by seven CNEs located downstream and upstream of SHOX, with proven enhancer capacity in chicken limbs. CNVs of the downstream CNEs have been reported in many idiopathic short stature (ISS) cases, however, only recently have a few CNVs of the upstream enhancers been identified. Here, we set out to provide insight into: (i) the cis-regulatory role of these upstream CNEs in human cells, (ii) the prevalence of upstream CNVs in ISS, and (iii) the chromatin architecture of the SHOX cis-regulatory landscape in chicken and human cells. Firstly, luciferase assays in human U2OS cells, and 4C-seq both in chicken limb buds and human U2OS cells, demonstrated cis-regulatory enhancer capacities of the upstream CNEs. Secondly, CNVs of these upstream CNEs were found in three of 501 ISS patients. Finally, our 4C-seq interaction map of the SHOX region reveals a cis-regulatory domain spanning more than 1 Mb and harbouring putative new cis-regulatory elements.

Regulation of gene transcription by Polycomb proteins

The Polycomb group (PcG) of proteins defines a subset of factors that physically associate and function to maintain the positional identity of cells from the embryo to adult stages. PcG has long been considered a paradigmatic model for epigenetic maintenance of gene transcription programs. Despite intensive research efforts to unveil the molecular mechanisms of action of PcG proteins, several fundamental questions remain unresolved: How many different PcG complexes exist in mammalian cells? How are PcG complexes targeted to specific loci? How does PcG regulate transcription? In this review, we discuss the diversity of PcG complexes in mammalian cells, examine newly identified modes of recruitment to chromatin, and highlight the latest insights into the molecular mechanisms underlying the function of PcGs in transcription regulation and three-dimensional chromatin conformation.

The long noncoding RNA Gm15055 represses Hoxa gene expression by recruiting PRC2 to the gene cluster

The Hox genes encode transcription factors that determine embryonic pattern formation. In embryonic stem cells, the Hox genes are silenced by PRC2. Recent studies have reported a role for long noncoding RNAs in PRC2 recruitment in vertebrates. However, little is known about how PRC2 is recruited to the Hox genes in ESCs. Here, we used stable knockdown and knockout strategies to characterize the function of the long noncoding RNA Gm15055 in the regulation of Hoxa genes in mouse ESCs. We found that Gm15055 is highly expressed in mESCs and its expression is maintained by OCT4. Gm15055 represses Hoxa gene expression by recruiting PRC2 to the cluster and maintaining the H3K27me3 modification on Hoxa promoters. A chromosome conformation capture assay revealed the close physical association of the Gm15055 locus to multiple sites at the Hoxa gene cluster in mESCs, which may facilitate the in cis targeting of Gm15055 RNA to the Hoxa genes. Furthermore, an OCT4-responsive positive cis-regulatory element is found in the Gm15055 gene locus, which potentially regulates both Gm15055 itself and the Hoxa gene activation. This study suggests how PRC2 is recruited to the Hoxa locus in mESCs, and implies an elaborate mechanism for Hoxa gene regulation in mESCs.

Gene regulation during development in the light of topologically associating domains

During embryonic development, complex transcriptional programs govern the precision of gene expression. Many key developmental genes are regulated via cis-regulatory elements that are located far away in the linear genome. How sequences located hundreds of kilobases away from a promoter can influence its activity has been the subject of numerous speculations, which all underline the importance of the 3D-organization of the genome. The recent advent of chromosome conformation capture techniques has put into focus the subdivision of the genome into topologically associating domains (TADs). TADs may influence regulatory activities on multiple levels. The relative invariance of TAD limits across cell types suggests that they may form fixed structural domains that could facilitate and/or confine long-range regulatory interactions. However, most recent studies suggest that interactions within TADs are more variable and dynamic than initially described. Hence, different models are emerging regarding how TADs shape the complex 3D conformations, and thereafter influence the networks of cis-interactions that govern gene expression during development. For further resources related to this article, please visit the WIREs website.

Nanoscale spatial organization of the HoxD gene cluster in distinct transcriptional states

Chromatin condensation plays an important role in the regulation of gene expression. Recently, it was shown that the transcriptional activation of Hoxd genes during vertebrate digit development involves modifications in 3D interactions within and around the HoxD gene cluster. This reorganization follows a global transition from one set of regulatory contacts to another, between two topologically associating domains (TADs) located on either side of the HoxD locus. Here, we use 3D DNA FISH to assess the spatial organization of chromatin at and around the HoxD gene cluster and report that although the two TADs are tightly associated, they appear as spatially distinct units. We measured the relative position of genes within the cluster and found that they segregate over long distances, suggesting that a physical elongation of the HoxD cluster can occur. We analyzed this possibility by super-resolution imaging (STORM) and found that tissues with distinct transcriptional activity exhibit differing degrees of elongation. We also observed that the morphological change of the HoxD cluster in developing digits is associated with its position at the boundary between the two TADs. Such variations in the fine-scale architecture of the gene cluster suggest causal links among its spatial configuration, transcriptional activation, and the flanking chromatin context.

Demystifying the secret mission of enhancers: linking distal regulatory elements to target genes

Enhancers are short regulatory sequences bound by sequence-specific transcription factors and play a major role in the spatiotemporal specificity of gene expression patterns in development and disease. While it is now possible to identify enhancer regions genomewide in both cultured cells and primary tissues using epigenomic approaches, it has been more challenging to develop methods to understand the function of individual enhancers because enhancers are located far from the gene(s) that they regulate. However, it is essential to identify target genes of enhancers not only so that we can understand the role of enhancers in disease but also because this information will assist in the development of future therapeutic options. After reviewing models of enhancer function, we discuss recent methods for identifying target genes of enhancers. First, we describe chromatin structure-based approaches for directly mapping interactions between enhancers and promoters. Second, we describe the use of correlation-based approaches to link enhancer state with the activity of nearby promoters and/or gene expression. Third, we describe how to test the function of specific enhancers experimentally by perturbing enhancer–target relationships using high-throughput reporter assays and genome editing. Finally, we conclude by discussing as yet unanswered questions concerning how enhancers function, how target genes can be identified, and how to distinguish direct from indirect changes in gene expression mediated by individual enhancers.

Polycomb repressive complex PRC1 spatially constrains the mouse embryonic stem cell genome

The Polycomb repressive complexes PRC1 and PRC2 maintain embryonic stem cell (ESC) pluripotency by silencing lineage-specifying developmental regulator genes. Emerging evidence suggests that Polycomb complexes act through controlling spatial genome organization. We show that PRC1 functions as a master regulator of mouse ESC genome architecture by organizing genes in three-dimensional interaction networks. The strongest spatial network is composed of the four Hox gene clusters and early developmental transcription factor genes, the majority of which contact poised enhancers. Removal of Polycomb repression leads to disruption of promoter-promoter contacts in the Hox gene network. In contrast, promoter-enhancer contacts are maintained in the absence of Polycomb repression, with accompanying widespread acquisition of active chromatin signatures at network enhancers and pronounced transcriptional upregulation of network genes. Thus, PRC1 physically constrains developmental transcription factor genes and their enhancers in a silenced but poised spatial network. We propose that the selective release of genes from this spatial network underlies cell fate specification during early embryonic development.

Detection of Short-Range DNA Interactions in Mammalian Cells Using High-Resolution Circular Chromosome Conformation Capture Coupled to Deep Sequencing

DNA interactions shape the genome to physically and functionally connect regulatory elements to their target genes. Studying these interactions is crucial to understanding the molecular mechanisms that regulate gene expression. In this chapter, we present a protocol for high-resolution circular chromosome conformation capture coupled to deep sequencing. This methodology allows to investigate short-range DNA interactions (<100 kbp) and to obtain high-resolution DNA interaction maps of loci. It is a powerful tool to explore how regulatory elements and genes are connected together.

Exploring CTCF and cohesin related chromatin architecture at HOXA gene cluster in primary human fibroblasts

Spatial expression patterns of homeobox (HOX) genes delineate positional identity of primary fibroblasts from different topographic sites. The molecular mechanism underlying the establishing or maintaining of HOX gene expression pattern remains an attractive developmental issue to be addressed. Our previous work suggested a critical role of CTCF/cohesin-mediated higher- order chromatin structure in RA-induced HOXA activation in human teratocarcinoma NT2/D1 cells. This study investigated the recruitment of CTCF and cohesin, and the higher-order chromatin structure of the HOXA locus in fetal lung and adult foreskin fibroblasts, which display complementary HOXA gene expression patterns. Chromatin contacts between the CTCF-binding sites were observed with lower frequency in human foreskin fibroblasts. This observation is consistent with the lower level of cohesin recruitment and 5' HOXA gene expression in the same cells. We also showed that CTCF-binding site A56 (CBSA56) related chromatin structures exhibit the most notable changes in between the two types of cell, and hence may stand for one of the key CTCF-binding sites for cell-type specific chromatin structure organization. Together, these results imply that CTCF/cohesin coordinates HOXA cluster higher-order chromatin structure and expression during development, and provide insight into the relationship between cell-type specific chromatin organization and the spatial collinearity.

Epigenetic mechanisms and boundaries in the regulation of mammalian Hox clusters

Hox gene expression imparts segment identity to body structures along the anterior-posterior axis and is tightly governed by higher order chromatin mechanisms. Chromatin regulatory features of the homeotic complex are best defined in Drosophila melanogaster, where multiple cis-regulatory elements have been identified that ensure collinear Hox gene expression patterns in accordance with their genomic organization. Recent studies focused on delineating the epigenetic features of the vertebrate Hox clusters have helped reveal their dynamic chromatin organization and its impact on gene expression. Enrichment for the 'activating' H3K4me3 and 'repressive' H3K27me3 histone modifications is a particularly strong read-out for transcriptional status and correlates well with the evidence for chromatin loop domain structures and stage specific topological changes at these loci. However, it is not clear how such distinct domains are imposed and regulated independent of each other. Comparative analysis of the chromatin structure and organization of the homeotic gene clusters in fly and mammals is increasingly revealing the functional conservation of chromatin mediated mechanisms. Here we discuss the case for interspersed boundary elements existing within mammalian Hox clusters along with their possible roles and mechanisms of action. Recent studies suggest a role for factors other than the well characterized vertebrate boundary factor CTCF, such as the GAGA binding factor (GAF), in maintaining chromatin domains at the Hox loci. We also present data demonstrating how such regulatory elements may be involved in organizing higher order structure and demarcating active domains of gene expression at the mammalian Hox clusters.

One, Two, Three: Polycomb Proteins Hit All Dimensions of Gene Regulation

Polycomb group (PcG) proteins contribute to the formation and maintenance of a specific repressive chromatin state that prevents the expression of genes in a particular space and time. Polycomb repressive complexes (PRCs) consist of several PcG proteins with specific regulatory or catalytic properties. PRCs are recruited to thousands of target genes, and various recruitment factors, including DNA-binding proteins and non-coding RNAs, are involved in the targeting. PcG proteins contribute to a multitude of biological processes by altering chromatin features at different scales. PcG proteins mediate both biochemical modifications of histone tails and biophysical modifications (e.g., chromatin fiber compaction and three-dimensional (3D) chromatin conformation). Here, we review the role of PcG proteins in nuclear architecture, describing their impact on the structure of the chromatin fiber, on chromatin interactions, and on the spatial organization of the genome in nuclei. Although little is known about the role of plant PcG proteins in nuclear organization, much is known in the animal field, and we highlight similarities and differences in the roles of PcG proteins in 3D gene regulation in plants and animals.

Keeping the gate closed: functions of the polycomb repressive complex PRC2 in development

Plant ontogeny relies on the correct timing and sequence of transitions between individual developmental phases. These are specified by gene expression patterns that are established by the balanced action of activators and repressors. Polycomb repressive complexes (PRCs) represent an evolutionarily conserved system of epigenetic gene repression that governs the establishment and maintenance of cell, tissue and organ identity, contributing to the correct execution of the developmental programs. PRC2 is a four-subunit histone methyltransferase complex that catalyzes trimethylation of lysine 27 on histone H3 (H3K27me3), which contributes to the change of chromatin structure and long-lasting gene repression. Here, we review the composition and molecular function of the different known PRC2 complexes in plants, and focus on the role of PRC2 in mediating the establishment of different developmental phases and transitions between them.

Genome-wide binding and mechanistic analyses of Smchd1-mediated epigenetic regulation

Structural maintenance of chromosomes flexible hinge domain containing 1 (Smchd1) is an epigenetic repressor with described roles in X inactivation and genomic imprinting, but Smchd1 is also critically involved in the pathogenesis of facioscapulohumeral dystrophy. The underlying molecular mechanism by which Smchd1 functions in these instances remains unknown. Our genome-wide transcriptional and epigenetic analyses show that Smchd1 binds cis-regulatory elements, many of which coincide with CCCTC-binding factor (Ctcf) binding sites, for example, the clustered protocadherin (Pcdh) genes, where we show Smchd1 and Ctcf act in opposing ways. We provide biochemical and biophysical evidence that Smchd1-chromatin interactions are established through the homodimeric hinge domain of Smchd1 and, intriguingly, that the hinge domain also has the capacity to bind DNA and RNA. Our results suggest Smchd1 imparts epigenetic regulation via physical association with chromatin, which may antagonize Ctcf-facilitated chromatin interactions, resulting in coordinated transcriptional control.

Defining the Human Brain Proteome Using Transcriptomics and Antibody-Based Profiling with a Focus on the Cerebral Cortex

The mammalian brain is a complex organ composed of many specialized cells, harboring sets of both common, widely distributed, as well as specialized and discretely localized proteins. Here we focus on the human brain, utilizing transcriptomics and public available Human Protein Atlas (HPA) data to analyze brain-enriched (frontal cortex) polyadenylated messenger RNA and long non-coding RNA and generate a genome-wide draft of global and cellular expression patterns of the brain. Based on transcriptomics analysis of altogether 27 tissues, we have estimated that approximately 3% (n=571) of all protein coding genes and 13% (n=87) of the long non-coding genes expressed in the human brain are enriched, having at least five times higher expression levels in brain as compared to any of the other analyzed peripheral tissues. Based on gene ontology analysis and detailed annotation using antibody-based tissue micro array analysis of the corresponding proteins, we found the majority of brain-enriched protein coding genes to be expressed in astrocytes, oligodendrocytes or in neurons with molecular properties linked to synaptic transmission and brain development. Detailed analysis of the transcripts and the genetic landscape of brain-enriched coding and non-coding genes revealed brain-enriched splice variants. Several clusters of neighboring brain-enriched genes were also identified, suggesting regulation of gene expression on the chromatin level. This multi-angle approach uncovered the brain-enriched transcriptome and linked genes to cell types and functions, providing novel insights into the molecular foundation of this highly specialized organ.

Evolutionary comparison reveals that diverging CTCF sites are signatures of ancestral topological associating domains borders

Increasing evidence in the last years indicates that the vast amount of regulatory information contained in mammalian genomes is organized in precise 3D chromatin structures. However, the impact of this spatial chromatin organization on gene expression and its degree of evolutionary conservation is still poorly understood. The Six homeobox genes are essential developmental regulators organized in gene clusters conserved during evolution. Here, we reveal that the Six clusters share a deeply evolutionarily conserved 3D chromatin organization that predates the Cambrian explosion. This chromatin architecture generates two largely independent regulatory landscapes (RLs) contained in two adjacent topological associating domains (TADs). By disrupting the conserved TAD border in one of the zebrafish Six clusters, we demonstrate that this border is critical for preventing competition between promoters and enhancers located in separated RLs, thereby generating different expression patterns in genes located in close genomic proximity. Moreover, evolutionary comparison of Six-associated TAD borders reveals the presence of CCCTC-binding factor (CTCF) sites with diverging orientations in all studied deuterostomes. Genome-wide examination of mammalian HiC data reveals that this conserved CTCF configuration is a general signature of TAD borders, underscoring that common organizational principles underlie TAD compartmentalization in deuterostome evolution.

The role of chromosome domains in shaping the functional genome

The genome must be highly compacted to fit within eukaryotic nuclei but must be accessible to the transcriptional machinery to allow appropriate expression of genes in different cell types and throughout developmental pathways. A growing body of work has shown that the genome, analogously to proteins, forms an ordered, hierarchical structure that closely correlates and may even be causally linked with regulation of functions such as transcription. This review describes our current understanding of how these functional genomic "secondary and tertiary structures" form a blueprint for global nuclear architecture and the potential they hold for understanding and manipulating genomic regulation.

Excessive versus physiologically relevant levels of retinoic acid in embryonic stem cell differentiation

Over the past two decades, embryonic stem cells (ESCs) have been established as a valuable system to study the complex molecular events that underlie the collinear activation of Hox genes during development. When ESCs are induced to differentiate in response to retinoic acid (RA), Hox genes are transcriptionally activated in their chromosomal order, with the most 3' Hox genes activated first, sequentially followed by more 5' Hox genes. In contrast to the low levels of RA detected during gastrulation (∼33 nM), a time when Hox genes are induced during embryonic development, high levels of RA are used to study Hox gene activation in ESCs in vitro (1-10 µM). This compelled us to compare RA-induced ESC differentiation in vitro with Hox gene activation in vivo. In this study, we show that treatment of ESCs for 2 days with RA best mimics activation of Hox genes during embryonic development. Furthermore, we show that defects in Hox gene expression known to occur in embryos lacking the histone acetyltransferase MOZ (also called MYST3 or KAT6A) were masked in Moz-deficient ESCs when excessive RA (0.5-5 µM) was used. The role of MOZ in Hox gene activation was only evident when ESCs were differentiated at low concentrations of RA, namely 20 nM, which is similar to RA levels in vivo. Our results demonstrate that using RA at physiologically relevant levels to study the activation of Hox genes, more accurately reflects the molecular events during the early phase of Hox gene activation in vivo.

Nuclear matrix and structural and functional compartmentalization of the eucaryotic cell nucleus

Becoming popular at the end of the 20th century, the concept of the nuclear matrix implies the existence of a nuclear skeleton that organizes functional elements in the cell nucleus. This review presents a critical analysis of the results obtained in the study of nuclear matrix in the light of current views on the organization of the cell nucleus. Numerous studies of nuclear matrix have failed to provide evidence of the existence of such a structure. Moreover, the existence of a filamentous structure that supports the nuclear compartmentalization appears to be unnecessary, since this function is performed by the folded genome itself.

MOZ and BMI1 play opposing roles during Hox gene activation in ES cells and in body segment identity specification in vivo

Hox genes underlie the specification of body segment identity in the anterior-posterior axis. They are activated during gastrulation and undergo a dynamic shift from a transcriptionally repressed to an active chromatin state in a sequence that reflects their chromosomal location. Nevertheless, the precise role of chromatin modifying complexes during the initial activation phase remains unclear. In the current study, we examined the role of chromatin regulators during Hox gene activation. Using embryonic stem cell lines lacking the transcriptional activator MOZ and the polycomb-family repressor BMI1, we showed that MOZ and BMI1, respectively, promoted and repressed Hox genes during the shift from the transcriptionally repressed to the active state. Strikingly however, MOZ but not BMI1 was required to regulate Hox mRNA levels after the initial activation phase. To determine the interaction of MOZ and BMI1 in vivo, we interrogated their role in regulating Hox genes and body segment identity using Moz;Bmi1 double deficient mice. We found that the homeotic transformations and shifts in Hox gene expression boundaries observed in single Moz and Bmi1 mutant mice were rescued to a wild type identity in Moz;Bmi1 double knockout animals. Together, our findings establish that MOZ and BMI1 play opposing roles during the onset of Hox gene expression in the ES cell model and during body segment identity specification in vivo. We propose that chromatin-modifying complexes have a previously unappreciated role during the initiation phase of Hox gene expression, which is critical for the correct specification of body segment identity.

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

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.

Unraveling the 3D genome: genomics tools for multiscale exploration

A decade of rapid method development has begun to yield exciting insights into the 3D architecture of the metazoan genome and the roles it may play in regulating transcription. Here we review core methods and new tools in the modern genomicist's toolbox at three length scales, ranging from single base pairs to megabase-scale chromosomal domains, and discuss the emerging picture of the 3D genome that these tools have revealed. Blind spots remain, especially at intermediate length scales spanning a few nucleosomes, but thanks in part to new technologies that permit targeted alteration of chromatin states and time-resolved studies, the next decade holds great promise for hypothesis-driven research into the mechanisms that drive genome architecture and transcriptional regulation.

[The intimate genome… in three dimensions]

Over the past decade, techniques based on chromosome conformation capture (3C) have accelerated our understanding of eukaryote's nuclear architecture. Coupled to high throughput sequencing and bioinformatics they have unveiled different organizational levels of the genome at an unprecedented scale. Initially performed using large populations of cells, a new variant of these techniques can be applied to single cell. Although it can be shown that chromosome folding varies from one cell to the other, their overall organization into topologically associating domains is conserved between cells of the same population. Interestingly, the predicted chromosome structures reveal that regions engaged in trans-chromosomal interactions are preferentially localized at the surface of the chromosome territory. These results confirm and extend previous observations on individual loci therefore highlighting the power of 3C based techniques.

CTCF establishes discrete functional chromatin domains at the Hox clusters during differentiation

Polycomb and Trithorax group proteins encode the epigenetic memory of cellular positional identity by establishing inheritable domains of repressive and active chromatin within the Hox clusters. Here we demonstrate that the CCCTC-binding factor (CTCF) functions to insulate these adjacent yet antagonistic chromatin domains during embryonic stem cell differentiation into cervical motor neurons. Deletion of CTCF binding sites within the Hox clusters results in the expansion of active chromatin into the repressive domain. CTCF functions as an insulator by organizing Hox clusters into spatially disjoint domains. Ablation of CTCF binding disrupts topological boundaries such that caudal Hox genes leave the repressed domain and become subject to transcriptional activation. Hence, CTCF is required to insulate facultative heterochromatin from impinging euchromatin to produce discrete positional identities.

YAP regulates the expression of Hoxa1 and Hoxc13 in mouse and human oral and skin epithelial tissues

Yes-associated protein (YAP) is a Hippo signaling transcriptional coactivator that plays pivotal roles in stem cell proliferation, organ size control, and tumor development. The downstream targets of YAP have been shown to be highly context dependent. In this study, we used the embryonic mouse tooth germ as a tool to search for the downstream targets of YAP in ectoderm-derived tissues. Yap deficiency in the dental epithelium resulted in a small tooth germ with reduced epithelial cell proliferation. We compared the gene expression profiles of embryonic day 14.5 (E14.5) Yap conditional knockout and YAP transgenic mouse tooth germs using transcriptome sequencing (RNA-Seq) and further confirmed the differentially expressed genes using real-time PCR and in situ hybridization. We found that YAP regulates the expression of Hoxa1 and Hoxc13 in oral and dental epithelial tissues as well as in the epidermis of skin during embryonic and adult stages. Sphere formation assay suggested that Hoxa1 and Hoxc13 are functionally involved in YAP-regulated epithelial progenitor cell proliferation, and chromatin immunoprecipitation (ChIP) assay implies that YAP may regulate Hoxa1 and Hoxc13 expression through TEAD transcription factors. These results provide mechanistic insights into abnormal YAP activities in mice and humans.

Chromosome 7 gain and DNA hypermethylation at the HOXA10 locus are associated with expression of a stem cell related HOX-signature in glioblastoma

Background

HOX genes are a family of developmental genes that are expressed neither in the developing forebrain nor in the normal brain. Aberrant expression of a HOX-gene dominated stem-cell signature in glioblastoma has been linked with increased resistance to chemo-radiotherapy and sustained proliferation of glioma initiating cells. Here we describe the epigenetic and genetic alterations and their interactions associated with the expression of this signature in glioblastoma.

Results

We observe prominent hypermethylation of the HOXA locus 7p15.2 in glioblastoma in contrast to non-tumoral brain. Hypermethylation is associated with a gain of chromosome 7, a hallmark of glioblastoma, and may compensate for tumor-driven enhanced gene dosage as a rescue mechanism by preventing undue gene expression. We identify the CpG island of the HOXA10 alternative promoter that appears to escape hypermethylation in the HOX-high glioblastoma. An additive effect of gene copy gain at 7p15.2 and DNA methylation at key regulatory CpGs in HOXA10 is significantly associated with HOX-signature expression. Additionally, we show concordance between methylation status and presence of active or inactive chromatin marks in glioblastoma-derived spheres that are HOX-high or HOX-low, respectively.

Conclusions

Based on these findings, we propose co-evolution and interaction between gene copy gain, associated with a gain of chromosome 7, and additional epigenetic alterations as key mechanisms triggering a coordinated, but inappropriate, HOX transcriptional program in glioblastoma.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0583-7) contains supplementary material, which is available to authorized users.

Epigenomic analysis of the HOX gene loci reveals mechanisms that may control canonical expression patterns in AML and normal hematopoietic cells

HOX genes are highly expressed in many acute myeloid leukemia (AML) samples, but the patterns of expression and associated regulatory mechanisms are not clearly understood. We analyzed RNA sequencing data from 179 primary AML samples and normal hematopoietic cells to understand the range of expression patterns in normal versus leukemic cells. HOX expression in AML was restricted to specific genes in the HOXA or HOXB loci, and was highly correlated with recurrent cytogenetic abnormalities. However, the majority of samples expressed a canonical set of HOXA and HOXB genes that was nearly identical to the expression signature of normal hematopoietic stem/progenitor cells. Transcriptional profiles at the HOX loci were similar between normal cells and AML samples, and involved bidirectional transcription at the center of each gene cluster. Epigenetic analysis of a subset of AML samples also identified common regions of chromatin accessibility in AML samples and normal CD34(+) cells that displayed differences in methylation depending on HOX expression patterns. These data provide an integrated epigenetic view of the HOX gene loci in primary AML samples, and suggest that HOX expression in most AML samples represents a normal stem cell program that is controlled by epigenetic mechanisms at specific regulatory elements.

A variably occupied CTCF binding site in the ultrabithorax gene in the Drosophila bithorax complex

Although the majority of genomic binding sites for the insulator protein CCCTC-binding factor (CTCF) are constitutively occupied, a subset show variable occupancy. Such variable sites provide an opportunity to assess context-specific CTCF functions in gene regulation. Here, we have identified a variably occupied CTCF site in the Drosophila Ultrabithorax (Ubx) gene. This site is occupied in tissues where Ubx is active (third thoracic leg imaginal disc) but is not bound in tissues where the Ubx gene is repressed (first thoracic leg imaginal disc). Using chromatin conformation capture, we show that this site preferentially interacts with the Ubx promoter region in the active state. The site lies close to Ubx enhancer elements and is also close to the locations of several gypsy transposon insertions that disrupt Ubx expression, leading to the bx mutant phenotype. gypsy insertions carry the Su(Hw)-dependent gypsy insulator and were found to affect both CTCF binding at the variable site and the chromatin topology. This suggests that insertion of the gypsy insulator in this region interferes with CTCF function and supports a model for the normal function of the variable CTCF site as a chromatin loop facilitator, promoting interaction between Ubx enhancers and the Ubx transcription start site.

Deep conservation of wrist and digit enhancers in fish

There is no obvious morphological counterpart of the autopod (wrist/ankle and digits) in living fishes. Comparative molecular data may provide insight into understanding both the homology of elements and the evolutionary developmental mechanisms behind the fin to limb transition. In mouse limbs the autopod is built by a "late" phase of Hoxd and Hoxa gene expression, orchestrated by a set of enhancers located at the 5' end of each cluster. Despite a detailed mechanistic understanding of mouse limb development, interpretation of Hox expression patterns and their regulation in fish has spawned multiple hypotheses as to the origin and function of "autopod" enhancers throughout evolution. Using phylogenetic footprinting, epigenetic profiling, and transgenic reporters, we have identified and functionally characterized hoxD and hoxA enhancers in the genomes of zebrafish and the spotted gar, Lepisosteus oculatus, a fish lacking the whole genome duplication of teleosts. Gar and zebrafish "autopod" enhancers drive expression in the distal portion of developing zebrafish pectoral fins, and respond to the same functional cues as their murine orthologs. Moreover, gar enhancers drive reporter gene expression in both the wrist and digits of mouse embryos in patterns that are nearly indistinguishable from their murine counterparts. These functional genomic data support the hypothesis that the distal radials of bony fish are homologous to the wrist and/or digits of tetrapods.

Spatial genome organization: contrasting views from chromosome conformation capture and fluorescence in situ hybridization

Although important for gene regulation, most studies of genome organization use either fluorescence in situ hybridization (FISH) or chromosome conformation capture (3C) methods. FISH directly visualizes the spatial relationship of sequences but is usually applied to a few loci at a time. The frequency at which sequences are ligated together by formaldehyde cross-linking can be measured genome-wide by 3C methods, with higher frequencies thought to reflect shorter distances. FISH and 3C should therefore give the same views of genome organization, but this has not been tested extensively. We investigated the murine HoxD locus with 3C carbon copy (5C) and FISH in different developmental and activity states and in the presence or absence of epigenetic regulators. We identified situations in which the two data sets are concordant but found other conditions under which chromatin topographies extrapolated from 5C or FISH data are not compatible. We suggest that products captured by 3C do not always reflect spatial proximity, with ligation occurring between sequences located hundreds of nanometers apart, influenced by nuclear environment and chromatin composition. We conclude that results obtained at high resolution with either 3C methods or FISH alone must be interpreted with caution and that views about genome organization should be validated by independent methods.