I-CeuI fragment analysis of the Shigella species: evidence for large-scale chromosome rearrangement in S. dysenteriae and S. flexneri

I-CeuI fragments of four Shigella species were analyzed to investigate their taxonomic distance from Escherichia coli and to collect substantiated evidence of their genetic relatedness because their ribosomal RNA sequences and similarity values of their chromosomal DNA/DNA hybridization had proved their taxonomic identity. I-CeuI digestion of genomic DNAs yielded seven fragments in every species, indicating that all the Shigella species contained seven sets of ribosome RNA operons. To determine the fragment identities, seven genes were selected from each I-CeuI fragment of E. coli strain K-12 and used as hybridization probes. Among the four Shigella species, S. boydii and S. sonnei showed hybridization patterns similar to those observed for E. coli strains; each gene probe hybridized to the I-CeuI fragments with sizes similar to that of the corresponding E. coli fragment. In contrast, S. dysenteriae and S. flexneri showed distinct patterns; rcsF and rbsR genes that located on different I-CeuI fragments in E. coli, fragments D and E, were found to co-locate on a fragment. Further analysis using an additional three genes that located on fragment D in K-12 revealed that some chromosome rearrangements involving the fragments corresponding to fragments D and E of K-12 took place in S. dysenteriae and S. flexneri.     

Ribosomal RNA operon anti-termination. Function of leader and spacer region box B-box A sequences and their conservation in diverse micro-organisms

All Escherichia coli rrn operons show a common motif in which anti-terminator box B-box A sequences occur twice, first in the leader and again in the 16 S-23 S spacer. In this study we have analyzed several aspects of rrn anti-termination by leader and spacer anti-terminator sequences. Using DNA synthesis and a plasmid test system, we incorporated random changes into the leader anti-terminator region and examined these mutations for their ability to read through a strong terminator. We also examined anti-termination by synthetic box A and by rrn spacer region sequences. Information derived from these experiments was used to search the rrn sequences of other micro-organisms for possible anti-termination features. Our principal conclusions were that: (1) box A was sufficient for terminator readthrough; (2) we could show no positive requirement for box B in our test system; (3) many of the negative anti-terminator mutations caused a promoter up-effect in the absence of a terminator; (4) the search of rrn operons from other micro-organisms revealed that anti-terminator-like box B-box A sequences exist in leader and spacer regions of both eubacteria and archaebacteria. The frequent occurrence of this pattern suggested that the E. coli rrn anti-termination motif is widespread in nature and has been conserved in microbial evolution.     

Gene conservation and loss in the mutS-rpoS genomic region of pathogenic Escherichia coli

The extent and nature of DNA polymorphism in the mutS-rpoS region of the Escherichia coli genome were assessed in 21 strains of enteropathogenic E. coli (EPEC) and enterohemorrhagic E. coli (EHEC) and in 6 strains originally isolated from natural populations. The intervening region between mutS and rpoS was amplified by long-range PCR, and the resulting amplicons varied substantially in length (7.8 to 14.2 kb) among pathogenic groups. Restriction maps based on five enzymes and sequence analysis showed that strains of the EPEC 1, EPEC 2, and EHEC 2 groups have a long mutS-rpoS region composed of a approximately 6.0-kb DNA segment found in strain K-12 and a novel DNA segment ( approximately 2.9 kb) located at the 3' end of rpoS. The novel segment contains three genes (yclC, pad1, and slyA) that occur in E. coli O157:H7 and related strains but are not found in K-12 or members of the ECOR group A. Phylogenetic analysis of the common sequences indicates that the long intergenic region is ancestral and at least two separate deletion events gave rise to the shorter regions characteristic of the E. coli O157:H7 and K-12 lineages.     

Chromosome of the enterohemorrhagic Escherichia coli O157:H7; comparative analysis with K-12 MG1655 revealed the acquisition of a large amount of foreign DNAs.

A complete Xba I and Bln I cleavage map was constructed for the chromosome of an enterohemorrhagic Escherichia coli (EHEC) O157:H7 strain isolated from an outbreak in Sakai City, Japan, in 1996. A comparative chromosome analysis with E. coli K-12 strain MG1655 was made. The EHEC chromosome was approximately 5600 kb in length, 1 Mb larger than that of MG1655. Despite the marked difference in chromosome length, the location and direction of seven rRNA operons of the EHEC strain were similar to those for MG1655. Overall organization of genes common in both strains is also highly conserved. Chromosome expansion was observed throughout the EHEC chromosome, albeit in an uneven manner. A large portion of the chromosome enlargement was observed in the region surrounding the replication terminus, particularly in a segment containing the terA locus. Sample sequencing of 3627 random shotgun clones suggested the presence of approximately 1550 kb strain-specific DNAs on the EHEC chromosome, most of which are likely to be of foreign origin.     

Construction and Initial Characterization of Escherichia coli Strains with Few or No Intact Chromosomal rRNA Operons

The Escherichia coli genome carries seven rRNA (rrn) operons, each containing three rRNA genes. The presence of multiple operons has been an obstacle to many studies of rRNA because the effect of mutations in one operon is diluted by the six remaining wild-type copies. To create a tool useful for manipulating rRNA, we sequentially inactivated from one to all seven of these operons with deletions spanning the 16S and 23S rRNA genes. In the final strain, carrying no intact rRNA operon on the chromosome, rRNA molecules were expressed from a multicopy plasmid containing a single rRNA operon (prrn). Characterization of these rrn deletion strains revealed that deletion of two operons was required to observe a reduction in the growth rate and rRNA/protein ratio. When the number of deletions was extended from three to six, the decrease in the growth rate was slightly more than the decrease in the rRNA/protein ratio, suggesting that ribosome efficiency was reduced. This reduction was most pronounced in the Delta7 prrn strain, in which the growth rate, unlike the rRNA/protein ratio, was not completely restored to wild-type levels by a cloned rRNA operon. The decreases in growth rate and rRNA/protein ratio were surprisingly moderate in the rrn deletion strains; the presence of even a single operon on the chromosome was able to produce as much as 56% of wild-type levels of rRNA. We discuss possible applications of these strains in rRNA studies.     

Regulation of the Escherichia coli rrnB P2 promoter

The seven rRNA operons in Escherichia coli each contain two promoters, rrn P1 and rrn P2. Most previous studies have focused on the rrn P1 promoters. Here we report a systematic analysis of the activity and regulation of the rrnB P2 promoter in order to define the intrinsic properties of rrn P2 promoters and to understand better their contributions to rRNA synthesis when they are in their natural setting downstream of rrn P1 promoters. In contrast to the conclusions reached in some previous studies, we find that rrnB P2 is regulated: it displays clear responses to amino acid availability (stringent control), rRNA gene dose (feedback control), and changes in growth rate (growth rate-dependent control). Stringent control of rrnB P2 requires the alarmone ppGpp, but growth rate-dependent control of rrnB P2 does not require ppGpp. The rrnB P2 core promoter sequence (-37 to +7) is sufficient to serve as the target for growth rate-dependent regulation.     

Loss of the spacer loop sequence from the rrnB operon in the Escherichia coli K-12 subline that bears the relA1 mutation.

A polymorphism affecting the spacer region of the rrnB rRNA operon is described. Strains from a major Escherichia coli K-12 subbranch are missing a 106-nucleotide portion of the rrnB 16S-to-23S spacer, and a 20-nucleotide sequence is found in its place. We have called this mutant operon rrnB2. The rrnB2 spacer was most probably derived from either rrnC or rrnE. This alteration of rrnB may have occurred by a recombinational exchange or by gene conversion. In the genealogy of E. coli K-12 strains, the appearance of rrnB2 is associated with the spontaneous occurrence of the first relaxed mutation, but attempts to show a selective relationship between the two mutational events have had negative results. The sequences of the rrnG and rrnC 16S-to-23S spacers have also been determined and their comparisons to the other rrn operons encoding tRNAGlu2 are presented.     

rRNA operon multiplicity in Escherichia coli and the physiological implications of rrn inactivation.

Here we present evidence that only five of the seven rRNA operons present in Escherichia coli are necessary to support near-optimal growth on complex media. Seven rrn operons are necessary, however, for rapid adaptation to nutrient and temperature changes, suggesting it is the ability to adapt quickly to changing environmental conditions that has provided the selective pressure for the persistence of seven rrn operons in E. coli. We have also found that one consequence of rrn operon inactivation is a miscoordination of the concentrations of initiation factor IF3 and ribosomes.     

Tn10-mediated inversions fuse uridine phosphorylase (udp) and rRNA genes of Escherichia coli.

Two strains carrying metE::Tn10 insertions (upstream of the udp gene) were used to isolate mutants of Escherichia coli overexpressing udp. These strains differ in their gene order; one contains an inversion between the rrnD and rrnE rRNA operons. Selection was based on the ability of overexpressed Udp to complement thymine auxotrophy. Chromosomal rearrangements that connect the udp gene and promoters of different rrn operons were obtained by this selection. Seven of 14 independent mutants selected in one of the initial strains contained similar inversions of the metE-rrnD segment of the chromosome (about 12% of its length). Another mutant contained traces of a more complicated event, inversion between rrnB and rrnG operons, which was followed by reinversion of the segment between metE and the hybrid rrnG/B operon. Similar inversions (udp-rrn) in a strain already carrying an rrnE-rrnD inversion flip the chromosomal segment between metE and rrnD/E in the opposite direction. In this case, inversions are also accompanied by duplications of the chromosomal region between the rrnA and hybrid udp-rrnD/E operons. PCR amplification with a set of oligonucleotides from the rrn, Tn5, and met genes was used for more detailed mapping. Amplified fragments of the rearranged chromosomes connecting rrnD sequences and insertion elements were sequenced, and inversion endpoints were established.     

Nucleotide sequence of the mannitol (mtl) operon in Escherichia coli

The nucleotide sequence of the known portions of the mannitol operon in Escherichia coli (mtlOPAD) has been determined. Both the operator-promoter region and the intercistronic region between the mtlA and mtlD genes (encoding the mannitol-specific Enzyme II of the phosphotransferase system and mannitol-1-phosphate dehydrogenase, respectively) show parallels with corresponding regions of the glucitol (gut) operon, but neither the mtlA nor the mtlD gene products show obvious homology with the corresponding gene products of the glucitol operon. Five potential cyclic AMP receptor protein binding sites were identified in the mtlOP region, all showing near identity with the consensus sequence. Four regions of dyad symmetry (four to seven bases in length), serving as potential repressor binding sites, overlap with the potential cyclic AMP receptor protein binding sites. Repetitive extragenic palindromic (REP) sequences, forming stem-loop structures in the intercistronic region between mtlA and mtlD and following the mtlD gene were identified. Probable terminator sequences were not found in any of these three regulatory regions. Mannitol-1-phosphate dehydrogenase exhibits two overlapping, potential NAD+ binding sites near the N-terminus of the protein. Computer techniques were used to analyse the mtlD gene and its product.     

Comparative and functional analysis of the rRNA-operons and their tRNA gene complement in different lactic acid bacteria

The complete genome sequences of the lactic acid bacteria (LAB), Lactobacillus plantarum, Lactococcus lactis, and Lactobacillus johnsonii were used to compare location, sequence, organisation, and regulation of the ribosomal RNA (rrn) operons. All rrn operons of the examined LAB diverge from the origin of replication, which is compatible with their efficient expression. All operons show a common organisation of 5'-16S-23S-5S-3' structure, but differ in the number, location and specificity of the tRNA genes. In the 16S-23S intergenic spacer region, two of the five rrn operons of Lb. plantarum and three of the six of Lb. johnsonii contain tRNA-ala and tRNA-ile genes, while L. lactis has a tRNA-ala gene in all six operons. The number of tRNA genes following the 5S rRNA gene ranges up to 14, 16, and 21 for L. lactis, Lb. johnsonii and Lb. plantarum, respectively. The tRNA gene complements are similar to each other and to those of other bacteria. Micro-heterogeneity was found within the rRNA structural genes and spacer regions of each strain. In the rrn operon promoter regions of Lb. plantarum and L. lactis marked differences were found, while the promoter regions of Lb. johnsonii showed a similar tandem promoter structure in all operons. The rrn promoters of L. lactis show either a single or a tandem promoter structure. All promoters of Lb. plantarum contain two or three -10 and -35 regions, of which either zero to two were followed by an UP-element. The Lb. plantarum rrnA, rrnB, and rrnC promoter regions display similarity to the rrn promoter structure of Esherichia coli. Differences in regulation between the five Lb. plantarum promoters were studied using a low copy promoter-probe plasmid. Taking copy number and growth rate into account, a differential expression over time was shown. Although all five Lb. plantarum rrn promoters are significantly different, this study shows that their activity was very similar under the circumstances tested. An active promoter was also identified within the Lb. plantarum rrnC operon preceding a cluster of 17 tRNA genes.     

Genomic peculiarity of coding sequences and metabolic potential of probiotic Escherichia coli strain Nissle 1917 inferred from raw genome data

Probiotic Escherichia coli strain Nissle 1917 (O6:K5:H1) is a commensal E. coli isolate that has a long tradition in medicine for the treatment of various intestinal disorders in humans. To elucidate the molecular basis of its probiotic nature, we started sequencing the genome of this organism with a whole-genome shotgun approach. A 7.8-fold coverage of the genomic sequence has been generated and is now in the finishing stage. To exploit the genome data as early as possible and to generate hypotheses for functional studies, the unfinished sequencing data were analyzed in this work using a new method [Sun, J., Zeng, A.P., 2004. IdentiCS--identification of coding sequence and in silico reconstruction of the metabolic network directly from unannotated low-coverage bacterial genome sequence. BMC Bioinformatics 5, 112] which is particularly suitable for the prediction of coding sequences (CDSs) from unannotated genome sequence. The CDSs predicted for E. coli Nissle 1917 were compared with those of all five other sequenced E. coli strains (E. coli K-12 MG1655, E. coli K-12 W3110, E. coli CFT073, EHEC O157:H7 EDL933 and EHEC O157:H7 Sakai) published to date. Five thousand one hundred and ninety-two CDSs were predicted for E. coli Nissle 1917, of which 1065 were assigned with enzyme EC numbers. The comparison of all predicted CDSs of E. coli Nissle 1917 to the other E. coli strains revealed 108 CDSs specific for this isolate. They are organized as four big genome islands and many other smaller gene clusters. Based on CDSs with EC numbers for enzymes, the potential metabolic network of Nissle 1917 was reconstructed and compared to those of the other five E. coli strains. Overall, the comparative genomic analysis sheds light on the genomic peculiarity of the probiotic E. coli strain Nissle 1917 and is helpful for designing further functional studies long before the sequencing project is completely finished.