A Dual Barcoding Approach to Bacterial Strain Nomenclature: Genomic Taxonomy of Klebsiella pneumoniae Strains.

Sublineages (SLs) within microbial species can differ widely in their ecology and pathogenicity, and their precise definition is important in basic research and for industrial or public health applications. Widely accepted strategies to define SLs are currently missing, which confuses communication in population biology and epidemiological surveillance. Here, we propose a broadly applicable genomic classification and nomenclature approach for bacterial strains, using the prominent public health threat Klebsiella pneumoniae as a model. Based on a 629-gene core genome multilocus sequence typing (cgMLST) scheme, we devised a dual barcoding system that combines multilevel single linkage (MLSL) clustering and life identification numbers (LINs). Phylogenetic and clustering analyses of >7,000 genome sequences captured population structure discontinuities, which were used to guide the definition of 10 infraspecific genetic dissimilarity thresholds. The widely used 7-gene multilocus sequence typing (MLST) nomenclature was mapped onto MLSL SLs (threshold: 190 allelic mismatches) and clonal group (threshold: 43) identifiers for backwards nomenclature compatibility. The taxonomy is publicly accessible through a community-curated platform (https://bigsdb.pasteur.fr/klebsiella), which also enables external users' genomic sequences identification. The proposed strain taxonomy combines two phylogenetically informative barcode systems that provide full stability (LIN codes) and nomenclatural continuity with previous nomenclature (MLSL). This species-specific dual barcoding strategy for the genomic taxonomy of microbial strains is broadly applicable and should contribute to unify global and cross-sector collaborative knowledge on the emergence and microevolution of bacterial pathogens.

Multidrug-resistant Corynebacterium diphtheriae in people with travel history from West Africa to France, March to September 2023.

We describe 10 unlinked cases of Corynebacterium diphtheriae infection (nine cutaneous, one respiratory) in France in 2023 in persons travelling from Guinea, Mali, Senegal, Niger or Nigeria and Central African Republic. Four isolates were toxigenic. Seven genomically unrelated isolates were multidrug-resistant, including a toxigenic respiratory isolate with high-level resistance to macrolides and beta-lactams. The high rates of resistance, including against first-line agents, call for further microbiological investigations to guide clinical management and public health response in ongoing West African outbreaks.

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.

Multidrug-resistant toxigenic Corynebacterium diphtheriae sublineage 453 with two novel resistance genomic islands.

Antimicrobial therapy is important for case management of diphtheria, but knowledge on the emergence of multidrug-resistance in Corynebacterium diphtheriae is scarce. We report on the genomic features of two multidrug-resistant toxigenic isolates sampled from wounds in France 3 years apart. Both isolates were resistant to spiramycin, clindamycin, tetracycline, kanamycin and trimethoprim-sulfamethoxazole. Genes ermX, cmx, aph(3')-Ib, aph(6)-Id, aph(3')-Ic, aadA1, dfrA15, sul1, cmlA, cmlR and tet(33) were clustered in two genomic islands, one consisting of two transposons and one integron, the other being flanked by two IS6100 insertion sequences. One isolate additionally presented mutations in gyrA and rpoB and was resistant to ciprofloxacin and rifampicin. Both isolates belonged to sublineage 453 (SL453), together with 25 isolates from 11 other countries (https://bigsdb.pasteur.fr/diphtheria/). SL453 is a cosmopolitan toxigenic sublineage of C. diphtheriae, a subset of which acquired multidrug resistance. Even though penicillin, amoxicillin and erythromycin, recommended as the first line in the treatment of diphtheria, remain active, surveillance of diphtheria should consider the risk of dissemination of multidrug-resistant strains and their genetic elements.

Corynebacteria of the diphtheriae Species Complex in Companion Animals: Clinical and Microbiological Characterization of 64 Cases from France.

Corynebacteria of the diphtheriae species complex (CdSC) can cause diphtheria in humans and have been reported from companion animals. We aimed to describe animal infection cases caused by CdSC isolates. A total of 18,308 animals (dogs, cats, horses, and small mammals) with rhinitis, dermatitis, nonhealing wounds, and otitis were sampled in metropolitan France (August 2019 to August 2021). Data on symptoms, age, breed, and the administrative region of origin were collected. Cultured bacteria were analyzed for tox gene presence, production of the diphtheria toxin, and antimicrobial susceptibility and were genotyped by multilocus sequence typing. Corynebacterium ulcerans was identified in 51 cases, 24 of which were toxigenic. Rhinitis was the most frequent presentation (18/51). Eleven cases (6 cats, 4 dogs, and 1 rat) were monoinfections. Large-breed dogs, especially German shepherds (9 of 28 dogs; P < 0.00001), were overrepresented. C. ulcerans isolates were susceptible to all tested antibiotics. tox-positive C. diphtheriae was identified in 2 horses. Last, 11 infections cases (9 dogs and 2 cats; mostly chronic otitis and 2 sores) had tox-negative C. rouxii, a recently defined species. C. rouxii and C. diphtheriae isolates were susceptible to most antibiotics tested, and almost all of these infections were polymicrobial. Monoinfections with C. ulcerans point toward a primary pathogenic potential to animals. C. ulcerans represents an important zoonotic risk, and C. rouxii may represent a novel zoonotic agent. This case series provides novel clinical and microbiological data on CdSC infections and underlines the need for management of animals and their human contacts. IMPORTANCE We report on the occurrence and clinical and microbiological characteristics of infections caused by members of the CdSC in companion animals. This is the first study based on the systematic analysis of a very large animal cohort (18,308 samples), which provides data on the frequency of CdSC isolates in various types of clinical samples from animals. Awareness of this zoonotic bacterial group remains low among veterinarians and veterinary laboratories, among which it is often considered commensal in animals. We suggest that in the case of CdSC detection in animals, the veterinary laboratories should be encouraged to send the samples to a reference laboratory for analysis of the presence of the tox gene. This work is relevant to the development of guidelines in the case of CdSC infections in animals and underlines their public health relevance given the zoonotic transmission risk.

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.

Ongoing diphtheria outbreak in Yemen: a cross-sectional and genomic epidemiology study.

BACKGROUND: An outbreak of diphtheria, declared in Yemen in October, 2017, is ongoing. We did a cross-sectional study to investigate the epidemiological, clinical, and microbiological features of the outbreak. METHODS: Probable cases of diphtheria that were defined clinically and recorded through a weekly electronic diseases early warning system (from 2017, week 22, to 2020, week 17) were used to identify trends of the outbreak (we divided the epidemic into three time periods: May 29, 2017, to June 10, 2018; June 11, 2018, to June 3, 2019; and June 4, 2019, to April 26, 2020). We used the line list of diphtheria reports for governorate-level descriptions. Vaccination coverage was estimated using the 2017 and 2018 annual reports by the national Expanded Programme on Immunization. To confirm cases biologically, Corynebacterium diphtheriae was isolated and identified from throat swabs using standard microbiological culture and identification procedures. We assessed differences in the temporal and geographical distributions of cases, including between different age groups. For in-depth microbiological analysis, tox gene and species-specific rpoB real-time PCR, Illumina genomic sequencing, antimicrobial susceptibility analysis (disk diffusion, E-test), and the Elek diphtheria toxin production test were done on confirmed cases. We used genomic data for phylogenetic analyses and to estimate the nucleotide substitution rate. FINDINGS: The Yemen diphtheria outbreak affected almost all governorates (provinces), with 5701 probable cases and 330 deaths recorded up to April 26, 2020. We collected clinical data for 888 probable cases with throat swab samples referred for biological confirmation, and genomic data for 42 positive cases, corresponding to 43 isolates (two isolates from one culture were included due to distinct colony morphologies). The median age of patients was 12 years (range 0·2-80). The proportion of cases in children aged 0-4 years was reduced during the second time period, after a vaccination campaign, compared with the first period (19% [95% CI 18-21] in the first period vs 14% [12-15] in the second period, p<0·0001). Among 43 tested isolates, 39 (91%) produced the diphtheria toxin and two had low level (0·25 mg/L) antimicrobial resistance to penicillin. We identified six C diphtheriae phylogenetic sublineages, four of which are genetically related to isolates from Saudi Arabia, Eritrea, and Somalia. Inter-sublineage genomic variations in genes associated with antimicrobial resistance, iron acquisition, and adhesion were observed. The predominant sublineage (30 [70%] of 43 isolates) was resistant to trimethoprim and was associated with unique genomic features, more frequent neck swelling (p=0·0029) and a younger age of patients (p=0·060) compared with the other sublineages. Its evolutionary rate was estimated at 1·67 × 10-6 substitutions per site per year, placing its most recent common ancestor in 2015, and indicating silent circulation of C diphtheriae in Yemen before the outbreak was declared. INTERPRETATION: In the Yemen outbreak, C diphtheriae shows high phylogenetic, genomic, and phenotypic variation. Laboratory capacity and real-time microbiological monitoring of diphtheria outbreaks need to be scaled up to inform case management and transmission control of diphtheria. Catch-up vaccination might have provided some protection to the targeted population (children aged 0-4 years). FUNDING: National Centre of the Public Health Laboratories (Yemen), Institut Pasteur, and the French Government Investissement d'Avenir Programme. TRANSLATION: For the Arabic translation of the abstract see Supplementary Materials section.

Corynebacterium rouxii sp. nov., a novel member of the diphtheriae species complex.

A group of six clinical isolates previously identified as Corynebacterium diphtheriae biovar Belfanti, isolated from human cutaneous or peritoneum infections and from one dog, were characterized by genomic sequencing, biochemical analysis and MALDI-TOF mass spectrometry. The six isolates were negative for the diphtheria toxin gene. Phylogenetic analyses showed that the six isolates (including FRC0190T) are clearly demarcated from C. diphtheriae, Corynebacterium belfantii, Corynebacterium ulcerans and Corynebacterium pseudotuberculosis. The average nucleotide identity of FRC0190T with C. diphtheriae NCTC11397T was 92.6%, and was 91.8% with C. belfantii FRC0043T. C. diphtheriae subsp. lausannense strain CHUV2995T appeared to be a later heterotypic synonym of C. belfantii (ANI, 99.3%). Phenotyping data revealed an atypical negative or heterogeneous intermediate maltose fermentation reaction for the six isolates. MALDI-TOF mass spectrometry differentiated the new group from the other Corynebacterium taxa by the presence of specific spectral peaks. rpoB sequences showed identity to atypical, maltose-negative C. diphtheriae biovar Belfanti isolates previously described from two cats in the USA. We propose the name Corynebacterium rouxii sp. nov. for the novel group, with FRC0190T (= CIP 111752T = DSM 110354T) as type strain.

Population genomics and antimicrobial resistance in Corynebacterium diphtheriae.

BACKGROUND: Corynebacterium diphtheriae, the agent of diphtheria, is a genetically diverse bacterial species. Although antimicrobial resistance has emerged against several drugs including first-line penicillin, the genomic determinants and population dynamics of resistance are largely unknown for this neglected human pathogen. METHODS: Here, we analyzed the associations of antimicrobial susceptibility phenotypes, diphtheria toxin production, and genomic features in C. diphtheriae. We used 247 strains collected over several decades in multiple world regions, including the 163 clinical isolates collected prospectively from 2008 to 2017 in France mainland and overseas territories. RESULTS: Phylogenetic analysis revealed multiple deep-branching sublineages, grouped into a Mitis lineage strongly associated with diphtheria toxin production and a largely toxin gene-negative Gravis lineage with few toxin-producing isolates including the 1990s ex-Soviet Union outbreak strain. The distribution of susceptibility phenotypes allowed proposing ecological cutoffs for most of the 19 agents tested, thereby defining acquired antimicrobial resistance. Penicillin resistance was found in 17.2% of prospective isolates. Seventeen (10.4%) prospective isolates were multidrug-resistant (≥ 3 antimicrobial categories), including four isolates resistant to penicillin and macrolides. Homologous recombination was frequent (r/m = 5), and horizontal gene transfer contributed to the emergence of antimicrobial resistance in multiple sublineages. Genome-wide association mapping uncovered genetic factors of resistance, including an accessory penicillin-binding protein (PBP2m) located in diverse genomic contexts. Gene pbp2m is widespread in other Corynebacterium species, and its expression in C. glutamicum demonstrated its effect against several beta-lactams. A novel 73-kb C. diphtheriae multiresistance plasmid was discovered. CONCLUSIONS: This work uncovers the dynamics of antimicrobial resistance in C. diphtheriae in the context of phylogenetic structure, biovar, and diphtheria toxin production and provides a blueprint to analyze re-emerging diphtheria.

Klebsiella pneumoniae carriage in low-income countries: antimicrobial resistance, genomic diversity and risk factors.

Background Klebsiella pneumoniae (hereafter, Kp) is a major public health threat responsible for high levels of multidrug resistant (MDR) human infections. Besides, Kp also causes severe infections in the community, especially in Asia and Africa. Although most Kp infections are caused by endogenous intestinal carriage, little is known about the prevalence and microbiological characteristics of Kp in asymptomatic human carriage, and attached risk factors including environmental sources exposure. Methods Here, 911 pregnant women from communities in Madagascar, Cambodia, and Senegal were screened for gut colonization by Kp. Characteristics of Kp strains (antimicrobial susceptibility, genomic diversity, virulence, and resistance genes) were defined, and associated risk factors were investigated. Results Kp carriage rate was 55.9%, and Kp populations were highly heterogeneous (6 phylogroups, 325 sequence types, Simpson index 99.6%). One third of Kp isolates had acquired antimicrobial resistance genes. MDR-Kp (11.7% to 39.7%) and extended spectrum beta-lactamase (ESBL)-producing Kp (0.7% to 14.7%) varied among countries. Isolates with virulence genes were detected (14.5%). Environmental exposure factors including food, animal contacts, or hospitalization of household members were associated with carriage of Kp, antimicrobial resistance and hypervirulence. However, risk factors were country-specific and Kp subpopulation-specific. Conclusion This large-scale multicenter study uncovers the huge diversity of Kp in human gut carriage, demonstrates that antimicrobial resistance is widespread in communities of three low-income countries, and underlines the challenges posed by Kp colonization to the control of antimicrobial resistance.