Genome analysis of Escherichia coli promoter sequences evidences that DNA static curvature plays a more important role in gene transcription than has previously been anticipated

Single-molecule visualization of twin-supercoiled domains generated during transcription

Transcription-coupled supercoiling of DNA is a key factor in chromosome compaction and the regulation of genetic processes in all domains of life. It has become common knowledge that, during transcription, the DNA-dependent RNA polymerase (RNAP) induces positive supercoiling ahead of it (downstream) and negative supercoils in its wake (upstream), as rotation of RNAP around the DNA axis upon tracking its helical groove gets constrained due to drag on its RNA transcript. Here, we experimentally validate this so-called twin-supercoiled-domain model with in vitro real-time visualization at the single-molecule scale. Upon binding to the promoter site on a supercoiled DNA molecule, RNAP merges all DNA supercoils into one large pinned plectoneme with RNAP residing at its apex. Transcription by RNAP in real time demonstrates that up- and downstream supercoils are generated simultaneously and in equal portions, in agreement with the twin-supercoiled-domain model. Experiments carried out in the presence of RNases A and H, revealed that an additional viscous drag of the RNA transcript is not necessary for the RNAP to induce supercoils. The latter results contrast the current consensus and simulations on the origin of the twin-supercoiled domains, pointing at an additional mechanistic cause underlying supercoil generation by RNAP in transcription.

Manipulation of topoisomerase expression inhibits cell division but not growth and reveals a distinctive promoter structure in Synechocystis

In cyanobacteria DNA supercoiling varies over the diurnal cycle and is integrated with temporal programs of transcription and replication. We manipulated DNA supercoiling in Synechocystis sp. PCC 6803 by CRISPRi-based knockdown of gyrase subunits and overexpression of topoisomerase I (TopoI). Cell division was blocked but cell growth continued in all strains. The small endogenous plasmids were only transiently relaxed, then became strongly supercoiled in the TopoI overexpression strain. Transcript abundances showed a pronounced 5'/3' gradient along transcription units, incl. the rRNA genes, in the gyrase knockdown strains. These observations are consistent with the basic tenets of the homeostasis and twin-domain models of supercoiling in bacteria. TopoI induction initially led to downregulation of G+C-rich and upregulation of A+T-rich genes. The transcriptional response quickly bifurcated into six groups which overlap with diurnally co-expressed gene groups. Each group shows distinct deviations from a common core promoter structure, where helically phased A-tracts are in phase with the transcription start site. Together, our data show that major co-expression groups (regulons) in Synechocystis all respond differentially to DNA supercoiling, and suggest to re-evaluate the long-standing question of the role of A-tracts in bacterial promoters.

Spatiotemporal Coupling of DNA Supercoiling and Genomic Sequence Organization-A Timing Chain for the Bacterial Growth Cycle?

In this article we describe the bacterial growth cycle as a closed, self-reproducing, or autopoietic circuit, reestablishing the physiological state of stationary cells initially inoculated in the growth medium. In batch culture, this process of self-reproduction is associated with the gradual decline in available metabolic energy and corresponding change in the physiological state of the population as a function of "travelled distance" along the autopoietic path. We argue that this directional alteration of cell physiology is both reflected in and supported by sequential gene expression along the chromosomal OriC-Ter axis. We propose that during the E. coli growth cycle, the spatiotemporal order of gene expression is established by coupling the temporal gradient of supercoiling energy to the spatial gradient of DNA thermodynamic stability along the chromosomal OriC-Ter axis.

DNA sequence encodes the position of DNA supercoils

The three-dimensional organization of DNA is increasingly understood to play a decisive role in vital cellular processes. Many studies focus on the role of DNA-packaging proteins, crowding, and confinement in arranging chromatin, but structural information might also be directly encoded in bare DNA itself. Here, we visualize plectonemes (extended intertwined DNA structures formed upon supercoiling) on individual DNA molecules. Remarkably, our experiments show that the DNA sequence directly encodes the structure of supercoiled DNA by pinning plectonemes at specific sequences. We develop a physical model that predicts that sequence-dependent intrinsic curvature is the key determinant of pinning strength and demonstrate this simple model provides very good agreement with the data. Analysis of several prokaryotic genomes indicates that plectonemes localize directly upstream of promoters, which we experimentally confirm for selected promotor sequences. Our findings reveal a hidden code in the genome that helps to spatially organize the chromosomal DNA.

RAIL: A new tool for defining bacterial promoter regions

In this issue of the Journal of Bacteriology, Hustmyer and colleagues describe a new method for rapidly generating reporter libraries (Hustmyer citation). This RAIL technique (Rapid Arbitrary PCR Insertion Libraries) uses arbitrary PCR and isothermal DNA assembly to insert random fragments of promoter regions into reporter plasmids, resulting in libraries that can be screened to identify regions required for gene expression. This technique will likely be useful for a number of different genetic applications.

Distinguishing between productive and abortive promoters using a random forest classifier in Mycoplasma pneumoniae

Distinguishing between promoter-like sequences in bacteria that belong to true or abortive promoters, or to those that do not initiate transcription at all, is one of the important challenges in transcriptomics. To address this problem, we have studied the genome-reduced bacterium Mycoplasma pneumoniae, for which the RNAs associated with transcriptional start sites have been recently experimentally identified. We determined the contribution to transcription events of different genomic features: the –10, extended –10 and –35 boxes, the UP element, the bases surrounding the –10 box and the nearest-neighbor free energy of the promoter region. Using a random forest classifier and the aforementioned features transformed into scores, we could distinguish between true, abortive promoters and non-promoters with good –10 box sequences. The methods used in this characterization of promoters can be extended to other bacteria and have important applications for promoter design in bacterial genome engineering.

A repetitive DNA element regulates expression of the Helicobacter pylori sialic acid binding adhesin by a rheostat-like mechanism

During persistent infection, optimal expression of bacterial factors is required to match the ever-changing host environment. The gastric pathogen Helicobacter pylori has a large set of simple sequence repeats (SSR), which constitute contingency loci. Through a slipped strand mispairing mechanism, the SSRs generate heterogeneous populations that facilitate adaptation. Here, we present a model that explains, in molecular terms, how an intergenically located T-tract, via slipped strand mispairing, operates with a rheostat-like function, to fine-tune activity of the promoter that drives expression of the sialic acid binding adhesin, SabA. Using T-tract variants, in an isogenic strain background, we show that the length of the T-tract generates multiphasic output from the sabA promoter. Consequently, this alters the H. pylori binding to sialyl-Lewis x receptors on gastric mucosa. Fragment length analysis of post-infection isolated clones shows that the T-tract length is a highly variable feature in H. pylori. This mirrors the host-pathogen interplay, where the bacterium generates a set of clones from which the best-fit phenotypes are selected in the host. In silico and functional in vitro analyzes revealed that the length of the T-tract affects the local DNA structure and thereby binding of the RNA polymerase, through shifting of the axial alignment between the core promoter and UP-like elements. We identified additional genes in H. pylori, with T- or A-tracts positioned similar to that of sabA, and show that variations in the tract length likewise acted as rheostats to modulate cognate promoter output. Thus, we propose that this generally applicable mechanism, mediated by promoter-proximal SSRs, provides an alternative mechanism for transcriptional regulation in bacteria, such as H. pylori, which possesses a limited repertoire of classical trans-acting regulatory factors.

Structural properties of prokaryotic promoter regions correlate with functional features

The structural properties of the DNA molecule are known to play a critical role in transcription. In this paper, the structural profiles of promoter regions were studied within the context of their diversity and their function for eleven prokaryotic species; Escherichia coli, Klebsiella pneumoniae, Salmonella Typhimurium, Pseudomonas auroginosa, Geobacter sulfurreducens Helicobacter pylori, Chlamydophila pneumoniae, Synechocystis sp., Synechoccocus elongates, Bacillus anthracis, and the archaea Sulfolobus solfataricus. The main anchor point for these promoter regions were transcription start sites identified through high-throughput experiments or collected within large curated databases. Prokaryotic promoter regions were found to be less stable and less flexible than the genomic mean across all studied species. However, direct comparison between species revealed differences in their structural profiles that can not solely be explained by the difference in genomic GC content. In addition, comparison with functional data revealed that there are patterns in the promoter structural profiles that can be linked to specific functional loci, such as sigma factor regulation or transcription factor binding. Interestingly, a novel structural element clearly visible near the transcription start site was found in genes associated with essential cellular functions and growth in several species. Our analyses reveals the great diversity in promoter structural profiles both between and within prokaryotic species. We observed relationships between structural diversity and functional features that are interesting prospects for further research to yet uncharacterized functional loci defined by DNA structural properties.

Role of DNA sequence based structural features of promoters in transcription initiation and gene expression

Regulatory information for transcription initiation is present in a stretch of genomic DNA, called the promoter region that is located upstream of the transcription start site (TSS) of the gene. The promoter region interacts with different transcription factors and RNA polymerase to initiate transcription and contains short stretches of transcription factor binding sites (TFBSs), as well as structurally unique elements. Recent experimental and computational analyses of promoter sequences show that they often have non-B-DNA structural motifs, as well as some conserved structural properties, such as stability, bendability, nucleosome positioning preference and curvature, across a class of organisms. Here, we briefly describe these structural features, the differences observed in various organisms and their possible role in regulation of gene expression.

Assessing diversity of DNA structure-related sequence features in prokaryotic genomes

Prokaryotic genomes are diverse in terms of their nucleotide and oligonucleotide composition as well as presence of various sequence features that can affect physical properties of the DNA molecule. We present a survey of local sequence patterns which have a potential to promote non-canonical DNA conformations (i.e. different from standard B-DNA double helix) and interpret the results in terms of relationships with organisms' habitats, phylogenetic classifications, and other characteristics. Our present work differs from earlier similar surveys not only by investigating a wider range of sequence patterns in a large number of genomes but also by using a more realistic null model to assess significant deviations. Our results show that simple sequence repeats and Z-DNA-promoting patterns are generally suppressed in prokaryotic genomes, whereas palindromes and inverted repeats are over-represented. Representation of patterns that promote Z-DNA and intrinsic DNA curvature increases with increasing optimal growth temperature (OGT), and decreases with increasing oxygen requirement. Additionally, representations of close direct repeats, palindromes and inverted repeats exhibit clear negative trends with increasing OGT. The observed relationships with environmental characteristics, particularly OGT, suggest possible evolutionary scenarios of structural adaptation of DNA to particular environmental niches.

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 Escherichia coli. Fourier analysis of promoter sequences for supercoiling-sensitive genes reveals the tendency in selection of sequences with helical periodicities close to 10nt for relaxation-induced genes and to 11nt 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 the 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 50nt indicate their role in regulation of supercoiling-sensitive genes.

DNA structural properties in the classification of genomic transcription regulation elements

It has been long known that DNA molecules encode information at various levels. The most basic level comprises the base sequence itself and is primarily important for the encoding of proteins and direct base recognition by DNA-binding proteins. A more elusive level consists of the local structural properties of the DNA molecule wherein the DNA sequence only plays an indirect supportive role. These properties are nevertheless an important factor in a large number of biomolecular processes and can be considered as informative signals for the presence of a variety of genomic features. Several recent studies have unequivocally shown the benefit of relying on such DNA properties for modeling and predicting genomic features as diverse as transcription start sites, transcription factor binding sites, or nucleosome occupancy. This review is meant to provide an overview of the key aspects of these DNA conformational and physicochemical properties. To illustrate their potential added value compared to relying solely on the nucleotide sequence in genomics studies, we discuss their application in research on transcription regulation mechanisms as representative cases.

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.

Genome-scale computational analysis of DNA curvature and repeats in Arabidopsis and rice uncovers plant-specific genomic properties

Background

Due to its overarching role in genome function, sequence-dependent DNA curvature continues to attract great attention. The DNA double helix is not a rigid cylinder, but presents both curvature and flexibility in different regions, depending on the sequence. More in depth knowledge of the various orders of complexity of genomic DNA structure has allowed the design of sophisticated bioinformatics tools for its analysis and manipulation, which, in turn, have yielded a better understanding of the genome itself. Curved DNA is involved in many biologically important processes, such as transcription initiation and termination, recombination, DNA replication, and nucleosome positioning. CpG islands and tandem repeats also play significant roles in the dynamics and evolution of genomes.

Results

In this study, we analyzed the relationship between these three structural features within rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana) genomes. A genome-scale prediction of curvature distribution in rice and Arabidopsis indicated that most of the chromosomes of both genomes have maximal chromosomal DNA curvature adjacent to the centromeric region. By analyzing tandem repeats across the genome, we found that frequencies of repeats are higher in regions adjacent to those with high curvature value. Further analysis of CpG islands shows a clear interdependence between curvature value, repeat frequencies and CpG islands. Each CpG island appears in a local minimal curvature region, and CpG islands usually do not appear in the centromere or regions with high repeat frequency. A statistical evaluation demonstrates the significance and non-randomness of these features.

Conclusions

This study represents the first systematic genome-scale analysis of DNA curvature, CpG islands and tandem repeats at the DNA sequence level in plant genomes, and finds that not all of the chromosomes in plants follow the same rules common to other eukaryote organisms, suggesting that some of these genomic properties might be considered as specific to plants.

A temperature-induced narrow DNA curvature range sustains the maximum activity of a bacterial promoter in vitro

Among the molecular strategies bacteria have set up to quickly match their transcriptional program to new environments, changes in sequence-mediated DNA curvature play a crucial role. Bacterial promoters, especially those of mesophilic bacteria, are in general preceded by a curved region. The marked thermosensitivity of curved DNA stretches allows bacteria to rapidly sense outer temperature variations and affects transcription by favoring the binding of activators or repressors. Curved DNA is also able to influence the transcriptional activity of a bacterial promoter directly, without the involvement of trans-acting regulators. This study attempts to quantitatively analyze the role of DNA curvature in thermoregulated gene expression using a real-time in vitro transcription model system based on a specific fluorescence molecular beacon. By analyzing the temperature-dependent expression of a reporter gene in a construct carrying a progressively decreasing bent sequence upstream from the promoter, we show that with a decrease in temperature a narrow curvature range accounts for a significant enhancement of promoter activity. This strengthens the view that DNA curvature-mediated regulation of gene expression is likely a strategy offering fine-tuning control possibilities and that, considering the widespread presence of curved sequences upstream from bacterial promoters, it may represent one of the most primitive forms of gene regulation.

Upstream curved sequences in E. coli are related to the regulation of transcription initiation

The advancement in Escherichia coli genome research has made the information regarding transcription start sites of many genes available. A study relying on the availability of transcription start locations was performed. The first question addressed was what an average DNA curvature profile upstream of genes would look like when these genes are aligned by transcription start sites in comparison to alignment by translation start sites. Since it was hypothesized that curvature plays a role in transcription regulation, the expectation was that curvature measurements relative to transcription starts, rather than translation, should strengthen the signal. Our study justified this expectation. The second question aimed to clarify the relation between DNA curvature and promoter strength. Through clustering based on DNA curvature profiles along promoter regions, a strong positive correlation between the promoter strength and the curved DNA was found. The third question dealt with dinucleotide periodicity in E. coli to see whether a periodicity pattern specific to promoter regions exists. Such unknown pattern might shed new light on transcription regulation mechanisms in E. coli. A sequence periodicity of about 11 bp is characteristic to the whole E. coli genome, and is especially well-expressed in intergenic regions. Here it was shown that regions of the size of about 100-150 bp centered 70-100 bp upstream to transcription starts carry hidden periodicity with a period of about 10.3 bp.

Transcriptional regulation of the capsular polysaccharide biosynthesis locus of streptococcus pneumoniae: a bioinformatic analysis

The polysaccharide capsule of Streptococcus pneumoniae is the main virulence factor, which makes the bacterium resistant to phagocytosis. Expression of capsular polysaccharide must be adjusted at different stages of pneumococcal infection, thus, their transcriptional regulation appears to be crucial. To get insight into the existence of regulatory mechanisms common to most serotypes, a bioinformatic analysis of the DNA region located upstream of the capsular locus was performed. With the exception of serotype 37, the capsular locus is located between dexB and aliA on the pneumococcal chromosome. Up to 26 different sequence organizations were found among pneumococci synthesizing their capsule through a Wzy-polymerase-dependent mechanism, mostly varying according to the presence/absence of distinct insertion elements. As a consequence, only ∼250 bp (including a 107 bp RUP_A element) was conserved in 86 sequences, although only a short (ca. 87 bp) region located immediately upstream of cpsA was strictly conserved in all the sequences analyzed. An exhaustive search for possible operator sequences was done. Interestingly, although the promoter region of serotype 3 isolates completely differs from that of other serotypes, most of the proteins proposed to regulate transcription in serotype 3 pneumococci were also predicted to function as possible regulators in non-serotype 3 S. pneumoniae isolates.

The DNA static curvature has a role in the regulation of the ompS1 porin gene in Salmonella enterica serovar Typhi

The DNA static curvature has been described to play a key role as a regulatory element in the transcription process of several bacterial genes. Here, the role of DNA curvature in the expression of the ompS1 porin gene in Salmonella enterica serovar Typhi is described. The web server mutacurve was used to predict mutations that diminished or restored the extent of DNA curvature in the 5' regulatory region of ompS1. Using these predictions, curvature was diminished by site-directed mutagenesis of only two residues, and curvature was restored by further mutagenesis of the same two residues. Lowering the extent of DNA curvature resulted in an increase in ompS1 expression and in the diminution of the affinity of the silencer proteins H-NS and StpA for the ompS1 5' regulatory region. These mutations were in a region shown not to contain the H-NS nucleation site, consistent with the notion that the effect on expression was due to changes in DNA structural topology.

The GC and window-averaged DNA curvature profile of secondary metabolite gene cluster in Aspergillus fumigatus genome

An immense variety of complex secondary metabolites is produced by filamentous fungi including Aspergillus fumigatus, a main inducer of invasive aspergillosis. The identification of fungal secondary metabolite gene cluster is essential for the characterization of fungal secondary metabolism in terms of genetics and biochemistry through recombinant technologies such as gene disruption and cloning. Most of the prediction methods for secondary metabolite gene cluster severely depend on homology searches. However, homology-based approach has intrinsic limitation to unknown or novel gene cluster. We analyzed the GC and window-averaged DNA curvature profile of 26 secondary metabolite gene clusters in the A. fumigatus genome to find out potential conserved features of secondary metabolite gene cluster. Fifteen secondary metabolite gene clusters showed a conserved pattern in window-averaged DNA curvature profile, that is, the DNA regions including secondary metabolic signature genes such as polyketide synthase, nonribosomal peptide synthase, and/or dimethylallyl tryptophan synthase consisted of window-averaged DNA curvature values lower than 0.18 and these DNA regions were at least 20 kb. Forty percent of secondary metabolite gene clusters with this conserved pattern were related to severe regulation by a transcription factor, LaeA. Our result could be used for identification of other fungal secondary metabolite gene clusters, especially for secondary metabolite gene cluster that is severely regulated by LaeA or other proteins with similar function to LaeA.

RutR is the uracil/thymine-sensing master regulator of a set of genes for synthesis and degradation of pyrimidines

Using the genomic SELEX, a total of six Escherichia coli DNA fragments have been identified, which formed complexes with transcription factor RutR. The RutR regulon was found to include a large number of genes encoding components for not only degradation of pyrimidines but also transport of glutamate, synthesis of glutamine, synthesis of pyrimidine nucleotides and arginine, and degradation of purines. DNase I footprinting indicated that RutR recognizes a palindromic sequence of TTGACCAnnTGGTCAA. The RutR box in P1 promoter of carAB encoding carbamoyl phosphate synthetase, a key enzyme of pyrimidine synthesis, overlaps with the PepA (CarP) repressor binding site, implying competition between RutR and PepA. Adding either uracil or thymine abolished RutR binding in vitro to the carAB P1 promoter. Accordingly, in the rutR-deletion mutant or in the presence of uracil, the activation in vivo of carAB P1 promoter was markedly reduced. Northern blot analysis of the RutR target genes indicated that RutR represses the Gad system genes involved in glutamate-dependent acid resistance and allantoin degradation. Altogether we propose that RutR is the pyrimidine sensor and the master regulator for a large set of the genes involved in the synthesis and degradation of pyrimidines.

The RovA regulons of Yersinia enterocolitica and Yersinia pestis are distinct: evidence that many RovA-regulated genes were acquired more recently than the core genome

RovA is a transcriptional activator of Yersinia invasin, an outer membrane protein involved in bacterial attachment and invasion across the intestinal epithelium. In Y. enterocolitica, a rovA mutant is attenuated for virulence compared with either wild-type or inv mutant strains, indicating that RovA may regulate additional virulence factors. Here, we used microarray analysis to define the RovA regulon. Curiously, there was little overlap between the RovA regulons of Y. enterocolitica and Y. pestis despite the fact that RovA itself is highly conserved between the two species. Some of these differences are explained by the observation that a number of RovA-regulated loci in Y. enterocolitica do not have orthologues in Y. pestis and vice versa, suggesting that RovA established regulatory control over genetic material acquired after the divergence of the species. Electromobility shift assays demonstrated that 15 of these RovA-regulated loci directly interact with RovA, and 11 of these promoters had similar affinity as observed for the inv promoter. H-NS and YmoA are believed to form a transcriptional repression complex on the inv promoter, and several studies indicate that RovA and H-NS have overlapping DNA binding sites. H-NS and YmoA regulated a subset of the RovA-regulated loci. Furthermore, H-NS directly bound to 14 of the 15 promoters bound by RovA. From these data, we hypothesize that RovA generally behaves as an anti-H-NS factor to alleviate transcriptional repression in Y. enterocolitica. A number of recent studies have presented data and a model suggesting that H-NS functions as a transcriptional silencer of horizontally acquired genes. This repression can be selectively relieved by regulators such as RovA, and the observation that nearly all RovA-activated genes are repressed by H-NS is consistent with this model.

Revisiting the relationship between compositional sequence complexity and periodicity

BACKGROUND

Given a big sequence fragment or a set of functionally related sequences we consider two problems of a sequence analysis associated with the given sequence(s). The first problem is to measure sequence complexity (repetitiveness, compactness) to estimate how informative the set as a whole is. Usually an obtained measure should be compared with an appropriate random background calculated using permutation of the given sequences. We propose a novel and effective approach for background information measurement instead of the usual sequence reshuffling. The second problem is to detect a periodic bias to determine if it is one of the set features. Sequence periodicity, when sometimes one has in mind hidden periodicity, is a very basic genomic property. The sequence period of 3, which is considered to characterize coding sequences, and period 10-11, which may be due to the alternation of hydrophobic and hydrophilic amino acids, DNA curvature, and bendability were discovered and described. Searching for periodical biases brought significant results in the study of sequence-dependent nucleosome positioning: nucleosomal sites carry hidden period of about 10.4 bases.

RESULTS

Calculated differences between genomic sequences and background showed high biological relevancy of the method that we proposed in this study. Our algorithm was applied to a few natural and artificial datasets. We constructed a simple "periodic" dataset by replacement of every tenth dinucleotide in each sequence of a trial set by the same dinucleotide "CC". We showed that the method reveals the introduced periodicity and that this periodical pattern carries higher information than in uninterrupted subsequences. An application of the method to the nucleosomal dataset revealed a weak pseudo-periodicity of 10.4 nucleotides confirming previous knowledge. An application of the method to Escherichia coli datasets revealed the well-known periodicity of 3bp as a genic attribute, a secondary genic period slightly larger than 11bp, and an intergenic period a bit smaller than 11bp.

CONCLUSIONS

We reported a novel compositional complexity-based method for sequence analysis. We found that the difference between the sequence complexity of a natural sequence and of background is especially high for a set consisting exclusively of coding sequences. Hidden periodicities were found with no need of any preliminary assumptions regarding a composition of periodic elements. We illustrated the power of the method by studying the sets with known weak periodic properties: a nucleosomal database and sets of different regions of E. coli. We showed that the method conveniently indicated all kinds of periodicity and related features in these sets of DNA sequences.

Blue–white selection of regulatory genes that affect the expression of dehalogenase IVa of Burkholderia cepacia MBA4

We have developed a method for rapid screening of genes that affected the expression of dehalogenase IVa of Burkholderia cepacia MBA4. The promoter region of the dehalogenase gene was used to drive the expression of a beta-galactosidase gene. A plasmid containing this reporter was first electroporated into MBA4, and a Tn5 containing suicidal plasmid was introduced subsequently. The use of electroporation was necessary because Escherichia coli mediated transconjugation was ineffective in plasmid-carrying MBA4. The number of integrants generated was directly proportional to the amount of plasmid DNA used. Integrants with an elevated beta-galactosidase activity were isolated. Mutants with a disruption in a putative iron-transporter gene and in a putative response regulator receiver gene were identified. The basal dehalogenase transcript levels of these mutants were higher than the wild type. These mutants also grow faster than the wild type in chloroacetate-containing medium. This methodology of isolating regulatory mutants is theoretically feasible and convenient for any kinds of bacteria.

A survey of genomic properties for the detection of regulatory polymorphisms

Advances in the computational identification of functional noncoding polymorphisms will aid in cataloging novel determinants of health and identifying genetic variants that explain human evolution. To date, however, the development and evaluation of such techniques has been limited by the availability of known regulatory polymorphisms. We have attempted to address this by assembling, from the literature, a computationally tractable set of regulatory polymorphisms within the ORegAnno database (http://www.oreganno.org). We have further used 104 regulatory single-nucleotide polymorphisms from this set and 951 polymorphisms of unknown function, from 2-kb and 152-bp noncoding upstream regions of genes, to investigate the discriminatory potential of 23 properties related to gene regulation and population genetics. Among the most important properties detected in this region are distance to transcription start site, local repetitive content, sequence conservation, minor and derived allele frequencies, and presence of a CpG island. We further used the entire set of properties to evaluate their collective performance in detecting regulatory polymorphisms. Using a 10-fold cross-validation approach, we were able to achieve a sensitivity and specificity of 0.82 and 0.71, respectively, and we show that this performance is strongly influenced by the distance to the transcription start site.

H-NS, the genome sentinel

Two recent reports have indicated that the H-NS protein in Salmonella enterica serovar Typhimurium has a key role in selectively silencing the transcription of large numbers of horizontally acquired AT-rich genes, including those that make up its major pathogenicity islands. Broadly similar conclusions have emerged from a study of H-NS binding to DNA in Escherichia coli. How do these findings affect our view of H-NS and its ability to influence bacterial evolution?

The distinctive signatures of promoter regions and operon junctions across prokaryotes

Here we show that regions upstream of first transcribed genes have oligonucleotide signatures that distinguish them from regions upstream of genes in the middle of operons. Databases of experimentally confirmed transcription units do not exist for most genomes. Thus, to expand the analyses into genomes with no experimentally confirmed data, we used genes conserved adjacent in evolutionarily distant genomes as representatives of genes inside operons. Likewise, we used divergently transcribed genes as representative examples of first transcribed genes. In model organisms, the trinucleotide signatures of regions upstream of these representative genes allow for operon predictions with accuracies close to those obtained with known operon data (0.8). Signature-based operon predictions have more similar phylogenetic profiles and higher proportions of genes in the same pathways than predicted transcription unit boundaries (TUBs). These results confirm that we are separating genes with related functions, as expected for operons, from genes not necessarily related, as expected for genes in different transcription units. We also test the quality of the predictions using microarray data in six genomes and show that the signature-predicted operons tend to have high correlations of expression. Oligonucleotide signatures should expand the number of tools available to identify operons even in poorly characterized genomes.

Involvement of DNA curvature in intergenic regions of prokaryotes

It is known that DNA curvature plays a certain role in gene regulation. The distribution of curved DNA in promoter regions is evolutionarily preserved, and it is mainly determined by temperature of habitat. However, very little is known on the distribution of DNA curvature in termination sites. Our main objective was to comprehensively analyze distribution of curved sequences upstream and downstream to the coding genes in prokaryotic genomes. We applied CURVATURE software to 170 complete prokaryotic genomes in a search for possible typical distribution of DNA curvature around starts and ends of genes. Performing cluster analyses and other statistical tests, we obtained novel results regarding various factors influencing curvature distribution in intergenic regions, such as growth temperature, A+T composition and genome size. We also analyzed intergenic regions between converging genes in 15 selected genomes. The results show that six genomes presented peaks of curvature excess larger than 3 SDs. Insufficient statistics did not allow us to draw further conclusion. Our hypothesis is that DNA curvature could affect transcription termination in many prokaryotes either directly, through contacts with RNA polymerase, or indirectly, via contacts with some regulatory proteins.