Previous Article | Next Article ![]()
Journal of Bacteriology, April 2003, p. 2131-2142, Vol. 185, No. 7
0021-9193/03/$08.00+0 DOI: 10.1128/JB.185.7.2131-2142.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Departments of Microbiology and Infectious Diseases,1 Biological Sciences,5 Biochemistry and Molecular Biology, University of Calgary, Calgary, T2N 4N1 Alberta, Canada,6 Cancer Research Center and University of Missouri, Columbia, Missouri 65201,3 Department of Microbiology, Harbin Medical University, Harbin 150086 ,4 Department of Microbiology, Peking University School of Basic Medical Sciences, Beijing 100083, China2
Received 6 May 2002/ Accepted 7 January 2003
|
|
|---|
|
|
|---|
Genomic diversification does occur among cells of the same clone in culture and presumably also in nature, including duplication and inversion of genomic segments, at frequencies as high as 10-5 to 10-3 (3, 13). However, most of these types of genomic changes are selected against in nature and therefore cannot be easily detected. In order to capture such genomic changes for the elucidation of the basic rules in genomic diversification during evolution, we examined serovar Typhimurium cultures that had been stocked for several decades in sealed agar stabs at room temperature, a condition allowing slow growth for some time until the nutrients become exhausted. This archival system of serovar Typhimurium cultures provides a unique chance to detect genetic and genomic changes. Assuming that the originally uniform medium would become distinct in different parts of the stab with the growth of the bacteria and their competition with one another, we would expect to see survivors with unique sets of spontaneous genomic changes selected in different parts of the stab in adaptation to the local specific ecological conditions. We aim at detecting the hypothesized genomic changes and documenting the genomic diversification in this system. We have a collection of spontaneous mutants of serovar Typhimurium (43), including those derived from strain LT2 (16) and those from strain LT7 (15, 16, 35). When a program of developing a chromosomal map of serovar Typhimurium was initiated a half century ago, it was recognized that auxotrophic mutants selected from the LT2 strain had low rates of reversion to prototrophy. However, auxotrophs selected from LT7 had high reversion rates, making it difficult to distinguish recombinants from revertants upon transduction with phage P22. Therefore, the resulting chromosomal map was constructed with LT2, not LT7.
The LT7 strain was initially designated a phenotypic mutator and later shown to be due to a mut defect (15, 35). This genotype would indicate a possibility that populations within the LT7 culture might have greater diversity than those within the LT2 culture. In a preliminary parallel study with LT2 and LT7 mutants, we found great uniformity with LT2 mutants; however, the archival LT7 mutants (in stabs at room temperature) showed much diversity in both colony morphology and genome structure. In contrast, the wild-type strain of LT7 that had been stocked at -70°C did not show variation in colony morphology or genome structure. Diversification of genome structure of the archived LT7 mutants included inversions, translocations, duplications, and point mutations. Interestingly, these genomic changes obviously occurred during archiving in the past four decades, and no further changes had been observed when the bacteria were grown in rich medium for up to nine successive single-colony isolations. We now report the detected genomic changes of serovar Typhimurium LT7 during the diversification process and discuss their possible roles in evolution.
|
|
|---|
kesander). Tn10 insertion mutants of serovar Typhimurium LT2 were obtained from numerous sources and are described elsewhere (2, 42); they were used for mapping the LT7 genome. Transfer of Tn10 insertions through bacteriophage P22-mediated transduction. A large number of Tn10 insertions into genes with known functions have been mapped on the genome of serovar Typhimurium LT2 (21). We transferred Tn10 insertions from serovar Typhimurium LT2 to serovar Typhimurium LT7 SGSC1417 by bacteriophage P22-mediated transduction to locate the same genes through homologous recombination. We made P22 lysates from a selected set of Tn10 insertion mutants of serovar Typhimurium LT2 by growing a 3-ml overnight culture in Luria-Bertani (LB) broth of these selected Tn10 mutants and inoculating these cultures with phage P22 at a multiplicity of infection of 1:100, followed by coincubation for 6 h. After destruction of living bacterial cells by two drops of chloroform and removal of the cell debris by centrifugation, the lysates, at 1011 PFU/ml, were ready for use in the transduction. For transferring the Tn10 insertions to serovar Typhimurium LT7, we spread 100 µl of an overnight culture of serovar Typhimurium LT7 and 20 µl of lysate onto an LB plate containing tetracycline. A colony was picked up and restreaked on another tetracycline plate for single-colony isolation. One colony from the second tetracycline plate was used for phenotype tests and mapping.
Enzymes and chemicals. I-CeuI, AvrII, and SpeI were purchased from New England BioLabs; XbaI and proteinase K were from Boehringer Mannheim. [32P]dCTP was from New England Nuclear. Most other chemicals were from the Sigma Chemical Co.
PFGE methods and genomic mapping.
Preparation of intact genomic DNA, endonuclease cleavage of DNA in agarose blocks, and separation of the DNA fragments by pulsed-field gel electrophoresis (PFGE) were as described previously (22, 29). PFGE was performed with the Bio-Rad CHEF Mapper or Bio-Rad CHEF DRII electrophoresis system. For PFGE, we normally use three cycles of conditions: the first for general separation at 30 s ramping to 90 s for 16 h at 6 V/cm and buffer temperature of 12°C, the second for expanding areas with closely packed small bands at 3 s ramping to 6 s, and the third for expanding areas with closely packed large bands at pulsing times based on the sizes of the bands. The total run times for the second and third cycles were usually 6 to 12 h, depending on the extent of the separation. Most runs were carried out at 120° angle. For very crowded areas of bands, a 150° angle was used. For determining the sizes of DNA fragments on the PFGE gel, we most often used only
ladder (New England Biolabs) as the size marker, but in many cases we also used bacterial genomic DNA cleaved with an endonuclease as markers. Among the ones we often used was serovar Typhimurium LT2 or serovar Typhi Ty2 DNA cleaved with XbaI, AvrII, or SpeI; sizes of these fragments had been determined previously (21, 24). Genomic mapping methods with I-CeuI were described (19) and further optimized (25). The techniques of double cleavage and end labeling were also described previously (21).
Long-range PCR procedure. For the preparation of templates, genomic DNA was extracted with the Promega Wizard genomic DNA purification kit (catalog no. A1120). We used the following 16 primers: primer 1, GGAACGTTGAAGACGACGAC; primer 2, GCAAGCTGCTTCCTGTTACC; primer 3, TCAATGCTGGAAAAGTCTTGC; primer 4, AACTGGTTCCTGGCAAAGTG; primer 5, CAGGCGCTCAGTAGTTGTTG; primer 6, CCCGTTTTACAGCGTTATGG; primer 7, AACCGAATGCAGGGATAGTC; primer 8, AAATGTCGGGACAAAAGTGC; primer 9, CATTTGGCTGCAAAACAGC; primer 10, CTGGGCGAATTCGATGATAC; primer 11, TTGGGAAACCGTATCCATTG; primer 12, AGACGGACATCGCCAATAAC; primer 13, GCGGTGGTTAGAGAAAGCAC; primer 14, CTATTTCGCCGGGAAGAATC; primer 15, ATGGGTATCATTGCGCTTTC; and primer 16, GCCGATGTCGTTAATCTGC. The locations of all of the primers except primers 1 and 2, which amplify a 0.5-kb segment of genomic DNA as a positive control of the conditions for the experiments, are shown on the map of LT7 (see Fig. 7B). The lower mix included 2 µl of deoxynucleoside triphosphates (10 mM), 4 µM forward primer (5 µl), 4 µM reverse primer (5 µl), and template (10 µl). After the mixture was heated at 90°C for 30 s, the upper mix (10x PCR buffer, 5 µl; 25 mM MgCl2, 3 µl; Tap polymerase, 0.25 µl; double-distilled H2O, 19.75 µl) was added. PCR was carried out for 30 cycles with denaturation at 96°C for 1 min, annealing at 57°C for 1 min, and extension at 72°C for 10 min. The PCR products were analyzed by gel electrophoresis with 1% agarose. We used a procedure similar to that described by Helm and Maloy (12).
![]() View larger version (30K): [in a new window] |
FIG. 7. Long-range PCR to confirm the genomic changes revealed by PFGE. (A) PCR products on an agarose gel, with the templates and primers shown on the right; (B) locations and orientations of the primers shown on the genome of LT7. HindIII-cleaved DNA was used as the DNA fragment size marker; negative controls of PCR products are not shown in this figure.
|
|
|
|---|
We suspected that there might be genomic changes among the bacteria with diverse colony sizes, which may be revealed through genomic analyses. For that, we needed first to make a genome map for LT7 and then locate the possible genomic changes on the map.
Endonuclease cleavages of serovar Typhimurium LT7 genomic DNA and the use of Tn10. The wild-type serovar Typhimurium LT7, SGSC1417, which had been stocked at -70°C and did not show diversification in colony size or genome structure (see below), was used to construct a genome map. Cleavages of genomic DNA of serovar Typhimurium LT7 with XbaI, I-CeuI, AvrII, and SpeI generated 25, 7, 13, and 35 fragments, respectively (Fig. 1). Taking advantage of the XbaI and AvrII cleavage sites within the Tn10 DNA sequence, we located a number of genes on the LT7 genome as landmarks through the Tn10 insertions transferred from LT2. Most of the Tn10 insertions were located at homologous sites in the genome of LT7. In the case of a Tn10 insertion mutant, an original XbaI fragment was usually missing because of the extra XbaI site in the Tn10. In the meantime, two new XbaI fragments would appear, adding up to the size of the missing XbaI fragment plus 9.3 kb, the size of Tn10; see the illustration in Fig. 3 in reference 22. AvrII works in a similar manner, except that it has two cleavage sites, one each in the IS10 sequences of Tn10 at the left and right ends. Therefore, an AvrII fragment will be split into three: two of the bacterial DNA (summing up to the wild-type fragment plus 2 kb, 1 kb each of the IS10 sequence from the end to the AvrII site) and one of the inner part of Tn10 between the two AvrII sites (ca. 7 kb). Unlike XbaI or AvrII, SpeI does not cut the Tn10 sequence. However, the SpeI fragments that have Tn10 insertions can be recognized by size increase (9.3 kb; see Fig. 1, lane 6). This feature is sometimes very advantageous for reliably assigning a Tn10 to a certain genomic location by the increased size of an SpeI fragment, whereas in the case of XbaI or AvrII, a Tn10 inserted a few kilobases to an end of a fragment may be difficult to locate on the genome.
![]() View larger version (49K): [in a new window] |
FIG. 1. PFGE patterns of genomic DNA of the wild serovar Typhimurium LT7 strain stored at -70°C cleaved with the endonucleases XbaI, AvrII, I-CeuI, and SpeI. Lanes: 1, XbaI cleavage of LT7 (3 small fragmentsW [6.5 kb], X [6.4 kb], and Y [1 kb]had run out of the gel); 2, I-CeuI cleavage of LT7; 3 and 4, AvrII cleavages of LT2 and LT7, respectively, for a comparison; two small fragments of LT7 had run out of the gel, including L (4 kb) and M (2 kb); 5, SpeI cleavage of LT7 (12 small fragments, X through II, ranging from 49 to 7 kb, had run out of the gel); 6, SpeI cleavage of LT7 serB::Tn10, showing the size shift of S112 with the Tn10 insertion (it is now 112 + 9, i.e., 121 kb).
|
![]() View larger version (23K): [in a new window] |
FIG. 3. Genomic translocation exemplified by strain 8111B. (A) PFGE gels of LT7 and 8111B. Lanes: 1, LT7 DNA cleaved by XbaI, with the fragments and their sizes indicated on the left of the PFGE gel; 2, 8111B DNA cleaved by XbaI, with deviations of the cleavage pattern from that of LT7 indicated on the right of the PFGE gel; 3, LT7 DNA cleaved by I-CeuI, with the fragments and their sizes indicated on the left of the PFGE gel; and 4, 8111B DNA cleaved by I-CeuI, with deviations of the cleavage pattern from that of LT7 indicated on the right of the PFGE gel. (B) Local comparison of LT7 and 8111B showing the translocation of I-CeuI D, which resulted in three hybrid rrn operons, disappearance of two XbaI fragments (F and E) and two SpeI fragments (K and A [not shown on the PFGE picture]), and appearance of two new fragments each from XbaI (198 and 478 kb) and SpeI (377 and 240 kb) digestions.
|
![]() View larger version (47K): [in a new window] |
FIG. 2. Genome map of serovar Typhimurium LT7. The seven rrn operons are shown inside the map.
|
Duplications. I-CeuI E is duplicated in several lines of the 8117 derivatives, including 8117C, 8117D, 8117E, and 8117F, and other strains (8615A, etc.). I-CeuI F is duplicated in 8618E1. These duplications also obviously occurred through homologous recombination between rrn operons (Fig. 4).
![]() View larger version (28K): [in a new window] |
FIG. 4. Genomic duplication exemplified by strain 8117C. (A) PFGE gels of LT7 and 8117C. Lanes: 1, LT7 DNA cleaved by XbaI; 2, 8117C DNA cleaved by XbaI; 3, LT7 DNA cleaved by I-CeuI; and 4, 8117C DNA cleaved by I-CeuI. (B) Local comparison of LT7 and 8117C showing the duplication of I-CeuI E, which resulted in two hybrid rrn operons, the disappearance of the XbaI fragment F and SpeI fragment K (not shown on the PFGE picture), and the appearance of one new fragment each from XbaI (449 kb) and SpeI (334 kb) digestions.
|
![]() ![]() View larger version (49K): [in a new window] |
FIG. 5. Genomic inversion exemplified by strain 8111D323. (A) PFGE gels of LT7 and 8111D323. Lanes: 1, LT7 DNA cleaved by XbaI; 2, 8111D323 DNA cleaved by XbaI; 3, LT7 DNA cleaved by AvrII; and 4, 8111D323 DNA cleaved by AvrII. (B) Local comparison of LT7 and 8111D323 showing the inversion of I-CeuI A, which resulted in two hybrid rrn operons, the disappearance of XbaI B and AvrII A, and appearance of two new fragments each from XbaI (411 and 322 kb) and AvrII (777 and 757 kb) digestions.
|
![]() View larger version (40K): [in a new window] |
FIG. 6. Point mutations as inferred from the creation or loss of endonuclease cleavage sites. (A) XbaI-cleaved genomic DNA of LT7 (lane 1) and 9059A (lane 2); (B) AvrII-cleaved genomic DNA of LT7 (lane 1) and 8608D (lane 2); (C) I-CeuI-cleaved genomic DNA of LT7 (lane 1) and 8111F231213 (lane 2). For these three strains, genomic translocations, duplications, and inversions had been ruled out.
|
Diversity of the global structure of the genome among the serovar Typhimurium LT7 mutants. Translocations, duplications, inversions, and point mutations all contributed to diversification of genome structure in the serovar Typhimurium LT7 mutants. Although many of the genomic changes reported above seemed to be a single event in a given strain, we could not rule out the possibility that they were the "net" result of multiple changes. For example, the translocations may be the final result of a number of inversions. A special case with strain 8608F2 clearly demonstrates that multiple genomic changes could occur in a given strain (Fig. 8). In 8608F2, multiple translocations resulted in an extraordinary genome structure. Figure 8 also presents some other examples of the diversifying genomes among the serovar Typhimurium LT7 mutants, including the I-CeuI F translocation in 9052F2, the point mutation in rrnD of 8111F23332, and I-CeuI F duplication in 8618E1. Note that in spite of the dramatic genomic changes, the genomic location of oriC changed not much.
![]() View larger version (24K): [in a new window] |
FIG. 8. Diverse genome structures among representative serovar Typhimurium LT7 mutants. (A) PFGE patterns of I-CeuI-cleaved genomic DNA of representative serovar Typhimurium LT7 mutants; (B) I-CeuI maps of these strains based on the data from panel A. In lane 3, the loss of the I-CeuI cleavage site in rrlD resulted in the disappearance of I-CeuI C; the disappearance of I-CeuI B and the appearance of I-CeuI B+C (1306 kb) are not obvious on this picture but have been confirmed in other experiments (data not shown). These maps, although constructed based initially on only I-CeuI data, were all confirmed by long-range PCR. The solid circle in I-CeuI C indicates the position of oriC.
|
|
View this table: [in a new window] |
TABLE 1. Summary of genomic changes on representative strains of S. enterica serovar Typhimurium
|
|
|
|---|
The diversity of genome structure among the archived LT7 strains that we have revealed in the present study is striking, considering the high degree of conservation in genome structure between serovar Typhimurium LT2 and E. coli K-12. Among nearly 200 strains examined, we found translocations, duplications, inversions, and point mutations. We did not find genomic insertions, which could be accounted for by the fact that the archived bacteria were in pure culture and therefore did not have the chance to acquire exogenous DNA. We did not find large-scale (12 kb or larger) deletions either, but we do not rule out the possibility that small deletions that are beyond the resolution of our mapping methods (ca. 12 kb) may have escaped our detection. We also considered the possible roles of IS200 in the genomic diversification. IS200 has six copies in serovar Typhimurium (44) and 25 copies in S. enterica serovar Typhi (40), so it may be a hot spot for homologous recombination. In serovar Typhi, an inversion around the replication terminus is mediated by a pair of IS200 (1). However, none of the rearrangements of the LT7 mutants detected in the present study was mediated by IS200. All copies of IS200 in serovar Typhimurium (6 copies) and in serovar Typhi (25 copies) are located in intergenic regions and are usually fairly stable in the genome, causing no genotypic or phenotypic changes. So we do not expect a major role of IS200 in the genomic diversification of LT7.
Translocations were frequent and always involved one of the three small I-CeuI fragmentsD, E, and Fin the present study. They could be the results of deletion and reinsertion, but it is more likely that they are the final products of multiple homologous recombination events because inversions seem to be more frequent among Salmonella spp. (28). Duplications also involve one of the three small I-CeuI fragments. In the present study, I-CeuI E and F are duplicated in 8117C and 8618E1, respectively; I-CeuI D is seen duplicated in other Salmonella strains (S.-L.L, unpublished data). Duplications may double the amount of certain gene products, which might be needed by the bacteria for adaptation.
Inversions often occur among enteric bacteria, presumably as a compensation mechanism to restore a physical balance of the chromosome when there is a major insertion or deletion on the genome (23, 26). In this work, however, inversion occurred when there was no detectable insertion or deletion. Probably, after archival storage for several decades under nonnatural conditions, a large number of genes have lost functions, but products of some other genes might be in greater demand and therefore have somehow obtained closer proximity to the origin of DNA replication (oriC) by inversion to reach a higher gene dosage.
Point mutations were inferred from the gain or loss of endonuclease cleavage sites. The number of nucleotides used in the cleavage sites by the four endonucleases, including 6 times 25 for XbaI, 6 times 13 for AvrII, 6 times 35 for SpeI, and 26 times 7 for I-CeuI, totals 620 bp, and the genome is 4,800 kb. Therefore, one base change in the cleavage sequences for the four endonucleases represents 7.7-kb base changes in the whole genome, although this number seems to be too large to be true and needs to be verified at least by sample sequencing. More startling, however, is the loss of an I-CeuI site in the 8111F derivatives (8111F and its descendants), including 8111F231213 (Fig. 6). I-CeuI is an intron-encoded endonuclease (10, 31, 32), which cleaves DNA within bacterial rrl genes and thus determines the copy number and genomic distribution of rrl genes (19). The I-CeuI cleavage site is highly conserved in most bacteria; a limited degeneracy is known, such as in Mycobacterium tuberculosis (4), but no recent mutation has been reported. If one base in the I-CeuI cleavage site (26 bp) represents 110 bp in the whole 23S rRNA gene (2.9 kb), frequent mutations may make 23S rRNA sequence of little value in phylogenetic studies, which is obviously not the case. We need to find out which base is changed in the 8111F strains that made one of the I-CeuI sites no longer cleavable and whether this 23S rRNA molecule is still functional with the mutation in the sequence.
One interesting question is where the genomic diversification events may have occurred. Based on our results, we speculate that slow growth in the stab during archiving at room temperature may have selected for some mutations that enhanced survival in the stab, which in turn led to genomic diversification, with the special genome structure that adapts best to the local ecological settings becoming fixed, although we do not completely rule out the possibility that some diversification events might have occurred in the first burst of luxuriant growth (40 years ago) before the period of stationary phase. In fact, it seems likely that cells were not in stationary phase in the stab but were growing very slowly, using as energy the bodies of their siblings. Initial heavy growth generated a large population that died over time and released nutrients that could be used by the surviving cells to support slow growth; ongoing death exceeded growth and the population was dying off very slowly over four decades. If growth (by partial replacement of dead cells) was occurring in the stabs, it would be hard to know whether the high frequency of mutant cells in an old stab reflects an increased mutation rate or differential reproduction of rare mutant types originally introduced into the stab at very low frequency. It is more likely that the observed increase in mutant frequency reflected selection rather than mutagenesis.
What guarantees the stability of the physical structure of the bacterial genome during evolution? Mechanisms responsible for such high degree of genomic conservation are still not fully understood. However, the methyl-directed mismatch repair (MMR) system may be a key factor, as the main genetic difference between the two commonly used serovar Typhimurium strains, LT7 and LT2, is that LT7 is a mutator with defects in mutL but LT2 is not and does not show genomic diversification (7). MMR systems repair DNA replication errors on newly synthesized DNA strands to guarantee the precise replication of the genomic DNA, whereby they inhibit recombination between nonidentical DNA sequences and safeguard genetic stability (11, 37). Genes involved in MMR, including mutH, mutL, mutS, etc., have been identified and characterized for their individual roles (5, 38, 41), with the MMR system in E. coli being the most extensively studied (36, 37). Cells having defects in MMR are mutators and have been reported for elevated spontaneous mutation rates (34, 35, 45). In the present study, LT7 as a mutator provided us with excellent chances to analyze genomic diversification; we are now sequencing genes in the MMR system, including mutL, for representative strains to investigate the relationships between the MMR defects and genomic diversification and to gain further insights into the molecular basis of genomic diversification. Whether or not MMR is involved in genomic diversification, two issues remain to be discussed and further explored. The first is how bacterial genome structure could be conserved if diversification took place all of the time, as in the case of the LT7 archival cultures. Working with stationary-phase cultures of E. coli, Finkel and Kolter found that populations of surviving cells were highly dynamic even after many months of incubation, and the diversification proceeded along different paths, resulting in the coexistence of multiple mutant forms (9). Many lines of evidence from researchers working with natural populations also indicate that genomic changes do occur frequently among bacteria, but most of them are lethal. As a result, only a very small proportion of the bacterial population undergoing genomic changes may survive under the specific natural selection pressure. Most genomic changes may be of small scale and may not drastically change the overall structure of the genome, as demonstrated by the distinct genomic differences and at the same time highly conserved overall genome structure of most of the more than 2,500 Salmonella species (18, 20, 21, 25, 28). Some lineages of bacteria of the same phylogenetic group, however, may have different genome structures, usually when they have large genomic insertions that may have disrupted the original genomic balance (23, 26, 27). Once a new balance is reached after genomic rearrangements, a new genome structure, very different as it may be from the majority of the bacteria in the same phylogenetic group, will be established and conserved (23). Combining our findings in previous and the current studies, we believe that the conservation of physical structure of bacterial chromosomes is guaranteed jointly by the requirement of a physical balance and by a genomic function that prohibits replication errors, such as the MMR system, although we need to experimentally evaluate the roles of MMR in the conservation of the bacterial genome structure.
The second issue is about the significance of genome diversification in bacterial evolution and, more specifically, speciation. In a model of bacterial speciation, the Adopt-Adapt Model (17, 30), we hypothesize that bacteria speciate by acquisition of exogenous DNA (the "adopt" step) and the ensuing adaptive genomic rearrangements (the "adapt" step). In fact, the bacterial speciation process involves a special kind of genomic diversification that ultimately leads to a new bacterial species. Both steps of the speciation process, adopt and adapt, would require the genome to be more tolerant to mismatches, i.e., for accepting the incoming DNA and for rearranging the chromosome. Further work is needed to establish whether or how the MMR system contributes to the genomic diversification events, as reported here, and to bacterial speciation.
The present study was supported by an operating grant from the Medical Research Council of Canada (MRC grant GOP-38106) to S.-L.L., an operating grant from the Natural Sciences and Engineering Research Council of Canada (NSERC; grant 216912-00) to S.-L.L., an operating grant from the Canadian Institutes for Health Research (CIHR; grant MOP-47817) to S.-L.L., an NSERC operating grant and a grant from the National Institutes of Health (AI34829) to K.E.S., a CIHR operating grant to R.N.J., and Cancer Research Center and University of Missouri Research Board grants to A.E. S.L.L. was also supported by the Alberta Cancer Board. W.-Q.L. is a summer student supported by Alberta Heritage Foundation for Medical Research (June-August 2001 and May-August 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»