Chromatin architecture of the human genome: gene-rich domains are enriched in open chromatin fibers

Gene-Specific H1 Eviction through a Transcriptional Activator→p300→NAP1→H1 Pathway

Linker histone H1 has been correlated with transcriptional inhibition, but the mechanistic basis of the inhibition and its reversal during gene activation has remained enigmatic. We report that H1-compacted chromatin, reconstituted in vitro, blocks transcription by abrogating core histone modifications by p300 but not activator and p300 binding. Transcription from H1-bound chromatin is elicited by the H1 chaperone NAP1, which is recruited in a gene-specific manner through direct interactions with activator-bound p300 that facilitate core histone acetylation (by p300) and concomitant eviction of H1 and H2A-H2B. An analysis in B cells confirms the strong dependency on NAP1-mediated H1 eviction for induction of the silent CD40 gene and further demonstrates that H1 eviction, seeded by activator-p300-NAP1-H1 interactions, is propagated over a CCCTC-binding factor (CTCF)-demarcated region through a distinct mechanism that also involves NAP1. Our results confirm direct transcriptional inhibition by H1 and establish a gene-specific H1 eviction mechanism through an activator→p300→NAP1→H1 pathway.

Chromosome territories and the global regulation of the genome

Spatial positioning is a fundamental principle governing nuclear processes. Chromatin is organized as a hierarchy from nucleosomes to Mbp chromatin domains (CD) or topologically associating domains (TADs) to higher level compartments culminating in chromosome territories (CT). Microscopic and sequencing techniques have substantiated chromatin organization as a critical factor regulating gene expression. For example, enhancers loop back to interact with their target genes almost exclusively within TADs, distally located coregulated genes reposition into common transcription factories upon activation, and Mbp CDs exhibit dynamic motion and configurational changes in vivo. A longstanding question in the nucleus field is whether an interactive nuclear matrix provides a direct link between structure and function. The findings of nonrandom radial positioning of CT within the nucleus suggest the possibility of preferential interaction patterns among populations of CT. Sequential labeling up to 10 CT followed by application of computer imaging and geometric graph mining algorithms revealed cell-type specific interchromosomal networks (ICN) of CT that are altered during the cell cycle, differentiation, and cancer progression. It is proposed that the ICN correlate with the global level of genome regulation. These approaches also demonstrated that the large scale 3-D topology of CT is specific for each CT. The cell-type specific proximity of certain chromosomal regions in normal cells may explain the propensity of distinct translocations in cancer subtypes. Understanding how genes are dysregulated upon disruption of the normal "wiring" of the nucleus by translocations, deletions, and amplifications that are hallmarks of cancer, should enable more targeted therapeutic strategies.

A common genomic code for chromatin architecture and recombination landscape

Recent findings established a link between DNA sequence composition and interphase chromatin architecture and explained the evolutionary conservation of TADs (Topologically Associated Domains) and LADs (Lamina Associated Domains) in mammals. This prompted us to analyse conformation capture and recombination rate data to study the relationship between chromatin architecture and recombination landscape of human and mouse genomes. The results reveal that: (1) low recombination domains and blocks of elevated linkage disequilibrium tend to coincide with TADs and isochores, indicating co-evolving regulatory elements and genes in insulated neighbourhoods; (2) double strand break (DSB) and recombination frequencies increase in the short loops of GC-rich TADs, whereas recombination cold spots are typical of LADs and (3) the binding and loading of proteins, which are critical for DSB and meiotic recombination (SPO11, DMC1, H3K4me3 and PRMD9) are higher in GC-rich TADs. One explanation for these observations is that the occurrence of DSB and recombination in meiotic cells are associated with compositional and epigenetic features (genomic code) that influence DNA stiffness/flexibility and appear to be similar to those guiding the chromatin architecture in the interphase nucleus of pre-leptotene cells.

Diverse Genome Topologies Characterize Dosage Compensation across Species

Dosage compensation is the process by which transcript levels of the X chromosome are equalized with those of autosomes. Although diverse mechanisms of dosage compensation have evolved across species, these mechanisms all involve distinguishing the X chromosome from autosomes. Because one chromosome is singled out from other chromosomes for precise regulation, dosage compensation serves as an important model for understanding how specific cis-elements are identified within the highly compacted 3D genome to co-regulate thousands of genes. Recently, multiple genomic approaches have provided key insights into the mechanisms of dosage compensation, extending what we have learned from classical genetic studies. In the future, newer genomic approaches that require little starting material show great promise to provide an understanding of the heterogeneity of dosage compensation between cells and how it functions in nonmodel organisms.

Characterizing the nuclease accessibility of DNA in human cells to map higher order structures of chromatin

Packaging of DNA into chromatin regulates DNA accessibility and consequently all DNA-dependent processes. The nucleosome is the basic packaging unit of DNA forming arrays that are suggested, by biochemical studies, to fold hierarchically into ordered higher-order structures of chromatin. This organization has been recently questioned using microscopy techniques, proposing an irregular structure. To address the principles of chromatin organization, we applied an in situ differential MNase-seq strategy and analyzed in silico the results of complete and partial digestions of human chromatin. We investigated whether different levels of chromatin packaging exist in the cell. We assessed the accessibility of chromatin within distinct domains of kb to Mb genomic regions, performed statistical analyses and computer modelling. We found no difference in MNase accessibility, suggesting no difference in fiber folding between domains of euchromatin and heterochromatin or between other sequence and epigenomic features of chromatin. Thus, our data suggests the absence of differentially organized domains of higher-order structures of chromatin. Moreover, we identified only local structural changes, with individual hyper-accessible nucleosomes surrounding regulatory elements, such as enhancers and transcription start sites. The regulatory sites per se are occupied with structurally altered nucleosomes, exhibiting increased MNase sensitivity. Our findings provide biochemical evidence that supports an irregular model of large-scale chromatin organization.

Physical principles of retroviral integration in the human genome

Certain retroviruses, including HIV, insert their DNA in a non-random fraction of the host genome via poorly understood selection mechanisms. Here, we develop a biophysical model for retroviral integration as stochastic and quasi-equilibrium topological reconnections between polymers. We discover that physical effects, such as DNA accessibility and elasticity, play important and universal roles in this process. Our simulations predict that integration is favoured within nucleosomal and flexible DNA, in line with experiments, and that these biases arise due to competing energy barriers associated with DNA deformations. By considering a long chromosomal region in human T-cells during interphase, we discover that at these larger scales integration sites are predominantly determined by chromatin accessibility. Finally, we propose and solve a reaction-diffusion problem that recapitulates the distribution of HIV hot-spots within T-cells. With few generic assumptions, our model can rationalise experimental observations and identifies previously unappreciated physical contributions to retroviral integration site selection.

H3K9me3-heterochromatin loss at protein-coding genes enables developmental lineage specification

Gene silencing by chromatin compaction is integral to establishing and maintaining cell fates. Trimethylated histone 3 lysine 9 (H3K9me3)-marked heterochromatin is reduced in embryonic stem cells compared to differentiated cells. However, the establishment and dynamics of closed regions of chromatin at protein-coding genes, in embryologic development, remain elusive. We developed an antibody-independent method to isolate and map compacted heterochromatin from low-cell number samples. We discovered high levels of compacted heterochromatin, H3K9me3-decorated, at protein-coding genes in early, uncommitted cells at the germ-layer stage, undergoing profound rearrangements and reduction upon differentiation, concomitant with cell type-specific gene expression. Perturbation of the three H3K9me3-related methyltransferases revealed a pivotal role for H3K9me3 heterochromatin during lineage commitment at the onset of organogenesis and for lineage fidelity maintenance.

Stem cell fate in cancer growth, progression and therapy resistance

Although we have come a long way in our understanding of the signals that drive cancer growth, and how these signals can be targeted, effective control of this disease remains a key scientific and medical challenge. The therapy resistance and relapse that are commonly seen are driven in large part by the inherent heterogeneity within cancers that allows drugs to effectively eliminate some, but not all, malignant cells. Here, we focus on the fundamental drivers of this heterogeneity by examining emerging evidence that shows that these traits are often controlled by the disruption of normal cell fate and aberrant adoption of stem cell signals. We discuss how undifferentiated cells are preferentially primed for transformation and often serve as the cell of origin for cancers. We also consider evidence showing that activation of stem cell programmes in cancers can lead to progression, therapy resistance and metastatic growth and that targeting these attributes may enable better control over a difficult disease.

Polymer Simulations of Heteromorphic Chromatin Predict the 3D Folding of Complex Genomic Loci

Chromatin folded into 3D macromolecular structures is often analyzed by chromosome conformation capture (3C) and fluorescence in situ hybridization (FISH) techniques, but these frequently provide contradictory results. Chromatin can be modeled as a simple polymer composed of a connected chain of units. By embedding data for epigenetic marks (H3K27ac), chromatin accessibility (assay for transposase-accessible chromatin using sequencing [ATAC-seq]), and structural anchors (CCCTC-binding factor [CTCF]), we developed a highly predictive heteromorphic polymer (HiP-HoP) model, where the chromatin fiber varied along its length; combined with diffusing protein bridges and loop extrusion, this model predicted the 3D organization of genomic loci at a population and single-cell level. The model was validated at several gene loci, including the complex Pax6 gene, and was able to determine locus conformations across cell types with varying levels of transcriptional activity and explain different mechanisms of enhancer use. Minimal a priori knowledge of epigenetic marks is sufficient to recapitulate complex genomic loci in 3D and enable predictions of chromatin folding paths.

•

HiP-HoP: highly predictive heteromorphic polymer model to analyze chromatin structure

•

Polymer simulations use widely available epigenetic and protein binding data as input

•

Validate HiP-HoP model at complex loci using 3D FISH and Capture-C

•

Simulations uncover striking conformational variability in chromatin fiber folding

Nuclear cytometry and chromatin organization

The nuclear-targeting chemical probe, for the detection and quantification of DNA within cells, has been a mainstay of cytometry-from the colorimetric Feulgen stain to smart fluorescent agents with tuned functionality. The level of nuclear structure and function at which the probe aims to readout, or indeed at which a DNA-targeted drug acts, is shadowed by a wide range of detection modalities and analytical methods. These methods are invariably limited in terms of the resolution attainable versus the volume occupied by targeted chromatin structures. The scalar challenge arises from the need to understand the extent and different levels of compaction of genomic DNA and how such structures can be re-modeled, reported, or even perturbed by both probes and drugs. Nuclear cytometry can report on the complex levels of chromatin order, disorder, disassembly, and even active disruption by probes and drugs. Nuclear probes can report defining features of clinical and therapeutic interest as in NETosis and other cell death processes. New cytometric approaches continue to bridge the scalar challenges of analyzing chromatin organization. Advances in super-resolution microscopy address the resolution and depth of analysis issues in cellular systems. Typical of recent insights into chromatin organization enabled by exploiting a DNA interacting probe is ChromEM tomography (ChromEMT). ChromEMT uses the unique properties of the anthraquinone-based cytometric dye DRAQ5™ to reveal that local and global 3D chromatin structures effect differences in compaction. The focus of this review is nuclear and chromatin cytometry, with linked reference to DNA targeting probes and drugs as exemplified by the anthracenediones.

A promoter interaction map for cardiovascular disease genetics

Over 500 genetic loci have been associated with risk of cardiovascular diseases (CVDs); however, most loci are located in gene-distal non-coding regions and their target genes are not known. Here, we generated high-resolution promoter capture Hi-C (PCHi-C) maps in human induced pluripotent stem cells (iPSCs) and iPSC-derived cardiomyocytes (CMs) to provide a resource for identifying and prioritizing the functional targets of CVD associations. We validate these maps by demonstrating that promoters preferentially contact distal sequences enriched for tissue-specific transcription factor motifs and are enriched for chromatin marks that correlate with dynamic changes in gene expression. Using the CM PCHi-C map, we linked 1999 CVD-associated SNPs to 347 target genes. Remarkably, more than 90% of SNP-target gene interactions did not involve the nearest gene, while 40% of SNPs interacted with at least two genes, demonstrating the importance of considering long-range chromatin interactions when interpreting functional targets of disease loci.

Epigenetics and epitranscriptomics in temporal patterning of cortical neural progenitor competence

Yoon et al. review epigenetic and epitranscriptomic mechanisms that regulate the lineage specification of neural progenitor cells in the developing brain.

During embryonic brain development, neural progenitor/stem cells (NPCs) sequentially give rise to different subtypes of neurons and glia via a highly orchestrated process. To accomplish the ordered generation of distinct progenies, NPCs go through multistep transitions of their developmental competence. The molecular mechanisms driving precise temporal coordination of these transitions remains enigmatic. Epigenetic regulation, including changes in chromatin structures, DNA methylation, and histone modifications, has been extensively investigated in the context of cortical neurogenesis. Recent studies of chemical modifications on RNA, termed epitranscriptomics, have also revealed their critical roles in neural development. In this review, we discuss advances in understanding molecular regulation of the sequential lineage specification of NPCs in the embryonic mammalian brain with a focus on epigenetic and epitranscriptomic mechanisms. In particular, the discovery of lineage-specific gene transcripts undergoing rapid turnover in NPCs suggests that NPC developmental fate competence is determined much earlier, before the final cell division, and is more tightly controlled than previously appreciated. We discuss how multiple regulatory systems work in harmony to coordinate NPC behavior and summarize recent findings in the context of a model of epigenetic and transcriptional prepatterning to explain NPC developmental competence.

ChromoTrace: Computational reconstruction of 3D chromosome configurations for super-resolution microscopy

The 3D structure of chromatin plays a key role in genome function, including gene expression, DNA replication, chromosome segregation, and DNA repair. Furthermore the location of genomic loci within the nucleus, especially relative to each other and nuclear structures such as the nuclear envelope and nuclear bodies strongly correlates with aspects of function such as gene expression. Therefore, determining the 3D position of the 6 billion DNA base pairs in each of the 23 chromosomes inside the nucleus of a human cell is a central challenge of biology. Recent advances of super-resolution microscopy in principle enable the mapping of specific molecular features with nanometer precision inside cells. Combined with highly specific, sensitive and multiplexed fluorescence labeling of DNA sequences this opens up the possibility of mapping the 3D path of the genome sequence in situ. Here we develop computational methodologies to reconstruct the sequence configuration of all human chromosomes in the nucleus from a super-resolution image of a set of fluorescent in situ probes hybridized to the genome in a cell. To test our approach, we develop a method for the simulation of DNA in an idealized human nucleus. Our reconstruction method, ChromoTrace, uses suffix trees to assign a known linear ordering of in situ probes on the genome to an unknown set of 3D in-situ probe positions in the nucleus from super-resolved images using the known genomic probe spacing as a set of physical distance constraints between probes. We find that ChromoTrace can assign the 3D positions of the majority of loci with high accuracy and reasonable sensitivity to specific genome sequences. By simulating appropriate spatial resolution, label multiplexing and noise scenarios we assess our algorithms performance. Our study shows that it is feasible to achieve genome-wide reconstruction of the 3D DNA path based on super-resolution microscopy images.

The 3D structure of DNA in the nucleus is known to be important for many aspects of DNA function, such as how gene expression is regulated. However, current techniques to localise or determine 3D DNA structure are often indirect. The advent of super-resolution microscopy, at a resolution of 20 nm or better can directly visualize fluorescent probes bound to specific DNA in the nucleus. However it is not trivial to associate how many specific stretches of DNA lie relative to each other, making reliable and precise 3D mapping of large stretches of the genome difficult. Here, we propose a method that leverages the fact that we know the sequence of the genome and the resolution of the super-resolution microscope. Our method, ChromoTrace, uses a computer science data structure, suffix trees, that allow one to simultaneous search the entire genome for specific sub-sequences. To show that our method works, we build a simulation scheme for simulating DNA as ensembles of polymer chains in a nucleus and explore the sensitivity of our method to different types of error. ChromoTrace can robustly and accurately reconstruct 3D paths in our simulations.

Genomic and Proteomic Resolution of Heterochromatin and Its Restriction of Alternate Fate Genes

Heterochromatin is integral to cell identity maintenance by impeding the activation of genes for alternate cell fates. Heterochromatic regions are associated with histone 3 lysine 9 trimethylation (H3K9me3) or H3K27me3, but these modifications are also found in euchromatic regions that permit transcription. We discovered that resistance to sonication is a reliable indicator of the heterochromatin state, and we developed a biophysical method (gradient-seq) to discriminate subtypes of H3K9me3 and H3K27me3 domains in sonication-resistant heterochromatin (srHC) versus euchromatin. These classifications are more accurate than the histone marks alone in predicting transcriptional silence and resistance of alternate fate genes to activation during direct cell conversion. Our proteomics of H3K9me3-marked srHC and functional screens revealed diverse proteins, including RBMX and RBMXL1, that impede gene induction during cellular reprogramming. Isolation of srHC with gradient-seq provides a genome-wide map of chromatin structure, elucidating subtypes of repressed domains that are uniquely predictive of diverse other chromatin properties.

Regulation of chromatin folding by conformational variations of nucleosome linker DNA

Linker DNA conformational variability has been proposed to direct nucleosome array folding into more or less compact chromatin fibers but direct experimental evidence for such models are lacking. Here, we tested this hypothesis by designing nucleosome arrays with A-tracts at specific locations in the nucleosome linkers to induce inward (AT-IN) and outward (AT-OUT) bending of the linker DNA. Using electron microscopy and analytical centrifugation techniques, we observed spontaneous folding of AT-IN nucleosome arrays into highly compact structures, comparable to those induced by linker histone H1. In contrast, AT-OUT nucleosome arrays formed less compact structures with decreased nucleosome interactions similar to wild-type nucleosome arrays. Adding linker histone H1 further increased compaction of the A-tract arrays while maintaining structural differences between them. Furthermore, restriction nuclease digestion revealed a strongly reduced accessibility of nucleosome linkers in the compact AT-IN arrays. Electron microscopy analysis and 3D computational Monte Carlo simulations are consistent with a profound zigzag linker DNA configuration and closer nucleosome proximity in the AT-IN arrays due to inward linker DNA bending. We propose that the evolutionary preferred positioning of A-tracts in DNA linkers may control chromatin higher-order folding and thus influence cellular processes such as gene expression, transcription and DNA repair.

3D genome structure modeling by Lorentzian objective function

The 3D structure of the genome plays a vital role in biological processes such as gene interaction, gene regulation, DNA replication and genome methylation. Advanced chromosomal conformation capture techniques, such as Hi-C and tethered conformation capture, can generate chromosomal contact data that can be used to computationally reconstruct 3D structures of the genome. We developed a novel restraint-based method that is capable of reconstructing 3D genome structures utilizing both intra-and inter-chromosomal contact data. Our method was robust to noise and performed well in comparison with a panel of existing methods on a controlled simulated data set. On a real Hi-C data set of the human genome, our method produced chromosome and genome structures that are consistent with 3D FISH data and known knowledge about the human chromosome and genome, such as, chromosome territories and the cluster of small chromosomes in the nucleus center with the exception of the chromosome 18. The tool and experimental data are available at https://missouri.box.com/v/LorDG.

Effects of short indels on protein structure and function in human genomes

Insertions and deletions (indels) represent the second most common type of genetic variations in human genomes. Indels can be deleterious and contribute to disease susceptibility as recent genome sequencing projects revealed a large number of indels in various cancer types. In this study, we investigated the possible effects of small coding indels on protein structure and function, and the baseline characteristics of indels in 2504 individuals of 26 populations from the 1000 Genomes Project. We found that each population has a distinct pattern in genes with small indels. Frameshift (FS) indels are enriched in olfactory receptor activity while non-frameshift (NFS) indels are enriched in transcription-related proteins. Structural analysis of NFS indels revealed that they predominantly adopt coil or disordered conformations, especially in proteins with transcription-related NFS indels. These results suggest that the annotated coding indels from the 1000 Genomes Project, while contributing to genetic variations and phenotypic diversity, generally do not affect the core protein structures and have no deleterious effect on essential biological processes. In addition, we found that a number of reference genome annotations might need to be updated due to the high prevalence of annotated homozygous indels in the general population.

cGAS surveillance of micronuclei links genome instability to innate immunity

DNA is strictly compartmentalized within the nucleus to prevent autoimmunity; despite this, cyclic GMP-AMP synthase (cGAS), a cytosolic sensor of double-stranded DNA, is activated in autoinflammatory disorders and by DNA damage. Precisely how cellular DNA gains access to the cytoplasm remains to be determined. Here, we report that cGAS localizes to micronuclei arising from genome instability in a mouse model of monogenic autoinflammation, after exogenous DNA damage and spontaneously in human cancer cells. Such micronuclei occur after mis-segregation of DNA during cell division and consist of chromatin surrounded by its own nuclear membrane. Breakdown of the micronuclear envelope, a process associated with chromothripsis, leads to rapid accumulation of cGAS, providing a mechanism by which self-DNA becomes exposed to the cytosol. cGAS is activated by chromatin, and consistent with a mitotic origin, micronuclei formation and the proinflammatory response following DNA damage are cell-cycle dependent. By combining live-cell laser microdissection with single cell transcriptomics, we establish that interferon-stimulated gene expression is induced in micronucleated cells. We therefore conclude that micronuclei represent an important source of immunostimulatory DNA. As micronuclei formed from lagging chromosomes also activate this pathway, recognition of micronuclei by cGAS may act as a cell-intrinsic immune surveillance mechanism that detects a range of neoplasia-inducing processes.

Nucleosome-free DNA regions differentially affect distant communication in chromatin

Communication between distantly spaced genomic regions is one of the key features of gene regulation in eukaryotes. Chromatin per se can stimulate efficient enhancer-promoter communication (EPC); however, the role of chromatin structure and dynamics in this process remains poorly understood. Here we show that nucleosome spacing and the presence of nucleosome-free DNA regions can modulate chromatin structure/dynamics and, in turn, affect the rate of EPC in vitro and in silico. Increasing the length of internucleosomal linker DNA from 25 to 60 bp results in more efficient EPC. The presence of longer nucleosome-free DNA regions can positively or negatively affect the rate of EPC, depending upon the length and location of the DNA region within the chromatin fiber. Thus the presence of histone-free DNA regions can differentially affect the efficiency of EPC, suggesting that gene regulation over a distance could be modulated by changes in the length of internucleosomal DNA spacers.

Towards quantitative analysis of gene regulation by enhancers

Enhancers are regulatory DNA sequences that can activate transcription over large distances. Recent studies have revealed the widespread role of distant activation in eukaryotic gene regulation and in the development of various human diseases, including cancer. Here we review recent progress in the field, focusing on new experimental and computational approaches that quantify the role of chromatin structure and dynamics during enhancer-promoter interactions in vitro and in vivo.

SAF-A Regulates Interphase Chromosome Structure through Oligomerization with Chromatin-Associated RNAs

Higher eukaryotic chromosomes are organized into topologically constrained functional domains; however, the molecular mechanisms required to sustain these complex interphase chromatin structures are unknown. A stable matrix underpinning nuclear organization was hypothesized, but the idea was abandoned as more dynamic models of chromatin behavior became prevalent. Here, we report that scaffold attachment factor A (SAF-A), originally identified as a structural nuclear protein, interacts with chromatin-associated RNAs (caRNAs) via its RGG domain to regulate human interphase chromatin structures in a transcription-dependent manner. Mechanistically, this is dependent on SAF-A's AAA⁺ ATPase domain, which mediates cycles of protein oligomerization with caRNAs, in response to ATP binding and hydrolysis. SAF-A oligomerization decompacts large-scale chromatin structure while SAF-A loss or monomerization promotes aberrant chromosome folding and accumulation of genome damage. Our results show that SAF-A and caRNAs form a dynamic, transcriptionally responsive chromatin mesh that organizes large-scale chromosome structures and protects the genome from instability.

Genomic properties of chromosomal bands are linked to evolutionary rearrangements and new centromere formation in primates

Chromosomal rearrangements in humans are largely related to pathological conditions, and phenotypic effects are also linked to alterations in the expression profile following nuclear relocation of genes between functionally different compartments, generally occupying the periphery or the inner part of the cell nuclei. On the other hand, during evolution, chromosomal rearrangements may occur apparently without damaging phenotypic effects and are visible in currently phylogenetically related species. To increase our insight into chromosomal reorganisation in the cell nucleus, we analysed 18 chromosomal regions endowed with different genomic properties in cell lines derived from eight primate species covering the entire evolutionary tree. We show that homologous loci, in spite of their evolutionary relocation along the chromosomes, generally remain localised to the same functional compartment of the cell nuclei. We conclude that evolutionarily successful chromosomal rearrangements are those that leave the nuclear position of the regions involved unchanged. On the contrary, in pathological situations, the effect typically observed is on gene structure alteration or gene nuclear reposition. Moreover, our data indicate that new centromere formation could potentially occur everywhere in the chromosomes, but only those emerging in very GC-poor/gene-poor regions, generally located in the nuclear periphery, have a high probability of being retained through evolution. This suggests that, in the cell nucleus of related species, evolutionary chromosomal reshufflings or new centromere formation does not alter the functionality of the regions involved or the interactions between different loci, thus preserving the expression pattern of orthologous genes.

Targeted Degradation of CTCF Decouples Local Insulation of Chromosome Domains from Genomic Compartmentalization

The molecular mechanisms underlying folding of mammalian chromosomes remain poorly understood. The transcription factor CTCF is a candidate regulator of chromosomal structure. Using the auxin-inducible degron system in mouse embryonic stem cells, we show that CTCF is absolutely and dose-dependently required for looping between CTCF target sites and insulation of topologically associating domains (TADs). Restoring CTCF reinstates proper architecture on altered chromosomes, indicating a powerful instructive function for CTCF in chromatin folding. CTCF remains essential for TAD organization in non-dividing cells. Surprisingly, active and inactive genome compartments remain properly segregated upon CTCF depletion, revealing that compartmentalization of mammalian chromosomes emerges independently of proper insulation of TADs. Furthermore, our data support that CTCF mediates transcriptional insulator function through enhancer blocking but not as a direct barrier to heterochromatin spreading. Beyond defining the functions of CTCF in chromosome folding, these results provide new fundamental insights into the rules governing mammalian genome organization.

How does chromatin package DNA within nucleus and regulate gene expression?

The human body is made up of 60 trillion cells, each cell containing 2 millions of genomic DNA in its nucleus. How is this genomic deoxyribonucleic acid [DNA] organised into nuclei? Around 1880, W. Flemming discovered a nuclear substance that was clearly visible on staining under primitive light microscopes and named it 'chromatin'; this is now thought to be the basic unit of genomic DNA organization. Since long before DNA was known to carry genetic information, chromatin has fascinated biologists. DNA has a negatively charged phosphate backbone that produces electrostatic repulsion between adjacent DNA regions, making it difficult for DNA to fold upon itself. In this article, we will try to shed light on how does chromatin package DNA within nucleus and regulate gene expression?

Anatomy of Mammalian Replication Domains

Genetic information is faithfully copied by DNA replication through many rounds of cell division. In mammals, DNA is replicated in Mb-sized chromosomal units called “replication domains.” While genome-wide maps in multiple cell types and disease states have uncovered both dynamic and static properties of replication domains, we are still in the process of understanding the mechanisms that give rise to these properties. A better understanding of the molecular basis of replication domain regulation will bring new insights into chromosome structure and function.

Transcription-dependent radial distribution of TCF7L2 regulated genes in chromosome territories

Human chromosomes occupy distinct territories in the interphase nucleus. Such chromosome territories (CTs) are positioned according to gene density. Gene-rich CTs are generally located in the center of the nucleus, while gene-poor CTs are positioned more towards the nuclear periphery. However, the association between gene expression levels and the radial positioning of genes within the CT is still under debate. In the present study, we performed three-dimensional fluorescence in situ hybridization experiments in the colorectal cancer cell lines DLD-1 and LoVo using whole chromosome painting probes for chromosomes 8 and 11 and BAC clones targeting four genes with different expression levels assessed by gene expression arrays and RT-PCR. Our results confirmed that the two over-expressed genes, MYC on chromosome 8 and CCND1 on chromosome 11, are located significantly further away from the center of the CT compared to under-expressed genes on the same chromosomes, i.e., DLC1 and SCN3B. When CCND1 expression was reduced after silencing the major transcription factor of the WNT/β-catenin signaling pathway, TCF7L2, the gene was repositioned and mostly detected in the interior of the CT. Thus, we suggest a non-random distribution in which over-expressed genes are located more towards the periphery of the respective CTs.

An Isochore Framework Underlies Chromatin Architecture

A recent investigation showed the existence of correlations between the architectural features of mammalian interphase chromosomes and the compositional properties of isochores. This result prompted us to compare maps of the Topologically Associating Domains (TADs) and of the Lamina Associated Domains (LADs) with the corresponding isochore maps of mouse and human chromosomes. This approach revealed that: 1) TADs and LADs correspond to isochores, i.e., isochores are the genomic units that underlie chromatin domains; 2) the conservation of TADs and LADs in mammalian genomes is explained by the evolutionary conservation of isochores; 3) chromatin domains corresponding to GC-poor isochores (e.g., LADs) show not only self-interactions but also intrachromosomal interactions with other domains also corresponding to GC-poor isochores even if located far away; in contrast, chromatin domains corresponding to GC-rich isochores (e.g., TADs) show more localized chromosomal interactions, many of which are inter-chromosomal. In conclusion, this investigation establishes a link between DNA sequences and chromatin architecture, explains the evolutionary conservation of TADs and LADs and provides new information on the spatial distribution of GC-poor/gene-poor and GC-rich/gene-rich chromosomal regions in the interphase nucleus.

Herpesvirus Latency: On the Importance of Positioning Oneself

The nucleus is composed of multiple compartments and domains, which directly or indirectly influence many cellular processes including gene expression, RNA splicing and maturation, protein post-translational modifications, and chromosome segregation. Nuclear-replicating viruses, especially herpesviruses, have co-evolved with the cell, adopting strategies to counteract and eventually hijack this hostile environment for their own benefit. This allows them to persist in the host for the entire life of an individual and to ensure their maintenance in the target species. Herpesviruses establish latency in dividing or postmitotic cells from which they can efficiently reactivate after sometimes years of a seemingly dormant state. Therefore, herpesviruses circumvent the threat of permanent silencing by reactivating their dormant genomes just enough to escape extinction, but not too much to avoid life-threatening damage to the host. In addition, herpesviruses that establish latency in dividing cells must adopt strategies to maintain their genomes in the daughter cells to avoid extinction by dilution of their genomes following multiple cell divisions. From a biochemical point of view, reactivation and maintenance of viral genomes in dividing cells occur successfully because the viral genomes interact with the nuclear architecture in a way that allows the genomes to be transmitted faithfully and to benefit from the nuclear micro-environments that allow reactivation following specific stimuli. Therefore, spatial positioning of the viral genomes within the nucleus is likely to be essential for the success of the latent infection and, beyond that, for the maintenance of herpesviruses in their respective hosts.

Variable chromatin structure revealed by in situ spatially correlated DNA cleavage mapping

Chromatin structure at the length scale encompassing local nucleosome-nucleosome interactions is thought to play a crucial role in regulating transcription and access to DNA. However, this secondary structure of chromatin remains poorly understood compared with the primary structure of single nucleosomes or the tertiary structure of long-range looping interactions. Here we report the first genome-wide map of chromatin conformation in human cells at the 1-3 nucleosome (50-500 bp) scale, obtained using ionizing radiation-induced spatially correlated cleavage of DNA with sequencing (RICC-seq) to identify DNA-DNA contacts that are spatially proximal. Unbiased analysis of RICC-seq signal reveals regional enrichment of DNA fragments characteristic of alternating rather than adjacent nucleosome interactions in tri-nucleosome units, particularly in H3K9me3-marked heterochromatin. We infer differences in the likelihood of nucleosome-nucleosome contacts among open chromatin, H3K27me3-marked, and H3K9me3-marked repressed chromatin regions. After calibrating RICC-seq signal to three-dimensional distances, we show that compact two-start helical fibre structures with stacked alternating nucleosomes are consistent with RICC-seq fragmentation patterns from H3K9me3-marked chromatin, while non-compact structures and solenoid structures are consistent with open chromatin. Our data support a model of chromatin architecture in intact interphase nuclei consistent with variable longitudinal compaction of two-start helical fibres.

Effects of DNA supercoiling on chromatin architecture

Disruptions in chromatin structure are necessary for the regulation of eukaryotic genomes, from remodelling of nucleosomes at the base pair level through to large-scale chromatin domains that are hundreds of kilobases in size. RNA polymerase is a powerful motor which, prevented from turning with the tight helical pitch of the DNA, generates over-wound DNA ahead of itself and under-wound DNA behind. Mounting evidence supports a central role for transcription-dependent DNA supercoiling in disrupting chromatin structure at all scales. This supercoiling changes the properties of the DNA helix in a manner that substantially alters the binding specificity of DNA binding proteins and complexes, including nucleosomes, polymerases, topoisomerases and transcription factors. For example, transient over-wound DNA destabilises nucleosome core particles ahead of a transcribing polymerase, whereas under-wound DNA facilitates pre-initiation complex formation, transcription factor binding and nucleosome core particle association behind the transcribing polymerase. Importantly, DNA supercoiling can also dissipate through DNA, even in a chromatinised context, to influence both local elements and large chromatin domains. We propose a model in which changes in unconstrained DNA supercoiling influences higher levels of chromatin organisation through the additive effects of DNA supercoiling on both DNA-protein and DNA-nucleosome interactions. This model links small-scale changes in DNA and chromatin to the higher-order fibre and large-scale chromatin structures, providing a mechanism relating gene regulation to chromatin architecture in vivo.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

Bidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.

Results

Here, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.

Conclusions

Our analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3125-0) contains supplementary material, which is available to authorized users.

Epigenomics and the structure of the living genome

Eukaryotic genomes are packaged into an extensively folded state known as chromatin. Analysis of the structure of eukaryotic chromosomes has been revolutionized by development of a suite of genome-wide measurement technologies, collectively termed “epigenomics.” We review major advances in epigenomic analysis of eukaryotic genomes, covering aspects of genome folding at scales ranging from whole chromosome folding down to nucleotide-resolution assays that provide structural insights into protein-DNA interactions. We then briefly outline several challenges remaining and highlight new developments such as single-cell epigenomic assays that will help provide us with a high-resolution structural understanding of eukaryotic genomes.

Effects of DNA supercoiling on chromatin architecture

Disruptions in chromatin structure are necessary for the regulation of eukaryotic genomes, from remodelling of nucleosomes at the base pair level through to large-scale chromatin domains that are hundreds of kilobases in size. RNA polymerase is a powerful motor which, prevented from turning with the tight helical pitch of the DNA, generates over-wound DNA ahead of itself and under-wound DNA behind. Mounting evidence supports a central role for transcription-dependent DNA supercoiling in disrupting chromatin structure at all scales. This supercoiling changes the properties of the DNA helix in a manner that substantially alters the binding specificity of DNA binding proteins and complexes, including nucleosomes, polymerases, topoisomerases and transcription factors. For example, transient over-wound DNA destabilises nucleosome core particles ahead of a transcribing polymerase, whereas under-wound DNA facilitates pre-initiation complex formation, transcription factor binding and nucleosome core particle association behind the transcribing polymerase. Importantly, DNA supercoiling can also dissipate through DNA, even in a chromatinised context, to influence both local elements and large chromatin domains. We propose a model in which changes in unconstrained DNA supercoiling influences higher levels of chromatin organisation through the additive effects of DNA supercoiling on both DNA-protein and DNA-nucleosome interactions. This model links small-scale changes in DNA and chromatin to the higher-order fibre and large-scale chromatin structures, providing a mechanism relating gene regulation to chromatin architecture in vivo.

Chromatin structure outside and inside the nucleus

The structure of the 30-nm chromatin fiber has provided, over the years, an important reference in chromatin studies. Originally derived from electron microscopic studies of soluble chromatin fibers released by restriction digestion, the gross structural features of such fragments have been supported by biophysical methods such as low angle X-ray and neutron scattering, sedimentation, light scattering, and electric dichroism. Electron microscopy and sedimentation velocity measurements demonstrated that reconstituted chromatin fibers, prepared from repeating arrays of high affinity nucleosome positioning sequences, retain the same overall features as observed for native chromatin fibers. It had been suggested that the 30 nm fiber might be the form assumed in vivo by transcriptionally silent chromatin, but individual gene or genome-wide studies of chromatin released from nuclei do not reveal any such simple correlation. Furthermore, even though the 30 nm fiber has been thought to represent an intermediate in the hierarchical folding of DNA into chromosomes, most analyses of chromatin folding within the nucleus do not detect any regular extended compact structures. However, there are important exceptions in chicken erythroid cell nuclei as well as in transcribed regions that form extended loops. Localized domains within the nucleus, either at the surface of chromosome domains or constrained as a specialized kind of constitutive heterochromatin by specific DNA binding proteins, may adopt 30 nm fiber-like structures.

Impact of the X Chromosome and sex on regulatory variation

The X Chromosome, with its unique mode of inheritance, contributes to differences between the sexes at a molecular level, including sex-specific gene expression and sex-specific impact of genetic variation. Improving our understanding of these differences offers to elucidate the molecular mechanisms underlying sex-specific traits and diseases. However, to date, most studies have either ignored the X Chromosome or had insufficient power to test for the sex-specific impact of genetic variation. By analyzing whole blood transcriptomes of 922 individuals, we have conducted the first large-scale, genome-wide analysis of the impact of both sex and genetic variation on patterns of gene expression, including comparison between the X Chromosome and autosomes. We identified a depletion of expression quantitative trait loci (eQTL) on the X Chromosome, especially among genes under high selective constraint. In contrast, we discovered an enrichment of sex-specific regulatory variants on the X Chromosome. To resolve the molecular mechanisms underlying such effects, we generated chromatin accessibility data through ATAC-sequencing to connect sex-specific chromatin accessibility to sex-specific patterns of expression and regulatory variation. As sex-specific regulatory variants discovered in our study can inform sex differences in heritable disease prevalence, we integrated our data with genome-wide association study data for multiple immune traits identifying several traits with significant sex biases in genetic susceptibilities. Together, our study provides genome-wide insight into how genetic variation, the X Chromosome, and sex shape human gene regulation and disease.

Nucleosomal arrays self-assemble into supramolecular globular structures lacking 30-nm fibers

The existence of a 30‐nm fiber as a basic folding unit for DNA packaging has remained a topic of active discussion. Here, we characterize the supramolecular structures formed by reversible Mg²⁺‐dependent self‐association of linear 12‐mer nucleosomal arrays using microscopy and physicochemical approaches. These reconstituted chromatin structures, which we call “oligomers”, are globular throughout all stages of cooperative assembly and range in size from ~50 nm to a maximum diameter of ~1,000 nm. The nucleosomal arrays were packaged within the oligomers as interdigitated 10‐nm fibers, rather than folded 30‐nm structures. Linker DNA was freely accessible to micrococcal nuclease, although the oligomers remained partially intact after linker DNA digestion. The organization of chromosomal fibers in human nuclei in situ was stabilized by 1 mM MgCl₂, but became disrupted in the absence of MgCl₂, conditions that also dissociated the oligomers in vitro. These results indicate that a 10‐nm array of nucleosomes has the intrinsic ability to self‐assemble into large chromatin globules stabilized by nucleosome–nucleosome interactions, and suggest that the oligomers are a good in vitro model for investigating the structure and organization of interphase chromosomes.

The "lnc" between 3D chromatin structure and X chromosome inactivation

The long non-coding RNA Xist directs a remarkable instance of developmentally regulated, epigenetic change known as X Chromosome Inactivation (XCI). By spreading in cis across the X chromosome from which it is expressed, Xist RNA facilitates the creation of a heritably silent, heterochromatic nuclear territory that displays a three-dimensional structure distinct from that of the active X chromosome. How Xist RNA attaches to and propagates across a chromosome and its influence over the three-dimensional (3D) structure of the inactive X are aspects of XCI that have remained largely unclear. Here, we discuss studies that have made significant contributions towards answering these open questions.

The folding landscape of the epigenome

The role of the spatial organization of chromatin in gene regulation is a long-standing but still open question. Experimentally it has been shown that the genome is segmented into epigenomic chromatin domains that are organized into hierarchical sub-nuclear spatial compartments. However, whether this non-random spatial organization only reflects or indeed contributes-and how-to the regulation of genome function remains to be elucidated. To address this question, we recently proposed a quantitative description of the folding properties of the fly genome as a function of its epigenomic landscape using a polymer model with epigenomic-driven attractions. We propose in this article, to characterize more deeply the physical properties of the 3D epigenome folding. Using an efficient lattice version of the original block copolymer model, we study the structural and dynamical properties of chromatin and show that the size of epigenomic domains and asymmetries in sizes and in interaction strengths play a critical role in the chromatin organization. Finally, we discuss the biological implications of our findings. In particular, our predictions are quantitatively compatible with experimental data and suggest a different mean of self-interaction in euchromatin versus heterochromatin domains.

A new oxygen modification cyclooctaoxygen binds to nucleic acids as sodium crown complex

BACKGROUND

Oxygen exists in two gaseous and six solid allotropic modifications. An additional allotropic modification of oxygen, the cyclooctaoxygen, was predicted to exist in 1990.

METHODS

Cyclooctaoxygen sodium was synthesized in vitro from atmospheric oxygen, or catalase effect-generated oxygen, under catalysis of cytosine nucleosides and either ninhydrin or eukaryotic low-molecular weight RNA. Thin-layer chromatographic mobility shift assays were applied on specific nucleic acids and the cyclooctaoxygen sodium complex.

RESULTS

We report the first synthesis and characterization of cyclooctaoxygen as its sodium crown complex, isolated in the form of three cytosine nucleoside hydrochloride complexes. The cationic cyclooctaoxygen sodium complex is shown to bind to nucleic acids (RNA and DNA), to associate with single-stranded DNA and spermine phosphate, and to be essentially non-toxic to cultured mammalian cells at 0.1-1.0mM concentration.

CONCLUSIONS

We postulate that cyclooctaoxygen is formed in most eukaryotic cells in vivo from dihydrogen peroxide in a catalase reaction catalyzed by cytidine and RNA. A molecular biological model is deduced for a first epigenetic shell of eukaryotic in vivo DNA. This model incorporates an epigenetic explanation for the interactions of the essential micronutrient selenium (as selenite) with eukaryotic in vivo DNA.

GENERAL SIGNIFICANCE

Since the sperminium phosphate/cyclooctaoxygen sodium complex is calculated to cover the active regions (2.6%) of bovine lymphocyte interphase genome, and 12.4% of murine enterocyte mitotic chromatin, we propose that the sperminium phosphate/cyclooctaoxygen sodium complex coverage of nucleic acids is essential to eukaryotic gene regulation and promoted proto-eukaryotic evolution.

Two Populations of Viral Minichromosomes Are Present in a Geminivirus-Infected Plant Showing Symptom Remission (Recovery)

Geminiviruses are important plant pathogens characterized by circular, single-stranded DNA (ssDNA) genomes. However, in the nuclei of infected cells, viral double-stranded DNA (dsDNA) associates with host histones to form a minichromosome. In phloem-limited geminiviruses, the characterization of viral minichromosomes is hindered by the low concentration of recovered complexes due to the small number of infected cells. Nevertheless, geminiviruses are both inducers and targets of the host posttranscriptional gene silencing (PTGS) and transcriptional gene silencing (TGS) machinery. We have previously characterized a "recovery" phenomenon observed in pepper plants infected with pepper golden mosaic virus (PepGMV) that is associated with a reduction of viral DNA and RNA levels, the presence of virus-related siRNAs, and an increase in the levels of viral DNA methylation. Initial micrococcal nuclease-based assays pinpointed the presence of different viral chromatin complexes in symptomatic and recovered tissues. Using the pepper-PepGMV system, we developed a methodology to obtain a viral minichromosome-enriched fraction that does not disturb the basic chromatin structural integrity, as evaluated by the detection of core histones. Using this procedure, we have further characterized two populations of viral minichromosomes in PepGMV-infected plants. After further purification using sucrose gradient sedimentation, we also observed that minichromosomes isolated from symptomatic tissue showed a relaxed conformation (based on their sedimentation rate), are associated with a chromatin activation marker (H3K4me3), and present a low level of DNA methylation. The minichromosome population obtained from recovered tissue, on the other hand, sedimented as a compact structure, is associated with a chromatin-repressive marker (H3K9me2), and presents a high level of DNA methylation.

IMPORTANCE

Viral minichromosomes have been reported in several animal and plant models. However, in the case of geminiviruses, there has been some recent discussion about the importance of this structure and the significance of the epigenetic modifications that it can undergo during the infective cycle. Major problems in this type of studies are the low concentration of these complexes in an infected plant and the asynchronicity of infected cells along the process; therefore, the complexes isolated in a given moment usually represent a mixture of cells at different infection stages. The recovery process observed in PepGMV-infected plants and the isolation procedure described here provide two distinct populations of minichromosomes that will allow a more precise characterization of the modifications of viral DNA and its host proteins associated along the infective cycle. This structure could be also an interesting model to study several processes involving plant chromatin.

Two Populations of Viral Minichromosomes Are Present in a Geminivirus-Infected Plant Showing Symptom Remission (Recovery)

Geminiviruses are important plant pathogens characterized by circular, single-stranded DNA (ssDNA) genomes. However, in the nuclei of infected cells, viral double-stranded DNA (dsDNA) associates with host histones to form a minichromosome. In phloem-limited geminiviruses, the characterization of viral minichromosomes is hindered by the low concentration of recovered complexes due to the small number of infected cells. Nevertheless, geminiviruses are both inducers and targets of the host posttranscriptional gene silencing (PTGS) and transcriptional gene silencing (TGS) machinery. We have previously characterized a "recovery" phenomenon observed in pepper plants infected with pepper golden mosaic virus (PepGMV) that is associated with a reduction of viral DNA and RNA levels, the presence of virus-related siRNAs, and an increase in the levels of viral DNA methylation. Initial micrococcal nuclease-based assays pinpointed the presence of different viral chromatin complexes in symptomatic and recovered tissues. Using the pepper-PepGMV system, we developed a methodology to obtain a viral minichromosome-enriched fraction that does not disturb the basic chromatin structural integrity, as evaluated by the detection of core histones. Using this procedure, we have further characterized two populations of viral minichromosomes in PepGMV-infected plants. After further purification using sucrose gradient sedimentation, we also observed that minichromosomes isolated from symptomatic tissue showed a relaxed conformation (based on their sedimentation rate), are associated with a chromatin activation marker (H3K4me3), and present a low level of DNA methylation. The minichromosome population obtained from recovered tissue, on the other hand, sedimented as a compact structure, is associated with a chromatin-repressive marker (H3K9me2), and presents a high level of DNA methylation.

IMPORTANCE

Viral minichromosomes have been reported in several animal and plant models. However, in the case of geminiviruses, there has been some recent discussion about the importance of this structure and the significance of the epigenetic modifications that it can undergo during the infective cycle. Major problems in this type of studies are the low concentration of these complexes in an infected plant and the asynchronicity of infected cells along the process; therefore, the complexes isolated in a given moment usually represent a mixture of cells at different infection stages. The recovery process observed in PepGMV-infected plants and the isolation procedure described here provide two distinct populations of minichromosomes that will allow a more precise characterization of the modifications of viral DNA and its host proteins associated along the infective cycle. This structure could be also an interesting model to study several processes involving plant chromatin.

Stable Chromosome Condensation Revealed by Chromosome Conformation Capture

Chemical cross-linking and DNA sequencing have revealed regions of intra-chromosomal interaction, referred to as topologically associating domains (TADs), interspersed with regions of little or no interaction, in interphase nuclei. We find that TADs and the regions between them correspond with the bands and interbands of polytene chromosomes of Drosophila. We further establish the conservation of TADs between polytene and diploid cells of Drosophila. From direct measurements on light micrographs of polytene chromosomes, we then deduce the states of chromatin folding in the diploid cell nucleus. Two states of folding, fully extended fibers containing regulatory regions and promoters, and fibers condensed up to 10-fold containing coding regions of active genes, constitute the euchromatin of the nuclear interior. Chromatin fibers condensed up to 30-fold, containing coding regions of inactive genes, represent the heterochromatin of the nuclear periphery. A convergence of molecular analysis with direct observation thus reveals the architecture of interphase chromosomes.

Polycomb-mediated chromatin compaction weathers the STORM

A recent super-resolution imaging study by Boettiger et al. elegantly demonstrates that three epigenetically defined, and functionally disparate, chromatin states have distinct folding characteristics in Drosophila nuclei.

Large-scale probabilistic 3D organization of human chromosome territories

There is growing evidence that chromosome territories (CT) have a probabilistic non-random arrangement within the cell nucleus of mammalian cells including radial positioning and preferred patterns of interchromosomal interactions that are cell-type specific. While it is generally assumed that the three-dimensional (3D) arrangement of genes within the CT is linked to genomic regulation, the degree of non-random organization of individual CT remains unclear. As a first step to elucidating the global 3D organization (topology) of individual CT, we performed multi-color fluorescence in situ hybridization using six probes extending across each chromosome in human WI38 lung fibroblasts. Six CT were selected ranging in size and gene density (1, 4, 12, 17, 18 and X). In-house computational geometric algorithms were applied to measure the 3D distances between every combination of probes and to elucidate data-mined structural patterns. Our findings demonstrate a high degree of non-random arrangement of individual CT that vary from chromosome to chromosome and display distinct changes during the cell cycle. Application of a classic, well-defined data mining and pattern recognition approach termed the 'k-means' generated 3D models for the best fit arrangement of each chromosome. These predicted models correlated well with the detailed distance measurements and analysis. We propose that the unique 3D topology of each CT and characteristic changes during the cell cycle provide the structural framework for the global gene expression programs of the individual chromosomes.

Genome Organization and Chromosome Architecture

How the same DNA sequences can function in the three-dimensional architecture of interphase nucleus, fold in the very compact structure of metaphase chromosomes, and go precisely back to the original interphase architecture in the following cell cycle remains an unresolved question to this day. The solution to this question presented here rests on the correlations that were found to hold between the isochore organization of the genome and the architecture of chromosomes from interphase to metaphase. The key points are the following: (1) The transition from the looped domains and subdomains of interphase chromatin to the 30-nm fiber loops of early prophase chromosomes goes through their unfolding into an extended chromatin structure (probably a 10-nm "beads-on-a-string" structure); (2) the architectural proteins of interphase chromatin, such as CTCF and cohesin subunits, are retained in mitosis and are part of the discontinuous protein scaffold of mitotic chromosomes; and (3) the conservation of the link between architectural proteins and their binding sites on DNA through the cell cycle explains the reversibility of the interphase to mitosis process and the "mitotic memory" of interphase architecture.

Nuclear positioning rather than contraction controls ordered rearrangements of immunoglobulin loci

Progenitor-B cells recombine their immunoglobulin (Ig) loci to create unique antigen receptors. Despite a common recombination machinery, the Ig heavy and Ig light chain loci rearrange in a stepwise manner. We studied pre-pro-B cells and Rag(-/-) progenitor-B cells to determine whether Ig locus contraction or nuclear positioning is decisive for stepwise rearrangements. We found that both Ig loci were contracted in pro-B and pre-B cells. Igh relocated from the nuclear lamina to central domains only at the pro-B cell stage, whereas, Igκ remained sequestered at the lamina, and only at the pre-B cell stage located to central nuclear domains. Finally, in vitro induced re-positioning of Ig alleles away from the nuclear periphery increased germline transcription of Ig loci in pre-pro-B cells. Thus, Ig locus contraction juxtaposes genomically distant elements to mediate efficient recombination, however, sequential positioning of Ig loci away from the nuclear periphery determines stage-specific accessibility of Ig loci.

Hierarchical block matrices as efficient representations of chromosome topologies and their application for 3C data integration

MOTIVATION

Recent advancements in molecular methods have made it possible to capture physical contacts between multiple chromatin fragments. The resulting association matrices provide a noisy estimate for average spatial proximity that can be used to gain insights into the genome organization inside the nucleus. However, extracting topological information from these data is challenging and their integration across resolutions is still poorly addressed. Recent findings suggest that a hierarchical approach could be advantageous for addressing these challenges.

RESULTS

We present an algorithmic framework, which is based on hierarchical block matrices (HBMs), for topological analysis and integration of chromosome conformation capture (3C) data. We first describe chromoHBM, an algorithm that compresses high-throughput 3C (HiT-3C) data into topological features that are efficiently summarized with an HBM representation. We suggest that instead of directly combining HiT-3C datasets across resolutions, which is a difficult task, we can integrate their HBM representations, and describe chromoHBM-3C, an algorithm which merges HBMs. Since three-dimensional (3D) reconstruction can also benefit from topological information, we further present chromoHBM-3D, an algorithm which exploits the HBM representation in order to gradually introduce topological constraints to the reconstruction process. We evaluate our approach in light of previous image microscopy findings and epigenetic data, and show that it can relate multiple spatial scales and provide a more complete view of the 3D genome architecture.

AVAILABILITY AND IMPLEMENTATION

The presented algorithms are available from: https://github.com/yolish/hbm

CONTACT

ys388@cam.ac.uk or pl219@cam.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.