# Whole Genome Sequencing of Drug-Resistant Vibrio cholerae Serotype Ogawa from an Outbreak in Khyber Pakhtunkhwa

**Authors:** Aftab Ali, Momin Khan, Taj Ali Khan, Sajjad Ahmad, Noor Rahman, Aiman Waheed, Taane G. Clark

PMC · DOI: 10.3390/pathogens15010039 · 2025-12-29

## TL;DR

This study uses whole genome sequencing to analyze drug-resistant Vibrio cholerae from a cholera outbreak in Pakistan, revealing high antibiotic resistance and genetic links to international strains.

## Contribution

The study provides genomic insights into multidrug-resistant Vibrio cholerae isolates from a recent outbreak in Pakistan and identifies resistance determinants.

## Key findings

- V. cholerae isolates showed high resistance to ampicillin, sulfamethoxazole/trimethoprim, erythromycin, and tetracycline.
- Genomic analysis revealed multiple resistance determinants including mutations in gyrA, parC, and genes like tetA and blaSHV.
- The isolates were genetically related to both historical Pakistani strains and recent international isolates.

## Abstract

Background: Cholera, caused by Vibrio cholerae, remains endemic in many developing countries, including Pakistan. The extensive use of antibiotics has led to the emergence of antimicrobial resistance in V. cholerae, limiting available treatment options. In this study, we performed molecular characterisation of antibiotic-resistant V. cholerae serotype Ogawa isolates from a recent cholera outbreak in Khyber Pakhtunkhwa, Pakistan. Methodology: Suspected cholera stool samples were collected from hospitalised patients at various district hospitals of Khyber Pakhtunkhwa Province (KPK), Pakistan. The samples were transported to the Public Health Reference Microbiology Laboratory at Khyber Medical University, Peshawar. V. cholerae were identified based on colonial morphology, Gram staining, and biochemical tests using EPI 10E. For serotype identification, monovalent antisera were used. Antibiotic susceptibility testing (AST) was performed using CLSI M45 and EUCAST guidelines. DNA was extracted from pure colonies of multidrug-resistant (MDR) V. cholerae and subjected to whole-genome sequencing (WGS) for genomic characterisation using an Illumina MiSeq platform. Results: Of the 350 active diarrheal cases investigated, 70 were confirmed as V. cholerae. The outbreak was initially reported in Dir and was subsequently followed by a high incidence of cholera in the Peshawar district of KPK. All strains belong to the Ogawa serotype, which shows high antibiotic resistance, particularly to ampicillin (n = 62, 88.57%), Sulfamethoxazole/Trimethoprim (n = 60, 85.71%), Erythromycin (n = 59, 84.29%), and Tetracycline (n = 53, 75.71%). The lowest resistance was against Meropenem (n = 1, 1.4%), followed by amikacin (n = 7, 10.0%) and levofloxacin (n = 13, 18.57%). Furthermore, 34 (48.57%) of the isolates were MDR, while 13 (18.57%) were extensively drug-resistant. Six samples were selected for whole-genome sequencing. The selection of six V. cholerae samples for WGS was based on their drug resistance pattern and origin of isolation. At the genomic level, all sequenced V. cholerae strains harboured multiple antimicrobial resistance determinants. Quinolone resistance was associated with mutations and genes in gyrA, gyrB, parC, and parE; resistance to sulfamethoxazole–trimethoprim with folA, folP, and dfr; tetracycline resistance with tetA and tet35; chloramphenicol resistance with catB and S10p; and aminoglycoside resistance with hns, S12p, and gigB. In addition, β-lactam resistance was linked to the presence of efflux and β-lactamase genes, including blaSHV and mox-3. Mutations were identified in gyrA at positions S83I, S177A, and S202A, and in parC at positions S85L and I231V. Collectively, the presence of these resistance determinants likely enables V. cholerae to survive exposure to high concentrations of multiple antibiotics. Conclusions: Our V. cholerae isolates showed close genetic relatedness to previously sequenced strains from Pakistan (2010 and 2022), as well as to recently reported international strains from the USA, Australia, and China. These findings highlight both the long-term persistence of these lineages within Pakistan and their international dissemination, likely facilitated by globalisation.

## Linked entities

- **Genes:** GYRA (DNA GYRASE A) [NCBI Gene 820238], gyrB (DNA gyrase subunit B) [NCBI Gene 857440], CCL18 (C-C motif chemokine ligand 18) [NCBI Gene 6362], parE (DNA topoisomerase IV subunit B) [NCBI Gene 879897], folA (dihydrofolate reductase) [NCBI Gene 879233], folP (dihydropteroate synthase) [NCBI Gene 881715], DFR (dihydroflavonol 4-reductase) [NCBI Gene 544150], tet(A) (tetracycline efflux MFS transporter Tet(A)) [NCBI Gene 33941499], tet(35) (tetracycline efflux Na+/H+ antiporter family transporter Tet(35)) [NCBI Gene 1188639], TYRP1 (tyrosinase related protein 1) [NCBI Gene 7306], hns (histone-like protein Hns) [NCBI Gene 886187], gigB (anti-anti-sigma factor GigB) [NCBI Gene 9383709], bla SHV (class A extended-spectrum beta-lactamase SHV-2) [NCBI Gene 40101717], mox3 (putative monooxygenase) [NCBI Gene 18813168]
- **Chemicals:** ampicillin (PubChem CID 6249), Sulfamethoxazole/Trimethoprim (PubChem CID 358641), Erythromycin (PubChem CID 12560), Tetracycline (PubChem CID 54675776), Meropenem (PubChem CID 441130), amikacin (PubChem CID 37768), levofloxacin (PubChem CID 149096)
- **Diseases:** cholera (MONDO:0015766)
- **Species:** Vibrio cholerae (taxon 666), Mus musculus (taxon 10090)

## Full-text entities

- **Diseases:** diarrheal (MESH:D004403), Cholera (MESH:D002771)
- **Chemicals:** levofloxacin (MESH:D064704), Sulfamethoxazole (MESH:D013420), Quinolone (MESH:D015363), Tetracycline (MESH:D013752), amikacin (MESH:D000583), sulfamethoxazole-trimethoprim (MESH:D015662), beta-lactam (MESH:D047090), aminoglycoside (MESH:D000617), Trimethoprim (MESH:D014295), chloramphenicol (MESH:D002701), Meropenem (MESH:D000077731), ampicillin (MESH:D000667), S10p (-), Erythromycin (MESH:D004917)
- **Species:** Vibrio cholerae (species) [taxon 666], Homo sapiens (human, species) [taxon 9606]
- **Mutations:** S177A, S83I, I231V, S85L, S202A

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12845007/full.md

---
Source: https://tomesphere.com/paper/PMC12845007