Coexistence of different base periodicities in prokaryotic genomes as related to DNA curvature, supercoiling, and transcription

We analyzed the periodic patterns in E. coli promoters and compared the distributions of the corresponding patterns in promoters and in the complete genome to elucidate their function. Except the three-base periodicity, coincident with that in the coding regions and growing stronger in the region downstream from the transcriptions start (TS), all other salient periodicities are peaked upstream of TS. We found that helical periodicities with the lengths about B-helix pitch ~10.2-10.5 bp and A-helix pitch ~10.8-11.1 bp coexist in the genomic sequences. We mapped the distributions of stretches with A-, B-, and Z- like DNA periodicities onto E.coli genome. All three periodicities tend to concentrate within non-coding regions when their intensity becomes stronger and prevail in the promoter sequences. The comparison with available experimental data indicates that promoters with the most pronounced periodicities may be related to the supercoiling-sensitive genes.Comment: 23 pages, 6 figures, 2 table

Structural attributes of nucleotide sequences in promoter regions of supercoiling-sensitive genes: how to relate microarray expression data with genomic sequences

The level of supercoiling in the chromosome can affect gene expression. To clarify the basis of supercoiling sensitivity, we analyzed the structural features of nucleotide sequences in the vicinity of promoters for the genes with expression enhanced and decreased in response to loss of chromosomal supercoiling in E. coli. Fourier analysis of promoter sequences for supercoiling-sensitive genes reveals the tendency in selection of sequences with helical periodicities close to 10 nt for relaxation-induced genes and to 11 nt for relaxation-repressed genes. The helical periodicities in the subsets of promoters recognized by RNA polymerase with different sigma factors were also studied. A special procedure was developed for study of correlations between the intensities of periodicities in promoter sequences and the expression levels of corresponding genes. Significant correlations of expression with the AT content and with AT periodicities about 10, 11, and 50 nt indicate their role in regulation of supercoiling-sensitive genes.Comment: 38 pages, 12 figure

Quantumlike Chaos in the Frequency Distributions of the Bases A, C, G, T in Drosophila DNA

Continuous periodogram power spectral analyses of fractal fluctuations of frequency distributions of bases A, C, G, T in Drosophila DNA show that the power spectra follow the universal inverse power-law form of the statistical normal distribution. Inverse power-law form for power spectra of space-time fluctuations is generic to dynamical systems in nature and is identified as self-organized criticality. The author has developed a general systems theory, which provides universal quantification for observed self-organized criticality in terms of the statistical normal distribution. The long-range correlations intrinsic to self-organized criticality in macro-scale dynamical systems are a signature of quantumlike chaos. The fractal fluctuations self-organize to form an overall logarithmic spiral trajectory with the quasiperiodic Penrose tiling pattern for the internal structure. Power spectral analysis resolves such a spiral trajectory as an eddy continuum with embedded dominant wavebands. The dominant peak periodicities are functions of the golden mean. The observed fractal frequency distributions of the Drosophila DNA base sequences exhibit quasicrystalline structure with long-range spatial correlations or self-organized criticality. Modification of the DNA base sequence structure at any location may have significant noticeable effects on the function of the DNA molecule as a whole. The presence of non-coding introns may not be redundant, but serve to organize the effective functioning of the coding exons in the DNA molecule as a complete unit.Comment: 46 pages, 9 figure

A-tract clusters may facilitate DNA packaging in bacterial nucleoid

Molecular mechanisms of bacterial chromosome packaging are still unclear, as bacteria lack nucleosomes or other apparent basic elements of DNA compaction. Among the factors facilitating DNA condensation may be a propensity of the DNA molecule for folding due to its intrinsic curvature. As suggested previously, the sequence correlations in genome reflect such a propensity [Trifonov and Sussman (1980) Proc. Natl Acad. Sci. USA, 77, 3816–3820]. To further elaborate this concept, we analyzed positioning of A-tracts (the sequence motifs introducing the most pronounced DNA curvature) in the Escherichia coli genome. First, we observed that the A-tracts are over-represented and distributed ‘quasi-regularly’ throughout the genome, including both the coding and intergenic sequences. Second, there is a 10–12 bp periodicity in the A-tract positioning indicating that the A-tracts are phased with respect to the DNA helical repeat. Third, the phased A-tracts are organized in ∼100 bp long clusters. The latter feature was revealed with the help of a novel approach based on the Fourier series expansion of the A-tract distance autocorrelation function. Since the A-tracts introduce local bends of the DNA duplex and these bends accumulate when properly phased, the observed clusters would facilitate DNA looping. Also, such clusters may serve as binding sites for the nucleoid-associated proteins that have affinities for curved DNA (such as HU, H-NS, Hfq and CbpA). Therefore, we suggest that the ∼100 bp long clusters of the phased A-tracts constitute the ‘structural code’ for DNA compaction by providing the long-range intrinsic curvature and increasing stability of the DNA complexes with architectural proteins

Sequence periodicity of Escherichia coli is concentrated in intergenic regions

BACKGROUND: Sequence periodicity with a period close to the DNA helical repeat is a very basic genomic property. This genomic feature was demonstrated for many prokaryotic genomes. The Escherichia coli sequences display the period close to 11 base pairs. RESULTS: Here we demonstrate that practically only ApA/TpT dinucleotides contribute to overall dinucleotide periodicity in Escherichia coli. The noncoding sequences reveal this periodicity much more prominently compared to protein-coding sequences. The sequence periodicity of ApC/GpT, ApT and GpC dinucleotides along the Escherichia coli K-12 is found to be located as well mainly within the intergenic regions. CONCLUSIONS: The observed concentration of the dinucleotide sequence periodicity in the intergenic regions of E. coli suggests that the periodicity is a typical property of prokaryotic intergenic regions. We suppose that this preferential distribution of dinucleotide periodicity serves many biological functions; first of all, the regulation of transcription

Prediction of highly expressed genes in microbes based on chromatin accessibility

BACKGROUND: It is well known that gene expression is dependent on chromatin structure in eukaryotes and it is likely that chromatin can play a role in bacterial gene expression as well. Here, we use a nucleosomal position preference measure of anisotropic DNA flexibility to predict highly expressed genes in microbial genomes. We compare these predictions with those based on codon adaptation index (CAI) values, and also with experimental data for 6 different microbial genomes, with a particular interest in experimental data from Escherichia coli. Moreover, position preference is examined further in 328 sequenced microbial genomes. RESULTS: We find that absolute gene expression levels are correlated with the position preference in many microbial genomes. It is postulated that in these regions, the DNA may be more accessible to the transcriptional machinery. Moreover, ribosomal proteins and ribosomal RNA are encoded by DNA having significantly lower position preference values than other genes in fast-replicating microbes. CONCLUSION: This insight into DNA structure-dependent gene expression in microbes may be exploited for predicting the expression of non-translated genes such as non-coding RNAs that may not be predicted by any of the conventional codon usage bias approaches

Hierarchial Information Content, Unguistic Properties, and Protein-Binding Oligomers in Coding and NonCoding DNA Sequences

As much as 97 percent of the DNA in mammalian genomes apparently does not code for protein amino acid sequences. Some of the noncoding DNA is known to function in various gene regulatory roles. The remainder of the noncoding DNA consists mainly of introns the functions of which are largely unknown. In this study large (72,000 base pairs) concatenated sequences of human coding and intronic DNA were analyzed by means of information theoretic and linguistic FORTRAN algorithms on a Sun Sparc 1000 system. The aim was to determine the statistical and linguistic textures of the two categories of DNA as a means of developing a new line of evidence that might provide a basis for an empirical distinction between an intelligent-design origin and an evolutionary origin of genomes. Calculations were run on both the natural DNA sequences and on their randomized counterparts. Similar analyses were performed on the sequenced genome of Mycoplasma genitalium, which consists of 88% coding DNA and does not contain introns. The hierarchical information content of the human (concatenated) coding sequences examined in this study was 1.948-1.951 bitslnucleotide up to the dinucleotide level and 1.912-1.916 bits/nucleotide up to the pentanucleotide level. For intronic DNA the corresponding values were 1.905-1.947 and 1.876-1.901. The Shannon redundancies for the coding DNA sequences are 1.34-1.44% at the dinucleotide level and 2.70-2.83% at the pentanucleotide level. The corresponding values for intronic DNA are 1.36-3.68% and 304-4.84%. The linguistic vocabularies of coding and noncoding DNA sequences of comparable lengths show Significant differences in preferred (standard deviate 2: 3.0) oligomers and avoided (standard deviate s -3.0) oligomers. Intronic sequences exhibit marked modulo 2 periodicities in the spacing of pairs of mirror-symmetric oligomers whereas coding sequences do not show this periodicity. Mirror-complementary oligomers are less abundant than mirror-symmetric and tandemly repeating oligomers in both the coding and noncoding DNAs. Mirror-complementary oligomers occur with higher frequencies in intronic sequences compared to their randomized counterparts than in codonic sequences compared to their randomized counterparts. Coding sequences show marked periodicities modulo 3 in the spacing of tandemly repeating oligomers, whereas the intronic sequences examined in this study do not show this periodicity. The pattern of frequencies of protein-binding sequences in introns differs from that of coding DNA. It is concluded that significant statistical and linguistic differences exist between the coding and intronic DNA of the human genome. These results are consistent with the hypothesis that intronic DNA may playa variety of vital roles in the cell biology of development in multicellular organisms. It is plausible that these roles were present from the very beginning of the existence of organisms on the earth

Conservation of the links between gene transcription and chromosomal organization in the highly reduced genome of Buchnera aphidicola

<p>Abstract</p> <p>Background</p> <p>Genomic studies on bacteria have clearly shown the existence of chromosomal organization as regards, for example, to gene localization, order and orientation. Moreover, transcriptomic analyses have demonstrated that, in free-living bacteria, gene transcription levels and chromosomal organization are mutually influenced. We have explored the possible conservation of relationships between mRNA abundances and chromosomal organization in the highly reduced genome of <it>Buchnera aphidicola</it>, the primary endosymbiont of the aphids, and a close relative to <it>Escherichia coli</it>.</p> <p>Results</p> <p>Using an oligonucleotide-based microarray, we normalized the transcriptomic data by genomic DNA signals in order to have access to inter-gene comparison data. Our analysis showed that mRNA abundances, gene organization (operon) and gene essentiality are correlated in <it>Buchnera </it>(i.e., the most expressed genes are essential genes organized in operons) whereas no link between mRNA abundances and gene strand bias was found. The effect of <it>Buchnera </it>genome evolution on gene expression levels has also been analysed in order to assess the constraints imposed by the obligate symbiosis with aphids, underlining the importance of some gene sets for the survival of the two partners. Finally, our results show the existence of spatial periodic transcriptional patterns in the genome of <it>Buchnera</it>.</p> <p>Conclusion</p> <p>Despite an important reduction in its genome size and an apparent decay of its capacity for regulating transcription, this work reveals a significant correlation between mRNA abundances and chromosomal organization of the aphid-symbiont <it>Buchnera</it>.</p