Nucleosome Phasing in Tetrahymena Macronuclei1

Core-protected DNA can drive only 60% of the Tetrahymena thermophila macronuclear genome into duplexes in hybridization experiments. This core-protected DNA therefore contains only a subset of the genome complexity. We interpret this to mean that a large fraction, if not all, of the genome is phased with respect to nucleosome placement. Among the sequences present in total DNA and absent from core-protected DNA are most of the sequences containing N⁶-methyladenine (MeAde) residues, consistent with our previous demonstration that most of these residues lie in linker DNA. We show that these results are not due to artifacts resulting from the small size of the DNA driver, nor are they due to any sequence preferences exhibited by staphylococcal (staph) nuclease. This is the first evidence that nucleosome phasing may be a bulk genome characteristic.     

Enhancement of the nucleosomal pattern in sequences of lower complexity.

Intuitively, the complexity of a given DNA sequence is related to the number of various superimposed biological messages it contains. Here we assess the expectation that in nucleosome DNA sequences of lower linguistic complexity, the nucleosome DNA positioning pattern would be more pronounced than in those of higher linguistic complexity. The nucleosome DNA positioning pattern is one of the weakest (highly degenerate) sequence patterns. It has been extracted recently by specially designed multiple alignment procedures. We applied the most sensitive of these procedures to nearly equal subsets of a nucleosome database separated according to linguistic complexity. The pattern extracted from the subset of the simpler nucleosome sequences not only possesses all major attributes of the known nucleosomal pattern, but is substantially stronger with respect to amplitude in comparison with the total database. This result constitutes the first demonstration that a weak pattern can be significantly enhanced by selective treatment of a lower complexity subset of the sequence ensemble under consideration.     

Controls of Nucleosome Positioning in the Human Genome

Nucleosomes are important for gene regulation because their arrangement on the genome can control which proteins bind to DNA. Currently, few human nucleosomes are thought to be consistently positioned across cells; however, this has been difficult to assess due to the limited resolution of existing data. We performed paired-end sequencing of micrococcal nuclease-digested chromatin (MNase–seq) from seven lymphoblastoid cell lines and mapped over 3.6 billion MNase–seq fragments to the human genome to create the highest-resolution map of nucleosome occupancy to date in a human cell type. In contrast to previous results, we find that most nucleosomes have more consistent positioning than expected by chance and a substantial fraction (8.7%) of nucleosomes have moderate to strong positioning. In aggregate, nucleosome sequences have 10 bp periodic patterns in dinucleotide frequency and DNase I sensitivity; and, across cells, nucleosomes frequently have translational offsets that are multiples of 10 bp. We estimate that almost half of the genome contains regularly spaced arrays of nucleosomes, which are enriched in active chromatin domains. Single nucleotide polymorphisms that reduce DNase I sensitivity can disrupt the phasing of nucleosome arrays, which indicates that they often result from positioning against a barrier formed by other proteins. However, nucleosome arrays can also be created by DNA sequence alone. The most striking example is an array of over 400 nucleosomes on chromosome 12 that is created by tandem repetition of sequences with strong positioning properties. In summary, a large fraction of nucleosomes are consistently positioned—in some regions because they adopt favored sequence positions, and in other regions because they are forced into specific arrangements by chromatin remodeling or DNA binding proteins.     

Sequence motifs and free energies of selected natural and non-natural nucleosome positioning DNA sequences.

Our laboratories recently completed SELEX experiments to isolate DNA sequences that most-strongly favor or disfavor nucleosome formation and positioning, from the entire mouse genome or from even more diverse pools of chemically synthetic random sequence DNA. Here we directly compare these selected natural and non-natural sequences. We find that the strongest natural positioning sequences have affinities for histone binding and nucleosome formation that are sixfold or more lower than those possessed by many of the selected non-natural sequences. We conclude that even the highest-affinity sequence regions of eukaryotic genomes are not evolved for the highest affinity or nucleosome positioning power. Fourier transform calculations on the selected natural sequences reveal a special significance for nucleosome positioning of a motif consisting of approximately 10 bp periodic placement of TA dinucleotide steps. Contributions to histone binding and nucleosome formation from periodic TA steps are more significant than those from other periodic steps such as AA (=TT), CC (=GG) and more important than those from the other YR steps (CA (=TG) and CG), which are reported to have greater conformational flexibility in protein-DNA complexes even than TA. We report the development of improved procedures for measuring the free energies of even stronger positioning sequences that may be isolated in the future, and show that when the favorable free energy of histone-DNA interactions becomes sufficiently large, measurements based on the widely used exchange method become unreliable.     

Nucleosome positioning based on the sequence word composition

The DNA of all eukaryotic organisms is packaged into nucleosomes (a basic repeating unit of chromatin). A nucleosome consists of histone octamer wrapped by core DNA and linker histone H1 associated with linker DNA. It has profound effects on all DNA-dependent processes by affecting sequence accessibility. Understanding the factors that influence nucleosome positioning has great help to the study of genomic control mechanism. Among many determinants, the inherent DNA sequence has been suggested to have a dominant role in nucleosome positioning in vivo. Here, we used the method of minimum redundancy maximum relevance (mRMR) feature selection and the nearest neighbor algorithm (NNA) combined with the incremental feature selection (IFS) method to identify the most important sequence features that either favor or inhibit nucleosome positioning. We analyzed the words of 53,021 nucleosome DNA sequences and 50,299 linker DNA sequences of Saccharomyces cerevisiae. 32 important features were abstracted from 5,460 features, and the overall prediction accuracy through jackknife cross-validation test was 76.5%. Our results support that sequence-dependent DNA flexibility plays an important role in positioning nucleosome core particles and that genome sequence facilitates the rapid nucleosome reassembly instead of nucleosome depletion. Besides, our results suggest that there exist some additional features playing a considerable role in discriminating nucleosome forming and inhibiting sequences. These results confirmed that the underlying DNA sequence plays a major role in nucleosome positioning.     

Sequence Signatures of Nucleosome Positioning in Caenorhabditis elegans

Our recent investigation in the protist Trichomonas vaginalis suggested a DNA sequence periodicity with a unit length of 120.9 nt, which represents a sequence signature for nucleosome positioning. We now extended our observation in higher eukaryotes and identified a similar periodicity of 175 nt in length in Caenorhabditis elegans. In the process of defining the sequence compositional characteristics, we found that the 10.5-nt periodicity, the sequence signature of DNA double helix, may not be sufficient for cross-nucleosome positioning but provides essential guiding rails to facilitate positioning. We further dissected nucleosome-protected sequences and identified a strong positive purine (AG) gradient from the 5′-end to the 3′-end, and also learnt that the nucleosome-enriched regions are GC-rich as compared to the nucleosome-free sequences as purine content is positively correlated with GC content. Sequence characterization allowed us to develop a hidden Markov model (HMM) algorithm for decoding nucleosome positioning computationally, and based on a set of training data from the fifth chromosome of C. elegans, our algorithm predicted 60%-70% of the well-positioned nucleosomes, which is 15%-20% higher than random positioning. We concluded that nucleosomes are not randomly positioned on DNA sequences and yet bind to different genome regions with variable stability, well-positioned nucleosomes leave sequence signatures on DNA, and statistical positioning of nucleosomes across genome can be decoded computationally based on these sequence signatures.     

Effects of DNA sequence and histone-histone interactions on nucleosome placement

Using competitive reconstitution, we have refined the parameters for the binding of histone octamers to artificial nucleosome-positioning sequences of the form: (A/T3nn(G/C)3nn. We find that the optimal period between flexible segments is approximately 10.1 base-pairs, supporting the view that the DNA on the nucleosome surface is overwound. The strongest requirement for flexible DNA is near the protein dyad. However, we see no indication of changes in DNA helical repeat in this region. Using a series of repetitive sequences, we confirm that neither all A/T-rich nor all G/C-rich regions are identical in promoting nucleosome formation. Surprisingly, A/T-rich segments containing the TpA step, subject to purine-purine clash in the minor groove, favor nucleosome formation over sequences lacking this step. Short tracts of adenine residues are found to position on the histone surface like other A/T-rich regions, in the manner predicted by the direction of their sequence-directed bends as determined by electrophoretic methods. Tracts containing five adenine residues are extremely aniostropic in their flexibility and are strongly detrimental to nucleosome formation when positioned for major groove compression. Longer adenine tracts are found to position near the ends of the nucleosomal DNA. However, other positions may be occupied by an A12 tract, with only a minor penalty in the free energy of nucleosome formation. Overall, reconstituted nucleosome positions are translationally degenerate, suggesting a weak dependence on DNA flexibility for nucleosome positioning. Dinucleosomal reconstitutions on tandem dimers of the 5 S RNA gene of Lytechinus variegatus demonstrate a weak phasing dependence for the interaction between nucleosomes. This interaction is maximal for the 202 base-pair repeat and suggests a co-operative mechanism for the formation of ordered nucleosomal arrays based on a combination of DNA flexibility and nucleosome-nucleosome interactions.     

High‐resolution analysis of DNA synthesis start sites and nucleosome architecture at efficient mammalian replication origins

DNA replication origins are poorly characterized genomic regions that are essential to recruit and position the initiation complex to start DNA synthesis. Despite the lack of specific replicator sequences, initiation of replication does not occur at random sites in the mammalian genome. This has lead to the view that DNA accessibility could be a major determinant of mammalian origins. Here, we performed a high‐resolution analysis of nucleosome architecture and initiation sites along several origins of different genomic location and firing efficiencies. We found that mammalian origins are highly variable in nucleosome conformation and initiation patterns. Strikingly, initiation sites at efficient CpG island‐associated origins always occur at positions of high‐nucleosome occupancy. Origin recognition complex (ORC) binding sites, however, occur at adjacent but distinct positions marked by labile nucleosomes. We also found that initiation profiles mirror nucleosome architecture, both at endogenous origins and at a transgene in a heterologous system. Our studies provide a unique insight into the relationship between chromatin structure and initiation sites in the mammalian genome that has direct implications for how the replication programme can be accommodated to diverse epigenetic scenarios.     

Nucleosome structure and conformational changes.

We have used a variety of chemical probes to measure the accessibility of DNA on the surface of the nucleosome. We review these results, and describe new experiments which show that T4 phage DNA can form complexes with the core histones, possessing the properties of normal nucleosomes. Since T4 DNA is largely occupied by glucose residues in the major groove, this suggests (as did earlier probe experiments) that the major groove is not filled with histone amino acid side chains. We also report results of recent measurements which appear to show that only a few strong charge interactions are involved in the attachment of the terminal 20 nucleotide pairs at each end of nucleosome core DNA. We speculate on the possible functional significance of the accessibility of DNA revealed by all of these experiments. We have also examined conformational changes induced in nucleosomes at high ionic strength (0.5-0.7M NaCl). The frictional coefficient is found to undergo a small increase in this region, not consistent with models in which the nucleosome is completely unfolded, but possibly reflecting the dissociation of terminal DNA from the nucleosome surface.     

Nucleosomes affect local transformation efficiency

Genetic transformation is a natural process during which foreign DNA enters a cell and integrates into the genome. Apart from its relevance for horizontal gene transfer in nature, transformation is also the cornerstone of today''s recombinant gene technology. Despite its importance, relatively little is known about the factors that determine transformation efficiency. We hypothesize that differences in DNA accessibility associated with nucleosome positioning may affect local transformation efficiency. We investigated the landscape of transformation efficiency at various positions in the Saccharomyces cerevisiae genome and correlated these measurements with nucleosome positioning. We find that transformation efficiency shows a highly significant inverse correlation with relative nucleosome density. This correlation was lost when the nucleosome pattern, but not the underlying sequence was changed. Together, our results demonstrate a novel role for nucleosomes and also allow researchers to predict transformation efficiency of a target region and select spots in the genome that are likely to yield higher transformation efficiency.     

Telomeric nucleosomes are intrinsically mobile

Nucleosomes are no longer considered only static basic units that package eukaryotic DNA but they emerge as dynamic players in all chromosomal processes. Regulatory proteins can gain access to recognition sequences hidden by the histone octamer through the action of ATP-dependent chromatin remodeling complexes that cause nucleosome sliding. In addition, it is known that nucleosomes are able to spontaneously reposition along the DNA due to intrinsic dynamic properties, but it is not clear yet to what extent sequence-dependent dynamic properties contribute to nucleosome repositioning. Here, we study mobility of nucleosomes formed on telomeric sequences as a function of temperature and ionic strength. We find that telomeric nucleosomes are highly intrinsically mobile under physiological conditions, whereas nucleosomes formed on an average DNA sequence mostly remain in the initial position. This indicates that DNA sequence affects not only the thermodynamic stability and the positioning of nucleosomes but also their dynamic properties. Moreover, our findings suggest that the high mobility of telomeric nucleosomes may be relevant to the dynamics of telomeric chromatin.     

Molecular components of the archaeal nucleosome

Here we describe the organization of the archaeal nucleosome, in which four archaeal histones are circumscribed by approximately 80 bp of DNA. Through a combination of sequence comparisons, 3D structural studies, site-directed mutagenesis and assays for DNA binding, we have assigned functions to most of the individual residues in the histone fold of the representative archaeal histone rHMfB. By SELEX selection, the sequences of DNA molecules that are most readily bound and wrapped by rHMfB into archaeal nucleosomes in vitro have been identified, and these define DNA structures that position archaeal nucleosome assembly.     

DNA sequence alterations affect nucleosome array formation of the chicken ovalbumin gene

The role of the large amount (more than half of the genome) of noncoding DNA in higher organisms is not well understood. DNA evolved to function in the context of chromatin, and the possibility exists that some of the noncoding DNA serves to influence chromatin structure and function. In this age of genomics and bioinformatics, genomic DNA sequences are being searched for informational content beyond the known genetic code. The discovery that period-10 non-T, A/T, G (VWG) triplets are among the most abundant motifs in human genomic DNA suggests that they may serve some function in higher organisms. In this paper, we provide direct evidence that the regular oscillation of period-10 VWG that occurs in the chicken ovalbumin gene sequence with a dinucleosome-like period facilitates nucleosome array formation. Using a linker histone-dependent in vitro chromatin assembly system that spontaneously aligns nucleosomes into a physiological array, we show that nucleosomes tend to avoid DNA regions with low period-10 VWG counts. This avoidance leads to the formation of an array with a nucleosome repeat equal to half the period value of the oscillation in period-10 VWG, as determined by Fourier analysis. Two different half-period deletions in the wild-type DNA sequence altered the nucleosome array, as predicted computationally. In contrast, a full-period deletion had an insignificant effect on the nucleosome array formed, also consistent with the prediction. An inversion mutation, with no DNA sequences deleted, again altered the nucleosome array formed, as predicted computationally. Hence, a VWG dinucleosome signal is plausible.