30nm染色质纤维结构与表观遗传调控

真核生物的DNA以染色质形式通过逐级折叠压缩形成高级结构存在于细胞核中。染色质高级结构直接参与了真核基因的转录调控和其它与DNA相关的生物学事件,因此研究染色质高级结构对了解表观遗传学分子机制有着至关重要的作用。近些年,研究者们针对30 nm染色质高级结构提出了两个模型:螺线管模型和Zig-Zag模型。2014年,我们利用体外染色质组装体系重建了30 nm染色质纤维,运用高精度冷冻电镜技术得到了分辨率为11?的30 nm染色质纤维的精细结构,提出了30 nm染色质高级结构的左手双螺旋Zig-Zag模型。本文综述了30 nm染色质纤维结构研究方面的相关进展,并对30 nm染色质高级结构的表观遗传调控机理以及单分子成像和操纵技术在研究30 nm染色质高级结构中潜在的应用作出讨论和展望。     

HiCRes: a computational method to estimate and predict the genomic resolution of Hi-C libraries

Three-dimensional (3D) conformation of the chromatin is crucial to stringently regulate gene expression patterns and DNA replication in a cell-type specific manner. Hi-C is a key technique for measuring 3D chromatin interactions genome wide. Estimating and predicting the resolution of a library is an essential step in any Hi-C experimental design. Here, we present the mathematical concepts to estimate the resolution of a dataset and predict whether deeper sequencing would enhance the resolution. We have developed HiCRes, a docker pipeline, by applying these concepts to several Hi-C libraries.     

Uncovering the statistical physics of 3D chromosomal organization using data-driven modeling

In recent years, much effort has been devoted to understanding the three-dimensional (3D) organization of the genome and how genomic structure mediates nuclear function. The development of experimental techniques that combine DNA proximity ligation with high-throughput sequencing, such as Hi-C, have substantially improved our knowledge about chromatin organization. Numerous experimental advancements, not only utilizing DNA proximity ligation but also high-resolution genome imaging (DNA tracing), have required theoretical modeling to determine the structural ensembles consistent with such data. These 3D polymer models of the genome provide an understanding of the physical mechanisms governing genome architecture. Here, we present an overview of the recent advances in modeling the ensemble of 3D chromosomal structures by employing the maximum entropy approach combined with polymer physics. Particularly, we discuss the minimal chromatin model (MiChroM) along with the “maximum entropy genomic annotations from biomarkers associated with structural ensembles” (MEGABASE) model, which have been remarkably successful in the accurate modeling of chromosomes consistent with both Hi-C and DNA-tracing data.     

Structural organization of rat hepatocyte chromatin as visualized in thin frozen sections selectively stained for DNA

We studied the structure of rat hepatocyte chromatin in situ using thin frozen sections selectively stained for DNA after aldehyde fixation. Our results indicate that intranucleolar chromatin is arranged into three different organization levels, confirming the observations on Epon-embedded chromatin. These are: completely extended DNA filaments, with a thickness of approximately 3 nm, clustered in loose, roundish agglomerates, very long fibers with a thickness ranging from 15 to 35 nm and compact chromatin clumps. Both the fibers and the chromatin clumps frequently appeared to be composed of nucleosome-like particles. In the extranucleolar chromatin, agglomerates of extended DNA filaments and long fibers were never visualized. In contrast to data from Epon-embedded chromatin, we noticed that in frozen sections neither the nucleolar nor the extranucleolar compact chromatin appear to be organized into discrete, 20 to 30 nm superordered fibers.     

The diameters of frozen-hydrated chromatin fibers increase with DNA linker length: evidence in support of variable diameter models for chromatin

The diameters of chromatin fibers from Thyone briareus (sea cucumber) sperm (DNA linker length, n = 87 bp) and Necturus maculosus (mudpuppy) erythrocytes (n = 48 bp) were investigated. Soluble fibers were frozen into vitrified aqueous solutions of physiological ionic strength (124 mM), imaged by cryo-EM, and measured interactively using quantitative computer image-processing techniques. Frozen-hydrated Thyone and Necturus fibers had significantly different mean diameters of 43.5 nm (SD = 4.2 nm; SEM = 0.61 nm) and 32.0 nm (SD = 3.0 nm; SEM = 0.36 nm), respectively. Evaluation of previously published EM data shows that the diameters of chromatin from a large number of sources are proportional to linker length. In addition, the inherent variability in fiber diameter suggests a relationship between fiber structure and the heterogeneity of linker length. The cryo-EM data were in quantitative agreement with space-filling double-helical crossed-linker models of Thyone and Necturus chromatin. The data, however, do not support solenoid or twisted-ribbon models for chromatin that specify a constant 30 nm diameter. To reconcile the concept of solenoidal packing with the data, we propose a variable-diameter solid-solenoid model with a fiber diameter that increases with linker length. In principle, each of the variable diameter models for chromatin can be reconciled with local variations in linker length.     

Visualization of G1 chromosomes: a folded, twisted, supercoiled chromonema model of interphase chromatid structure

We have used light microscopy and serial thin-section electron microscopy to visualize intermediates of chromosome decondensation during G1 progression in synchronized CHO cells. In early G1, tightly coiled 100-130-nm "chromonema" fibers are visualized within partially decondensed chromatin masses. Progression from early to middle G1 is accompanied by a progressive uncoiling and straightening of these chromonema fibers. Further decondensation in later G1 and early S phase results in predominantly 60-80-nm chromonema fibers that can be traced up to 2-3 microns in length as discrete fibers. Abrupt transitions in diameter from 100-130 to 60-80 nm along individual fibers are suggestive of coiling of the 60-80-nm chromonema fibers to form the thicker 100-130-nm chromonema fiber. Local unfolding of these chromonema fibers, corresponding to DNA regions tens to hundreds of kilobases in length, reveal more loosely folded and extended 30-nm chromatin fibers. Kinks and supercoils appear as prominent features at all observed levels of folding. These results are inconsistent with prevailing models of chromosome structure and, instead, suggest a folded chromonema model of chromosome structure.     

The sequencing bias relaxed characteristics of Hi-C derived data and implications for chromatin 3D modeling

The 3D chromatin structure modeling by chromatin interactions derived from Hi-C experiments is significantly challenged by the intrinsic sequencing biases in these experiments. Conventional modeling methods only focus on the bias among different chromatin regions within the same experiment but neglect the bias arising from different experimental sequencing depth. We now show that the regional interaction bias is tightly coupled with the sequencing depth, and we further identify a chromatin structure parameter as the inherent characteristics of Hi-C derived data for chromatin regions. Then we present an approach for chromatin structure prediction capable of relaxing both kinds of sequencing biases by using this identified parameter. This method is validated by intra and inter cell-line comparisons among various chromatin regions for four human cell-lines (K562, GM12878, IMR90 and H1hESC), which shows that the openness of chromatin region is well correlated with chromatin function. This method has been executed by an automatic pipeline (AutoChrom3D) and thus can be conveniently used.     

Pulsed-field gel electrophoresis analysis of higher-order chromatin structures of Zea mays. Highly methylated DNA in the 50 kb chromatin structure

We have investigated the presence of higher-order chromatin structures in different maize tissues. Taking advantage of the pulsed-field gel electrophoresis technique to analyse large DNA fragments from intact nuclei and cells, we have determined the size distribution of the high-molecular-weight DNA fragments obtained from chromatin degradation by endogenous nucleases in isolated nuclei. Chromatin digestion leads to the appearance of stable DNA fragments of about 50 kb in all the tissues examined, suggesting the folding of DNA in higher-order chromatin domain structures. It has been reported that such chromatin domains are formed by loops of the 30 nm fibres anchored to the nuclear matrix by a complex set of proteins, including DNA topoisomerase II. Treatment of maize protoplasts with the calcium ionophore A23187 and the antitumour drug VM-26, which specifically inhibit the religation of the cleaved DNA in the topoisomerase II reaction, also produces the 50 kb structure. Analysis of the DNA contained in the 50 kb chromatin structure shows a higher degree of methylation than in bulk maize chromosomal DNA. The role of methylated DNA in the chromatin folding is discussed.     

Radial density distribution of chromatin: evidence that chromatin fibers have solid centers

Fiber diameter, radial distribution of density, and radius of gyration were determined from scanning transmission electron microscopy (STEM) of unstained, frozen-dried chromatin fibers. Chromatin fibers isolated under physiological conditions (ionic strength, 124 mM) from Thyone briareus sperm (DNA linker length, n = 87 bp) and Necturus maculosus erythrocytes (n = 48 bp) were analyzed by objective image-processing techniques. The mean outer diameters were determined to be 38.0 nm (SD = 3.7 nm; SEM = 0.36 nm) and 31.2 nm (SD = 3.6 nm; SEM = 0.32 nm) for Thyone and Necturus, respectively. These data are inconsistent with the twisted-ribbon and solenoid models, which predict constant diameters of approximately 30 nm, independent of DNA linker length. Calculated radial density distributions of chromatin exhibited relatively uniform density with no central hole, although the 4-nm hole in tobacco mosaic virus (TMV) from the same micrographs was visualized clearly. The existence of density at the center of chromatin fibers is in strong disagreement with the hollow-solenoid and hollow-twisted-ribbon models, which predict central holes of 16 and 9 nm for chromatin of 38 and 31 nm diameter, respectively. The cross-sectional radii of gyration were calculated from the radial density distributions and found to be 13.6 nm for Thyone and 11.1 nm for Necturus, in good agreement with x-ray and neutron scattering. The STEM data do not support the solenoid or twisted-ribbon models for chromatin fiber structure. They do, however, support the double-helical crossed-linker models, which exhibit a strong dependence of fiber diameter upon DNA linker length and have linker DNA at the center.     

Chromatin fibers are left-handed double helices with diameter and mass per unit length that depend on linker length

Four classes of models have been proposed for the internal structure of eukaryotic chromosome fibers--the solenoid, twisted-ribbon, crossed-linker, and superbead models. We have collected electron image and x-ray scattering data from nuclei, and isolated chromatin fibers of seven different tissues to distinguish between these models. The fiber diameters are related to the linker lengths by the equation: D(N) = 19.3 + 0.23 N, where D(N) is the external diameter (nm) and N is the linker length (base pairs). The number of nucleosomes per unit length of the fibers is also related to linker length. Detailed studies were done on the highly regular chromatin from erythrocytes of Necturus (mud puppy) and sperm of Thyone (sea cucumber). Necturus chromatin fibers (N = 48 bp) have diameters of 31 nm and have 7.5 +/- 1 nucleosomes per 10 nm along the axis. Thyone chromatin fibers (N = 87 bp) have diameters of 39 nm and have 12 +/- 2 nucleosomes per 10 nm along the axis. Fourier transforms of electron micrographs of Necturus fibers showed left-handed helical symmetry with a pitch of 25.8 +/- 0.8 nm and pitch angle of 32 +/- 3 degrees, consistent with a double helix. Comparable conclusions were drawn from the Thyone data. The data do not support the solenoid, twisted-ribbon, or supranucleosomal particle models. The data do support two crossed-linker models having left-handed double-helical symmetry and conserved nucleosome interactions.     

Microrheology for Hi-C Data Reveals the Spectrum of the Dynamic 3D Genome Organization

The one-dimensional information of genomic DNA is hierarchically packed inside the eukaryotic cell nucleus and organized in a three-dimensional (3D) space. Genome-wide chromosome conformation capture (Hi-C) methods have uncovered the 3D genome organization and revealed multiscale chromatin domains of compartments and topologically associating domains (TADs). Moreover, single-nucleosome live-cell imaging experiments have revealed the dynamic organization of chromatin domains caused by stochastic thermal fluctuations. However, the mechanism underlying the dynamic regulation of such hierarchical and structural chromatin units within the microscale thermal medium remains unclear. Microrheology is a way to measure dynamic viscoelastic properties coupling between thermal microenvironment and mechanical response. Here, we propose a new, to our knowledge, microrheology for Hi-C data to analyze the dynamic compliance property as a measure of rigidness and flexibility of genomic regions along with the time evolution. Our method allows the conversion of an Hi-C matrix into the spectrum of the dynamic rheological property along the genomic coordinate of a single chromosome. To demonstrate the power of the technique, we analyzed Hi-C data during the neural differentiation of mouse embryonic stem cells. We found that TAD boundaries behave as more rigid nodes than the intra-TAD regions. The spectrum clearly shows the dynamic viscoelasticity of chromatin domain formation at different timescales. Furthermore, we characterized the appearance of synchronous and liquid-like intercompartment interactions in differentiated cells. Together, our microrheology data derived from Hi-C data provide physical insights into the dynamics of the 3D genome organization.     

ASHIC: hierarchical Bayesian modeling of diploid chromatin contacts and structures

The recently developed Hi-C technique has been widely applied to map genome-wide chromatin interactions. However, current methods for analyzing diploid Hi-C data cannot fully distinguish between homologous chromosomes. Consequently, the existing diploid Hi-C analyses are based on sparse and inaccurate allele-specific contact matrices, which might lead to incorrect modeling of diploid genome architecture. Here we present ASHIC, a hierarchical Bayesian framework to model allele-specific chromatin organizations in diploid genomes. We developed two models under the Bayesian framework: the Poisson-multinomial (ASHIC-PM) model and the zero-inflated Poisson-multinomial (ASHIC-ZIPM) model. The proposed ASHIC methods impute allele-specific contact maps from diploid Hi-C data and simultaneously infer allelic 3D structures. Through simulation studies, we demonstrated that ASHIC methods outperformed existing approaches, especially under low coverage and low SNP density conditions. Additionally, in the analyses of diploid Hi-C datasets in mouse and human, our ASHIC-ZIPM method produced fine-resolution diploid chromatin maps and 3D structures and provided insights into the allelic chromatin organizations and functions. To summarize, our work provides a statistically rigorous framework for investigating fine-scale allele-specific chromatin conformations. The ASHIC software is publicly available at https://github.com/wmalab/ASHIC.     

Order and disorder in 30 nm chromatin fibers.

The tendency of DNA to form fibers upon condensation with counterions is reviewed. It is shown that chromatin fibers may acquire a relatively constant diameter of about 30 nm simply as an optimal size achieved upon neutralization of DNA, without requiring a repetitive internal structure. Thus the size of chromatin fibers would not be determined by any specific spatial interaction between DNA and histones. The driving force for the formation of fibers in chromatin would be similar to that found in proteins when they acquire a compact globular shape.