Whole-Genome Sequencing of Corynebacterium diphtheriae Isolates Recovered from an Inner-City Population Demonstrates the Predominance of a Single Molecular Strain

Phenotypic and molecular characterization and complete genome sequence of a Corynebacterium diphtheriae strain isolated from cutaneous infection in an immunized individual

Diphtheria is an infectious disease potentially fatal that constitutes a threat to global health security, with possible local and systemic manifestations that result mainly from the production of diphtheria toxin (DT). In the present work, we report a case of infection by Corynebacterium diphtheriae in a cutaneous lesion of a fully immunized individual and provided an analysis of the complete genome of the isolate. The clinical isolate was first identified by MALDI-TOF Mass Spectrometry. The commercial strip system and mPCR performed phenotypic and genotypic characterization, respectively. The antimicrobial susceptibility profile was determined by the disk diffusion method. Additionally, genomic DNA was sequenced and analyzed for species confirmation and sequence type (ST) determination. Detection of resistance and virulence genes was performed by comparisons against ResFinder and VFDB databases. The isolate was identified as a nontoxigenic C. diphtheriae biovar Gravis strain. Its genome presented a size of 2.46 Mbp and a G + C content of 53.5%. Ribosomal Multilocus Sequence Typing (rMLST) allowed the confirmation of species as C. diphtheriae with 100% identity. DDH in silico corroborated this identification. Moreover, MLST analyses revealed that the isolate belongs to ST-536. No resistance genes were predicted or mutations detected in antimicrobial-related genes. On the other hand, virulence genes, mostly involved in iron uptake and adherence, were found. Presently, we provided sufficient clinical data regarding the C. diphtheriae cutaneous infection in addition to the phenotypic and genomic data of the isolate. Our results indicate a possible circulation of ST-536 in Brazil, causing cutaneous infection. Considering that cases of C. diphtheriae infections, as well as diphtheria outbreaks, have still been reported in several regions of the world, studies focusing on taxonomic analyzes and predictions of resistance genes may help to improve the diagnosis and to monitor the propagation of resistant clones. In addition, they can contribute to understanding the association between variation in genetic factors and resistance to antimicrobials.

Genomic characterization of cocirculating Corynebacterium diphtheriae and non-diphtheritic Corynebacterium species among forcibly displaced Myanmar nationals, 2017-2019

Respiratory diphtheria is a serious infection caused by toxigenic Corynebacterium diphtheriae, and disease transmission mainly occurs through respiratory droplets. Between 2017 and 2019, a large diphtheria outbreak among forcibly displaced Myanmar nationals densely settled in Bangladesh was investigated. Here we utilized whole-genome sequencing (WGS) to characterize recovered isolates of C. diphtheriae and two co-circulating non-diphtheritic Corynebacterium (NDC) species - C. pseudodiphtheriticum and C. propinquum. C. diphtheriae isolates recovered from all 53 positive cases in this study were identified as toxigenic biovar mitis, exhibiting intermediate resistance to penicillin, and formed four phylogenetic clusters circulating among multiple refugee camps. Additional sequenced isolates collected from two patients showed co-colonization with non-toxigenic C. diphtheriae biovar gravis, one of which exhibited decreased susceptibility to the first-line antibiotics and harboured a novel 23-kb multidrug resistance plasmid. Results of phylogenetic reconstruction and virulence-related gene contents of the recovered NDC isolates indicated they were likely commensal organisms, though 80.4 %(45/56) were not susceptible to erythromycin, and most showed high minimum inhibition concentrations against azithromycin. These results demonstrate the high resolution with which WGS can aid molecular investigation of diphtheria outbreaks, through the quantification of bacterial genetic relatedness, as well as the detection of virulence factors and antibiotic resistance markers among case isolates.

Genomic Epidemiology of Corynebacterium diphtheriae in New Caledonia

An increasing number of isolations of Corynebacterium diphtheriae has been observed in recent years in the archipelago of New Caledonia. We aimed to analyze the clinical and microbiological features of samples with C. diphtheriae. All C. diphtheriae isolates identified in New Caledonia from May 2015 to May 2019 were included. For each case, a retrospective consultation of the patient files was conducted. Antimicrobial susceptibility phenotypes, tox gene and diphtheria toxin expression, biovar, and the genomic sequence were determined. Core genome multilocus sequence typing (cgMLST), 7-gene MLST, and search of genes of interest were performed from genomic assemblies. Fifty-eight isolates were included, with a median age of patients of 28 years (range: 9 days to 78 years). Cutaneous origin accounted for 51 of 58 (87.9%) isolates, and C. diphtheriae was associated with Staphylococcus aureus and/or Streptococcus pyogenes in three-quarters of cases. Half of cases came either from the main city Noumea (24%, 14/58) or from the sparsely populated island of Lifou (26%, 15/58). Six tox-positive isolates were identified, associated with recent travel to Vanuatu; 5 of these cases were linked and cgMLST confirmed recent transmission. Two cases of endocarditis in young female patients with a history of rheumatic fever involved tox-negative isolates. The 58 isolates were mostly susceptible to commonly used antibiotics. In particular, no isolate was resistant to the first-line molecules amoxicillin or erythromycin. Resistance to tetracycline was found in a genomic cluster of 17 (29%) isolates, 16 of which carried the tetO gene. There were 13 cgMLST sublineages, most of which were also observed in the neighboring country Australia. Cutaneous infections may harbor nontoxigenic C. diphtheriae isolates, which circulate largely silently in nonspecific wounds. The possible introduction of tox-positive strains from a neighboring island illustrates that diphtheria surveillance should be maintained in New Caledonia, and that immunization in neighboring islands must be improved. Genomic sequencing uncovers how genotypes circulate locally and across neighboring countries. IMPORTANCE The analysis of C. diphtheriae from the tropical archipelago of New Caledonia revealed a high genetic diversity with sublineages that may be linked to Polynesia, Australia, or metropolitan France. Genomic typing allowed confirming or excluding suspected transmission events among cases and contacts. A highly prevalent tetracycline-resistant sublineage harboring the tetO gene was uncovered. Toxigenic isolates were observed from patients returning from Vanuatu, showing the importance of improving vaccination coverage in settings where it is insufficient. This study also illustrates the importance for diphtheria surveillance of the inclusion of isolates from cutaneous sources in addition to respiratory cases, in order to provide a more complete epidemiological picture of the diversity and transmission of C. diphtheriae.

Genomic Analysis of Corynebacterium diphtheriae Strains Isolated in the Years 2007-2022 with a Report on the Identification of the First Non-Toxigenic tox Gene-Bearing Strain in Poland

Infections caused by non-toxigenic Corynebacterium diphtheriae have been reported every year in Poland since 2004, with the ST8 biovar gravis strains being most commonly isolated. This study analyzed thirty strains isolated between 2017 and 2022 and six previously isolated strains. All the strains were characterized using classic methods in terms of species, biovar level, and diphtheria toxin production, as well as by means of whole genome sequencing. The phylogenetic relationship based on SNP analysis was determined. The number of C. diphtheriae infections has been rising in Poland every year with a maximum of 22 cases in the year 2019. Since 2022, only the non-toxigenic gravis ST8 (most common) and mitis ST439 (less common) strains have been isolated. An analysis of the genomes of the ST8 strains showed that they had many potential virulence factors, such as adhesins and iron-uptake systems. The situation rapidly changed in 2022 and strains from different STs were isolated (ST32, 40, and 819). The ST40 biovar mitis strain was found to be non-toxigenic tox gene-bearing (NTTB), with the tox gene inactivated due to a single nucleotide deletion. Such strains were previously isolated in Belarus. The sudden appearance of new C. diphtheriae strains with different STs and the isolation of the first NTTB strain in Poland indicate that C. diphtheriae should be classified as a pathogen of special public health concern.

Endemic erythromycin resistant Corynebacterium diphtheriae in Vietnam in the 1990s

Diphtheria is a potentially fatal respiratory disease caused by toxigenic forms of the Gram-positive bacterium Corynebacterium diphtheriae. Despite the availability of treatments (antitoxin and antimicrobials) and effective vaccines, the disease still occurs sporadically in low-income countries and in higher income where use of diphtheria vaccine is inconsistent. Diphtheria was highly endemic in Vietnam in the 1990s; here, we aimed to provide some historical context to the circulation of erythromycin resistant organisms in Vietnam during this period. After recovering 54 C. diphtheriae isolated from clinical cases of diphtheria in Ho Chi Minh City between 1992 and 1998 we conducted whole genome sequencing and analysis. Our data outlined substantial genetic diversity among the isolates, illustrated by seven distinct Sequence Types (STs), but punctuated by the sustained circulation of ST67 and ST209. With the exception of one isolate, all sequences contained the tox gene, which was classically located on a corynebacteriophage. All erythromycin resistant isolates, accounting for 13 % of organisms in this study, harboured a novel 18 kb erm(X)-carrying plasmid, which exhibited limited sequence homology to previously described resistance plasmids in C. diphtheriae. Our study provides historic context for the circulation of antimicrobial resistant C. diphtheriae in Vietnam; these data provide a framework for the current trajectory in global antimicrobial resistance trends.

Genomic Epidemiology and Strain Taxonomy of Corynebacterium diphtheriae

Corynebacterium diphtheriae is highly transmissible and can cause large diphtheria outbreaks where vaccination coverage is insufficient. Sporadic cases or small clusters are observed in high-vaccination settings. The phylogeography and short timescale evolution of C. diphtheriae are not well understood, in part due to a lack of harmonized analytical approaches of genomic surveillance and strain tracking. We combined 1,305 genes with highly reproducible allele calls into a core genome multilocus sequence typing (cgMLST) scheme. We analyzed cgMLST gene diversity among 602 isolates from sporadic clinical cases, small clusters, or large outbreaks. We defined sublineages based on the phylogenetic structure within C. diphtheriae and strains based on the highest number of cgMLST mismatches within documented outbreaks. We performed time-scaled phylogenetic analyses of major sublineages. The cgMLST scheme showed high allele call rate in C. diphtheriae and the closely related species C. belfantii and C. rouxii. We demonstrate its utility to delineate epidemiological case clusters and outbreaks using a 25 mismatches threshold and reveal a number of cryptic transmission chains, most of which are geographically restricted to one or a few adjacent countries. Subcultures of the vaccine strain PW8 differed by up to 20 cgMLST mismatches. Phylogenetic analyses revealed a short-timescale evolutionary gain or loss of the diphtheria toxin and biovar-associated genes. We devised a genomic taxonomy of strains and deeper sublineages (defined using a 500-cgMLST-mismatch threshold), currently comprising 151 sublineages, only a few of which are geographically widespread based on current sampling. The cgMLST genotyping tool and nomenclature was made publicly accessible (https://bigsdb.pasteur.fr/diphtheria). Standardized genome-scale strain genotyping will help tracing transmission and geographic spread of C. diphtheriae. The unified genomic taxonomy of C. diphtheriae strains provides a common language for studies of ecology, evolution, and virulence heterogeneity among C. diphtheriae sublineages.

Phenotypic and Genotypic Correlates of Penicillin Susceptibility in Nontoxigenic Corynebacterium diphtheriae, British Columbia, Canada, 2015-2018

In 2015, the Clinical and Laboratory Standards Institute (CLSI) updated its breakpoints for penicillin susceptibility in Corynebacterium species from <1 mg/L to <0.12 mg/L. We assessed the effect of this change on C. diphtheriae susceptibility reported at an inner city, tertiary care center in Vancouver, British Columbia, Canada, during 2015–2018 and performed whole-genome sequencing to investigate phenotypic and genotypic resistance to penicillin. We identified 44/45 isolates that were intermediately susceptible to penicillin by the 2015 breakpoint, despite meeting previous CLSI criteria for susceptibility. Sequencing did not reveal β-lactam resistance genes. Multilocus sequence typing revealed a notable predominance of sequence type 76. Overall, we saw no evidence of penicillin nonsusceptibility at the phenotypic or genotypic level in C. diphtheriae isolates from our institution. The 2015 CLSI breakpoint change could cause misclassification of penicillin susceptibility in C. diphtheriae isolates, potentially leading to suboptimal antimicrobial treatment selection.

Genomic epidemiology of nontoxigenic Corynebacterium diphtheriae from King County, Washington State, USA between July 2018 and May 2019

Between July 2018 and May 2019, Corynebacterium diphtheriae was isolated from eight patients with non-respiratory infections, seven of whom experienced homelessness and had stayed at shelters in King County, WA, USA. All isolates were microbiologically identified as nontoxigenic C. diphtheriae biovar mitis. Whole-genome sequencing confirmed that all case isolates were genetically related, associated with sequence type 445 and differing by fewer than 24 single-nucleotide polymorphisms (SNPs). Compared to publicly available C. diphtheriae genomic data, these WA isolates formed a discrete cluster with SNP variation consistent with previously reported outbreaks. Virulence-related gene content variation within the highly related WA cluster isolates was also observed. These results indicated that genome characterization can readily support epidemiology of nontoxigenic C. diphtheriae.