Genomic evolution and rearrangement of CTX-Φ prophage elements in Vibrio cholerae during the 2018-2022 cholera outbreaks in The Democratic Republic of Congo

Leonid M. Irenge, Jérôme Ambroise, Bertrand Bearzatto, Jean-François Durant, Louisette K. Wimba, Jean-Luc Gala
doi: https://doi.org/10.1101/2024.03.30.24305082

Abstract

Between 2018 and 2022, we conducted whole-genome sequencing and phylogenomic analysis on 173 Vibrio cholerae bacteria (172 V. cholerae O1 and 1 V. cholerae non-O1/non-O139 isolates) from cholera patients across four provinces in the Democratic Republic of Congo (North-Kivu, South-Kivu, Tanganyika, and Kasai Oriental). All V. cholerae O1 isolates were classified into the AFR10d and AFR10e sublineages of AFR10 lineage, originating from the third wave of the seventh El Tor cholera pandemic (7PET). Compared to the strains analysed between 2014 and 2017, we observed only limited genetic changes in the core genome of both sublineages. However, AFR10e expanded across the four provinces, whereas AFR10d appeared to become extinct by the end of 2020. The V. cholerae isolates from 2022 exhibited significant rearrangement in the CTX prophage and associated phage satellites. Notably, this included the tandem repeat of a novel environmental satellite phage RS1 downstream the ctxB toxin gene of the CTX-Φ-3 prophage on the large chromosome, and two tandemly arrayed copies of the pre-CTX-Φ prophage precursor on the small chromosome. In conclusion, while the core genome of V. cholerae O1 AFR10d and AFR10e sublineages showed minimal changes, significant alterations in the content and organisation of elements associated with the CTX-Φ and pre-CTX prophages were identified in AFR10e V. cholerae O1 isolated in 2022.

Introduction

Cholera, a severe watery diarrheal disease caused by Vibrio cholerae serogroups O1 and/or O139 (1, 2) remains a significant healthcare concern in developing countries (1, 3, 4). It is estimated to cause up to 3 million cases of cholera and approximately 95,000 deaths annually (5, 6). Over the past decade, Sub-Saharan Africa has become the most afflicted by cholera (7, 8), now referred to as the “new cholera homeland” (9), with countries in this region reporting 83% of global cases between 2000 and 2015 (10). The DRC is among the leading countries in terms of cholera cases, contributing 5-14% of the global case count annually (5, 11, 12). In 2020, DRC reported 19,789 cases, ranking second only to Yemen in terms of numbers (13). Recent data from 2023 indicate a deepening cholera crisis in the country, with over 41,000 cases and 314 deaths reported (14).

Cholera has traditionally been most prevalent in the Great Lakes Region (GLR) (16), but recent reports show it has spread to most DRC provinces, including areas well beyond the GLR (17). For example, between 2018-2022, Kasai Oriental province in the central region of DRC reported 2,176 cholera cases (18). The initial introduction of V. cholerae O1 strains into the DRC occurred during the first wave of P7ET, known as the T5 introduction (1970–1972) (19, 20). A subsequent introduction, in the early nineties, part of the P7ET’s third wave and designated as the T10 introduction event, brought a strain from southern Asia (21), which has persisted in the GRL through successive outbreaks (22, 23). Previous studies, including our own, have documented the division of the T10 (AFR10) lineage into two sublineages, AFR10d, and AFR10e, and identified a regional focus of endemic disease derived from the T10 introduction (22).

In this study, we conducted whole-genome sequencing and phylogenetic analyses on V. cholerae isolates collected from cholera cases in the North-Kivu, South-Kivu, and Tanganyika between 2018 and 2022, as well as from cases in Kasai Oriental specifically observed in 2020. By comparing these genomes with data from V. cholerae isolates collected between 2014 and 2017, we aim to assess the genetic evolution of V. cholerae in the DRC, and explore potential genetic links between isolates from the GLR and the central Kasai Oriental province.

Material and methods

Ethical clearance

The study protocols received approval from the Ethical Review Board (ERB) of the ISTM/Bukavu, DRC (ISTM-BUKAVU/CRPS/CIES/ML/0016/2023). The ERB granted a waiver of informed consent, determining that the study complied with international guidelines applicable to research during severe outbreaks. To protect confidentiality, samples were analysed anonymously.

Phenotypic analysis

The isolation of V. cholerae bacteria and subsequent phenotypical tests and antimicrobial susceptibility assays followed protocols previously described (22). Rectal swabs were collected from patients exhibiting clinical symptoms of cholera, notably severe watery diarrhoea. These samples were incubated in alkaline peptone water for 6-8 hours, then streaked onto thiosulfate-citrate-bile salts-sucrose (TCBS) agar Petri dishes and incubated at 37°C overnight. Colonies appearing large yellow, flattened with opaque centres and translucent peripheries were further cultured on nutritive alkaline agar.

For phenotypical confirmation of V. cholerae standard microbiological tests were carried out (24), alongside immunological tests to determine the serotypes of the isolates. Antimicrobial susceptibility was assessed using the diffusion method on Mueller-Hinton agar plates, using antibiotic disks. Bacterial isolates were prepared at a density equivalent to 0.5 to 0.6 McFarland standard, and the zone of inhibition around antibiotics discs were measured after 16–20□hours of incubation.

Minimum inhibitory concentration (MIC) for ciprofloxacin was determined with MIC test strips (Liofilchem, Roseto degli Abruzzi, Italy) following the manufacturer’s guidelines. The interpretation of MIC values inhibition zone diameters was based on the EUCAST clinical breakpoint values (version 13.0, https://www.eucast.org/clinical_breakpoints).

Whole-genome sequencing

Genomic DNA (gDNA) was isolated from all V. cholerae isolates in our collection using an EZ1 Advanced XL Biorobot and the Tissue DNA Kit (QIAGEN, Hilden, Germany), employing the Bacterial Card according as per the manufacturer’s instructions. Quantification of the isolated DNA was performed using a Qubit® fluorometer (Thermo Fisher Scientific, Eugene, OR, USA).

For whole-genome sequencing (WGS) of V. cholerae, we initially used Illumina short-reads sequencing. Short-read libraries were prepared from 70 ng DNA using an Illumina DNA Prep kit (Illumina, San Diego, CA, USA), and then sequenced using paired-end (2 x 300 bp) reads on a MiSeq platform (Illumina, San Diego, CA).

To investigate the genomic organisation within the prophage V. cholerae CTX-Φ, WGS was also carried out using the long-reads sequencing on the GridION platform (Oxford Nanopore Technologies, UK) for a subsect of 20 samples, selected across different collection years. Long-reads libraries were prepared from 400 ng of genomic DNA using the rapid barcoding kit 96 V14 (SQK-RBK114) and sequenced on an R10.4.1 flow cell run in a GridION device over a 72-hour period.

Bioinformatics

Raw sequencing data from all V. cholerae isolates were submitted to the European Nucleotide Archive (ENA) under accession number PRJEB66194 (http://www.ebi.ac.uk/ena).

For Illumina short-reads, our bioinformatics pipeline started with quality checking via FastQC v.0.11.9 (25) followed by trimmomatic v.0.39 (26). De novo assembly of paired-end reads for each V. cholerae isolate was performed using Spades v.3.13.0 (27), generating a draft genome sequence. Assembly quality was assessed using QUAST 5.0.2 (28) and completeness was verified with the BUSCO v.4.1.1 algorithm, using the bacteria_odb10 dataset (29).

SNP-based phylogenetic analysis of V. cholerae genomes, including those from our previous (22) and other recent studies, was conducted using kSNP3.0, using the default settings, except for a k-mer size of 19 and the ‘-core’ option. Within identified clusters, representative genomes were selected, and a simplified phylogenomic tree was built using the R package ggtree (31). In-silico Multi-Locus Sequence Typing (MLST) was performed on all genomes by using the screen-Blast-mlst function of the custom Pathogenomics R package (https://github.com/UCL-CTMA/Pathogenomics), requiring 100% sequence identity and coverage for allele type assignment. All genomes were screened for virulence factor genes as per the VFDB database ((http://www.mgc.ac.cn/VFs) (32), with an identity cut-off set at 80%. Genotypic AMR profiles were assessed using AMRFinderPlus v3.10.24 and ARIBA v2.14.6.

Genomic organisation of the CTX-Φ prophage I, V. Cholerae O1 was determined using de novo assembled genomes from GridION long reads with canu v.2.2 (32). Representative sequences of the CTX-Φ prophage were localised using Blastn v.2.12.0 (33) in these contiguous genomes. Finally, the graphical representation of prophage organisation was created using the Gviz v.1.21.1 package in Bioconductor.

Results

V. cholerae serotyping and antimicrobial susceptibility

In this study, 173 clinical isolates of V. cholerae were evaluated. The phenotypic characterisation data and antimicrobial susceptibility profiles are presented in Table 1 (Table near here1).

View this table:
Table 1. Antimicrobial susceptibility patterns and antimicrobial resistance (AMR) genes of V. cholerae isolates from eastern DRC.

AMP: Ampicillin; CHL: Chloramphenicol; CIP: Ciprofloxacin; DOX: Doxycycline; ERY: Erythromycin; PXB: Polymyxin B; SXT: Sulfamethoxazole-Trimethoprim; TET: Tetracycline. Antimicrobial susceptibility profile of isolates is expressed as the percentage of resistant isolates of V. cholerae. AMR genes are in in italic.

Of these isolates, 172 tested positive for agglutination with the O1 antiserum. Among these, 145 isolates agglutinated with the Ogawa and 27 with the Inaba antisera. All 172 V. cholerae isolates O1 were found to be susceptible to ampicillin, chloramphenicol, and tetracycline. However, they showed resistance to cotrimoxazole and polymyxin B, the latter being indicative of the El Tor biotype of V. cholerae (34). While all Inaba isolates (n = 27) were susceptible to ciprofloxacin, the Ogawa isolates (n = 146) demonstrated reduced susceptibility to this antibiotic. The non-O1/non-O139 isolate was susceptible to all antimicrobial agents tested including polymyxin B.

Phylogenetics of V. cholerae O1 isolates in the DRC provinces under study

Phylogenetic analyses revealed that all 172 V. cholerae O1 isolates belonged to the AFR10 lineage (comprising AFR10d and AFR10e sublineages), aligning with previous reports from the same region (22, 23). Consistent with an earlier study, DRC AFR10 V. cholerae O1 isolates were divided into two sublineages: AFR10e and AFR10d, associated with sequence types ST69 and ST515, respectively as shown in Fig. 1 (insert figure 1 near here).

Figure 1.
Figure 1.

Phylogenomic analysis of V. cholerae O1 isolates (2014–2022) from four investigated DRC provinces.

The 7PET V. cholerae O1 biotype EL Tor N19691 was used as an outgroup to root the tree. The scale bar represents substitutions per variable site in the core genome. V. cholerae O1 isolates from DRC are in black, while isolates from other African countries and the isolate N19691 are in grey.

All isolates from the Kasai Oriental (n = 9) and Tanganyika provinces (n = 2) characterised in this study belonged to the AFR10e sublineage (Fig.2) (insert figure 2 near here). From the 62 V. cholerae O1 isolates characterised in South-Kivu province between 2018 and 2022, the vast majority (n = 60) belonged to the AFR10e lineage, whereas only two isolates (collected in 2018) belonged to the AFR10d sublineage. Similarly, in North Kivu, the incidence of AFR10d declined dramatically by the end of 2019, with only one AFR10d isolate identified in 2020, and none reported since. During the 2018-2022 period, there were 73 isolates from the AFR10e and 25 from the AFR10d sublineage. Therefore, in contrast to our observations during the 2014-2017 period, the AFR10e sublineage significantly outpaced its AFR10d counterpart in both North-Kivu and South-Kivu provinces (see Fig. 2).

Figure 2.
Figure 2.

Temporal prevalence of AFR10d and AFR10e V. cholerae O1 isolates in the four DRC provinces under investigation. For the Kasai-Oriental and Tanganyika provinces, strains were collected only in 2020 and 2022, respectively.

AMR genes

The analysis of antimicrobial (AMR) genes revealed that all DRC AFR10d and AFR10e isolates, consistent with previous studies, harboured the SXT/R391 integrative conjugative element ICEVchBan5, which carries genes associated with reduced susceptibility to cotrimoxazole (35) (Table 1). Similar to findings in all isolates of the AFR10 sublineage (23), both AFR10d and AFR10e isolates exhibited the S83I substitution in the gyrA protein.

Additionally, all AFR10e isolates presented the S85L substitution in the parC protein, a combination associated with reduced susceptibility to quinolones in V. cholerae (23). Notably, no AFR10d isolate in this study was found to carry both the S83I substitution in gyrA and the S85L substitution in parC. This observation was also true for clinical V. cholerae isolates from Tanganyika (n = 3), and Kasai Oriental (n = 9), all phylogenetically belonging to the AFR10e sublineage, and carrying both the S83I gyrA, and the S85L parC substitutions. Moreover, all V. cholerae O1 isolates characterised in this study possessed the catB9, and floR AMR genes, which confer resistance to chloramphenicol, the beta-lactamase varG gene, and the sulI gene associated with sulphonamide resistance. Surprisingly, the presence of these AMR genes did not compromise the susceptibility of AFR10 isolates to ampicillin, chloramphenicol, and tetracycline.

Virulence genes

AFR10d and AFR10e isolates analysed in DRC between 2018 and 2022 possessed the same virulence genes as those characterised in the earlier 2014-2017 study (22). These include:

  1. Virulence genes associated with the CTX-Φ-3 prophage and its phage satellites (36, 37), as previously described in all AFR10 isolates (21);

  2. The complete Vibrio pathogenicity island-1 (VPI-1) (38), the Vibrio Seventh Pandemic Island-I (VSP-I), as well as the Vibrio Seventh Pandemic Island (VSP-II) (39). Notably, VSP-II exhibits a deletion spanning from ORF VC_0495 to VC_0512.

Furthermore, no genetic changes were identified in the virulence genes compared to those in the previous study (22).

CTX-Φ prophage organization in O1 isolates

Long reads sequencing of 20 out of 172 V. cholerae O1 isolates using the GridION ONT platform, followed by assembly with Canu, generated contiguous genomes for each isolate. This approach facilitated the comprehensive characterisation of repetitive elements within the CTX-Φ prophage.

Isolates collected between 2018 and 2020 (Fig. 3A) (insert figure 3A near here). harboured a complete CTX-Φ-3 prophage, which includes the RS2 satellite phage and the core CTX-Φ prophage carrying a classical ctxB gene. This prophage was integrated into the large chromosome between the loci VC1465 and the rtxA gene of the RTX toxin locus. In most of these isolates, the CTX-Φ-3 prophage was flanked upstream by one canonical RS1 satellite phage (RS1ca), aligning with observations of several CTX-Φ prophages in atypical El Tor V. cholerae strains (36, 37).

Figure 3A.
Figure 3A.

Representative organisation of the CTX-Φ prophage and its satellite in V. cholerae isolates collected between 2018 and 2020 in the DRC.

Bloc arrows represent genes present in each element of the CTX-Φ-3-phage and indicate the direction of transcription. RS1ca and RS2ca represent RS1 and RS2 canonical satellite phages.

Isolates collected in 2022 displayed significant changes in the organisation of the CTX-Φ prophage array and its satellites on both the large and the small chromosomes. On the large chromosome, one RS1ca satellite phage flanked the CTX-Φ-3 prophage upstream, while two tandem copies of an environmental RS1, named RS1env-DRC, flanked the CTX-3 Φ prophage downstream of the ctxB toxin gene (Fig. 3B) (insert figure 3B near here). The RS1env-DRC variant is characterised by canonical rstA and rstB genes, an rstC gene with a single amino acid mismatch (L64P) compared to the canonical rstC gene in CTX-Φ prophage, and an rstR gene coding for a 86-aa long protein that matches perfectly (100%) with the rstR characterised in the environmental Mozambican V. cholerae IB1617 isolate (40).

Figure 3B.
Figure 3B.

Representative organisation of the CTX-Φ prophage and its satellites in both the large and small chromosomes in V. cholerae isolates collected in 2022 in the DRC.

Canonical sequences for each gene are shown in grey, while environmental sequences are shown in blue.

Another notable finding was observed in DRC V. cholerae O1 bacteria isolated from 2022, specifically on the small chromosome with the presence of two tandemly arrayed copies of a pre-CTX-Φ prophage inserted between VCA0569 and VCA0570 loci (37) (Fig. 3B). This arrangement resembles the configuration of the pre-CTX-Φ prophage on the small chromosome of the V. cholerae O139 VC1374 (41). Despite sharing a similar pre-CTX-Φ prophage configuration and an identical rstR gene (named rstRZHJ) with the V. cholerae O139 VC1374 isolate, the 2022-DRC V. cholerae O1 isolates exhibited several amino acid substitutions in the rstA, rstB and cep proteins compared to both the V. cholerae O139 VC1374, and to V. cholerae non-O1/non-O139 VCE232 (42). These variants in the rstA, rstB, and cep proteins have been designated as rstADRC, rstBDRC, and cepDRC, respectively, to distinguish them from their counterparts in the environment.

Consequently, V. cholerae O1 isolates collected in 2022 demonstrated a unique configuration of the CTX-Φ prophage array and its satellites on the large chromosome (RS1ca - RS2ca - psh - cep - pIIIctx - ace - zot - ctxA - ctxBclassical - RS1env-DRC- RS1env-DRC), and on the small chromosome (a tandem arrangement of rstRZHJ - rstADRC - rstBDRC - psh - cepDRC - pIIIctx - ace - zot). Notably, such a hybrid configuration has not been previously reported in 7PET V. cholerae O1 from the third wave.

Characterization of the non-O1/non-O139 V. cholerae isolate

In this study, a single non-O1/non-O139 V. cholerae bacterium, isolated from a patient residing the city of Fizi, on the west bank of Lake Tanganyika, was characterized. This isolate exhibited a novel allele of the housekeeping pyrC gene, which led to designation as ST 362 by the curator of the MLST database. The combination of this allele with existing alleles for the other housekeeping genes constitutive of the MLST typing scheme (adk = 96, gyrB = 41; mdh = 109; metE = 237; pntA = 136; purM = 1; pyrC = 362) resulted in a new MLST profile, assigned the number 1577.

Similar to other non-O1/O139 V. cholerae isolates previously identified in the DRC (22), this ST1577 isolate lacked most PAI genes characteristic of the 7PET. Nonetheless, it harboured several virulence genes, including those related to Type VI Secretion System (T6SS) - VCA0109, VCA0122, vgrG.2 and vgrG.3 -, the rtx cluster - rtxA, rtxB, rtxC, rtxD genes -, and the virulence-associated (vas) operon.

Discussion

The objective of this study was to phylogenetically characterise V. cholerae isolates associated with cholera outbreaks in DRC from 2018 to 2022, building upon a previous study that analysed isolates from the eastern DRC between 2014 and 2017 (22). Our results are consistent with previous reports that V. cholerae O1 isolates from the third wave and T10 transmission (AFR10) remain the major clade responsible for cholera outbreaks in the eastern part of the country (13, 22, 23). This evidence marks the sustained prevalence of a single clade as unprecedented when compared to the establishment of other V. cholerae O1 clades in neighbouring African countries (43–45).

Furthermore, our findings are in line with recent studies that show an AFR10d decline and the dominance of the AFR10e sublineage (23). In particular, we demonstrate that V. cholerae O1 isolates from Kasaï Oriental province, the first from central DRC to be molecularly characterised, belong to the AFR10e sublineage. These isolates are genetically linked to cholera outbreaks in the GLR provinces, highlighting the geographical spread and persistence of this clade within the DRC.

Despite geographic isolation of the Kasai Oriental province, which lacks direct natural linkage to the GLR through river or lake, the identification of AFR10e isolates in this area suggests potential transmission routes between the GLR and central regions. Human movement, including train travel, could facilitate such transmission. However, AFR10e V. cholerae bacteria isolated in 2022 exhibited significant genetic changes, in particular in the CTX-F prophage configuration. Although rearrangements in CTX-F prophage have often been linked to specific waves of V. cholerae O1 (36), our results show that all analysed V. cholerae O1 isolates still belong to the 10th transmission event.

Noteworthy, an upstream RS1ca satellite and two downstream tandem copies of RS1env-DRC, a novel environmental RS1 satellite, flanked the CTX-F-3 prophage on the large chromosome. Additionally, two tandemly arranged pre-CTX-Fenv prophages containing the rstRZHJ sequence, similar to those described in a V. cholerae O139 isolate (41), were discovered on the small chromosome of these DRC V. cholerae O1 isolates. This rearrangement results in the co-existence of a CTX-3-F prophage, pre-CTX-F prophage, and RS1 satellites with both clinical and environmental origins within the same bacterial cell. Considering the crucial role of RS1 and pre-CTX-F proteins in regulating the replication, morphogenesis, and the virulence of the CTX-F phage (46), these observed genomic rearrangements may significantly influence the evolutionary dynamics of V. cholerae, including replication, excision, and the horizontal transfer of CTX-F prophage elements.

The discovery of these genetic changes in AFR10e isolates collected in 2022 from various provinces across the DRC not only underscores their expansion into new territories but also suggests an enhanced capacity for environmental adaptation. This adaptability may provide them with a competitive advantage over V. cholerae O1 isolates from previous years (2014–2021), potentially facilitating their rapid spread and dominance.

In conclusion, our study not only confirms the continued dominance of AFR10e sublineage in cholera outbreaks within the DRC but also highlights potential transmission routes and reveals significant genetic alterations within the CTX-Φ prophage. These findings underscore the critical need for ongoing monitoring of V. cholerae O1 to better predict and mitigate future outbreaks, emphasising the importance of understanding genetic evolution in the effectiveness of public health responses.

Data Availability

All the data produced in the present work are contained in the manuscript

https://www.eucast.org/clinical_breakpoints

http://www.mgc.ac.cn/VFs

https://github.com/UCL-CTMA/Pathogenomics

Declaration of interest statement

No potential conflicts of interest were reported by the authors.

Funding details

This study was funded by the Belgian Cooperation Agency of the ARES (Académie de Recherche et d’Enseignement Supérieur) under grant COOP-CONV-20-022. The funder had no role in the study design, collection, analysis and interpretation of data, preparation of the manuscript, or the decision to submit the paper for publication.

Acknowledgments

We extend our gratitude to the technical teams of C.D.R.M/AFIA-TELE (South-Kivu Province) and AMI-LABO (North-Kivu province) for their invaluable assistance in the sampling and processing of faecal samples from patients suspected of having cholera, which facilitated the isolation of V. cholerae. Special thanks are due to the laboratories and authors who submitted sequences to the GISAID EpiFlu and the NCBI databases.

Further acknowledgment goes to the EMBL-EBI team (Genome Campus, Hinxton, Cambridgeshire, CB10 1SD, UK) for providing ENA accession numbers, thus enhancing the accessibility of our data. We are particularly grateful to Dr. Sophie Octavia, the curator of the Vibrio cholerae MLST database, for her expertise and assistance in assigning the new MLST profile.

References

  1. 1.↵
    Clemens JD, Nair GB, Ahmed T, et al. Cholera. Lancet. 2017;390(10101):1539–49.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Kanungo S, Azman AS, Ramamurthy T, Deen J, Dutta S. Cholera. Lancet. 2022;399:1429–1440.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Lippi D, Gotuzzo E, Caini S. Cholera. Microbiol Spectr. 2016;4(4).
  4. 4.↵
    Affairs OUNOftCoH. Latin America & The Caribbean Situation Update (12-18 December 2022) OCHA; 2022 [Available from: https://reliefweb.int/report/haiti/latin-america-caribbean-weekly-situation-update-12-18-december-2022-19-december-2022.
  5. 5.↵
    Ali M, Nelson AR, Lopez AL, et al. Updated global burden of cholera in endemic countries. PLoS Negl Trop Dis. 2015;9:e0003832.
    OpenUrlCrossRefPubMed
  6. 6.↵
    Lessler J, Moore SM, Luquero FJ, et al. Mapping the burden of cholera in sub-Saharan Africa and implications for control: an analysis of data across geographical scales. Lancet. 2018;391:1908–1915.
    OpenUrlPubMed
  7. 7.↵
    Mengel MA, Delrieu I, Heyerdahl L, et al. Cholera outbreaks in Africa. Curr Top Microbiol Immunol. 2014;379:117–144.
    OpenUrlPubMed
  8. 8.↵
    Zheng Q, Luquero FJ, Ciglenecki I, et al. Cholera outbreaks in sub-Saharan Africa during 2010-2019: a descriptive analysis. Int J Infect Dis. 2022;122:215–221.
    OpenUrl
  9. 9.↵
    Gaffga NH, Tauxe RV, Mintz ED. Cholera: a new homeland in Africa? Am J Trop Med Hyg. 2007;77:705–713.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    Machado A, Amorim E, Bordalo AA. Major Stressors Favoring Cholera Trigger and Dissemination in Guinea-Bissau (West Africa). Int J Environ Res Public Health. 2021;18.
  11. 11.↵
    Ingelbeen B, Hendrickx D, Miwanda B, van der Sande MAB, Mossoko M, Vochten H, et al. Recurrent Cholera Outbreaks, Democratic Republic of the Congo, 2008-2017. Emerg Infect Dis. 2019;25:856–864.
    OpenUrlCrossRefPubMed
  12. 12.↵
    (WHO) WHO. Disease Outbreak News ; Cholera - Global situation 2022 [Available from: https://www.who.int/emergencies/disease-outbreak-news/item/2022-DON426.
  13. 13.↵
    Alam MT, Mavian C, Paisie TK, et al. Emergence and Evolutionary Response of Vibrio cholerae to Novel Bacteriophage, Democratic Republic of the Congo. Emerg Infect Dis. 2022;28:2482–2490.
    OpenUrl
  14. 14.↵
    Organization WH. A new resolve to eliminate cholera in DRC. 2023.
  15. 15.
    Region WHO-A. Weekly Bulletin on Outbreaks and Other Emergencies-Week 25: 12 June-18 June 2023.
  16. 16.↵
    Bompangue D, Giraudoux P, Piarroux M, Mutombo G, et al. Cholera epidemics, war and disasters around Goma and Lake Kivu: an eight-year survey. PLoS Negl Trop Dis. 2009;3:e436.
    OpenUrlCrossRefPubMed
  17. 17.↵
    Africa W-ROf. Weekly bulletin on outbreaks and other emergencies—week 1 2018 2018 [Available from: http://apps.who.int/iris/bitstream/10665/259809/1/OEW1-2018.pdf.
  18. 18.↵
    Division Provinciale de la Santé du Kasai Oriental RC. Nombre de cas de choléra dans la province du Kasai Oriental. Ministère de la Santé; 2020.
  19. 19.↵
    Ramamurthy T, Mutreja A, Weill FX, et al. Corrigendum: Revisiting the Global Epidemiology of Cholera in Conjunction With the Genomics of Vibrio cholerae. Front Public Health. 2019;7:237.
    OpenUrl
  20. 20.↵
    Mutreja A, Kim DW, Thomson NR, Connor TR, Lee JH, Kariuki S, et al. Evidence for several waves of global transmission in the seventh cholera pandemic. Nature. 2011;477:462–465.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    Weill FX, Domman D, Njamkepo E, et al. Genomic history of the seventh pandemic of cholera in Africa. Science. 2017;358:785–789.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    Irenge LM, Ambroise J, Mitangala P, et al. Genomic analysis of pathogenic isolates of Vibrio cholerae from eastern Democratic Republic of the Congo (2014-2017). PLoS Negl Trop Dis. 2020;14:e0007642.
    OpenUrlCrossRef
  23. 23.↵
    Hounmanou YMG, Njamkepo E, Rauzier J, et al. Genomic Microevolution of Vibrio cholerae O1, Lake Tanganyika Basin, Africa. Emerg Infect Dis. 2023;29:149–153.
    OpenUrl
  24. 24.↵
    Ramamurthy T, Das B, Chakraborty S, et al. Diagnostic techniques for rapid detection of Vibrio cholerae O1/O139. Vaccine. 2020;38 Suppl 1:A73–A82.
    OpenUrl
  25. 25.↵
    Andrews S. FastQC: a quality control tool for high throughput sequence data. Babraham Bioinformatics, Babraham Institute, Cambridge, United Kingdom; 2010.
  26. 26.↵
    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–2120.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    Bankevich A, Nurk S, Antipov D, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. Journal of computational biology. 2012;19:455–477.
    OpenUrlCrossRefPubMed
  28. 28.↵
    Gurevich A, Saveliev V, Vyahhi N, et al. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29:1072–1075.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    Seppey M, Manni M, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness. Gene prediction: methods and protocols. 2019:227–245.
  30. 30.
    Gardner SN, Hall BG. When whole-genome alignments just won’t work: kSNP v2 software for alignment-free SNP discovery and phylogenetics of hundreds of microbial genomes. PloS one. 2013;8:e81760.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Yu G. Using ggtree to visualize data on treelJlike structures. Current protocols in bioinformatics. 2020;69:e96.
    OpenUrlCrossRef
  32. 32.↵
    Chen L, Yang J, Yu J, et al. VFDB: a reference database for bacterial virulence factors. Nucleic acids research. 2005;33:D325–D328.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    Camacho C, Coulouris G, Avagyan V, et al. BLAST+: architecture and applications. BMC bioinformatics. 2009;10:1–9.
    OpenUrlCrossRefPubMed
  34. 34.↵
    Matson JS, Yoo HJ, Hakansson K, et al. Polymyxin B resistance in El Tor Vibrio cholerae requires lipid acylation catalyzed by MsbB. J Bacteriol. 2010;192:2044–2052.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    Wozniak RA, Fouts DE, Spagnoletti M, et al. Comparative ICE genomics: insights into the evolution of the SXT/R391 family of ICEs. PLoS Genet. 2009;5:e1000786.
    OpenUrlCrossRefPubMed
  36. 36.↵
    Kim EJ, Lee D, Moon SH, et al. CTX prophages in Vibrio cholerae O1 strains. J Microbiol Biotechnol. 2014;24:725–731.
    OpenUrl
  37. 37.↵
    Pant A, Das B, Bhadra RK. CTX phage of Vibrio cholerae: Genomics and applications. Vaccine. 2020;38:A7–A12.
    OpenUrl
  38. 38.↵
    Kumar A, Das B, Kumar N. Vibrio Pathogenicity Island-1: The Master Determinant of Cholera Pathogenesis. Front Cell Infect Microbiol. 2020;10:561296.
    OpenUrl
  39. 39.↵
    Dziejman M, Balon E, Boyd D, et al. Comparative genomic analysis of Vibrio cholerae: genes that correlate with cholera endemic and pandemic disease. Proc Natl Acad Sci U S A. 2002;99:1556–1561.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    Choi SY, Lee JH, Kim EJ, Lee HR, et al. Classical RS1 and environmental RS1 elements in Vibrio cholerae O1 El Tor strains harbouring a tandem repeat of CTX prophage: revisiting Mozambique in 2005. J Med Microbiol. 2010;59:302–308.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    Li X, Zhao L, Gao H, et al. A novel pre-CTX prophage in the Vibrio cholerae serogroup O139 strain. Infect Genet Evol. 2020;81:104238.
    OpenUrl
  42. 42.↵
    Verma J, Bag S, Saha B, et al. Genomic plasticity associated with antimicrobial resistance in Vibrio cholerae. Proc Natl Acad Sci U S A. 2019;116:6226–6231.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    Mwaba J, Debes AK, et al. Three transmission events of Vibrio cholerae O1 into Lusaka, Zambia. BMC Infect Dis. 2021;21:570.
    OpenUrl
  44. 44.
    Bwire G, Sack DA, Almeida M, et al. Molecular characterization of Vibrio cholerae responsible for cholera epidemics in Uganda by PCR, MLVA and WGS. PLoS Negl Trop Dis. 2018;12:e0006492.
    OpenUrlCrossRef
  45. 45.↵
    Hounmanou YMG, Leekitcharoenphon P, Kudirkiene E, et al. Genomic insights into Vibrio cholerae O1 responsible for cholera epidemics in Tanzania between 1993 and 2017. PLoS Negl Trop Dis. 2019;13:e0007934.
    OpenUrl
  46. 46.↵
    Davis BM, Kimsey HH, Kane AV, et al. A satellite phage-encoded antirepressor induces repressor aggregation and cholera toxin gene transfer. EMBO J. 2002;21:4240–4249.
    OpenUrlAbstract/FREE Full Text
Back to top
Posted April 01, 2024.
Download PDF
Data/Code
Email
Genomic evolution and rearrangement of CTX-Φ prophage elements in Vibrio cholerae during the 2018-2022 cholera outbreaks in The Democratic Republic of Congo
Leonid M. Irenge, Jérôme Ambroise, Bertrand Bearzatto, Jean-François Durant, Louisette K. Wimba, Jean-Luc Gala
medRxiv 2024.03.30.24305082; doi: https://doi.org/10.1101/2024.03.30.24305082
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area