Short palindromic repetitive DNA elements in enterobacteria: a survey

Non-serogroup O1/O139 agglutinable Vibrio cholerae: a phylogenetically and genealogically neglected yet emerging potential pathogen of clinical relevance

Somatic antigen agglutinable type-1/139 Vibrio cholerae (SAAT-1/139-Vc) members or O1/O139 V. cholerae have been described by various investigators as pathogenic due to their increasing virulence potential and production of choleragen. Reported cholera outbreak cases around the world have been associated with these choleragenic V. cholerae with high case fatality affecting various human and animals. These virulent Vibrio members have shown genealogical and phylogenetic relationship with the avirulent somatic antigen non-agglutinable strains of 1/139 V. cholerae (SANAS-1/139- Vc) or O1/O139 non-agglutinating V. cholerae (O1/O139-NAG-Vc). Reports on implication of O1/O139-NAGVc members in most sporadic cholera/cholera-like cases of diarrhea, production of cholera toxin and transmission via consumption and/or contact with contaminated water/seafood are currently on the rise. Some reported sporadic cases of cholera outbreaks and observed change in nature has also been tracable to these non-agglutinable Vibrio members (O1/O139-NAGVc) yet there is a sustained paucity of research interest on the non-agglutinable V. cholerae members. The emergence of fulminating extraintestinal and systemic vibriosis is another aspect of SANAS-1/139- Vc implication which has received low attention in terms of research driven interest. This review addresses the need to appraise and continually expand research based studies on the somatic antigen non-serogroup agglutinable type-1/139 V. cholerae members which are currently prevalent in studies of water bodies, fruits/vegetables, foods and terrestrial environment. Our opinion is amassed from interest in integrated surveillance studies, management/control of cholera outbreaks as well as diarrhea and other disease-related cases both in the rural, suburban and urban metropolis.

TnpAREP and REP sequences dissemination in bacterial genomes: REP recognition determinants

REP, diverse palindromic DNA sequences found at high copy number in many bacterial genomes, have been attributed important roles in cell physiology but their dissemination mechanisms are poorly understood. They might represent non-autonomous transposable elements mobilizable by TnpAREP, the first prokaryotic domesticated transposase associated with REP. TnpAREP, fundamentally different from classical transposases, are members of the HuH superfamily and closely related to the transposases of the IS200/IS605 family. We previously showed that Escherichia coli TnpAREP processes cognate single stranded REP in vitro and that this activity requires the integrity of the REP structure, in particular imperfect palindromes interrupted by a bulge and preceded by a conserved DNA motif. A second group of REPs rather carry perfect palindromes, raising questions about how the latter are recognized by their cognate TnpAREP. To get insight into the importance of REP structural and sequence determinants in these two groups, we developed an in vitro activity assay coupled to a mutational analysis for three different TnpAREP/REP duos via a SELEX approach. We also tackled the question of how the cleavage site is selected. This study revealed that two TnpAREP groups have co-evolved with their cognate REPs and use different strategies to recognize their REP substrates.

Neutralism versus selectionism: Chargaff's second parity rule, revisited

Of Chargaff's four "rules" on DNA base frequencies, the functional interpretation of his second parity rule (PR2) is the most contentious. Thermophile base compositions (GC%) were taken by Galtier and Lobry (1997) as favoring Sueoka's neutral PR2 hypothesis over Forsdyke's selective PR2 hypothesis, namely that mutations improving local within-species recombination efficiency had generated a genome-wide potential for the strands of duplex DNA to separate and initiate recombination through the "kissing" of the tips of stem-loops. However, following Chargaff's GC rule, base composition mainly reflects a species-specific, genome-wide, evolutionary pressure. GC% could not have consistently followed the dictates of temperature, since it plays fundamental roles in both sustaining species integrity and, through primarily neutral genome-wide mutation, fostering speciation. Evidence for a local within-species recombination-initiating role of base order was obtained with a novel technology that masked the contribution of base composition to nucleic acid folding energy. Forsdyke's results were consistent with his PR2 hypothesis, appeared to resolve some root problems in biology and provided a theoretical underpinning for alignment-free taxonomic analyses using relative oligonucleotide frequencies (k-mer analysis). Moreover, consistent with Chargaff's cluster rule, discovery of the thermoadaptive role of the "purine-loading" of open reading frames made less tenable the Galtier-Lobry anti-selectionist arguments.

Enterobacterial Small Mobile Sequences Carry Open Reading Frames and are Found Intragenically–-Evolutionary Implications for Formation of New Peptides

Intergenic repeat units of 127-bp (RU-1) and 168-bp (RU-2), as well as a newly-found class of 103-bp (RU-3), represent small mobile sequences in enterobacterial genomes present in multiple intergenic regions. These repeat sequences display similarities to eukaryotic miniature inverted-repeat transposable elements (MITE). The RU mobile elements have not been reported to encode amino acid sequences. An in silico approach was used to scan genomes for location of repeat units. RU sequences are found to have open reading frames, which are present in annotated gene loci whereby the RU amino acid sequence is maintained. Gene loci that display repeat units include those that encode large proteins which are part of super families that carry conserved domains and those that carry predicted motifs such as signal peptide sequences and transmembrane domains. A putative exported protein in Y. pestis and a phylogenetically conserved putative inner membrane protein in Salmonella species represent some of the more interesting constructs. We hypothesize that a major outcome of RU open reading frame fusions is the evolutionary emergence of new proteins.

DNA repeat sequences: diversity and versatility of functions

Although discovered decades ago, the molecular identification, the diversity and versatility of functions, and the evolutionary origin of repeat DNA sequences (REPs) containing palindromic units in prokaryotes are now bringing attention to a wide range of biological scientists. A brief account of the current state of the repeat DNA sequences is presented here.

Conformational diversity of single-stranded DNA from bacterial repetitive extragenic palindromes: Implications for the DNA recognition elements of transposases

Repetitive extragenic palindrome (REP)—associated tyrosine transposase enzymes (RAYTs) bind REP DNA domains and catalyze their cleavage. Genomic sequence analyses identify potential noncoding REP sequences associated with RAYT-encoding genes. To probe the conformational space of potential RAYT DNA binding domains, we report here spectroscopic and calorimetric measurements that detect and partially characterize the solution conformational heterogeneity of REP oligonucleotides from six bacterial species. Our data reveal most of these REP oligonucleotides adopt multiple conformations, suggesting that RAYTs confront a landscape of potential DNA substrates in dynamic equilibrium that could be selected, enriched, and/or induced via differential binding. Thus, the transposase-bound DNA motif may not be the predominant conformation of the isolated REP domain. Intriguingly, for several REPs, the circular dichroism spectra suggest guanine tetraplexes as potential alternative or additional RAYT recognition elements, an observation consistent with these REP domains being highly nonrandom, with tetraplex-favoring 5′-G and 3′-C-rich segments. In fact, the conformational heterogeneity of REP domains detected and reported here, including the formation of noncanonical DNA secondary structures, may reflect a general feature required for recognition by RAYT transposases. Based on our biophysical data, we propose guanine tetraplexes as an additional DNA recognition element for binding by RAYT transposase enzymes. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 585–596, 2015.

Translatomics combined with transcriptomics and proteomics reveals novel functional, recently evolved orphan genes in Escherichia coli O157:H7 (EHEC)

Background

Genomes of E. coli, including that of the human pathogen Escherichia coli O157:H7 (EHEC) EDL933, still harbor undetected protein-coding genes which, apparently, have escaped annotation due to their small size and non-essential function. To find such genes, global gene expression of EHEC EDL933 was examined, using strand-specific RNAseq (transcriptome), ribosomal footprinting (translatome) and mass spectrometry (proteome).

Results

Using the above methods, 72 short, non-annotated protein-coding genes were detected. All of these showed signals in the ribosomal footprinting assay indicating mRNA translation. Seven were verified by mass spectrometry. Fifty-seven genes are annotated in other enterobacteriaceae, mainly as hypothetical genes; the remaining 15 genes constitute novel discoveries. In addition, protein structure and function were predicted computationally and compared between EHEC-encoded proteins and 100-times randomly shuffled proteins. Based on this comparison, 61 of the 72 novel proteins exhibit predicted structural and functional features similar to those of annotated proteins. Many of the novel genes show differential transcription when grown under eleven diverse growth conditions suggesting environmental regulation. Three genes were found to confer a phenotype in previous studies, e.g., decreased cattle colonization.

Conclusions

These findings demonstrate that ribosomal footprinting can be used to detect novel protein coding genes, contributing to the growing body of evidence that hypothetical genes are not annotation artifacts and opening an additional way to study their functionality. All 72 genes are taxonomically restricted and, therefore, appear to have evolved relatively recently de novo.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2456-1) contains supplementary material, which is available to authorized users.

Intrastrand triplex DNA repeats in bacteria: a source of genomic instability

Repetitive nucleic acid sequences are often prone to form secondary structures distinct from B-DNA. Prominent examples of such structures are DNA triplexes. We observed that certain intrastrand triplex motifs are highly conserved and abundant in prokaryotic genomes. A systematic search of 5246 different prokaryotic plasmids and genomes for intrastrand triplex motifs was conducted and the results summarized in the ITxF database available online at http://bioinformatics.uni-konstanz.de/utils/ITxF/. Next we investigated biophysical and biochemical properties of a particular G/C-rich triplex motif (TM) that occurs in many copies in more than 260 bacterial genomes by CD and nuclear magnetic resonance spectroscopy as well as in vivo footprinting techniques. A characterization of putative properties and functions of these unusually frequent nucleic acid motifs demonstrated that the occurrence of the TM is associated with a high degree of genomic instability. TM-containing genomic loci are significantly more rearranged among closely related Escherichia coli strains compared to control sites. In addition, we found very high frequencies of TM motifs in certain Enterobacteria and Cyanobacteria that were previously described as genetically highly diverse. In conclusion we link intrastrand triplex motifs with the induction of genomic instability. We speculate that the observed instability might be an adaptive feature of these genomes that creates variation for natural selection to act upon.

The Origin of Mutants Under Selection: How Natural Selection Mimics Mutagenesis (Adaptive Mutation)

Selection detects mutants but does not cause mutations. Contrary to this dictum, Cairns and Foster plated a leaky lac mutant of Escherichia coli on lactose medium and saw revertant (Lac(+)) colonies accumulate with time above a nongrowing lawn. This result suggested that bacteria might mutagenize their own genome when growth is blocked. However, this conclusion is suspect in the light of recent evidence that revertant colonies are initiated by preexisting cells with multiple copies the conjugative F'lac plasmid, which carries the lac mutation. Some plated cells have multiple copies of the simple F'lac plasmid. This provides sufficient LacZ activity to support plasmid replication but not cell division. In nongrowing cells, repeated plasmid replication increases the likelihood of a reversion event. Reversion to lac(+) triggers exponential cell growth leading to a stable Lac(+) revertant colony. In 10% of these plated cells, the high-copy plasmid includes an internal tandem lac duplication, which provides even more LacZ activity—sufficient to support slow growth and formation of an unstable Lac(+) colony. Cells with multiple copies of the F'lac plasmid have an increased mutation rate, because the plasmid encodes the error-prone (mutagenic) DNA polymerase, DinB. Without DinB, unstable and stable Lac(+) revertant types form in equal numbers and both types arise with no mutagenesis. Amplification and selection are central to behavior of the Cairns-Foster system, whereas mutagenesis is a system-specific side effect or artifact caused by coamplification of dinB with lac. Study of this system has revealed several broadly applicable principles. In all populations, gene duplications are frequent stable genetic polymorphisms, common near-neutral mutant alleles can gain a positive phenotype when amplified under selection, and natural selection can operate without cell division when variability is generated by overreplication of local genome subregions.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

GTAG- and CGTC-tagged palindromic DNA repeats in prokaryotes

BACKGROUND

REPs (Repetitive Extragenic Palindromes) are small (20-40 bp) palindromic repeats found in high copies in some prokaryotic genomes, hypothesized to play a role in DNA supercoiling, transcription termination, mRNA stabilization.

RESULTS

We have monitored a large number of REP elements in prokaryotic genomes, and found that most can be sorted into two large DNA super-families, as they feature at one end unpaired motifs fitting either the GTAG or the CGTC consensus. Tagged REPs have been identified in >80 species in 8 different phyla. GTAG and CGTC repeats reside predominantly in microorganisms of the gamma and alpha division of Proteobacteria, respectively. However, the identification of members of both super- families in deeper branching phyla such Cyanobacteria and Planctomycetes supports the notion that REPs are old components of the bacterial chromosome. On the basis of sequence content and overall structure, GTAG and CGTC repeats have been assigned to 24 and 4 families, respectively. Of these, some are species-specific, others reside in multiple species, and several organisms contain different REP types. In many families, most units are close to each other in opposite orientation, and may potentially fold into larger secondary structures. In different REP-rich genomes the repeats are predominantly located between unidirectionally and convergently transcribed ORFs. REPs are predominantly located downstream from coding regions, and many are plausibly transcribed and function as RNA elements. REPs located inside genes have been identified in several species. Many lie within replication and global genome repair genes. It has been hypothesized that GTAG REPs are miniature transposons mobilized by specific transposases known as RAYTs (REP associated tyrosine transposases). RAYT genes are flanked either by GTAG repeats or by long terminal inverted repeats (TIRs) unrelated to GTAG repeats. Moderately abundant families of TIRs have been identified in multiple species.

CONCLUSIONS

CGTC REPs apparently lack a dedicated transposase. Future work will clarify whether these elements may be mobilized by RAYTs or other transposases, and assess if de-novo formation of either GTAG or CGTC repeats type still occurs.

An integrated toolkit for accurate prediction and analysis of cis-regulatory motifs at a genome scale

MOTIVATION

We present an integrated toolkit, BoBro2.0, for prediction and analysis of cis-regulatory motifs. This toolkit can (i) reliably identify statistically significant cis-regulatory motifs at a genome scale; (ii) accurately scan for all motif instances of a query motif in specified genomic regions using a novel method for P-value estimation; (iii) provide highly reliable comparisons and clustering of identified motifs, which takes into consideration the weak signals from the flanking regions of the motifs; and (iv) analyze co-occurring motifs in the regulatory regions.

RESULTS

We have carried out systematic comparisons between motif predictions using BoBro2.0 and the MEME package. The comparison results on Escherichia coli K12 genome and the human genome show that BoBro2.0 can identify the statistically significant motifs at a genome scale more efficiently, identify motif instances more accurately and get more reliable motif clusters than MEME. In addition, BoBro2.0 provides correlational analyses among the identified motifs to facilitate the inference of joint regulation relationships of transcription factors.

AVAILABILITY

The source code of the program is freely available for noncommercial uses at http://code.google.com/p/bobro/.

CONTACT

xyn@bmb.uga.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Breaking and joining single-stranded DNA: the HUH endonuclease superfamily

HUH endonucleases are numerous and widespread in all three domains of life. The major function of these enzymes is processing a range of mobile genetic elements by catalysing cleavage and rejoining of single-stranded DNA using an active-site Tyr residue to make a transient 5'-phosphotyrosine bond with the DNA substrate. These enzymes have a key role in rolling-circle replication of plasmids and bacteriophages, in plasmid transfer, in the replication of several eukaryotic viruses and in various types of transposition. They have also been appropriated for cellular processes such as intron homing and the processing of bacterial repeated extragenic palindromes. Here, we provide an overview of these fascinating enzymes and their functions, using well-characterized examples of Rep proteins, relaxases and transposases, and we explore the molecular mechanisms used in their diverse activities.

Bacterial helicases in post-transcriptional control

Among the five superfamilies of helicases involved in RNA and DNA metabolism, superfamily 2 and superfamily 5 include bacterial RNA-helicases. These enzymes have been shown to be involved in ribosome biogenesis and post-transcriptional gene regulation. Here, we focus on bacterial regulatory mechanisms that are mediated by RNA helicases belonging to superfamily 2, which includes DEAD-box and DEAH-box helicases. Some of these helicases are part of bacterial degradosomes and were shown to unwind RNA duplexes. We will review examples where these enzymes have been implicated in translatability and metabolic stability of bacterial transcripts. This article is part of a Special Issue entitled: The Biology of RNA helicases - Modulation for life.

Multiple pathways of duplication formation with and without recombination (RecA) in Salmonella enterica

Duplications are often attributed to "unequal recombination" between separated, directly repeated sequence elements (>100 bp), events that leave a recombinant element at the duplication junction. However, in the bacterial chromosome, duplications form at high rates (10(-3)-10(-5)/cell/division) even without recombination (RecA). Here we describe 1800 spontaneous lac duplications trapped nonselectively on the low-copy F'(128) plasmid, where lac is flanked by direct repeats of the transposable element IS3 (1258 bp) and by numerous quasipalindromic REP elements (30 bp). Duplications form at a high rate (10(-4)/cell/division) that is reduced only about 11-fold in the absence of RecA. With and without RecA, most duplications arise by recombination between IS3 elements (97%). Formation of these duplications is stimulated by IS3 transposase (Tnp) and plasmid transfer functions (TraI). Three duplication pathways are proposed. First, plasmid dimers form at a high rate stimulated by RecA and are then modified by deletions between IS3 elements (resolution) that leave a monomeric plasmid with an IS3-flanked lac duplication. Second, without RecA, duplications occur by single-strand annealing of DNA ends generated in different sister chromosomes after transposase nicks DNA near participating IS3 elements. The absence of RecA may stimulate annealing by allowing chromosome breaks to persist. Third, a minority of lac duplications (3%) have short (0-36 bp) junction sequences (SJ), some of which are located within REP elements. These duplication types form without RecA, Tnp, or Tra by a pathway in which the palindromic junctions of a tandem inversion duplication (TID) may stimulate deletions that leave the final duplication.

Noncoding RNAs binding to the nucleoid protein HU in Escherichia coli

Some unidentified RNA molecules, together with the nucleoid protein HU, were suggested to be involved in the nucleoid structure of Escherichia coli. HU is a conserved protein known for its role in binding to DNA and maintaining negative supercoils in the latter. HU also binds to a few RNAs, but the full spectrum of its binding targets in the cell is not known. To understand any interaction of HU with RNA in the nucleoid structure, we immunoprecipitated potential HU-RNA complexes from cells and examined bound RNAs by hybridization to whole-genome tiling arrays. We identified associations between HU and 10 new intragenic and intergenic noncoding RNAs (ncRNAs), 2 of which are homologous to the annotated bacterial interspersed mosaic elements (BIMEs) and boxC DNA repeat elements. We confirmed direct binding of HU to BIME RNA in vitro. We also studied the nucleoid shape of HU and two of the ncRNA mutants (nc1 and nc5) by transmission electron microscopy and showed that both HU and the two ncRNAs play a role in nucleoid morphology. We propose that at least two of the ncRNA species complex with HU and help the formation or maintenance of the architecture of the E. coli chromosome. We also observed binding of HU with rRNA and tRNA segments, a few small RNAs, and a distinct small set of mRNAs, although the significance, if any, of these associations is not known.

The processing of repetitive extragenic palindromes: the structure of a repetitive extragenic palindrome bound to its associated nuclease

Extragenic sequences in genomes, such as microRNA and CRISPR, are vital players in the cell. Repetitive extragenic palindromic sequences (REPs) are a class of extragenic sequences, which form nucleotide stem-loop structures. REPs are found in many bacterial species at a high copy number and are important in regulation of certain bacterial functions, such as Integration Host Factor recruitment and mRNA turnover. Although a new clade of putative transposases (RAYTs or TnpA_REP) is often associated with an increase in these repeats, it is not clear how these proteins might have directed amplification of REPs. We report here the structure to 2.6 Å of TnpA_REP from Escherichia coli MG1655 bound to a REP. Sequence analysis showed that TnpA_REP is highly related to the IS200/IS605 family, but in contrast to IS200/IS605 transposases, TnpA_REP is a monomer, is auto-inhibited and is active only in manganese. These features suggest that, relative to IS200/IS605 transposases, it has evolved a different mechanism for the movement of discrete segments of DNA and has been severely down-regulated, perhaps to prevent REPs from sweeping through genomes.

Structuring the bacterial genome: Y1-transposases associated with REP-BIME sequences

REPs are highly repeated intergenic palindromic sequences often clustered into structures called BIMEs including two individual REPs separated by short linker of variable length. They play a variety of key roles in the cell. REPs also resemble the sub-terminal hairpins of the atypical IS200/605 family of insertion sequences which encode Y1 transposases (TnpA(IS200/IS605)). These belong to the HUH endonuclease family, carry a single catalytic tyrosine (Y) and promote single strand transposition. Recently, a new clade of Y1 transposases (TnpA(REP)) was found associated with REP/BIME in structures called REPtrons. It has been suggested that TnpA(REP) is responsible for REP/BIME proliferation over genomes. We analysed and compared REP distribution and REPtron structure in numerous available E. coli and Shigella strains. Phylogenetic analysis clearly indicated that tnpA(REP) was acquired early in the species radiation and was lost later in some strains. To understand REP/BIME behaviour within the host genome, we also studied E. coli K12 TnpA(REP) activity in vitro and demonstrated that it catalyses cleavage and recombination of BIMEs. While TnpA(REP) shared the same general organization and similar catalytic characteristics with TnpA(IS200/IS605) transposases, it exhibited distinct properties potentially important in the creation of BIME variability and in their amplification. TnpA(REP) may therefore be one of the first examples of transposase domestication in prokaryotes.

Impact of small repeat sequences on bacterial genome evolution

Intergenic regions of prokaryotic genomes carry multiple copies of terminal inverted repeat (TIR) sequences, the nonautonomous miniature inverted-repeat transposable element (MITE). In addition, there are the repetitive extragenic palindromic (REP) sequences that fold into a small stem loop rich in G–C bonding. And the clustered regularly interspaced short palindromic repeats (CRISPRs) display similar small stem loops but are an integral part of a complex genetic element. Other classes of repeats such as the REP2 element do not have TIRs but show other signatures. With the current availability of a large number of whole-genome sequences, many new repeat elements have been discovered. These sequences display diverse properties. Some show an intimate linkage to integrons, and at least one encodes a small RNA. Many repeats are found fused with chromosomal open reading frames, and some are located within protein coding sequences. Small repeat units appear to work hand in hand with the transcriptional and/or post-transcriptional apparatus of the cell. Functionally, they are multifaceted, and this can range from the control of gene expression, the facilitation of host/pathogen interactions, or stimulation of the mammalian immune system. The CRISPR complex displays dramatic functions such as an acquired immune system that defends against invading viruses and plasmids. Evolutionarily, mobile repeat elements may have influenced a cycle of active versus inactive genes in ancestral organisms, and some repeats are concentrated in regions of the chromosome where there is significant genomic plasticity. Changes in the abundance of genomic repeats during the evolution of an organism may have resulted in a benefit to the cell or posed a disadvantage, and some present day species may reflect a purification process. The diverse structure, eclectic functions, and evolutionary aspects of repeat elements are described.

Within-genome evolution of REPINs: a new family of miniature mobile DNA in bacteria

Repetitive sequences are a conserved feature of many bacterial genomes. While first reported almost thirty years ago, and frequently exploited for genotyping purposes, little is known about their origin, maintenance, or processes affecting the dynamics of within-genome evolution. Here, beginning with analysis of the diversity and abundance of short oligonucleotide sequences in the genome of Pseudomonas fluorescens SBW25, we show that over-represented short sequences define three distinct groups (GI, GII, and GIII) of repetitive extragenic palindromic (REP) sequences. Patterns of REP distribution suggest that closely linked REP sequences form a functional replicative unit: REP doublets are over-represented, randomly distributed in extragenic space, and more highly conserved than singlets. In addition, doublets are organized as inverted repeats, which together with intervening spacer sequences are predicted to form hairpin structures in ssDNA or mRNA. We refer to these newly defined entities as REPINs (REP doublets forming hairpins) and identify short reads from population sequencing that reveal putative transposition intermediates. The proximal relationship between GI, GII, and GIII REPINs and specific REP-associated tyrosine transposases (RAYTs), combined with features of the putative transposition intermediate, suggests a mechanism for within-genome dissemination. Analysis of the distribution of REPs in a range of RAYT–containing bacterial genomes, including Escherichia coli K-12 and Nostoc punctiforme, show that REPINs are a widely distributed, but hitherto unrecognized, family of miniature non-autonomous mobile DNA.

DNA sequences that copy themselves throughout genomes, and make no specific contribution to reproductive success, are by definition “selfish.” Such DNA is a feature of the genomes of all organisms and evident by virtue of its repetitive nature. In bacteria the predominant repetitive sequences are short (∼20 bp), extragenic, and palindromic. These so-called REP sequences may occur many hundreds of times per genome, but their origins and means of dissemination have been a longstanding mystery. We show that REPs are components of higher-order replicative entities termed REPINs, which are themselves thought to be derived from REP sequences that flanked an ancestral autonomous selfish element. In this ancestral state the REP sequences were likely to have been critical for the movement of the selfish element, but were devoid of any capacity to replicate independently. REPINs, on the other hand, have evolved to have a life of their own, albeit one that exploits—even enslaves—a genetic element upon which their existence depends. REPINs are the ultimate non-autonomous, super-streamlined, selfish element and are widespread among bacteria.

Reductive divergence of enterobacterial repetitive intergenic consensus sequences among Gammaproteobacteria genomes

Enterobacterial repetitive intergenic consensus (ERIC) sequence is a transcription-modulating, nonautonomous, miniature inverted-repeat transposable element. Its origin and the mechanism of highly varying incidences, limited to Enterobacteriaceae and Vibrionaceae, have not been identified. In this study, distribution and divergence of ERICs along bacterial taxonomie units were analyzed. ERICs were found among five families of gammaproteobacteria, with the copy numbers varying with exponential increments. The variability was explained by genus (45%) and species (36%) affiliations, indicating that copy numbers are specific to sub-family taxa. ERICs were interspersed in genomes with considerable divergences. Locations of ERICs in a genome appeared to be strongly conserved in a strain, moderately in a species or a genus, and weakly in a family. ERICs in different species of a genus were from the identical population of sequences while ERICs in different genera of a family were nearly identical. However, ERICs in different families formed distinct monophylectic groups, implying vertical transmission of diverging population of sequences. In spite of large difference in copy numbers, overall intra-genome evolutionary distances among ERICs were similar among different species, except for a few genomes. The exceptions substantiated hypotheses of genetic drifts and horizontal gene transfers of mobility capacity. Therefore, the confined, variable distribution of ERIC could be explained as a two-step evolution: introduction and proliferation of ERIC in one of the progenitors of gammaproteobacteria, followed by vertical transmission under negative selection. Deterioration of sequences and reduction in copy number were concluded to be the predominant patterns in the evolution of ERIC loci.

Identification, variation and transcription of pneumococcal repeat sequences

Background

Small interspersed repeats are commonly found in many bacterial chromosomes. Two families of repeats (BOX and RUP) have previously been identified in the genome of Streptococcus pneumoniae, a nasopharyngeal commensal and respiratory pathogen of humans. However, little is known about the role they play in pneumococcal genetics.

Results

Analysis of the genome of S. pneumoniae ATCC 700669 revealed the presence of a third repeat family, which we have named SPRITE. All three repeats are present at a reduced density in the genome of the closely related species S. mitis. However, they are almost entirely absent from all other streptococci, although a set of elements related to the pneumococcal BOX repeat was identified in the zoonotic pathogen S. suis. In conjunction with information regarding their distribution within the pneumococcal chromosome, this suggests that it is unlikely that these repeats are specialised sequences performing a particular role for the host, but rather that they constitute parasitic elements. However, comparing insertion sites between pneumococcal sequences indicates that they appear to transpose at a much lower rate than IS elements. Some large BOX elements in S. pneumoniae were found to encode open reading frames on both strands of the genome, whilst another was found to form a composite RNA structure with two T box riboswitches. In multiple cases, such BOX elements were demonstrated as being expressed using directional RNA-seq and RT-PCR.

Conclusions

BOX, RUP and SPRITE repeats appear to have proliferated extensively throughout the pneumococcal chromosome during the species' past, but novel insertions are currently occurring at a relatively slow rate. Through their extensive secondary structures, they seem likely to affect the expression of genes with which they are co-transcribed. Software for annotation of these repeats is freely available from ftp://ftp.sanger.ac.uk/pub/pathogens/strep_repeats/.

Fur and the novel regulator YqjI control transcription of the ferric reductase gene yqjH in Escherichia coli

Iron acquisition in aerobic habitats is complicated by the low solubility of ferric hydroxides. Siderophores that bind ferric iron with high affinity are used to mobilize iron. The reduction of ferric iron to the ferrous form can be coupled to the release of iron from siderophores. Iron is also stored intracellularly as a ferric mineral in proteins, such as ferritin, and must be reduced during release. In Escherichia coli, the yqjH gene encodes a putative ferric siderophore reductase that is also part of the Fur regulon. Here we show that YqjH has ferric reductase activity and is required for iron homeostasis in E. coli. Divergently transcribed from yqjH is the yqjI gene, which encodes a novel member of the winged-helix family of transcriptional regulators and also contains an N-terminal extension similar to the Ni(2+)-binding C-terminal tail of SlyD. Deletion of yqjI leads to constitutive high-level activity of the yqjH and yqjI promoters. Purified YqjI binds inverted repeat target sequences within the yqjH and yqjI promoters. We also observed that YqjI-dependent transcriptional repression is reduced when cells are exposed to elevated nickel levels, resulting in increased expression of yqjH and yqjI. YqjI binding to nickel or iron reduces YqjI DNA-binding activity in vitro. Furthermore, we found that elevated nickel stress levels disrupt iron homeostasis in E. coli and that deletion of yqjH increases nickel toxicity. Our results suggest that the YqjI protein controls expression of yqjH to help maintain iron homeostasis under conditions (such as elevated cellular nickel levels) that disrupt iron metabolism.

Interspersed DNA repeats bcr1-bcr18 of Bacillus cereus group bacteria form three distinct groups with different evolutionary and functional patterns

Many short (<400 bp) interspersed sequence repeats exist in bacteria, yet little is known about their origins, mode of generation, or possible function. Here, we present a comprehensive analysis of 18 different previously identified repeated DNA elements, bcr1-bcr18 (Økstad OA, Hegna I, Lindback T, Rishovd AL, Kolstø AB. 1999. Genome organization is not conserved between Bacillus cereus and Bacillus subtilis. Microbiology. 145:621-631.; Tourasse NJ, Helgason E, Økstad OA, Hegna IK, Kolstø AB. 2006. The Bacillus cereus group: novel aspects of population structure and genome dynamics. J Appl Microbiol. 101:579-593.), in 36 sequenced genomes from the Bacillus cereus group of bacteria. This group consists of genetically closely related species with variable pathogenic specificity toward different hosts and includes among others B. anthracis, B. cereus, and B. thuringiensis. The B. cereus group repeat elements could be classified into three categories with different properties: Group A elements (bcr1-bcr3) exhibited highly variable copy numbers ranging from 4 to 116 copies per strain, showed a nonconserved chromosomal distribution pattern between strains, and displayed several features characteristic of mobile elements. Group B repeats (bcr4-bcr6) were present in 0-10 copies per strain and were associated with strain-specific genes and disruptions of genome synteny, implying a possible contribution to genome rearrangements and/or horizontal gene transfer events. bcr5, in particular, was associated with large gene clusters showing resemblance to integrons. In agreement with their potentially mobile nature or involvement in horizontal transfers, the sequences of the repeats from Groups A and B (bcr1-bcr6) followed a phylogeny different from that of the host strains. Conversely, repeats from Group C (bcr7-bcr18) had a conserved chromosomal location and orthologous gene neighbors in the investigated B. cereus group genomes, and their phylogeny matched that of the host chromosome. Several of the group C repeats exhibited a conserved secondary structure or had parts of the structure conserved, possibly indicating functional RNAs. Accordingly, five of the repeats in group C overlapped regions encoding previously characterized riboswitches. Similarly, other group C repeats could represent novel riboswitches, encode small RNAs, and/or constitute other types of regulatory elements with specific biological functions. The current analysis suggests that the multitude of repeat elements identified in the B. cereus group promote genome dynamics and plasticity and could contribute to the flexible and adaptive life style of these bacteria.

The tandem inversion duplication in Salmonella enterica: selection drives unstable precursors to final mutation types

During growth under selection, mutant types appear that are rare in unselected populations. Stress-induced mechanisms may cause these structures or selection may favor a series of standard events that modify common preexisting structures. One such mutation is the short junction (SJ) duplication with long repeats separated by short sequence elements: AB*(CD)*(CD)*E (* = a few bases). Another mutation type, described here, is the tandem inversion duplication (TID), where two copies of a parent sequence flank an inverse-order segment: AB(CD)(E'D'C'B')(CD)E. Both duplication types can amplify by unequal exchanges between direct repeats (CD), and both are rare in unselected cultures but common after prolonged selection for amplification. The observed TID junctions are asymmetric (aTIDs) and may arise from a symmetrical precursor (sTID)-ABCDE(E'D'C'B'A')ABCDE-when sequential deletions remove each palindromic junction. Alternatively, one deletion can remove both sTID junctions to generate an SJ duplication. It is proposed that sTID structures form frequently under all growth conditions, but are usually lost due to their instability and fitness cost. Selection for increased copy number helps retain the sTID and favors deletions that remodel junctions, improve fitness, and allow higher amplification. Growth improves with each step in formation of an SJ or aTID amplification, allowing selection to favor completion of the mutation process.

Identification and characterization of repetitive extragenic palindromes (REP)-associated tyrosine transposases: implications for REP evolution and dynamics in bacterial genomes

Background

Bacterial repetitive extragenic palindromes (REPs) compose a distinct group of genomic repeats. They usually occur in high abundance (>100 copies/genome) and are often arranged in composite repetitive structures - bacterial interspersed mosaic elements (BIMEs). In BIMEs, regularly spaced REPs are present in alternating orientations. BIMEs and REPs have been shown to serve as binding sites for several proteins and suggested to play role in chromosome organization and transcription termination. Their origins are, at present, unknown.

Results

In this report, we describe a novel class of putative transposases related to IS200/IS605 transposase family and we demonstrate that they are obligately associated with bacterial REPs. Open reading frames coding for these REP-associated tyrosine transposases (RAYTs) are always flanked by two REPs in inverted orientation and thus constitute a unit reminiscent of typical transposable elements. Besides conserved residues involved in catalysis of DNA cleavage, RAYTs carry characteristic structural motifs that are absent in typical IS200/IS605 transposases. DNA sequences flanking rayt genes are in one third of examined cases arranged in modular BIMEs. RAYTs and their flanking REPs apparently coevolve with each other. The rayt genes themselves are subject to rapid evolution, substantially exceeding the substitution rate of neighboring genes. Strong correlation was found between the presence of a particular rayt in a genome and the abundance of its cognate REPs.

Conclusions

In light of our findings, we propose that RAYTs are responsible for establishment of REPs and BIMEs in bacterial genomes, as well as for their exceptional dynamics and species-specifity. Conversely, we suggest that BIMEs are in fact a special type of nonautonomous transposable elements, mobilizable by RAYTs.

Duplication frequency in a population of Salmonella enterica rapidly approaches steady state with or without recombination

Tandem duplications are among the most common mutation events. The high loss rate of duplication suggested that the frequency of duplications in a bacterial population (1/1000) might reflect a steady state dictated by relative rates of formation (k(F)) and loss (k(L)). This possibility was tested for three genetic loci. Without homologous recombination (RecA), duplication loss rate dropped essentially to zero, but formation rate decreased only slightly and a steady state was still reached rapidly. Under all conditions, steady state was reached faster than predicted by formation and loss rates alone. A major factor in determining steady state proved to be the fitness cost, which can exceed 40% for some genomic regions. Depending on the region tested, duplications reached 40-98% of the steady-state frequency within 30 generations-approximately the growth required for a single cell to produce a saturated overnight culture or form a large colony on solid medium (10(9) cells). Long-term bacterial populations are stably polymorphic for duplications of every region of their genome. These polymorphisms contribute to rapid genetic adaptation by providing frequent preexisting mutations that are beneficial whenever imposed selection favors increases in some gene activity. While the reported results were obtained with the bacterium Salmonella enterica, the genetic implications seem likely to be of broader biological relevance.

Small RNA-induced differential degradation of the polycistronic mRNA iscRSUA

Most polycistronic genes are expressed in a single transcript, in which each cistron produces a fixed amount of protein. In this report, we show the first example of differential degradation of a polycistronic gene induced by a small regulatory RNA (sRNA). Our data show that the iron-responsive sRNA, RyhB, binds to the second cistron of the polycistronic mRNA, iscRSUA, which encodes the necessary machinery for biosynthesis of Fe-S clusters, and promotes the cleavage of the downstream iscSUA transcript. This cleavage gives rise to the remaining 5'-section of the transcript encoding IscR, a transcriptional regulator responsible for activation and repression of several genes depending on the cellular Fe-S level. Our data indicate that the iscR transcript is stable and that translation is active. The stability of the iscR transcript depends on a 111-nucleotide long non-translated RNA section located between iscR and iscS, which forms a strong repetitive extragenic palindromic secondary structure and may protect against ribonucleases degradation. This novel regulation shows how sRNAs and mRNA structures can work together to modulate the transcriptional response to a specific stress.

Intergenic regions of Borrelia plasmids contain phylogenetically conserved RNA secondary structure motifs

BACKGROUND

Borrelia species are unusual in that they contain a large number of linear and circular plasmids. Many of these plasmids have long intergenic regions. These regions have many fragmented genes, repeated sequences and appear to be in a state of flux, but they may serve as reservoirs for evolutionary change and/or maintain stable motifs such as small RNA genes.

RESULTS

In an in silico study, intergenic regions of Borrelia plasmids were scanned for phylogenetically conserved stem loop structures that may represent functional units at the RNA level. Five repeat sequences were found that could fold into stable RNA-type stem loop structures, three of which are closely linked to protein genes, one of which is a member of the Borrelia lipoprotein_1 super family genes and another is the complement regulator-acquiring surface protein_1 (CRASP-1) family. Modeled secondary structures of repeat sequences display numerous base-pair compensatory changes in stem regions, including C-G-->A-U transversions when orthologous sequences are compared. Base-pair compensatory changes constitute strong evidence for phylogenetic conservation of secondary structure.

CONCLUSION

Intergenic regions of Borrelia species carry evolutionarily stable RNA secondary structure motifs. Of major interest is that some motifs are associated with protein genes that show large sequence variability. The cell may conserve these RNA motifs whereas allow a large flux in amino acid sequence, possibly to create new virulence factors but with associated RNA motifs intact.

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.

Fluorescent in situ hybridization and flow cytometry as tools to evaluate the treatments for the control of slime-forming enterobacteria in paper mills

Slime formation is a serious problem nowadays in the paper industry. Some enterobacteria are associated with the formation of slime deposits in paper and board mills. Detection and characterization of slime forming bacteria, belonging to the genus Enterobacter, Raoultella, and Klebsiella have been achieved by fluorescence in situ hybridization (FISH), using one probe based on the enterobacterial repetitive intergenic consensus sequence and other two rRNA targeted oligonucleotide probes. The effects of three kinds of antimicrobiological products (biocides, dispersants, and enzymes) on these enterobacterial cells were analyzed by flow cytometry (FC). Biocides B: utrol 1009 and 1072 were the most effective microbiocides against all enterobacterial cells analyzed, reaching 90% of dead bacteria after 24 h. However, the enzymatic treatment (Buzyme) was not equally efficient on enterobacteria and its microbiocide capacity varied depending on the type of microorganism. FISH and FC were effective tools to detect important slime forming enterobacteria and to select specific treatments to control microbial problems in the paper industry.

Systematic identification of stem-loop containing sequence families in bacterial genomes

BACKGROUND

Analysis of non-coding sequences in several bacterial genomes brought to the identification of families of repeated sequences, able to fold as secondary structures. These sequences have often been claimed to be transcribed and fulfill a functional role. A previous systematic analysis of a representative set of 40 bacterial genomes produced a large collection of sequences, potentially able to fold as stem-loop structures (SLS). Computational analysis of these sequences was carried out by searching for families of repetitive nucleic acid elements sharing a common secondary structure.

RESULTS

The initial clustering procedure identified clusters of similar sequences in 29 genomes, corresponding to about 1% of the whole population. Sequences selected in this way have a substantially higher aptitude to fold into a stable secondary structure than the initial set. Removal of redundancies and regrouping of the selected sequences resulted in a final set of 92 families, defined by HMM analysis. 25 of them include all well-known SLS containing repeats and others reported in literature, but not analyzed in detail. The remaining 67 families have not been previously described. Two thirds of the families share a common predicted secondary structure and are located within intergenic regions.

CONCLUSION

Systematic analysis of 40 bacterial genomes revealed a large number of repeated sequence families, including known and novel ones. Their predicted structure and genomic location suggest that, even in compact bacterial genomes, a relatively large fraction of the genome consists of non-protein-coding sequences, possibly functioning at the RNA level.

Exploring the evolution of the Bacillus cereus group repeat element bcr1 by comparative genome analysis of closely related strains

bcr1 is a chromosomal approximately 155 bp repeated element found uniquely and ubiquitously in the Bacillus cereus group of Gram-positive bacteria; it exhibits several features characteristic of mobile elements, including a variable distribution pattern between strains. Here, highly similar bcr1 elements in non-conserved genomic loci are identified in a set of closely related B. cereus and Bacillus thuringiensis strains near the Bacillus anthracis phylogenetic cluster. It is also shown that bcr1 may be present on small RNA transcripts in the 100-400 bp size range. In silico folding of bcr1 at the RNA level indicated that transcripts may form a double-hairpin-like structure predicted to have high structural stability. A functional role of bcr1 at the RNA level is supported by multiple cases of G-U base-pairing, and compensatory mutations maintaining structural stability of the RNA fold. In silico folding at the DNA level produced similar predicted structures, with the potential to form a cruciform structure at open DNA complexes. The predicted structural stability was greater for bcr1 elements showing high sequence identities to bcr1 elements in non-conserved chromosomal loci in other strains, relative to other bcr1 copies. bcr1 mobility could thus be dependent on the formation of a stable DNA or RNA intermediate. Furthermore, bcr1 elements potentially encoding structurally stable and less stable transcripts were phylogenetically intermixed, indicating that loss of bcr1 mobility may have occurred multiple times during evolution. Repeated elements with similar features in other bacteria have been shown to provide functions such as mRNA stabilization, transcription termination and/or promoter function. Similarly, bcr1 may constitute a mobile element which occasionally gains a function when it enters an appropriate chromosomal locus.

Small mobile sequences in bacteria display diverse structure/function motifs

Small repeat sequences in bacterial genomes, which represent non-autonomous mobile elements, have close similarities to archaeon and eukaryotic miniature inverted repeat transposable elements. These repeat elements are found in both intergenic and intragenic chromosomal regions, and contain an array of diverse motifs. These can include DNA sequences containing an integration host factor binding site and a proposed DNA methyltransferase recognition site, transcribed RNA secondary structural motifs, which are involved in mRNA regulation, and translated open reading frames found fused to other open reading frames. Some bacterial mobile element fusions are in evolutionarily conserved protein and RNA genes. Others might represent or lead to creation of new protein genes. Here we review the remarkable properties of these small bacterial mobile elements in the context of possible beneficial roles resulting from random insertions into the genome.

Structural organization and functional properties of miniature DNA insertion sequences in yersiniae

YPALs (Yersinia palindromic sequences) are miniature DNA insertions scattered along the chromosomes of yersiniae. The spread of these intergenic repeats likely occurred via transposition, as suggested by the presence of target site duplications at their termini and the identification of syntenic chromosomal regions which differ in the presence/absence of YPAL DNA among Yersinia strains. YPALs tend to be inserted closely downstream from the stop codon of flanking genes, and many YPAL targets overlap rho-independent transcriptional terminator-like sequences. This peculiar pattern of insertion supports the hypothesis that most of these repeats are cotranscribed with upstream sequences into mRNAs. YPAL RNAs fold into stable hairpins which may modulate mRNA decay. Accordingly, we found that YPAL-positive transcripts accumulate in Yersinia enterocolitica cells at significantly higher levels than homologous transcripts lacking YPAL sequences in their 3' untranslated region.

Stem-loop structures in prokaryotic genomes

BACKGROUND

Prediction of secondary structures in the expressed sequences of bacterial genomes allows to investigate spontaneous folding of the corresponding RNA. This is particularly relevant in untranslated mRNA regions, where base pairing is less affected by interactions with the translation machinery. Relatively large stem-loops significantly contribute to the formation of more complex secondary structures, often important for the activity of sequence elements controlling gene expression.

RESULTS

Systematic analysis of the distribution of stem-loop structures (SLSs) in 40 wholly-sequenced bacterial genomes is presented. SLSs were searched as stems measuring at least 12 bp, bordering loops 5 to 100 nt in length. G-U pairing in the stems was allowed. SLSs found in natural genomes are constantly more numerous and stable than those expected to randomly form in sequences of comparable size and composition. The large majority of SLSs fall within protein-coding regions but enrichment of specific, non random, SLS sub-populations of higher stability was observed within the intergenic regions of the chromosomes of several species. In low-GC firmicutes, most higher stability intergenic SLSs resemble canonical rho-independent transcriptional terminators, but very frequently feature at the 5'-end an additional A-rich stretch complementary to the 3' uridines. In all species, a clearly biased SLS distribution was observed within the intergenic space, with most concentrating at the 3'-end side of flanking CDSs. Some intergenic SLS regions are members of novel repeated sequence families.

CONCLUSION

In depth analysis of SLS features and distribution in 40 different bacterial genomes showed the presence of non random populations of such structures in all species. Many of these structures are plausibly transcribed, and might be involved in the control of transcription termination, or might serve as RNA elements which can enhance either the stability or the turnover of cotranscribed mRNAs. Three previously undescribed families of repeated sequences were found in Yersiniae, Bordetellae and Enterococci.

Enterobacterial repetitive intergenic consensus sequence repeats in yersiniae: genomic organization and functional properties

Genome-wide analyses carried out in silico revealed that the DNA repeats called enterobacterial repetitive intergenic consensus sequences (ERICs), which are present in several Enterobacteriaceae, are overrepresented in yersiniae. From the alignment of DNA regions from the wholly sequenced Yersinia enterocolitica 8081 and Yersinia pestis CO92 strains, we could establish that ERICs are miniature mobile elements whose insertion leads to duplication of the dinucleotide TA. ERICs feature long terminal inverted repeats (TIRs) and can fold as RNA into hairpin structures. The proximity to coding regions suggests that most Y. enterocolitica ERICs are cotranscribed with flanking genes. Elements which either overlap or are located next to stop codons are preferentially inserted in the same (or B) orientation. In contrast, ERICs located far apart from open reading frames are inserted in the opposite (or A) orientation. The expression of genes cotranscribed with A- and B-oriented ERICs has been monitored in vivo. In mRNAs spanning B-oriented ERICs, upstream gene transcripts accumulated at lower levels than downstream gene transcripts. This difference was abolished by treating cells with chloramphenicol. We hypothesize that folding of B-oriented elements is impeded by translating ribosomes. Consequently, upstream RNA degradation is triggered by the unmasking of a site for the RNase E located in the right-hand TIR of ERIC. A-oriented ERICs may act in contrast as upstream RNA stabilizers or may have other functions. The hypothesis that ERICs act as regulatory RNA elements is supported by analyses carried out in Yersinia strains which either lack ERIC sequences or carry alternatively oriented ERICs at specific loci.

Characterization of the molecular mechanisms involved in the differential production of erythrose-4-phosphate dehydrogenase, 3-phosphoglycerate kinase and class II fructose-1,6-bisphosphate aldolase in Escherichia coli

A gapA-pgk gene tandem coding the glyceraldehyde 3-phosphate dehydrogenase and 3-phosphoglycerate kinase, is most frequently found in bacteria. However, in Enterobacteriaceae, gapA is replaced by an epd open reading frame (ORF) coding an erythrose-4-phosphate dehydrogenase and an fbaA ORF coding the class II fructose-1,6-bisphosphate aldolase follows pgk. Although epd expression is very low in Escherichia coli, we show that, in the presence of glucose, the 3 epd, pgk and fbaA ORFs are efficiently cotranscribed from promoter epd P0. Conservation of promoter epd P0 is likely due to its important role in modulation of the metabolic flux during glycolysis and gluconeogenesis. As a consequence, we found that the epd translation initiation region and ORF have been adapted in order to limit epd translation and to create an efficient RNase E entry site. We also show that fbaA is cotranscribed with pgk, from promoter epd P0 or an internal pgk P1 promoter of the extended -10 class. The differential expression of pgk and fbaA also depends upon an RNase E segmentation process, leading to individual mRNAs with different stabilities. The secondary structures of the RNA regions containing the RNase E sites were experimentally determined which brings important information on the structural features of RNase E ectopic sites.

Palindromic repetitive DNA elements with coding potential in Methanocaldococcus jannaschii

We have identified 141 novel palindromic repetitive elements in the genome of euryarchaeon Methanocaldococcus jannaschii. The total length of these elements is 14.3kb, which corresponds to 0.9% of the total genomic sequence and 6.3% of all extragenic regions. The elements can be divided into three groups (MJRE1-3) based on the sequence similarity. The low sequence identity within each of the groups suggests rather old origin of these elements in M. jannaschii. Three MJRE2 elements were located within the protein coding regions without disrupting the coding potential of the host genes, indicating that insertion of repeats might be a widespread mechanism to enhance sequence diversity in coding regions.

Systematic determination of the mosaic structure of bacterial genomes: species backbone versus strain-specific loops

Background

Public databases now contain multitude of complete bacterial genomes, including several genomes of the same species. The available data offers new opportunities to address questions about bacterial genome evolution, a task that requires reliable fine comparison data of closely related genomes. Recent analyses have shown, using pairwise whole genome alignments, that it is possible to segment bacterial genomes into a common conserved backbone and strain-specific sequences called loops.

Results

Here, we generalize this approach and propose a strategy that allows systematic and non-biased genome segmentation based on multiple genome alignments. Segmentation analyses, as applied to 13 different bacterial species, confirmed the feasibility of our approach to discern the 'mosaic' organization of bacterial genomes. Segmentation results are available through a Web interface permitting functional analysis, extraction and visualization of the backbone/loops structure of documented genomes. To illustrate the potential of this approach, we performed a precise analysis of the mosaic organization of three E. coli strains and functional characterization of the loops.

Conclusion

The segmentation results including the backbone/loops structure of 13 bacterial species genomes are new and available for use by the scientific community at the URL: .

The bzd gene cluster, coding for anaerobic benzoate catabolism, in Azoarcus sp. strain CIB

We report here that the bzd genes for anaerobic benzoate degradation in Azoarcus sp. strain CIB are organized as two transcriptional units, i.e., a benzoate-inducible catabolic operon, bzdNOPQMSTUVWXYZA, and a gene, bzdR, encoding a putative transcriptional regulator. The last gene of the catabolic operon, bzdA, has been expressed in Escherichia coli and encodes the benzoate-coenzyme A (CoA) ligase that catalyzes the first step in the benzoate degradation pathway. The BzdA enzyme is able to activate a wider range of aromatic compounds than that reported for other previously characterized benzoate-CoA ligases. The reduction of benzoyl-CoA to a nonaromatic cyclic intermediate is carried out by a benzoyl-CoA reductase (bzdNOPQ gene products) detected in Azoarcus sp. strain CIB extracts. The bzdW, bzdX, and bzdY gene products show significant similarity to the hydratase, dehydrogenase, and ring-cleavage hydrolase that act sequentially on the product of the benzoyl-CoA reductase in the benzoate catabolic pathway of Thauera aromatica. Benzoate-CoA ligase assays and transcriptional analyses based on lacZ-reporter fusions revealed that benzoate degradation in Azoarcus sp. strain CIB is subject to carbon catabolite repression by some organic acids, indicating the existence of a physiological control that connects the expression of the bzd genes to the metabolic status of the cell.

PriA is essential for viability of the Escherichia coli topoisomerase IV parE10(Ts) mutant

The parE10(Ts) mutation, which renders Escherichia coli thermosensitive for growth by inactivation of the essential E. coli topoisomerase topo IV, is lethal at all temperatures when PriA, the main replication restart protein, is absent. This lethality is suppressed by the activation of a PriA-independent replication restart pathway (dnaC809 mutation). This result suggests that topo IV acts prior to full-chromosome replication completion.

The RNA degradosome and poly(A) polymerase of Escherichia coli are required in vivo for the degradation of small mRNA decay intermediates containing REP-stabilizers

In Escherichia coli, REP-stabilizers are structural elements in polycistronic messages that protect 5'-proximal cistrons from 3'-->5' exonucleolytic degradation. The stabilization of a protected cistron can be an important determinant in the level of gene expression. Our results suggest that RNase E, an endoribonuclease, initiates the degradation of REP-stabilized mRNA. However, subsequent degradation of mRNA fragments containing a REP-stabilizer poses a special challenge to the mRNA degradation machinery. Two enzymes, the DEAD-box RNA helicase, RhlB and poly(A) polymerase (PAP) are required to facilitate the degradation of REP-stabilizers by polynucleotide phosphorylase (PNPase). This is the first in vivo evidence that these enzymes are required for the degradation of REP-stabilizers. Furthermore, our results show that REP degradation by RhlB and PNPase requires their association with RNase E as components of the RNA degradosome, thus providing the first in vivo evidence that this ribonucleolytic multienzyme complex is involved in the degradation of structured mRNA fragments.

Transposases are responsible for the target specificity of IS1397 and ISKpn1 for two different types of palindromic units (PUs)

Insertion sequences (IS)1397 and ISKpn1, found in Escherichia coli and Klebsiella pneumoniae, respectively, are IS3 family members that insert specifically into short palindromic repeated sequences (palindromic units or PUs). In this paper, we first show that although PUs are naturally absent from extrachromosomal elements, both ISs are able to transpose from the chromosome or from a plasmid into PUs artificially introduced into target plasmids. We also show that ISKpn1 target specificity is restricted to K.pneumoniae Z1 PU type, whereas IS1397 target specificity is less stringent since the IS targets the three E.coli Y, Z1 and Z2 PU types indifferently. Experiments of transposition of both ISs driven by both transposases demonstrate that the inverted repeats flanking the ISs are not responsible for this target specificity, which is entirely due to the transposase itself. Implications on ISs evolution are presented.

Polarisation of prokaryotic chromosomes

In many prokaryotes, asymmetrical mutational or selective pressures have caused compositional skews between complementary strands of replication arms, especially sensitive in the distribution of guanine and cytosine. In Escherichia coli, most of the guanine/cytosine skew is caused by mutation rates differing on leading and lagging strands, but contribution of skewed functionally important guanine-rich motifs (Chi and Rag sites), which control chromosome repair or positioning, is noticeable. Interference between replication and gene expression plays a minor role. The situation may be different in other bacteria. Studies of chromosome processing and bacterial taxonomy might profit from consideration of chromosome polarisation.

Palindromic unit-independent transposition of IS1397 in Yersinia pestis

Palindromic units (PUs) are intergenic repeated sequences scattered over the chromosomes of Escherichia coli and several other enterobacteria. In the latter, IS1397, an E. coli insertion sequence specific to PUs, transposes into PUs with sequences close to the E. coli consensus. Reasons for this insertion specificity can relate to either a direct recognition of the target (by its sequence or its structure) by the transposase or an interaction between a specific host protein and the PU target DNA sequence. In this study, we show that for Yersinia pestis, a species deprived of PUs, IS1397 can transpose onto its chromosome, with transpositional hot spots. Our results are in favor of a direct recognition of target DNA by IS1397 transposase.

Proliferation and deterioration of Rickettsia palindromic elements

It has been suggested that Rickettsia Palindromic Elements (RPEs) have evolved as selfish DNA that mediate protein sequence evolution by being targeted to genes that code for RNA and proteins. Here, we have examined the phylogenetic depth of two RPEs that are located close to the genes encoding elongation factors Tu (tuf) and G (fus) in Rickettsia. An exceptional organization of the elongation factor genes was found in all 11 species examined, with complete or partial RPEs identified downstream of the tuf gene (RPE-tuf) in six species and of the fus gene (RPE-fus) in 10 species. A phylogenetic reconstruction shows that both RPE-tuf and RPE-fus have evolved in a manner that is consistent with the expected species divergence. The analysis provides evidence for independent loss of RPE-tuf in several species, possibly mediated by short repetitive sequences flanking the site of excision. The remaining RPE-tuf sequences evolve as neutral sequences in different stages of deterioration. Likewise, highly fragmented remnants of the RPE-fus sequence were identified in two species. This suggests that genome-specific differences in the content of RPEs are the result of recent loss rather than recent proliferation.

Molecular characterization of long direct repeat (LDR) sequences expressing a stable mRNA encoding for a 35-amino-acid cell-killing peptide and a cis-encoded small antisense RNA in Escherichia coli

Genome sequence analyses of Escherichia coli K-12 revealed four copies of long repetitive elements. These sequences are designated as long direct repeat (LDR) sequences. Three of the repeats (LDR-A, -B, -C), each approximately 500 bp in length, are located as tandem repeats at 27.4 min on the genetic map. Another copy (LDR-D), 450 bp in length and nearly identical to LDR-A, -B and -C, is located at 79.7 min, a position that is directly opposite the position of LDR-A, -B and -C. In this study, we demonstrate that LDR-D encodes a 35-amino-acid peptide, LdrD, the overexpression of which causes rapid cell killing and nucleoid condensation of the host cell. Northern blot and primer extension analysis showed constitutive transcription of a stable mRNA (approximately 370 nucleotides) encoding LdrD and an unstable cis-encoded antisense RNA (approximately 60 nucleotides), which functions as a trans-acting regulator of ldrD translation. We propose that LDR encodes a toxin-antitoxin module. LDR-homologous sequences are not pre-sent on any known plasmids but are conserved in Salmonella and other enterobacterial species.