Assessing the Genetic Diversity of Austrian Corynebacterium diphtheriae Clinical Isolates, 2011 to 2019

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.

Analysis of the HbpA Protein from Corynebacterium diphtheriae Clinical Isolates and Identification of a Putative Hemoglobin-Binding Site on HbpA

The Corynebacterium diphtheriae hemoglobin-binding protein HbpA is critical for the acquisition of iron from the hemoglobin-haptoglobin complex (Hb-Hp). Previous studies using C. diphtheriae strain 1737 showed that large aggregates formed by HbpA are associated with iron transport activity and enhanced binding to Hb-Hp; however, specific regions within HbpA required for Hb-Hp binding or iron uptake have not been identified. In this study, we characterized two clinical isolates from Austria, designated 07-18 and 09-15, which express HbpA proteins that share only 53% and 44% sequence identity, respectively, to the strain 1737 HbpA protein. The HbpA proteins expressed by the Austrian strains had functional and structural properties similar to those of the HbpA protein in strain 1737 despite the limited sequence similarity. These shared characteristics between the HbpA proteins included similar cellular localization, aggregate formation, and Hb and Hb-Hp binding. Additionally, the Austrian strains were able to acquire iron from Hb and Hb-Hp, and deletion of the hbpA gene from these two clinical isolates reduced their ability to use Hb-Hp as an iron source. A sequence comparison between the HbpA proteins from 1737 and the Austrian strains assisted in the identification of a putative Hb-binding site that shared similar characteristics with the Hb-binding regions in Staphylococcus aureus NEAT domains. Amino acid substitutions within this conserved Hb-binding region significantly reduced Hb and Hb-Hp binding and diminished the hemin-iron uptake function of HbpA. These findings represent important advances in our understanding of the interaction of HbpA with human hemoproteins. IMPORTANCE Hemoglobin (Hb) is the primary source of iron in humans, and the acquisition of hemin-iron from Hb is critical for many bacterial pathogens to infect and survive in the human host. In this study, we have examined the C. diphtheriae Hb-binding protein HbpA in two clinical isolates and show that these proteins, despite limited sequence similarity, are functionally equivalent to the previously described HbpA protein in strain 1737. A sequence comparison between these three strains led to the identification of a conserved Hb-binding site, which will further our understanding of how this novel protein functions in hemin-iron transport and, more generally, will expand our knowledge on how Hb interacts with proteins.

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.