Nucleosome Distribution Research Articles

Compaction of Single-Molecule Megabase-Long Chromatin under the Influence of Macromolecular Crowding

The megabase-sized length of chromatin is highly relevant to the state of chromatin in vivo, where it is subject to a highly crowded environment and is organized in topologically associating domains of similar dimension. We developed an in vitro experimental chromatin model system reconstituted from T4 DNA (approximately 166 kbp) and histone octamers and studied the monomolecular compaction of this megabase-sized chromatin fiber under the influence of macromolecular crowding. We used single-molecule fluorescence microscopy and observed compaction in aqueous solutions containing poly(ethylene glycol) in the presence of monovalent (Na+ and K+) and divalent (Mg2+) cations. Both DNA and chromatin demonstrated compaction under comparable conditions in the presence of poly(ethylene glycol) and Na+ or Mg2+ salt. However, the mechanism of the compaction changed from a first-order phase transition for DNA to a continuous folding for megabase-sized chromatin fibers. A more efficient and pronounced chromatin compaction was observed in the presence of Na+ compared to K+. A flow-stretching technique to unfold DNA and chromatin coils was used to gain further insight into the morphology of partially folded chromatin fibers. The results revealed a distribution of partially folded chromatin fibers. This variability is likely the result of the heterogeneous distribution of nucleosomes on the DNA chain. The packaging of DNA in the form of chromatin in the crowded nuclear environment appears essential to ensure gradual conformational changes of DNA.

Chromatin Affinity Purification (ChAP) from Arabidopsis thaliana Rosette Leaves Using in vivo Biotinylation System.

Chromatin Affinity Purification (ChAP) is widely used to study chromatin architecture and protein complexes interacting with DNA. Here we present an efficient method for ChAP from Arabidopsis thaliana rosette leaves, in which in vivo biotinylation system is used. The chromatin is digested by Micrococcal Nuclease (MNase), hence the distribution of nucleosomes is also achieved. The in vivo biotinylation system was initially developed for Drosophila melanogaster ( Mito et al., 2005 ), but the presented protocol has been developed specifically for Arabidopsis thaliana ( Sura et al., 2017 ).

Drosophila Model for the Analysis of Genesis of LIM-kinase 1-Dependent Williams-Beuren Syndrome Cognitive Phenotypes: INDELs, Transposable Elements of the Tc1/Mariner Superfamily and MicroRNAs.

Genomic disorders, the syndromes with multiple manifestations, may occur sporadically due to unequal recombination in chromosomal regions with specific architecture. Therefore, each patient may carry an individual structural variant of DNA sequence (SV) with small insertions and deletions (INDELs) sometimes less than 10 bp. The transposable elements of the Tc1/mariner superfamily are often associated with hotspots for homologous recombination involved in human genetic disorders, such as Williams Beuren Syndromes (WBS) with LIM-kinase 1-dependent cognitive defects. The Drosophila melanogaster mutant agnts3 has unusual architecture of the agnostic locus harboring LIMK1: it is a hotspot of chromosome breaks, ectopic contacts, underreplication, and recombination. Here, we present the analysis of LIMK1-containing locus sequencing data in agnts3 and three D. melanogaster wild-type strains—Canton-S, Berlin, and Oregon-R. We found multiple strain-specific SVs, namely, single base changes and small INDEls. The specific feature of agnts3 is 28 bp A/T-rich insertion in intron 1 of LIMK1 and the insertion of mobile S-element from Tc1/mariner superfamily residing ~460 bp downstream LIMK1 3′UTR. Neither of SVs leads to amino acid substitutions in agnts3 LIMK1. However, they apparently affect the nucleosome distribution, non-canonical DNA structure formation and transcriptional factors binding. Interestingly, the overall expression of miRNAs including the biomarkers for human neurological diseases, is drastically reduced in agnts3 relative to the wild-type strains. Thus, LIMK1 DNA structure per se, as well as the pronounced changes in total miRNAs profile, probably lead to LIMK1 dysregulation and complex behavioral dysfunctions observed in agnts3 making this mutant a simple plausible Drosophila model for WBS.

A comparison of nucleosome organization in Drosophila cell lines.

Changes in the distribution of nucleosomes along the genome influence chromatin structure and impact gene expression by modulating the accessibility of DNA to transcriptional machinery. However, the role of genome-wide nucleosome positioning in gene expression and in maintaining differentiated cell states remains poorly understood. Drosophila melanogaster cell lines represent distinct tissue types and exhibit cell-type specific gene expression profiles. They thus could provide a useful tool for investigating cell-type specific nucleosome organization of an organism’s genome. To evaluate this possibility, we compared genome-wide nucleosome positioning and occupancy in five different Drosophila tissue-specific cell lines, and in reconstituted chromatin, and then tested for correlations between nucleosome positioning, transcription factor binding motifs, and gene expression. Nucleosomes in all cell lines were positioned in accordance with previously known DNA-nucleosome interactions, with helically repeating A/T di-nucleotide pairs arranged within nucleosomal DNAs and AT-rich pentamers generally excluded from nucleosomal DNA. Nucleosome organization in all cell lines differed markedly from in vitro reconstituted chromatin, with highly expressed genes showing strong nucleosome organization around transcriptional start sites. Importantly, comparative analysis identified genomic regions that exhibited cell line-specific nucleosome enrichment or depletion. Further analysis of these regions identified 91 out of 16,384 possible heptamer sequences that showed differential nucleosomal occupation between cell lines, and 49 of the heptamers matched one or more known transcription factor binding sites. These results demonstrate that there is differential nucleosome positioning between these Drosophila cell lines and therefore identify a system that could be used to investigate the functional significance of differential nucleosomal positioning in cell type specification.

Genome-scale identification of nucleosome organization by using 1000 porcine oocytes at different developmental stages.

The nucleosome is the basic structural unit of chromosomes, and its occupancy and distribution in promoters are crucial for the regulation of gene expression. During the growth process of porcine oocytes, the “growing” oocytes (SF) have a much higher transcriptional activity than the “fully grown” oocytes (BF). However, the chromosome status of the two kinds of oocytes remains poorly understood. In this study, we profiled the nucleosome distributions of SF and BF with as few as 1000 oocytes. By comparing the altered regions, we found that SF tended toward nucleosome loss and more open chromosome architecture than BF did. BF had decreased nucleosome occupancy in the coding region and increased nucleosome occupancy in the promoter compared to SF. The nucleosome occupancy of SF was higher than that of BF in the GC-poor regions, but lower than that of BF in the GC-rich regions. The nucleosome distribution around the transcriptional start site (TSS) of all the genes of the two samples was basically the same, but the nucleosome occupancy around the TSS of SF was lower than that of BF. GO functional annotation of genes with different nucleosome occupancy in promoter showed the genes were mainly involved in cell, cellular process, and metabolic process biological process. The results of this study revealed the dynamic reorganization of porcine oocytes in different developmental stages and the critical role of nucleosome arrangement during the oocyte growth process.

Genome-wide analysis of chromatin structures in Trypanosoma brucei using high-resolution MNase-ChIP-seq

Specific DNA-protein interactions are the basis for many important cellular mechanisms like the regulation of gene expression or replication. Knowledge about the precise genomic locations of DNA-protein interactions is important because it provides insight into the regulation of these processes. Recently, we have adapted an approach that combines micrococcal nuclease (MNase) digestion of chromatin with chromatin immunoprecipitation in Trypanosoma brucei. Here, we describe in detail how this method can be used to map the genome-wide distribution of nucleosomes or other DNA-binding proteins at high resolution in T. brucei.

Regulated large-scale nucleosome density patterns and precise nucleosome positioning correlate with V(D)J recombination

We show that the physical distribution of nucleosomes at antigen receptor loci is subject to regulated cell type-specific and lineage-specific positioning and correlates with the accessibility of these gene segments to recombination. At the Ig heavy chain locus (IgH), a nucleosome in pro-B cells is generally positioned over each IgH variable (VH) coding segment, directly adjacent to the recombination signal sequence (RSS), placing the RSS in a position accessible to the recombination activating gene (RAG) recombinase. These changes result in establishment of a specific chromatin organization at the RSS that facilitates accessibility of the genomic DNA for the RAG recombinase. In contrast, in mouse embryonic fibroblasts the coding segment is depleted of nucleosomes, which instead cover the RSS, thereby rendering it inaccessible. Pro-T cells exhibit a pattern intermediate between pro-B cells and mouse embryonic fibroblasts. We also find large-scale variations of nucleosome density over hundreds of kilobases, delineating chromosomal domains within IgH, in a cell type-dependent manner. These findings suggest that developmentally regulated changes in nucleosome location and occupancy, in addition to the known chromatin modifications, play a fundamental role in regulating V(D)J recombination. Nucleosome positioning-which has previously been observed to vary locally at individual enhancers and promoters-may be a more general mechanism by which cells can regulate the accessibility of the genome during development, at scales ranging from several hundred base pairs to many kilobases.

Epigenetics in male reproduction: effect of paternal diet on sperm quality and offspring health.

Epigenetic inheritance and its underlying molecular mechanisms are among the most intriguing areas of current biological and medical research. To date, studies have shown that both female and male germline development follow distinct paths of epigenetic events and both oocyte and sperm possess their own unique epigenomes. Fertilizing male and female germ cells deliver not only their haploid genomes but also their epigenomes, which contain the code for preimplantation and postimplantation reprogramming and embryonal development. For example, in spermatozoa, DNA methylation profile, DNA-associated proteins, protamine 1:protamine 2 ratio, nucleosome distribution pattern, histone modifications and other properties make up a unique epigenetic landscape. However, epigenetic factors and mechanisms possess certain plasticity and are affected by environmental conditions. Paternal and maternal lifestyle, including physical activity, nutrition and exposure to hazardous substances, can alter the epigenome and, moreover, can affect the health of their children. In male reproductive health, data are emerging on epigenetically mediated effects of a man's diet on sperm quality, for example through phytochemicals, minerals and vitamins, and nutritional support for subfertile men is already being used. In addition, studies in animal models and human epidemiological data point toward a transgenerational effect of the paternally contributed sperm epigenome on offspring health.

Enzymatic and chemical mapping of nucleosome distribution in purified micro- and macronuclei of the ciliated model organism, Tetrahymena thermophila.

Genomic distribution of the nucleosome, the basic unit of chromatin, contains important epigenetic information. To map nucleosome distribution in structurally and functionally differentiated micronucleus (MIC) and macronucleus (MAC) of the ciliate Tetrahymena thermophila, we have purified MIC and MAC and performed micrococcal nuclease (MNase) digestion as well as hydroxyl radical cleavage. Different factors that may affect MNase digestion were examined, to optimize mono-nucleosome production. Mono-nucleosome purity was further improved by ultracentrifugation in a sucrose gradient. As MNase concentration increased, nucleosomal DNA sizes in MIC and MAC converged on 147 bp, as expected for the nucleosome core particle. Both MNase digestion and hydroxyl radical cleavage consistently showed a nucleosome repeat length of ~200 bp in MAC of Tetrahymena, supporting ~50 bp of linker DNA. Our work has systematically tested methods currently available for mapping nucleosome distribution in Tetrahymena, and provided a solid foundation for future epigenetic studies in this ciliated model organism.

Dissecting relative contributions of cis- and trans-determinants to nucleosome distribution by comparing Tetrahymena macronuclear and micronuclear chromatin.

The ciliate protozoan Tetrahymena thermophila contains two types of structurally and functionally differentiated nuclei: the transcriptionally active somatic macronucleus (MAC) and the transcriptionally silent germ-line micronucleus (MIC). Here, we demonstrate that MAC features well-positioned nucleosomes downstream of transcription start sites and flanking splice sites. Transcription-associated trans-determinants promote nucleosome positioning in MAC. By contrast, nucleosomes in MIC are dramatically delocalized. Nucleosome occupancy in MAC and MIC are nonetheless highly correlated with each other, as well as with in vitro reconstitution and predictions based upon DNA sequence features, revealing unexpectedly strong contributions from cis-determinants. In particular, well-positioned nucleosomes are often matched with GC content oscillations. As many nucleosomes are coordinately accommodated by both cis- and trans-determinants, we propose that their distribution is shaped by the impact of these nucleosomes on the mutational and transcriptional landscape, and driven by evolutionary selection.

Changes in Gene Expression and Estrogen Receptor Cistrome in Mouse Liver Upon Acute E2 Treatment.

Transcriptional regulation by the estrogen receptor-α (ER) has been investigated mainly in breast cancer cell lines, but estrogens such as 17β-estradiol (E2) exert numerous extrareproductive effects, particularly in the liver, where E2 exhibits both protective metabolic and deleterious thrombotic actions. To analyze the direct and early transcriptional effects of estrogens in the liver, we determined the E2-sensitive transcriptome and ER cistrome in mice after acute administration of E2 or placebo. These analyses revealed the early induction of genes involved in lipid metabolism, which fits with the crucial role of ER in the prevention of liver steatosis. Characterization of the chromatin state of ER binding sites (BSs) in mice expressing or not ER demonstrated that ER is not required per se for the establishment and/or maintenance of chromatin modifications at the majority of its BSs. This is presumably a consequence of a strong overlap between ER and hepatocyte nuclear factor 4α BSs. In contrast, 40% of the BSs of the pioneer factor forkhead box protein a (Foxa2) were dependent upon ER expression, and ER expression also affected the distribution of nucleosomes harboring dimethylated lysine 4 of Histone H3 around Foxa2 BSs. We finally show that, in addition to a network of liver-specific transcription factors including CCAAT/enhancer-binding protein and hepatocyte nuclear factor 4α, ER might be required for proper Foxa2 function in this tissue.

A one-dimensional statistical mechanics model for nucleosome positioning on genomic DNA

The first level of folding of DNA in eukaryotes is provided by the so-called ‘10 nm chromatin fibre’, where DNA wraps around histone proteins (∼10 nm in size) to form nucleosomes, which go on to create a zig-zagging bead-on-a-string structure. In this work we present a one-dimensional statistical mechanics model to study nucleosome positioning within one such 10 nm fibre. We focus on the case of genomic sheep DNA, and we start from effective potentials valid at infinite dilution and determined from high-resolution in vitro salt dialysis experiments. We study positioning within a polynucleosome chain, and compare the results for genomic DNA to that obtained in the simplest case of homogeneous DNA, where the problem can be mapped to a Tonks gas []. First, we consider the simple, analytically solvable, case where nucleosomes are assumed to be point-like. Then, we perform numerical simulations to gauge the effect of their finite size on the nucleosomal distribution probabilities. Finally we compare nucleosome distributions and simulated nuclease digestion patterns for the two cases (homogeneous and sheep DNA), thereby providing testable predictions of the effect of sequence on experimentally observable quantities in experiments on polynucleosome chromatin fibres reconstituted in vitro.

Major Determinants of Nucleosome Organization

The stereotypical nucleosome organization at the gene promoters, a nucleosome-free region (NFR) flanked by regularly spaced nucleosomes [1-5], is common among different species, e.g. yeast, fly, mouse, human. We study the major determinants of nucleosome organization in vivo, and the common principles that can explain the similarities in nucleosome distribution among different organisms. We show that in yeast, the DNA sequence alone is not able to generate either the proper NFRs, or the nucleosome phasing. Transcription, chromatin remodeling factors, and transcription factors play key roles in establishing the distribution of the nucleosomes [6-9]. We analyze the roles of remodelers Isw1, Isw2, Chd1 and RSC in chromatin organization. We find that in vivo, Isw2 and RSC are responsible for creating proper NFRs and positioning the neighboring nucleosomes at the correct locations, while Isw1 and Chd1 act on gene bodies to generate regular arrays of nucleosomes.

Hierarchical regulation of the genome: global changes in nucleosome organization potentiate genome response.

Nucleosome occupancy is critically important in regulating access to the eukaryotic genome. Few studies in human cells have measured genome-wide nucleosome distributions at high temporal resolution during a response to a common stimulus. We measured nucleosome distributions at high temporal resolution following Kaposi's-sarcoma-associated herpesvirus (KSHV) reactivation using our newly developed mTSS-seq technology, which maps nucleosome distribution at the transcription start sites (TSS) of all human genes. Nucleosomes underwent widespread changes in organization 24 hours after KSHV reactivation and returned to their basal nucleosomal architecture 48 hours after KSHV reactivation. The widespread changes consisted of an indiscriminate remodeling event resulting in the loss of nucleosome rotational phasing signals. Additionally, one in six TSSs in the human genome possessed nucleosomes that are translationally remodeled. 72% of the loci with translationally remodeled nucleosomes have nucleosomes that moved to positions encoded by the underlying DNA sequence. Finally we demonstrated that these widespread alterations in nucleosomal architecture potentiated regulatory factor binding. These descriptions of nucleosomal architecture changes provide a new framework for understanding the role of chromatin in the genomic response, and have allowed us to propose a hierarchical model for chromatin-based regulation of genome response.

Comprehensive nucleosome mapping of the human genome in cancer progression.

Altered chromatin structure is a hallmark of cancer, and inappropriate regulation of chromatin structure may represent the origin of transformation. Important studies have mapped human nucleosome distributions genome wide, but the role of chromatin structure in cancer progression has not been addressed. We developed a MNase-Transcription Start Site Sequence Capture method (mTSS-seq) to map the nucleosome distribution at human transcription start sites genome-wide in primary human lung and colon adenocarcinoma tissue. Here, we confirm that nucleosome redistribution is an early, widespread event in lung (LAC) and colon (CRC) adenocarcinoma. These altered nucleosome architectures are consistent between LAC and CRC patient samples indicating that they may serve as important early adenocarcinoma markers. We demonstrate that the nucleosome alterations are driven by the underlying DNA sequence and potentiate transcription factor binding. We conclude that DNA-directed nucleosome redistributions are widespread early in cancer progression. We have proposed an entirely new hierarchical model for chromatin-mediated genome regulation.

Using deformation energy to analyze nucleosome positioning in genomes

By modulating the accessibility of genomic regions to regulatory proteins, nucleosome positioning plays important roles in cellular processes. Although intensive efforts have been made, the rules for determining nucleosome positioning are far from satisfaction yet. In this study, we developed a biophysical model to predict nucleosomal sequences based on the deformation energy of DNA sequences, and validated it against the experimentally determined nucleosome positions in the Saccharomyces cerevisiae genome, achieving very high success rates. Furthermore, using the deformation energy model, we analyzed the distribution of nucleosomes around the following three types of DNA functional sites: (1) double strand break (DSB), (2) single nucleotide polymorphism (SNP), and (3) origin of replication (ORI). We have found from the analyzed energy spectra that a remarkable “trough” or “valley” occurs around each of these functional sites, implying a depletion of nucleosome density, fully in accordance with experimental observations. These findings indicate that the deformation energy may play a key role for accurately predicting nucleosome positions, and that it can also provide a quantitative physical approach for in-depth understanding the mechanism of nucleosome positioning.

DPNuc: Identifying Nucleosome Positions Based on the Dirichlet Process Mixture Model.

Nucleosomes and the free linker DNA between them assemble the chromatin. Nucleosome positioning plays an important role in gene transcription regulation, DNA replication and repair, alternative splicing, and so on. With the rapid development of ChIP-seq, it is possible to computationally detect the positions of nucleosomes on chromosomes. However, existing methods cannot provide accurate and detailed information about the detected nucleosomes, especially for the nucleosomes with complex configurations where overlaps and noise exist. Meanwhile, they usually require some prior knowledge of nucleosomes as input, such as the size or the number of the unknown nucleosomes, which may significantly influence the detection results. In this paper, we propose a novel approach DPNuc for identifying nucleosome positions based on the Dirichlet process mixture model. In our method, Markov chain Monte Carlo (MCMC) simulations are employed to determine the mixture model with no need of prior knowledge about nucleosomes. Compared with three existing methods, our approach can provide more detailed information of the detected nucleosomes and can more reasonably reveal the real configurations of the chromosomes; especially, our approach performs better in the complex overlapping situations. By mapping the detected nucleosomes to a synthetic benchmark nucleosome map and two existing benchmark nucleosome maps, it is shown that our approach achieves a better performance in identifying nucleosome positions and gets a higher F-score. Finally, we show that our approach can more reliably detect the size distribution of nucleosomes.

Decoupling of divergent gene regulation by sequence-specific DNA binding factors

Divergent gene pairs (DGPs) are abundant in eukaryotic genomes. Since two genes in a DGP potentially share the same regulatory sequence, one might expect that they should be co-regulated. However, an inspection of yeast DGPs containing cell-cycle or stress response genes revealed that most DGPs are differentially-regulated. The mechanism underlying DGP differential regulation is not understood. Here, we showed that co- versus differential regulation cannot be explained by genetic features including promoter length, binding site orientation, TATA elements, nucleosome distribution, or presence of non-coding RNAs. Using time-lapse fluorescence microscopy, we carried out an in-depth study of a differentially regulated DGP, PFK26-MOB1. We found that their differential regulation is mainly achieved through two DNA-binding factors, Tbf1 and Mcm1. Similar to ‘enhancer-blocking insulators’ in higher eukaryotes, these factors shield the proximal promoter from the action of more distant transcription regulators. We confirmed the blockage function of Tbf1 using synthetic promoters. We further presented evidence that the blockage mechanism is widely used among genome-wide DGPs. Besides elucidating the DGP regulatory mechanism, our work revealed a novel class of insulators in yeast.

N6-Methyldeoxyadenosine Marks Active Transcription Start Sites in Chlamydomonas

N(6)-methyldeoxyadenosine (6mA or m(6)A) is a DNA modification preserved in prokaryotes to eukaryotes. It is widespread in bacteria and functions in DNA mismatch repair, chromosome segregation, and virulence regulation. In contrast, the distribution and function of 6mA in eukaryotes have been unclear. Here, we present a comprehensive analysis of the 6mA landscape in the genome of Chlamydomonas using new sequencing approaches. We identified the 6mA modification in 84% of genes in Chlamydomonas. We found that 6mA mainly locates at ApT dinucleotides around transcription start sites (TSS) with a bimodal distribution and appears to mark active genes. A periodic pattern of 6mA deposition was also observed at base resolution, which is associated with nucleosome distribution near the TSS, suggesting a possible role in nucleosome positioning. The new genome-wide mapping of 6mA and its unique distribution in the Chlamydomonas genome suggest potential regulatory roles of 6mA in gene expression in eukaryotic organisms.

DNA topoisomerase I affects polycomb group protein-mediated epigenetic regulation and plant development by altering nucleosome distribution in Arabidopsis.

It has been perplexing that DNA topoisomerases, enzymes that release DNA supercoils, play specific roles in development. In this study, using a floral stem cell model in Arabidopsis thaliana, we uncovered a role for TOPOISOMERASE1α (TOP1α) in Polycomb Group (PcG) protein-mediated histone 3 lysine 27 trimethylation (H3K27me3) at, and transcriptional repression of, the stem cell maintenance gene WUSCHEL (WUS). We demonstrated that H3K27me3 deposition at other PcG targets also requires TOP1α. Intriguingly, the repression of some, as well as the expression of many, PcG target genes requires TOP1α. The mechanism that unifies the opposing effects of TOP1α appears to lie in its role in decreasing nucleosome density, which probably allows the binding of factors that either recruit PcG, as we demonstrated for AGAMOUS at the WUS locus, or counteract PcG-mediated regulation. Although TOP1α reduces nucleosome density at all genes, the lack of a 5' nucleosome-free region is a feature that distinguishes PcG targets from nontargets and may condition the requirement for TOP1α for their expression. This study uncovers a connection between TOP1α and PcG, which explains the specific developmental functions of TOP1α.