Previous Article | Next Article ![]()
Journal of Bacteriology, October 2004, p. 6437-6442, Vol. 186, No. 19
0021-9193/04/$08.00+0 DOI: 10.1128/JB.186.19.6437-6442.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
R. Frederick,2 D. W
och,1,
P. Portier,1,¶ M. Blot,1,|| and J. Adams1,2*
Plasticité et Expression des Génomes Microbiens, CNRS FRE2383, Université Joseph-Fourier, Grenoble, France,1 Department of Molecular, Cellular and Developmental Biology, University of Michigan, Ann Arbor, Michigan2
Received 2 April 2004/ Accepted 29 June 2004
|
|
|---|
5-kb range but was bordered in all but two cases by sequences belonging to the preferred consensus target sequence for IS5, YTAR. The ubiquity of such deletions in independent stab cultures and the increase in their frequency over time argue that they have a selective advantage. It is speculated that the loss of the crl locus is responsible for the selective advantage of the deletions. |
|
|---|
40 h. Since the generation time of E. coli in a rich medium under ideal conditions is
40 min, we must conclude that this organism spends the vast majority of its life cycle under conditions of starvation. It is now known that E. coli cultures entering stationary phase undergo many complex physiological and biochemical changes (17, 28). Evidence is now accumulating that cultures maintained in stationary phase continue to undergo genetic changes. Kolter and colleagues (41, 42) have shown that mutants carrying loss-of-function alleles of both lrp and rpoS predominate in stationary-phase cultures, while Naas et al. (26, 27) have described extensive genetic variations arising from insertion sequence (IS) element transposition in stab cultures of E. coli maintained for up to 30 years. In a series of communications, Eisenstark and colleagues (13, 20, 30, 37) showed that archival cultures of Salmonella have accumulated variations in genome size and gene order as well as variations at the gene level. For this communication, we used comparative gene hybridization to systematically scan for changes in the genome of E. coli K-12 maintained in long-term stab cultures.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bacterial strains used for this study
|
Molecular methods. Standard molecular biology procedures were used (35). The digested total DNAs of the bacteria were hybridized with different digoxigenin-labeled probes. The IS5 and IS30 probes were prepared as previously described (26). IS911 DNA was obtained by PCR amplification with the pOF139 plasmid DNA as a template (32) and with primers IS911-1 (5'-GAAGTGGCACACTGAATTTGG) and IS911-2 (5'-TGAAGTGGTCAACAAAAACTGC). The other DNA probes were obtained by PCR amplification with pairs of internal primers for the following genes: for lcpA, primers LcpA1 (5'-AACGAACTGAACGAAGCGGCG) and LcpA2 (5'-AATCTCCTGAATGCGGTCGGC); for dinB, primers DinB1 (5'-AAATTATGCCCACATCTCACC) and DinB2 (5'-CGCGGTGGCGATTAGATCAGC); for yafO, primers YafO1 (5'-CAACCGGTTGCGGTTCTTTC) and YafO2 (5'-GGTTATGCTTGCGTCCACCG); for pepD, primers PepD1 (5'-CGTTCAGGGAGGTTTCAACC) and PepD2 (5'-CACGAAAGATCCTATCCAGC); for crl, primers Crl1 (5'-GCAACCCGAAGAGCAGATTG) and Crl2 (5'-CCGTTAACTTCACCGGCTCG); for phoE, primers PhoE1 (5'-TGCTTCTGCACGCTTGCCTG) and PhoE2 (5'-CGCCAGTAAAGATGGCGACC); for proB, primers ProB1 (5'-GTTGAACTTGTTCGCGAGTG) and ProB2 (5'-CCGAAAATCCAGCGTTTACGG); and for ykfC, primers YkfC1 (5'-AGCTTGCCACATGGGCAGCC) and YkfC2 (5'-GTACGCCACCGCGGGTTTCC). The nucleotide sequence of the region flanking the proB-IS30 deletion was determined with the purified PCR product obtained with an internal IS5 primer (5'-GACTGAGTCAGCCGAGAAG) in combination with primer PepD1, PepD3 (5'-GCTGGATAGGATCTTTCGTG), or YafO1.
Comparative genomic hybridization. Genomic DNAs were isolated by standard procedures from overnight cultures grown in LB medium (35). After digestion with AluI, the DNAs were labeled by random priming with 33P-labeled nucleotides as described previously. The genomic DNAs were then hybridized to E. coli DNA array high-density nylon filters (Sigma-Genosys) by the procedures recommended by the manufacturer. The filters were then exposed overnight on a phosphorimaging screen. The filters were stripped by standard procedures as recommended by the manufacturer. Each high-density array was probed a maximum of six times. The arrays were analyzed with Array Gauge, v. 1.21, software (Fuji Medical Systems, Stamford, Conn.). Intensity measurements at each spot for each exposure were then obtained and were normalized by use of the following equation: normalized spot intensity = (spot intensity reading) x (10,000/average intensity of the reference spots). For each comparison, a total of four measured spot intensities were used. Modifying the work of Riehle et al. (34), we used a regression model, W1 = mW2 + b, to fit the strain W3110 data to obtain a total blot intensity-corrected replicate-averaged value for each gene by the equation W = (W1 + mW2 + b)/2, where W1 and W2 represent different exposures of one strain. In the same way, intensity-corrected replicate-averaged values were calculated for strains A59 and Y34. A residual was then calculated by fitting the data to a second regression model, y = mW + b, and determining the difference between the predicted and measured spot intensities. The residuals were sorted according to their gene order in E. coli K-12, and moving averages were calculated based on a window size of 20. Confidence intervals were calculated as ±3 standard deviations, determined as previously described (34).
|
|
|---|
Both A59 and Y34 were compared to a culture of the ancestral strain of W3110, which had been stored at 80°C. The comparative genomic hybridization intensities for W3110 and Y34 indicated that two large deletions occurred in strain Y34. Deletion I spans
19 kb and is located between 5.4 and 5.9 min on the E. coli K-12 map, while deletion II spans
23 kb and is located between 72.1 and 72.5 min on the map (Fig. 1). The genes located within these two deletions are listed in Table 2. Comparative genomic hybridization intensities for W3110 and the intermediate culture A59 indicated that A59 does not possess any deletions (data not shown).
![]() View larger version (15K): [in a new window] |
FIG. 1. Comparative genomic hybridization of strain W3110 and strain Y34. Genes are plotted by their coordinates according to the E. coli K-12 genome map (6). Confidence intervals were calculated as described in the text.
|
|
View this table: [in a new window] |
TABLE 2. Deletions detected in strain Y34
|
An analysis of the 22 strains used to inoculate stab Y (including A59 and A83) revealed that none were proline auxotrophs. We therefore concluded that deletion I arose during the 5 years of maintenance of this stab culture. To determine the dynamics of the proline auxotroph phenotype in stab cultures, we tested a total of 360 clones isolated from a series of 12 independent stab cultures maintained for 2, 42, or 183 days. Figure 2 shows that the frequency of proline auxotrophs increased steadily, to an average frequency of 0.22 after 183 days. In addition, we tested a total of 201 clones isolated from two additional mixed stab cultures maintained for 5 years (stab Y) (Table 1) and 4.5 years (stab Z) (Table 1). Among these 201 clones, the frequency of proline auxotrophs reached 0.21 and 0.40 after 4.5 and 5 years of starvation, respectively.
![]() View larger version (8K): [in a new window] |
FIG. 2. Increase in average frequency of proline auxotrophs with age of stab culture. Error bars indicate standard errors.
|
![]() View larger version (19K): [in a new window] |
FIG. 3. The type I deletions can be grouped into the following 12 categories: B52 (D22), H122 (H148), F10, G113 (G92 and G95), Y54 (Y34, Y42, Y50, and Y59), E26 (HB101, W2914, and W2915), Z22 (Z27, Z28, and Z57), Z14 (Z29), Y39 (Y4, Y5, Y15, Y24, Y26, Y35, Y45, Y46, Y49, Y58, and Y60), Y52 (Y51 and Y53), Z20 (Z24), and E41. The deletions are displayed as thick lines. Comparisons of the deletion lengths distinguish variable and common deleted regions. At the opposite end from the IS element, the last open reading frame in the common deleted region is crl. The IS30 element is inserted inside an inactivated single copy of IS911. The distances are given (in kilobases) according to the published E. coli K-12 genome (6).
|
![]() View larger version (49K): [in a new window] |
FIG. 4. Sequences of the flanking regions of type 1 deletions. The sequences at the left boundary for H122, B52, F10, G113, Y54, E26, HB101, Y39, Y52, Z20, and E41 matched the consensus preferred target sequence, YTAR, for IS5. The sequences for Y14 and Z52 differed from this consensus sequence by one base. The end of IS5A is underlined and shown in lowercase, and the target sequences of the IS5 insertion are shown in bold. The nucleotide positions on the chromosome are given according the E. coli K-12 genomic sequence (6).
|
(gpt-proA)62 or proA2 (7). In addition, this strain has been used to identify the function of the crl gene by complementation (4, 29). We therefore characterized this deletion by sequencing the region encompassing the junction. The results (Fig. 4) revealed that the deletion present in HB101 is identical to the E26 deletion. In the HB101 lineage (5), strain W2915 was the first ancestral clone to be described as a proline auxotroph, but it has been suggested (M. Berlyn, personal communication) that this phenotype was already present in W2914, the direct ancestor of W2915. We verified that both of these strains indeed carry the same deletion. |
|
|---|
Porwillok et al. (30) have described archival strains of Salmonella enterica that have apparently suffered deletion as well as the amplification of genetic material. Although not common, there have been earlier reports of selection for deletions arising in long-term culture. Riehle et al. (34) used comparative genomic hybridization and observed two deletions in six lines of E. coli B evolving for
2,000 generations under high-temperature stress. Modi et al. (24) described the selection of a plasmid deletion arising in a long-term glucose-limited chemostat culture under conditions in which the so-called hunger response and the starvation response overlapped (14). Similarly, Adams et al. (1) and Dunham et al. (12) observed gross chromosomal changes and rearrangements in evolving yeast populations, some of which were large deletions. The evolutionary significance of deletions and gene loss has been revealed by comparative evolutionary studies, in which the inactivation of nonfunctional or inessential genes often precedes gene loss. This pattern is seen most dramatically in the evolution of parasitic or symbiotic prokaryotic genomes, which show a significant reduction in size (23, 25).
Involvement of IS5 in the generation of type I deletions. The location of an IS5A sequence at the right-hand boundary of all of the type I deletions seen points to the involvement of this element in the generation of the deletions. First described in the case of IS1 (33), several other IS elements have been shown to generate deletions that terminate precisely at one end of the IS element. The mechanism of the formation of such deletions and the transposition of the element appear to share a common pathway (15). Deletions are thought to be generated by a two-step process in which the insertion of a new IS copy near a preexisting IS is followed by homologous recombination between the two copies. IS5 transposes nonrandomly, and a preferred consensus target sequence, YTAR, has been described (21). All but two of the deletions possessed the sequence TTAR, belonging to the preferred consensus sequence YTAR, at the left boundary. The two deletions that did not possess the TTAR sequence at the left boundary had a four-base sequence at the left boundary that differed by only one base from the YTAR consensus. It is therefore hard to escape the conclusion that the IS5 transposition machinery is involved in the generation of these deletions.
Many reports have shown that repeated elements in the genome can often be recombinogenic, facilitating the formation of both duplications and deletions (3, 10, 11, 16, 18, 19). In a recent report, Dunham et al. (12) showed that deletions and major chromosomal rearrangements arising in Saccharomyces cerevisiae were all bordered by sigma, tau, or delta repeated elements. Gene duplications, by providing regions of adjacent homology, are intrinsically unstable and can result in frequent amplifications and deletions of duplicated regions via unequal recombination events (2).
Selective advantage of type I deletions.
Four lines of evidence point to the selective advantage of type I deletions in stab cultures. (i) A type I deletion was observed both in Y34, a clone isolated from a stab culture that was maintained for 5 years, and in several other stab cultures. (ii) The frequency of type I deletions increased steadily over time in short-term stab cultures (Fig. 2). (iii) The variety of different type I deletions seen, which were often characteristic of the particular stab culture, indicates that such deletions arose repeatedly during the evolution of the cultures. Thus, their ubiquity cannot be explained by the occurrence of a unique event occurring in a clone containing a selectively favored mutation ("hitchhiking"). (iv) A type I deletion is present in the widely used laboratory strain of E. coli K-12, HB101, as well as in other strains of the HB101 lineage, such as W2914 and W2915. In addition, a deletion in the same region occurs in strain
711 (8).
Deletion of the crl locus. It is noteworthy that the right boundaries of type I deletions vary over a 5-kb range but that in all 12 categories, the crl gene was deleted and appeared to define the limit of the smallest deletion. We therefore speculate that the selective advantage of type I deletions resides in the loss of the crl gene. In all but two cases, the left boundary of the deletion was characterized by the existence of the YTAR consensus preferred target sequence for IS5 transposition. YTAR sequences are observed frequently in the E. coli genome, occurring every 50 to 100 bases, and can be found adjacent to the left boundary of the deletions. Thus, the absence of sequences of the YTAR family cannot be explained to invoke the lower size limit of type I deletions. Although we cannot rule out mechanical constraints for the lower size limit of type I deletions, the available information for the crl gene implicates it as a likely candidate. First identified as a transcriptional activator of csgAB, the operon encoding the surface protein curli (4, 29), the cytoplasmic Crl factor also modulates the expression of rpoS during starvation (31). Mutations in rpoS that confer a selective advantage during starvation in liquid cultures have already been described (41), suggesting that other mutations modulating the rpoS regulon may be advantageous under starvation conditions.
The dinB locus is maintained by selection.
While the lower size limit of type I deletions can be speculated to be defined by the crl locus, the upper size limit of these deletions appears to be defined by the presence of the dinB locus, which encodes DNA polymerase IV. As with the lower limit defined by the crl locus, the paucity of YTAR sequences in the genome cannot explain this upper limit, as they occur frequently in the E. coli genome, on an average of once every 50 to 100 bases. Similarly, it is unlikely that mechanical constraints define this upper limit, as Brooke and Valvano (8) have described a large IS5-mediated deletion of the crl-dinB-lpcA region in the E. coli K-12 derivative
711. However, it is noteworthy that DNA polymerase IV is involved in SOS-induced untargeted mutagenesis (38, 39) and in the appearance of mutations in cells that are not growing (22). Furthermore, the loss of dinB activity results in a marked reduction in fitness compared with dinB+ strains (40). It is therefore tempting to suggest that the upper size limit of type I deletions is determined by dinB activity and the generation of adaptive genetic variation in stab cultures undergoing low rates of cell turnover. The selective advantage of mutator strains in continuous cultures has likewise been shown to be due to an increase in the rate of generation of genetic variations among mutator strains (9).
In conclusion, the use of comparative genomic hybridization has allowed us to identify large-scale deletions occurring in cultures evolving under long-term conditions of extreme starvation. The analysis of one category of deletion has further allowed us to identify loci that we speculate are both selectively lost and selectively maintained in such cultures.
This work was supported in part by the Centre National de la Recherche Scientifique, the Université Joseph FourierGrenoble 1, a national program of the Bureau des Ressources Génétiques (1999-2000), and NIH grant AI55756.
This paper is dedicated to the memory of Michel Blot, whose contributions to the microbial evolution community will long be remembered. ![]()
Present address: Institut des Sciences du Végétal, CNRS UPR2355, Centre National de la Recherche Scientifique, 91 198 Gif-sur-Yvette, France. ![]()
Present address: Institute of Environmental Sciences, Jagiellonian University, 30-387 Cracow, Poland. ![]()
¶ Present address: Laboratoire d'Ecologie Microbienne, CNRS UMR5557, Université Claude-Bernard Lyon 1, Villeurbanne, France. ![]()
|| Deceased (8 September 2002). ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»