Comparative Genomics of Vibrio cholerae from Haiti, Asia, and Africa - Vol. 17 No. 11 - November 2011 - Emerging Infectious Disease journal - CDC

Volume 17, Number 11—November 2011
THEME ISSUE
CHOLERA IN HAITI

Research

Comparative Genomics of Vibrio cholerae from Haiti, Asia, and Africa

Article Contents

• Methods
• Results
• Discussion
• Acknowledgments
• References
• Figure 1
• Figure 2
• Figure 3
• Table 1
• Table 2
• Table A1
• Technical Appendix 1  [440 KB - 3 pages]
• Technical Appendix 2  [1.89 MB - 3 pages]
• Suggested Citation

Aleisha R. Reimer, Gary Van Domselaar, Steven Stroika, Matthew Walker, Heather Kent, Cheryl Tarr, Deborah Talkington, Lori Rowe, Melissa Olsen-Rasmussen, Michael Frace, Scott Sammons, Georges Anicet Dahourou, Jacques Boncy, Anthony M. Smith, Philip Mabon, Aaron Petkau, Morag Graham, Matthew W. Gilmour, Peter Gerner-Smidt , and the V. cholerae Outbreak Genomics Task Force

Author affiliations: Public Health Agency of Canada, Winnipeg, Manitoba, Canada (A.R. Reimer, G. Van Domselaar, M. Walker, H. Kent, P. Mabon, A. Petkau, M. Graham, M.W. Gilmour); Centers for Disease Control and Prevention, Atlanta, Georgia, USA (S. Stroika, C. Tarr, D. Talkington, L. Rowe, M. Olsen-Rasmussen, M. Frace, S. Sammons, P. Gerner-Smidt); Centers for Disease Control and Prevention, Port-au-Prince, Haiti (G.A. Dahourou); Ministry of Public Health and Population, Port-au-Prince (J. Boncy); National Institute for Communicable Diseases, Johannesburg, South Africa (A.M. Smith)

Suggested citation for this article

Abstract

Cholera was absent from the island of Hispaniola at least a century before an outbreak that began in Haiti in the fall of 2010. Pulsed-field gel electrophoresis (PFGE) analysis of clinical isolates from the Haiti outbreak and recent global travelers returning to the United States showed indistinguishable PFGE fingerprints. To better explore the genetic ancestry of the Haiti outbreak strain, we acquired 23 whole-genome Vibrio cholerae sequences: 9 isolates obtained in Haiti or the Dominican Republic, 12 PFGE pattern-matched isolates linked to Asia or Africa, and 2 nonmatched outliers from the Western Hemisphere. Phylogenies for whole-genome sequences and core genome single-nucleotide polymorphisms showed that the Haiti outbreak strain is genetically related to strains originating in India and Cameroon. However, because no identical genetic match was found among sequenced contemporary isolates, a definitive genetic origin for the outbreak in Haiti remains speculative.

The current (seventh) cholera pandemic was caused by serogroup O1 El Tor biotypes of Vibrio cholerae. This biotype first emerged on the Indonesian island of Sulawesi in 1961, then subsequently spread throughout Asia and Africa, where endemic and epidemic disease persists today (1,2). Seventh cholera pandemic biotypes were introduced into Peru in 1991 and subsequently spread across South and Central America, but these biotypes never reached the island of Hispaniola. Recent endemic and epidemic cases in Asia and Africa are increasingly attributed to genetically atypical El Tor variants that share characteristics of classical and El Tor strains (1,3,4).
After the 2010 earthquake in Haiti, an outbreak of cholera emerged that resulted in >385,000 infections and 5,800 deaths as of July 7, 2011 (5). The outbreak strain quickly spread to the neighboring Dominican Republic and globally as travelers returned home from affected regions (6,7). Concurrent cholera cases in the United States, linked by travel to cholera-endemic regions in Asia and Africa, were identified by national surveillance activities of PulseNet USA (Centers for Disease Control and Prevention [CDC], Atlanta, GA, USA.)
Serotyping, biotyping, and pulsed-field gel electrophoresis (PFGE) fingerprinting investigations suggested that the travel-associated cases could be genetically related to the Haiti outbreak strain (8). Because of the historical absence of cholera in Haiti before the 2010 earthquake, speculation abounds that the outbreak strain was imported into Haiti. Although clonality of the Haiti outbreak strain has been inferred by phenotypic characterization and genotypic subtyping, thereby supporting a single foreign source hypothesis (6,8), definitive evidence, e.g., by whole-genome sequencing for the genetic ancestry of the Haitian strain, is lacking.
Preliminary comparative analysis of whole-genome sequences from two 2010 Haiti outbreak isolates with genomes from historical cholera cases resulted in speculation that the outbreak originated in southern Asia (9). However, this study lacked recent, globally distributed cholera case isolates and particularly lacked studied genomes from Africa, to which cholera is endemic. We selected contemporary V. cholerae isolates from clinical infections, attributed to geographically distinct locations and sharing PFGE fingerprints with Haiti outbreak strains, from the PulseNet USA database for comparative whole-genome analysis. Although detailed epidemiologic investigations are essential for unequivocally attributing geographic origin(s) and means of cholera introduction into Haiti, genome sequences of these 23 contemporary isolates showed details related to genetic content and diversity that were otherwise missed with lower-resolution PFGE subtyping, thereby providing useful genetic ancestry information for interpreting the outbreak in Haiti.

Methods

Patients and Isolates

V. cholerae isolates and travel histories from cholera case-patients in the United States were referred to CDC. A strain from an outbreak in Cameroon in 2010, isolated from a specimen received at CDC, and an isolate from South Africa likely linked to an outbreak in Zimbabwe in 2009 were also included in this study (10). Isolates C6706 and 3569–08 were acquired during the outbreak in Latin America in 1991 and from the US Gulf Coast in 2008, respectively. All strains were characterized as V. cholerae O1 on the basis of standard biochemical, cholera toxin, and serologic testing performed as described (11,12). PFGE was performed according to the PulseNet standardized protocol with restriction enzymes SfiI and NotI; PFGE patterns were designated by using BioNumerics version 5.10 (Applied Maths Inc., Sint-Martins-Latem, Belgium) and compared by unweighted pair group method with arithmetic mean analysis (DICE coefficient 1.5% tolerance and optimization). Strain designations and other information are shown in Table 1.

Whole-Genome Data Acquisition, Assembly, and Annotation

Single-end pyrosequencing reads (GS FLX-Titanium; Roche Diagnostics, Indianapolis, IN, USA) and single-end 36-bp or 76-bp Illumina reads (GAIIe sequencer; Illumina, San Diego, CA, USA) were acquired and yielded >99% genome coverage and 32× and 240× average coverage depths, respectively (Table 2). Pyrosequencing reads were first assembled de novo by using Newbler version 2.5.3 (Roche Diagnostics). To correct potential base-calling errors attributed to homopolymers, Illumina GAIIe reads (average 14 million reads/genome) were mapped to the Newbler contigs by using CLC Genomics Workbench version 4.5 (www.clcbio.com/index.php?id=1042) and yielded an average combined coverage depth of 270× per genome.

Both chromosomes of Haiti outbreak isolate 2010EL-1786 were sequenced to full closure by using PCR and Sanger sequence-based bridging of contigs and a fosmid library of templates. Optical mapping also supported the contig ordering derived for 2010EL-1786. For all remaining isolates, Illumina-supplemented, homopolymer-corrected, Newbler-assembled contigs were prepared as pseudogenomes by first linking contigs with a linker sequence containing stop codons in all 6 translation reading frames. These high-coverage pseudogenomes were used for downstream analyses. Identification of coding sequences was achieved by using Glimmer3 (14). Genome annotation was achieved by using an automated, in-house, modified version of GenDB version 2.2 (15) and manual curation for regions of interest.

Comparative Genomics

Whole-Genome Alignment and Core Genome Phylogeny

Whole-genome alignments of all study isolates and 5 available reference V. cholerae genomes (Table 1) were performed by using Progressive Mauve (16) and visualized by using PhyML 3.0 (17). To determine vertical inheritance patterns, study genomes were analyzed with historical V. cholerae genomes (isolates M66–2, MJ-1236, CIRS101, and N16961) by using phylogenetic analysis of high-quality single-nucleotide polymorphisms (hqSNPs) contained in core genes. Coding region predictions were analyzed by using parallelized BLASTn (http://blast.ncbi.nlm.gov/Blast.cgi) to identify highly similar orthologs in all strains. Highly similar orthologs were defined as those containing a high-scoring segment pair >400 bp and identity >97%. Each orthologous loci set was multiply aligned by using ClustalW (18). Multiple alignments were manually inspected to remove erroneously aligned regions; indel-associated SNPs and loci containing >30 SNPs were also excluded. Each SNP column from each multiple nucleotide alignment was analyzed for hqSNPs, defined as those containing no gaps or ambiguous basecalls, and having an adjusted quality score >90 (of a maximum score of 93). A total of 4,376 hqSNPs were identified from 632 orthologous loci and extracted from the alignments to prepare a compressed pseudoalignment composed of hqSNPs (Technical Appendix 1 ). This pseudoalignment was used to build a maximum-likelihood phylogenetic tree by using PhyML 3.0 (17). Branch confidences were estimated by using the approximate likelihood-ratio test (19).

BLAST Atlases

A circular BLAST atlas was generated for each chromosome by using Haiti isolate 2010EL-1786 as mapping reference. Glimmer3 was used to predict coding sequences contained on pseudogenomes for the remaining isolates sequenced in this study and for 4 available genomes (14). Reference isolate 2010EL-1786 was mapped against the resulting translated coding sequences by using BLASTx with a percentage identity cutoff value of 70% and an expected cutoff value of 1 × 10^–10 for high-scoring segment pairs >100 aa. The results were visualized by using GView (20). Sequence accession numbers are shown in Table 1.

Results

SfiI and NotI PFGE Patterns of Recent Global Cholera Isolates

Nine V. cholerae isolates directly associated with the outbreak on Hispaniola were examined, 7 of which had indistinguishable SfiI and NotI PFGE patterns designated PulseNet USA patterns KZGS12.0088 and KZGN11.0092, respectively (Table 1). Also sequenced were a hemolytic variant and a nonhemolytic variant that harbored a minor variation of the main Haiti outbreak PFGE pattern and were derived from an isolate from 1 patient in Haiti (Table 1). Twelve contemporary V. cholerae isolates from global sources with matched PFGE fingerprints were also sequenced. Infections for these 12 contemporary isolates originated (by documented patient travel) from regions of Pakistan, India, or Nepal. Two additional isolates were from patients in outbreaks in Cameroon and South Africa likely connected to the cholera outbreak in Zimbabwe in 2009 (21). Although all sequenced clinical isolates were serogroup O1, Inaba and Ogawa serotypes were observed among PFGE pattern-matched isolates (Table 1). All strains were biotype El Tor and all produced cholera toxin.

Phylogenetics of Strains

Figure 1

Figure 1. Whole-genome neighbor-joining tree of Vibrio cholerae isolate from cholera outbreak in Haiti, fall 2010; concurrent clinical isolates with pulsed-field gel electrophoresis pattern-matched combinations; reference isolates sequenced in this study; and available...

Haiti outbreak isolates and 12 global PFGE pattern-matched V. cholerae isolates belong to phylogroup 1 of the seventh pandemic clade. The phylogenetic tree based on whole-genome sequencing showed clustering of the 9 Hispaniola isolates (8 from Haiti and 1 related isolate from the Dominican Republic) with 12 other PFGE pattern-matched isolates. All 21 isolates were in 1 cluster relative to non-PFGE–pattern-matched outliers (Figure 1). When compared with historical reference genomes, the closest ancestors for Haiti genome sequences (2010–2011; derived herein) were isolates CIRS101 from Dhaka, Bangladesh (2002) and MJ-1236 from Matlab, Bangladesh (1994). These data confirm the genetic relatedness also inferred by PFGE subtyping and further support inclusion of the Haiti outbreak isolates in phylogroup 1 of the seventh pandemic clade (Figure 1). The whole-genome sequencing dataset showed that additional underlying genetic diversity was present across PFGE pattern-matched isolates (including 9 isolated from Hispaniola) not observed by PFGE subtyping.

Common Mobile Elements and Genes of Haiti Outbreak Strain and PFGE Pattern-matched Isolates

Figure 2

Figure 2. BLAST (http://blast.ncbi.nlm.gov/Blast.cgi) atlas of predicted protein homologies mapped against the closed genome of Haiti Vibrio cholerae outbreak type isolate 2010EL-1786, fall 2010. Full color saturation represents 100% sequence homology, and...

V. cholerae macrodiversity is commonly attributed to presence or absence of mobile genetic elements (22). The contiguous genome derived for Haiti isolate 2010EL-1786 was used as the outbreak type strain and harbored 2 circular chromosomes of 3.03 Mbp (chromosome I) and 1.05 Mbp (chromosome II), which encoded 2,920 and 1,051 predicted coding sequences, respectively. Pairwise comparisons of all coding sequences from each study genome with all coding sequences from reference isolate 2010EL-1786 (all vs. all comparison) showed congruent gene content and low overall diversity on larger chromosome I (Figure 2). One noteworthy exception was the absence of Vibrio pathogenicity island 1 in the 2005 isolate 3582–05 from Pakistan. This island contains essential cholera virulence factors, including the tcp gene cluster, which encodes toxin-coregulated pilus involved in V. cholerae colonization of the human intestine and necessary for horizontal transfer of the cholera toxin bacteriophage. This finding was the only macroscopic difference observed between isolate 3582–05 and PFGE matches. All Haiti outbreak and PFGE pattern-matched isolates contain an integrated conjugative element belonging to the SXT/R391 family (SXT-ICE) that carries genes conferring antimicrobial drug resistance. No macroscopic differences were observed in SXT-ICE among Haiti outbreak and PFGE pattern-matched isolates (Figure 2; Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 1, panel A).
Smaller chromosome II was more content variable and divergent across study strains. These findings were largely attributable to the hypervariable superintegron region, an ≈120-kb gene capture system predominantly encoding hypothetical proteins (Figure 2; Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 1, panel B) (13). Gene polymorphisms observed in the 9 sequenced isolates from Hispaniola also localized primarily within the superintegron region. Despite these observed differences, no major deletions in the superintegron were observed among PFGE pattern-matched isolates (Figure 2; Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 1). Thus, phylogeny derived from V. cholerae whole-genome sequencing (Figure 1) showed genetic diversity within PFGE pattern-matched isolates. However, binary (present or absent) gene content assessment failed to pinpoint extensive contiguous diversity outside the superintegron region.

Shared Ancestry between Isolates from Haiti, India, and Cameroon

A core genome phylogeny was constructed on the basis of 4,376 hqSNPs found within 632 orthologous core genes (0.81 Mbp) that were universally present in all 27 study and reference genomes (Technical Appendix 1 ; Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 2). Among 9 sequences from Hispaniola isolates, 0–2 SNPs were observed (Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 2). Hispaniola isolates differed from PFGE pattern-matched genomes from other locations by 4–25 SNPs, and genomes with nonmatched PFGE patterns differed from the outbreak isolates by 13–3,361 SNPs. Notably, phylogeny based on hqSNPs showed clustering of the Haiti strain with 3 epidemiologically unrelated clinical isolates, which represented isolates from 2 travelers from the United States to India in 2009 and a patient in Cameroon in 2010. Isolates 2009V-1085 (India, 2009), 2009V-1096 (India, 2009), and 2010EL-1749 (Cameroon, 2010) were most related to the Haiti isolates. These 3 isolates had 4–7 core hqSNPs when compared with the outbreak strain, and the derived sequence for a 2008 clinical isolate from Nepal differed from outbreak isolates by 7–8 core hqSNPs (Technical Appendix 1 ; Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 2).
Conversely, historical isolates (1970–2005) from Pakistan, Bangladesh, the US Gulf Coast, and South America, and recent clinical isolates (2009–2010) from cases linked to Pakistan or South Africa independently clustered away from Haiti outbreak isolates (Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 2). Clade analysis of outbreak isolates and highly related isolates 2009V-1085, 2009V-1096, and 2010EL-1749, identified 25 hqSNPs in 24 conserved loci that distinguish members of this clade (Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 3; Table A1). Resulting distances suggest that the outbreak isolates have a closer genetic relationship with 2009V-1085 and 2009V-1096 from India (7–10 hqSNPs) than with 2010EL-1749 from Cameroon (10–13 hqSNPs).

Comparison of Haiti Outbreak Genomes

Across the 18 described hypervariable V. cholerae mobile genetic elements sequences (representing >300 kb of the total genome), no macroscopic differences were observed among the 9 Hispaniola isolate sequences (Figure 2; Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 1), and as stated, only 2 hqSNPs were identified in the core genome. Pairwise alignment of the complete genome of study reference 2010EL-1786 with available genome data for 2 sequenced Haiti 2010 outbreak isolates, designated H1 and H2 (9), showed only 3 polymorphisms across the entire genome. However, because the available H1 and H2 consensus sequences contain ambiguous basecalls, these nucleotides were excluded from our comparative analyses. Nonetheless, these data confirm the clonal nature of the Haiti outbreak strain.

Structural and Alleleic Profiles of Isolates Carrying a Hybrid Cholera Toxin Prophage

Figure 3

Figure 3. Genetic structure of cholera toxin (CTX) prophage and associated elements in Haiti cholera outbreak Vibrio cholerae isolate 2010EL-1786, fall 2010. The toxin-linked cryptic (TLC) element is not drawn to scale. Black...

Structure and allelic profiles of the CTXϕ prophage have been used for V. cholerae lineage analysis (23). Chromosome I of Haiti isolate 2010EL-1786 harbors 1 hybrid CTXϕ characterized by a 1-nt variant of the classical ctxB allele (ctxB-7) and El Tor rstR flanked by a toxin-linked cryptic element and El Tor–type RS1 element with an intact rstC locus (Figure 3). The SNP at ctxB codon19 results in replacement of the classical cholera toxin B histidine residue with asparagine, and this ctxB-7 allele was observed among all Hispaniola isolates (Table 1). Five of the 12 PFGE pattern-matched isolates from other locations (2008–2010) also shared this variant ctxB allele. The remaining 7 PFGE pattern-matched isolates encoded classical ctxB alleles.

Discussion

Public health investigators use PFGE, the current standard technique for subtyping most bacterial enteric pathogens, to link patients infected with a particular pathogen to a specific infection source(s) by fingerprint matching to pathogens isolated from environmental samples. Whole-genome sequencing has recently emerged as an enhanced laboratory tool for high-resolution analysis of microbial diversity and has been successfully used to investigate bacterial disease outbreaks (24–26). Because whole-genome sequencing can provide pathogen genetic fingerprints at single-nucleotide resolution, it should revolutionize the diagnosis, surveillance, and control of microbial diseases.
For molecular epidemiologic investigations using whole-genome sequencing, an expansive number of isolates from an outbreak would ideally be selected to ensure broad coverage for possible genotype variants within that population that might otherwise be masked with lower-resolution typing methods. In addition, outlier isolates from different locations that are indistinguishable or related by several diverse subtyping methods should also be subjected to whole-genome sequencing to contextualize the diversity seen within the outbreak population and to find other clonal relationships In this study, a temporal and geographic distribution of outbreak isolates was selected to confirm clonality of the outbreak strain and to gain insight into the microevolution of V. cholerae during an outbreak. Additionally, minor PFGE and nonhemolytic variants observed among outbreak isolates were also sequenced to confirm their clonal relationships with isolates exhibiting the main outbreak pattern and phenotype.
The PulseNet USA database substantially contributed to this work by identifying genetically related (using PFGE typing) and epidemiologically relevant isolates for whole-genome sequencing analyses. Notably, one 2008 isolate from a traveler from the United States to Nepal was identified and included in this study, although we acknowledge that the evolutionary relationship of the Haiti strain to strain(s) circulating in Nepal during 2010 may not be ideally represented by this 2008 isolate. Microbial evolution will have occurred during 2008–2010, and global travel may have introduced additional strains into Nepal in the interim, such that the 2008 isolate from Nepal may differ substantially from a strain circulating in Nepal in 2010, the suggested progenitor of the outbreak strain. Unfortunately, 2010 isolates from Nepal were not available for analysis.
Also identified in the PulseNet USA database was 1 PFGE pattern-matched isolate from western Africa. The close genetic relationship of this isolate from Cameroon to the Haiti strain suggests that a potential link between western Africa and the Haiti outbreak cannot be ignored. Further studies on additional isolates from western Africa are required to confirm or refute this possibility. Similarity of whole-genome sequences for Haiti isolates, PFGE pattern-matched isolates, and other seventh pandemic strains confirmed the clonal nature of the 2010–2011 cholera outbreak strain and the close genetic relationships for the studied strains initially suggested by PFGE subtyping (Figure 1). Previous V. cholerae studies have reported that seventh pandemic strains are clonal, sharing near identical gene content on a highly related genome backbone but containing variable mobile genetic elements or gene cassettes (27). Despite dynamic horizontal gene transfer (22), we identified only a few nucleotide differences among mobile sequences of the 9 sequenced 2010–2011 outbreak-related Hispaniola isolates and the 12 recent PFGE pattern-matched clinical isolates (Figure 2).
Extensive recombination in V. cholerae genomes may confound evolutionary relationship analyses as strains and lineages undergo reassortment (1). However, base substitutions acquired horizontally as recombination segments generally occur with localized density (28). Although we cannot guarantee that recombinant segments were absent from the core genome phylogeny (Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 2), the even spatial and genome-wide distribution of core genome hqSNPs suggests that they were vertically inherited. We have derived a useful phylogenetic approximation of isolate relatedness on the basis of hqSNPs, which supports shared ancestry for the Haiti outbreak isolates and 12 recent clinical isolates sharing PFGE patterns (Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 2). Sequenced isolates from India and Cameroon (2009–2010) were shown to be the closest genetic relatives among the non-Hispaniola isolates (isolated in 1991–2010; this study) and 4 other available reference V. cholerae genomes (isolated in 1937–2002). The ctxB allele variant (ctxB-7) of the Haiti strain (and its genetic relatives) was first observed among isolates from a cholera outbreak in Orissa, India, in 2007 (29), but the ctxB-7 allele has since also been observed in isolates from southern Asia and more recently from western Africa (8,30).
The genetic makeup of the Haiti outbreak strain will likely have substantial public health implications for Haiti and other susceptible locations. Our reasoning is that the atypical O1 El Tor V. cholerae strains (CIRS101 and CIRS101-like variants) have already emerged as the predominant clone causing cholera in Asia and Africa and have displaced prototypical O1 El Tor strains (3,4,29). Unfortunately, atypical O1 El Tor V. cholerae strains appear to have retained the relative environmental fitness of their prototypical O1 El Tor ancestors while acquiring enhanced virulence traits, such as classical or hybrid CTX prophage and SXT-ICE (4). Thus, with higher relative fitness and virulent and antimicrobial drug–resistant phenotypes, the Haiti outbreak strain harbors infectivity and ecologic persistence advantages over other seventh pandemic strains. Consequently, the Haiti outbreak strain (or its genetic ancestor) may easily replace current El Tor V. cholerae strains circulating in the Western Hemisphere to become endemic (like other atypical El Tor strains) and will likely cause future outbreaks. Such dire predictions warrant enhanced epidemiologic surveillance and renewed priorities aimed at cholera prevention.
Absence of cholera in Haiti over the past century; the clonal nature of the outbreak strain; and a massive influx of international travelers, aid workers, and supplies after the 2010 earthquake suggest an outside infection source for the 2010–2011 outbreak. Our core genome phylogeny (Technical Appendix 2  [PDF - 1.89 MB - 3 pages] Figure 2) suggests that the Haiti outbreak strain most likely derived from an ancestor related to isolates from within or near the Indian subcontinent. However, concurrent identification of a 2010 isolate from Cameroon as a close genetic relative of the Haiti outbreak strain illustrates that whole-genome sequencing on such a relatively small number (n = 27) of V. cholerae isolates is insufficient to exclude other plausible ancestral geographic locations.
Our study results are consistent with recent findings of Chin et al. (9), who concluded that two 2010 Haiti outbreak isolates shared ancestry with variant O1 El Tor strains isolated in Bangladesh in 2002 and 2008 and a more distant relationship with an isolate from an outbreak in Latin American in 1991. The vertical inheritance pattern of hqSNPs in our study provide unequivocal genetic evidence for introduction of the outbreak strain into Haiti from an external source as opposed to local aquatic emergence. However, the specific geographic source and mode of entry of the outbreak strain into Haiti cannot be proven by microbiological investigations. Only large-scale epidemiologic studies and microbiological data can provide conclusive evidence of how cholera was introduced into Haiti. This whole-genome sequencing study provides expanded evidence that variant O1 El Tor V. cholerae appeared in Haiti by importation and has generated a whole-genome sequencing dataset for future study.

Ms Reimer is a biologist at the Public Health Agency of Canada in Winnipeg, Manitoba, Canada, and project leader for comparative genomics projects of foodborne and waterborne bacterial pathogens. Her research interests are application of bacterial genomics to disease surveillance, outbreak response, and public health priorities.

Acknowledgments

This study was supported in part by Transformational Medical Technologies Program Contract B1042551 from the Department of Defense Chemical and Biological Defense Program through the Defense Threat Reduction Agency.

References

1. Feng  L, Reeves  PR, Lan  R, Ren  Y, Gao  C, Zhou  Z, A recalibrated molecular clock and independent origins for the cholera pandemic clones. PLoS ONE. 2008;3:e4053. DOIPubMed
2. Safa  A, Nair  GB, Kong  RY. Evolution of new variants of Vibrio cholerae O1. Trends Microbiol. 2010;18:46–54. DOIPubMed
3. Nguyen  BM, Lee  JH, Cuong  NT, Choi  SY, Hien  NT, Anh  DD, Cholera outbreaks caused by an altered Vibrio cholerae O1 El Tor biotype strain producing classical cholera toxin B in Vietnam in 2007 to 2008. J Clin Microbiol. 2009;47:1568–71. DOIPubMed
4. Grim  CJ, Hasan  NA, Taviani  E, Haley  B, Chun  J, Brettin  TS, Genome sequence of hybrid Vibrio cholerae O1 MJ-1236, B-33, and CIRS101 and comparative genomics with V. cholerae. J Bacteriol. 2010;192:3524–33. DOIPubMed
5. Pan American Health Organization. Cholera and post-earthquake response in Haiti, 2011 Jul 25 [cited 2011 Aug 6]. http://www.who.int/hac/crises/hti/sitreps/haiti_health_cluster_bulletin_25july2011.pdf .
6. Centers for Disease Control and Prevention. Update on cholera—Haiti, Dominican Republic, and Florida, 2010. MMWR Morb Mortal Wkly Rep. 2010;59:1637–41.PubMed
7. Gilmour  MW, Martel-Laferrière  V, Lévesque  S, Gaudreau  C, Bekal  S, Nadon  C, Vibrio cholerae in traveler from Haiti to Canada. Emerg Infect Dis. 2011;17:1124–5.PubMed
8. Talkington  D, Bopp  C, Tarr  C, Parsons  MB, Dahourou  G, Freeman  M, Characterization of toxigenic Vibrio cholerae from Haiti, 2010–2011. Emerg Infect Dis. 2011;17:2122–9.
9. Chin  CS, Sorenson  J, Harris  JB, Robins  WP, Charles  RC, Jean-Charles  RR, The origin of the Haitian cholera outbreak strain. N Engl J Med. 2011;364:33–42. DOIPubMed
10. Archer  BN, Cengimbo  A, De Jong  GM, Keddy  KH, Smith  AM, Sooka  A, Cholera outbreak in South Africa: preliminary descriptive epidemiology on laboratory-confirmed cases, 15 November 2008 to 30 April 2009. Communicable Diseases Surveillance Bulletin. 2009;7:3–8.
11. Sakazaki  R, Shimada  T. Serovars of Vibrio cholerae identified during 1970–1975. Jpn J Med Sci Biol. 1977;30:279–82.PubMed
12. Sakazaki  R, Tamura  K, Gomez  CZ, Sen  R. Serological studies on the cholera group of vibrios. Jpn J Med Sci Biol. 1970;23:13–20.PubMed
13. Heidelberg  JF, Eisen  JA, Nelson  WC, Clayton  RA, Gwinn  ML, Dodson  RJ, DNA sequence of both chromosomes of the cholera pathogen Vibrio cholerae. Nature. 2000;406:477–83. DOIPubMed
14. Delcher  AL, Bratke  KA, Powers  EC, Salzberg  SL. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics. 2007;23:673–9. DOIPubMed
15. Meyer  F, Goesmann  A, McHardy  AC, Bartels  D, Bekel  T, Clausen  J, GenDB–an open source genome annotation system for prokaryote genomes. Nucleic Acids Res. 2003;31:2187–95. DOIPubMed
16. Darling  AE, Mau  B, Perna  NT. Progressive Mauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE. 2010;5:e11147. DOIPubMed
17. Guindon  S, Delsuc  F, Dufayard  JF, Gascuel  O. Estimating maximum likelihood phylogenies with PhyML. Methods Mol Biol. 2009;537:113–37. DOIPubMed
18. Thompson  JD, Higgins  DG, Gibson  TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80. DOIPubMed
19. Anisimova  M, Gascuel  O. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol. 2006;55:539–52. DOIPubMed
20. Petkau  A, Stuart-Edwards  M, Stothard  P, Van Domselaar  G. Interactive microbial genome visualization with GView. Bioinformatics. 2010;26:3125–6. DOIPubMed
21. Mason  PR. Zimbabwe experiences the worst epidemic of cholera in Africa. J Infect Dev Ctries. 2009;3:148–51. DOIPubMed
22. Chun  J, Grim  CJ, Hasan  NA, Je  HL, Seon  YC, Haley  BJ, Comparative genomics reveals mechanism for short-term and long-term clonal transitions in pandemic Vibrio cholerae. Proc Natl Acad Sci U S A. 2009;106:15442–7. DOIPubMed
23. Olsvik  O, Wahlberg  J, Petterson  B, Uhlen  M, Popovic  T, Wachsmuth  IK, 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. 1993;31:22–5.PubMed
24. Gilmour  MW, Graham  M, Van Domselaar  G, Tyler  S, Kent  H, Trout-Yakel  KM, High-throughput genome sequencing of two Listeria monocytogenes clinical isolates during a large foodborne outbreak. BMC Genomics. 2010;11:120.PubMed
25. Gardy  JL, Johnston  JC, Ho Sui  SJ, Cook  VJ, Shah  L, Brodkin  E, Whole genome sequencing and social network analysis of a tuberculosis outbreak. N Engl J Med. 2011;364:730–9. DOIPubMed
26. Taboada  EN, van Belkum  A, Yuki  N, Acedillo  RR, Godschalk  PC, Koga  M, Comparative genomic analysis of Campylobacter jejuni associated with Guillain-Barré and Miller Fisher syndromes: neuropathogenic and enteritis-associated isolates can share high levels of genomic similarity. BMC Genomics. 2007;8:359. DOIPubMed
27. Cho  YJ, Yi  H, Lee  JH, Kim  DW, Chun  J. Genomic evolution of Vibrio cholerae. Curr Opin Microbiol. 2010;13:646–51. DOIPubMed
28. Croucher  NJ, Harris  SR, Fraser  C, Quail  MA, Burton  J, van der Linden  M, Rapid pneumococcal evolution in response to clinical interventions. Science. 2011;331:430–4. DOIPubMed
29. Goel  AK, Jain  M, Kumar  P, Bhadauria  S, Kmboj  DV, Singh  L. A new variant of Vibrio cholerae O1 El Tor causing cholera in India. J Infect. 2008;57:280–1. DOIPubMed
30. Quilici  ML, Massenet  D, Gake  B, Bwalki  B, Olson  DM. Vibrio cholerae O1 variant with reduced susceptibility to ciprofloxacin, western Africa. Emerg Infect Dis. 2010;16:1804–5.PubMed

Figures

• Figure 1. Whole-genome neighbor-joining tree of Vibrio cholerae isolate from cholera outbreak in Haiti, fall 2010; concurrent clinical isolates with pulsed-field gel electrophoresis pattern-matched combinations; reference isolates sequenced in this study;...
• Figure 2. BLAST (http://blast.ncbi.nlm.gov/Blast.cgi) atlas of predicted protein homologies mapped against the closed genome of Haiti Vibrio cholerae outbreak type isolate 2010EL-1786, fall 2010. Full color saturation represents 100% sequence homology,...
• Figure 3. Genetic structure of cholera toxin (CTX) prophage and associated elements in Haiti cholera outbreak Vibrio cholerae isolate 2010EL-1786, fall 2010. The toxin-linked cryptic (TLC) element is not drawn to...

Tables

• Table 1. Characteristics of Vibrio cholerae isolates from Haiti, Asia, Africa, and the United States
• Table 2. Next-generation sequence average coverage and number of mapped reads for Vibrio cholerae isolates from Haiti, Asia, and Africa
• Table A1. Conserved open reading frames and hqSNPs used to refine evolutionary relationship between highly related outbreak Vibrio cholerae isolates and isolates 2009V-1085 (India), 2009V-1096 (India), and 2010EL-1749 (Cameroon)

Technical Appendices

• Technical Appendix 1.    [440 KB - 3 pages]
• Technical Appendix 2.    [1.89 MB - 3 pages]

Suggested citation for this article: Reimer AR, Van Domselaar G, Stroika S, Walker M, Kent H, Tarr C, et al. Comparative genomics of Vibrio cholerae from Haiti, Asia, and Africa. Emerg Infect Dis [serial on the Internet]. 2011 Nov [date cited]. http://dx.doi.org/10.3201/eid1711.110794

DOI: 10.3201/eid1711.110794

¹Members of the V. cholerae Outbreak Genomics Task Force are Arunmozhi Balajee, Shanna Bolcen, Cheryl A. Bopp, John Besser, Ifeoma Ezeoke, Patricia Fields, Molly Freeman, Lori Gladney, Dhwani Govil, Michael S. Humphrys, Maria Sjölund-Karlsson, Karen H. Keddy, Elizabeth Neuhaus, Michele M. Parsons, Efrain Ribot, Maryann Turnsek, Shaun Tyler, Jean M. Whichard, Anne Whitney, and the authors.

Comparative Genomics of Vibrio cholerae from Haiti, Asia, and Africa - Vol. 17 No. 11 - November 2011 - Emerging Infectious Disease journal - CDC