IAI FigSearch
Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Byun, R.
Right arrow Articles by Reeves, P. R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Byun, R.
Right arrow Articles by Reeves, P. R.

 Previous Article  |  Next Article 

Infection and Immunity, March 1999, p. 1116-1124, Vol. 67, No. 3
0019-9567/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.

Evolutionary Relationships of Pathogenic Clones of Vibrio cholerae by Sequence Analysis of Four Housekeeping Genes

Roy Byun, Liam D. H. Elbourne, Ruiting Lan, and Peter R. Reeves*

Department of Microbiology, The University of Sydney, Sydney, New South Wales 2006, Australia

Received 17 July 1998/Returned for modification 6 October 1998/Accepted 10 December 1998


    ABSTRACT
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Studies of the Vibrio cholerae population, using molecular typing techniques, have shown the existence of several pathogenic clones, mainly sixth-pandemic, seventh-pandemic, and U.S. Gulf Coast clones. However, the relationship of the pathogenic clones to environmental V. cholerae isolates remains unclear. A previous study to determine the phylogeny of V. cholerae by sequencing the asd (aspartate semialdehyde dehydrogenase) gene of V. cholerae showed that the sixth-pandemic, seventh-pandemic, and U.S. Gulf Coast clones had very different asd sequences which fell into separate lineages in the V. cholerae population. As gene trees drawn from a single gene may not reflect the true topology of the population, we sequenced the mdh (malate dehydrogenase) and hlyA (hemolysin A) genes from representatives of environmental and clinical isolates of V. cholerae and found that the mdh and hlyA sequences from the three pathogenic clones were identical, except for the previously reported 11-bp deletion in hlyA in the sixth-pandemic clone. Identical sequences were obtained, despite average nucleotide differences in the mdh and hlyA genes of 1.52 and 3.25%, respectively, among all the isolates, suggesting that the three pathogenic clones are closely related. To extend these observations, segments of the recA and dnaE genes were sequenced from a selection of the pathogenic isolates, where the sequences were either identical or substantially different between the clones. The results show that the three pathogenic clones are very closely related and that there has been a high level of recombination in their evolution.


    INTRODUCTION
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Vibrio cholerae is a gram-negative bacterium which comprises part of the autochthonous microflora of aquatic environments, often found in close association with a variety of algae and crustaceans (13, 14, 22, 24). Of medical importance, however, is that certain members of the species have evolved mechanisms to become pathogenic to humans, with the potential to cause the severe life-threatening diarrheal disease cholera. A characteristic of the disease is its ability to emerge as explosive outbreaks in human populations. Since epidemiological records of cholera were initiated, the outbreaks have been divided into seven pandemics, with the fifth, sixth, and seventh pandemics caused by strains which carry the O1 antigen (39).

At the end of 1992, a strain of V. cholerae with a novel antigen emerged as a major cause of cholera around the Bay of Bengal in India and Bangladesh (11, 41). Prior to this, non-O1 V. cholerae strains were known to be responsible only for sporadic cases of gastroenteritis and for extraintestinal infections (36). The new form of the antigen was designated O139, and the strain is known as V. cholerae O139 Bengal, after its initial appearance in the Indian subcontinent. V. cholerae O139 Bengal rapidly spread through the immunologically naive populations in neighboring Asian countries. Genetic studies indicate it is closely related to the O1 seventh-pandemic clone and presumably arose by lateral transfer of genes for lipopolysaccharide biosynthesis from an O139 strain to an organism of the seventh-pandemic clone (4, 5).

The techniques traditionally used to assess the relationship between V. cholerae isolates were based mainly on the biochemical characteristics. Recently, various molecular biology-based techniques have been used to study the relationships among clinical and environmental isolates. They include multilocus enzyme electrophoresis (MLEE) (10, 16, 45), pulse-field gel electrophoresis (PFGE) (9), ribotyping (25, 40), and randomly amplified polymorphic DNA (RAPD) (43, 46), which have differentiated isolates of the V. cholerae population into different electrophoretic types (ETs or zymovars), PFGE types, ribotypes, and RAPD fingerprint types, respectively. Application of these molecular epidemiological techniques has shown the existence within the V. cholerae population of several pathogenic clones, primarily isolates considered to be remnants of the sixth pandemic, isolates from the seventh pandemic, and isolates from the U.S. Gulf Coast region of North America.

Another molecular technique which has been used to study the relationships of the pathogenic clones to the predominantly nontoxigenic, environmental isolates of V. cholerae is comparative nucleotide sequence analysis. This technique provides particularly valuable data for population genetic studies, aimed at determining the genetic structures of populations of bacteria and understanding the evolutionary processes that affect rates of nucleotide and amino acid substitutions. Karaolis et al. (26) analyzed the sequence variation in the asd gene from 45 isolates of V. cholerae. No variation was found within the sixth-pandemic, seventh-pandemic, or U.S. Gulf Coast clones, but the asd sequences of the three clones were not closely related.

A single locus may not be representative of a given genome, and as MLEE and other data are discordant with the conclusion from the asd sequences, we sequenced the mdh (malate dehydrogenase [MDH]) gene and a segment of the hlyA (hemolysin) gene from 32 isolates of V. cholerae. In contrast to findings for the asd gene, we found no variation in the mdh gene and hlyA gene (except for the 11-bp deletion in the sixth-pandemic clone) within or between isolates of the pathogenic clones, suggesting that the clones are very closely related. This observation was supported by the limited sequencing of segments of the recA and dnaE genes, which also showed that the level of recombination is high for V. cholerae.


    MATERIALS AND METHODS
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Bacterial isolates. In this study, we examined a total of 33 V. cholerae isolates, comprising 13 O1 clinical isolates of the sixth-pandemic clone (M642, M644, M648, M967, and 569B), along with isolates from pre-seventh-pandemic outbreaks (M543, M640, M645, and M802), from the seventh-pandemic outbreak (M663 and M793), and from outbreaks in the U.S. Gulf Coast region (M794 and M796); 2 O1 nontoxigenic environmental isolates (M535 and M536); 2 O139 Bengal isolates (M539 and M831); and 16 environmental, nontoxigenic, non-O1, non-O139 isolates (M548, M549, M550, M551, M552, M553, M554, M555, M556, M557, M558, M559, M560, M561, M562, and M563). These isolates are from diverse geographical locations and were previously described by Karaolis et al. (26). Strain M967 (#75) is a sixth-pandemic isolate from Japan (1921), and 569B is a remnant of the sixth pandemic from India (1940). Vibrio mimicus M547 was selected for use as an outgroup.

DNA methods. Isolates were stored at -70°C and subcultured onto nutrient agar, from which a single colony was selected and chromosomal DNA was extracted as previously described (3). PCR was performed in reaction mixtures containing 50 mM KCl, 10 mM Tris-HCl (pH 9.0), 1.5 mM MgCl2, bovine serum albumin (200 µg/ml), 100 µM each deoxyribonucleoside triphosphate, 0.3 µM each primer, purified chromosomal DNA (~10 ng/µl), and Taq polymerase (0.02 U/µl). Amplification was performed in an FTS-960 thermal cycler (Corbett Research) with the following program: denaturation at 94°C for 2 min, followed by 35 cycles of 94°C for 15 s, 50 to 60°C for 15 s, and 72°C for 30 s, and a final cycle of 72°C for 10 min. Amplified products were resolved by 1% agarose gel electrophoresis with 0.5× Tris-borate-EDTA as the running buffer and visualized by ethidium bromide staining followed by UV transillumination. PCR primers used in this study are listed in Table 1. PCR amplicons were purified by using the Promega Wizard PCR purification system and sequenced by the dye terminator method at the Sydney University and Prince Alfred Hospital Macromolecular Analysis Centre, using a model 877 integrated thermal cycler and model 377 automated DNA sequencer (Applied Biosystems).

                              
View this table:
[in this window]
[in a new window]
  		TABLE 1.   PCR primers used in this study

Sequencing of the mdh gene of V. cholerae M793. Degenerate primers 502 and 457 were based on highly conserved regions in MDH of closely related species (Photobacterium spp. [accession no. P37226 {52}], Escherichia coli [accession no. M24777 {50}], and Salmonella enterica [accession no. P25077 {31}]) and used in a PCR with chromosomal DNA from strain M793. The expected fragment of 742 bp was excised from a 1% low-melting-temperature agarose gel with 1× Tris-acetate-EDTA as the running buffer, purified by using the Promega Wizard PCR purification system, and ligated into the cloning vector pGEM-T (Promega) for dye-labeled primer sequencing.

Primers 516 and 517, based on the partial mdh sequence of V. cholerae M793, were used in an inverse PCR (IPCR) (21) to amplify the flanking regions. The primers amplified a fragment of approximately 800 bp from the template DNA derived by digestion with PstI, which was cloned and sequenced as described above. Amplification of IPCR fragments from templates obtained from digestions with AccI, StyI, NcoI, SalI, SphI, XhoI, AflII, BclI, BglII, BssHII, EcoRI, SacI, and XbaI was unsuccessful.

Sequencing of the mdh, hlyA, recA, and dnaE genes from selected V. cholerae strains. Primers 540 and 541, based on the flanking sequences of the mdh gene of strain M793, were used to amplify and sequence a 1,039-bp fragment containing the entire 936-bp coding region of the mdh gene from purified chromosomal DNA. These primers failed to amplify from V. mimicus M547 and V. cholerae M552; however, PCR amplicons and partial mdh sequences were obtained from these isolates by using degenerate primers 584 and 585.

Primers 644 and 645, based on the hlyA sequence of V. cholerae 017 (accession no. Y00557 [1]), were used to amplify and sequence a 1,038-bp segment of the 2,226-bp coding sequence of the hlyA gene. Primers 644 and 645 failed to amplify from strains M547 and M552, despite the lower stringency of annealing conditions.

Primers 884 and 885, based on the recA sequence of V. cholerae 017 (accession no. X71969 [49]), were used to amplify and sequence a 1,041-bp segment of the 1,061-bp coding sequence of the recA gene. Primers 713 and 714, based on the dnaE sequence of V. cholerae C6706 (accession no. U30472 [19]), were used to amplify and sequence a 1,067-bp segment of the 3,477-bp coding region of dnaE.

Computer analysis of the sequences. The nucleotide sequences of the mdh, hlyA, recA, and dnaE genes were edited and assembled with the TED (20) and GAP4 (47) programs. Sequences were aligned with the CLUSTALW program, and phylogenetic analysis was performed with PHYLIP (17) and MULTICOMP (42). These programs are accessed through the Australian National Genomic Information Service at the University of Sydney.

Phylogenetic trees were constructed by the neighbor-joining method (44) for the mdh and hlyA genes (Fig. 2). The mdh gene tree was rooted with the partial mdh sequence from V. mimicus M547. The hlyA gene tree was rooted with the vmhA sequence of V. mimicus (accession no. U68271 [28]), which shows 76% nucleotide identity to hlyA from pathogenic isolates and is clearly the same gene with a different name.

Nucleotide sequence accession numbers. The GenBank accession numbers for the nucleotide sequences determined in this study are AF117833 to AF117883.


    RESULTS
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

Sequence of the V. cholerae mdh gene. At the time this study was initiated, no suitable housekeeping genes from V. cholerae were available in the databases. The mdh gene, coding for the metabolic enzyme MDH, was selected because it has been used previously for population studies, which would enable comparison of mdh variation between species (7). Also, the mdh gene trees constructed for E. coli and S. enterica isolates were shown to be congruent with phylogenetic relationships inferred from MLEE.

The mdh gene from strain M793 was sequenced by degenerate PCR and IPCR. The mdh gene shows high levels of similarity to the mdh genes of other bacteria: 72.33% identity to mdhA of a Photobacterium sp. (81.4% amino acid identity; accession no. P37226) and 72.44% identity to mdh of E. coli (79.8% amino acid identity; accession no. M24777). Surprisingly, the V. cholerae mdh gene shows only 70.5% identity to the mdh gene from a psychrophilic Vibrio sp. (37), which shows greater similarity to the mdhA gene of a Photobacterium sp. (74% identity). If the difference in the mdh genes of V. cholerae and the Vibrio sp. isolate reflects the overall genetic differences in their genomes, the taxonomical classification of these isolates requires reassessment.

The mdh gene of V. cholerae is 936 bp in length, coding for 311 amino acids. Compared to other bacterial mdh sequences, the V. cholerae MDH is one amino acid shorter than the 312-residue MDH of other, closely related species (37, 52).

Within the sequence obtained from the 0.8-kb IPCR fragment, a partial open reading frame was identified 327 bp upstream of mdh and in the opposite orientation (data not shown). The open reading frame identified shows 70.7% identity to the argR gene of E. coli (accession no. M17532 [30]). The argR gene is also found upstream of and in the opposite orientation from mdh in the E. coli and Haemophilus influenzae genomes, which suggests that the gene order around the mdh locus has been conserved in these bacteria.

Nucleotide sequence variation in V. cholerae mdh. The nucleotide sequence of the 936-bp coding region of mdh was determined for 32 isolates and a partial mdh sequence obtained from strain M552. Among the 32 complete mdh sequences, 16 variants were identified. Interestingly, 14 of the 15 pathogenic isolates used in this study have identical mdh sequences, the exception being the pre-seventh-pandemic strain M645. Other isolates with identical mdh sequences were environmental isolates M535 and M553, M557 and M558, and M559 and M560. For comparative analysis, the 16 unique mdh sequences were used.

There were 64 polymorphic sites within the mdh sequences analyzed (Fig. 1), with the majority of the substitutions occurring at the 3' end of the gene (codons 151 to 311). Similarly, the polymorphisms detected in mdh of E. coli and S. enterica were mostly within the 3' region of the gene (7). Of the 64 nucleotide substitutions, 57 occurred at the third base of a codon, 1 was at the second base, and 6 were at the first base. Four nonsynonymous substitutions, representing 1.28% of the 311 codons, were found. Of the 64 polymorphic sites, 36 were phylogenetically informative (at least two bases present in two or more of the 16 sequences).


View larger version (46K):
[in this window]
[in a new window]
  		FIG. 1.   Polymorphic sites within the mdh (A) and hlyA (B) genes of V. cholerae isolates. (A) Polymorphic sites within unique mdh sequences. M793 represents the sixth-pandemic, seventh-pandemic, and U.S. Gulf Coast clones, M535 also represents M553, M557, also represents M558, and M559 also represents M560. (B) Polymorphic sites within the 1,038-bp segment of the hlyA gene. M793 represents the seventh-pandemic and U.S. Gulf Coast clones, M642 represents the sixth-pandemic clone, and M555 also represents M559 and M560. Numbering of the polymorphic sites (vertical format) are from the first position of the sequence segment. The position within the codon for each polymorphic site is shown below the sequences. Asterisks indicate polymorphic sites which are informative.

The average pairwise difference for the 16 mdh sequence variants was 1.52%, with a maximum of 4.49% observed between strains M553 and M554 (Table 2). This is similar to the average pairwise difference for 21 unique sequences of asd (1.41%). Average nucleotide differences in mdh variants from E. coli and S. enterica populations were 1.1 and 4.5%, respectively, (7).

                              
View this table:
[in this window]
[in a new window]
  		TABLE 2.   Nucleotide differences between mdh and hlyA genes among isolatesa

The mdh sequences of the pathogenic isolates (except M645) were all identical, suggesting that the pathogenic clones are very closely related. The sequence is most similar to that of the environmental isolate M549, with seven synonymous substitutions and a pairwise difference of 0.75%. The pre-seventh-pandemic isolate M645 differed from the other pathogenic isolates at 1.18% of the nucleotide sites and was most similar to environmental isolates M559 and M560 (0.64% difference).

Partial mdh sequences from V. mimicus M547 and from strain M552 revealed that the average pairwise difference between the typical V. cholerae isolates and M547 is 10.52%, whereas M552 differed from the other V. cholerae isolates at 11.87% and from M547 at 10.45% of the nucleotide sites. The level of divergence of strain M552 from the other V. cholerae isolates and from V. mimicus suggests that it may represent a different species.

hlyA. The hlyA gene codes for the hemolysin which is traditionally used to differentiate between the two biotypes of V. cholerae. It was selected for this study because it is ubiquitous within V. cholerae (8), suggesting that it probably plays a role in the survival of V. cholerae in its natural environment. The loss of hemolytic function among isolates of the sixth pandemic and the gradual loss in isolates from the seventh pandemic would indicate that it is not essential for human pathogenesis (2).

The nucleotide sequence of a 1,038-bp segment of hlyA, extending from codons 12 to 358 and representing 46.6% of the 2,226-bp coding region of hlyA, was determined for 32 V. cholerae isolates. This segment could not be amplified from V. cholerae M552 and V. mimicus M547, even when less stringent conditions were used. From the 32 isolates studied, 17 different hlyA sequences were identified. Similar to the findings for the mdh gene, 14 of the 15 pathogenic isolates have the same hlyA sequence (except for the 11-bp deletion in sixth-pandemic isolates [1]). A different hlyA sequence is also shared by the environmental isolates M555, M559, and M560.

Analysis of the 17 hlyA sequence variants revealed 121 polymorphic nucleotide sites (Fig. 1): 84 at the third base of a codon, 16 at the second base, and 21 at the first base. We found 61 polymorphic sites to be phylogenetically informative, with detection of 24 nonsynonymous substitutions, which represents 6.94% of the 346 residues studied. The average pairwise percent difference within the V. cholerae isolates studied was 3.21%, with a maximum difference of 7.23% observed between isolates M549 and M554 (Table 2). The level of variation in hlyA is more than twice that observed in asd and mdh. The hlyA sequence from the pathogenic isolates is most closely related to the hlyA sequence of environmental O1 isolate M535, with a difference of 1.16% of the nucleotides.

Evidence for recombination. Application of the Stephens test for nonrandom clustering of polymorphic nucleotide sites (48) revealed no detectable cases of intragenic recombination over the 936-bp coding region of mdh. However, a significant partition of 60 bp (bases 354 to 414) (Fig. 1) supported by 16 sites was detected in the hlyA sequences, which separates all of the pathogenic strains (except M645) and the environmental isolates M535, M536, M548, M549, M550, M551, and M553 from the other isolates studied (P < 0.00001). Of the 16 sites, 4 were nonsynonymous, resulting in three amino acid substitutions. As several environmental isolates were affected by the recombination event, which is obscured by subsequent mutation within or proximal to the recombinant region, the recombination event probably occurred significantly before the emergence of the V. cholerae pathogenic clones. Omitting this region, the average pairwise difference between isolates falls to 2.25%, still higher than that for asd or mdh.

Phylogenetic analysis. Phylogenetic trees constructed from the mdh and hlyA sequences (Fig. 2), with and without the regions involved in recombination, show few examples of congruence between the two trees. Low bootstrap values were obtained for most of the nodes in the mdh gene tree, possibly due to recombinational events in mdh which were not detected by the Stephens test. The recombination in hlyA distributed the set of strains into two distinct clusters, which is evident even with the 60-bp segment omitted.


View larger version (17K):
[in this window]
[in a new window]
  		FIG. 2.   Phylogenetic tree for the mdh gene (A) and for the hlyA gene with (B) and without (C) the 60-bp region of recombination. The mdh tree was rooted with the partial mdh sequence of strain M547. The hlyA trees were rooted with the vmhA gene sequence of V. mimicus. Bootstrap values are percentages of 1,000 computer-generated trees and are shown at the nodes. Values of less than 50 are not shown.

Strain M554 is the most divergent V. cholerae strain, as was expected from the pairwise comparisons, although in the mdh tree, strains M557 and M558 cluster with strain M554, whereas in the hlyA tree, they cluster with the other isolates. The clinical pre-seventh-pandemic isolate M645 is found in different clusters than the other pathogenic isolates in both gene trees and has a significant sequence divergence from them in the mdh and hlyA genes, with 1.18 and 3.92% nucleotide differences, respectively.

Nucleotide sequence variation in recA and dnaE. There is no variation within the mdh and hlyA genes for 14 pathogenic isolates of V. cholerae, whereas three distinct sequences are found for asd. To extend these observations, two more genes were selected for sequencing. At the time these experiments were done, sequences of housekeeping genes from biosynthetic pathways more traditionally used for such studies were not available in the databases. Segments of the recA (coding for the RecA protein involved in homologous recombination) and dnaE (coding for the alpha  subunit of the DNA polymerase III holoenzyme) genes were selected because they encode proteins involved in housekeeping roles and therefore are not expected to be under diversifying selection. Sequences were obtained from four isolates of the sixth-pandemic clone (M642, M648, M967, and 569B), four pre-seventh-pandemic outbreak isolates (M645, M802, M543, and M640), two seventh-pandemic isolates (M793 and M663), two U.S. Gulf Coast isolates (M794 and M796), and two non-O1, non-O139 environmental isolates (M549 and M553) which are closely related to the pathogenic isolates in the mdh and hlyA genes.

recA. A 1,041-bp fragment of the recA gene, from positions 25 to 1065 (residues 9 to 354), representing 98.11% of the 1,061-bp coding region, was sequenced for the 14 selected V. cholerae isolates. All the pre-seventh-pandemic (except M645), seventh-pandemic, and U.S. Gulf Coast isolates had the same recA sequence, which is identical to the published recA sequence (accession no. U10162 [32]), in agreement with the differences noted for previously published recA sequences (accession no. X71969 [49] and X61384). The recA sequence of the sixth-pandemic clone differs from that of the seventh-pandemic and U.S. Gulf Coast clones at 48 nucleotide sites, or 4.59% the 1,041-bp segment (Table 3). Of the 48 substitutions, 3 were nonsynonymous. The recA sequence from isolates of the sixth-pandemic clone in this study, one of which is strain 569B, differed from the GenBank sequence of recA from strain 569B (accession no. L42384), which is identical to the recA sequence (U10162) from a U.S. Gulf Coast isolate. We believe that GenBank entry L42384 is recA of a seventh-pandemic strain.

                              
View this table:
[in this window]
[in a new window]
  		TABLE 3.   Nucleotide differences between recA and dnaE genes among isolatesa

The recA sequence difference between the sixth- and seventh-pandemic clones was greater than that observed between the clones in the asd locus. Visual inspection of the polymorphic sites (Fig. 3) within the recA sequences shows that the sixth-pandemic clone and the environmental isolate M549 differ substantially from the other recA sequences between bases 765 and 1011 at 15 sites, indicative of a recombination event. This conclusion is supported by statistical analysis of the data, using the Stephens test, which detected a significant partition, supported by the 15 sites, between (i) strain M549 and the sixth-pandemic isolates and (ii) the other V. cholerae isolates studied (P < 0.00001).


View larger version (35K):
[in this window]
[in a new window]
  		FIG. 3.   Polymorphic sites within the 1,041-bp fragment of recA (A) and the 1,067-bp fragment of dnaE (B) of selected V. cholerae isolates. (A) M793 represents the seventh-pandemic and U.S. Gulf Coast clones, and M642 represents the sixth-pandemic clone. (B) M793 represents the sixth- and seventh-pandemic clones, and M794 represents the U.S. Gulf Coast clone. See the legend to Fig. 1 for further details.

dnaE. A 1,067-bp fragment of the dnaE gene, from positions 1 to 1128 (residues 21 to 376), representing 30.69% of the 3,477-bp coding region, was sequenced for the 14 selected V. cholerae isolates. The sequences revealed that isolates of the sixth- and seventh-pandemic clones and from pre-seventh-pandemic outbreaks (except M645) all have the same dnaE sequence, which differs from that of the U.S. Gulf Coast clone at 21 synonymous sites, which represents 1.97% of the 1,067-bp sequence (Table 3). There was one site of conflict between the sequence from the seventh-pandemic clone and the published dnaE sequence of strain C6706 (accession no. U30472 [19]), at position 1027, which results in an amino acid substitution from Val to Ile at residue 343. However, as the dnaE sequences from 10 isolates in this study were identical, we are confident that our sequence is correct.

The dnaE sequence of the U.S. Gulf Coast clone was most similar to that of the environmental isolate M549, differing at only six nucleotide sites, which represents 0.56% of the region sequenced. Visual inspection of the polymorphic sites (Fig. 3) suggests that the U.S. Gulf Coast clone and M549 may have undergone recombination between bases 135 and 645, an inference supported by application of the Stephens test (P < 0.02594). The dnaE sequence from the sixth- and seventh-pandemic clones is most similar to that of environmental isolate M553, with a pairwise difference of 0.28%.


    DISCUSSION
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

The species V. cholerae shows considerable variation in ETs (16, 45), ribotypes (25, 29, 40), PFGE types (9), RAPD fingerprint types (43, 46), and insertion sequence fingerprint types (6). MLEE studies, which analyze the electrophoretic mobility differences in multiple housekeeping enzymes to determine strain relationships, have differentiated several ETs among pathogenic isolates, being remnants of the sixth-pandemic isolates seventh-pandemic (including O139 Bengal) isolates, isolates from the Latin American outbreak and from outbreaks in the U.S. Gulf Coast region (10), and clinical isolates from Australia (15) and the Amazon rainforest (12). The pathogenic isolates all exhibit similar electrophoretic profiles, differing at only 1 to 3 of the 12 to 16 loci analyzed, implying a closer relationship among these pathogenic ETs than to most of the ETs represented by the predominantly nonpathogenic, environmental isolates, which display diverse electrophoretic patterns and serotypes.

One method to confirm and clarify the phylogenies inferred from the MLEE studies is comparative nucleotide sequence analysis of housekeeping genes, which detects synonymous as well as nonsynonymous substitutions in chromosomally located genes for pairwise comparisons and construction of phylogenetic trees.

Nucleotide sequence variation in mdh, hlyA, recA, and dnaE. The mdh gene and a segment of the hlyA gene were analyzed for 32 V. cholerae isolates. Identical mdh and hylA sequences were found in isolates of the sixth-pandemic clone (except for an 11-bp deletion in hlyA) and isolates from pre-seventh-pandemic (except M645), seventh pandemic, O139 Bengal, and U.S. Gulf Coast outbreaks. The mdh and hlyA sequences from the pathogenic isolates were distinct from the sequences obtained from the environmental nontoxigenic, non-O1, non-O139 isolates, where the average pairwise nucleotide differences between variants were 1.52 and 3.21%, respectively. The pathogenic isolates (except M645) generally have identical recA and dnaE genes, although for recA, the sixth-pandemic clone differed substantially from the other pathogenic isolates, while for dnaE, the U.S. Gulf Coast clone had a sequence which differed substantially from that of the other pathogenic strains.

The frequent occurrence of identical housekeeping genes among the pathogenic isolates suggests that the pathogenic clones are indeed closely related, despite three asd sequence variants being present among these clones, as it seems highly unlikely that identical mdh, hlyA, recA, and dnaE genes could have been transferred into independent lineages, even by means such as hitchhiking with genes acquired as part of their adaptation to pathogenesis.

The pre-seventh-pandemic El Tor outbreak strains M543, M640, and M802, but not M645, can be considered at this stage to be precursors to the seventh-pandemic isolates, as for the five genes studied there are no differences between them and the seventh-pandemic isolates. Strain M645 appears to be an anomalous clinical isolate, as it differed from the other pathogenic isolates in all of the five genes studied and was located within different clusters than the other pathogenic isolates in the mdh and hlyA gene trees.

For each of the genes asd, mdh, hlyA, recA, and dnaE, there is no variation within the sixth-pandemic, seventh-pandemic, or U.S. Gulf Coast clones. These clones are clearly closely related, as for each pair they were identical at two or three of the five genes analyzed (Fig. 4). The significant nucleotide differences at the other loci must be due to recombination, as it seems inconceivable that mutation alone could give such levels of divergence while other genes did not diverge at all. The genes all encode proteins involved in housekeeping functions, with no reason to expect differences in the level of selection to account for the disparity in sequence variation in the different genes. Only the 11-bp deletion in the hlyA gene of the sixth-pandemic clone is attributed to mutation, and this mutation in the sixth-pandemic clone may well have been established by selection, as hlyA-negative forms appeared soon after the major expansion of the seventh-pandemic clone. Thus, among the 4,082-bp sequence in the four genes, there are no differences in the three pathogenic clones attributed to random genetic drift of neutral mutation. This finding indicates a much closer relationship than can be inferred from MLEE data.


View larger version (15K):
[in this window]
[in a new window]
  		FIG. 4.   Representation of the genetic identities of four housekeeping genes between the pathogenic clones of V. cholerae. Genes within the triangle are identical to the adjacent pathogenic clones. Genes outside the triangle are different between adjacent pathogenic clones; percent nucleotide differences are in parentheses.

MLEE studies show that the seventh-pandemic and U.S. Gulf Coast clones differ in the leucine aminopeptidase, DA1 (NADPH diaphorase), and NSE (carboxylesterase) loci, and the sixth- and seventh-pandemic clones differ in the 6-phosphogluconate dehydrogenase and glucose phosphate isomerase loci (16, 45). Whether the differences are due to recombination or arise by mutation cannot be determined by MLEE, but they are very obvious from the sequences. In light of our observations and conclusions from the sequence data, the difference in mobilities of the enzymes between the clones is most likely due to recombination, but this can be confirmed only by sequence analysis of the genes.

Recombination in V. cholerae. The recombination discussed above is not expected to be due to diversifying selection, suggesting a high level of recombination for V. cholerae. We applied the algorithm used by Maynard Smith et al. to detect levels of association between alleles at different loci (33) to the MLEE data set for 260 V. cholerae strains (45). The index of association (IA) is a generalized measure of linkage disequilibrium, with an expected value of zero if the association between loci is random. For all V. cholerae isolates, IA equals 1.57 with a standard error of 0.09, indicating there is a nonrandom distribution of alleles, which is evidence for a clonal population structure. However, this value is biased by the overrepresentation of isolates of the sixth and seventh pandemics of cholera, and when only ETs are considered, IA falls to -0.092 ± 0.17, a value consistent with a nonclonal or weakly clonal population.

A comparison of the phylogenetic trees for mdh, hlyA, and asd shows a lack of congruence between the trees, which is concordant with the conclusion from the statistical test on the isozyme data. The only exception to the lack of congruence are two pairs of strains, M559-M560 and M557-M558, each pair possessing either identical or very similar mdh, hlyA, and asd sequences.

The mdh and hlyA gene trees are more congruent with each other than with the asd gene tree, where the difference in phylogenies is due to recombination. For example, strain M555 has the same hlyA sequence as strains M559 and M560 and a similar mdh sequence (0.32% pairwise difference) but differs in its asd sequence from strains M559 and M560 by 4.16 and 4.06% of the nucleotides, respectively. In addition, strains M535 and M553 have identical mdh sequences and similar hlyA sequences (0.77% pairwise difference) but differ in the asd locus by 2.13%.

The asd locus has been affected by gene transfer among the pathogenic clones also, since these clonal lineages diverge significantly at this locus (26). Across the mdh and hlyA gene trees, the pathogenic clones cluster with the environmental isolates M535 and M548, but in the asd gene tree, it is only the seventh-pandemic clone which clusters with these two environmental isolates. The sixth-pandemic and U.S. Gulf Coast clones are found within different lineages in the asd gene tree, where the noncongruence in the genealogies of these pathogenic clones is most likely a result of independent gene transfers which occurred after their divergence from a common ancestor.

It is interesting that the asd locus has undergone two (or three) recombination events among the pathogenic isolates, whereas there has been only one event involving the whole gene in each of the dnaE and recA genes and no recombination in the mdh and hlyA genes in the pathogenic lineage. This is consistent with a high but random level of recombination involving large (greater than gene size) fragments. However, the fact that there is evidence for intragenic recombination in asd suggests that it is particularly subject to recombination. We have no explanation for the high rate of recombination observed for asd.

The emergence of cholera. Characteristics of the disease cholera, such as long-term immunity of the host, lack of a carrier state, and no known animal host, suggest that it was probably rare or nonexistent in the Paleolithic period, when the relative isolation in which the small hunter-gatherer societies existed could not have supported the continual propagation of such infectious diseases (18). Cholera most probably emerged after the Neolithic revolution, which occurred first in the Middle East some 10,000 years ago, where the invention and/or adoption of agricultural practices by nomadic groups enabled higher densities of humans to subsist. With the establishment of villages and their water supplies, the change from the nomadic to the sedentary lifestyle of human populations provided an opportunity for an environmental V. cholerae bacterium to acquire the necessary virulence mechanisms to survive and multiply in the specific niche of the intestines of humans. Diarrhea induced by extracellular proteins, mainly the cholera toxin, provided the means by which the organism could be released back to the environment, to await infection of the next host.

The genes encoding the cholera toxin comprise part of the genome of the phage encoding cholera toxin, CTXphi (51), which uses the toxin coregulated pilus (TCP) as a receptor, with the genes encoding the TCP being found within a potentially mobile element, the V. cholerae pathogenicity island (VPI) (27). It has been suggested that the adaptation to pathogenesis of V. cholerae involved a sequential process (35), initially requiring the expression of the TCP for CTXphi transduction. For the sixth- and seventh-pandemic clones, whether this process of acquisition occurred before or after their divergence from a common ancestor remains unclear, as the two clones differ in the chromosomal location and copy number of the CTX (cholera toxin) element; the sixth-pandemic clone containing two separate copies compared to the one to three tandem copies of the CTX element found in the seventh-pandemic clone (34). They also exhibit sequence divergence in the ctxB (38) and tcpA (23) genes, which may reflect diversifying selection pressures or indicate independent acquisitions of the CTX and VPI elements.

The emergence of a common ancestor of the present pathogenic clones of V. cholerae probably occurred relatively recently, as no variation was detected within the mdh gene or in the hlyA segment (except for the 11-bp deletion in the sixth-pandemic clone) of the pathogenic isolates over a 57-year time period. Similarly, no variation was detected in the recA and dnaE genes of the pathogenic clones, although recombination in the recA gene of the sixth-pandemic clone, in the dnaE gene of the U.S. Gulf Coast clone, and in the asd gene of these clonal lineages suggests that recombination is frequent in V. cholerae, higher than the mutation rate in these pathogenic clones.

The lack of mutational changes and the high frequency of recombination in the loci studied make it difficult for clear relationships to be determined for the pathogenic clones. From MLEE data, it appears that the U.S. Gulf Coast clone diverged before the emergence of the sixth- and seventh-pandemic clones, suggesting that the U.S. Gulf Coast isolates are remnants from one of the previous pandemics that swept across North America. The 11-bp deletion in hlyA of the sixth-pandemic clone and the fact that some of the characteristics which distinguish the classical and El Tor biotypes involve loss of function (e.g., Vogues-Proskauer reaction and hemagglutination) indicate that it diverged from a common ancestor with the seventh-pandemic clone which had these properties intact. The high rate of recombination and the existence of pathogenic strains from outbreaks between 1937 and 1954 which are very closely related to isolates of the seventh pandemic has major implications for our understanding of how new pandemics emerge. Recombination could be seen as a mechanism whereby recombinant phenotypes are generated from existing pathogenic isolates which, given the right selection pressures, emerge as new pandemics of cholera.


    ACKNOWLEDGMENTS

This project was supported by grants from the Australian Research Council and the National Health and Medical Research Council of Australia.


    FOOTNOTES

* Corresponding author. Mailing address: Department of Microbiology (GO8), The University of Sydney, Sydney, New South Wales 2006, Australia. Phone: (61)(2) 9351 2536. Fax: (61)(2) 9351 4571. E-mail: reeves{at}angis.usyd.edu.au.

Editor:   V. A. Fischetti


    REFERENCES
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References

1. Alm, R. A., U. H. Stroeher, and P. A. Manning. 1988. Extracellular proteins of Vibrio cholerae: nucleotide sequence of the structural gene (hlyA) for the haemolysin of the haemolytic El Tor strain 017 and characterization of the hlyA mutation in the non-haemolytic classical strain 569B. Mol. Microbiol. 2:481-488[Medline].
2. Barrett, T. J., and P. A. Blake. 1981. Epidemiological usefulness of changes in hemolytic activity of Vibrio cholerae biotype El Tor during the seventh pandemic. J. Clin. Microbiol. 13:126-129[Abstract/Free Full Text].
3. Bastin, D. A., L. K. Romana, and P. R. Reeves. 1991. Molecular cloning and expression in Escherichia coli K-12 of the rfb gene cluster determining the O antigen of an E. coli O111 strain. Mol. Microbiol. 5:2223-2231[Medline].
4. Berche, P., C. Poyart, E. Abachin, H. Lelievre, J. Vandepitte, A. Dodin, and J. M. Fournier. 1994. The novel epidemic strain O139 is closely related to the pandemic strain O1 of Vibrio cholerae. J. Infect. Dis. 170:701-704[Medline].
5. Bik, E. M., A. E. Bunschoten, R. D. Gouw, and F. R. Mooi. 1995. Genesis of the novel epidemic Vibrio cholerae O139 strain: evidence for horizontal transfer of genes involved in polysaccharide synthesis. EMBO J. 14:209-216[Medline].
6. Bik, E. M., R. D. Gouw, and F. R. Mooi. 1996. DNA fingerprinting of Vibrio cholerae strains with a novel insertion sequence element: a tool to identify epidemic strains. J. Clin. Microbiol. 34:1453-1461[Abstract].
7. Boyd, E. F., K. Nelson, F. S. Wang, T. S. Whittam, and R. K. Selander. 1994. Molecular genetic basis of allelic polymorphism in malate dehydrogenase (mdh) in natural populations of Escherichia coli and Salmonella enterica. Proc. Natl. Acad. Sci. USA 91:1280-1284[Abstract/Free Full Text].
8. Brown, M. H., and P. A. Manning. 1985. Haemolysin genes of Vibrio cholerae: presence of homologous DNA in non-haemolytic O1 and haemolytic non-O1 strains. FEMS Microbiol. Lett. 30:197-201.
9. Cameron, D. N., F. M. Khambaty, I. K. Wachsmuth, R. V. Tauxe, and T. J. Barrett. 1994. Molecular characterization of Vibrio cholerae O1 strains by pulsed-field gel electrophoresis. J. Clin. Microbiol. 32:1685-1690[Abstract/Free Full Text].
10. Chen, F., G. M. Evins, W. L. Cook, R. Almeida, N. Hargrett-Bean, and K. Wachsmuth. 1991. Genetic diversity among toxigenic and nontoxigenic Vibrio cholerae O1 isolated from the Western Hemisphere. Epidemiol. Infect. 107:225-233[Medline].
11. Cholera Working Group, International Centre for Diarrhoeal Diseases Research, Bangladesh. 1993. Large epidemic of cholera-like disease in Bangladesh caused by Vibrio cholerae O139 synonym Bengal. Lancet 342:387-390[Medline].
12. Coelho, A., J. R. Andrade, A. C. Vicente, and C. A. Salles. 1995. New variant of Vibrio cholerae O1 from clinical isolates in Amazonia. J. Clin. Microbiol. 33:114-118[Abstract].
13. Colwell, R. R., and A. Huq. 1994. Environmental reservoir of Vibrio cholerae. The causative agent of cholera. Ann. N. Y. Acad. Sci. 740:44-54[Medline].
14. Colwell, R. R., and W. M. Spira. 1992. The ecology of Vibrio cholerae, p. 107-127. In D. Barua, and W. B. Greenough III (ed.), Cholera. Plenum Medical Book Co., New York, N.Y.
15. Desmarchelier, P. M., H. Momen, and C. A. Salles. 1988. A zymovar analysis of Vibrio cholerae isolated in Australia. Trans. R. Soc. Trop. Med. Hyg. 82:914-917[Medline].
16. Evins, G. M., D. N. Cameron, J. G. Wells, K. D. Greene, T. Popovic, S. Giono-Cerezo, I. K. Wachsmuth, and R. V. Tauxe. 1995. The emerging diversity of the electrophoretic types of Vibrio cholerae in the Western Hemisphere. J. Infect. Dis. 172:173-179[Medline].
17. Felsenstein, J. 1993. PHYLIP package, version 3.5. University of Washington, Seattle, Wash. (http://evolution.genetics.washington.edu/phylip.html.)
18. Fenner, F. 1970. The effects of changing social organisation on the infectious diseases of man, p. 49-76. In S. V. Boyden (ed.), The impact of civilisation on the biology of man. Australian National University Press, Canberra, Australia.
19. Franco, A. A., P.-E. Yeh, J. A. Johnson, E. M. Barry, H. Guerra, R. Maurer, and J. G. Morris, Jr. 1996. Cloning and characterization of dnaE, encoding the catalytic subunit of replicative DNA polymerase III, from Vibrio cholerae strain C6706. Gene 175:281-283[Medline].
20. Gleeson, T. J., and R. Staden. 1991. An X Windows and UNIX implementation of our sequence analysis package. Comput. Appl. Biosci. 7:398[Free Full Text].
21. Hartl, D. L., and H. Ochman. 1994. Inverse polymerase chain reaction. Methods Mol. Biol. 31:187-196[Medline].
22. Huq, A., E. B. Small, P. A. West, M. I. Huq, R. Rahman, and R. R. Colwell. 1983. Ecological relationships between Vibrio cholerae and planktonic crustacean copepods. Appl. Environ. Microbiol. 45:275-283[Abstract/Free Full Text].
23. Iredell, J. R., and P. A. Manning. 1994. Biotype-specific tcpA genes in Vibrio cholerae FEMS Microbiol. Lett. 121:47-54[Medline].
24. Kaper, J., H. Lockman, R. R. Colwell, and S. W. Joseph. 1979. Ecology, serology, and enterotoxin production of Vibrio cholerae in Chesapeake Bay. Appl. Environ. Microbiol. 37:91-103[Abstract/Free Full Text].
25. Karaolis, D. K., R. Lan, and P. R. Reeves. 1994. Molecular evolution of the seventh-pandemic clone of Vibrio cholerae and its relationship to other pandemic and epidemic V. cholerae isolates. J. Bacteriol. 176:6199-6206[Abstract/Free Full Text].
26. Karaolis, D. K., R. Lan, and P. R. Reeves. 1995. The sixth and seventh cholera pandemics are due to independent clones separately derived from environmental, nontoxigenic, non-O1 Vibrio cholerae. J. Bacteriol. 177:3191-3198[Abstract/Free Full Text].
27. Karaolis, D. K. R., J. A. Johnson, C. C. Bailey, E. C. Boedeker, J. B. Kaper, and P. R. Reeves. 1998. A Vibrio cholerae pathogenicity island associated with epidemic and pandemic strains. Proc. Natl. Acad. Sci. USA 95:3134-3139[Abstract/Free Full Text].
28. Kim, G. T., J. Y. Lee, S. H. Huh, J. H. Yu, and I. S. Kong. 1997. Nucleotide sequence of the vmhA gene encoding hemolysin from Vibrio mimicus. Biochim. Biophys. Acta 1360:102-104[Medline].
29. Koblavi, S., F. Grimont, and P. A. Grimont. 1990. Clonal diversity of Vibrio cholerae O1 evidenced by rRNA gene restriction patterns. Res. Microbiol. 141:645-657[Medline].
30. Lim, D. B., J. D. Oppenheim, T. Eckhardt, and W. K. Maas. 1987. Nucleotide sequence of the argR gene of Escherichia coli K-12 and isolation of its product, the arginine repressor. Proc. Natl. Acad. Sci. USA 84:6697-6701[Abstract/Free Full Text].
31. Lu, C. D., and A. T. Abdelal. 1993. Complete sequence of the Salmonella typhimurium gene encoding malate dehydrogenase. Gene 123:143-144[Medline].
32. Margraf, R. L., A. I. Roca, and M. M. Cox. 1995. The deduced Vibrio cholerae RecA amino-acid sequence. Gene 152:135-136[Medline].
33. Maynard Smith, J., N. H. Smith, M. O'Rourke, and B. G. Spratt. 1993. How clonal are bacteria? Proc. Natl. Acad. Sci. USA 90:4384-4388[Abstract/Free Full Text].
34. Mekalanos, J. J. 1983. Duplication and amplification of toxin genes in Vibrio cholerae. Cell 35:253-263[Medline].
35. Mekalanos, J. J., E. J. Rubin, and M. K. Waldor. 1997. Cholera---molecular basis for emergence and pathogenesis. FEMS Immunol. Med. Microbiol. 18:241-248[Medline].
36. Morris, J. G., Jr. 1990. Non-O group 1 Vibrio cholerae: a look at the epidemiology of an occasional pathogen. Epidemiol. Rev. 12:179-191[Free Full Text].
37. Ohkuma, M., K. Ohtoko, N. Takada, T. Hamamoto, R. Usami, T. Kudo, and K. Horikoshi. 1996. Characterization of malate dehydrogenase from deep-sea psychrophilic Vibrio sp. strain no. 5710 and cloning of its gene. FEMS Microbiol. Lett. 137:247-252[Medline].
38. Olsvik, O., J. Wahlberg, B. Petterson, M. Uhlen, T. Popovic, I. K. Wachsmuth, and P. I. Fields. 1993. Use of automated sequencing of polymerase chain reaction-generated amplicons to identify three types of cholera toxin subunit B in Vibrio cholerae O1 strains. J. Clin. Microbiol. 31:22-25[Abstract/Free Full Text].
39. Pollitzer, R. 1959. History of the disease, p. 11-50. In R. Pollitzer (ed.), Cholera. World Health Organization, Geneva, Switzerland.
40. Popovic, T., C. Bopp, O. Olsvik, and K. Wachsmuth. 1993. Epidemiologic application of a standardized ribotype scheme for Vibrio cholerae O1. J. Clin. Microbiol. 31:2474-2482[Abstract/Free Full Text].
41. Ramamurthy, T., S. Garg, R. Sharma, S. K. Bhattacharya, G. B. Nair, T. Shimada, T. Takeda, T. Karasawa, H. Kurazano, A. Pal, and Y. Takeda. 1993. Emergence of a novel strain of Vibrio cholerae with epidemic potential in southern and eastern India. Lancet 341:703-704[Medline].
42. Reeves, P. R., L. Farnell, and R. Lan. 1994. MULTICOMP: a program for preparing sequence data for phylogenetic analysis. CABIOS 10:281-284[Abstract/Free Full Text].
43. Rivera, I. G., M. A. Chowdhury, A. Huq, D. Jacobs, M. T. Martins, and R. R. Colwell. 1995. Enterobacterial repetitive intergenic consensus sequences and the PCR to generate fingerprints of genomic DNAs from Vibrio cholerae O1, O139, and non-O1 strains. Appl. Environ. Microbiol. 61:2898-2904[Abstract].
44. Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406-425[Abstract].
45. Salles, C. A., and H. Momen. 1991. Identification of Vibrio cholerae by enzyme electrophoresis. Trans. R. Soc. Trop. Med. Hyg. 85:544-547[Medline].
46. Shangkuan, Y. H., C. M. Tsao, and H. C. Lin. 1997. Comparison of Vibrio cholerae O1 isolates by polymerase chain reaction fingerprinting and ribotyping. J. Med. Microbiol. 46:941-948[Abstract].
47. Staden, R. 1982. An interactive graphics program for comparing and aligning nucleic acid and amino acid sequences. Nucleic Acids Res. 10:2951-2961[Abstract/Free Full Text].
48. Stephens, J. C. 1985. Statistical methods of DNA sequence analysis: detection of intragenic recombination or gene conversion. Mol. Biol. Evol. 2:539-556[Abstract].
49. Stroeher, U. H., A. J. Lech, and P. A. Manning. 1994. Gene sequence of recA+ and construction of recA mutants of Vibrio cholerae. Mol. Gen. Genet. 244:295-302[Medline].
50. Vogel, R. F., K. D. Entian, and D. Mecke. 1987. Cloning and sequence of the mdh structural gene of Escherichia coli coding for malate dehydrogenase. Arch. Microbiol. 149:36-42[Medline].
51. Waldor, M. K., and J. J. Mekalanos. 1996. Lysogenic conversion by a filamentous phage encoding cholera toxin. Science 272:1910-1914[Abstract].
52. Welch, T. J., and D. H. Bartlett. 1997. Cloning, sequencing and overexpression of the gene encoding malate dehydrogenase from the deep-sea bacterium Photobacterium species strain SS9. Biochim. Biophys. Acta 1350:41-46[Medline].

Infection and Immunity, March 1999, p. 1116-1124, Vol. 67, No. 3
0019-9567/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.


This article has been cited by other articles:


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow