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.     

Evolutionary genetics of the isocitrate dehydrogenase gene (icd) in Escherichia coli and Salmonella enterica.

Sequences of the icd gene, encoding isocitrate dehydrogenase (IDH), were obtained for 33 strains representing the major phylogenetic lineages of Escherichia coli and Salmonella enterica. Evolutionary relationships of the strains based on variation in icd are generally similar to those previously obtained for several other housekeeping and for invasion genes, but the sequences of S. enterica subspecies V strains are unusual in being almost intermediate between those of the other S. enterica subspecies and E. coli. For S. enterica, the ratio of synonymous (silent) to nonsynonymous (replacement) nucleotide substitutions between pairs of strains was larger than comparable values for 12 other housekeeping and invasion genes, reflecting unusually strong purifying selection against amino acid replacement in the IDH enzyme. All amino acids involved in the catalytic activity and conformational changes of IDH are strictly conserved within and between species. In E. coli, the level of variation at the 3' end of the gene is elevated by the presence in some strains of a 165-bp replacement sequence supplied by the integration of either lambdoid phage 21 or defective prophage element e14. The 72 members of the E. coli Reference Collection (ECOR) and five additional E. coli strains were surveyed for the presence of phage 21 (as prophage) by PCR amplification of a phage 21-specific fragment in and adjacent to the host icd, and the sequence of the phage 21 segment extending from the 3' end of icd through the integrase gene (int) was determined in nine strains of E. coli. Phage 21 was found in 39% of E. coli strains, and its distribution among the ECOR strains is nonrandom. In two ECOR strains, the phage 21 int gene is interrupted by a 1,313-bp insertion element that has 99.3% nucleotide sequence identity with IS3411 of E. coli. The phylogenetic relationships of phage 21 strains derived from sequences of two different genomic regions were strongly incongruent, providing evidence of frequent recombination.     

Nucleotide sequence of aceK, the gene encoding isocitrate dehydrogenase kinase/phosphatase.

In Escherichia coli, the phosphorylation and dephosphorylation of isocitrate dehydrogenase (IDH) are catalyzed by a bifunctional protein kinase/phosphatase. We have determined the nucleotide sequence of aceK, the gene encoding IDH kinase/phosphatase. This gene consists of a single open reading frame of 1,734 base pairs preceded by a Shine-Dalgarno ribosome-binding site. Examination of the deduced amino acid sequence of IDH kinase/phosphatase revealed sequences which are similar to the consensus sequence for ATP-binding sites. This protein did not, however, exhibit the extensive sequence homologies which are typical of other protein kinases. Multiple copies of the REP family of repetitive extragenic elements were found within the intergenic region between aceA (encoding isocitrate lyase) and aceK. These elements have the potential for combining to form an exceptionally stable stem-loop structure (delta G = -54 kcal/mol [ca. -226 kJ/mol]) in the mRNA. This structure, which masks the ribosome-binding site and start codon for aceK, may contribute to the downshift in expression observed between aceA and aceK. Another potential stem-loop structure (delta G = -29 kcal/mol [ca. 121 kJ/mol]), unrelated to the REP sequences, was found within aceK.     

Genetic rearrangements in the tyrB-uvrA region of the enterobacterial chromosome: a potential cause for different class B acid phosphatase regulation in Salmonella enterica and Escherichia coli

Unlike in Escherichia coli, in Salmonella enterica production of class B acid phosphatase (AphA) was detectable also in cells growing in the presence of glucose. Characterization of the aphA locus from a S. enterica ser. typhi strain showed that the aphA determinant is very similar to the E. coli homolog, and that its chromosomal location between the highly conserved tyrB and uvrA genes is retained. However, the aphA flanking regions were found to be markedly different in the two species, either between tyrB and aphA or between aphA and uvrA. The differences in the aphA 5'-flanking region, which in S. enterica is considerably shorter than in E. coli (183 vs. 1121 bp) and includes potential promoter sequences not present in E. coli, could be responsible for the different regulation of class B acid phosphatase observed in the two species.     

Phylogeny and strain typing of Escherichia coli, inferred from variation at mononucleotide repeat loci

Multilocus sequencing of housekeeping genes has been used previously for bacterial strain typing and for inferring evolutionary relationships among strains of Escherichia coli. In this study, we used shorter intergenic sequences that contained simple sequence repeats (SSRs) of repeating mononucleotide motifs (mononucleotide repeats [MNRs]) to infer the phylogeny of pathogenic and commensal E. coli strains. Seven noncoding loci (four MNRs and three non-SSRs) were sequenced in 27 strains, including enterohemorrhagic (six isolates of O157:H7), enteropathogenic, enterotoxigenic, B, and K-12 strains. The four MNRs were also sequenced in 20 representative strains of the E. coli reference (ECOR) collection. Sequence polymorphism was significantly higher at the MNR loci, including the flanking sequences, indicating a higher mutation rate in the sequences flanking the MNR tracts. The four MNR loci were amplifiable by PCR in the standard ECOR A, B1, and D groups, but only one (yaiN) in the B2 group was amplified, which is consistent with previous studies that suggested that B2 is the most ancient group. High sequence compatibility was found between the four MNR loci, indicating that they are in the same clonal frame. The phylogenetic trees that were constructed from the sequence data were in good agreement with those of previous studies that used multilocus enzyme electrophoresis. The results demonstrate that MNR loci are useful for inferring phylogenetic relationships and provide much higher sequence variation than housekeeping genes. Therefore, the use of MNR loci for multilocus sequence typing should prove efficient for clinical diagnostics, epidemiology, and evolutionary study of bacteria.     

Enterobacterial repetitive intergenic consensus (ERIC) sequences in Escherichia coli: Evolution and implications for ERIC-PCR

Enterobacterial repetitive intergenic consensus (ERIC) sequences are 127-bp imperfect palindromes that occur in multiple copies in the genomes of enteric bacteria and vibrios. Here we investigate the distribution of these elements in the complete genome sequences of nine Escherichia coli (including Shigella species) strains. There is a significant tendency for copies to be adjacent to more highly expressed genes. There is considerable variation among strains with respect to the presence of an element in any particular intergenic region, but some copies appear to have been conserved since before the divergence of E. coli and Salmonella enterica. In comparisons of orthologous copies between these species, ERIC sequences are surprisingly conserved, implying that they have acquired some function, perhaps related to mRNA stability. The relationships among copies within E. coli are consistent with a master copy mode of generation. Insertion of new copies seems to occur at, and involve duplication of, the dinucleotide TA. Two classes of inserts of about 70 bp each occur at different specific sites within ERIC sequences; these inserts evolve independently of the ERIC sequences. The small number of ERIC sequences in E. coli genomes indicates that a widely used bacterial fingerprinting method using primers based on ERIC sequences (ERIC-PCR) does not rely on the presence of ERIC sequences.     

Phylogeny and Strain Typing of Escherichia coli, Inferred from Variation at Mononucleotide Repeat Loci

Multilocus sequencing of housekeeping genes has been used previously for bacterial strain typing and for inferring evolutionary relationships among strains of Escherichia coli. In this study, we used shorter intergenic sequences that contained simple sequence repeats (SSRs) of repeating mononucleotide motifs (mononucleotide repeats [MNRs]) to infer the phylogeny of pathogenic and commensal E. coli strains. Seven noncoding loci (four MNRs and three non-SSRs) were sequenced in 27 strains, including enterohemorrhagic (six isolates of O157:H7), enteropathogenic, enterotoxigenic, B, and K-12 strains. The four MNRs were also sequenced in 20 representative strains of the E. coli reference (ECOR) collection. Sequence polymorphism was significantly higher at the MNR loci, including the flanking sequences, indicating a higher mutation rate in the sequences flanking the MNR tracts. The four MNR loci were amplifiable by PCR in the standard ECOR A, B1, and D groups, but only one (yaiN) in the B2 group was amplified, which is consistent with previous studies that suggested that B2 is the most ancient group. High sequence compatibility was found between the four MNR loci, indicating that they are in the same clonal frame. The phylogenetic trees that were constructed from the sequence data were in good agreement with those of previous studies that used multilocus enzyme electrophoresis. The results demonstrate that MNR loci are useful for inferring phylogenetic relationships and provide much higher sequence variation than housekeeping genes. Therefore, the use of MNR loci for multilocus sequence typing should prove efficient for clinical diagnostics, epidemiology, and evolutionary study of bacteria.     

A novel IS element,IS621, of the IS110/IS492 family transposes to a specific site in repetitive extragenic palindromic sequences in Escherichia coli

An Escherichia coli strain, ECOR28, was found to have insertions of an identical sequence (1,279 bp in length) at 10 loci in its genome. This insertion sequence (named IS621) has one large open reading frame encoding a putative protein that is 326 amino acids in length. A computer-aided homology search using the DNA sequence as the query revealed that IS621 was homologous to the piv genes, encoding pilin gene invertase (PIV). A homology search using the amino acid sequence of the putative protein encoded by IS621 as the query revealed that the protein also has partial homology to transposases encoded by the IS110/IS492 family elements, which were known to have partial homology to PIV. This indicates that IS621 belongs to the IS110/IS492 family but is most closely related to the piv genes. In fact, a phylogenetic tree constructed on the basis of amino acid sequences of PIV proteins and transposases revealed that IS621 belongs to the piv gene group, which is distinct from the IS110/IS492 family elements, which form several groups. PIV proteins and transposases encoded by the IS110/IS492 family elements, including IS621, have four acidic amino acid residues, which are conserved at positions in their N-terminal regions. These residues may constitute a tetrad D-E(or D)-D-D motif as the catalytic center. Interestingly, IS621 was inserted at specific sites within repetitive extragenic palindromic (REP) sequences at 10 loci in the ECOR28 genome. IS621 may not recognize the entire REP sequence in transposition, but it recognizes a 15-bp sequence conserved in the REP sequences around the target site. There are several elements belonging to the IS110/IS492 family that also transpose to specific sites in the repeated sequences, as does IS621. IS621 does not have terminal inverted repeats like most of the IS110/IS492 family elements. The terminal sequences of IS621 have homology with the 26-bp inverted repeat sequences of pilin gene inversion sites that are recognized and used for inversion of pilin genes by PIV. This suggests that IS621 initiates transposition through recognition of their terminal regions and cleavage at the ends by a mechanism similar to that used for PIV to promote inversion at the pilin gene inversion sites.     

Nonrandom location of IS1 elements in the genomes of natural isolates of Escherichia coli

We have studied the spatial distribution of IS1 elements in the genomes of natural isolates comprising the ECOR reference collection of Escherichia coli. We find evidence for nonrandomness at three levels. Many pairs of IS1 elements are in much closer proximity (< 10 kb) than can be accounted for by chance. IS1 elements in close proximity were identified by long-range PCR amplification of the genomic sequence between them. Each amplified region was sequenced and its map location determined by database screening of DNA hybridization. Among the ECOR strains with at least two IS1 elements, 54% had one or more pairs of elements separated by < 10 kb. We propose that this type of clustering is a result of "local hopping," in which we assume that a significant proportion of tranposition events leads to the insertion of a daughter IS element in the vicinity of the parental element. A second level of nonrandomness is found in strains with a modest number of IS1 elements that are mapped through the use of inverse PCR to amplify flanking genomic sequences: in these strains, the insertion sites tend to be clustered over a smaller region of chromosome than would be expected by chance. A third level of nonrandomness is observed in the composite distribution of IS elements across strains: among 20 mapped IS1 elements, none were found in the region of 48-77 minutes, a significant gap. One region of the E. coli chromosome, at 98 min, had a cluster of IS1 elements in seven ECOR strains of diverse phylogenetic origin. We deduce from sequence analysis that this pattern of distribution is a result of initial insertion in the most recent common ancestor of these strains and therefore not a hot spot of insertion. Analysis using long- range PCR with primers for IS2 and IS3 also yielded pairs of elements in close proximity, suggesting that these elements may also occasionally transpose by local hopping.      

The nucleotide sequence of the mitochondrial DNA genome of an abundant petite mutant of Saccharomyces cerevisiae carrying the ori1 replication origin

We have determined the 903 bp nucleotide sequence of the mitochondrial DNA genome of a Saccharomyces cerevisiae petite mutant BB5. This petite, containing the 265 nucleotide ori1 region, is representative of a class of petites arising at exceptionally high frequency within the population of spontaneous petites derived from a particular mit- strain Mb12. The DNA sequences of both the ori1 region and the flanking intergenic regions have been compared to those of the corresponding regions of mtDNA in a previously reported petite strain, a1/1R/1 of Bernardi's laboratory, that has a similar (880 bp) repeat unit. The BB5 petite genome carries a canonical ori1 sequence that is identical in both petite mtDNAs, but the flanking intergenic sequences show significant differences between the two petite strains. The divergence is considered to arise from differences in the sequences flanking ori1 in the respective parent strains.     

Integration host factor binds to a unique class of complex repetitive extragenic DNA sequences in Escherichia coli

Interspersed repeated DNA sequences are characteristic features of both prokaryotic and eukaryotic genomes. REP sequences are defined as conserved repetitive extragenic palindromic sequences and are found in Escherichia coli, Salmonella typhimurium and other closely related enteric bacteria. These REP sequences may participate in the folding of the bacterial chromosome. In this work we describe a unique class of 28 conserved complex REP clusters, about 100bp long, in which two inverted REPs are separated by a singular integration host factor (IHF) recognition sequence. We term these sequences RIP (for repetitive IHF-binding palindromic) elements and demonstrate that IHF binds to them specifically. It is estimated that there are about 70 RIP elements in E. coli. Our analysis shows that the RIP elements are evenly distributed around the bacterial chromosome. The possible function of the RIP element is discussed.     

A novel human nonviral retroposon derived from an endogenous retrovirus.

In a human genome, we found dispersed repetitive sequences homologous to part of a human endogenous retrovirus termed HERV-K which resembled mouse mammary tumor virus. For elucidation of their structure and organization, we cloned some of these sequences from a human gene library. The sequence common to the cloned DNA was ca. 630 base-pairs (bp) in length with an A-rich tail at the 3' end and was found to be a SINE (short interspersed repeated sequence) type nonviral retroposon. In this retroposon, the 5' end had multiple copies of a 40 bp direct repeat very rich in GC content and about the next 510 nucleotides were homologous to the 3' long terminal repeat and its upstream flanking region of the HERV-K genome. This retroposon was thus given the name, SINE-R element since most of it derived from a retrovirus. SINE-R elements were present at 4,000 to 5,000 copies per haploid human genome. The nucleotide sequence was ca. 90% homologous among the cloned elements.     

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.     

Joint Distribution of Insertion Elements Is4 and Is 5 in Natural Isolates of ESCHERICHIA COLI

A reference collection of natural isolates of Escherichia coli has been studied in order to determine the distribution, abundance and joint occurrence of DNA insertion elements IS4 and IS5. Among these isolates, 36% were found to contain IS4 and 30% were found to contain IS5. Among strains containing IS4 the mean number of copies per strain was 4.4 +/- 0.8; the comparable figure for IS5 was 3.7 +/- 1.0. Although the presence of the elements among the isolates was independent, among those isolates containing both IS4 and IS5, there was a significant negative correlation in the number of copies of the elements. The reference collection was also studied for the presence of the DNA sequences flanking the single copy of IS4 in the chromosome of E. coli K12. Homologous sequences were found in only 26% of the isolates. The sequences flanking the IS4 invariably occur together, and their presence is significantly correlated with the presence of IS4. In eight of the strains that carry these flanking sequences, an IS4 is located between them, and the sequences are present at the homologous position as in the K12 strain. We suggest that IS4 and its flanking sequences share a common mechanism of dissemination, such as plasmids, and we present evidence that they are included in a much larger transposable element.