- Genome Report
- Open Access
Prevailing clone (ST69) of Vibrio cholerae O139 in India over 10 years
Gut Pathogens volume 9, Article number: 60 (2017)
Vibrio cholerae is responsible for the cause of severe life-threatening infection known as cholera. The study aimed to analyze the genetic make-up of V. cholerae O139 isolates from India and compare its phylogeny with the global strains. The genome data revealed that all isolates were of same sequence type (ST69) which belongs to seventh pandemic clone, with same virulence gene profile and, antimicrobial resistance gene profile except for two isolates. No known CRISPR repeats were identified in any of these isolates. Three different phages were identified among the isolates. All the isolates were found to harbour int SXT and seventh pandemic-specific gene (VC2346). Results from this study enhance our understanding on the persistence of ST69 V. cholerae O139 over 20 years.
Historically, only the O1 serogroup of V. cholerae has been the causative agent of epidemic cholera. In late 1992, in India and Bangladesh, a novel serogroup designated V. cholerae O139 emerged, replacing V. cholerae O1, and gave rise to major cholera outbreaks . Based on the whole genome phylogenies of V. cholerae, the relation of isolates responsible for global cholera could be understood . Nearly 141 draft genomes of V. cholerae O139 are accessible in the NCBI genome projects. Evolution and adaptation mechanism of this pathogen can be demonstrated based on genomic diversity of clinical isolates from various geographical regions and timescale . The current study aims to perform and report the whole genome shotgun sequences of V. cholerae O139 isolates for the first time in India. The genomes reported here will help to better understand the evolution and endemicity of cholera in India.
Clinical strains of V. cholerae O139 obtained from the past 20 years (FC1105-2005, FC1225-2001, FC1341-2002, FC1384-2000, FC1817-1994, FC1877-1995, FC2271-1997, FC2273-1998, FC3611a-1999 and FC3611b-1997) were revived from the repository at the Department of Clinical Microbiology, Christian Medical College, Vellore, India. Isolates were identified with standard biochemical methods and were confirmed by agglutination with O139 antisera raised in-house. The genomic DNA was extracted using Qiagen automated DNA extraction method with QIAsymphony DSP Virus/Pathogen Mini Kit (QIAsymphony, Qiagen, Germany). All 10 V. cholerae O139 was confirmed based on the 16S rRNA gene, by amplification, sequencing, and BLAST against the NCBI database.
Next generation sequencing
The whole genome shotgun sequencing for the V. cholerae O139 isolates was performed using Ion Torrent (PGM, Life Technologies) with 400-bp read chemistry (Life Technologies, Carlsbad, CA). Genomic DNA libraries were prepared using Ion Plus Fragment Library Kit; (Life Technologies) according to manufacturer’s instructions. AMPure beads were used to purify the genomic libraries and their concentrations were determined using the Qubit 3.0 fluorimeter (Invitrogen, Merelbeke, Belgium). Emulsion PCR was performed on pooled libraries (Ion One Touch Hi-Q 400 Template Kit v2 DL Kit; Life Technologies), and template-positive Ion Sphere particles were enriched using Dynabeads Myone streptavidin C1 beads. Finally, pooled samples were loaded on Ion 318V2 chip for sequencing.
Genome assembly and annotation
The generated whole genome data were assembled de novo using Assembler SPAdes v.220.127.116.11 embedded in Torrent suite server v.5.0.5. The genome sequence was annotated using PATRIC (Pathosystems Resource Integration Center), the bacterial bioinformatics database and analysis resource (http://www.patricbrc.org) , and the NCBI Prokaryotic Genome Automatic Annotation Pipeline (PGAAP) (http://www.ncbi.nlm.nih.gov/genomes/static/Pipeline.html).
Downstream genome analysis
The whole genome data were analyzed using open access tools at Center for Genomic Epidemiology (CGE) web-based server. Sequence types for the study isolates were determined using MLST 1.8 (https://cge.cbs.dtu.dk//services/MLST/), and antimicrobial resistance were identified using ResFinder 2.1 (https://cge.cbs.dtu.dk//services/ResFinder/) with the 90% threshold for identity and with 60% of minimum length coverage. Genome data was screened for virulence genes through PATRIC database. The presence of plasmids was analyzed using PlasmidFinder 1.3 (https://cge.cbs.dtu.dk//services/PlasmidFinder/) with 95% threshold for identity. The study isolates were further screened for the presence of prophage sequences within the genome using PHAST (PHAge Search Tool) (http://phast.wishartlab.com) .
The web-based MyDbFinder 1.0 (https://cge.cbs.dtu.dk/services/MyDbFinder/) was used, in silico to determine the SXT elements (int SXT), specific integrase genes of class 1 integron (intI gene) and specific gene (VC2346) of the seventh pandemic strains with a selected threshold equal to 98% identity as previously described .
The phylogenetic analysis was performed based on MLST housekeeping genes (adk, gyrB, mdh, metE, pntA, purM, pyrC) using the MEGA7 software . The evolutionary history was inferred using the Neighbor-Joining method. The evolutionary distances were computed using the Maximum Composite Likelihood method. goEBURST analysis for the V. cholerae O139 strains was performed using PHYLOViZ 1.1 software  to identify the relation between study isolates and the global strains with sequence types of V. cholerae (ST168-KC895170.1, ST169-KC895206.1, ST170-KJ657169.1, ST182-KC895159.1, ST429-KC895162.1, N16961-AE003852.1, M2140_ST70-CP013315.1, M802_ST71-LT907989.1, M1616_ST73-KC895170.1, M796_ST75-DQ316973.1).
Results and discussion
Genome features, virulence and antimicrobial resistance determinants
Assembly of the raw reads of the isolates ranged from 69 to 196X coverage (Table 1). NGS revealed that all isolates belong to same sequence type ST69. Similarly in a study by Siriphap et al. , most of the seventh pandemic isolates belonged to ST69 and all the clinical strains of O1 and O139 were highly conserved to ST69. Literature study reveals that the Bay of Bengal as a major hub connecting the spread of cholera around the globe and has continually being evolved in South Asia, later spread around the globe in at least three independent but overlapping waves. Besides V. cholerae serogroup O1 is widespread throughout the globe, whereas serogroup O139 is almost solely confined to Asia .The probable transmission of V. cholerae O139 sequence types (ST69, ST187, ST161, ST162 and ST163) that reported elsewhere from South Asia was depicted in Fig. 1.
ResFinder 2.1 returned same antimicrobial resistance genes for all except two isolates (Table 1). catB9 gene responsible for chloramphenicol resistance was found common in all isolates. strA and strB (streptomycin resistance), floR (florfenicol/chloramphenicol resistance), and sul2 (sulphonamide resistance) genes were identified in all except FC1225 and FC1341 strains. The occurrence of these resistance genes were previously been reported in V. cholerae O139 serogroup. Number of antimicrobial resistance genes identified in these genomes as per comprehensive antibiotic resistance database, CARD and antibiotic resistance genes database, ARDB database matching were ~ 17 and ~ 6 genes. PlasmidFinder 1.3 resulted negative for all 10 isolates. Additionally, all 10 isolates were found to possess virulence genes including ctxA, ctxB, rtx, hlyA, zot, ace, tcpT, higB, doc, parE.
Repeat regions and insertional elements
The clusters of regularly interspaced short palindromic repeats (CRISPR) and spacer sequences in the genome were screened using CRISPR finder (http://crispr.u-psud.fr/Server/)  and PATRIC database which resulted in no CRISPR regions in these genomes (Table 1).
However, ISfinder (https://www-is.biotoul.fr/blast.php)  resulted various IS elements that were listed in Table 1. Most of the IS elements such as IS3-ISVch4, IS200-IS1004, ISAs1-IS1358 (O139), IS481, Tn3-ISShfr9, IS66-ISPsy43, IS5-ISVch8, IS5-ISVch5, IS3-ISVch4, Tn3-ISPsy42 were previously reported from V. cholerae genomes. Of these, ISAs1-IS1358 was reported to be from V. cholerae O139. Prophage screening results showed that the isolates FC1105, FC1225, FC1341, FC1877, FC3611a and FC3611b were found to have V. cholerae phage CTX (NC_015209). The phage K139 (NC_003313) was detected in FC1817, whereas VFJ (NC_021562) was seen in FC2271 and FC2273 respectively.
All study isolates lacked the intI gene but were possessed the SXT element (int SXT) with 99.9% identity to the reference sequence, which was mostly known to be, associated with the sulfamethoxazole and trimethoprim resistant strains . All isolates harboured the seventh pandemic-specific gene (VC2346), suggesting that they belong to the same clonal linage and might originate from a common ancestor of the seventh pandemic strain.
The phylogenetic analysis based on MLST housekeeping genes revealed the relation between ST69 Indian strains and global strains (Fig. 2). All the study isolates are closely related to ST429 originated from M802 (ST71) global isolates. Also, M2140 (ST70), M796 (ST75), ST168, M1616 (ST73), ST169, ST170 and ST182 belong to the same clonal complex with ST71 as the founder type. goEBURST analysis of V. cholerae confirmed that ST69 belong to the clonal complex CC0 with the above-mentioned sequence types (Fig. 3).
The genome sequences of V. cholerae O139 strains have been deposited at DDBJ/ENA/Genbank under the accession numbers NGQQ00000000–NGQZ00000000. The version described in this manuscript is version one.
To the best of our knowledge, this is the first report of ST69 V. cholerae O139 in India. Though there are few reports of clinical V. cholerae isolates from India, our study adds new dimension to the earlier finding to show that the sequence type of V. cholerae O139 continues to persist in India for over 20 years is of ST69. This implies that clonal expansion of V. cholerae O1 seventh pandemic has resulted in ST69 strains. The isolates exhibited same virulence profile but different antimicrobial resistance gene profile. This study helps in understanding the epidemiology of V. cholerae and genomic relationship between these isolates with the strains from China and Thailand. Also a continuous surveillance of V. cholerae O139 serogroup can provide more insights into the evolution of V. cholerae.
Jesudason MV, John TJ. The Vellore vibrio watch. Lancet. 1996;347:1493–4.
Mutreja A, Kim DW, Thomson NR, Connor TR, Lee JH, Kariuki S, et al. Evidence for several waves of global transmission in the seventh cholera pandemic. Nature. 2011;477:462–5.
Yi Y, Lu N, Liu F, Li J, Zhang R, Jia L. Genome sequence and comparative analysis of a Vibrio cholerae O139 strain E306 isolated from a cholera case in China. Gut Pathog. 2014;6:3.
Wattam AR, Abraham D, Dalay O, Disz TL, Driscoll T, Gabbard JL, et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Res. 2014;42:D581–91.
Zhou Y, Liang Y, Lynch K, Dennis JJ, Wishart DS. PHAST: a fast phage search too. Nucleic Acids Res. 2011;39:W347–52.
Siriphap A, Leekitcharoenphon P, Kaas RS, Theethakaew C, Aarestrup FM, Sutheinkul O, Hendriksen RS. Characterization and genetic variation of Vibrio cholerae isolated from clinical and environmental sources in Thailand. PLoS ONE. 2017;12:e0169324.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9.
Francisco AP, Vaz C, Monteiro PT, Melo-Cristino J, Ramirez M, Carriço JA. PHYLOViZ: phylogenetic inference and data visualization for sequence based typing methods. BMC Bioinform. 2012;13:87.
Moore S, Thomson N, Mutreja A, Piarroux R. Widespread epidemic cholera caused by a restricted subset of Vibrio cholerae clones. Clin Microbiol Infect. 2014;20:373–9.
Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:W52–7.
Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006;34:D32–6.
DRNK, MSDP, SA and BV designed the study. MSDP and TS performed the laboratory work. DRNK and MSDP analysed the data and wrote the manuscript. All authors read and approved the final manuscript.
The authors gratefully acknowledge the Institutional Review Board of the Christian Medical College, Vellore (83-i/11/13) for approving the study and providing lab space and facilities.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
Ethics approval and consent to participate
The study was funded by Indian Council of Medical Research (ICMR), New Delhi (Ref. No: AMR/TF/55/13ECDII dated 23/10/2013).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.