Epigenetic regulatory layers in the 3D nucleus.

Nearly 7 decades have elapsed since Francis Crick introduced the central dogma of molecular biology, as part of his ideas on protein synthesis, setting the fundamental rules of sequence information transfer from DNA to RNAs and proteins. We have since learned that gene expression is finely tuned in time and space, due to the activities of RNAs and proteins on regulatory DNA elements, and through cell-type-specific three-dimensional conformations of the genome. Here, we review major advances in genome biology and discuss a set of ideas on gene regulation and highlight how various biomolecular assemblies lead to the formation of structural and regulatory features within the nucleus, with roles in transcriptional control. We conclude by suggesting further developments that will help capture the complex, dynamic, and often spatially restricted events that govern gene expression in mammalian cells.

Three-dimensional genome architecture: players and mechanisms.

The different cell types of an organism share the same DNA, but during cell differentiation their genomes undergo diverse structural and organizational changes that affect gene expression and other cellular functions. These can range from large-scale folding of whole chromosomes or of smaller genomic regions, to the re-organization of local interactions between enhancers and promoters, mediated by the binding of transcription factors and chromatin looping. The higher-order organization of chromatin is also influenced by the specificity of the contacts that it makes with nuclear structures such as the lamina. Sophisticated methods for mapping chromatin contacts are generating genome-wide data that provide deep insights into the formation of chromatin interactions, and into their roles in the organization and function of the eukaryotic cell nucleus.

Cell-type specialization is encoded by specific chromatin topologies.

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.

R-Loops Enhance Polycomb Repression at a Subset of Developmental Regulator Genes.

R-loops are three-stranded nucleic acid structures that form during transcription, especially over unmethylated CpG-rich promoters of active genes. In mouse embryonic stem cells (mESCs), CpG-rich developmental regulator genes are repressed by the Polycomb complexes PRC1 and PRC2. Here, we show that R-loops form at a subset of Polycomb target genes, and we investigate their contribution to Polycomb repression. At R-loop-positive genes, R-loop removal leads to decreased PRC1 and PRC2 recruitment and Pol II activation into a productive elongation state, accompanied by gene derepression at nascent and processed transcript levels. Stable removal of PRC2 derepresses R-loop-negative genes, as expected, but does not affect R-loops, PRC1 recruitment, or transcriptional repression of R-loop-positive genes. Our results highlight that Polycomb repression does not occur via one mechanism but consists of different layers of repression, some of which are gene specific. We uncover that one such mechanism is mediated by an interplay between R-loops and RING1B recruitment.

Multiplex-GAM: genome-wide identification of chromatin contacts yields insights overlooked by Hi-C.

Technology for measuring 3D genome topology is increasingly important for studying gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of genome architecture mapping (GAM), a ligation-free technique that maps chromatin contacts genome-wide. We perform a detailed comparison of multiplex-GAM and Hi-C using mouse embryonic stem cells. When examining the strongest contacts detected by either method, we find that only one-third of these are shared. The strongest contacts specifically found in GAM often involve 'active' regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain 'inactive' regions. Our work shows that active genomic regions are involved in extensive complex contacts that are currently underestimated in ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies.

Methods for mapping 3D chromosome architecture.

Determining how chromosomes are positioned and folded within the nucleus is critical to understanding the role of chromatin topology in gene regulation. Several methods are available for studying chromosome architecture, each with different strengths and limitations. Established imaging approaches and proximity ligation-based chromosome conformation capture (3C) techniques (such as DNA-FISH and Hi-C, respectively) have revealed the existence of chromosome territories, functional nuclear landmarks (such as splicing speckles and the nuclear lamina) and topologically associating domains. Improvements to these methods and the recent development of ligation-free approaches, including GAM, SPRITE and ChIA-Drop, are now helping to uncover new aspects of 3D genome topology that confirm the nucleus to be a complex, highly organized organelle.

LifeTime and improving European healthcare through cell-based interceptive medicine.

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade.

A Finer Print Than TADs: PRC1-Mediated Domains.

Polycomb proteins are well-known epigenetic repressors with unexplained roles in chromatin folding. In this issue of Molecular Cell, Kundu et al. (2017) investigate the structures of PRC1-mediated domains in stem cells and probe their changes upon differentiation and in PRC knockouts.

Comparison of the Hi-C, GAM and SPRITE methods using polymer models of chromatin.

Hi-C, split-pool recognition of interactions by tag extension (SPRITE) and genome architecture mapping (GAM) are powerful technologies utilized to probe chromatin interactions genome wide, but how faithfully they capture three-dimensional (3D) contacts and how they perform relative to each other is unclear, as no benchmark exists. Here, we compare these methods in silico in a simplified, yet controlled, framework against known 3D structures of polymer models of murine and human loci, which can recapitulate Hi-C, GAM and SPRITE experiments and multiplexed fluorescence in situ hybridization (FISH) single-molecule conformations. We find that in silico Hi-C, GAM and SPRITE bulk data are faithful to the reference 3D structures whereas single-cell data reflect strong variability among single molecules. The minimal number of cells required in replicate experiments to return statistically similar contacts is different across the technologies, being lowest in SPRITE and highest in GAM under the same conditions. Noise-to-signal levels follow an inverse power law with detection efficiency and grow with genomic distance differently among the three methods, being lowest in GAM for genomic separations >1 Mb.

Complex multi-enhancer contacts captured by genome architecture mapping.

The organization of the genome in the nucleus and the interactions of genes with their regulatory elements are key features of transcriptional control and their disruption can cause disease. Here we report a genome-wide method, genome architecture mapping (GAM), for measuring chromatin contacts and other features of three-dimensional chromatin topology on the basis of sequencing DNA from a large collection of thin nuclear sections. We apply GAM to mouse embryonic stem cells and identify enrichment for specific interactions between active genes and enhancers across very large genomic distances using a mathematical model termed SLICE (statistical inference of co-segregation). GAM also reveals an abundance of three-way contacts across the genome, especially between regions that are highly transcribed or contain super-enhancers, providing a level of insight into genome architecture that, owing to the technical limitations of current technologies, has previously remained unattainable. Furthermore, GAM highlights a role for gene-expression-specific contacts in organizing the genome in mammalian nuclei.

GAMIBHEAR: whole-genome haplotype reconstruction from Genome Architecture Mapping data.

MOTIVATION: Genome Architecture Mapping (GAM) was recently introduced as a digestion- and ligation-free method to detect chromatin conformation. Orthogonal to existing approaches based on chromatin conformation capture (3C), GAM's ability to capture both inter- and intra-chromosomal contacts from low amounts of input data makes it particularly well suited for allele-specific analyses in a clinical setting. Allele-specific analyses are powerful tools to investigate the effects of genetic variants on many cellular phenotypes including chromatin conformation, but require the haplotypes of the individuals under study to be known a priori. So far, however, no algorithm exists for haplotype reconstruction and phasing of genetic variants from GAM data, hindering the allele-specific analysis of chromatin contact points in non-model organisms or individuals with unknown haplotypes. RESULTS: We present GAMIBHEAR, a tool for accurate haplotype reconstruction from GAM data. GAMIBHEAR aggregates allelic co-observation frequencies from GAM data and employs a GAM-specific probabilistic model of haplotype capture to optimize phasing accuracy. Using a hybrid mouse embryonic stem cell line with known haplotype structure as a benchmark dataset, we assess correctness and completeness of the reconstructed haplotypes, and demonstrate the power of GAMIBHEAR to infer accurate genome-wide haplotypes from GAM data. AVAILABILITY AND IMPLEMENTATION: GAMIBHEAR is available as an R package under the open-source GPL-2 license at https://bitbucket.org/schwarzlab/gamibhear. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Formation of new chromatin domains determines pathogenicity of genomic duplications.

Chromosome conformation capture methods have identified subchromosomal structures of higher-order chromatin interactions called topologically associated domains (TADs) that are separated from each other by boundary regions. By subdividing the genome into discrete regulatory units, TADs restrict the contacts that enhancers establish with their target genes. However, the mechanisms that underlie partitioning of the genome into TADs remain poorly understood. Here we show by chromosome conformation capture (capture Hi-C and 4C-seq methods) that genomic duplications in patient cells and genetically modified mice can result in the formation of new chromatin domains (neo-TADs) and that this process determines their molecular pathology. Duplications of non-coding DNA within the mouse Sox9 TAD (intra-TAD) that cause female to male sex reversal in humans, showed increased contact of the duplicated regions within the TAD, but no change in the overall TAD structure. In contrast, overlapping duplications that extended over the next boundary into the neighbouring TAD (inter-TAD), resulted in the formation of a new chromatin domain (neo-TAD) that was isolated from the rest of the genome. As a consequence of this insulation, inter-TAD duplications had no phenotypic effect. However, incorporation of the next flanking gene, Kcnj2, in the neo-TAD resulted in ectopic contacts of Kcnj2 with the duplicated part of the Sox9 regulatory region, consecutive misexpression of Kcnj2, and a limb malformation phenotype. Our findings provide evidence that TADs are genomic regulatory units with a high degree of internal stability that can be sculptured by structural genomic variations. This process is important for the interpretation of copy number variations, as these variations are routinely detected in diagnostic tests for genetic disease and cancer. This finding also has relevance in an evolutionary setting because copy-number differences are thought to have a crucial role in the evolution of genome complexity.

Altered nuclear architecture in blood cells from Huntington's disease patients.

BACKGROUND: Cell nuclear architecture has been explored in cancer and laminopathies but not in neurodegenerative disorders. Huntington's disease (HD) is a neurodegenerative disorder that leads to neuronal death. Chromosome-wide changes in gene expression have been reported in HD, not only in the brain but also in peripheral blood cells, but whether this translates into nuclear and chromosome architecture alterations has not yet been studied. METHODS: We investigate nuclear structure and chromosome organization in HD blood cells using fluorescence in situ hybridization in ultrathin cryosections (cryoFISH), coupled with machine learning image analysis to evaluate size, distribution, and morphology of nuclei and chromosomes. Four chromosomes were analyzed based on up- or downregulation of gene expression in HD. RESULTS: We show that blood cells from HD patients display increased nuclear size and filamentary shape, increased size of gene-rich chromosome 19, decreased filamentary shape of gene-rich chromosome 22, and a more radially centralized position for chromosome 19, whereas chromosomes 4 and 5 do not show detectable differences. CONCLUSIONS: We identify gross changes in nuclear architecture and chromosome organization associated with HD in blood. This adds a new layer of information onto disrupting mechanisms in HD and increases the potential of using blood to survey HD.

Potential use of macroalgae Gracilaria gracilis in diets for European seabass (Dicentrarchus labrax): Health benefits from a sustainable source.

Seaweeds still possess a large undisclosed potential, mainly due to their constituent's richness, which may have several uses for society. In aquaculture, they may play a role as an ecological sustainable aquafeed supplement to increase overall health and fight pathogenic outbreaks. This study aimed to evaluate the general health modulation that the inclusion of Gracilaria gracilis could accomplish in the diet of Dicentrarchus labrax. Dried algae at 2.5% and 5% and algal extract at 0.35% inclusion levels were supplemented to seabass diet to evaluate possible growth, haematological, immunological, antioxidant, metabolic, and intestinal morphological modulations. The supplementations did not impact growth or feed utilization, and barely affected the haematological profile and some metabolic parameters. Nevertheless, it caused a marked outcome on lysozyme, some oxidative stress biomarkers, and intestine morphology, suggesting beneficial consequences from the algal inclusion. Dried algae powder, with a 2.5% inclusion, boosted immune response, with higher plasmatic lysozyme and intestinal acid goblet cells and protected against oxidative damages by improved enzymatic and non-enzymatic responses. Thus, we provide evidence that dietary seaweed application may be a path towards a more sustainable aquaculture industry.

Inference of chromosome 3D structures from GAM data by a physics computational approach.

The combination of modelling and experimental advances can provide deep insights for understanding chromatin 3D organization and ultimately its underlying mechanisms. In particular, models of polymer physics can help comprehend the complexity of genomic contact maps, as those emerging from technologies such as Hi-C, GAM or SPRITE. Here we discuss a method to reconstruct 3D structures from Genome Architecture Mapping (GAM) data, based on PRISMR, a computational approach introduced to find the minimal polymer model best describing Hi-C input data from only polymer physics. After recapitulating the PRISMR procedure, we describe how we extended it for treating GAM data. We successfully test the method on a 6 Mb region around the Sox9 gene and, at a lower resolution, on the whole chromosome 7 in mouse embryonic stem cells. The PRISMR derived 3D structures from GAM co-segregation data are finally validated against independent Hi-C contact maps. The method results to be versatile and robust, hinting that it can be similarly applied to different experimental data, such as SPRITE or microscopy distance data.

RNA polymerase II primes Polycomb-repressed developmental genes throughout terminal neuronal differentiation.

Polycomb repression in mouse embryonic stem cells (ESCs) is tightly associated with promoter co-occupancy of RNA polymerase II (RNAPII) which is thought to prime genes for activation during early development. However, it is unknown whether RNAPII poising is a general feature of Polycomb repression, or is lost during differentiation. Here, we map the genome-wide occupancy of RNAPII and Polycomb from pluripotent ESCs to non-dividing functional dopaminergic neurons. We find that poised RNAPII complexes are ubiquitously present at Polycomb-repressed genes at all stages of neuronal differentiation. We observe both loss and acquisition of RNAPII and Polycomb at specific groups of genes reflecting their silencing or activation. Strikingly, RNAPII remains poised at transcription factor genes which are silenced in neurons through Polycomb repression, and have major roles in specifying other, non-neuronal lineages. We conclude that RNAPII poising is intrinsically associated with Polycomb repression throughout differentiation. Our work suggests that the tight interplay between RNAPII poising and Polycomb repression not only instructs promoter state transitions, but also may enable promoter plasticity in differentiated cells.

Short RNAs are transcribed from repressed polycomb target genes and interact with polycomb repressive complex-2.

Polycomb proteins maintain cell identity by repressing the expression of developmental regulators specific for other cell types. Polycomb repressive complex-2 (PRC2) catalyzes trimethylation of histone H3 lysine-27 (H3K27me3). Although repressed, PRC2 targets are generally associated with the transcriptional initiation marker H3K4me3, but the significance of this remains unclear. Here, we identify a class of short RNAs, approximately 50-200 nucleotides in length, transcribed from the 5' end of polycomb target genes in primary T cells and embryonic stem cells. Short RNA transcription is associated with RNA polymerase II and H3K4me3, occurs in the absence of mRNA transcription, and is independent of polycomb activity. Short RNAs form stem-loop structures resembling PRC2 binding sites in Xist, interact with PRC2 through SUZ12, cause gene repression in cis, and are depleted from polycomb target genes activated during cell differentiation. We propose that short RNAs play a role in the association of PRC2 with its target genes.