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.     

Sequence Structure of Hidden 10.4-base Repeat in the Nucleosomes of C. elegans

Abstract

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.     

Unmethylated and methylated CpG dinucleotides distinctively regulate the physical properties of DNA

In eukaryotic cells, DNA has to bend significantly to pack inside the nucleus. Physical properties of DNA such as bending flexibility and curvature are expected to affect DNA packaging and partially determine the nucleosome positioning patterns inside a cell. DNA CpG methylation, the most common epigenetic modification found in DNA, is known to affect the physical properties of DNA. However, its detailed role in nucleosome formation is less well‐established. In this study, we evaluated the effect of defined CpG patterns (unmethylated and methylated) on DNA structure and their respective nucleosome‐forming ability. Our results suggest that the addition of CpG dinucleotides, either as a (CG)_n stretch or (CGX₈)_n repeats at 10 bp intervals, lead to reduced hydrodynamic radius and decreased nucleosome‐forming ability of DNA. This effect is more predominant for a DNA stretch ((CG)₅) located in the middle of a DNA fragment. Methylation of CpG sites, surprisingly, seems to reduce the difference in DNA structure and nucleosome‐forming ability among DNA constructs with different CpG patterns. Our results suggest that unmethylated and methylated CpG patterns can play very different roles in regulating the physical properties of DNA. CpG methylation seems to reduce the DNA conformational variations affiliated with defined CpG patterns. Our results can have significant bearings in understanding the nucleosome positioning pattern in living organisms modulated by DNA sequences and epigenetic features. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 517–524, 2014.     

The DNA Sequence-dependence of Nucleosome Positioning in vivo and in vitro

Abstract

The contribution of histone-DNA interactions to nucleosome positioning in vivo is currently a matter of debate. We argue here that certain nucleosome positions, often in promoter regions, in yeast may be, at least in part, specified by the DNA sequence. In contrast other positions may be poorly specified. Positioning thus has both statistical and DNA-determined components. We further argue that the relative affinity of the octamer for different DNA sequences can vary and therefore the interaction of histones with the DNA is a ‘tunable’ property.     

A determining influence for CpG dinucleotides on nucleosome positioning in vitro

DNA sequence information that directs the translational positioning of nucleosomes can be attenuated by cytosine methylation when a short run of CpG dinucleotides is located close to the dyad axis of the nucleosome. Here, we show that point mutations introduced to re-pattern methylation at the (CpG)₃ element in the chicken β^A-globin promoter sequence themselves strongly influenced nucleosome formation in reconstituted chromatin. The disruptive effect of cytosine methylation on nucleosome formation was found to be determined by the sequence context of CpG dinucleotides, not just their location in the positioning sequence. Additional mutations indicated that methylation can also promote the occupation of certain nucleosome positions. DNase I analysis demonstrated that these genetic and epigenetic modifications altered the structural characteristics of the (CpG)₃ element. Our findings support a proposal that the intrinsic structural properties of the DNA at the −1.5 site, as occupied by (CpG)₃ in the nucleosome studied, can be decisive for nucleosome formation and stability, and that changes in anisotropic DNA bending or flexibility at this site explain why nucleosome positioning can be exquisitely sensitive to genetic and epigenetic modification of the DNA sequence.     

Nucleosome DNA Bendability Matrix (C. elegans)

Abstract

An original signal extraction procedure is applied to database of 146 base nucleosome core DNA sequences from C. elegans (S. M. Johnson et al. Genome Research 16, 1505–1516, 2006). The positional preferences of various dinucleotides within the 10.4 base nucleosome DNA repeat are calculated, resulting in derivation of the nucleosome DNA bendability matrix of 16x10 elements. A simplified one-line presentation of the matrix (“consensus” repeat) is…A(TTTCCGGAAA)T…. All 6 chromosomes of C. elegans conform to the bendability pattern. The strongest affinity to their respective positions is displayed by dinucleotides AT and CG, separated within the repeat by 5 bases. The derived pattern makes a basis for sequence-directed mapping of nucleosome positions in the genome of C. elegans. As the first complete matrix of bendability available the pattern may serve for iterative calculations of the species-specific matrices of bendability applicable to other genomic sequences.     

Sequence-directed Mapping of Nucleosome Positions

Abstract

A nucleosome DNA sequence probe is designed that combines recently derived RR/YY counter-phase and AA/TT in-phase periodical patterns. A simple nucleosome mapping procedure is introduced for prediction of the nucleosome positions in the sequence of interest, to serve as a guide for experimental studies of the chromatin structure.     

Creation of a SINE enriched library for the isolation of polymorphic (AGC)n microsatellite markers in the bovine genome

Trinucleotide (AGC)_n microsatellites are found as 3′ tails of the artiodactyl short interspersed nuclear element (SINE) A-dimer. We describe a polymerase chain reaction (PCR)-based method for the construction of a plasmid library enriched for SINE (AGC)_n microsatellites. By amplifying Sau3AI inserts with a conserved SINE primer and a flanking vector primer, a 35-fold enrichment of (AGC)_n microsatellites over a conventional genomic library was obtained. The SINE primer was used for both sequencing of AGC-containing inserts and analysis of polymorphism. Twenty-three unique reverse primers were synthesized and used on bovine genomic DNA, 21 producing PCR products of expected size. Five polymorphic (AGC)_n microsatellites with 2–4 alleles each were characterized. Allele sizes differed by a 3 bp motif and lacked the stutter bands associated with dinucleotide repeats. A tendency of increased polymorphism for longer AGC repeat arrays was observed. High stringency selection for positive clones containing eight or more AGC repeats can thus facilitate the isolation of polymorphic (AGC)_n microsatellites, Enrichment for (AGC), microsatellites by SINE-vector PCR can be applied to other bovidae species, such as sheep or goat, containing the artiodactyl SINE elements.     

A partial digest approach to restriction site mapping

We present a new, practical algorithm to resolve the experimental data in restriction site analysis, which is a common technique for mapping DNA. Specifically, we assert that multiple digestions with a single restriction enzyme can provide sufficient information to identify the positions of the restriction sites with high probability. The motivation for the new approach comes from combinatorial results on the number of mutually homeometric sets in one dimension, where two sets ofn points are homeometric if the multiset ofn(n−1)/2 distances they determine are the same. Since experimental data contain errors, we propose algorithms for reconstructing sets from noisy interpoint distances, including the possibility of missing fragments. We analyse the performance of these algorithms under a reasonable probability distribution, establishing a relative error limit ofr=Θ(1/n ²) beyond which our technique becomes infeasible. Through simulations, we establish that our technique is robust enough to reconstruct data with relative errors of up to 7.0% in the measured fragment lengths for typical problems, which appears sufficient for certain biological applications.     

香菇基因组中EST-SSR的构成和分布

从真菌基因组计划网站(FGP)和NCBI网站数据库下载了符合条件的总长度为8.1×106bp的11150条香菇的EST(包括10条cDNA)序列,通过SSRhunter1.3软件结合手工查找,从中发现2.83%即316条EST含有一共469个SSR,平均每17.3kb出现一个EST-SSR。在所有EST-SSR中,三碱基和六碱基SSR出现最多,分别占EST-SSR总数的38.00%和20.00%。出现较多的基序为(A)n、(T)n、(GA)n、(AG)n、(TGA)n、(GAT)n和(TCTTT)n,占所有EST-SSR的35.39%。     

Analysis of Chromatin Structure in the Control Regions of the <Emphasis Type="Italic">Chlamydomonas HSP70A</Emphasis> and <Emphasis Type="Italic">RBCS2</Emphasis> Genes

We have used DNaseI and micrococcal nuclease sensitivity assays to determine the chromatin structures in the control regions of the Chlamydomonas reinhardtii HSP70A and RBCS2 genes. Both genes appear to be organized into nucleosome arrays, which exhibit shorter nucleosome repeat lengths than bulk chromatin. In HSP70A we have identified up to four confined DNaseI hypersensitive sites, three of them localize to the promoter region, a fourth one to the fourth intron. Three hypersensitive sites map close to putative heat shock elements, one close to a CCAAT-box. All hypersensitive sites are located to internucleosomal linkers. Alternative nucleosome positions at half-nucleosomal phasing were constitutively detected in the HSP70A promoter region, indicating local chromatin remodelling. Upon heat shock, dramatic changes in the nucleosome structure of HSP70A were detected that particularly affected the promoter, but also a region within the fourth intron. In contrast, light induction entailed no change in HSP70A chromatin. In the RBCS2 control region we identified a strong DNaseI hypersensitive site that maps close to a CCAAT-box. This site forms the boundary of a nucleosome array with a region of ~700 bp apparently devoid of nucleosomes. This study demonstrates that chromatin structure may be determined readily at fairly high resolution in Chlamydomonas, suggesting this organism as a well-suited model for studying the role of chromatin structure on gene expression in photosynthetic eukaryotes.     

基于遗传片段分析系统的转录起始位点分析技术:从预测到结果评估

【目的】以遗传片段分析仪内标法替代传统放射性标记引物延伸技术进行样本转录起始位点(TSS)分析,并弥补引物延伸技术应用于未知样本缺乏前期预测和后期评估环节,形成一套基于遗传片段分析仪内标法分析未知样品TSS的完整技术方案。【方法】以粘球菌Myxococcus DK1622来源的双拷贝Gro ELs基因为素材;首先从预测出发,利用数据库进行启动子和转录起始位点预测;其次,根据预测结果设计合成荧光标记引物进行靶标m RNA的反转录;再次,应用遗传片段分析技术内标法鉴定分析粘球菌来源的双拷贝Gro ELs基因转录起始位点(TSS)及其丰度;最后,应用正态分布理论进行鉴定结果评估。【结果】明确了转录起始位点的数量、转录丰度及最可能的TSS位点:粘球菌DK1622基因组中Gro EL1拷贝存在1个启动子,TSS位点为TSS_(286);Gro EL2拷贝存在2个启动子,TSS位点分别为TSS_(548)和TSS_(502),其中TSS_(548)转录丰度是TSS_(502)的13.8倍,Gro EL1的TSS_(286)丰度是gro EL2的TSS_(548)丰度的14.3倍。【结论】预测结果指明了实验设计的范围,遗传片段分析仪内标检测法替代传统放射性标记法使实验更加简便、安全、自动、准确,正态分布理论进一步评估了实验结果的可信度,三者接合形成了完善的转录起始位点鉴定技术方案。     

Structure-based Analysis of DNA Sequence Patterns Guiding Nucleosome Positioning in vitro

Abstract

Recent studies of genome-wide nucleosomal organization suggest that the DNA sequence is one of the major determinants of nucleosome positioning. Although the search for underlying patterns encoded in nucleosomal DNA has been going on for about 30 years, our knowledge of these patterns still remains limited. Based on our evaluations of DNA deformation energy, we developed new scoring functions to predict nucleosome positioning. There are three principal differences between our approach and earlier studies: (i) we assume that the length of nucleosomal DNA varies from 146 to 147 bp; (ii) we consider the anisotropic flexibility of pyrimidine-purine (YR) dimeric steps in the context of their neighbors (e.g., YYRR versus RYRY); (iii) we postulate that alternating AT-rich and GC-rich motifs reflect sequence-dependent interactions between histone arginines and DNA in the minor groove. Using these functions, we analyzed 20 nucleosome positions mapped in vitro at single nucleotide resolution (including clones 601, 603, 605, the pGUB plasmid, chicken β-globin and three 5S rDNA genes). We predicted 15 of the 20 positions with 1-bp precision, and two positions with 2-bp precision. The predicted position of the ‘601’ nucleosome (i.e., the optimum of the computed score) deviates from the experimentally determined unique position by no more than 1 bp—an accuracy exceeding that of earlier predictions.

Our analysis reveals a clear heterogeneity of the nucleosomal sequences which can be divided into two groups based on the positioning ‘rules’ they follow. The sequences of one group are enriched by highly deformable YR/YYRR motifs at the minor-groove bending sites SHL ±3.5 and ±5.5, which is similar to the α-satellite sequence used in most crystallized nucleosomes. Apparently, the positioning of these nucleosomes is determined by the interactions between histones H2A/H2B and the terminal parts of nucleosomal DNA. In the other group (that includes the ‘601’ clone) the same YR/YYRR motifs occur predominantly at the sites SHL ±1.5. The interaction between the H3/H4 tetramer and the central part of the nucleosomal DNA is likely to be responsible for the positioning of nucleosomes of this group, and the DNA trajectory in these nucleosomes may differ in detail from the published structures.

Thus, from the stereochemical perspective, the in vitro nucleosomes studied here follow either an X-ray-like pattern (with strong deformations in the terminal parts of nucleosomal DNA), or an alternative pattern (with the deformations occurring predominantly in the central part of the nucleosomal DNA). The results presented here may be useful for genome-wide classification of nucleosomes, linking together structural and thermodynamic characteristics of nucleosomes with the underlying DNA sequence patterns guiding their positions.     

Structural models for the complex of chemotaxis inhibitory protein of Staphylococcus aureus with the C5a receptor

The study presents structural models for the complex of the chemotaxis inhibitory protein of Staphylococcus aureus, CHIPS, and receptor for anaphylotoxin C5a, C5aR. The models are based on the recently found NMR structure of the complex between CHIPS fragment 31-121 and C5aR fragment 7-28, as well as on previous results of molecular modeling of C5aR. Simple and straightforward modeling procedure selected low-energy conformations of the C5aR fragment 8-41 that simultaneously fit the NMR structure of the C5aR 10-18 fragment and properly orient the NMR structure of CHIPS_31-121 relative to C5aR. Extensive repacking of the side chains of CHIPS_31-121 and C5aR_8-41 predicted specific residue-residue interactions on the interface between CHIPS and C5aR. Many of these interactions were rationalized with experimental data obtained by site-directed mutagenesis of CHIPS and C5aR. The models correctly showed that CHIPS binds only to the first binding site of C5a to C5aR not competing with C5a fragment 59-74, which binds the second binding site of C5aR. The models also predict that two elements of CHIPS, fragments 48-58 and 97-111, may be used as structural templates for potential inhibitors of C5a.     

Nucleosome Positioning,Nucleosome Spacing and the Nucleosome Code

Abstract

Nucleosome positioning has been the subject of intense study for many years. The properties of micrococcal nuclease, the enzyme central to these studies, are discussed. The various methods used to determine nucleosome positions in vitro and in vivo are reviewed critically. These include the traditional low resolution method of indirect end-labelling, high resolution methods such as primer extension, monomer extension and nucleosome sequencing, and the high throughput methods for genome-wide analysis (microarray hybridisation and parallel sequencing). It is established that low resolution mapping yields an averaged chromatin structure, whereas high resolution mapping reveals the weighted superposition of all the chromatin states in a cell population. Mapping studies suggest that yeast DNA contains information specifying the positions of nucleosomes and that this code is made use of by the cell. It is proposed that the positioning code facilitates nucleosome spacing by encoding information for multiple alternative overlapping nucleosomal arrays. Such a code might facilitate the shunting of nucleosomes from one array to another by ATP-dependent chromatin remodelling machines.