Sequence structure of hidden 10.4-base repeat in the nucleosomes of C. elegans

By measuring prevailing distances between YY, YR, RR, and RY dinucleotides in the large database of the nucleosome DNA fragments from C. elegans, the consensus sequence structure of the nucleosome DNA repeat of C. elegans was reconstructed: (YYYYYRRRRR)n. An actual period was estimated to be 10.4 bases. The pattern is fully consistent with the nucleosome DNA patterns of other eukaryotes, as established earlier, and, thus, the YYYYYRRRRR repeat can be considered as consensus nucleosome DNA sequence repeat across eukaryotic species. Similar distance analysis for [A, T] dinucleotides suggested the related pattern (TTTYTARAAA)n where the TT and AA dinucleotides display rather out of phase behavior, contrary to the "AA or TT" in-phase periodicity, considered in some publications. A weak 5-base periodicity in the distribution of TA dinucleotides was detected.     

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.     

Universal full-length nucleosome mapping sequence probe

For the computational sequence-directed mapping of the nucleosomes, the knowledge of the nucleosome positioning motifs – 10–11 base long sequences – and respective matrices of bendability, is not sufficient, since there is no justified way to fuse these motifs in one continuous nucleosome DNA sequence. Discovery of the strong nucleosome (SN) DNA sequences, with visible sequence periodicity allows derivation of the full-length nucleosome DNA bendability pattern as matrix or consensus sequence. The SN sequences of three species (A. thaliana, C. elegans, and H. sapiens) are aligned (512 sequences for each species), and long (115 dinucleotides) matrices of bendability derived for the species. The matrices have strong common property – alternation of runs of purine–purine (RR) and pyrimidine–pyrimidine (YY) dinucleotides, with average period 10.4 bases. On this basis the universal [R,Y] consensus of the nucleosome DNA sequence is derived, with exactly defined positions of respective penta- and hexamers RRRRR, RRRRRR, YYYYY, and YYYYYY.     

Signals in chicken beta-globin DNA influence chromatin assembly in vitro.

We have confirmed the result that chicken beta-globin gene chromatin, which possesses the characteristics of active chromatin in erythroid cells, has shortened internucleosome spacings compared with bulk chromatin or that of the ovalbumin gene, which is inactive. To understand how the short (approximately 180-bp) nucleosome repeat arises specifically on beta-globin DNA, we have studied chromatin assembly of cloned chicken beta-globin DNA in a defined in vitro system. With chicken erythrocyte core histones and linker histone H5 as the only cellular components, a cloned 6.2-kb chicken beta-globin DNA fragment assembled into chromatin possessing a regular 180 +/- 5-bp repeat, very similar to what is observed in erythroid cells. A 2-kb DNA subfragment containing the beta A gene and promoter region, but lacking the downstream intergenic region between the beta A and epsilon genes, failed to generate a regular nucleosome array in vitro, suggesting that the intergenic region facilitates linker histone-induced nucleosome alignment. When the beta A gene was placed on a plasmid that contained a known chromatin-organizing signal, nucleosome alignment with a 180-bp periodicity was restored, whereas nucleosomes on flanking plasmid sequences possessed a 210-bp spacing periodicity. Our results suggest that the shortened 180-bp nucleosome spacing periodicity observed in erythroid cells is encoded in the beta-globin DNA sequence and that nucleosome alignment by linker histones is facilitated by sequences in the beta A-epsilon intergenic region.     

Periodicity in prokaryotic and eukaryotic genomes identified by power spectrum analysis

We used a power spectrum method to identify periodic patterns in nucleotide sequence, and characterized nucleotide sequences that confer periodicities to prokaryotic and eukaryotic genomes and genomes. A 10-bp periodicity was prevalent in hyperthermophilic bacteria and archaebacteria, and an 11-bp periodicity was prevalent in eubacteria. The 10-bp periodicity was also prevalent in the eukaryotes such as the worm Caenorhabditis elegans. Additionally, in the worm genome, a 68-bp periodicity in chromosome I, a 59-bp periodicity in chromosome II, and a 94-bp periodicity in chromosome III were found. In human chromosomes 21 and 22, approximately 167- or 84-bp periodicity was detected along the entire length of these chromosomes. Because the 167-bp is identical to the length of DNA that forms two complete helical turns in nucleosome organization, we speculated that the respective sequences may correspond to arrays of a special compact form of nucleosomes clustered in specific regions of the human chromosomes. This periodic element contained a high frequency of TGG. TGG-rich sequences are known to form a specific subset of folded DNA structures, and therefore, the sequences might have potential to form specific higher order structures related to the clustered occurrence of a specific form of the speculated nucleosomes.     

Surfeit locus gene homologs are widely distributed in invertebrate genomes.

The mouse Surfeit locus contains six sequence-unrelated genes (Surf-1 to -6) arranged in the tightest gene cluster so far described for mammals. The organization and juxtaposition of five of the Surfeit genes (Surf-1 to -5) are conserved between mammals and birds, and this may reflect a functional or regulatory requirement for the gene clustering. We have undertaken an evolutionary study to determine whether the Surfeit genes are conserved and clustered in invertebrate genomes. Drosophila melanogaster and Caenorhabditis elegans homologs of the mouse Surf-4 gene, which encodes an integral membrane protein associated with the endoplasmic reticulum, have been isolated. The amino acid sequences of the Drosophila and C. elegans homologs are highly conserved in comparison with the mouse Surf-4 protein. In particular, a dilysine motif implicated in endoplasmic reticulum localization of the mouse protein is conserved in the invertebrate homologs. We show that the Drosophila Surf-4 gene, which is transcribed from a TATA-less promoter, is not closely associated with other Drosophila Surfeit gene homologs but rather is located upstream from sequences encoding a homolog of a yeast seryl-tRNA synthetase protein. There are at least two closely linked Surf-3/rpL7a genes or highly polymorphic alleles of a single Surf-3/rpL7a gene in the C. elegans genome. The chromosomal locations of the C. elegans Surf-1, Surf-3/rpL7a, and Surf-4 genes have been determined. In D. melanogaster the Surf-3/rpL7a, Surf-4, and Surf-5 gene homologs and in C. elegans the Surf-1, Surf-3/rpL7a, Surf-4, and Surf-5 gene homologs are located on completely different chromosomes, suggesting that any requirement for the tight clustering of the genes in the Surfeit locus is restricted to vertebrate lineages.     

Correlation between dinucleotide periodicities and nucleosome positioning on mouse satellite DNA

Previous studies demonstrated 16 well-defined nucleosome locations (A-P) on a tandemly repeated prototype 234 base pair (bp) mouse satellite repeat unit. We have aligned the A-P fragments to search for DNA sequence elements that might contribute to nucleosome placement at these positions. Our results demonstrate a strikingly regular, uninterrupted, periodic pattern for the AA dinucleotide occurrences along the entire length of the aligned fragments. The periodicity of the AA occurrences is about 9.7 bp. The pattern exhibits a local minimum at position 74, near the nucleosome dyad axis of symmetry. Other dinucleotides--including AC: GT, CA: TG, and CC: GG--are also placed periodically, but their patterns of occurrence are less regular and less frequent than AA. The calculated spacings between consecutive preferred nucleosome locations on mouse satellite DNA are nearly identical, corresponding to multiples of 9.7 bp. The correlation between the periodicity of dinucleotide occurrences and the average spacing of nucleosome positions suggests that the preferred nucleosome locations recur at intervals that may correspond to the DNA helical repeat in the mouse satellite nucleosomes, and that the histone octamers sample (or slip along) the duplex in steps of 9.7 bp during nucleosome formation on mouse satellite DNA.     

Nucleosome formation potential of eukaryotic DNA: calculation and promoters analysis.

MOTIVATION: A rapid growth in the number of genes with known sequences calls for developing automated tools for their classification and analysis. It became clear that nucleosome packaging of eukaryotic DNA is very important for gene functioning. Automated computer tools for characterization of nucleosome packaging density could be useful for studying of gene regulation and genome annotation. RESULTS: A program for constructing nucleosome formation potential profiles of eukaryotic DNA sequences was developed. Nucleosome packaging density was analyzed for different functional types of human promoters. It was found that in promoters of tissue-specific genes, the nucleosome formation potential was essentially higher than in genes expressed in many tissues, or housekeeping genes. Hence, capability of nucleosome positioning in the promoter region may serve as a factor regulating gene expression. AVAILABILITY: The program for nucleosome sites recognition is included into the GeneExpress system; section 'DNA Nucleosomal Organization', http://wwwmgs.bionet.nsc.ru/mgs/programs/recon/.     

Modified curved DNA that could allow local DNA underwinding at the nucleosomal pseudodyad fails to position a nucleosome in vivo.

In competitive in vitro reconstitution experiments synthetic DNA composed of tandem repeats of the repetitive sequence (A/T)3NN(G/C)3NN, specifically the 20 bp 'TG sequence' (5'-TCGGTGTTAGAGCCTGTAAC-3'), was reported to associate with the histone octamer with an affinity higher than that of nucleosomally derived DNA. However, at least two groups have independently shown that tandem repeats of the TG sequence do not accommodate a stably positioned nucleosome in vivo. It was suggested that the anisotropic flexibility of the TG sequence, governed by a 10 bp sequence periodicity, is incompatible with the required underwinding of the DNA helix at the nucleosome pseudodyad while maintaining a bending preference that can be accommodated in the remainder of the nucleosome. Here we test this hypothesis directly by studying the in vivo nucleosomal structure of modified TG sequences designed to accommodate underwinding at the pseudodyad. We show that these modifications are not sufficient to allow stable incorporation of the TG sequence repeat into a nucleosome in vivo, but do note invasion from one end of the TG heptamer of a translationally random but rotationally constrained nucleosome. We discuss possible reasons for the absence of nucleosomes from the TG sequence in vivo.     

技术与方法体外组装含有组蛋白变体H2A.Z及H3.3的核小体

核小体是构成真核生物染色质的基本结构单位,组蛋白变体H2A.Z及H3.3对染色质结构及基因转录过程发挥着重要的调控作用。体内研究核小体及染色质结构受到诸多因素限制,体外重构含有H2A.Z及H3.3的核小体结构是研究与组蛋白变体相关基因表达调控的重要方法之一。实验表达纯化了6种组蛋白,在复性的过程中装配了含有H2A.Z和H3.3的组蛋白八聚体。基于DNA序列10bp周期性及序列模体设计了3条易于形成核小体的DNA序列,通过PCR大量扩增的方法,回收了标记Cy3荧光分子的目的DNA序列。采用盐透析法体外组装了含有H2A.Z和H3.3的核小体结构,利用荧光标记、EB染色及考马斯亮蓝染色检测了含有组蛋白变体的核小体形成效率及形成过程的吉布斯自由能变化。结果发现,设计的3条DNA序列可以有效地组装形成含有组蛋白电梯的核小体结构,而且随着组蛋白八聚体与DNA比例的增加,核小体的形成效率显著提高;采用Cy3荧光标记可以灵敏且定量地计算组装过程的吉布斯自由能。该方法的建立对研究组蛋白变体相关的结构生物学及转录调控等具有一定的意义。     

Improved alignment of nucleosome DNA sequences using a mixture model

DNA sequences that are present in nucleosomes have a preferential approximately 10 bp periodicity of certain dinucleotide signals, but the overall sequence similarity of the nucleosomal DNA is weak, and traditional multiple sequence alignment tools fail to yield meaningful alignments. We develop a mixture model that characterizes the known dinucleotide periodicity probabilistically to improve the alignment of nucleosomal DNAs. We assume that a periodic dinucleotide signal of any type emits according to a probability distribution around a series of 'hot spots' that are equally spaced along nucleosomal DNA with 10 bp period, but with a 1 bp phase shift across the middle of the nucleosome. We model the three statistically most significant dinucleotide signals, AA/TT, GC and TA, simultaneously, while allowing phase shifts between the signals. The alignment is obtained by maximizing the likelihood of both Watson and Crick strands simultaneously. The resulting alignment of 177 chicken nucleosomal DNA sequences revealed that all 10 distinct dinucleotides are periodic, however, with only two distinct phases and varying intensity. By Fourier analysis, we show that our new alignment has enhanced periodicity and sequence identity compared with center alignment. The significance of the nucleosomal DNA sequence alignment is evaluated by comparing it with that obtained using the same model on non-nucleosomal sequences.     

Artificial nucleosome positioning sequences tested in yeast minichromosomes: a strong rotational setting is not sufficient to position nucleosomes in vivo.

DNA sequences that support bending around the histone octamer ('rotational setting') are considered to be a major determinant of nucleosome positions. TG5 is an artificial positioning sequence containing 100 bp of an (A/T)3NN(G/C)3NN motif repeated with a 10 bp period. It provides a strong rotational setting and is superior to natural sequences in nucleosome formation in vitro [Shrader, T.E. and Crothers, D.M. (1989) Proc. Natl. Acad. Sci. USA, 86, 7418-7422]. To investigate the contribution of the rotational setting to nucleosome positioning in vivo, TG sequences were inserted in a nucleosome, at the edge of a nucleosome and in a nuclease sensitive region of yeast minichromosomes and the chromatin structures were analysed. In none of the constructs were TG sequences folded in a positioned nucleosome, demonstrating that the rotational setting played a subordinate role in the rough positioning in vivo. The rotational setting might fine tune the positions. Positioned nucleosomes were found overlapping the ends of TG, indicating that a discontinuity of the 10 bp periodicity of (A/T)3 and (G/C)3 near the centre of a nucleosome might be favourable for positioning and serve as a translational signal.     

Predicting nucleosome positioning in genomes: physical and bioinformatic approaches

In eukaryotic genomes, nucleosomes are responsible for packaging DNA and controlling gene expression. For this reason, an increasing interest is arising on computational methods capable of predicting the nucleosome positioning along genomes. In this review we describe and compare bioinformatic and physical approaches adopted to predict nucleosome occupancy along genomes. Computational analyses attempt at decoding the experimental nucleosome maps of genomes in terms of certain dinucleotide step periodicity observed along DNA. Such investigations show that highly significant information about the occurrence of a nucleosome along DNA is intrinsic in certain features of the sequence suggesting that DNA of eukaryotic genomes encodes nucleosome organization. Besides the bioinformatic approaches, physical models were proposed based on the sequence dependent conformational features of the DNA chain, which govern the free energy needed to transform recurrent DNA tracts along the genome into the nucleosomal shape.     

Functional elements and domains inferred from sequence comparisons of a heat shock gene in two nematodes

Caenorhabditis elegans and Caenorhabditis briggsae are two closely related nematode species that are nearly identical morphologically. Interspecific cross-hybridizing DNA appears to be restricted primarily to coding regions. We compared portions of the hsp-3 homologs, two grp 78-like genes, from C. elegans and C. briggsae and detected regions of DNA identity in the coding region, the 5' flanking DNAs, and the introns. The hsp-3 homologs share approximately 98% and 93% identity at the amino acid and nucleotide levels, respectively. Using the nucleotide substitution rate at the silent third position of the codons, we have estimated a lower limit for the date of divergence between C. elegans and C. briggsae to be approximately 23-32 million years ago. The 5' flanking DNAs and one of the introns contain elements that are highly conserved between C. elegans and C. briggsae. Some of the regions of nucleotide identity in the 5' flanking DNAs correspond to previously detected identities including viral enhancer sequences, a heat shock element, and an element present in the regulatory regions of mammalian grp78 and grp94 genes. We propose that a comparison of C. elegans and C. briggsae sequences will be useful in the detection of potential regulatory and structural elements.