Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 28, Number 11—November 2022

Genomic Epidemiology of Vibrio cholerae O139, Zhejiang Province, China, 1994–2018

Yun Luo, Julian Ye, Michael Payne, Dalong Hu, Jianmin Jiang1, and Ruiting Lan1Comments to Author 
Author affiliations: Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, China (Y. Luo, J. Ye, J. Jiang); University of New South Wales, Sydney, New South Wales, Australia (Y. Luo, M. Payne, D. Hu, R. Lan); Key Laboratory of Vaccine, Hangzhou (J. Jiang)

Cite This Article


Cholera caused by Vibrio cholerae O139 was first reported in Bangladesh and India in 1992. To determine the genomic epidemiology and origins of O139 in China, we sequenced 104 O139 isolates collected from Zhejiang Province, China, during 1994–2018 and compared them with 57 O139 genomes from other countries in Asia. Most Zhejiang isolates fell into 3 clusters (C1–C3), which probably originated in India (C1) and Thailand (C2 and C3) during the early 1990s. Different clusters harbored different antimicrobial resistance genes and IncA/C plasmids. The integrative and conjugative elements carried by Zhejiang isolates were of a new type, differing from ICEVchInd4 and SXTMO10 by single-nucleotide polymorphisms and presence of genes. Quinolone resistance–conferring mutations S85L in parC and S83I in gyrA occurred in 71.2% of the Zhejiang isolates. The ctxB copy number differed among the 3 clusters. Our findings provided new insights for prevention and control of O139 cholera .

Cholera is an acute watery diarrheal disease that has caused 7 global pandemics since 1817. The current, ongoing seventh pandemic started in 1961 and continues today (1). The causative agent of cholera is Vibrio cholerae; serogroups O1 and O139 cause epidemic- and pandemic-level disease. Serogroup O139 first caused an outbreak in Bangladesh and India in 1992 (2,3). However, the epidemic O139 clone was later found to be a derivative of a seventh pandemic O1 strain, having had its O1 gene cluster replaced with an O139 O antigen gene cluster (4) and therefore genetically belonging to the seventh pandemic clone and sharing the same sequence type (5).

V. cholerae O139 spread to China and was reported in Xinjiang Uygur Autonomous Region (6) and Guangdong Province (7) in 1993 and in Jiangxi Province and the cities of Beijing and Shanghai in 1994 (8). Most studies on O139 in China have focused on virulence and resistance-gene profiles, cholera toxin (CTX) types (7,8), and plasmid carriage (9). The genomic epidemiology of O139 in China and the phylogenetic relationship of isolates from China to isolates from other countries in Asia are still unknown. A study of 9 O139 isolates suggested that O139 reached China soon after outbreaks in India in early 1990s and became dominant a few years later (10).

Antimicrobial therapy (in addition to rehydration therapy) plays a vital role in the management of cholera patients (11). In a previous study of 340 O139 isolates collect in China during 1993–2009, resistances to streptomycin, trimethoprim/sulfamethoxazole, and polymyxin B were found in isolates from early years (12). IncA/C conjugative plasmids can effectively mobilize genes associated with resistance to different classes of antibiotics, including β-lactams, aminoglycosides, chloramphenicol, folate-pathway inhibitors, quinolones, and tetracycline (13). IncA/C plasmids are widely present in Enterobacterales but not common in V. cholerae populations (14), although they have been found in the seventh pandemic V. cholerae lineage (11).

In this study, we sequenced the genomes of 104 V. cholerae O139 isolates collected from Zhejiang Province, China, during 1994–2018. Comparative genomic and phylogenetic analyses revealed the genetic characteristics of V. cholerae O139 isolates in Zhejiang and their evolutionary relationships to isolates from countries in Asia. We also analyzed the virulence and antimicrobial resistance (AMR) gene profiles and the distribution of IncA/C plasmids to elucidate the evolution of virulence and AMR.



We recovered 104 V. cholerae O139 isolates collected during 1994–2018 from the Zhejiang Provincial Center for Disease Control and Prevention (Zhejiang CDC). We downloaded 57 public V. cholerae genomes from countries in Asia (Appendix 1 Table 1) and 133 publicly available V. cholerae genomes from China from the European Nucleotide Archive database ( and identified them by searching for the O139 O-antigen–specific wbf gene using BLASTN version 2.9.0 (15).

Genome Sequencing

We performed whole-genome sequencing by using the Illumina Hiseq X-ten sequencing platform with TruePrepTM DNA Library Prep Kit version 2 and 150-bp paired-end sequencing (Illumina, We checked all input read sets for contamination by using kraken2 with a threshold of 10% for non–V. cholerae reads (16). We submitted genome sequences obtained in this study as raw reads under the National Center for Biotechnology Information’s Sequence Read Archive database (Bioproject no. PRJNA643344).

Single-Nucleotide polymorphism Calling and Phylogenetic Analyses

We identified single-nucleotide polymorphisms (SNPs) by using a section of the SaRTree (17) pipeline. We removed Superintegron sequences on the small chromosome and all recombinant SNPs. The reference genome sequence (GenBank accession no. GCF_900324445.1) was from Bangladesh strain 4295STDY6534216, isolated in 2014 (18). We allocated SNPs to each branch of the tree by using the SaRTree pipeline (17). We performed phylogenetic analysis by constructing a maximum-likelihood tree using IQ-Tree version 2.0.4 (19) under default parameters (transversion model with AG = CT and empirical base frequencies) with 1,000 bootstrap replicates.

Antimicrobial-Resistance and Virulence Genes

For all genomes, we predicted AMR genes by using ABRicate ( with the AMRFinderPlus gene database (20), plasmids by using PlasmidFinder (21), and virulence genes by using a customized database of 67 virulence genes (Appendix 1 Table 2). We applied a cutoff of percentage nucleotide identity at 80% for virulence genes and plasmids and at 60% for resistance genes. We used k-mer alignment (22) to map raw reads against all these genes. As criteria for gene presence, we used a combination of minimum identity and coverage thresholds from ABRicate or the ratio of the gene depth to the average depth of housekeeping genes >20% from KMA. We used CNVnator version 0.4.1 (23) with default settings and a bin size of 100 bp to calculate copy numbers of ctxB genes.


Whole-Genome Sequencing of V. cholerae O139 Isolates from Zhejiang Province

Figure 1

Distribution of Vibrio cholerae O139 isolates, by clusters and year of isolation, Zhejiang Province, China, 1994–2018. Bar sections represent isolate numbers in different clusters in each year (Appendix 1 Table 4; Appendix 2 Figure 2).

Figure 1. Distribution of Vibrio cholerae O139 isolates, by clusters and year of isolation, Zhejiang Province, China, 1994–2018. Bar sections represent isolate numbers in different clusters in each year (...

We recovered and sequenced 104 V. cholerae O139 isolates collected during 1994 to 2018 by Zhejiang CDC (Figure 1; Appendix 1 Table 1). For comparison with other isolates from China, we included 9 published isolates from Shanghai (10). We used another 133 publicly available O139 genomes from China without metadata only to infer phylogenetic relationships with Zhejiang isolates. For international comparison, we included 48 publicly available O139 genomes from India (19), Bangladesh (19), and Thailand (10). The earliest isolates were collected in 1983 in Bangladesh; other isolates were collected during 1992–2014 (Appendix 2 Figure 1).

Phylogenetic Analysis of O139 Isolates from China and Worldwide

Figure 2

Maximum-likelihood phylogenetic tree of 161 Vibrio cholerae O139 (sequence type 69) isolates from Zhejiang Province, China, 1994–2018, and isolates from outside of China. The tree was rooted using the seventh pandemic O1 strain N16961 as an outgroup. Lineage 1 (L1) and lineage 2 (L2) are demarcated with red dashed lines and pink dashed-line boxes. The 3 clusters (C1, C2, and C3) are collapsed to reduce figure size (Appendix 2 Figure 2). Key branches are marked with a branch number followed in brackets by the number of single nucleotide polymorphisms that supported the branch. The colored solid squares at the end of isolate names indicate the location of isolation of the isolates. GenBank accession numbers were used as isolate names for O139 isolates not from Zhejiang Province. SH, Shanghai; ZJ, Zhejiang.

Figure 2. Maximum-likelihood phylogenetic tree of 161 Vibrio choleraeO139 (sequence type 69) isolates from Zhejiang Province, China, 1994–2018, and isolates from outside of China. The tree was rooted using the...

We identified 629 SNPs from our 104 Zhejiang O139 isolates, 9 Shanghai isolates, and 48 isolates from other countries; 501 SNPs were on chromosome I and 128 SNPs were on chromosome II. We constructed a phylogenetic tree of the 161 isolates, using N16961 as the outgroup (Figure 2; Appendix 2 Figure 2). We further identified branch-supporting SNPs (Appendix 1 Table 3; Appendix 2 Figure 3). The tree can be divided into 2 distinctive linages, defined as lineages 1 (L1) and 2 (L2). L1 contained 11 isolates from this study and 15 isolates from Bangladesh and India. L2 contained 92 isolates from this study and 8 Shanghai isolates.

The isolates from China grouped together as 3 clusters, which were each supported by unique SNPs (Appendix 2 Figures 2, 3). Cluster 1 (C1) consisted of 11 Zhejiang isolates collected during 1996–2000 supported by 3 SNPs on branch 233. Cluster 2 (C2) consisted of 59 isolates from this study and 4 Shanghai isolates collected during 1998–2008. C2 was supported by 1 SNP on branch 49. Cluster 3 (C3) consisted of 28 isolates from this study and 4 isolates from Shanghai collected during 1998–2018. C3 was supported by 1 SNP (branch 45). C1 was located within L1, whereas C2 and C3 were located within L2. The closest ancestral isolate of C1 was an isolate from India. C2 and C3 grouped together as L2, and all isolates in L2 originated in China. The L2 node was supported by 3 SNPs (branch 177), and the closest ancestral isolate of L2 was from Thailand.

Six isolates from Zhejiang fell outside of C2 and C3 and were referred to as outliers. These 6 isolates were obtained from patients in the 1990s. One isolate (V01) was isolated in 1994 from the first clinical case-patient in Zhejiang but shared a common ancestor with 1 isolate from Shanghai (isolated in 1994) and was sibling to L2.

We also constructed a phylogenetic tree of 349 genomes that included the 133 isolates from China without metadata (Appendix 2 Figure 4). The 2 lineages (L1 and L2) and 3 clusters were preserved on this tree. Most isolates from China fell into C2 (92/133 [69.2%]), 9 isolates fell into C3, and only 2 isolates fell into C1.

Genetic Elements and Virulence Genes

We searched the 104 genomes from this study by using ABRicate with our custom database of 67 known V. cholerae virulence genes. We further confirmed by reads mapping any genomes with a negative result for any of these virulence genes. The presence of a gene was a combined result of ABRicate searches of assembled genomes and reads mapping. Four core CTX phage genes (ace, zot, ctxA, and ctxB) and the repeat sequences were all present in 94 genomes. Three genomes (V31, V32, and V33) were all negative for these genes. The repeat sequences were not well assembled and on different contigs, whereas ctxAB, zot, and ace of 89/104 isolates were on the same contig.

All 104 genomes contained 19 Vibrio pathogenicity island genes, 18 genes on 2 Vibrio seventh pandemic islands, and 19 type VI secretion system–related genes. All but 1 genome contained the intact repeats-in-toxin gene cluster. (Appendix 1 Table 2). We found Vibrio pathogenicity island genes on the same contigs in 99/104 genomes, Vibrio seventh pandemic island I genes on the same contigs in 103/104 genomes, and Vibrio seventh pandemic island II genes on the same contigs in 104/104 genomes.

Because a strain may contain multiple copies of the CTX phage (24), we used CNVnator to estimate the number of copies of the ctxB gene in the 104 isolates by mapping reads to V. cholerae seventh pandemic reference genome N16961. The ctxB gene copies differed in the 3 clusters; on average, C2 had 4.4 copies/isolate, C3 had 1.2 copies/isolate, and C1 had 1.3 copies/isolate. A total of 68 isolates carried multiple copies of ctxB (range 4–22 copies) (Appendix 1 Table 2), whereas 38 isolates carried only 1 copy of ctxB. Most C2 isolates (81.4% [48/59]) carried >2 copies (average 5 copies) of ctxB. In contrast, only 3 isolates (10.7% [3/28]) in C3 carried >1 copy (average 2.6 copies), and 6 isolates (54.5% [6/11]) in C1 carried 2 copies. Six outlier isolates carried multiple copies of ctxB (range 2–22 copies). All 101 ctxB–positive isolates contained ctxB genotype 3.

Antimicrobial-Resistance Genes and Resistance Mutations of the O139 Isolates

We found the chromosomally encoded resistance genes varG and catB9 in all isolates. floR, dfr18, sul2, aph(3′′)-lb, and aph(6)-ld were in most of the isolates, including international isolates, and 8 AMR genes detected only in cluster 3 were present at low frequencies (3.13%–15.63%) (Table). We found blaTEM-1, catA2, aac(3)-lld, aadA2, aph(3′)-la, mph(E), msr(E), sul, dfrA12, tet(M), and tet(Y) only in isolates from China, including Shanghai isolates, whereas blaCMY-2, blaOXA-1, catB3, aac(6′)-lb-cr5, aadA3, ere(A), mph(A), mph(F), qnrA1, qnrA7, aar-3, dfrA27, dfrA32, tet(A), and tet(D) were in Zhejiang isolates only (Appendix 1 Table 4). C1 carried only AMR genes common to all isolates, aph(3′)-Ia and sul1 were common to C2 and C3, aac(3)II, aadA2, tet(D), mph(E), msr(E), blaTEM-1, and catA2 were more common in C2, and tet(M), mph(A), and dfrA12 were more common in C3.

We also searched these isolates for quinolone-resistance mutations. Seventy-four isolates (74/104 [71.2%]) harbored mutations Ser85Leu in parC and Ser83Ile in gyrA. Ten isolates had a mutation in Asp87 of gyrA, of which 6 had Asp87Tyr, 3 had Asp87Gly, and 1 had Asp87Asn.

Association of Plasmids and Integrative and Conjugative Elements with AMR Genes

We analyzed the integrative and conjugative elements (ICEs) carried by our isolates and compared them with the 2 known ICE variants in O139 (SXTMO10 and ICEVchInd4) (25). SNP and phylogenetic analyses found that all O139 ICEs were closely related (difference of 0–13 SNPs) (Appendix 2 Figure 5). Our ICEs differed from ICEVchInd4 by 0–7 SNPs and from SXTMO10 by 10–13 SNPs. However, most of our isolates carried the 4 SXTMO10 genes, including dfrA18 that are absent in ICEVchInd4. The AMR genes present in our isolates, dfrA18, floR, aph(3′′)-Ib, aph(6)-Id, and sul2, were probably carried by the ICE.

Eighty-three isolates carried an IncA/C plasmid. We found 2 IncA/C subtypes (IncA/C2_1_JN157804 type, belonging to plasmid pNDM-KN-lineage, and IncA/C_1_FJ705807 type, belonging to pRA1-lineage) (Appendix 1 Table 4). None of the C1 isolates contained an IncA/C plasmid. All except 2 C2 isolates contained an IncA/C2 plasmid (Appendix 1 Table 4). k-mer alignment analysis indicated that these C2 isolates carried a plasmid nearly identical to the known V. cholerae plasmid pVC1447 of 160 kb (9). pVC1447 is known to carry aadA, sul1, tetD, blaTEM, catA2, mph(E), tet(R), mel, qacEdelta1, and folP genes (9). The last 4 AMR genes were not found in any of our C2 isolates. Most C2 isolates carried aadA2, sul1, tetD, catA2, mph(E), msr(E), and blaTEM. Seven C2 isolates lost >1 of the AMR genes. Isolate V29 and V30 lost sul1, tetD, mph(E), and msr(E) genes. The 2 C2 isolates without the pVC1447-like plasmid did not contain any of the pVC1447 AMR genes. All except 4 C3 isolates carried an IncA/C plasmid; 17 had the IncA/C_1_FJ705807 replicon type, and 5 had the IncA/C2_1_JN157804 replicon type. We further determined that the IncA/C_1_FJ705807 type plasmid is a novel plasmid that was most closely related to Aeromonas veronii plasmid p158496 (26), whereas the IncA/C2_1_JN157804 type plasmid was most similar to V. cholerae O139 pVC211 (GenBank accession no. KY399978.1). The p158496-like plasmid in the 17 C3 isolates shared an average nucleotide identity of 97.7% and length coverage of 82.89% with the 158 kb p158496 and probably carried aadA2, tet(D), tet(M), mph(A), dfrA12, and sul1 genes. However, more than half of these C3 isolates lost the tet(D) gene. Two outlier isolates (V17 and V18) also carried the p158496-like plasmid. The pVC211-like plasmid in the 5 C3 isolates shared an average nucleotide identity of 99.09% and length coverage of 92.46% with the 148 kb pVC211 plasmid and probably carried aadA16, tet(A), mph(A), dfrA27, qnrA1, and arr-3 resistance genes. Some isolates had further loss and gain of AMR genes.


The first V. cholerae O139 isolate in Zhejiang Province was reported in September 1994, which was 16 months after the first O139 case reported in China (6). Phylogenetic analysis grouped Zhejiang isolates into 2 independent lineages (L1 and L2) and 3 clusters (C1, C2, and C3). The origin of C1 was probably India and the origin of L2 (C2 and C3) was probably Thailand. However, considerable uncertainty exists, as it does with L1, the sister clade of C1, which contained both India and Bangladesh isolates. Similarly, L2, which contained C2 and C3 of isolates from China only, shared a most recent common ancestor with isolates from India, Bangladesh, and Thailand. More isolates from the other countries in Asia would be required to resolve the origins of the clusters in China.

Other isolates in China also fell into the 3 Zhejiang clusters, suggesting that these clusters were circulating across China. However, because the publicly available O139 genomes from other parts of China contained no location metadata, we cannot infer whether O139 reached Zhejiang first and then spread to other parts of China or vice versa.

All isolates in this study were ctxB genotype 3. However, a study of isolates from south China found that a small proportion of ctxB genotypes 1 and 5 in isolates from the 1990s, although >90% of the isolates were ctxB genotype 3 (7). Most C2 isolates carried multiple copies of ctxB, suggesting that the cluster carried multiple copies of the CTX phage. The number of CTX carried by O139 may vary (24). In our study, we observed that the variation in the number of CTX carried occurred along lineages. The higher number of ctx copies might lead to greater toxin production, potentially affecting disease outcomes.

On the basis of the presence of AMR genes and resistant mutations, we determined that the evolution of resistance to antimicrobials changed substantially over time. Tetracycline resistance genes tet(M) and tet(Y) were present only in isolates in China, and tet(A) and tet(D) were only present in Zhejiang isolates. tet(M) was found in 65.6% of C3 isolates; some C3 isolates carried both tet(D) and tet(M). Previous studies found that O139 isolates from 1991–2013 in Thailand and from 1997 in India were susceptible to tetracycline (27,28), suggesting that the earlier O139 isolates in Asia did not carry the tet genes. Because tetracyclines were overused in China (29), it is not surprising that C2 and C3 isolates acquired tetracycline-resistance genes, and these events probably occurred in China.

mph(A) was present in 68.8% of C3 isolates only. mph(A) conferring azithromycin resistance is plasmid-borne and rarely found in V. cholerae (30). Azithromycin was first used in clinical treatment in 1988 (31). In our study, mph(A) was first identified in an isolate in 1998, only 10 years after azithromycin was first used for treatment. The high percentage of mph(A) in O139 C3 isolates in this study is concerning. However, because cholera was relatively infrequent in China, the acquisition of such resistance may not be attributable to selection pressure from clinical antimicrobial treatment.

C2 and C3 shared 2 nonsynonymous mutations, 1 each in the genes encoding for penicillin-binding protein 2 and a lytic murein transglycosylase with affinity to β-lactam antibiotic resistance (32). These genomic changes were previously reported in Shanghai O139 isolates and were attributed to the increasing usage of β-lactam antibiotics (10). These mutations were present only in isolates originating in China in L2 and may have evolved in China.

Plasmid analysis found that C2 and C3 isolates acquired different IncA/C plasmids. C2 carried a known plasmid (pVC1447), whereas C3 isolates acquired 2 different IncA/C plasmids, A. veronii plasmid p158496-like and V. cholerae pVC211-like. These plasmids were probably the carriers of the new AMR genes in different clusters and contributed to the differences of AMR gene profiles between clusters.

Although O139 spread to China in 1993 (6), our earliest isolates in Zhejiang were from 1994 and did not belong to any of the 3 clusters. Five more unclustered isolates were from 1996–1998, all of which belong to L2. Therefore, in early years of the O139 epidemic, multiple independent introductions of O139 cholera to Zhejiang directly from other countries or indirectly from other parts of China had probably occurred. However, the 3 clusters flourished at different times were successively replaced during 1994–2018. C1 was found in 1996 and persisted until 2000, C2 during 1998–2008, and C3 during 1998–2018. The earliest isolate from both C2 and C3 were 1998, suggesting that C2 and C3 were imported to China at similar times or a single importation of the most recent common ancestor of C2 and C3 had occurred from which the 2 clusters diverged in China. C2 became a dominant population in Zhejiang during 2000–2008 and then C3 took over from 2009, replacing the other clusters. Therefore, Zhejiang experienced O139 cholera in 3 waves caused by 3 clusters, each lasting up to a decade.

The epidemiologic pattern uncovered raises many interesting questions, most notably regarding what advantage did subsequent clusters have over their predecessor. C2 and C3 carried more resistance genes than C1. Although we have no AMR phenotypic data, difference in their AMR gene profiles suggests that AMR may have been the driver that caused C1 to be replaced by C2 and C3. C2 overall carried more copies of ctxB, suggesting that it may produce more CTX toxin than C3. Nearly 50% of the isolates from other parts of China from the unpublished genomes belonged to C2, suggesting that C2 was quite prevalent and more successful than C3. The increased number of copies of ctxB probably contributed to its success in replacing C1. However, this explanation does not account for why C2 was subsequently replaced by C3 in Zhejiang. Again, C3 acquired resistance to several additional AMR genes that may explain its fitness advantage over C2, given that tet(M) and mph(A) were only detected in C3 and dfrA12 was mainly present in C3. In addition, the AMR genes present at low frequency in different C3 isolates may have also collectively contributed to C3’s fitness. However, 63% of C2 isolates simultaneously carried blaTEM1, catA2, aac(3)-lld, aadA2, aph(3′)-la, sul1, mph(E), and msr(E) genes, a pattern not present in other clusters.

The persistence of each cluster for many years in Zhejiang is also intriguing. The clusters possibly were circulating in other parts of China and spread to Zhejiang. Most of the other isolates in China fell into C2, and Shanghai isolates were shown to be ancestral to some Zhejiang isolates within C2 and C3 (Appendix 2 Figure 2), supporting this hypothesis. Isolates may have also been continuously imported from other countries. However, we have no isolates from other countries of corresponding years to examine this hypothesis. Another possibility is that O139 has spread to the environment in Zhejiang, where it has established itself as a local reservoir. However, our extensive sampling of river waters over 2 years in 2 cities in Zhejiang only found non-O1/non-O139 isolates (33), although the sampling done in that study had only 2 years overlap with the isolation years of O139 isolates from those cities. Thus, it is less likely that these O139 cases were from local environmental reservoirs. A recent study of cholera in Africa also found repeated importation rather than local environmental reservoirs as the source of the seventh pandemic cholera during cholera resurgence over a 40-year period (11).

This study describes the possible origin, evolution, and spread of O139 cholera in a single province, Zhejiang. Further studies are required to expand this analysis to the national level. Most V. cholerae O139 isolates in Zhejiang grouped into 3 major clusters, which were probably derived from multiple independent importation events directly or indirectly from other countries in Asia and prevailed over the period 1994–2018, with one cluster replacing another sequentially. Variations in AMR gene content or resistance mutations suggest that acquisition of AMR probably has played a role in the succession of the V. cholerae O139 clusters in Zhejiang.

Ms. Luo is a PhD candidate in the School of Biotechnology and Biomolecular Sciences, University of New South Wales. Her primary research interests are the genomic epidemiology of enteric pathogens and control of infectious diseases.



We thank the staff of Zhejiang CDC for their help with the isolation and conservation of the V. cholerae isolates used in this study. We also thank Shichang Xia for his conceptualization and support for this project. Y.L. is a PhD student supported by Australian Government Research Training Program Scholarship.

Author contributions: Y.L., J.J., and R.L. conceived of the study. J.Y. curated the data and resources. Y.L. did the data analysis and wrote the draft of the manuscript. M.P. and D.H. developed the methodology. M.P., D.L., and R.L. verified the underlying data. M.P. and R.L. reviewed and edited the manuscript. All authors reviewed and approved the final version of the manuscript. All authors had full access to all the data in this study. R.L. and J.J. had final responsibility for the decision to submit for publication.



  1. Clemens  JD, Nair  GB, Ahmed  T, Qadri  F, Holmgren  J. Cholera. Lancet. 2017;390:153949. DOIPubMedGoogle Scholar
  2. Ramamurthy  T, Garg  S, Sharma  R, Bhattacharya  SK, Nair  GB, Shimada  T, et al. Emergence of novel strain of Vibrio cholerae with epidemic potential in southern and eastern India. Lancet. 1993;341:7034. DOIPubMedGoogle Scholar
  3. Choleraworkinggroupinternation ; Cholera Working Group, International Centre for Diarrhoeal Diseases Research, Bangladesh. Large epidemic of cholera-like disease in Bangladesh caused by Vibrio cholerae O139 synonym Bengal. Lancet. 1993;342:38790. DOIPubMedGoogle Scholar
  4. Stroeher  UH, Jedani  KE, Dredge  BK, Morona  R, Brown  MH, Karageorgos  LE, et al. Genetic rearrangements in the rfb regions of Vibrio cholerae O1 and O139. Proc Natl Acad Sci U S A. 1995;92:103748. DOIPubMedGoogle Scholar
  5. Karaolis  DK, Lan  R, Reeves  PR. The sixth and seventh cholera pandemics are due to independent clones separately derived from environmental, nontoxigenic, non-O1 Vibrio cholerae. J Bacteriol. 1995;177:31918. DOIPubMedGoogle Scholar
  6. Liu  J, Gao  S, Gao  T, Qi  G, Cao  X, Duan  G, et al. Diarrhea outbreak by Vibrio cholerae O139 in Keping Xinjiang [in Chinese]. Dis Surveill. 1993;8:2389.
  7. Li  BS, Xiao  Y, Wang  DC, Tan  HL, Ke  BX, He  DM, et al. Genetic relatedness of selected clinical Vibrio cholerae O139 isolates from the southern coastal area of China over a 20-year period. Epidemiol Infect. 2016;144:267987. DOIPubMedGoogle Scholar
  8. Zhang  P, Zhou  H, Diao  B, Li  F, Du  P, Li  J, et al. A molecular surveillance reveals the prevalence of Vibrio cholerae O139 isolates in China from 1993 to 2012. J Clin Microbiol. 2014;52:114652. DOIPubMedGoogle Scholar
  9. Wang  R, Liu  H, Zhao  X, Li  J, Wan  K. IncA/C plasmids conferring high azithromycin resistance in vibrio cholerae. Int J Antimicrob Agents. 2018;51:1404. DOIPubMedGoogle Scholar
  10. Hu  D, Yin  Z, Yuan  C, Yang  P, Qian  C, Wei  Y, et al. Changing molecular epidemiology of Vibrio cholerae outbreaks in Shanghai, China. mSystems. 2019;4:e0056119. DOIPubMedGoogle Scholar
  11. Weill  F-X, Domman  D, Njamkepo  E, Tarr  C, Rauzier  J, Fawal  N, et al. Genomic history of the seventh pandemic of cholera in Africa. Science. 2017;358:7859. DOIPubMedGoogle Scholar
  12. Yu  L, Zhou  Y, Wang  R, Lou  J, Zhang  L, Li  J, et al. Multiple antibiotic resistance of Vibrio cholerae serogroup O139 in China from 1993 to 2009. PLoS One. 2012;7:e38633. DOIPubMedGoogle Scholar
  13. Carraro  N, Matteau  D, Burrus  V, Rodrigue  S. Unraveling the regulatory network of IncA/C plasmid mobilization: When genomic islands hijack conjugative elements. Mob Genet Elements. 2015;5:15. DOIPubMedGoogle Scholar
  14. Carraro  N, Rivard  N, Ceccarelli  D, Colwell  RR, Burrus  V. IncA/C conjugative plasmids mobilize a new family of multidrug resistance islands in clinical Vibrio cholerae non-O1/non-O139 isolates from Haiti. MBio. 2016;7:e0050916. DOIPubMedGoogle Scholar
  15. Bhumiratana  A, Siriphap  A, Khamsuwan  N, Borthong  J, Chonsin  K, Sutheinkul  O. O serogroup-specific touchdown-multiplex polymerase chain reaction for detection and identification of Vibrio cholerae O1, O139, and non-O1/non-O139. Biochem Res Int. 2014;2014:295421. DOIPubMedGoogle Scholar
  16. Wood  DE, Salzberg  SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15:R46. DOIPubMedGoogle Scholar
  17. Hu  D, Liu  B, Wang  L, Reeves  PR. Living trees: high-quality reproducible and reusable construction of bacterial phylogenetic trees. Mol Biol Evol. 2020;37:56375. DOIPubMedGoogle Scholar
  18. Dorman  MJ, Domman  D, Uddin  MI, Sharmin  S, Afrad  MH, Begum  YA, et al. High quality reference genomes for toxigenic and non-toxigenic Vibrio cholerae serogroup O139. Sci Rep. 2019;9:5865. DOIPubMedGoogle Scholar
  19. Minh  BQ, Schmidt  HA, Chernomor  O, Schrempf  D, Woodhams  MD, von Haeseler  A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37:15304. DOIPubMedGoogle Scholar
  20. Feldgarden  M, Brover  V, Haft  DH, Prasad  AB, Slotta  DJ, Tolstoy  I, et al. Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates. Antimicrob Agents Chemother. 2019;63:e0048319. DOIPubMedGoogle Scholar
  21. Carattoli  A, Zankari  E, García-Fernández  A, Voldby Larsen  M, Lund  O, Villa  L, et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58:3895903. DOIPubMedGoogle Scholar
  22. Clausen  PTLC, Aarestrup  FM, Lund  O. Rapid and precise alignment of raw reads against redundant databases with KMA. BMC Bioinformatics. 2018;19:307. DOIPubMedGoogle Scholar
  23. Abyzov  A, Urban  AE, Snyder  M, Gerstein  M. CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing. Genome Res. 2011;21:97484. DOIPubMedGoogle Scholar
  24. Sharma  C, Nair  GB, Mukhopadhyay  AK, Bhattacharya  SK, Ghosh  RK, Ghosh  A. Molecular characterization of Vibrio cholerae O1 biotype El Tor strains isolated between 1992 and 1995 in Calcutta, India: evidence for the emergence of a new clone of the El Tor biotype. J Infect Dis. 1997;175:113441. DOIPubMedGoogle Scholar
  25. Wozniak  RAF, Fouts  DE, Spagnoletti  M, Colombo  MM, Ceccarelli  D, Garriss  G, et al. Comparative ICE genomics: insights into the evolution of the SXT/R391 family of ICEs. PLoS Genet. 2009;5:e1000786. DOIPubMedGoogle Scholar
  26. Lovero  KG, Mota-Bravo  L, Rasko  D. Closed genome sequence of an environmental Aeromonas veronii strain from California, United States, with an IncA/C plasmid carrying an extended-spectrum β-lactamase gene, blaVEB-3. Microbiol Resour Announc. 2022;11:e0103321. DOIPubMedGoogle Scholar
  27. Siriphap  A, Leekitcharoenphon  P, Kaas  RS, Theethakaew  C, Aarestrup  FM, Sutheinkul  O, et al. Characterization and genetic variation of Vibrio cholerae isolated from clinical and environmental sources in Thailand. PLoS One. 2017;12:e0169324. DOIPubMedGoogle Scholar
  28. Vijayalakshmi  N, Rao  RS, Badrinath  S. Minimum inhibitory concentration (MIC) of some antibiotics against Vibrio cholerae O139 isolates from Pondicherry. Epidemiol Infect. 1997;119:258. DOIPubMedGoogle Scholar
  29. Zhu  Y-G, Johnson  TA, Su  J-Q, Qiao  M, Guo  G-X, Stedtfeld  RD, et al. Diverse and abundant antibiotic resistance genes in Chinese swine farms. Proc Natl Acad Sci U S A. 2013;110:343540. DOIPubMedGoogle Scholar
  30. Mohanraj  RS, Samanta  P, Mukhopadhyay  AK, Mandal  J. Haitian-like genetic traits with creeping MIC of Azithromycin in Vibrio cholerae O1 isolates from Puducherry, India. J Med Microbiol. 2020;69:3728. DOIPubMedGoogle Scholar
  31. David  G. Antimicrobial drugs: chronicle of a twentieth century medical triumph. First edition. Oxford: Oxford University Press; 2008.
  32. Lamers  RP, Nguyen  UT, Nguyen  Y, Buensuceso  RN, Burrows  LL. Loss of membrane-bound lytic transglycosylases increases outer membrane permeability and β-lactam sensitivity in Pseudomonas aeruginosa. MicrobiologyOpen. 2015;4:87995. DOIPubMedGoogle Scholar
  33. Luo  Y, Wang  H, Liang  J, Qian  H, Ye  J, Chen  L, et al. Population structure and multidrug resistance of non-O1/non-O139 Vibrio cholerae in freshwater rivers in Zhejiang, China. Microb Ecol. 2021;82:31933. DOIPubMedGoogle Scholar




Cite This Article

DOI: 10.3201/eid2811.212066

Original Publication Date: October 13, 2022

1These authors contributed equally to this article.

Table of Contents – Volume 28, Number 11—November 2022

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Ruiting Lan, School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW 2052, Australia

Send To

10000 character(s) remaining.


Page created: September 28, 2022
Page updated: October 21, 2022
Page reviewed: October 21, 2022
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.