# Genetic approach toward linkage of Iran 2012–2016 cholera outbreaks with 7th pandemic Vibrio cholerae

Fatemeh Jalalizadeh, Elisabeth Njamkepo, François-Xavier Weill, Forough Goodarzi, Marjan Rahnamaye-Farzami, Roghieh Sabourian, Bita Bakhshi

► To cite this version:

Fatemeh Jalalizadeh, Elisabeth Njamkepo, François-Xavier Weill, Forough Goodarzi, Marjan Rahnamaye-Farzami, et al.. Genetic approach toward linkage of Iran 2012-2016 cholera outbreaks with 7th pandemic Vibrio cholerae. BMC Microbiology, 2024, 24, pp.33. 10.1186/s12866-024-03185-9. hal-04613488

HAL Id: hal-04613488

https://hal.science/hal-04613488

Submitted on 17 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Genetic approach toward linkage of Iran 2012–2016 cholera outbreaks with 7th pandemic Vibrio cholerae

Fatemeh Jalalizadeh, Elisabeth Njamkepo, François-Xavier Weill, Forough Goodarzi, Marjan Rahnamaye-Farzami, Roghieh Sabourian, Bita Bakhshi

► To cite this version:

Fatemeh Jalalizadeh, Elisabeth Njamkepo, François-Xavier Weill, Forough Goodarzi, Marjan Rahnamaye-Farzami, et al.. Genetic approach toward linkage of Iran 2012-2016 cholera outbreaks with 7th pandemic Vibrio cholerae. BMC Microbiology, 2024, 24, pp.33. 10.1186/s12866-024-03185-9. hal-04613488

HAL Id: hal-04613488

https://hal.science/hal-04613488v1

Submitted on 17 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

# Genetic approach toward linkage of Iran 2012–2016 cholera outbreaks with 7th pandemic *Vibrio cholerae*

Fatemeh Jalalizadeh¹, Elisabeth Njamkepo², François-Xavier Weill², Forough Goodarzi¹, Marjan Rahnamaye-Farzami³, Roghieh Sabourian³ and Bita Bakhshi¹*

## Abstract

*Vibrio cholerae*, as a natural inhabitant of the marine environment is among the world-leading causes of diarrheal diseases. The present study aimed to investigate the genetic relatedness of Iran 2012–2016 *V. cholerae* outbreaks with 7th pandemic cholera and to further characterize the non-ST69/non-ST75 sequence types strains by whole-genome sequencing (WGS).

Twenty *V. cholerae* isolates related to 2012, 2013, 2015 and 2016 cholera outbreaks were studied by two genotyping methods – Pulsed-field Gel Electrophoresis (PFGE) and Multi-locus Sequence Typing (MLST)-and by antimicrobial susceptibility testing. Seven sequence types (STs) and sixteen pulsotypes were detected. Sequence type 69 was the most abundant ST confirming that most (65%, 13/20) of the studied isolates collected in Iran between 2012 and 2016 belonged to the 7th pandemic clone. All these ST69 isolates (except two) exhibited similar pulsotypes. ST75 was the second most abundant ST. It was identified in 2015 and 2016. ST438, ST178, ST579 and STs of 983 and 984 (as newfound STs) each were only detected in one isolate. All strains collected in 2016 appeared as distinct STs and pulsotypes indicative of probable different originations. All ST69 strains were resistant to nalidixic acid. Moreover, resistance to nalidixic acid, trimethoprim-sulfamethoxazole and tetracycline was only observed in strains of ST69. These properties propose the ST69 as a unique genotype derived from a separate lineage with distinct resistance properties. The circulation of *V. cholerae* ST69 and its traits in recent years in Iran proposes the 7th pandemic strains as the ongoing causes of cholera outbreaks in this country, although the role of ST75 as the probable upcoming dominant ST should not be ignored.

Genomic analysis of non-ST69/non-ST75 strains in this study showed ST579 is the most similar ST type to 7th pandemic sequence types, due to the presence of wild type-El Tor sequences of tcpA and VC-1319, VC-1320, VC-1577, VC-1578 genes (responsible for polymyxin resistance in El Tor biotype), the traits of rstC of RS1 phage in one strain of this ST type and the presence of VPI-1 and VSP-I islands in ST579 and ST178 strains. In silico analysis showed no significant presence of resistance genes/cassettes/plasmids within non-ST69/non-ST75 strains genomes. Overall, these data indicate the higher susceptibility of *V. cholerae* non-ST69/non-ST75 strains in comparison with more ubiquitous and more circulating ST69 and ST75 strains.

In conclusion, the occurrence of small outbreaks and sporadic cholera cases due to *V. cholerae* ST69 in recent years in Iran shows the 7th pandemic strains as the persistent causes of cholera outbreaks in this country, although the role

<footnote>*Correspondence:
Bita Bakhshi
b.bakhshi@modares.ac.ir
Full list of author information is available at the end of the article</footnote>
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

# Genetic approach toward linkage of Iran 2012–2016 cholera outbreaks with 7th pandemic *Vibrio cholerae*

Fatemeh Jalalizadeh¹, Elisabeth Njamkepo², François-Xavier Weill², Forough Goodarzi¹, Marjan Rahnamaye-Farzami³, Roghieh Sabourian³ and Bita Bakhshi¹*

## Abstract

*Vibrio cholerae*, as a natural inhabitant of the marine environment is among the world-leading causes of diarrheal diseases. The present study aimed to investigate the genetic relatedness of Iran 2012–2016 *V. cholerae* outbreaks with 7th pandemic cholera and to further characterize the non-ST69/non-ST75 sequence types strains by whole-genome sequencing (WGS).

Twenty *V. cholerae* isolates related to 2012, 2013, 2015 and 2016 cholera outbreaks were studied by two genotyping methods – Pulsed-field Gel Electrophoresis (PFGE) and Multi-locus Sequence Typing (MLST)-and by antimicrobial susceptibility testing. Seven sequence types (STs) and sixteen pulsotypes were detected. Sequence type 69 was the most abundant ST confirming that most (65%, 13/20) of the studied isolates collected in Iran between 2012 and 2016 belonged to the 7th pandemic clone. All these ST69 isolates (except two) exhibited similar pulsotypes. ST75 was the second most abundant ST. It was identified in 2015 and 2016. ST438, ST178, ST579 and STs of 983 and 984 (as newfound STs) each were only detected in one isolate. All strains collected in 2016 appeared as distinct STs and pulsotypes indicative of probable different originations. All ST69 strains were resistant to nalidixic acid. Moreover, resistance to nalidixic acid, trimethoprim-sulfamethoxazole and tetracycline was only observed in strains of ST69. These properties propose the ST69 as a unique genotype derived from a separate lineage with distinct resistance properties. The circulation of *V. cholerae* ST69 and its traits in recent years in Iran proposes the 7th pandemic strains as the ongoing causes of cholera outbreaks in this country, although the role of ST75 as the probable upcoming dominant ST should not be ignored.

Genomic analysis of non-ST69/non-ST75 strains in this study showed ST579 is the most similar ST type to 7th pandemic sequence types, due to the presence of wild type-El Tor sequences of tcpA and VC-1319, VC-1320, VC-1577, VC-1578 genes (responsible for polymyxin resistance in El Tor biotype), the traits of rstC of RS1 phage in one strain of this ST type and the presence of VPI-1 and VSP-I islands in ST579 and ST178 strains. In silico analysis showed no significant presence of resistance genes/cassettes/plasmids within non-ST69/non-ST75 strains genomes. Overall, these data indicate the higher susceptibility of *V. cholerae* non-ST69/non-ST75 strains in comparison with more ubiquitous and more circulating ST69 and ST75 strains.

In conclusion, the occurrence of small outbreaks and sporadic cholera cases due to *V. cholerae* ST69 in recent years in Iran shows the 7th pandemic strains as the persistent causes of cholera outbreaks in this country, although the role

<footnote>*Correspondence:
Bita Bakhshi
b.bakhshi@modares.ac.ir
Full list of author information is available at the end of the article</footnote>
of ST75 as the second most contributed ST should not be ignored. The occurrence of non-ST69/non-ST75 sequence types with some virulence factors characteristics in border provinces in recent years is noteworthy, and further studies together with surveillance efforts are expected to determine their likely route of transport.

**Keywords** Vibrio cholerae, PFGE, MLST, Cholera, Iran, WGS

**Introduction**

*Vibrio cholerae*, as a bacterial water and foodborne pathogen, causes significant morbidity and mortality worldwide. Cholera symptoms include dehydrating diarrhea, vomiting and abdominal cramps [1, 2]. *V. cholerae* is responsible for seven important pandemics of cholera since 1817 [3]. According to a World Health Organization (WHO) report, it has been evaluated that millions of cholera cases and tens of thousands of mortalities happen each year around the world [4]. Epidemiological and source tracking investigation of cholera outbreaks in human societies often requires typing of strains to identify the prevalence and relatedness of epidemic and endemic isolates. All these events emphasize on typing investigations to determine the source of outbreak and improve public health, and monitoring vaccination programs [5-7].

Different typing methods have been developed for outbreak surveillance of *V. cholerae* strains. The conventional phenotypic epidemiological tools including serotyping, phage typing and antibiogram are inadequate to provide standardized description of clonal definition of *V. cholerae* isolates involved in outbreaks. Different genotyping tools have also been developed for subspecies-level classification of *V. cholerae*. These methods include serotyping, multilocus sequence typing (MLST) [8], multilocus variable number of tandem repeats (VNTR) analysis (MLVA) [9], single-nucleotide polymorphism (SNP)- based approaches and Core genome MLST (cgMLST) in dealing with large data sets [10]. Despite some important limitations to each, these methods are progressively used individually or in combination, to provide the best approach to bacterial genetic relatedness.

Among genotyping methods, Pulsed-field Gel Electrophoresis (PFGE) has been selected as the gold standard in molecular epidemiology of small-scale outbreaks [11]. The power of discrimination and reproducibility of PFGE has made it a widely useful method for bacterial species genotyping [12]. Some limitations of PFGE include identical band patterns are not always guaranteeing homologous genetic content and the inability of gel-based methods to resolve similar sized bands which could be resolved in combination with other genotyping methods.

Multi-locus Sequence Typing (MLST) as a sequence-based method, investigates variations in seven housekeeping genes and provides information on population

structure and evolution of bacterial species as well as their epidemiological significance. Furthermore, MLST provides comparable, portable and reproducible genotypic data which can be readily exchanged among laboratories [13-16].

Core genome MLST (cgMLST) which relies on whole-genome sequencing (WGS), although providing a superior scheme and more phylogenetic resolution than traditional MLST [10], there is a need to continue with MLST amplicon sequencing to allow some level of genomic epidemiology where access to second or third generation sequencing is limited, especially in dealing with large data sets in low-income and middle-income countries.

In global surveillance of pandemic-causing pathogens, such as *V. cholerae*, it is important to use an efficient and easy to use genotyping tool, with the potential to be applied to all *V. cholerae* strains around the world.

Combining PFGE and MLST analysis can provide higher resolution in characterization of *V. cholerae* isolates and their relationship to each other. In the present study, we investigated the genetic relatedness of *V. cholerae* strains collected from 2012-2016 cholera outbreaks in Iran using PFGE and MLST analysis. Moreover, the non-ST69/non-ST75 sequence types strains were further characterized using whole-genome sequencing (WGS).

**Materials and methods**

**Ethics statement**

This study protocol was approved by the Ethics Committee of Tarbiat Modares University (Code: IR.TMU.REC.1397.092). All patients completed a questionnaire and gave informed consent prior to the study.

**Bacterial strains**

We studied 20 *V. cholerae* strains from stool sample of cholera patients in Iran in 2012 (n=5), 2013 (n=5), 2015 (n=6) and 2016 (n=4). These *V. cholerae* strains were selected from our collection containing 11 strains isolated in 2012 (~20% of total cholera cases in Iran this year), 33 in 2013 (~13% of total cholera cases in Iran this year), 6 in 2015 (~7% of total cholera cases in Iran this year) and 4 in 2016 (no data on total cholera cases). Due to high level of clonality among strains of 2012 and 2013 [17, 18], five representatives from each year were randomly selected and studied.
of ST75 as the second most contributed ST should not be ignored. The occurrence of non-ST69/non-ST75 sequence types with some virulence factors characteristics in border provinces in recent years is noteworthy, and further studies together with surveillance efforts are expected to determine their likely route of transport.

**Keywords** Vibrio cholerae, PFGE, MLST, Cholera, Iran, WGS

**Introduction**

*Vibrio cholerae*, as a bacterial water and foodborne pathogen, causes significant morbidity and mortality worldwide. Cholera symptoms include dehydrating diarrhea, vomiting and abdominal cramps [1, 2]. *V. cholerae* is responsible for seven important pandemics of cholera since 1817 [3]. According to a World Health Organization (WHO) report, it has been evaluated that millions of cholera cases and tens of thousands of mortalities happen each year around the world [4]. Epidemiological and source tracking investigation of cholera outbreaks in human societies often requires typing of strains to identify the prevalence and relatedness of epidemic and endemic isolates. All these events emphasize on typing investigations to determine the source of outbreak and improve public health, and monitoring vaccination programs [5-7].

Different typing methods have been developed for outbreak surveillance of *V. cholerae* strains. The conventional phenotypic epidemiological tools including serotyping, phage typing and antibiogram are inadequate to provide standardized description of clonal definition of *V. cholerae* isolates involved in outbreaks. Different genotyping tools have also been developed for subspecies-level classification of *V. cholerae*. These methods include serotyping, multilocus sequence typing (MLST) [8], multilocus variable number of tandem repeats (VNTR) analysis (MLVA) [9], single-nucleotide polymorphism (SNP)- based approaches and Core genome MLST (cgMLST) in dealing with large data sets [10]. Despite some important limitations to each, these methods are progressively used individually or in combination, to provide the best approach to bacterial genetic relatedness.

Among genotyping methods, Pulsed-field Gel Electrophoresis (PFGE) has been selected as the gold standard in molecular epidemiology of small-scale outbreaks [11]. The power of discrimination and reproducibility of PFGE has made it a widely useful method for bacterial species genotyping [12]. Some limitations of PFGE include identical band patterns are not always guaranteeing homogeneous genetic content and the inability of gel-based methods to resolve similar sized bands which could be resolved in combination with other genotyping methods.

Multi-locus Sequence Typing (MLST) as a sequence-based method, investigates variations in seven housekeeping genes and provides information on population

structure and evolution of bacterial species as well as their epidemiological significance. Furthermore, MLST provides comparable, portable and reproducible genotypic data which can be readily exchanged among laboratories [13-16].

Core genome MLST (cgMLST) which relies on whole-genome sequencing (WGS), although providing a superior scheme and more phylogenetic resolution than traditional MLST [10], there is a need to continue with MLST amplicon sequencing to allow some level of genomic epidemiology where access to second or third generation sequencing is limited, especially in dealing with large data sets in low-income and middle-income countries.

In global surveillance of pandemic-causing pathogens, such as *V. cholerae*, it is important to use an efficient and easy to use genotyping tool, with the potential to be applied to all *V. cholerae* strains around the world.

Combining PFGE and MLST analysis can provide higher resolution in characterization of *V. cholerae* isolates and their relationship to each other. In the present study, we investigated the genetic relatedness of *V. cholerae* strains collected from 2012-2016 cholera outbreaks in Iran using PFGE and MLST analysis. Moreover, the non-ST69/non-ST75 sequence types strains were further characterized using whole-genome sequencing (WGS).

**Materials and methods**

**Ethics statement**

This study protocol was approved by the Ethics Committee of Tarbiat Modares University (Code: IR.TMU.REC.1397.092). All patients completed a questionnaire and gave informed consent prior to the study.

**Bacterial strains**

We studied 20 *V. cholerae* strains from stool sample of cholera patients in Iran in 2012 (n=5), 2013 (n=5), 2015 (n=6) and 2016 (n=4). These *V. cholerae* strains were selected from our collection containing 11 strains isolated in 2012 (~20% of total cholera cases in Iran this year), 33 in 2013 (~13% of total cholera cases in Iran this year), 6 in 2015 (~7% of total cholera cases in Iran this year) and 4 in 2016 (no data on total cholera cases). Due to high level of clonality among strains of 2012 and 2013 [17, 18], five representatives from each year were randomly selected and studied.
The isolates were confirmed by biochemical and molecular methods. After enrichment by alkaline peptone water, yellow colonies on thiosulfate citrate bile salts sucrose (TCBS agar) were isolated and tested for oxidize, Kligler iron slant agar reaction, arginine and esculin hydrolysis, Lysine and ornithine decarboxylase activity and growth in 0% NaCl [19, 20]. All twenty strains were characterized by polyclonal O1 and Ogawa and Inaba specific antisera (Bahar Afshan, Iran). Biotyping of isolates was performed by the Voges-Proskauer (VP) test and the sheep red blood cells Hemagglutination assay [21]. Molecular identification of *V. cholerae* isolates was performed by PCR with primers for targeting the 16S-23S rRNA intergenic region as previously described [22].

## Antimicrobial susceptibility testing

Antibiotic susceptibility was tested by disk diffusion method according to Clinical and Laboratory Standards Institute (CLSI-M45) guidelines [23, 24] with following antibiotics: ampicillin (10 µg), cefotaxime (30 µg), cefotazidime (30 µg), gentamicin (10 µg), tetracycline (30 µg), ciprofloxacin (5 µg), trimethoprim-sulfamethoxazole (1.23/23.75 µg), chloramphenicol (30 µg), azithromycin, erythromycin and nalidixic acid. All antimicrobial disks were purchased from Mast Group (Mast Group Ltd, UK). The minimum inhibitory concentrations (MICs) of ampicillin were determined by the agar dilution method. The MIC E-test strip (Liofilchem, Italy) was used for determination of MIC of azithromycin and erythromycin for all strains.

## Pulsed-field gel electrophoresis

PFGE was performed on 10 isolates from 2015 and 2016, based on the Standard Operating Procedure for PulseNet PFGE of *Vibrio cholerae* and *Vibrio parahaemolyticus* [25]. Briefly, Proteinase K (Bioneer, South Korea) and NotI restriction enzyme (Thermo Fisher Scientific, US) was used for digestion of the genomic DNA of bacterial cells. PFGE was performed with a CHEF Mapper system (BioRad) with a two-block program under the following conditions: The block 1 ranged from 2–10 s for 13 h and block 2 from 20–25 s for 6 h at 6.0 V/cm. *Salmonella* serotype Branderup H9812 was utilized as DNA molecular mass size marker. PFGE pattern of isolates 2012 and 2013 were available from our previous researches [17, 18]. PFGE banding patterns were analyzed by BioNumerics software (version 7.6.3; Applied Maths, Belgium). The dendrogram was generated using Dice similarity index and UPGMA method. The strains with more than > 30% similarity fell into major clusters, then types and subtypes were assigned according to criteria presented by Tenover et al. [26].

## Multi-locus Sequence Typing (MLST)

Genomic DNA was extracted with AccuPrep genomic extraction kit (Bioneer, Korea) for all amplification tests. Seven housekeeping genes (*adk*, *gyrB*, *metE*, *mdh*, *pntA*, *purM*, *pyrC*) were used for MLST analysis [8]. Amplification of genes was performed by primers that are shown in Table 1. Each PCR was performed with an initial denaturation at 95°C for 3 min, denaturing at 95°C for 30 s, extension at 72°C for 1 min and a final extension at 72°C for 2 min repeated for 34 cycles. Annealing temperatures were adjusted according to melting temperature of each primer sets. PCR reactions were carried out in 25 µl total reaction volume containing 8 µl PCR-Master mix (Taq DNA Polymerase 2×Master Mix RED, 1.5 mM MgCl2) (Ampliqon, Denmark), 2 µl template DNA, 1 µl of each of the primers (10 pmol), and 13 µl sterilized distilled water. The DNA purification and sanger sequencing were performed by Microsynth (Switzerland). Sequences were trimmed, placed in the right direction and aligned at the same start and stop positions by CLC Genomic Workbench (version 12.0; Qiagen) and BioNumerics softwares (version 7.6.3; Applied Maths, Belgium). The alleles were uploaded to PubMLST and assigning the alleles and sequence type of strains were defined according to PubMLST website (https://pubmlst.org/organisms/vibrio-cholerae).

## Real-time PCR for toxin gene detection

Real-time PCR was used for detection of cholera toxin (*ctxA/ctxB* genes) and toxin co-regulated pilus (*tcpA* gene). Primer sequences are described in Table 1. The SYBR green Premix EX Taq mixture (Takara, Japan) was in a total volume of 20µL containing 0.8µL of each primer pair, 2µL of DNA sample, and 6.4µL of distilled water. RT- PCR cycling conditions were as follows: 95°C for 5 min and 40 cycles of 95°C for 15s, 60°C for 43 s.

## Whole-genome sequencing (WGS) and analysis of selected strains

Genomic DNA was extracted using AccuPrep genomic extraction kit (Bioneer, Korea). Whole genome sequencing was performed using the NovaSeq 6000 and HiSeq 2500 instruments (Illumina, San Diego, California, USA). DNA libraries were prepared using a Nextera DNA Flex Library Preparation Kit (Illumina), followed by 2×150 bp paired-end sequencing runs with a median coverage of 171 X (95–628). Short-read sequence data were submitted in the European Nucleotide Archive (ENA) (http://www.ebi.ac.uk/ena) and the accession numbers are listed in Table 2.

Short reads were first cleaned using FqCleanER version 3.0 (https://gitlab.pasteur.fr/GIPhy/fqCleanER)
The isolates were confirmed by biochemical and molecular methods. After enrichment by alkaline peptone water, yellow colonies on thiosulfate citrate bile salts sucrose (TCBS agar) were isolated and tested for oxidize, Kligler iron slant agar reaction, arginine and esculin hydrolysis, Lysine and ornithine decarboxylase activity and growth in 0% NaCl [19, 20]. All twenty strains were characterized by polyclonal O1 and Ogawa and Inaba specific antisera (Bahar Afshan, Iran). Biotyping of isolates was performed by the Voges-Proskauer (VP) test and the sheep red blood cells Hemagglutination assay [21]. Molecular identification of *V. cholerae* isolates was performed by PCR with primers for targeting the 16S-23S rRNA intergenic region as previously described [22].

## Antimicrobial susceptibility testing

Antibiotic susceptibility was tested by disk diffusion method according to Clinical and Laboratory Standards Institute (CLSI-M45) guidelines [23, 24] with following antibiotics: ampicillin (10 µg), cefotaxime (30 µg), cefotazidime (30 µg), gentamicin (10 µg), tetracycline (30 µg), ciprofloxacin (5 µg), trimethoprim-sulfamethoxazole (1.23/23.75 µg), chloramphenicol (30 µg), azithromycin, erythromycin and nalidixic acid. All antimicrobial disks were purchased from Mast Group (Mast Group Ltd, UK). The minimum inhibitory concentrations (MICs) of ampicillin were determined by the agar dilution method. The MIC E-test strip (Liofilchem, Italy) was used for determination of MIC of azithromycin and erythromycin for all strains.

## Pulsed-field gel electrophoresis

PFGE was performed on 10 isolates from 2015 and 2016, based on the Standard Operating Procedure for PulseNet PFGE of *Vibrio cholerae* and *Vibrio parahaemolyticus* [25]. Briefly, Proteinase K (Bioneer, South Korea) and NotI restriction enzyme (Thermo Fisher Scientific, US) was used for digestion of the genomic DNA of bacterial cells. PFGE was performed with a CHEF Mapper system (BioRad) with a two-block program under the following conditions: The block 1 ranged from 2–10 s for 13 h and block 2 from 20–25 s for 6 h at 6.0 V/cm. *Salmonella* serotype Branderup H9812 was utilized as DNA molecular mass size marker. PFGE pattern of isolates 2012 and 2013 were available from our previous researches [17, 18]. PFGE banding patterns were analyzed by BioNumerics software (version 7.6.3; Applied Maths, Belgium). The dendrogram was generated using Dice similarity index and UPGMA method. The strains with more than > 30% similarity fell into major clusters, then types and subtypes were assigned according to criteria presented by Tenover et al. [26].

## Multi-locus Sequence Typing (MLST)

Genomic DNA was extracted with AccuPrep genomic extraction kit (Bioneer, Korea) for all amplification tests. Seven housekeeping genes (*adk*, *gyrB*, *metE*, *mdh*, *pntA*, *purM*, *pyrC*) were used for MLST analysis [8]. Amplification of genes was performed by primers that are shown in Table 1. Each PCR was performed with an initial denaturation at 95°C for 3 min, denaturing at 95°C for 30 s, extension at 72°C for 1 min and a final extension at 72°C for 2 min repeated for 34 cycles. Annealing temperatures were adjusted according to melting temperature of each primer sets. PCR reactions were carried out in 25 µl total reaction volume containing 8 µl PCR-Master mix (Taq DNA Polymerase 2×Master Mix RED, 1.5 mM MgCl2) (Ampliqon, Denmark), 2 µl template DNA, 1 µl of each of the primers (10 pmol), and 13 µl sterilized distilled water. The DNA purification and sanger sequencing were performed by Microsynth (Switzerland). Sequences were trimmed, placed in the right direction and aligned at the same start and stop positions by CLC Genomic Workbench (version 12.0; Qiagen) and BioNumerics softwares (version 7.6.3; Applied Maths, Belgium). The alleles were uploaded to PubMLST and assigning the alleles and sequence type of strains were defined according to PubMLST website (https://pubmlst.org/organisms/vibrio-cholerae).

## Real-time PCR for toxin gene detection

Real-time PCR was used for detection of cholera toxin (*ctxA/ctxB* genes) and toxin co-regulated pilus (*tcpA* gene). Primer sequences are described in Table 1. The SYBR green Premix EX Taq mixture (Takara, Japan) was in a total volume of 20µL containing 0.8µL of each primer pair, 2µL of DNA sample, and 6.4µL of distilled water. RT- PCR cycling conditions were as follows: 95°C for 5 min and 40 cycles of 95°C for 15s, 60°C for 43 s.

## Whole-genome sequencing (WGS) and analysis of selected strains

Genomic DNA was extracted using AccuPrep genomic extraction kit (Bioneer, Korea). Whole genome sequencing was performed using the NovaSeq 6000 and HiSeq 2500 instruments (Illumina, San Diego, California, USA). DNA libraries were prepared using a Nextera DNA Flex Library Preparation Kit (Illumina), followed by 2×150 bp paired-end sequencing runs with a median coverage of 171 X (95–628). Short-read sequence data were submitted in the European Nucleotide Archive (ENA) (http://www.ebi.ac.uk/ena) and the accession numbers are listed in Table 2.

Short reads were first cleaned using FqCleanER version 3.0 (https://gitlab.pasteur.fr/GIPhy/fqCleanER)
<table><caption>**Table 1** Primers used in this study</caption><thead><tr><th>Gene</th><th>Gene product</th><th>Oligonucleotide sequence</th><th>Product size (bp)</th><th>Annealing</th><th>Reference</th></tr></thead><tbody><tr><td>adk</td><td>adenylate kinase</td><td>CATCATTCTTCTCGGTGCTC<br>AGTGCCGTCAAACTTCAGGTA</td><td>592</td><td>59°C</td><td>[8]</td></tr><tr><td>gyrB</td><td>DNA gyrase subunit B</td><td>GTACGTTTCTGGCCTAGTGC<br>GGGTCTTTTTCCTGACAATC</td><td>749</td><td>60°C</td><td>[8]</td></tr><tr><td>metE</td><td>methionine synthase</td><td>CGGGTGACTTTGCTTGGT<br>CAGATCGACTGGGCTGTG</td><td>827</td><td>58°C</td><td>[8]</td></tr><tr><td>mdh</td><td>malate dehydrogenase</td><td>ATGAAAGTCGCTGTTATTGG<br>GCCGCTTGGCCCATAGAAAG</td><td>692</td><td>60°C</td><td>[8]</td></tr><tr><td>pntA</td><td>pyridine nucleotide transhydrogenase</td><td>GGCCAGCCCAAAATCCT<br>GATATTGCCGTCTTTTTCTT</td><td>758</td><td>59°C</td><td>[8]</td></tr><tr><td>purM</td><td>phosphoribosylformylglycinamide cyclo-ligase</td><td>GGTGTCGATATTGATGCAGG<br>GGAATGTTTTCCCAGAAGCC</td><td>734</td><td>58.4°C</td><td>[8]</td></tr><tr><td>pyrC</td><td>Dihydroorotase</td><td>ATCATGCCTAACACGGTTCC<br>TTCAAACTTCGGCATA</td><td>726</td><td>56°C</td><td>[8]</td></tr><tr><td>ctxA/B</td><td>Cholera toxin</td><td>TATGCCAAGAGGACAGAGTGAG<br>AACATATCCATCATCGTGCCTAAC</td><td>115</td><td>55°C</td><td>[27]</td></tr><tr><td>tcpA</td><td>Toxin co-regulated pilus subunit A</td><td>GGTCAGCCTTGGTAAGGTT<br>CAAGATCAGCGACAGCAG</td><td>225</td><td>50°C</td><td>This study</td></tr></tbody></table>

to eliminate adaptor sequences [28], correct sequencing errors [29], and discard low-quality short-reads. Assemblies were generated de novo by SPAdes version 3.15.0 [30] with default settings.

Genome assemblies were aligned with the reference genome of Vibrio cholerae O1 El Tor N16961 using Blast software (https://blast.ncbi.nlm.nih.gov). Different genetic markers were analyzed against reference sequences of the O-antigen biosynthetic gene clusters (O-AGCs), CTX prophage, the ctxB gene, the toxin co-regulated pilus (TCP) genes, the Vibrio pathogenicity islands 1/2, and the Vibrio seventh pandemic islands I/II (Table 3).

The presence and type of genomic islands (GIs), acquired antibiotic resistance genes (ARGs) and plasmids were determined with VCGIDB (http://leb.snu.ac.kr/vcgidb/index) [31], ResFinder v4.0.1 (https://cge.cbs.dtu.dk/services/ResFinder/) [32] and PlasmidFinder v1.3 (https://cge.cbs.dtu.dk/services/PlasmidFinder/) [33].

## Additional genomic data

In July of 2021, only five genomic sequences from ST579 (n=2) and ST178 (n=3) strains were available via PubMLST (https://pubmlst.org/organisms/vibrio-cholerae). Fasta sequence files from available sequence data were downloaded and included in the analyses (Table 2).

## Results

### Identification and confirmation of V. cholerae isolates

Eighteen V. cholerae isolates were confirmed as O1 serogroup. One was O2 and the remaining isolate was O7. According to VP and the sheep red blood cells hemagglutination tests, 19 out of twenty isolates were identified as El Tor biotype and only one isolate (V. cholerae O2) showed phenotypic characteristics of classical biotype.

### Antimicrobial susceptibility testing

The antimicrobial susceptibility showed that all the isolates were susceptible to cefotaxime, ceftazidime, gentamicin, chloramphenicol and azithromycin. Resistance to trimethoprim-sulfamethoxazole, tetracycline and nalidixic acid was seen in 40% (8/20), 40% (8/20) and 65% (13/20) of isolates, respectively (Table 4). All strains except one (ID:391) were susceptible to ampicillin as determined by MIC method (Table 4).

All ST69 strains were resistant to nalidixic acid. Moreover, resistance to nalidixic acid, trimethoprim-sulfamethoxazole and tetracycline was only observed in ST69 strains.

### MLST analysis

All amplicon sequences were submitted and assigned by PubMLST. In total, seven sequence types (STs) were identified among 20 V. cholerae isolates in this study. The ST69 was the most abundant ST and was identified in 13 (65%) isolates (from 2012, 2013 and 2015). This
<table><caption>**Table 1** Primers used in this study</caption><thead><tr><th>Gene</th><th>Gene product</th><th>Oligonucleotide sequence</th><th>Product size (bp)</th><th>Annealing</th><th>Reference</th></tr></thead><tbody><tr><td>adk</td><td>adenylate kinase</td><td>CATCATTCTTCTCGGTGCTC<br>AGTGCCGTCAAACTTCAGGTA</td><td>592</td><td>59°C</td><td>[8]</td></tr><tr><td>gyrB</td><td>DNA gyrase subunit B</td><td>GTACGTTTCTGGCCTAGTGC<br>GGGTCTTTTTCCTGACAATC</td><td>749</td><td>60°C</td><td>[8]</td></tr><tr><td>metE</td><td>methionine synthase</td><td>CGGGTGACTTTGCTTGGT<br>CAGATCGACTGGGCTGTG</td><td>827</td><td>58°C</td><td>[8]</td></tr><tr><td>mdh</td><td>malate dehydrogenase</td><td>ATGAAAGTCGCTGTTATTGG<br>GCCGCTTGGCCCATAGAAAG</td><td>692</td><td>60°C</td><td>[8]</td></tr><tr><td>pntA</td><td>pyridine nucleotide transhydrogenase</td><td>GGCCAGCCCAAAATCCT<br>GATATTGCCGTCTTTTTCTT</td><td>758</td><td>59°C</td><td>[8]</td></tr><tr><td>purM</td><td>phosphoribosylformylglycinamide cyclo-ligase</td><td>GGTGTCGATATTGATGCAGG<br>GGAATGTTTTCCCAGAAGCC</td><td>734</td><td>58.4°C</td><td>[8]</td></tr><tr><td>pyrC</td><td>Dihydroorotase</td><td>ATCATGCCTAACACGGTTCC<br>TTCAAACTTCGGCATA</td><td>726</td><td>56°C</td><td>[8]</td></tr><tr><td>ctxA/B</td><td>Cholera toxin</td><td>TATGCCAAGAGGACAGAGTGAG<br>AACATATCCATCATCGTGCCTAAC</td><td>115</td><td>55°C</td><td>[27]</td></tr><tr><td>tcpA</td><td>Toxin co-regulated pilus subunit A</td><td>GGTCAGCCTTGGTAAGGTT<br>CAAGATCAGCGACAGCAG</td><td>225</td><td>50°C</td><td>This study</td></tr></tbody></table>

to eliminate adaptor sequences [28], correct sequencing errors [29], and discard low-quality short-reads. Assemblies were generated de novo by SPAdes version 3.15.0 [30] with default settings.

Genome assemblies were aligned with the reference genome of Vibrio cholerae O1 El Tor N16961 using Blast software (https://blast.ncbi.nlm.nih.gov). Different genetic markers were analyzed against reference sequences of the O-antigen biosynthetic gene clusters (O-AGCs), CTX prophage, the *ctxB* gene, the toxin co-regulated pilus (TCP) genes, the Vibrio pathogenicity islands 1/2, and the Vibrio seventh pandemic islands I/II (Table 3).

The presence and type of genomic islands (GIs), acquired antibiotic resistance genes (ARGs) and plasmids were determined with VCGIDB (http://leb.snu.ac.kr/vcgidb/index) [31], ResFinder v4.0.1 (https://cge.cbs.dtu.dk/services/ResFinder/) [32] and PlasmidFinder v1.3 (https://cge.cbs.dtu.dk/services/PlasmidFinder/) [33].

## Additional genomic data

In July of 2021, only five genomic sequences from ST579 (n=2) and ST178 (n=3) strains were available via PubMLST (https://pubmlst.org/organisms/vibrio-cholerae). Fasta sequence files from available sequence data were downloaded and included in the analyses (Table 2).

## Results

### Identification and confirmation of V. cholerae isolates

Eighteen V. cholerae isolates were confirmed as O1 serogroup. One was O2 and the remaining isolate was O7. According to VP and the sheep red blood cells hemagglutination tests, 19 out of twenty isolates were identified as El Tor biotype and only one isolate (V. cholerae O2) showed phenotypic characteristics of classical biotype.

### Antimicrobial susceptibility testing

The antimicrobial susceptibility showed that all the isolates were susceptible to cefotaxime, ceftazidime, gentamicin, chloramphenicol and azithromycin. Resistance to trimethoprim-sulfamethoxazole, tetracycline and nalidixic acid was seen in 40% (8/20), 40% (8/20) and 65% (13/20) of isolates, respectively (Table 4). All strains except one (ID:391) were susceptible to ampicillin as determined by MIC method (Table 4).

All ST69 strains were resistant to nalidixic acid. Moreover, resistance to nalidixic acid, trimethoprim-sulfamethoxazole and tetracycline was only observed in ST69 strains.

### MLST analysis

All amplicon sequences were submitted and assigned by PubMLST. In total, seven sequence types (STs) were identified among 20 V. cholerae isolates in this study. The ST69 was the most abundant ST and was identified in 13 (65%) isolates (from 2012, 2013 and 2015). This
<table><caption>**Table 2** The non-ST69/non-ST75 strains characterized via genomic analysis</caption><thead><tr><th>Isolate ID</th><th>Genome ID<br>(accession no)/<br>Database</th><th>Genomic<br>coverage<br>(X)</th><th>Country of<br>isolation</th><th>Year of<br>isolation</th><th>ST</th><th>Source of<br>isolation</th><th>O-serogroup<br>(in silico)</th><th>AMR<br>phenotype</th><th>AMR gene</th><th>Plasmids</th><th>Genomic<br>Islands</th><th>ctxA/B</th><th>tcpA</th><th>Reference</th></tr></thead><tbody><tr><td>195</td><td>2016_1<br>s(ERR6323226)/<br>(EBI-ENA)</td><td>195</td><td>Iran</td><td>2016</td><td>579</td><td>Human</td><td>O1</td><td>Pan suscep-<br>tible</td><td>qnrC4</td><td>/</td><td>VPI-1</td><td>/</td><td>wt</td><td>This study</td></tr><tr><td>2134</td><td>2134/PubMLST /</td><td>India</td><td>2017</td><td>579</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>qnrC4</td><td>/</td><td>VPI-1</td><td>/</td><td>wt</td><td>PubMLST</td></tr><tr><td>2106</td><td>2106/PubMLST /</td><td>Morocco</td><td>2015</td><td>579</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>qnrC4</td><td>/</td><td>VPI-1</td><td>/</td><td>wt</td><td>PubMLST</td></tr><tr><td>295</td><td>CNRVC170253<br>(ERR6323228)/<br>(EBI-ENA)</td><td>171</td><td>Iran</td><td>2016</td><td>178</td><td>Human</td><td>O1</td><td>Pan suscep-<br>tible</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>This study</td></tr><tr><td>1884</td><td>1884/PubMLST /</td><td>Russia</td><td>2004</td><td>178</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>PubMLST</td></tr><tr><td>1402</td><td>1402/PubMLST /</td><td>Ukraine</td><td>2016</td><td>178</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>PubMLST</td></tr><tr><td>1403</td><td>1403/PubMLST /</td><td>Ukraine</td><td>2017</td><td>178</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>PubMLST</td></tr><tr><td>491</td><td>2012_4<br>(ERR6323225)/<br>(EBI-ENA)</td><td>628</td><td>Iran</td><td>2012</td><td>438</td><td>Human</td><td>O2</td><td>Pan suscep-<br>tible</td><td>/</td><td>/</td><td>/</td><td>/</td><td>/</td><td>This study</td></tr><tr><td>391</td><td>CNRVC170247<br>(ERR6323227)/<br>(EBI-ENA)</td><td>126</td><td>Iran</td><td>2012</td><td>983</td><td>Human</td><td>O7</td><td>Ampicillin</td><td>bla<sub>CARB-<br>/catB9</sub></td><td>/</td><td>/</td><td>/</td><td>/</td><td>This study</td></tr><tr><td>495</td><td>CNRVC170254<br>(ERR6323229)/<br>(EBI-ENA)</td><td>95</td><td>Iran</td><td>2016</td><td>984</td><td>Human</td><td>O1</td><td>Pan suscep-<br>tible</td><td>/</td><td>/</td><td>/</td><td>/</td><td>/</td><td>This study</td></tr></tbody></table>

Pan susceptible: Susceptible to all antibiotics

wt: Wild-type V. cholerae O1 El Tor Inaba strain N16961

(/): Not detected
<table><caption>**Table 2** The non-ST69/non-ST75 strains characterized via genomic analysis</caption><thead><tr><th>Isolate ID</th><th>Genome ID<br>(accession no)/<br>Database</th><th>Genomic<br>coverage<br>(X)</th><th>Country of<br>isolation</th><th>Year of<br>isolation</th><th>ST</th><th>Source of<br>isolation</th><th>O-serogroup<br>(in silico)</th><th>AMR<br>phenotype</th><th>AMR gene</th><th>Plasmids</th><th>Genomic<br>Islands</th><th>ctxA/B</th><th>tcpA</th><th>Reference</th></tr></thead><tbody><tr><td>195</td><td>2016_1<br>s(ERR6323226)/<br>(EBI-ENA)</td><td>195</td><td>Iran</td><td>2016</td><td>579</td><td>Human</td><td>O1</td><td>Pan suscep-<br>tible</td><td>qnrC4</td><td>/</td><td>VPI-1</td><td>/</td><td>wt</td><td>This study</td></tr><tr><td>2134</td><td>2134/PubMLST /</td><td>India</td><td>2017</td><td>579</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>qnrC4</td><td>/</td><td>VPI-1</td><td>/</td><td>wt</td><td>PubMLST</td></tr><tr><td>2106</td><td>2106/PubMLST /</td><td>Morocco</td><td>2015</td><td>579</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>qnrC4</td><td>/</td><td>VPI-1</td><td>/</td><td>wt</td><td>PubMLST</td></tr><tr><td>295</td><td>CNRVC170253<br>(ERR6323228)/<br>(EBI-ENA)</td><td>171</td><td>Iran</td><td>2016</td><td>178</td><td>Human</td><td>O1</td><td>Pan suscep-<br>tible</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>This study</td></tr><tr><td>1884</td><td>1884/PubMLST /</td><td>Russia</td><td>2004</td><td>178</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>PubMLST</td></tr><tr><td>1402</td><td>1402/PubMLST /</td><td>Ukraine</td><td>2016</td><td>178</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>PubMLST</td></tr><tr><td>1403</td><td>1403/PubMLST /</td><td>Ukraine</td><td>2017</td><td>178</td><td>Human</td><td>O1</td><td>No data</td><td>No data</td><td>/</td><td>/</td><td>VSP-1</td><td>/</td><td>/</td><td>PubMLST</td></tr><tr><td>491</td><td>2012_4<br>(ERR6323225)/<br>(EBI-ENA)</td><td>628</td><td>Iran</td><td>2012</td><td>438</td><td>Human</td><td>O2</td><td>Pan suscep-<br>tible</td><td>/</td><td>/</td><td>/</td><td>/</td><td>/</td><td>This study</td></tr><tr><td>391</td><td>CNRVC170247<br>(ERR6323227)/<br>(EBI-ENA)</td><td>126</td><td>Iran</td><td>2012</td><td>983</td><td>Human</td><td>O7</td><td>Ampicillin</td><td>bla<sub>CARB-<br>/catB9</sub></td><td>/</td><td>/</td><td>/</td><td>/</td><td>This study</td></tr><tr><td>495</td><td>CNRVC170254<br>(ERR6323229)/<br>(EBI-ENA)</td><td>95</td><td>Iran</td><td>2016</td><td>984</td><td>Human</td><td>O1</td><td>Pan suscep-<br>tible</td><td>/</td><td>/</td><td>/</td><td>/</td><td>/</td><td>This study</td></tr></tbody></table>

Pan susceptible: Susceptible to all antibiotics

wt: Wild-type V. cholerae O1 El Tor Inaba strain N16961

(/): Not detected
<table><caption>**Table 3** The reference DNA sequences used in this study</caption><thead><tr><th>Genomic sequence analyzed</th><th>Accession numbers (GenBank/DDBJ)</th></tr></thead><tbody><tr><td>Vibrio cholerae O1 El Tor N16961</td><td>AE003852.1 and AE003853.1</td></tr><tr><td>CTX prophage</td><td>CTX-1, AE003852, coordinates 15666967–1573281<br>CTX-2, CP001486, coordinates 852233–858550</td></tr><tr><td>ctxB gene</td><td>(AE003852, coordinates 1566967–1567341)</td></tr><tr><td>toxin co-regulated pilus (TCP) genes</td><td>(AE003852, coordinates 890449–891123)</td></tr><tr><td>Vibrio pathogenicity island 1</td><td>(VPI-1, AE003852, coordinates 873242–914124)</td></tr><tr><td>Vibrio pathogenicity island 2</td><td>(VPI-2, AE003852, coordinates 1896092–1952861)</td></tr><tr><td>Vibrio seventh pandemic island I</td><td>(VSP-I, AE003852, coordinates 175343–189380)</td></tr><tr><td>Vibrio seventh pandemic island II</td><td>(VSP-II, AE003852, coordinates 523156–550021)</td></tr><tr><td>O-AGC</td><td>rfb from O1 (LC594800), O2 (LC594901) and O7(LC594961)</td></tr></tbody></table>

<table><caption>**Table 4** Antimicrobial resistance phenotypes of studied *V. cholerae* isolates</caption><thead><tr><th>Isolate ID</th><th>Isolation Province</th><th>ST</th><th>Year of Isolation</th><th>GM</th><th>CHL</th><th>AZM</th><th>ERY</th><th>AMP</th><th>TET</th><th>NAL</th><th>SXT</th><th>CIP</th><th>CTX</th><th>CZA</th></tr></thead><tbody><tr><td>191</td><td>Baluchistan</td><td>69</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>391</td><td>Baluchistan</td><td>983</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>I</td><td>R</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>491</td><td>Baluchistan</td><td>438</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>891</td><td>Baluchistan</td><td>69</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>991</td><td>Baluchistan</td><td>69</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>192</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>292</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>392</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>492</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>592</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>194</td><td>Khuzestan</td><td>75</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>894</td><td>Khorasan</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>I</td><td>S</td><td>S</td></tr><tr><td>994</td><td>Khorasan</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>104</td><td>Tehran</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>I</td><td>S</td><td>S</td></tr><tr><td>114</td><td>Qazvin</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>124</td><td>Qom</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>I</td><td>S</td><td>S</td></tr><tr><td>195</td><td>Bushehr</td><td>579</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>295</td><td>Bushehr</td><td>178</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>I</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>395</td><td>Khuzestan</td><td>75</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>495</td><td>Khuzestan</td><td>984</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr></tbody></table>

Abbreviations: GM gentamicin, CHL chloramphenicol, AZM azithromycin, ERY erythromycin, AMP ampicillin, TET tetracycline, NAL nalidixic acid, SXT trimethoprim-sulfamethoxazole, CIP ciprofloxacin, CTX cefotaxime, CZA ceftazidime, S susceptible, R resistant

ST is the predominant ST among seventh pandemic *V. cholerae* El Tor isolates, which has also been reported from at least 37 countries (https://pubmlst.org/organisms/vibrio-cholerae). ST75, was the second most prevalent ST identified in 10% of our isolates (from 2015 and 2016). The remaining 5 STs, including 2 new-found STs (983 and 984), each were detected in a single isolate. Based on MLST data, the genetic relatedness of these seven STs was assessed by a minimum spanning tree (Fig. 1). Accordingly, ST75 is the most closely

related ST to ST69, and ST579 is a single-locus variant (SLV) of ST75.

Members of all STs were defined as O1/El Tor serogroup/biotype, except for ST983 (O7/El Tor) and ST438 (O2/Classical). The STs, their allele types and geographical location of *V. cholerae* isolates is summarized in Table 5 and Fig. 2. As depicted in Fig. 2, ST69 strains were more widely distributed around the country, while ST75 strains were restricted to South West Iran. Moreover, strains from 2012 and 2013 were isolated from South
<table><caption>**Table 3** The reference DNA sequences used in this study</caption><thead><tr><th>Genomic sequence analyzed</th><th>Accession numbers (GenBank/DDBJ)</th></tr></thead><tbody><tr><td>Vibrio cholerae O1 El Tor N16961</td><td>AE003852.1 and AE003853.1</td></tr><tr><td>CTX prophage</td><td>CTX-1, AE003852, coordinates 15666967–1573281<br>CTX-2, CP001486, coordinates 852233–858550</td></tr><tr><td>ctxB gene</td><td>(AE003852, coordinates 1566967–1567341)</td></tr><tr><td>toxin co-regulated pilus (TCP) genes</td><td>(AE003852, coordinates 890449–891123)</td></tr><tr><td>Vibrio pathogenicity island 1</td><td>(VPI-1, AE003852, coordinates 873242–914124)</td></tr><tr><td>Vibrio pathogenicity island 2</td><td>(VPI-2, AE003852, coordinates 1896092–1952861)</td></tr><tr><td>Vibrio seventh pandemic island I</td><td>(VSP-I, AE003852, coordinates 175343–189380)</td></tr><tr><td>Vibrio seventh pandemic island II</td><td>(VSP-II, AE003852, coordinates 523156–550021)</td></tr><tr><td>O-AGC</td><td>rfb from O1 (LC594800), O2 (LC594901) and O7(LC594961)</td></tr></tbody></table>

<table><caption>**Table 4** Antimicrobial resistance phenotypes of studied *V. cholerae* isolates</caption><thead><tr><th>Isolate ID</th><th>Isolation Province</th><th>ST</th><th>Year of Isolation</th><th>GM</th><th>CHL</th><th>AZM</th><th>ERY</th><th>AMP</th><th>TET</th><th>NAL</th><th>SXT</th><th>CIP</th><th>CTX</th><th>CZA</th></tr></thead><tbody><tr><td>191</td><td>Baluchistan</td><td>69</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>391</td><td>Baluchistan</td><td>983</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>I</td><td>R</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>491</td><td>Baluchistan</td><td>438</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>891</td><td>Baluchistan</td><td>69</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>991</td><td>Baluchistan</td><td>69</td><td>2012</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>192</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>292</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>392</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>492</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>592</td><td>Baluchistan</td><td>69</td><td>2013</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>R</td><td>R</td><td>S</td><td>S</td><td>S</td></tr><tr><td>194</td><td>Khuzestan</td><td>75</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>894</td><td>Khorasan</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>I</td><td>S</td><td>S</td></tr><tr><td>994</td><td>Khorasan</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>104</td><td>Tehran</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>I</td><td>S</td><td>S</td></tr><tr><td>114</td><td>Qazvin</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>124</td><td>Qom</td><td>69</td><td>2015</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>R</td><td>S</td><td>I</td><td>S</td><td>S</td></tr><tr><td>195</td><td>Bushehr</td><td>579</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>295</td><td>Bushehr</td><td>178</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>I</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>395</td><td>Khuzestan</td><td>75</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr><tr><td>495</td><td>Khuzestan</td><td>984</td><td>2016</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td><td>S</td></tr></tbody></table>

Abbreviations: GM gentamicin, CHL chloramphenicol, AZM azithromycin, ERY erythromycin, AMP ampicillin, TET tetracycline, NAL nalidixic acid, SXT trimethoprim-sulfamethoxazole, CIP ciprofloxacin, CTX cefotaxime, CZA ceftazidime, S susceptible, R resistant

ST is the predominant ST among seventh pandemic *V. cholerae* El Tor isolates, which has also been reported from at least 37 countries (https://pubmlst.org/organisms/vibrio-cholerae). ST75, was the second most prevalent ST identified in 10% of our isolates (from 2015 and 2016). The remaining 5 STs, including 2 new-found STs (983 and 984), each were detected in a single isolate. Based on MLST data, the genetic relatedness of these seven STs was assessed by a minimum spanning tree (Fig. 1). Accordingly, ST75 is the most closely

related ST to ST69, and ST579 is a single-locus variant (SLV) of ST75.

Members of all STs were defined as O1/El Tor serogroup/biotype, except for ST983 (O7/El Tor) and ST438 (O2/Classical). The STs, their allele types and geographical location of *V. cholerae* isolates is summarized in Table 5 and Fig. 2. As depicted in Fig. 2, ST69 strains were more widely distributed around the country, while ST75 strains were restricted to South West Iran. Moreover, strains from 2012 and 2013 were isolated from South
<figure><img src="image_3.png" /><figcaption>**Fig. 1** The minimum spanning tree based on STs of V. cholerae strains. The size of nodes represents STs numbers and distance between nodes represents the locus differences. The digits next to the lines between the two nodes represent the number of different loci</figcaption></figure>

<table><caption>**Table 5** Allele types and ST numbers of *V. cholerae* strains in relation to their virulence genes content</caption><thead><tr><th>Isolate ID</th><th>Year of Isolation</th><th>adk</th><th>gyrB</th><th>mdh</th><th>metE</th><th>pntA</th><th>purM</th><th>pyrC</th><th>ST</th><th>Serogroup</th><th>ctxA/B</th><th>tcpA</th></tr></thead><tbody><tr><td>191</td><td>2012</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>391</td><td>2012</td><td>1</td><td>153</td><td>68</td><td>243</td><td>2</td><td>43</td><td>1</td><td>983</td><td>O7</td><td>-</td><td>-</td></tr><tr><td>491</td><td>2012</td><td>2</td><td>44</td><td>11</td><td>64</td><td>3</td><td>8</td><td>43</td><td>438</td><td>O2</td><td>-</td><td>-</td></tr><tr><td>891</td><td>2012</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>991</td><td>2012</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>192</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>292</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>392</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>492</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>592</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>194</td><td>2015</td><td>7</td><td>2</td><td>4</td><td>37</td><td>12</td><td>1</td><td>38</td><td>75</td><td>O1</td><td>-</td><td>+</td></tr><tr><td>894</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>994</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>104</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>114</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>124</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>195</td><td>2016</td><td>7</td><td>2</td><td>4</td><td>13</td><td>12</td><td>1</td><td>38</td><td>579</td><td>O1</td><td>-</td><td>+</td></tr><tr><td>295</td><td>2016</td><td>1</td><td>31</td><td>56</td><td>99</td><td>70</td><td>48</td><td>13</td><td>178</td><td>O1</td><td>-</td><td>-</td></tr><tr><td>395</td><td>2016</td><td>7</td><td>2</td><td>4</td><td>37</td><td>12</td><td>1</td><td>38</td><td>75</td><td>O1</td><td>-</td><td>+</td></tr><tr><td>495</td><td>2016</td><td>2</td><td>38</td><td>161</td><td>230</td><td>69</td><td>1</td><td>259</td><td>984</td><td>O1</td><td>-</td><td>-</td></tr></tbody></table>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** The minimum spanning tree based on STs of V. cholerae strains. The size of nodes represents STs numbers and distance between nodes represents the locus differences. The digits next to the lines between the two nodes represent the number of different loci</figcaption></figure>

<table><caption>**Table 5** Allele types and ST numbers of V. cholerae strains in relation to their virulence genes content</caption><thead><tr><th>Isolate ID</th><th>Year of Isolation</th><th>adk</th><th>gyrB</th><th>mdh</th><th>metE</th><th>pntA</th><th>purM</th><th>pyrC</th><th>ST</th><th>Serogroup</th><th>ctxA/B</th><th>tcpA</th></tr></thead><tbody><tr><td>191</td><td>2012</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>391</td><td>2012</td><td>1</td><td>153</td><td>68</td><td>243</td><td>2</td><td>43</td><td>1</td><td>983</td><td>O7</td><td>-</td><td>-</td></tr><tr><td>491</td><td>2012</td><td>2</td><td>44</td><td>11</td><td>64</td><td>3</td><td>8</td><td>43</td><td>438</td><td>O2</td><td>-</td><td>-</td></tr><tr><td>891</td><td>2012</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>991</td><td>2012</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>192</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>292</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>392</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>492</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>592</td><td>2013</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>194</td><td>2015</td><td>7</td><td>2</td><td>4</td><td>37</td><td>12</td><td>1</td><td>38</td><td>75</td><td>O1</td><td>-</td><td>+</td></tr><tr><td>894</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>994</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>104</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>114</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>124</td><td>2015</td><td>7</td><td>11</td><td>4</td><td>37</td><td>12</td><td>1</td><td>20</td><td>69</td><td>O1</td><td>+</td><td>+</td></tr><tr><td>195</td><td>2016</td><td>7</td><td>2</td><td>4</td><td>13</td><td>12</td><td>1</td><td>38</td><td>579</td><td>O1</td><td>-</td><td>+</td></tr><tr><td>295</td><td>2016</td><td>1</td><td>31</td><td>56</td><td>99</td><td>70</td><td>48</td><td>13</td><td>178</td><td>O1</td><td>-</td><td>-</td></tr><tr><td>395</td><td>2016</td><td>7</td><td>2</td><td>4</td><td>37</td><td>12</td><td>1</td><td>38</td><td>75</td><td>O1</td><td>-</td><td>+</td></tr><tr><td>495</td><td>2016</td><td>2</td><td>38</td><td>161</td><td>230</td><td>69</td><td>1</td><td>259</td><td>984</td><td>O1</td><td>-</td><td>-</td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>**Fig. 2** The distribution of V. cholerae strains in this study according to the sequence types and year of isolation. Sequence types detected in each province are shown by ST number. Members of each ST are defined within parenthesis. Isolates belonging to each year are depicted by a distinct color. As depicted in picture, strains of non-ST69/non-ST75 were mainly isolated from border provinces in south of Iran. Free map was obtained from Wikimedia commons (https://commons.wikimedia.org/wiki/File:Iran_location_map.svg)</figcaption></figure>

East Iran, the Baluchistan province which shares a common border with Afghanistan and proposes their probable import from this neighbor country. It should be noted that cholera cases in border provinces in south of Iran were solely of non-ST69/non-ST75 (Fig. 2), raising the hypothesis of their importing from our neighboring countries, all of those countries encounter cholera epidemics annually. A phylogenetic tree of V. cholerae isolates with ST 178, 438, 983, 984 and 579 from PubMLST is shown in Fig. 3. The tree is reflecting the relatedness of non-ST69/non-ST75 according to MLST data. Accordingly, the ST983 and ST984 are the two most closely related STs, which together with ST438 comprising a distinct cluster from ST178 and ST579.

## Detection of toxin genes

Real-time PCR showed the presence of *ctxA/B* gene in 13/20 (65%) strains (Table 5). In fact, all ST69 strains harbored the *ctxA/B* genes while, all non-ST69 strains lacked the toxin genes. The *tcpA* gene was detected in 16/20 (75%) strains (strains of ST69, ST75 and ST579 genotypes).

## PFGE analysis

Seven different pulsotype letters (A-G) were assigned to 20 V. cholerae strains which were analyzed based on PFGE method. Eleven isolates (from 2012, 2013 and 2015) showed identical or very similar pulsotypes, which was considered as the dominant pulsotype (A) of this study. All the isolates with pulsotype A uniformly belonged to ST69 (Fig. 4).

## Whole-genome sequence analyses of ST178, ST983, ST984, ST438 and ST579 strains

Except for well-known ST69 and ST75 types, little information is available on genomic data and gene content of strains belonging to other sequence types. Therefore, non-ST69/non-ST75 strains in this study were subjected to WGS to make a conceptual analysis and comparison (Table 2).

WGS analysis of non-ST69/non-ST75 strains including 5 from Iran (this study) and 5 from PubMLST showed that none of them harbored *ctxA* or *ctxB* genes, while a wild type-El Tor *tcpA* was present only in strains of ST579 (Table 2). The O1-antigen gene (*rfbO1*) was
<figure><img src="image_4.png" /><figcaption>**Fig. 2** The distribution of V. cholerae strains in this study according to the sequence types and year of isolation. Sequence types detected in each province are shown by ST number. Members of each ST are defined within parenthesis. Isolates belonging to each year are depicted by a distinct color. As depicted in picture, strains of non-ST69/non-ST75 were mainly isolated from border provinces in south of Iran. Free map was obtained from Wikimedia commons (https://commons.wikimedia.org/wiki/File:Iran_location_map.svg)</figcaption></figure>

East Iran, the Baluchistan province which shares a common border with Afghanistan and proposes their probable import from this neighbor country. It should be noted that cholera cases in border provinces in south of Iran were solely of non-ST69/non-ST75 (Fig. 2), raising the hypothesis of their importing from our neighboring countries, all of those countries encounter cholera epidemics annually. A phylogenetic tree of V. cholerae isolates with ST 178, 438, 983, 984 and 579 from PubMLST is shown in Fig. 3. The tree is reflecting the relatedness of non-ST69/non-ST75 according to MLST data. Accordingly, the ST983 and ST984 are the two most closely related STs, which together with ST438 comprising a distinct cluster from ST178 and ST579.

## Detection of toxin genes

Real-time PCR showed the presence of *ctxA/B* gene in 13/20 (65%) strains (Table 5). In fact, all ST69 strains harbored the *ctxA/B* genes while, all non-ST69 strains lacked the toxin genes. The *tcpA* gene was detected in 16/20 (75%) strains (strains of ST69, ST75 and ST579 genotypes).

## PFGE analysis

Seven different pulsotype letters (A-G) were assigned to 20 V. cholerae strains which were analyzed based on PFGE method. Eleven isolates (from 2012, 2013 and 2015) showed identical or very similar pulsotypes, which was considered as the dominant pulsotype (A) of this study. All the isolates with pulsotype A uniformly belonged to ST69 (Fig. 4).

## Whole-genome sequence analyses of ST178, ST983, ST984, ST438 and ST579 strains

Except for well-known ST69 and ST75 types, little information is available on genomic data and gene content of strains belonging to other sequence types. Therefore, non-ST69/non-ST75 strains in this study were subjected to WGS to make a conceptual analysis and comparison (Table 2).

WGS analysis of non-ST69/non-ST75 strains including 5 from Iran (this study) and 5 from PubMLST showed that none of them harbored *ctxA* or *ctxB* genes, while a wild type-El Tor *tcpA* was present only in strains of ST579 (Table 2). The O1-antigen gene (*rfbO1*) was
<figure><img src="image_5.png" /><figcaption>**Fig. 3** A phylogenetic tree of all V. cholerae isolates with ST178, ST438, ST983, ST984 and ST579 sequence types which have been deposited in PubMLST (https://pubmlst.org/organisms/vibrio-cholerae). The tree was constructed based on MLST profiles using MEGA software (version 11) based on the maximum composite likelihood and UPGMA. * V. cholerae isolates with WGS data used in this study</figcaption></figure>

present in all strains of ST579, ST178 and ST984 indica-
tive of O1 serogroup specificity.

Wild type-El Tor sequences of VC-1319, VC-1320,
VC-1577, VC-1578 genes (responsible for polymyxin
resistance in El Tor biotype) was detected in all strains of
ST579 (from Iran, India and Morocco). These genes are
supposed to be altered in Classical strains which restore
their susceptibility to polymyxin B [34, 35].

Moreover, the *rstC* gene, the signature of RS1 phage,
was also present in one genome of ST579 (ID:2106
from pubMLST). This means that ST579 genomes are
prone for acquisition of CTX page and its satellite RS1,
although its stability is not well guaranteed.

The beta-lactamase *bla*~CARB-7~ gene was detected in an
ampicillin-resistant ST983 strain from Iran (ID:391).

The *catB9* gene was also found in this ST983 strain,
which was phenotypically susceptible to chlorampheni-
col. As indicated in previous studies, this gene is a silent
chloramphenicol acetyl transferase gene within a super
integron of *V. cholerae* which is not considered as indic-
ative of resistance when identified in whole-genome
sequences of *V. cholerae* El Tor [36, 37]. The *qnrVC4*
gene (Quinolone resistance pentapeptide repeat protein
QnrVC4) was detected in all three ST579 strains. This

gene has also been found in quinolone-susceptible strains
of ST75 (SLV of ST579) in South Africa [38]. This gene
different alleles of which are located in a chromosomal
super-integron, is not necessarily related to resistance as
indicated here and in several previous studies [39, 40].

Discussion

As an endemic disease, cholera outbreaks are annually
reported from Iran; However, cross-border cases from
neighbor countries including Pakistan, Afghanistan
and Iraq, are also observed in some provinces. To our
knowledge, approximately 400 cholera cases have been
reported from Iran during 2012–2016 (mainly May–
November) [41, 42]. This prompted us to undertake the
genetic characterization of recent cholera outbreaks and
better understanding of endemic cholera.

The significant contribution of current study was the
investigation of prevalent STs in accordance with related
pulsotypes in Iran during 2012–2016. The 7th pandemic
clone, ST69, was revealed as the predominant circulat-
ing clone in Iran during 2012–2016 and all strains in
this ST appeared as pulsotype A or F and their subtypes
(A1 to A6) (F1 and F2). Subtypes of each pulsotype are
indicative of 1–2 genetic event, i.e., a point mutation or
<figure><img src="image_5.png" /><figcaption>**Fig. 3** A phylogenetic tree of all V. cholerae isolates with ST178, ST438, ST983, ST984 and ST579 sequence types which have been deposited in PubMLST (https://pubmlst.org/organisms/vibrio-cholerae). The tree was constructed based on MLST profiles using MEGA software (version 11) based on the maximum composite likelihood and UPGMA. * V. cholerae isolates with WGS data used in this study</figcaption></figure>

present in all strains of ST579, ST178 and ST984 indica-
tive of O1 serogroup specificity.

Wild type-El Tor sequences of VC-1319, VC-1320,
VC-1577, VC-1578 genes (responsible for polymyxin
resistance in El Tor biotype) was detected in all strains of
ST579 (from Iran, India and Morocco). These genes are
supposed to be altered in Classical strains which restore
their susceptibility to polymyxin B [34, 35].

Moreover, the *rstC* gene, the signature of RS1 phage,
was also present in one genome of ST579 (ID:2106
from pubMLST). This means that ST579 genomes are
prone for acquisition of CTX page and its satellite RS1,
although its stability is not well guaranteed.

The beta-lactamase *bla*~CARB-7~ gene was detected in an
ampicillin-resistant ST983 strain from Iran (ID:391).

The *catB9* gene was also found in this ST983 strain,
which was phenotypically susceptible to chlorampheni-
col. As indicated in previous studies, this gene is a silent
chloramphenicol acetyl transferase gene within a super
integron of *V. cholerae* which is not considered as indic-
ative of resistance when identified in whole-genome
sequences of *V. cholerae* El Tor [36, 37]. The *qnrVC4*
gene (Quinolone resistance pentapeptide repeat protein
QnrVC4) was detected in all three ST579 strains. This

gene has also been found in quinolone-susceptible strains
of ST75 (SLV of ST579) in South Africa [38]. This gene
different alleles of which are located in a chromosomal
super-integron, is not necessarily related to resistance as
indicated here and in several previous studies [39, 40].

Discussion

As an endemic disease, cholera outbreaks are annually
reported from Iran; However, cross-border cases from
neighbor countries including Pakistan, Afghanistan
and Iraq, are also observed in some provinces. To our
knowledge, approximately 400 cholera cases have been
reported from Iran during 2012–2016 (mainly May–
November) [41, 42]. This prompted us to undertake the
genetic characterization of recent cholera outbreaks and
better understanding of endemic cholera.

The significant contribution of current study was the
investigation of prevalent STs in accordance with related
pulsotypes in Iran during 2012–2016. The 7th pandemic
clone, ST69, was revealed as the predominant circulat-
ing clone in Iran during 2012–2016 and all strains in
this ST appeared as pulsotype A or F and their subtypes
(A1 to A6) (F1 and F2). Subtypes of each pulsotype are
indicative of 1–2 genetic event, i.e., a point mutation or
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Pulsed-field Gel Electrophoresis (PFGE) patterns of V. cholerae strains in this study. Each strain is shown in relation to year of isolation, pulsotype name and related sequence type (ST). V. cholerae strains are classified to three clusters A, B and C based on a similarity cutoff of 30%. Moreover, each pulsotype is assigned a distinct letter (A to G) and subtypes of each pulsotype are defined with numerical subscripts. * Full-length blots/gels associated to PFGE patterns of V. cholerae strains in current study are presented in Supplementary Fig. 1</figcaption></figure>

an insertion or deletion of DNA and are considered as minor differences from predominant pulsotype [26], consistent with ST69 sub-lineages [10].

Interestingly, all isolates of the present study with pulsotype A, uniformly fell in ST69. Moreover, the two ST75 isolates together with a single locus variant of ST75 (ST579) were assigned to a separate pulsotype letter (pulsotype E). This may validate PFGE as a mirror of MLST typing and a powerful epidemiological tool in *V. cholerae* investigations. The global circulation of ST69 might lead to some genetic events (mutations, deletion, insertion) within the genome content of *V. cholerae* strains which may be reflected in minor differences in PFGE patterns. By careful consideration of studies reporting PFGE patterns of ST69 strains, these minor variations in PFGE patterns of ST69 strains could be clearly elicited [43]. This means that strains belonging to ST69 might show minor differences in pulsotype, despite of identical ST type. ST69 as cholera 7th pandemic strain, has been reported from Asia, Africa, America and Europe [43–46].

ST75, as the second most prevalent ST in our study, has also been identified in seven countries including Iran, China, Russia, USA, Ukraine and Thailand and South

Africa (according to pubMLST database), but its frequency did not display considerable changes among years except for 2009, which affected Thailand and became more prevalent than the ST69 clone afterward (Fig. 5). Subsequently, related ST75 strains have emerged in several countries and is now somehow widespread in Asia and America continents [8, 47, 48]. In China, ST75 was the most prevalent sequence type among non-7th pandemic clone strains in Zhejiang Province, its emergence and potential spread draw significant attention due its probable threat to public health [47]. Moreover, ST75 recently emerged and became more prevalent than the pandemic clone in South Africa during 2018–2020 [38]. The occurrence of two isolates of ST75 in 2015 and 2016 is noteworthy, indicating the probable beginning of sporadic cholera due to ST75 in Iran; however, it should be considered that they may be endemic to the region but not previously detected by the surveillance system. The advent of ST75 in provinces bordering the Persian Gulf, like some other countries which were affected by US indigenous Gulf Coast-like strains, can enhance the probability of ST75 overtaking in similar regions in global view.
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Pulsed-field Gel Electrophoresis (PFGE) patterns of V. cholerae strains in this study. Each strain is shown in relation to year of isolation, pulsotype name and related sequence type (ST). V. cholerae strains are classified to three clusters A, B and C based on a similarity cutoff of 30%. Moreover, each pulsotype is assigned a distinct letter (A to G) and subtypes of each pulsotype are defined with numerical subscripts. * Full-length blots/gels associated to PFGE patterns of V. cholerae strains in current study are presented in Supplementary Fig. 1</figcaption></figure>

an insertion or deletion of DNA and are considered as minor differences from predominant pulsotype [26], consistent with ST69 sub-lineages [10].

Interestingly, all isolates of the present study with pulsotype A, uniformly fell in ST69. Moreover, the two ST75 isolates together with a single locus variant of ST75 (ST579) were assigned to a separate pulsotype letter (pulsotype E). This may validate PFGE as a mirror of MLST typing and a powerful epidemiological tool in *V. cholerae* investigations. The global circulation of ST69 might lead to some genetic events (mutations, deletion, insertion) within the genome content of *V. cholerae* strains which may be reflected in minor differences in PFGE patterns. By careful consideration of studies reporting PFGE patterns of ST69 strains, these minor variations in PFGE patterns of ST69 strains could be clearly elicited [43]. This means that strains belonging to ST69 might show minor differences in pulsotype, despite of identical ST type. ST69 as cholera 7th pandemic strain, has been reported from Asia, Africa, America and Europe [43–46].

ST75, as the second most prevalent ST in our study, has also been identified in seven countries including Iran, China, Russia, USA, Ukraine and Thailand and South

Africa (according to pubMLST database), but its frequency did not display considerable changes among years except for 2009, which affected Thailand and became more prevalent than the ST69 clone afterward (Fig. 5). Subsequently, related ST75 strains have emerged in several countries and is now somehow widespread in Asia and America continents [8, 47, 48]. In China, ST75 was the most prevalent sequence type among non-7th pandemic clone strains in Zhejiang Province, its emergence and potential spread draw significant attention due its probable threat to public health [47]. Moreover, ST75 recently emerged and became more prevalent than the pandemic clone in South Africa during 2018–2020 [38]. The occurrence of two isolates of ST75 in 2015 and 2016 is noteworthy, indicating the probable beginning of sporadic cholera due to ST75 in Iran; however, it should be considered that they may be endemic to the region but not previously detected by the surveillance system. The advent of ST75 in provinces bordering the Persian Gulf, like some other countries which were affected by US indigenous Gulf Coast-like strains, can enhance the probability of ST75 overtaking in similar regions in global view.
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Geographic distribution of ST75, ST178, ST579 and ST438 around the world (PubMLST website, July 2021). A) According to country and number of isolates. B) According to identification year and number of isolates. * No value indicates that the year of isolation of V. cholerae strains is unknown</figcaption></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

ST438 which was detected among our isolates of 2012, was previously confirmed in three countries, Hungary, Tunisia and the US. This ST is, for the first time, reported from an Asian country (Fig. 5). The

single Iranian strain (ID: 491) belonging to this ST type, appeared as classical biotype according to phenotypic analyses. Due to little available data on sequence types of non-O1/non-O139 strains and considerable diversity
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Geographic distribution of ST75, ST178, ST579 and ST438 around the world (PubMLST website, July 2021). A) According to country and number of isolates. B) According to identification year and number of isolates. * No value indicates that the year of isolation of V. cholerae strains is unknown</figcaption></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

ST438 which was detected among our isolates of 2012, was previously confirmed in three countries, Hungary, Tunisia and the US. This ST is, for the first time, reported from an Asian country (Fig. 5). The

single Iranian strain (ID: 491) belonging to this ST type, appeared as classical biotype according to phenotypic analyses. Due to little available data on sequence types of non-O1/non-O139 strains and considerable diversity
among their population, careful consideration should
be undertaken in interpretations.

Large-scale cholera outbreaks have decreased in recent years in Iran due to improvements in hygienic level and lifestyle of people in this country, although small outbreaks and sporadic cholera cases mainly of non-ST69/non-ST75 strains is still an ongoing concern in bordering provinces. Genomic analysis of non-ST69/ non-ST75 strains in this study showed i) the presence of wild type sequence of *tcpA* in ST579 strains, ii) the wild type-El Tor sequence of VC-1319, VC-1320, VC-1577, VC-1578 genes (responsible for polymyxin resistance in El Tor biotype) in strains of ST579, iii) the *rstC* gene, indicative of RS1 phage in one of the ST579 strains, and iv) the presence of VPI-1 and VSP-I islands in ST579 and ST178 strains. These data, all together, indicate that ST579 has more virulence and pathogenic potential compared with other sequence types under study (ST178, ST983, ST984 and ST438). In silico analysis of resistance genes and cassettes showed that i) none of the known antibiotic resistance genes or cassettes are present among the non-ST69/non-ST75 strains under study, ii) SXT constin responsible for resistance to several antimicrobials was not detected, and iii) no plasmid was detected among the isolates. Overall, these data indicate the high susceptibility of *V. cholerae* non-ST69/non-ST75 strains in comparison with more ubiquitous and more globally ST69 and ST75 strains.

In conclusion, the occurrence of *V. cholerae* strains of non-ST69/non-ST75 sequence types with some traits of virulence factors in recent years in Iran is noteworthy. Extensive studies together with surveillance efforts are expected to determine their likely route of transport. The occurrence of these STs, although still susceptible to current antimicrobial agents, is important because they may gradually change in antibiotic resistance gene content. Moreover, the circulation of *V. cholerae* ST69 in recent years in Iran shows the 7th pandemic strains as the persistent causes of cholera outbreaks in this country, although the role of ST75 as the second most contributed ST should not be ignored.

**Supplementary Information**

The online version contains supplementary material available at https://doi.org/10.1186/s12866-024-03185-9.

**Additional file 1.**

**Acknowledgements**

Not Applicable.

**Authors' contributions**

FJ was collected the data and performed the analysis. EN and FXW were contributed the analysis. MRF and RS were collected bacterial isolates. FG was contributed editing the manuscript. BB was supervised the whole project.

**Funding**

This research was supported by Iran National Science Foundation (INSF) grant number 96000794. The funding body played no role in the design of the study and collection, analysis, interpretation of data, and in writing the manuscript.

**Availability of data and materials**

The confirmed DNA sequences of MLST analysis were deposited in GenBank under accession numbers MN026296-MN026329. Short-read sequence data from WGS were deposited in European Nucleotide Archive (https://www.ebi.ac.uk/) under accession numbers ERR6323225- ERR6323229.

All datasets of the current study are available within article or can be obtained from corresponding with no restriction.

**Declarations**

**Ethical approval and consent to participate**

This study protocol was approved by the Ethics Committee of Tarbiat Modares University (Code: IR.TMU.REC.1397.092). all methods were carried out in accordance with relevant guidelines and regulations. All patients completed a questionnaire and gave informed consent prior to the study.

**Consent for publication**

Not applicable.

**Competing interests**

The authors declare no competing interests.

**Author details**

¹Department of Bacteriology, Faculty of Medical Sciences, Tarbiat Modares University, Tehran, Iran. ²Institut Pasteur, Université Paris Cité, Paris, France.

³Research Center of Health Reference Laboratories, Tehran, Iran.

Received: 7 May 2023 Accepted: 8 January 2024

Published online: 22 January 2024

**References**

1. Pierce NF, Sack RB, Mitra RC, Banwell JG, Brigham KL, Fedson DS, et al. Replacement of water and electrolyte losses in cholera by an oral glucose—electrolyte solution. Ann Intern Med. 1969;70(6):1173–81.

2. Hirschhorn N, Kinzie JL, Sachar DB, Northrup RS, Taylor JO, Ahmad SZ, et al. Decrease in net stool output in cholera during intestinal perfusion with glucose-containing solutions. N Engl J Med. 1968;279(4):176–81.

3. Lacey SW. Cholera: calamitous past, ominous future. Clin Infect Dis. 1995;20(5):1409–19.

4. Ali M, Lopez AL, You YA, Kim YE, Sah B, Maskery B, et al. The global burden of cholera. Bull World Health Organ. 2012;90(3):209–18.

5. Foley SL, Lynne AM, Nayak R. Molecular typing methodologies for microbial source tracking and epidemiological investigations of gram-negative bacterial foodborne pathogens. Infect Genet Evol. 2009;9(4):430–40.

6. Olive DM, Bean P. Principles and applications of methods for DNA-based typing of microbial organisms. J Clin Microbiol. 1999;37(6):1661–9.

7. Rahaman M, Islam T, Colwell RR, Alam M. Molecular tools in understanding the evolution of Vibrio cholerae. Front Microbiol. 2015;6:1040.

8. Octavia S, Salim A, Kurniawan J, Lam C, Leung Q, Ahsan S, et al. Population structure and evolution of non-O1/non-O139 Vibrio cholerae by multilocus sequence typing. PLoS ONE. 2013;8(6):e65342.

9. Olsen JS, Aarskaug T, Skogan G, Fykse EM, Ellingsen AB, Blatny JM. Evaluation of a highly discriminating multiplex multi-locus variable-number of tandem-repeats (MLVA) analysis for Vibrio cholerae. J Microbiol Methods. 2009;78(3):271–85.

10. Liang KY, Orata FD, Islam MT, Nasreen T, Alam M, Tarr CL, et al. A Vibrio cholerae core genome multilocus sequence typing scheme to facilitate the epidemiological study of cholera. J Bacteriol. 2020;202(24):e00086-e120.

11. Goering RV. Pulsed field gel electrophoresis: a review of application and interpretation in the molecular epidemiology of infectious disease. Infect Genet Evol. 2010;10(7):866–75.
among their population, careful consideration should
be undertaken in interpretations.

Large-scale cholera outbreaks have decreased in recent years in Iran due to improvements in hygienic level and lifestyle of people in this country, although small outbreaks and sporadic cholera cases mainly of non-ST69/non-ST75 strains is still an ongoing concern in bordering provinces. Genomic analysis of non-ST69/ non-ST75 strains in this study showed i) the presence of wild type sequence of *tcpA* in ST579 strains, ii) the wild type-El Tor sequence of VC-1319, VC-1320, VC-1577, VC-1578 genes (responsible for polymyxin resistance in El Tor biotype) in strains of ST579, iii) the *rstC* gene, indicative of RS1 phage in one of the ST579 strains, and iv) the presence of VPI-1 and VSP-I islands in ST579 and ST178 strains. These data, all together, indicate that ST579 has more virulence and pathogenic potential compared with other sequence types under study (ST178, ST983, ST984 and ST438). In silico analysis of resistance genes and cassettes showed that i) none of the known antibiotic resistance genes or cassettes are present among the non-ST69/non-ST75 strains under study, ii) SXT constin responsible for resistance to several antimicrobials was not detected, and iii) no plasmid was detected among the isolates. Overall, these data indicate the high susceptibility of *V. cholerae* non-ST69/non-ST75 strains in comparison with more ubiquitous and more globally ST69 and ST75 strains.

In conclusion, the occurrence of *V. cholerae* strains of non-ST69/non-ST75 sequence types with some traits of virulence factors in recent years in Iran is noteworthy. Extensive studies together with surveillance efforts are expected to determine their likely route of transport. The occurrence of these STs, although still susceptible to current antimicrobial agents, is important because they may gradually change in antibiotic resistance gene content. Moreover, the circulation of *V. cholerae* ST69 in recent years in Iran shows the 7th pandemic strains as the persistent causes of cholera outbreaks in this country, although the role of ST75 as the second most contributed ST should not be ignored.

**Supplementary Information**

The online version contains supplementary material available at https://doi.org/10.1186/s12866-024-03185-9.

**Additional file 1.**

**Acknowledgements**

Not Applicable.

**Authors' contributions**

FJ was collected the data and performed the analysis. EN and FXW were contributed the analysis. MRF and RS were collected bacterial isolates. FG was contributed editing the manuscript. BB was supervised the whole project.

**Funding**

This research was supported by Iran National Science Foundation (INSF) grant number 96000794. The funding body played no role in the design of the study and collection, analysis, interpretation of data, and in writing the manuscript.

**Availability of data and materials**

The confirmed DNA sequences of MLST analysis were deposited in GenBank under accession numbers MN026296-MN026329. Short-read sequence data from WGS were deposited in European Nucleotide Archive (https://www.ebi.ac.uk/) under accession numbers ERR6323225- ERR6323229.

All datasets of the current study are available within article or can be obtained from corresponding with no restriction.

**Declarations**

**Ethical approval and consent to participate**

This study protocol was approved by the Ethics Committee of Tarbiat Modares University (Code: IR.TMU.REC.1397.092). all methods were carried out in accordance with relevant guidelines and regulations. All patients completed a questionnaire and gave informed consent prior to the study.

**Consent for publication**

Not applicable.

**Competing interests**

The authors declare no competing interests.

**Author details**

¹Department of Bacteriology, Faculty of Medical Sciences, Tarbiat Modares University, Tehran, Iran. ²Institut Pasteur, Université Paris Cité, Paris, France.

³Research Center of Health Reference Laboratories, Tehran, Iran.

Received: 7 May 2023 Accepted: 8 January 2024

Published online: 22 January 2024

**References**

1. Pierce NF, Sack RB, Mitra RC, Banwell JG, Brigham KL, Fedson DS, et al. Replacement of water and electrolyte losses in cholera by an oral glucose—electrolyte solution. Ann Intern Med. 1969;70(6):1173–81.

2. Hirschhorn N, Kinzie JL, Sachar DB, Northrup RS, Taylor JO, Ahmad SZ, et al. Decrease in net stool output in cholera during intestinal perfusion with glucose-containing solutions. N Engl J Med. 1968;279(4):176–81.

3. Lacey SW. Cholera: calamitous past, ominous future. Clin Infect Dis. 1995;20(5):1409–19.

4. Ali M, Lopez AL, You YA, Kim YE, Sah B, Maskery B, et al. The global burden of cholera. Bull World Health Organ. 2012;90(3):209–18.

5. Foley SL, Lynne AM, Nayak R. Molecular typing methodologies for microbial source tracking and epidemiological investigations of gram-negative bacterial foodborne pathogens. Infect Genet Evol. 2009;9(4):430–40.

6. Olive DM, Bean P. Principles and applications of methods for DNA-based typing of microbial organisms. J Clin Microbiol. 1999;37(6):1661–9.

7. Rahaman M, Islam T, Colwell RR, Alam M. Molecular tools in understanding the evolution of Vibrio cholerae. Front Microbiol. 2015;6:1040.

8. Octavia S, Salim A, Kurniawan J, Lam C, Leung Q, Ahsan S, et al. Population structure and evolution of non-O1/non-O139 Vibrio cholerae by multilocus sequence typing. PLoS ONE. 2013;8(6):e65342.

9. Olsen JS, Aarskaug T, Skogan G, Fykse EM, Ellingsen AB, Blatny JM. Evaluation of a highly discriminating multiplex multi-locus variable-number of tandem-repeats (MLVA) analysis for Vibrio cholerae. J Microbiol Methods. 2009;78(3):271–85.

10. Liang KY, Orata FD, Islam MT, Nasreen T, Alam M, Tarr CL, et al. A Vibrio cholerae core genome multilocus sequence typing scheme to facilitate the epidemiological study of cholera. J Bacteriol. 2020;202(24):e00086-e120.

11. Goering RV. Pulsed field gel electrophoresis: a review of application and interpretation in the molecular epidemiology of infectious disease. Infect Genet Evol. 2010;10(7):866–75.
12. Lukinmaa S, Nakari UM, Eklund M, Siitonen A. Application of molecular genetic methods in diagnostics and epidemiology of food-borne bacterial pathogens. *Apnis*. 2004;112(11-12):908-29.

13. Teh CSJ, Chua KH, Thong KL. Genetic variation analysis of Vibrio cholerae using multilocus sequencing typing and multi-virulence locus sequencing typing. *Infect Genet Evol*. 2011;11(5):1121-8.

14. Pavón ABI, Maiden MC. Multilocus sequence typing. *Molecular Epidemiology of Microorganisms*: Springer; 2009. p. 129-40.

15. Chan M-S, Maiden MC, Spratt BG. Database-driven multi locus sequence typing (MLST) of bacterial pathogens. *Bioinformatics*. 2001;17(11):1077-83.

16. Urwin R, Maiden MC. Multi-locus sequence typing: a tool for global epidemiology. *Trends Microbiol*. 2003;11(10):479-87.

17. Bakhshi B, Mahmoudi-Aznaveh A, Salimi-Khorashad A. Clonal dissemi- nation of a single Vibrio cholerae O1 biotype El Tor strain in Sistan-Bal- uchestan province of Iran during 2013. *Curr Microbiol*. 2015;71(2):163-9.

18. Bakhshi B, Boustanshenas M, Mahmoudi-aznaveh A. Emergence of V ibrio cholerae O 1 classical biotype in 2012 in Iran. *Lett Appl Microbiol*. 2014;58(2):145-9.

19. Stavric S, Bachanan B. The isolation and identification of V. cholerae O1 and non-O1 from foods. Government of Canada, Health Protection Branch, Ottawa, Polyscience Publication, MFLP-72, Canada. 1995.

20. Choopun N, Louis V, Huq A, Colwell RR. Simple procedure for rapid iden- tification of Vibrio cholerae from the aquatic environment. *Appl Environ Microbiol*. 2002;68(2):995-8.

21. Control CFD, Prevention. Laboratory methods for the diagnosis of Vibrio cholerae. Atlanta, GA. 1994;500:38-67.

22. Chun J, Rivera IN, Colwell RR. Analysis of 16S-23S rRNA intergenic spacer of Vibrio cholerae and Vibrio mimicus for detection of these species. *Gene Probes*: Springer; 2002. p. 171-8.

23. CLSI C. Methods for antimicrobial dilution and disk susceptibility testing of infrequently isolated or fastidious bacteria; approved guideline. Clinical and Laboratory Standards Institute (M45eA2). 2010.

24. Weinstein MP, Limbago B, Patel J, Mathers A, Campeau S, Mazzulli T, et al. M100 performance standards for antimicrobial susceptibility testing. *CLSI*; 2018.

25. vibrio-pfge-protocol [Available from: https://www.cdc.gov/pulsenet/pdf/ vibrio-pfge-protocol-508c.pdf].

26. Tenover F, Arbeit R, Goering R, the Molecular Typing Working Group of the Society for Healthcare Epidemiology of America. How to select and interpret molecular strain typing methods for epidemiological studies of bacterial infections: a review for healthcare epidemiologists. *Infect Control Hosp Epidemiol*. 1997;18(6):426-39.

27. Amin Marashi SM, Bakhshi B, Imani Fooladi AA, Tavakoli A, Sharifnia A, Pourshafie MR. Quantitative expression of cholera toxin mRNA in Vibrio cholerae isolates with different CTX cassette arrangements. *J Med Microbiol*. 2012;61(8):1071-3.

28. Criscuolo A. AlienTrimmer User Guide.

29. Liu Y, Schröder J, Schmidt B. Musket: a multistage k-mer spectrum- based error corrector for Illumina sequence data. *Bioinformatics*. 2013;29(3):308-15.

30. Prjibelski A, Antipov D, Meleshko D, Lapidus A, Korobeynikov A. Using SPAdes de novo assembler. *Curr Protoc Bioinformatics*. 2020;70(1):e102.

31. Hur Y, Chalita M, Ha S-m, Baek J, Chun J. VCGIDB: a database and web resource for the genomic islands from Vibrio Cholerae. *Pathogens*. 2019;8(4):261.

32. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. *J Antimi- crob Chemother*. 2012;67(11):2640-4.

33. 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(7):3895-903.

34. Herrera CM, Crofts AA, Henderson JC, Pingali SC, Davies BW, Trent MS. The Vibrio cholerae VprA-VprB two-component system controls virulence through endotoxin modification. *MBio*. 2014;5(6):e02283-e2314.

35. Matson JS, Livny J, DiRita VJ. A putative Vibrio cholerae two-component system controls a conserved periplasmic protein in response to the antimicrobial peptide polymyxin B. *PLoS ONE*. 2017;12(10):e0186199.

36. Rowe-Magnus DA, Guerout AM, Mazel D. Bacterial resistance evolu- tion by recruitment of super-integron gene cassettes. *Mol Microbiol*. 2002;43(6):1657-69.

37. 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(6364):785-9.

38. Smith AM, Weill F-X, Njamkepo E, Ngomane HM, Ramalwa N, Sekwadi P, et al. Emergence of Vibrio cholerae O1 sequence type 75, South Africa, 2018–2020. *Emerg Infect Dis*. 2021;27(11):2927.

39. Fonseca ÉL, dos Santos FF, Vieira VV, Vicente AC. New qnr gene cassettes associated with superintegron repeats in Vibrio cholerae O1. *Emerg Infect Dis*. 2008;14(7):1129.

40. Zhou Y-Y, Ma L-Y, Yu L, Lu X, Liang W-L, Kan B, et al. Quinolone Resistance Genes and Their Contribution to Resistance in Vibrio cholerae Serogroup O139. *Antibiotics*. 2023;12(2):416.

41. Hajia M, Sohrabi A. In silico characteristics for re-emerging possibil- ity of Vibrio cholerae genotypes in Iran. *New Microbes New Infect*. 2019;31:100577.

42. Mafi M, Goya MM, Hajia M. A five-year study on the epidemiological approaches to cholera in Iran. *Caspian J Intern Med*. 2016;7(3):162.

43. Liao F, Mo Z, Chen M, Pang B, Fu X, Xu W, et al. Comparison and evalu- ation of the molecular typing methods for toxigenic Vibrio cholerae in Southwest China. *Front Microbiol*. 2018;9:905.

44. Vilela FP, Falcão JP. Analysis of the antimicrobial resistance gene frequency in whole-genome sequenced Vibrio from Latin American countries. *J Med Microbiol*. 2021;70(9):001428.

45. Breurec S, FranckT, Njamkepo E, Mbecko J-R, Rauzier J, Sanke-Waigana H, et al. Seventh pandemic vibrio cholerae O1 sublineages, Central African Republic. *Emerg Infect Dis*. 2021;27(1):262.

46. Greig DR, Schaefer U, Octavia S, Hunter E, Chattaway MA, Dallman TJ, et al. Evaluation of whole-genome sequencing for identification and typing of Vibrio cholerae. *J Clin Microbiol*. 2018;56(11):10-128. https://doi.org/10. 1128/jcm.00831-18.

47. Luo Y, Octavia S, Jin D, Ye J, Miao Z, Jiang T, et al. US Gulf-like toxigenic O1 Vibrio cholerae causing sporadic cholera outbreaks in China. *J Infect*. 2016;72(5):564-72.

48. Okada K, Roobthaisong A, Swaddiwudhipong W, Hamada S, Chantaroj S. Vibrio cholerae O1 isolate with novel genetic background, Thailand- Myanmar. *Emerg Infect Dis*. 2013;19(6):1015.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in pub- lished maps and institutional affiliations.
12. Lukinmaa S, Nakari UM, Eklund M, Siitonen A. Application of molecular genetic methods in diagnostics and epidemiology of food-borne bacterial pathogens. *Apnis*. 2004;112(11-12):908-29.

13. Teh CSJ, Chua KH, Thong KL. Genetic variation analysis of Vibrio cholerae using multilocus sequencing typing and multi-virulence locus sequencing typing. *Infect Genet Evol*. 2011;11(5):1121-8.

14. Pavón ABI, Maiden MC. Multilocus sequence typing. *Molecular Epidemiology of Microorganisms*: Springer; 2009. p. 129-40.

15. Chan M-S, Maiden MC, Spratt BG. Database-driven multi locus sequence typing (MLST) of bacterial pathogens. *Bioinformatics*. 2001;17(11):1077-83.

16. Urwin R, Maiden MC. Multi-locus sequence typing: a tool for global epidemiology. *Trends Microbiol*. 2003;11(10):479-87.

17. Bakhshi B, Mahmoudi-Aznaveh A, Salimi-Khorashad A. Clonal dissemi- nation of a single Vibrio cholerae O1 biotype El Tor strain in Sistan-Bal- uchestan province of Iran during 2013. *Curr Microbiol*. 2015;71(2):163-9.

18. Bakhshi B, Boustanshenas M, Mahmoudi-aznaveh A. Emergence of V ibrio cholerae O 1 classical biotype in 2012 in Iran. *Lett Appl Microbiol*. 2014;58(2):145-9.

19. Stavric S, Bachanan B. The isolation and identification of V. cholerae O1 and non-O1 from foods. Government of Canada, Health Protection Branch, Ottawa, Polyscience Publication, MFLP-72, Canada. 1995.

20. Choopun N, Louis V, Huq A, Colwell RR. Simple procedure for rapid iden- tification of Vibrio cholerae from the aquatic environment. *Appl Environ Microbiol*. 2002;68(2):995-8.

21. Control CFD, Prevention. Laboratory methods for the diagnosis of Vibrio cholerae. Atlanta, GA. 1994;500:38-67.

22. Chun J, Rivera IN, Colwell RR. Analysis of 16S-23S rRNA intergenic spacer of Vibrio cholerae and Vibrio mimicus for detection of these species. *Gene Probes*: Springer; 2002. p. 171-8.

23. CLSI C. Methods for antimicrobial dilution and disk susceptibility testing of infrequently isolated or fastidious bacteria; approved guideline. *Clinical and Laboratory Standards Institute (M45eA2)*. 2010.

24. Weinstein MP, Limbago B, Patel J, Mathers A, Campeau S, Mazzulli T, et al. M100 performance standards for antimicrobial susceptibility testing. *CLSI*; 2018.

25. vibrio-pfge-protocol [Available from: https://www.cdc.gov/pulsenet/pdf/ vibrio-pfge-protocol-508c.pdf].

26. Tenover F, Arbeit R, Goering R, the Molecular Typing Working Group of the Society for Healthcare Epidemiology of America. How to select and interpret molecular strain typing methods for epidemiological studies of bacterial infections: a review for healthcare epidemiologists. *Infect Control Hosp Epidemiol*. 1997;18(6):426-39.

27. Amin Marashi SM, Bakhshi B, Imani Fooladi AA, Tavakoli A, Sharifnia A, Pourshafie MR. Quantitative expression of cholera toxin mRNA in Vibrio cholerae isolates with different CTX cassette arrangements. *J Med Microbiol*. 2012;61(8):1071-3.

28. Criscuolo A. AlienTrimmer User Guide.

29. Liu Y, Schröder J, Schmidt B. Musket: a multistage k-mer spectrum- based error corrector for Illumina sequence data. *Bioinformatics*. 2013;29(3):308-15.

30. Prjibelski A, Antipov D, Meleshko D, Lapidus A, Korobeynikov A. Using SPAdes de novo assembler. *Curr Protoc Bioinformatics*. 2020;70(1):e102.

31. Hur Y, Chalita M, Ha S-m, Baek J, Chun J. VCGIDB: a database and web resource for the genomic islands from Vibrio Cholerae. *Pathogens*. 2019;8(4):261.

32. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. *J Antimi- crob Chemother*. 2012;67(11):2640-4.

33. 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(7):3895-903.

34. Herrera CM, Crofts AA, Henderson JC, Pingali SC, Davies BW, Trent MS. The Vibrio cholerae VprA-VprB two-component system controls virulence through endotoxin modification. *MBio*. 2014;5(6):e02283-e2314.

35. Matson JS, Livny J, DiRita VJ. A putative Vibrio cholerae two-component system controls a conserved periplasmic protein in response to the antimicrobial peptide polymyxin B. *PLoS ONE*. 2017;12(10):e0186199.

36. Rowe-Magnus DA, Guerout AM, Mazel D. Bacterial resistance evolu- tion by recruitment of super-integron gene cassettes. *Mol Microbiol*. 2002;43(6):1657-69.

37. 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(6364):785-9.

38. Smith AM, Weill F-X, Njamkepo E, Ngomane HM, Ramalwa N, Sekwadi P, et al. Emergence of Vibrio cholerae O1 sequence type 75, South Africa, 2018–2020. *Emerg Infect Dis*. 2021;27(11):2927.

39. Fonseca ÉL, dos Santos FF, Vieira VV, Vicente AC. New qnr gene cassettes associated with superintegron repeats in Vibrio cholerae O1. *Emerg Infect Dis*. 2008;14(7):1129.

40. Zhou Y-Y, Ma L-Y, Yu L, Lu X, Liang W-L, Kan B, et al. Quinolone Resistance Genes and Their Contribution to Resistance in Vibrio cholerae Serogroup O139. *Antibiotics*. 2023;12(2):416.

41. Hajia M, Sohrabi A. In silico characteristics for re-emerging possibil- ity of Vibrio cholerae genotypes in Iran. *New Microbes New Infect*. 2019;31:100577.

42. Mafi M, Goya MM, Hajia M. A five-year study on the epidemiological approaches to cholera in Iran. *Caspian J Intern Med*. 2016;7(3):162.

43. Liao F, Mo Z, Chen M, Pang B, Fu X, Xu W, et al. Comparison and evalu- ation of the molecular typing methods for toxigenic Vibrio cholerae in Southwest China. *Front Microbiol*. 2018;9:905.

44. Vilela FP, Falcão JP. Analysis of the antimicrobial resistance gene frequency in whole-genome sequenced Vibrio from Latin American countries. *J Med Microbiol*. 2021;70(9):001428.

45. Breurec S, Franck T, Njamkepo E, Mbecko J-R, Rauzier J, Sanke-Waigana H, et al. Seventh pandemic vibrio cholerae O1 sublineages, Central African Republic. *Emerg Infect Dis*. 2021;27(1):262.

46. Greig DR, Schaefer U, Octavia S, Hunter E, Chattaway MA, Dallman TJ, et al. Evaluation of whole-genome sequencing for identification and typing of Vibrio cholerae. *J Clin Microbiol*. 2018;56(11):10-128. https://doi.org/10. 1128/jcm.00831-18.

47. Luo Y, Octavia S, Jin D, Ye J, Miao Z, Jiang T, et al. US Gulf-like toxigenic O1 Vibrio cholerae causing sporadic cholera outbreaks in China. *J Infect*. 2016;72(5):564-72.

48. Okada K, Roobthaisong A, Swaddiwudhipong W, Hamada S, Chantaroj S. Vibrio cholerae O1 isolate with novel genetic background, Thailand- Myanmar. *Emerg Infect Dis*. 2013;19(6):1015.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.