Backward compatibility of whole genome sequencing data with MLVA typing using a new MLVAtype shiny application for Vibrio cholerae

Vibrio cholerae O1 associated with recent endemic cholera shows temporal changes in serotype, genotype, and drug-resistance patterns in Bangladesh

BACKGROUND

Despite the advancement in our understanding of cholera and its etiological agent, Vibrio cholerae, the prevention and treatment of the disease are often hindered due to rapid changes in drug response pattern, serotype, and the major genomic islands namely, the CTX-prophage, and related genetic characteristics. In the present study, V. cholerae (n = 172) associated with endemic cholera in Dhaka during the years 2015-2021 were analyzed for major phenotypic and genetic characteristics, including drug resistance patterns.

RESULTS

Results revealed that the V. cholerae strains belonged to serogroup O1 biotype El Tor carrying El Tor -specific genes rtxC, tcpA El Tor, and hlyA El Tor, but possessed classical-biotype cholera toxin. Serotypes of V. cholerae strains differed temporally in predominance with Inaba during 2015-2017, and again in 2020-2021, while Ogawa was the predominant serotype in 2018-2019. Also, ctxB1 was predominant in V. cholerae associated with cholera during 2015-2017, while ctxB7 was predominant in 2018, and in the subsequent years, as observed until 2021. V. cholerae strains differed in their antibiotic resistance pattern with a majority (97%) being multi-drug resistant (MDR) and belonging to six sub-groups. Notably, one of these MDR strains was resistant to eleven of the eighteen antibiotics tested, with resistance to fourth-generation cephalosporin (cefepime), and aztreonam. This extreme drug resistant (XDR) strain carried resistance-related genes namely, extended-spectrum β-lactamases (ESBL), blaOXA-1 and blaPER-3.

CONCLUSION

The observed temporal switching of serotypes, as well as the ctxB genotype, and the emergence of MDR/XDR V. cholerae and their association with endemic cholera in Dhaka underscore the need for routine monitoring of the pathogen for proper patient management.

Multilocus Variable-Number Tandem-Repeat Analysis of Enterohemorrhagic Escherichia coli Serogroups O157, O26, and O111 Based on a De Novo Look-Up Table Constructed by Regression Analysis

Multilocus variable-number tandem-repeat analysis (MLVA) is a widely accepted molecular typing tool for enterohemorrhagic Escherichia coli (EHEC). However, ensuring the accuracy of MLVA data among multiple laboratories remains difficult. We developed a method of constructing adjusted look-up tables, which are necessary for MLVA profiling, at each laboratory using a regression analysis based on electrophoresis data from 24 in-house reference strains. On performing MLVA against 51 EHEC O157 isolates, the repeat numbers of 46 isolates were determined accurately using the look-up table with a 99% prediction interval, an outcome superior to that when using a 95% prediction interval. For the remaining five isolates, although the electrophoresis size fell outside the look-up table, we were able to predict the repeat number accurately by extrapolation or the nearest values of the look-up table. Our approach provides more accurate results than a nonadjusted conventional look-up table for calibrating MLVA profiles.

Genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan

Tularemia is a highly dangerous zoonotic infection due to the bacteria Francisella tularensis. Low genetic diversity promoted the use of polymorphic tandem repeats (MLVA) as first-line assay for genetic description. Whole genome sequencing (WGS) is becoming increasingly accessible, opening the perspective of a time when WGS might become the universal genotyping assay. The main goal of this study was to describe F. tularensis strains circulating in Kazakhstan based on WGS data and develop a MLVA assay compatible with in vitro and in silico analysis. In vitro MLVA genotyping and WGS were performed for the vaccine strain and for 38 strains isolated in Kazakhstan from natural water bodies, ticks, rodents, carnivores, and from one migratory bird, an Isabellina wheatear captured in a rodent burrow. The two genotyping approaches were congruent and allowed to attribute all strains to two F. tularensis holarctica lineages, B.4 and B.12. The seven tandem repeats polymorphic in the investigated strain collection could be typed in a single multiplex PCR assay. Identical MLVA genotypes were produced by in vitro and in silico analysis, demonstrating full compatibility between the two approaches. The strains from Kazakhstan were compared to all publicly available WGS data of worldwide origin by whole genome SNP (wgSNP) analysis. Genotypes differing at a single SNP position were collected within a time interval of more than fifty years, from locations separated from each other by more than one thousand kilometers, supporting a role for migratory birds in the worldwide spread of the bacteria.

Genotyping of Francisella tularensis has become a routine practice in epidemiology. Despite rapidly accumulating knowledge, the phylogeography of the pathogen is still poorly understood and discussions about geographic and temporal origins continue. One important reason is the poor characterization of the pathogen in many tularemia-endemic countries. This article describes the genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan using tandem repeat polymorphisms as well as whole genome sequencing. Thirty-nine strains were analyzed and two lineages were identified, namely B.4 and B.12. The study demonstrates a wider distribution of genotype B.4 in Asia, and identified a more basal branching point in this subclade. The obtained data support the Asian origin hypothesis for F. tularensis. The finding of identical genotypes in strains separated in time by decades and a thousand-kilometers geographic distance, confirms the ability of the bacteria for long-term preservation and fast long distances spread. The isolation of F. tularensis subsp. holarctica from the bird species Isabellina wheatear allows speculating about a major contribution of birds to the phylogeography of the pathogen. A genotyping protocol was developed utilizing seven polymorphic tandem repeats, two of which were identified within the framework of this work. The in vitro and in silico results are identical when using sequencing reads of 300 base-pairs or more.

Genomic analysis of pathogenic isolates of Vibrio cholerae from eastern Democratic Republic of the Congo (2014-2017)

BACKGROUND

Over the past recent years, Vibrio cholerae has been associated with outbreaks in sub-Saharan Africa, notably in Democratic Republic of the Congo (DRC). This study aimed to determine the genetic relatedness of isolates responsible for cholera outbreaks in eastern DRC between 2014 and 2017, and their potential spread to bordering countries.

METHODS/PRINCIPAL FINDINGS

Phenotypic analysis and whole genome sequencing (WGS) were carried out on 78 clinical isolates of V. cholerae associated with cholera in eastern provinces of DRC between 2014 and 2017. SNP-based phylogenomic data show that most isolates (73/78) were V. cholerae O1 biotype El Tor with CTX-3 type prophage. They fell within the third transmission wave of the current seventh pandemic El Tor (7PET) lineage and were contained in the introduction event (T)10 in East Africa. These isolates clustered in two sub-clades corresponding to Multiple Locus Sequence Types (MLST) profiles ST69 and the newly assigned ST515, the latter displaying a higher genetic diversity. Both sub-clades showed a distinct geographic clustering, with ST69 isolates mostly restricted to Lake Tanganyika basin and phylogenetically related to V. cholerae isolates associated with cholera outbreaks in western Tanzania, whereas ST515 isolates were disseminated along the Albertine Rift and closely related to isolates in South Sudan, Uganda, Tanzania and Zambia. Other V. cholerae isolates (5/78) were non-O1/non-O139 without any CTX prophage and no phylogenetic relationship with already characterized non-O1/non-O139 isolates.

CONCLUSIONS/SIGNIFICANCE

Current data confirm the association of both DRC O1 7PET (T)10 sub-clades ST69 and ST515 with recurrent outbreaks in eastern DRC and at regional level over the past 10 years. Interestingly, while ST69 is predominantly a locally endemic sequence type, ST515 became adaptable enough to expand across DRC neighboring countries.