<figure><img src="image_1.png" /></figure>

HAL
open science

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

Gabriele Arcari, Melanie Hennart, Edgar Badell, Sylvain Brisse

► To cite this version:

Gabriele Arcari, Melanie Hennart, Edgar Badell, Sylvain Brisse. Multidrug-resistant toxigenic Corynebacterium diphtheriae sublineage 453 with two novel resistance genomic islands. Microbial Genomics, 2023, 9 (1), pp.mgen000923. 10.1099/mgen.0.000923. hal-05127229

HAL Id: hal-05127229

https://hal.science/hal-05127229v1

Submitted on 30 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Multidrug-resistant toxigenic *Corynebacterium diphtheriae* sublineage 453 with two novel resistance genomic islands

Gabriele Arcari¹,², Mélanie Hennart¹,³, Edgar Badell¹,⁴ and Sylvain Brisse¹,⁴,*

## Abstract

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')-lb, 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://bigdb.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.

## DATA SUMMARY

The complete genome sequence data for isolates FRC0137 and FRC0375 were deposited in the European Nucleotide Archive and are available at accessions GCA_902808445.2 and GCA_902808935.2, respectively.

## INTRODUCTION

Toxigenic *Corynebacterium diphtheriae* is responsible for diphtheria, a potentially severe infection of the upper respiratory tract that may include sore throat, neck swelling and distant toxic signs. The pathognomonic symptom is the presence of a pseudomembrane, which adheres on the tonsils, oropharynx and pharynx. If untreated, the disease can potentially cause death from suffocation or due to distal toxic damage [1].

The incidence of diphtheria has dropped significantly following the Global immunization initiative [2]. Nevertheless, outbreaks of diphtheria can still occur when vaccination of populations is sub-optimal, as observed in the ex-Soviet Union countries in the 1990s [3] and more recently in the Rohingya [4, 5] and Yemeni populations [6].

Diphtheria therapy may include the timely administration of diphtheria anti-toxin (DAT), but this product is in short supply [7]. Antimicrobial treatment is central in clinical case management and for contacts' prophylaxis, and typically comprises oral or parenteral penicillin, amoxicillin or erythromycin. Even though antimicrobial resistance is more frequent in some other

Received 01 June 2022; Accepted 03 November 2022; Published 27 January 2023

**Author affiliations:** ¹Institut Pasteur, Université Paris Cité, Biodiversity and Epidemiology of Bacterial Pathogens, Paris, France; ²Department of Molecular Medicine, Sapienza Università di Roma, Rome, Italy; ³Collège doctoral, Sorbonne Université, F-75005 Paris, France; ⁴Institut Pasteur, National Reference Center for Corynebacteria of the Diphtheriae Complex, Paris, France.

*Correspondence: Sylvain Brisse, sylvain.brisse@pasteur.fr

**Keywords:** *Corynebacterium diphtheriae*; diphtheria; toxigenic strains; multidrug-resistance; phylogeography; genomic epidemiology.

**Abbreviations:** AMR, antimicrobial resistance; DAT, diphtheria anti-toxin; GNAT, GCN5-related N-acetyltransferase family; MDR, multidrug-resistant; NRC-CCD, National Reference Center for Corynebacteria of the diphtheriae complex; NTTB, non-toxigenic tox gene-bearing; ONT, Oxford Nanopore Technologies; RGI, resistance genomic island; SL, sublineage.

All sequencing data were deposited in the European Nucleotide Archive and are available within Bioproject PRJEB22103, under biosamples SAMEA6566531 (FRC0137) and SAMEA104208685 (FRC0375).

000923 © 2023 The Authors

<figure><img src="image_3.png" /></figure>
**Impact Statement**

The genetic elements that carry antimicrobial resistance genes in *Corynebacterium diphtheriae* are not well documented and, until now, phylogenetic sublineages of toxigenic strains are considered as being geographically restricted.

In this study we uncovered the existence of a cosmopolitan toxigenic sublineage, a subset of which presented multiple antimicrobial resistance genes clustered in two genomic islands. This work contributes to a better understanding of the dynamics that underlie the spread of antimicrobial resistance in toxigenic *C. diphtheriae*.

*Corynebacterium* species [1], *C. diphtheriae* isolates with reduced susceptibility towards penicillin, cefotaxime and erythromycin as well as tetracycline, sulphonamide or other agents have been reported [2, 8–10].

Yet, knowledge on the epidemiological links and dissemination of resistant *C. diphtheriae* is scarce. *C. diphtheriae* is phylogenetically structured into multiple sublineages, and most of them currently appear geographically restricted [6, 11, 12]. However, rather than true endemism, this may reflect both under-sampling and the lack of systematic genomic sequencing in microbiology laboratories. There are also large gaps in the understanding of the genomic mechanisms underlying antimicrobial resistance (AMR) in this species, and the physical organization of AMR genes within *C. diphtheriae* genomes is largely unknown except for early work on a few strains, e.g. [13]. Some recent studies have reported the presence of transposons [14] or of a plasmid [9] carrying a *pbp2m* transfer unit, sometimes associated with the *ermX* methyltransferase gene. Other members of the *Corynebacterium* genus harbour transposons, insertion sequences and plasmids that carry AMR genes [15–17].

In a previous study of the population structure and antimicrobial-resistance phenotypes of a large collection of *C. diphtheriae* [9], several multidrug-resistant (MDR) isolates were observed; MDR was defined as resistance to antimicrobial agents belonging to three or more categories – excluding fosfomycin to which all *Corynebacterium* are resistant [18]. Here, we investigate further two MDR isolates, FRC0137 and FRC0375. After long-read sequencing and genome-sequence completion, a focus was put on the genomic context of their resistance elements. Furthermore, we aimed to place the two isolates in their evolutionary context by analysing publicly available genomes that belong to the same *C. diphtheriae* sublineage (SL), SL453.

**METHODS**

Isolates FRC0137 and FRC0375 were sampled in 2012 and 2015, respectively. The identification took place at the National Reference Center for Corynebacteria of the *diphtheriae* complex (NRC-CCD) as previously described [9] using PCR assays and MALDI-TOF mass spectrometry (Bruker). The presence of the diphtheria toxin *tox* gene was determined by PCR [9] and confirmed by BLASTN of the DIP_RS12515 reference *tox* gene on the genomic sequences and a more recent real-time PCR assay [19]. The production of the toxin was assessed using the modified Elek test [20].

In addition, the isolates were characterized biochemically for pyrazinamidase, urease, nitrate reductase, for utilization of maltose and for glycogen fermentation. This allowed determining their biovar, based on the combination of nitrate reductase (positive in Mitis and Gravis, negative in Belfanti) and glycogen fermentation (positive in biovar Gravis only).

Phenotypic susceptibility to antimicrobial agents was tested by disc diffusion (BioRad, Marnes-la-Coquette, France) for the antibiotics listed in Table 1. When possible, the results were interpreted according to CA-SFM/EUCAST V.1.0 (January 2019) [21]. In other cases, the interpretative criteria published in Table III of the CA-SFM 2013 recommendations [22] were applied.

The strains were previously [9] sequenced with a NextSeq-500 instrument (Illumina, San Diego, CA, USA) at a minimum of 50×coverage depth; trimming was performed using AlienTrimmer v0.4.0 [23], redundant or over-represented reads were reduced using the khmer software v1.3 [24], sequencing errors were corrected using Musket v1.1 [25] and a *de novo* assembly was performed for each strain using SPAdes v3.12.0 [26]. The quality of the assembly generated by SPAdes was assessed using BUSCO [27]. Here, the two isolates were also subjected to Oxford Nanopore Technologies (ONT) sequencing, for which genomic DNA was extracted using a phenol-chloroform approach. Libraries were prepared using a 1D ligation sequencing kit (SQK-LSK-108) without fragmentation and sequenced using a MinION FLO-MIN-106 flow cell.

ONT and Illumina short reads were combined to generate a hybrid assembly using Unicycler v0.4.4 [28] (normal assembly mode, default parameters).

Genome annotation was performed using RAST [29]. The genomes were screened for resistance and plasmid replicon genes using the ResFinder 4.1 [30] and PlasmidFinder 2.1 [31] online tools at the https://cge.cbs.dtu.dk/services/ website. Synteny between the resistance regions of the FRC0137 isolate and homologous regions from other genomes was analysed by mapping BLASTN alignments with the Circos visualization tool [32]. A neighbor-joining tree was generated based on the concatenated nucleotide
<table><caption>**Table 1.** Antimicrobial resistance phenotypes and genotypes of isolates FRC0137 and FRC0375</caption><thead><tr><th>Antibiotic class</th><th>Antimicrobial agent</th><th>Diameter (mm) of FRC0137</th><th>Diameter (mm) of FRC0375</th><th>Breakpoint (EUCAST 2019) [21]</th><th>Resistance element (FRC0137)</th><th>Resistance element (FRC0375)</th></tr></thead><tbody><tr><td rowspan="5">β lactams</td><td>Penicillin (1 UI)</td><td>17</td><td>20</td><td>29</td><td>None detected</td><td>None detected</td></tr><tr><td>Amoxicillin</td><td>46</td><td>35</td><td>23</td><td>None detected</td><td>None detected</td></tr><tr><td>Oxacillin</td><td>30</td><td>26</td><td>20</td><td>None detected</td><td>None detected</td></tr><tr><td>Cefotaxime</td><td>37</td><td>28</td><td>26</td><td>None detected</td><td>None detected</td></tr><tr><td>Imipenem</td><td>57</td><td>44</td><td>24</td><td>None detected</td><td>None detected</td></tr><tr><td rowspan="3">Macrolides</td><td>Erythromycin</td><td>32</td><td>31</td><td>22</td><td>ermX</td><td>ermX</td></tr><tr><td>Clarithromycin</td><td>42</td><td>33</td><td>22</td><td>ermX</td><td>ermX</td></tr><tr><td>Azithromycin</td><td>21</td><td>27</td><td>22</td><td>ermX</td><td>ermX</td></tr><tr><td rowspan="2">Lincosamides</td><td>Spiramycin</td><td>22</td><td>11</td><td>24</td><td>ermX</td><td>ermX</td></tr><tr><td>Clindamycin</td><td>14</td><td>6</td><td>20</td><td>ermX</td><td>ermX</td></tr><tr><td>Streptogramin</td><td>Pristinamycin</td><td>42</td><td>40</td><td>22</td><td>None detected</td><td>None detected</td></tr><tr><td rowspan="2">Aminoglycosides</td><td>Kanamycin</td><td>6</td><td>6</td><td>17</td><td>aph(3')-Ib, aph(6)-Id, andA1</td><td>aph(3')-Ib, aph(6)-Id, andA1</td></tr><tr><td>Gentamicin</td><td>31</td><td>27</td><td>23</td><td>None detected</td><td>None detected</td></tr><tr><td>Rifamycin</td><td>Rifampicin</td><td>6</td><td>35</td><td>30</td><td>S44F in rpoB*</td><td>None detected</td></tr><tr><td>Tetracycline</td><td>Tetracycline</td><td>15</td><td>18</td><td>24</td><td>tet(33)</td><td>tet(33)</td></tr><tr><td>Quinolone</td><td>Ciprofloxacin</td><td>7</td><td>33</td><td>25</td><td>S90F and D94Y in gyrA†</td><td>None detected</td></tr><tr><td rowspan="2">Folate pathway inhibitors</td><td>Sulfonamide</td><td>6</td><td>6</td><td>17</td><td>sul1</td><td>sul1</td></tr><tr><td>Trimethoprim</td><td>6</td><td>6</td><td>16</td><td>dfrA15</td><td>dfrA15</td></tr><tr><td>Trimethoprim-Sulfamethoxazole</td><td></td><td>6</td><td>6</td><td>19</td><td>sul1, dfrA15</td><td>sul1, dfrA15</td></tr><tr><td>Fosfomycin</td><td>Fosfomycin</td><td>6</td><td>6</td><td>14</td><td>Intrinsic resistance</td><td>Intrinsic resistance</td></tr><tr><td>Glycopeptide</td><td>Vancomycin</td><td>27</td><td>19</td><td>17</td><td>None detected</td><td>None detected</td></tr><tr><td>Oxazolidinone</td><td>Linezolid</td><td>34</td><td>35</td><td>25</td><td>None detected</td><td>None detected</td></tr></tbody></table>

<footnote>*Same of *Corynebacterium urealyticum* in [38].</footnote>

<footnote>†Same position, different substitutions compared with C. urealyticum in [38].</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Schematic representation of the *Cd*-RGI-I and *Cd*-RGI-II resistance islands (a). The upper part represents a map of the approximately 20 kbp, which make up *Cd*-RGI-I of the FRC0137 isolate. Arrows indicate the genes and are colour coded as follows: red for antimicrobial resistance genes, yellow for ISs, green for integrase genes and purple for miscellaneous genes. Coloured horizontal bars below the genes indicate the several parts that compose the resistance genomic island: blue, the fragmented Tn5432; orange, Tn5564; and green, the integron. In the lower part, the synteny with several other isolates retrieved from the GenBank database is shown. Genes are colour-coded as in (a). Specifically, the Tn5432 alignment with a *Corynebacterium striatum* isolate is depicted in blue; the Tn5564 alignment with a *Corynebacterium singulare* is depicted in orange and the one with *Corynebacterium imitans* in yellow. The integron alignments' regions are drawn in green with the IncHI1B plasmid of *Klebsiella pneumoniae* and in aquamarine with the chromosome of *Acinetobacter toowneri*. (b) The upper part represents a map of the approximately 10 kbp, which makes up the *Cd*-RGI-II of isolate FRC0137. Arrows indicate the genes and are colour coded as in (a). In the lower part, is shown a synteny plot with *C. diphtheriae* CIP107572 isolated from Thailand in 1994. Beyond the loss of a small portion of the recombinase gene in FRC0137, these two islands have 100% sequence identity. Note that the genomic context in which the islands are located is different, indicating the horizontal transfer of the island.</figcaption></figure>

sequences of the 1305 genes of the cgMLST scheme [11] of SL453 isolates using the https://bigsdb.pasteur.fr/diphtheria/ platform, and the tree was visualized using iTOL [33].

# RESULTS

## Characteristics of the isolates and description of resistance genomic islands RGI-I and RGI-II

Isolates FRC0137 and FRC0375 were sampled from wound infections of two patients living in a single mainland France city. Strain FRC0137 was isolated in 2012 from a foot abscess in a patient with a recent hospitalization in South Asia. The wound was first treated by administration of amoxicillin-clavulanic acid for 7 days, and later on with vancomycin and ceftriaxone, which was followed by the resolution of the wound infection. Isolate FRC0375 was isolated in 2015 from an ankle ulceration in a patient who had also travelled to the same South Asian country but no earlier than 10 years before the isolation of the strain; in this case the infection was treated with an undefined macrolide and no further details on the outcome were available.

Phenotypic characterization showed that both isolates belonged to biovar Mitis (nitrate positive, glycogen negative) and produced the diphtheria toxin, as determined by Elek's test. Antimicrobial susceptibility testing showed that both isolates were resistant to clindamycin, spiramycin, kanamycin, tetracycline and trimethoprim-sulfamethoxazole; furthermore, FRC0137 was resistant to azithromycin, rifampicin and ciprofloxacin (Table 1). Both isolates were susceptible to erythromycin and clarithromycin even though they harboured the *ermX* gene (see below). They also appeared penicillin-resistant, as do most *C. diphtheriae* using the currently proposed breakpoint [9].

High-quality genome assemblies of FRC0137 and FRC0375 were obtained, with a size of 2549354 bp and of 2613747 bp, and G+C% content of 53.63 and 53.62%, respectively. The FRC0137 genome was assembled into two contigs whereas the genome sequence of FRC0375 was fully circularized. No plasmids were identified. Sequence analysis confirmed the biochemical characterization of toxigenicity and biovar: the two isolates carried an identical, non-disrupted copy of the *tox* gene (allele 2 in the BIGSdb database, https://bigsdb.pasteur.fr/diphtheria/) and were negative for the Gravis biovar-associated *spuA* gene region [9, 34].

The genomes of both isolates carried 12 antimicrobial resistance genes, which we found were grouped on two antimicrobial resistance islands. The first one, herein named *Cd*-RGI-I, had a length of 20769 bp in strain FRC0137 (Fig. 1a) and of 19328 bp in strain FRC0375. The difference in size is due to the loss of the remnants of an *ISCx1* in the latter isolate. The second resistance island was named *Cd*-RGI-II and was of size 8570 bp in both strains (Fig. 1b). *Cd*-RGI-I contained nine
antimicrobial resistance genes, plus a GCN5-related N-acetyltransferase family (GNAT) gene and a truncated version of the *qacE* gene. This resistance island may be divided into three distinct elements (Fig. 1a). The first element corresponds to the transposon Tn5432, which is composed by two IS1249 genes belonging to the IS256 family, and the macrolide resistance gene *ermX* as well as ISCx1. The latter, absent in the FRC0375 strain, may be considered as a genomic scar, since it is prematurely interrupted by a stop codon [17]. Transposon Tn5432 has already been described in other *Corynebacterium* species with some variation [15, 17].

The second element composing this resistance island is transposon Tn5564, which is bounded by an IS5564. It comprises the *cmx* gene coding for a chloramphenicol efflux pump, and an IS6 family member IS1628 gene. Tn5564 also includes the *aph(3')-Ib* and *aph(6)-Id* (formerly known as *strA* and *strB*) streptomycin resistance genes and the *aph(3')-Ic* aminoglycoside phosphotransferase, which was potentially imported with Tn5564 [15, 17].

The third element of *Cd*-RGI-I is composed by an integron flanked by two insertion sequences (IS6100 and IS1628), and comprises four gene cassettes corresponding to AMR genes (*dfrA15*, *cmlA*, *aadA1* and *sul1*), *ΔqacE* and a GNAT gene. Hence, this element may underlie the resistance to trimethoprim, sulfamethoxazole, chloramphenicol and aminoglycosides (*aad1*: streptomycin and spectinomycin). It was not possible to trace back this integron to any previously described *Corynebacterium* or Actinobacteria (GenBank database, last accessed date 7 December 2021). Rather, it may come from Gram-negative species, given that the best sequence matches (maximum 82% coverage and 99% identity) were observed with Gamma-proteobacteria such as *Escherichia coli* (LT985260.1), *Klebsiella pneumoniae* (CP030270.1, KX029332.1) and *Acinetobacter species* (CP046045.1, CP046596.1). In *Cd*-RGI-I, transposon Tn5432 might have been the first element to integrate into the chromosome of an ancestral *C. diphtheriae*, and may then have served as a 'hotspot' for the integration of the two other elements present in the region, as already proposed [35].

The second resistance island (*Cd*-RGI-II) is flanked by two IS6100 insertion sequences and may thus correspond to a transposon; it was identical (100% identity, 100% length coverage) in both isolates. It carries a recombinase gene followed by the *sul1* gene in a tail-to-tail configuration, by a deleted IS481, and by two other AMR genes: *cmlR* and *tet(33)*, potentially contributing to the observed resistance to chloramphenicol and tetracycline (Fig. 1b). We found that this island is present with the same composition but in a different genomic context in a phylogenetically distant (1050 different cgMLST alleles out of 1305 gene loci) *tox*-positive *C. diphtheriae* isolate CIP107572, which belongs to MLST sequence type ST258 and was isolated in Thailand in 1994 [9, 36, 37].

FRC0137 differed from FRC0375 by being resistant to two additional antimicrobial agents: rifampicin and ciprofloxacin. Consistently, in FRC0137 mutations in the genes encoding for the β-subunit of the RNA polymerase (*rpoB*) and for the gyrase subunit A (*gyrA*) were observed. FRC0137 displayed a deduced S444F alteration of the RpoB protein, previously described in a rifampicin-resistant *Corynebacterium urealyticum* [38]. However, this specific mutation did not affect the *rpoB* MLST allele, as it lied outside of the *rpoB* region used in MLST; three additional, synonymous changes were observed within the *rpoB* MLST template, leading to allele rpoB-13 (ST136) compared to FRC0375 *rpoB* sequence (allele 2, ST389).

The amino acid alterations S90F and D94Y deduced from the *gyrA* gene sequence probably underlie the observed ciprofloxacin resistance in FRC0137. *C. urealyticum* isolates displaying S90V and D94Y amino-acid substitutions were also resistant to quinolones [38], and both phenylalanine and valine have a hydrophobic side chain known to favour the disruption of the water-metal ion bridge [39].

## Evolution within sublineage SL453

To understand the evolutionary origins of multidrug-resistance in these isolates, we performed a comparative genomic analysis. According to cgMLST [11], they differed by 131 loci, and thus belong to the same sublineage, SL453 [11]. We found that 25 publicly available genomes also belonged to SL453. These were sampled from 1963 to 2016, nine of them from Malaysia [40], two from France [9, 41], five from Belarus [42], two from Spain [43] and one each from Australia, Guatemala, India, Romania, Russia, Thailand and the USA (Fig. 2b). One isolate had an unknown provenance.

We defined two branches within SL453 based on its phylogenetic structure (Fig. 2). Most SL453 isolates (24/27) carried the *tox* gene even though nine isolates of branch 2, associated with the 1990s outbreak in Belarus (Fig. 2), were defined as non-toxigenic toxin-bearing (NTTB) isolates as deduced from *tox* gene disruption and as previously reported for some of them [42]. Only three isolates of branch 2 harboured an intact copy of the *tox* gene, and three did not carry it. Branch 1 isolates were all hypothesized as being toxigenic, based on their complete *tox* gene, and had diverse geographical and temporal origins: the older isolate was from Thailand, 1996, while eight isolates were from a 2016 outbreak in Malaysia. The three other branch 1 isolates were from Europe, including FRC0137, FRC0375 and isolate CS1 from Spain. CS1 was isolated from a young patient having a travel history in Afghanistan [43] and belonged to the same outbreak cluster (i.e. less than 25 allelic mismatches [11] as FRC0375. Thus, the latter isolate from France is genetically closer from isolate CS1 from Spain, than from FRC0137. Notably, isolate CS1 lacks the *Cd*-RGI-II resistance island, suggesting an initial acquisition of this island in the ancestor of the three isolates, followed by its loss in the isolate from Spain.
<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Phylogeny of *Corynebacterium diphtheriae* sublineage 453, and distribution of isolates features (a). Midpoint-rooted neighbour-joining tree based on the concatenated nucleotide sequences of the 1305 genes of the cgMLST scheme of *C. diphtheriae*. Nodes are colour-coded as follows: empty stars indicate the absence of the *tox* gene, orange stars indicate non-toxigenic toxin bearing (NTTB) isolates, and red stars indicate *tox* positive, toxigenic isolates. Other data are colour-coded according to the keys or, for resistance genes: empty squares indicate the absence of the gene, filled squares indicate the presence of the gene, and a darker colour indicates isolates harbouring two copies of the specific gene (*sul*1 in this case). The scale indicates the number of nucleotide substitutions per site. (b) World map download from https://simplemaps.com/resources/svg-world displaying the country of origin of SL453 isolates. The image was edited using the OpenSource InkScape software.</figcaption></figure>

DISCUSSION

Toxigenic *C. diphtheriae* outbreaks have been reported in recent years [5, 6, 44–46], and occurred in countries with disrupted healthcare services, which do not achieve high vaccination coverage [1]. Antimicrobial treatment is an important addition to diphtheria antitoxin treatment, which is often unavailable and only targets the toxin effects. Therefore, antimicrobial resistance in toxigenic strains is a matter of public health concern.

Here we analysed two MDR toxigenic isolates and discovered two genomic islands that together carry 12 resistance genes and which were previously unreported in *C. diphtheriae*, to our knowledge. Even though most of the resistance genes they carried, except for *ermX*, may confer resistance to antimicrobial agents that are not recommended as he first line in the treatment of *C. diphtheriae* infection or colonization, the resistance profile of these two isolates has public health relevance. First, these genomic regions may in the future aggregate more AMR genes and disseminate them to other sublineages of *C. diphtheriae* or other *Corynebacterium* species. Second, these resistance genes may serve as markers for the yet untraced spread of MDR toxigenic *C. diphtheriae*. The fact that the resistance islands carry genes that confer resistance to agents that are not used against diphtheria, highlights the risk of long-term maintenance of resistance genes, either in the absence of antimicrobial selective pressure, or due to bystander selection of resistance in *C. diphtheriae* through exposure to antibiotics used to treat other infections [47].

The international spread of SL453 isolates and the evolutionary pattern of the RGI-I and RGI-II provide clues as to the scenario of their dissemination. The two isolates from France, and the one from Spain, were linked to South Asia. The acquisition of the resistance islands Cd-RGI-I and Cd-RGI-II might therefore have occurred in this world region. In fact, many MDR isolates previously described in *C. diphtheriae* were linked to this region [10, 12]. Acquisition of the RGIs in SL453 must have occurred after the evolutionary split from the lineage leading to the Malaysian outbreak strain, which does not carry them. We hypothesize that the common ancestor of the three RGI-I carrying isolates had acquired the two RGIs, with subsequent loss of RGI-II in the branch leading to isolate CS1. The three isolates have unexpected genetic relationships, with a close genetic identity of isolate CS1 with isolate FRC0137, which was isolated from a patient in the same French city as FRC0375. The lack of fit observed in this case between epidemiology and genomics illustrates the value of genomic epidemiology and underlines the need to improve the rate of *C. diphtheriae* infection reporting and microbiological characterization to better trace the links between sporadic infections.

Another isolate belonging to SL453, of ST136 (the same ST as FRC0137), was described in Canada in 2011, from a patient who returned from India with a toe infection [48]; and the antimicrobial susceptibility pattern of this isolate was similar to that of FRC0375. No other ST136 isolate was reported from India or elsewhere [12].

The FRC0137 isolate, isolated 3 years before the FRC0375 one, displayed a mutated *gyrA*, indicating that evolution towards quinolone resistance happened in the branch leading to FRC0137, potentially under quinolone selective pressure in Asia [49, 50].
Microevolution of *Corynebacterium* quinolone resistance in defined epidemiological contexts has been observed previously [6], underlying the high risk of emergence of resistance to this antimicrobial class during treatment in *C. diphtheriae* and other *Corynebacterium* species.

A previous study has shown that most sublineages harbouring the *tox* gene have a limited geographical distribution [9]. Our study, instead, uncovers the existence of a cosmopolitan toxigenic MDR clone. We concur with previous statements [11, 12] that *C. diphtheriae* is currently sampled at a largely insufficient rate to decipher its geographical spread with confidence, and it is to be expected that, as *C. diphtheriae* genomic epidemiology develops, more sublineages with broad spatial distribution will be described. Genomic surveillance of *C. diphtheriae* should be strengthened globally, in order to monitor the spread of sublineages of particular concern such as the MDR toxigenic SL453.

**Funding Information**

The National Reference Center for Corynebacteria of the diphtheriae complex receives support from Institut Pasteur and Public Health France (Santé Publique France, Saint Maurice, France). This work was also supported financially by the French Government's Investissement d'Avenir program Laboratoire d'Excellence "Integrative Biology of Emerging Infectious Diseases" (ANR-10-LABX-62-IBEID). M.H. was supported financially by the PhD grant "Codes4strains" from the European Joint Programme One Health, which has received funding from the European Union's Horizon 2020 Research and Innovation Programme under Grant Agreement No. 773830.

**Acknowledgements**

This work used the computational and storage services provided by the IT department at Institut Pasteur. We thank Nora Zidane (Institut Pasteur) for assistance with the genome sequencing of the isolates.

**Author contributions**

S.B. designed and coordinated the study, G.A. and M.H. performed the genomic analyses. E.B. performed the microbiological characterization of the *C. diphtheriae* isolates. S.B. and G.A. wrote the manuscript. All authors provided comments and approved the final manuscript version.

**Conflicts of interest**

The authors declare no conflict of interest.

**Ethical statement**

To conduct the research, we used bacterial strains, which are not considered human samples. Accordingly, this research was not considered human research and is out of the scope of the decree n° 2016-1537 of 16 November 2016 implementing law n° 2012-300 of 5 March 2012 on research involving human subjects. Therefore, no ethics approval was needed and no informed consent was required. All French bacteriological samples were collected in the frame of the French national diphtheria surveillance and are collected, coded, shipped, managed and analyzed according to the French National Reference Center protocols that received approval by French competent body CNIL (Commission Nationale Informatique et Libertés; approval number: 1474671). The research conformed to the principles of the Helsinki Declaration.

**References**

1. Sharma NC, Efstatiou A, Mokrousov I, Mutreja A, Das B, et al. *Diphtheria*. *Nat Rev Dis Primers* 2019;5:81.

2. Galazka A. The changing epidemiology of diphtheria in the vaccine era. *J Infect Dis* 2000;181 Suppl 1:S2-9.

3. Vitek CR, Wharton M. Diphtheria in the former Soviet Union: reemergence of a pandemic disease. *Emerg Infect Dis* 1998;4:539-550.

4. Matsuyama R, Akhmetzhanov AR, Endo A, Lee H, Yamaguchi T, et al. Uncertainty and sensitivity analysis of the basic reproduction number of diphtheria: a case study of a Rohingya refugee camp in Bangladesh, November-December 2017. *PeerJ* 2018;6:e4583.

5. Rahman MR, Islam K. Massive diphtheria outbreak among Rohingya refugees: lessons learnt. *J Travel Med* 2019;26.

6. Badell E, Alharazi A, Criscuolo A, Almoayed KAA, Lefrancq N, et al. Ongoing diphtheria outbreak in Yemen: a cross-sectional and genomic epidemiology study. *Lancet Microbe* 2021;2:e386-e396.

7. Wagner KS, Stickings P, White JM, Neal S, Crowcroft NS, et al. A review of the international issues surrounding the availability and demand for diphtheria antitoxin for therapeutic use. *Vaccine* 2009;28:14-20.

8. Benamrouche N, Hasnaoui S, Badell E, Guettou B, Lazri M, et al. Microbiological and molecular characterization of *Corynebacterium diphtheriae* isolated in Algeria between 1992 and 2015. *Clin Microbiol Infect* 2016;22:1005.

9. Hennart M, Panunzi LG, Rodrigues C, Gaday Q, Baines SL, et al. Population genomics and antimicrobial resistance in *Corynebacterium diphtheriae*. *Genome Med* 2020;12:107.

10. Husada D, Soegianto SDP, Kurniawati IS, Hendrata AP, Irawan E, et al. First-line antibiotic susceptibility pattern of toxigenic *Corynebacterium diphtheriae* in Indonesia. *BMC Infect Dis* 2019;19:1049.

11. Guglielmini J, Hennart M, Badell E, Toubiana J, Criscuolo A, et al. Genomic epidemiology and strain taxonomy of *Corynebacterium diphtheriae*. *J Clin Microbiol* 2021;59:e01581-21.

12. Will RC, Ramamurthy T, Sharma NC, Veeraraghavan B, Sangal L, et al. Spatiotemporal persistence of multiple, diverse clades and toxins of *Corynebacterium diphtheriae*. *Nat Commun* 2021;12:1500.

13. Barraud O, Badell E, Denis F, Guiso N, Ploy MC. Antimicrobial drug resistance in *Corynebacterium diphtheriae* mitis. *Emerg Infect Dis* 2011;17:2078-2080.

14. Forde BM, Henderson A, Playford EG, Looke D, Henderson BC, et al. Fatal respiratory diphtheria caused by β-lactam-resistant *Corynebacterium diphtheriae*. *Clin Infect Dis* 2021;73:e4531-e4538.

15. Tauch A, Kassing F, Kalinowski J, Pühler A. The *Corynebacterium xerosis* composite transposon Tn5432 consists of two identical insertion sequences, designated IS1249, flanking the erythromycin resistance gene ermCX. *Plasmid* 1995;34:119-131.

16. Tauch A, Zheng Z, Pühler A, Kalinowski J. *Corynebacterium striatum* chloramphenicol resistance transposon Tn5564: genetic organization and transposition in *Corynebacterium glutamicum*. *Plasmid* 1998;40:126-139.

17. Leyton B, Ramos JN, Baio PVP, Veras JFC, Souza C, et al. Treat me well or will resist: uptake of mobile genetic elements determine the resistome of *Corynebacterium striatum*. *Int J Mol Sci* 2021;22:7499.

18. Soriano F, Zapardiel J, Nieto E. Antimicrobial susceptibilities of *Corynebacterium* species and other non-spore-forming gram-positive bacilli to 18 antimicrobial agents. *Antimicrob Agents Chemother* 1995;39:208-214.

19. Badell E, Guillot S, Tulliez M, Pascal M, Panunzi LG, et al. Improved quadruplex real-time PCR assay for the diagnosis of diphtheria. *J Med Microbiol* 2019;68:1455-1465.
20. Engler KH, Glushkevich T, Mazurova IK, George RC, Efstratiou A. A modified Elek test for detection of toxigenic corynebacteria in the diagnostic laboratory. *J Clin Microbiol* 1997;35:495-498.

21. The European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters. Version 9.0; 2019. http://www.eucast.org

22. https://www.sfm-microbiologie.org/wp-content/uploads/2020/07/CASFM_2013.pdf.

23. Criscuolo A, Brisse S. AlienTrimmer: a tool to quickly and accurately trim off multiple short contaminant sequences from high-throughput sequencing reads. *Genomics* 2013;102:500-506.

24. Crusoe MR, Alameldin HF, Awad S, Boucher E, Caldwell A, et al. The khmer software package: enabling efficient nucleotide sequence analysis. *J 1000Res* 2015;4:900.

25. Liu Y, Schröder J, Schmidt B. Musket: a multistage k-mer spectrum-based error corrector for Illumina sequence data. *Bioinformatics* 2013;29:308-315.

26. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. *J Comput Biol* 2012;19:455-477.

27. Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. *Bioinformatics* 2015;31:3210-3212.

28. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycer: resolving bacterial genome assemblies from short and long sequencing reads. *PLoS Comput Biol* 2017;13:e1005595.

29. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, et al. The RAST Server: rapid annotations using subsystems technology. *BMC Genomics* 2008;9:75.

30. Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, et al. ResFinder 4.0 for predictions of phenotypes from genotypes. *J Antimicrob Chemother* 2020;75:3491-3500.

31. Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O, et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. *Antimicrob Agents Chemother* 2014;58:3895-3903.

32. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, et al. Circos: an information aesthetic for comparative genomics. *Genome Res* 2009;19:1639-1645.

33. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. *Nucleic Acids Res* 2019;47:W256-W259.

34. Santos AS, Ramos RT, Silva A, Hirata R Jr, Mattos-Guaraldi AL, et al. Searching whole genome sequences for biochemical identification features of emerging and reemerging pathogenic Corynebacterium species. *Funct Integr Genomics* 2018;18:593-610.

35. Nudel K, Zhao X, Basu S, Dong X, Hoffmann M, et al. Genomics of Corynebacterium striatum, an emerging multidrug-resistant pathogen of immunocompromised patients. *Clin Microbiol Infect* 2018;24:1016.

36. Grimont PAD, Grimont F, Efstratiou A, De Zoysa A, Mazurova I, et al. International nomenclature for Corynebacterium diphtheriae ribotypes. *Res Microbiol* 2004;155:162-166.

37. Paveenkittiporn W, Sripakdee S, Koobkratok O, Sangkitporn S, Kerdsin A. Molecular epidemiology and antimicrobial susceptibility of outbreak-associated Corynebacterium diphtheriae in Thailand, 2012. *Infect Genet Evol* 2019;75:104007.

38. Chapartegui-González I, Fernández-Martínez M, Rodríguez-Fernández A, Rocha DJP, Aguiar E, et al. Antimicrobial susceptibility and characterization of resistance mechanisms of Corynebacterium urealyticum clinical isolates. *Antibiotics* (Basel) 2020;9:404.

39. Aldred KJ, Kerns RJ, Osheroff N. Mechanism of quinolone action and resistance. *Biochemistry* 2014;53:1565-1574.

40. Ahmad N, Hii SYF, Mohd Khalid MKN, Abd Wahab MA, Hashim R, et al. First draft genome sequences of Malaysian clinical isolates of Corynebacterium diphtheriae. *Genome Announc* 2017;5:e01670-16.

41. Dazas M, Badell E, Carmi-Leroy A, Criscuolo A, Brisse S. Taxonomic status of Corynebacterium diphtheriae biovar Belfanti and proposal of Corynebacterium belfantii sp. nov. *Int J Syst Evol Microbiol* 2018;68:3826-3831.

42. Grosse-Kock S, Kolodkina V, Schwalbe EC, Blom J, Burkovski A, et al. Genomic analysis of endemic clones of toxigenic and nontoxigenic Corynebacterium diphtheriae in Belarus during and after the major epidemic in 1990s. *BMC Genomics* 2017;18:873.

43. Hoefer A, Pampaka D, Herrera-León S, Peiró S, Varona S, et al. Molecular and epidemiological characterization of toxigenic and nontoxigenic Corynebacterium diphtheriae, Corynebacterium belfantii, Corynebacterium rouxii, and Corynebacterium ulcerans isolates identified in Spain from 2014 to 2019. *J Clin Microbiol* 2021;59:e02410-20.

44. Mahomed S, Archary M, Mutevedzi P, Mahabeer Y, Govender P, et al. An isolated outbreak of diphtheria in South Africa, 2015. *Epidemiol Infect* 2017;145:2100-2108.

45. Lodeiro-Colatosti A, Reischl U, Holzmann T, Hernández-Pereira CE, Ríquez A, et al. Diphtheria outbreak in Amerindian Communities, Wonken, Venezuela, 2016-2017. *Emerg Infect Dis* 2018;24:1340-1344.

46. Strauss RA, Herrera-Leon L, Guillén AC, Castro JS, Lorenz E, et al. Molecular and epidemiologic characterization of the diphtheria outbreak in Venezuela. *Sci Rep* 2021;11:6378.

47. Tedijanto C, Olesen SW, Grad YH, Lipsitch M. Estimating the proportion of bystander selection for antibiotic resistance among potentially pathogenic bacterial flora. *Proc Natl Acad Sci U S A* 2018;115:E1198-8-E1195.

48. Mina NV, Burdz T, Wiebe D, Rai JS, Rahim T, et al. Canada's first case of a multidrug-resistant Corynebacterium diphtheriae strain, isolated from a skin abscess. *J Clin Microbiol* 2011;49:4003-4005.

49. Lundborg CS, Tamhankar AJ. Antibiotic residues in the environment of South East Asia. *BMJ* 2017;358:j2440.

50. Sierra JM, Martinez-Martinez L, Vázquez F, Giralt E, Vila J. Relationship between mutations in the gyrA gene and quinolone resistance in clinical isolates of Corynebacterium striatum and Corynebacterium amycolatum. *Antimicrob Agents Chemother* 2005;49:1714-1719.

Five reasons to publish your next article with a Microbiology Society journal

1. When you submit to our journals, you are supporting Society activities for your community.

2. Experience a fair, transparent process and critical, constructive review.

3. If you are at a Publish and Read institution, you'll enjoy the benefits of Open Access across our journal portfolio.

4. Author feedback says our Editors are 'thorough and fair' and 'patient and caring'.

5. Increase your reach and impact and share your research more widely.

Find out more and submit your article at microbiologyresearch.org.