Previous Article | Next Article ![]()
Infection and Immunity, March 1999, p. 1116-1124, Vol. 67, No. 3
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
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.
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.
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
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
![]()
ABSTRACT
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
![]()
INTRODUCTION
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
![]()
MATERIALS AND METHODS
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
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).
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 |
|---|
|
|
|---|
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).
|
|
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.
|
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
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.
|
|
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 |
|---|
|
|
|---|
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.
|
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.
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, CTX
(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
CTX
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 |
|---|
|
|
|---|
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 47. |
Staden, R.
1982.
An interactive graphics program for comparing and aligning nucleic acid and amino acid sequences.
Nucleic Acids Res.
10:2951-2961 |
| 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]. |
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»