A signal encoded in vertebrate DNA that influences nucleosome positioning and alignment

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.

DNaseI hypersensitivity at gene-poor, FSH dystrophy-linked 4q35.2

A subtelomeric region, 4q35.2, is implicated in facioscapulohumeral muscular dystrophy (FSHD), a dominant disease thought to involve local pathogenic changes in chromatin. FSHD patients have too few copies of a tandem 3.3-kb repeat (D4Z4) at 4q35.2. No phenotype is associated with having few copies of an almost identical repeat at 10q26.3. Standard expression analyses have not given definitive answers as to the genes involved. To investigate the pathogenic effects of short D4Z4 arrays on gene expression in the very gene-poor 4q35.2 and to find chromatin landmarks there for transcription control, unannotated genes and chromatin structure, we mapped DNaseI-hypersensitive (DH) sites in FSHD and control myoblasts. Using custom tiling arrays (DNase-chip), we found unexpectedly many DH sites in the two large gene deserts in this 4-Mb region. One site was seen preferentially in FSHD myoblasts. Several others were mapped >0.7 Mb from genes known to be active in the muscle lineage and were also observed in cultured fibroblasts, but not in lymphoid, myeloid or hepatic cells. Their selective occurrence in cells derived from mesoderm suggests functionality. Our findings indicate that the gene desert regions of 4q35.2 may have functional significance, possibly also to FSHD, despite their paucity of known genes.

Sequence information encoded in DNA that may influence long-range chromatin structure correlates with human chromosome functions

Little is known about the possible function of the bulk of the human genome. We have recently shown that long-range regular oscillation in the motif non-T, A/T, G (VWG) existing at ten-nucleotide multiples influences large-scale nucleosome array formation. In this work, we have determined the locations of all 100 kb regions that are predicted to form distinctive chromatin structures throughout each human chromosome (except Y). Using these data, we found that a significantly greater fraction of 300 kb sequences lacked annotated transcripts in genomic DNA regions > or = 300 kb that contained nearly continuous chromatin organizing signals than in control regions. We also found a relationship between the meiotic recombination frequency and the presence of strong VWG chromatin organizing signals. Large (> or = 300 kb) genomic DNA regions having low average recombination frequency are enriched in chromatin organizing signals. As additional controls, we show using chromosome 1 that the VWG motif signals are not enriched in randomly selected DNA regions having the mean size of the recombination coldspots, and that non-VWG motif sets do not generate signals that are enriched in recombination coldspots. We also show that tandemly repeated alpha satellite DNA contains strong VWG signals for the formation of distinctive nucleosome arrays, consistent with the low recombination activity of centromeres. Our correlations cannot be explained simply by variations in the GC content. Our findings suggest that a specific set of periodic DNA motifs encoded in genomic DNA, which provide signals for chromatin organization, influence human chromosome function.

Immunity related genes in dipterans share common enrichment of AT-rich motifs in their 5' regulatory regions that are potentially involved in nucleosome formation

Background

Understanding the transcriptional regulation mechanisms in response to environmental challenges is of fundamental importance in biology. Transcription factors associated to response elements and the chromatin structure had proven to play important roles in gene expression regulation. We have analyzed promoter regions of dipteran genes induced in response to immune challenge, in search for particular sequence patterns involved in their transcriptional regulation.

Results

5' upstream regions of D. melanogaster and A. gambiae immunity-induced genes and their corresponding orthologous genes in 11 non-melanogaster drosophilid species and Ae. aegypti share enrichment in AT-rich short motifs. AT-rich motifs are associated with nucleosome formation as predicted by two different algorithms. In A. gambiae and D. melanogaster, many immunity genes 5' upstream sequences also showed NFκB response elements, located within 500 bp from the transcription start site. In A. gambiae, the frequency of ATAA motif near the NFκB response elements was increased, suggesting a functional link between nucleosome formation/remodelling and NFκB regulation of transcription.

Conclusion

AT-rich motif enrichment in 5' upstream sequences in A. gambiae, Ae. aegypti and the Drosophila genus immunity genes suggests a particular pattern of nucleosome formation/chromatin organization. The co-occurrence of such motifs with the NFκB response elements suggests that these sequence signatures may be functionally involved in transcriptional activation during dipteran immune response. AT-rich motif enrichment in regulatory regions in this group of co-regulated genes could represent an evolutionary constrained signature in dipterans and perhaps other distantly species.

A high-resolution map of nucleosome positioning on a fission yeast centromere

A key element for defining the centromere identity is the incorporation of a specific histone H3, CENPA, known as Cnp1p in Schizosaccharomyces pombe. Previous studies have suggested that functional S. pombe centromeres lack regularly positioned nucleosomes and may involve chromatin remodeling as a key step of kinetochore assembly. We used tiling microarrays to show that nucleosomes are, in fact, positioned in regular intervals in the core of centromere 2, providing the first high-resolution map of regional centromere chromatin. Nucleosome locations are not disrupted by mutations in kinetochore protein genes cnp1, mis18, mis12, nuf2, mal2; overexpression of cnp1; or the deletion of ams2, which encodes a GATA-like factor participating in CENPA incorporation. Bioinformatics analysis of the centromere sequence indicates certain enriched motifs in linker regions between nucleosomes and reveals a sequence bias in nucleosome positioning. In addition, sequence analysis of nucleosome-free regions identifies novel binding sites of Ams2p. We conclude that centromeric nucleosome positions are stable and may be derived from the underlying DNA sequence.

Inferring natural selection on fine-scale chromatin organization in yeast

Despite its potential role in the evolution of complex phenotypes, the detection of negative (purifying) and positive selection on noncoding regulatory sequence has been elusive because of the inherent difficulty in predicting the functional consequences of mutations on noncoding sequence. Because the functioning of regulatory sequence depends upon both chromatin configuration and cis-regulatory factor binding, we investigate the idea that the functional conservation of regulatory regions should be associated with the conservation of sequence-dependent bending properties of DNA that determine its affinity for the nucleosome. Recent advances in the computational prediction of sequence-dependent affinity to nucleosomes provide an opportunity to distinguish between neutral and nonneutral evolution of fine-scale chromatin organization. Here, a statistical test is presented for detecting evolutionary conservation and/or adaptive evolution of nucleosome affinity from interspecies comparisons of DNA sequences. Local nucleosome affinities of homologous sequences were calculated using 2 recently published methods. A randomization test was applied to sites of mutation to evaluate the similarity of DNA-nucleosome affinity between several closely related species of Saccharomyces yeast. For most of the genes we analyzed, the conservation of local nucleosome affinity was detected at a few distinct locations in the upstream noncoding region. Our results also demonstrate that different patterns of chromatin evolution have shaped DNA-nucleosome interaction at the core promoters of TATA-containing and TATA-less genes and that elevated purifying selection has maintained low affinity for nucleosome in the core promoters of the latter group. Across the entire yeast genome, DNA-nucleosome interaction was also discovered to be significantly more conserved in TATA-less genes compared with TATA-containing genes.

Nucleosome positioning signals in genomic DNA

Although histones can form nucleosomes on virtually any genomic sequence, DNA sequences show considerable variability in their binding affinity. We have used DNA sequences of Saccharomyces cerevisiae whose nucleosome binding affinities have been experimentally determined (Yuan et al. 2005) to train a support vector machine to identify the nucleosome formation potential of any given sequence of DNA. The DNA sequences whose nucleosome formation potential are most accurately predicted are those that contain strong nucleosome forming or inhibiting signals and are found within nucleosome length stretches of genomic DNA with continuous nucleosome formation or inhibition signals. We have accurately predicted the experimentally determined nucleosome positions across a well-characterized promoter region of S. cerevisiae and identified strong periodicity within 199 center-aligned mononucleosomes studied recently (Segal et al. 2006) despite there being no periodicity information used to train the support vector machine. Our analysis suggests that only a subset of nucleosomes are likely to be positioned by intrinsic sequence signals. This observation is consistent with the available experimental data and is inconsistent with the proposal of a nucleosome positioning code. Finally, we show that intrinsic nucleosome positioning signals are both more inhibitory and more variable in promoter regions than in open reading frames in S. cerevisiae.

Analysis and development of the computer methods of nucleosome localization on DNA fragments with different AT-content

Trinucleotide parameter sets published previously were used for the development of the predictive method for the determining the nucleosome positions along the DNA. The choice of the type of parameter sets used depends upon AT-content of the fragment. Some limitations are imposed on these predictions due to the presence of A(n), T(n) tracts (in our case n>5 or =5) within the 145 bp fragment leading to the displacement or even the prohibition for the corresponding site to be occupied by nucleosomes. The predicted nucleosome positioning site with the large potential may influence on the choice of the proximal nucleosome positions with the weaker bending potentials as is revealed by the comparison with the micrococcal nuclease digestion map. Trinucleotide methods may be considered as advantageous in the comparison with the dinucleotide ones.

In silico approaches reveal the potential for DNA sequence-dependent histone octamer affinity to influence chromatin structure in vivo

Nucleosome positioning signals embedded within the DNA sequence have the potential to influence the detailed structure of the higher-order chromatin fibre. In two previous studies of long stretches of DNA, encompassing the chicken beta-globin and ovine beta-lactoglobulin genes, respectively, we mapped the relative affinity of every site for the core histone octamer. In both cases a periodic arrangement of the in vitro positioning sites suggests that they might influence the folding of a nucleosome chain into higher-order structure; this hypothesis was borne out in the case of the beta-lactoglobulin gene, where the distribution of the in vitro positioning sites is related to the positions nucleosomes actually occupy in sheep liver cells. Here, we have exploited the in vitro nucleosome positioning datasets to simulate nucleosomal organisation using in silico approaches. We use the high-resolution, quantitative positioning maps to define a one-dimensional positioning energy lattice, which can be populated with a defined number of nucleosomes. Monte Carlo techniques are employed to simulate the behaviour of the model at equilibrium to produce a set of configurations, which provide a probability-based occupancy map. Employing a variety of techniques we show that the occupancy maps are a sensitive function of the histone octamer density (nucleosome repeat length) and find that a minimal change in this property can produce dramatic localised changes in structure. Although simulations generally give rise to regular periodic nucleosomal arrangements, they often show octamer density-dependent discontinuities, which tend to co-localise with sequences that adopt distinctive chromatin structure in vivo. Furthermore, the overall organisation of simulated chromatin structures are more closely related to the situation in vivo than is the original in vitro positioning data, particularly at a nucleosome density corresponding to the in vivo state. Although our model is simplified, we argue that it provides a unique insight into the influence that DNA sequence can have in determining chromatin structure and could serve as a useful basis for the incorporation of other parameters.

Large-scale oscillation of structure-related DNA sequence features in human chromosome 21

Human chromosome 21 is the only chromosome in the human genome that exhibits oscillation of the (G+C) content of a cycle length of hundreds kilobases (kb) ( 500 kb near the right telomere). We aim at establishing the existence of a similar periodicity in structure-related sequence features in order to relate this (G+C)% oscillation to other biological phenomena. The following quantities are shown to oscillate with the same 500 kb periodicity in human chromosome 21: binding energy calculated by two sets of dinucleotide-based thermodynamic parameters, AA/TT and AAA/TTT bi- and tri-nucleotide density, 5'-TA-3' dinucleotide density, and signal for 10- or 11-base periodicity of AA/TT or AAA/TTT. These intrinsic quantities are related to structural features of the double helix of DNA molecules, such as base-pair binding, untwisting or unwinding, stiffness, and a putative tendency for nucleosome formation.

Aspects of large-scale chromatin structures in mouse liver nuclei can be predicted from the DNA sequence

The large amount of non-coding DNA present in mammalian genomes suggests that some of it may play a structural or functional role. We provide evidence that it is possible to predict computationally, from the DNA sequence, loci in mouse liver nuclei that possess distinctive nucleosome arrays. We tested the hypothesis that a 100 kb region of DNA possessing a strong, in-phase, dinucleosome period oscillation in the motif period-10 non-T, A/T, G, should generate a nucleosome array with a nucleosome repeat that is one-half of the dinucleosome oscillation period value, as computed by Fourier analysis of the sequence. Ten loci with short repeats, that would be readily distinguishable from the pervasive bulk repeat, were predicted computationally and then tested experimentally. We estimated experimentally that less than 20% of the chromatin in mouse liver nuclei has a nucleosome repeat length that is 15 bp, or more, shorter than the bulk repeat value of 195 +/- bp. All 10 computational predictions were confirmed experimentally with high statistical significance. Nucleosome repeats as short as 172 +/- 5 bp were observed for the first time in mouse liver chromatin. These findings may be useful for identifying distinctive chromatin structures computationally from the DNA sequence.

Long-range oscillation in a periodic DNA sequence motif may influence nucleosome array formation

We have experimentally examined the characteristics of nucleosome array formation in different regions of mouse liver chromatin, and have computationally analyzed the corresponding genomic DNA sequences. We have shown that the mouse adenosine deaminase (MADA) gene locus is packaged into an exceptionally regular nucleosome array with a shortened repeat, consistent with our computational prediction based on the DNA sequence. A survey of the mouse genome indicates that <10% of 70 kb windows possess a nucleosome-ordering signal, consisting of regular long-range oscillations in the period-10 triplet motif non-T, A/T, G (VWG), which is as strong as the signal in the MADA locus. A strong signal in the center of this locus, confirmed by in vitro chromatin assembly experiments, appears to cooperate with weaker, in-phase signals throughout the locus. In contrast, the mouse odorant receptor (MOR) locus, which lacks locus-wide signals, was representative of ∼40% of the mouse genomic DNA surveyed. Within this locus, nucleosome arrays were similar to those of bulk chromatin. Genomic DNA sequences which were computationally similar to MADA or MOR resulted in MADA- or MOR-like nucleosome ladders experimentally. Overall, we provide evidence that computationally predictable information in the DNA sequence may affect nucleosome array formation in animal tissue.

Wavelet Analysis of DNA Bending Profiles reveals Structural Constraints on the Evolution of Genomic Sequences

Analyses of genomic DNA sequences have shown in previous works that base pairs are correlated at large distances with scale-invariant statistical properties. We show in the present study that these correlations between nucleotides (letters) result in fact from long-range correlations (LRC) between sequence-dependent DNA structural elements (words) involved in the packaging of DNA in chromatin. Using the wavelet transform technique, we perform a comparative analysis of the DNA text and of the corresponding bending profiles generated with curvature tables based on nucleosome positioning data. This exploration through the optics of the so-called `wavelet transform microscope' reveals a characteristic scale of 100-200 bp that separates two regimes of different LRC. We focus here on the existence of LRC in the small-scale regime (≲ 200 bp). Analysis of genomes in the three kingdoms reveals that this regime is specifically associated to the presence of nucleosomes. Indeed, small scale LRC are observed in eukaryotic genomes and to a less extent in archaeal genomes, in contrast with their absence in eubacterial genomes. Similarly, this regime is observed in eukaryotic but not in bacterial viral DNA genomes. There is one exception for genomes of Poxviruses, the only animal DNA viruses that do not replicate in the cell nucleus and do not present small scale LRC. Furthermore, no small scale LRC are detected in the genomes of all examined RNA viruses, with one exception in the case of retroviruses. Altogether, these results strongly suggest that small-scale LRC are a signature of the nucleosomal structure. Finally, we discuss possible interpretations of these small-scale LRC in terms of the mechanisms that govern the positioning, the stability and the dynamics of the nucleosomes along the DNA chain. This paper is maily devoted to a pedagogical presentation of the theoretical concepts and physical methods which are well suited to perform a statistical analysis of genomic sequences. We review the results obtained with the so-called wavelet-based multifractal analysis when investigating the DNA sequences of various organisms in the three kingdoms. Some of these results have been announced in B. Audit et al. [1, 2].

Circle ligation of in vitro assembled chromatin indicates a highly flexible structure

Evidence is provided that some condensed linker histone-containing chromatin structures are highly flexible in solutions containing 2 mM Mg2+. Chromatin assembled in vitro +/- histone H5 on a 6.3 kb linear DNA fragment in 90 mM NaCl using the polyglutamic acid method sedimented fairly homogeneously. The H5-containing sample had s(20, w) values that were 58-69% greater than the sample lacking H5. Chromatin assembled on linear pUC19 plasmid DNA was treated with T4 DNA ligase in solutions containing 2 mM Mg2+ over a range of DNA concentrations. It was found that the intramolecular DNA ends of the chromatin could be joined together more efficiently than the intramolecular ends of the naked DNA at the higher DNA concentrations. This result could not be attributed to the effective reduction in DNA length by nucleosome formation. The chromatin structures formed did not have naked DNA tails extending from the ends as assessed by exonuclease III digestion. Chromatin assembled on DNA shortened by up to 420 bp gave very similar results, suggesting that the structure was a flexible one, rather than a rigid one having DNA ends that were fortuitously juxtaposed.

What positions nucleosomes?--A model

Here we propose a new determinant for localization of nucleosomes along genomic DNA, in addition to sequence-dependent features. The new specific class of chromatin scaling signals involves curved DNA. According to the observed positional distribution of DNA curvature, the new synchronizing signal occurs once per four nucleosomes on average. This new factor in nucleosome positioning should substantially influence the efficiency of biological reactions through regulatory factors microscopically and the entire chromatin structure through the 30 nm fiber structure macroscopically. Allocation of the new type of signals is found to be fixed evolutionarily although they could be shifted in accordance with the hierarchy of functional genomic structures.

Long-range correlations between DNA bending sites: relation to the structure and dynamics of nucleosomes

It has been established that the precise positioning of nucleosomes on genomic DNA can be achieved, at least for a minority of them, through sequence-dependent processes. However, to what extent DNA sequences play a role in the positioning of the major part of nucleosomes is still debated. The aim of the present study is to examine to what extent long-range correlations (LRC) are related to the presence of nucleosomes. Using the wavelet transform technique, we perform a comparative analysis of the DNA text and of the corresponding bending profiles generated with curvature tables based on nucleosome positioning data. The exploration of a number of eukaryotic and bacterial genomes through the optics of the so-called "wavelet transform microscope" reveals a characteristic scale of 100-200 bp that separates two regimes of different LRC. Here, we focus on the existence of LRC in the small-scale regime (10-200 bp) which are actually observed in eukaryotic genomes, in contrast to their absence in eubacterial genomes. Analysis of viral DNA genomes shows that, like their host's genomes, eukaryotic viruses present LRC but eubacterial viruses do not. There is one exception for genomes of poxviruses (Vaccinia and Melamoplus sanguinipes) which do not replicate in the cell nucleus and do not exhibit LRC. No small-scale LRC are detected in the genomes of all examined RNA viruses, with the exception of retroviruses. These results together with the observation of LRC between particular sequence motifs known to participate in the formation of nucleosomes (e.g. AA dinucleotides) strongly suggest that the 10-200 bp LRC are a signature of the sequence-dependence of nucleosome positioning. Finally, we discuss possible interpretations of these LRC in terms of the physical mechanisms that might govern the positioning and the dynamics of the nucleosomes along the DNA chain through cooperative processes.

Sequence patterns defining the 5' boundary of human genes

In studies of both short and relatively long human genomic DNA, we found a clustering of the consensus site for the transcription factor GCF at the 5' boundary of a subset of human genes. In studies of promoter regions with known transcription initiation site, the cluster of consensus GCF site appeared near the transcription initiation site and in some sequences it extended into the transcribed region defining the leader mRNA. We also found a detectable correlation between the 5' boundary of human genes and recognition motifs for other transcription factors that bind to GC-rich sequences. But in these cases, the correlation was not as general as the correlation observed for the consensus GCF site.

Radiolysis of nucleosome core DNA: a modelling approach

PURPOSE

To calculate the expected pattern of frank strand breaks (FSB) induced in the DNA of a nucleosome core particle (NCP) by the attack of radiolytic OH* radicals and to compare this with the experimental pattern of FSB induced by the in vitro irradiation of chicken NCP.

MATERIALS AND METHODS

The structure of NCP was obtained from the PDB databank. Using molecular modelling, the structure of a linear DNA fragment with the central 60 bp sequence of NCP-DNA was determined. The accessibility of the sugar attack sites to OH* radicals at each nucleotide of the linear fragment or of the NCP-DNA was calculated. The probability of reaction of OH* with each sugar reactive site was calculated using a Monte-Carlo method-based stochastic model.

RESULTS

The accessibility of attack sites in the NCP-DNA and the calculated probabilities of sugar attack are mainly modulated through masking by histones, and only slightly through bending. The periodicity of the calculated FSB probabilities is identical to that of the experimental FSB probabilities in chicken NCP (period 10.4 +/- 0.1 bp).

CONCLUSIONS

The good agreement between the calculated and experimental results recommends the presented modelling procedure as a tool in predicting the radiosensitivity of DNA in DNA-ligand complexes of known structure.

Syndromes associated with Homo sapiens pol II regulatory genes

The molecular basis of human characteristics is an intriguing but an unresolved problem. Human characteristics cover a broad spectrum, from the obvious to the abstract. Obvious characteristics may include morphological features such as height, shape, and facial form. Abstract characteristics may be hidden in processes that are controlled by hormones and the human brain. In this review we examine exaggerated characteristics presented as syndromes. Specifically, we focus on human genes that encode transcription factors to examine morphological, immunological, and hormonal anomalies that result from deletion, insertion, or mutation of genes that regulate transcription by RNA polymerase II (the Pol II genes). A close analysis of abnormal phenotypes can give clues into how sequence variations in regulatory genes and changes in transcriptional control may give rise to characteristics defined as complex traits.

Sequence-dependent DNA structure: dinucleotide conformational maps

We have used a computational model to calculate the potential energy surface for dinucleotide steps in double helical DNA as a function of the two principal degrees of freedom, slide and shift. By using a virtual bond to model the constraints imposed by the sugar-phosphate backbone, twist, roll, tilt and rise can be simultaneously optimised for any given values of slide and shift. Thus we have been able to construct complete conformational maps for all step types. For some steps, the maps agree well with experimental data from X-ray crystal structures, but other steps appear to be strongly perturbed by the effects of context (conformational coupling with the neighbouring steps). The optimised values of twist and roll show sequence-dependent variations consistent with the crystal structure data. The conformational maps allow us to construct adiabatic paths, and hence calculate the flexibility of each step with respect to slide and shift. Again the results agree well with the available experimental assignments of flexibility: YR steps, CA/TG and CG, are the most flexible and RR steps, such as AA, the least flexible.