Validating the significance of genomic properties of Chi sites from the distribution of all octamers in Escherichia coli

Biochemical characterization of RecBCD enzyme from an Antarctic Pseudomonas species and identification of its cognate Chi (χ) sequence

Pseudomonas syringae Lz4W RecBCD enzyme, RecBCD^Ps, is a trimeric protein complex comprised of RecC, RecB, and RecD subunits. RecBCD enzyme is essential for P. syringae growth at low temperature, and it protects cells from low temperature induced replication arrest. In this study, we show that the RecBCD^Ps enzyme displays distinct biochemical behaviors. Unlike E. coli RecBCD enzyme, the RecD subunit is indispensable for RecBCD^Ps function. The RecD motor activity is essential for the Chi-like fragments production in P. syringae, highlighting a distinct role for P. syringae RecD subunit in DNA repair and recombination process. Here, we demonstrate that the RecBCD^Ps enzyme recognizes a unique octameric DNA sequence, 5′-GCTGGCGC-3′ (Chi^Ps) that attenuates nuclease activity of the enzyme when it enters dsDNA from the 3′-end. We propose that the reduced translocation activities manifested by motor-defective mutants cause cold sensitivity in P. syrinage; emphasizing the importance of DNA processing and recombination functions in rescuing low temperature induced replication fork arrest.

Selection Effects on the Positioning of Genes and Gene Structures from the Interplay of Replication and Transcription in Bacterial Genomes

Bacterial chromosomes are partly shaped by the functional requirements for efficient replication, which lead to strand bias as commonly characterized by the excess of guanines over cytosines in the leading strand. Gene structures are also highly organized within bacterial genomes as a result of such functional constraints, displaying characteristic positioning and structuring along the genome. Here we analyze the gene structures in completely sequenced bacterial chromosomes to observe the positional constraints on gene orientation, length, and codon usage with regard to the positions of replication origin and terminus. Selection on these gene features is different in regions surrounding the terminus of replication from the rest of the genome, but the selection could be either positive or negative depending on the species, and these positional effects are partly attributed to the A-T enrichment near the terminus. Characteristic gene structuring relative to the position of replication origin and terminus is commonly observed among most bacterial species with circular chromosomes, and therefore we argue that the highly organized gene positioning as well as the strand bias should be considered for genomics studies of bacteria.

The DNA Exonucleases of Escherichia coli

DNA exonucleases, enzymes that hydrolyze phosphodiester bonds in DNA from a free end, play important cellular roles in DNA repair, genetic recombination and mutation avoidance in all organisms. This article reviews the structure, biochemistry, and biological functions of the 17 exonucleases currently identified in the bacterium Escherichia coli. These include the exonucleases associated with DNA polymerases I (polA), II (polB), and III (dnaQ/mutD); Exonucleases I (xonA/sbcB), III (xthA), IV, VII (xseAB), IX (xni/xgdG), and X (exoX); the RecBCD, RecJ, and RecE exonucleases; SbcCD endo/exonucleases; the DNA exonuclease activities of RNase T (rnt) and Endonuclease IV (nfo); and TatD. These enzymes are diverse in terms of substrate specificity and biochemical properties and have specialized biological roles. Most of these enzymes fall into structural families with characteristic sequence motifs, and members of many of these families can be found in all domains of life.

Homologous Recombination-Enzymes and Pathways

Homologous recombination is an ubiquitous process that shapes genomes and repairs DNA damage. The reaction is classically divided into three phases: presynaptic, synaptic, and postsynaptic. In Escherichia coli, the presynaptic phase involves either RecBCD or RecFOR proteins, which act on DNA double-stranded ends and DNA single-stranded gaps, respectively; the central synaptic steps are catalyzed by the ubiquitous DNA-binding protein RecA; and the postsynaptic phase involves either RuvABC or RecG proteins, which catalyze branch-migration and, in the case of RuvABC, the cleavage of Holliday junctions. Here, we review the biochemical properties of these molecular machines and analyze how, in light of these properties, the phenotypes of null mutants allow us to define their biological function(s). The consequences of point mutations on the biochemical properties of recombination enzymes and on cell phenotypes help refine the molecular mechanisms of action and the biological roles of recombination proteins. Given the high level of conservation of key proteins like RecA and the conservation of the principles of action of all recombination proteins, the deep knowledge acquired during decades of studies of homologous recombination in bacteria is the foundation of our present understanding of the processes that govern genome stability and evolution in all living organisms.

Measures of compositional strand bias related to replication machinery and its applications

The compositional asymmetry of complementary bases in nucleotide sequences implies the existence of a mutational or selectional bias in the two strands of the DNA duplex, which is commonly shaped by strand-specific mechanisms in transcription or replication. Such strand bias in genomes, frequently visualized by GC skew graphs, is used for the computational prediction of transcription start sites and replication origins, as well as for comparative evolutionary genomics studies. The use of measures of compositional strand bias in order to quantify the degree of strand asymmetry is crucial, as it is the basis for determining the applicability of compositional analysis and comparing the strength of the mutational bias in different biological machineries in various species. Here, we review the measures of strand bias that have been proposed to date, including the ∆GC skew, the B₁ index, the predictability score of linear discriminant analysis for gene orientation, the signal-to-noise ratio of the oligonucleotide bias, and the GC skew index. These measures have been predominantly designed for and applied to the analysis of replication-related mutational processes in prokaryotes, but we also give research examples in eukaryotes.

The processing of double-stranded DNA breaks for recombinational repair by helicase-nuclease complexes

Double-stranded DNA breaks are prepared for recombinational repair by nucleolytic digestion to form single-stranded DNA overhangs that are substrates for RecA/Rad51-mediated strand exchange. This processing can be achieved through the activities of multiple helicases and nucleases. In bacteria, the function is mainly provided by a stable multi-protein complex of which there are two structural classes; AddAB- and RecBCD-type enzymes. These helicase-nucleases are of special interest with respect to DNA helicase mechanism because they are exceptionally powerful DNA translocation motors, and because they serve as model systems for both single molecule studies and for understanding how DNA helicases can be coupled to other protein machinery. This review discusses recent developments in our understanding of the AddAB and RecBCD complexes, focussing on their distinctive strategies for processing DNA ends. We also discuss the extent to which bacterial DNA end resection mechanisms may parallel those used in eukaryotic cells.

Quantitative analysis of replication-related mutation and selection pressures in bacterial chromosomes and plasmids using generalised GC skew index

BACKGROUND

Due to their bi-directional replication machinery starting from a single finite origin, bacterial genomes show characteristic nucleotide compositional bias between the two replichores, which can be visualised through GC skew or (C-G)/(C+G). Although this polarisation is used for computational prediction of replication origins in many bacterial genomes, the degree of GC skew visibility varies widely among different species, necessitating a quantitative measurement of GC skew strength in order to provide confidence measures for GC skew-based predictions of replication origins.

RESULTS

Here we discuss a quantitative index for the measurement of GC skew strength, named the generalised GC skew index (gGCSI), which is applicable to genomes of any length, including bacterial chromosomes and plasmids. We demonstrate that gGCSI is independent of the window size and can thus be used to compare genomes with different sizes, such as bacterial chromosomes and plasmids. It can suggest the existence of different replication mechanisms in archaea and of rolling-circle replication in plasmids. Correlation of gGCSI values between plasmids and their corresponding host chromosomes suggests that within the same strain, these replicons have reproduced using the same replication machinery and thus exhibit similar strengths of replication strand skew.

CONCLUSIONS

gGCSI can be applied to genomes of any length and thus allows comparative study of replication-related mutation and selection pressures in genomes of different lengths such as bacterial chromosomes and plasmids. Using gGCSI, we showed that replication-related mutation or selection pressure is similar for replicons with similar machinery.

Entropic Profiler - detection of conservation in genomes using information theory

Background

In the last decades, with the successive availability of whole genome sequences, many research efforts have been made to mathematically model DNA. Entropic Profiles (EP) were proposed recently as a new measure of continuous entropy of genome sequences. EP represent local information plots related to DNA randomness and are based on information theory and statistical concepts. They express the weighed relative abundance of motifs for each position in genomes. Their study is very relevant because under or over-representation segments are often associated with significant biological meaning.

Findings

The Entropic Profiler application here presented is a new tool designed to detect and extract under and over-represented DNA segments in genomes by using EP. It allows its computation in a very efficient way by recurring to improved algorithms and data structures, which include modified suffix trees. Available through a web interface and as downloadable source code, it allows to study positions and to search for motifs inside the whole sequence or within a specified range. DNA sequences can be entered from different sources, including FASTA files, pre-loaded examples or resuming a previously saved work. Besides the EP value plots, p-values and z-scores for each motif are also computed, along with the Chaos Game Representation of the sequence.

Conclusion

EP are directly related with the statistical significance of motifs and can be considered as a new method to extract and classify significant regions in genomes and estimate local scales in DNA. The present implementation establishes an efficient and useful tool for whole genome analysis.

RecBCD enzyme and the repair of double-stranded DNA breaks

The RecBCD enzyme of Escherichia coli is a helicase-nuclease that initiates the repair of double-stranded DNA breaks by homologous recombination. It also degrades linear double-stranded DNA, protecting the bacteria from phages and extraneous chromosomal DNA. The RecBCD enzyme is, however, regulated by a cis-acting DNA sequence known as Chi (crossover hotspot instigator) that activates its recombination-promoting functions. Interaction with Chi causes an attenuation of the RecBCD enzyme's vigorous nuclease activity, switches the polarity of the attenuated nuclease activity to the 5' strand, changes the operation of its motor subunits, and instructs the enzyme to begin loading the RecA protein onto the resultant Chi-containing single-stranded DNA. This enzyme is a prototypical example of a molecular machine: the protein architecture incorporates several autonomous functional domains that interact with each other to produce a complex, sequence-regulated, DNA-processing machine. In this review, we discuss the biochemical mechanism of the RecBCD enzyme with particular emphasis on new developments relating to the enzyme's structure and DNA translocation mechanism.

Single molecule studies of homologous recombination

Single molecule methods offer an unprecedented opportunity to examine complex macromolecular reactions that are obfuscated by ensemble averaging. The application of single molecule techniques to study DNA processing enzymes has revealed new mechanistic details that are unobtainable from bulk biochemical studies. Homologous DNA recombination is a multi-step pathway that is facilitated by numerous enzymes that must precisely and rapidly manipulate diverse DNA substrates to repair potentially lethal breaks in the DNA duplex. In this review, we present an overview of single molecule assays that have been developed to study key aspects of homologous recombination and discuss the unique information gleaned from these experiments.

Identification of replication origins in prokaryotic genomes

The availability of hundreds of complete bacterial genomes has created new challenges and simultaneously opportunities for bioinformatics. In the area of statistical analysis of genomic sequences, the studies of nucleotide compositional bias and gene bias between strands and replichores paved way to the development of tools for prediction of bacterial replication origins. Only a few (about 20) origin regions for eubacteria and archaea have been proven experimentally. One reason for that may be that this is now considered as an essentially bioinformatics problem, where predictions are sufficiently reliable not to run labor-intensive experiments, unless specifically needed. Here we describe the main existing approaches to the identification of replication origin (oriC) and termination (terC) loci in prokaryotic chromosomes and characterize a number of computational tools based on various skew types and other types of evidence. We also classify the eubacterial and archaeal chromosomes by predictability of their replication origins using skew plots. Finally, we discuss possible combined approaches to the identification of the oriC sites that may be used to improve the prediction tools, in particular, the analysis of DnaA binding sites using the comparative genomic methods.

DNA double strand break repair and crossing over mediated by RuvABC resolvase and RecG translocase

DNA double-strand breaks threaten the stability of the genome, and yet are induced deliberately during meiosis in order to provoke homologous recombination and generate the crossovers needed to promote faithful chromosome transmission. Crossovers are secured via biased resolution of the double Holliday junction intermediates formed when both ends of the broken chromosome engage an intact homologue. To investigate whether the enzymes catalysing branch migration and resolution of Holliday junctions are directed to favour production of either crossover or noncrossover products, we engineered a genetic system based on DNA breakage induced by the I-SceI endonuclease to detect analogous exchanges in Escherichia coli where the enzymology of recombination is more fully understood. Analysis of the recombinants generated revealed approximately equal numbers of crossover and noncrossover products, regardless of whether repair is mediated via RecG, RuvABC, or the RusA resolvase. We conclude that there little or no control of crossing over at the level of Holliday junction resolution. Thus, if similar resolvases function during meiosis, additional factors must act to bias recombination in favour of crossover progeny.

Very small mobile repeated elements in cyanobacterial genomes

Mobile DNA elements play a major role in genome plasticity and other evolutionary processes, an insight gained primarily through the study of transposons and retrotransposons (generally approximately 1000 nt or longer). These elements spawn smaller parasitic versions (generally >100 nt) that propagate through proteins encoded by the full elements. Highly repeated sequences smaller than 100 nt have been described, but they are either nonmobile or their origins are not known. We have surveyed the genome of the multicellular cyanobacterium, Nostoc punctiforme, and its relatives for small dispersed repeat (SDR) sequences and have identified eight families in the range of from 21 to 27 nucleotides. Three of the families (SDR4, SDR5, and SDR6), despite little sequence similarity, share a common predicted secondary structure, a conclusion supported by patterns of compensatory mutations. The SDR elements are found in a diverse set of contexts, often embedded within tandemly repeated heptameric sequences or within minitransposons. One element (SDR5) is found exclusively within instances of an octamer, HIP1, that is highly over-represented in the genomes of many cyanobacteria. Two elements (SDR1 and SDR4) often are found within copies of themselves, producing complex nested insertions. An analysis of SDR elements within cyanobacterial genomes indicate that they are essentially confined to a coherent subgroup. The evidence indicates that some of the SDR elements, probably working through RNA intermediates, have been mobile in recent evolutionary time, making them perhaps the smallest known mobile elements.

Local Renyi entropic profiles of DNA sequences

Background

In a recent report the authors presented a new measure of continuous entropy for DNA sequences, which allows the estimation of their randomness level. The definition therein explored was based on the Rényi entropy of probability density estimation (pdf) using the Parzen's window method and applied to Chaos Game Representation/Universal Sequence Maps (CGR/USM). Subsequent work proposed a fractal pdf kernel as a more exact solution for the iterated map representation. This report extends the concepts of continuous entropy by defining DNA sequence entropic profiles using the new pdf estimations to refine the density estimation of motifs.

Results

The new methodology enables two results. On the one hand it shows that the entropic profiles are directly related with the statistical significance of motifs, allowing the study of under and over-representation of segments. On the other hand, by spanning the parameters of the kernel function it is possible to extract important information about the scale of each conserved DNA region. The computational applications, developed in Matlab m-code, the corresponding binary executables and additional material and examples are made publicly available at .

Conclusion

The ability to detect local conservation from a scale-independent representation of symbolic sequences is particularly relevant for biological applications where conserved motifs occur in multiple, overlapping scales, with significant future applications in the recognition of foreign genomic material and inference of motif structures.