Volume 28, Number 11—November 2022
Sequence-Based Identification of Metronidazole-Resistant Clostridioides difficile Isolates
The plasmid pCD-METRO confers metronidazole resistance in Clostridioides difficile. We showed high sequence similarity among pCD-METRO plasmids from different isolates and identified pCD-METRO and associated metronidazole-resistant isolates in clinical and veterinary reservoirs in the Americas. We recommend using PCR or genomic assays to detect pCD-METRO in metronidazole-resistant C. difficile.
Clostridioides difficile is a major cause of antibiotic-associated colitis (1). Antimicrobial drug–resistant infections are a global economic and healthcare burden (2). Resistance is generally low to commonly prescribed antimicrobial drugs used for primary C. difficile infections. However, high rates of metronidazole resistance have been observed for C. difficile isolates carrying the 7-kb plasmid pCD-METRO, in particular for isolates belonging to PCR ribotype (RT) 010 and RT020 (clade 1) and the epidemic strain RT027 (clade 2) (3) (Figure, panel A). This plasmid has been reported in C. difficile isolates from countries in Europe.
Since the discovery of pCD-METRO, we have implemented PCR that uses primers oBH-1 (5′-CCTCGTAGAATCCGGTGCAA-3′) and oBH-2 (5′-TATTTCCTTGCCGCTGAGGT-3′) for national sentinel surveillance and diagnostics of C. difficile infections in the Netherlands. The primers are specific for open reading frame (ORF) 6 of pCD-METRO (Figure, panel A). Since 2019, we have tested 3,257 isolates and identified 8 (0.25%) additional pCD-METRO–positive isolates; this percentage is consistent with previous findings (3). We have a total of 27 human and animal C. difficile isolates in our collection that are pCD-METRO–positive. Most of the isolates (22/27, 81%) belong to nontoxigenic PCR RT010, including isolate 1143 from Brazil. Isolate 1143 is one of 8 canine isolates that showed phenotypic resistance to metronidazole (MIC = 32 mg/L) by Etest on Brucella blood agar (BBA); the Etest was performed at the Universidade Federal do Rio de Janeiro in Brazil. The isolate from Brazil confirmed that pCD-METRO is present in C. difficile not only in Europe but also in South America. The 1143 isolate was not characterized further because it belonged to PCR RT010, in which pCD-METRO is most frequently observed. The high number of C. difficile RT010 isolates carrying pCD-METRO might be related to genomic background of the isolates (4) or sampling bias; a higher prevalence of metronidazole resistance has been observed among RT010 strains (3,5). Low-frequency horizontal gene transfer is more likely to occur after prolonged co-colonization of nontoxigenic C. difficile and pCD-METRO donor bacteria, and acquisition of the plasmid might occur from a source after metronidazole exposure. For example, dogs carry nontoxigenic C. difficile frequently and are often treated with metronidazole (6). The presence of pCD-METRO in toxigenic isolates might also be underestimated; antimicrobial susceptibility testing is not routinely performed, and plasmid carriage is not assessed, even when metronidazole treatment fails.
Among C. difficile isolates from the Netherlands, we identified a toxigenic pCD-METRO–positive isolate (LUMCMM20 0132, National Center for Biotechnology Information [NCBI] BioSample no. SAMN26573026) from a symptomatic patient with C. difficile infection. The isolate belonged to RT005, a ribotype not reported previously to carry pCD-METRO. RT005 accounts for ≈4% of C. difficile isolates in Europe (7) and shows a similar prevalence in the Netherlands. The patient did not respond to metronidazole treatment, and a metronidazole Etest on BBA, performed at Leiden University Medical Center, confirmed the isolate was metronidazole-resistant (MIC = 8 mg/L). In contrast, a plasmid-negative RT005 isolate obtained earlier from the same patient (LUMCMM20 0131, NCBI BioSample no. SAMN26573027) was metronidazole-susceptible (MIC = 0.125 mg/L), further suggesting acquired resistance after pCD-METRO acquisition. Illumina whole-genome sequencing (NCBI BioProject accession no. PRJNA814863) and analysis of draft genomes using Kbase (8) indicated LUMCMM20 0131 and LUMCMM20 0132 were highly homologous, had an average nucleotide identity (ANI) of >99.99%, and were categorized as sequence type 6, clade 1 (9). We performed phylogenomic analysis by using IQ-TREE (10) and Roary (11) to show the 2 patient isolates were distinct from the RT005 reference strain LL005 (ANI 99.91–99.92) and RT012 reference strain 630 (ANI 99.16–99.18) (Figure, panel B). Moreover, we identified only 1 single-nucleotide polymorphism (SNP) when we aligned reads from LUMCMM20 0132 in a reference assembly against the draft LUMCMM20 0131 genome (minimum coverage 10, minimum variant frequency 0.8). We revealed that differences in the 2 patient isolates were driven by pCD-METRO carriage in LUMCMM20 0132 in a pangenome analysis using Kbase (8). We identified the pCD-METRO contig in the draft genome by using a homology search, removed terminal repeats, and circularized the sequences by using Geneious R9.1 (https://www.geneious.com). The resulting plasmid sequence was 100% identical to the pCD-METRO reference sequence (GenBank accession no. OM972905) (Figure, Panel A), which likely explains the metronidazole-resistant phenotype.
Because the presence of pCD-METRO is rare, we identified pCD-METRO–positive isolates in public repositories. We queried a curated database of >661,000 assembled bacterial genomes (12) by using a compact bit-sliced signature index with a k-mer similarity threshold of 0.4. A total of 465 assemblies were returned, but only 1 C. difficile isolate had a close-hit of 0.99 k-mer similarity. The other hits had k-mer similarities of <0.49 and included different species. The C. difficile isolate containing a contig with sequence homology to pCD-METRO was V356 (NCBI BioSample no. SAMN08813897). V356 is a nontoxigenic sequence type 15 isolate cultured from an intensive care unit patient in the United States who was an asymptomatic C. difficile carrier; the isolate clustered with other nontoxigenic C. difficile genomes (13). The isolate was metronidazole-resistant (MIC = 16–24 mg/L) in an Etest on BBA medium (the test was performed at Brigham and Women’s Hospital at the time of identification). We assembled the whole-genome sequence of the isolate by using Kbase (8) and reconstructed the pCD-METRO plasmid from the draft genome sequence as described above. The plasmid had 2 SNPs compared with the pCD-METRO reference sequence: G5441C, resulting in a Q96E amino acid substitution in the ORF7 hydrolase protein, and C5904A upstream of ORF7 (Figure, panel A); other variants are described elsewhere (3). V356 extends the geographic range of pCD-METRO and associated plasmid-mediated metronidazole resistance to North America.
To facilitate homology-based identifications, we deposited a pCD-METRO sequence assembly (GenBank accession no. OM972905) for inclusion in databases of antimicrobial resistance and mobilization determinants, such as the Comprehensive Antibiotic Resistance Database (14) and PlasmidFinder (15). The deposited file also indicates the sequence variants described in this study.
SNPs in pCD-METRO have been reported in ORF1, the ORF6–ORF7 intergenic region, ORF7, and ORF8, but not in the region that contains ORF2–6; major deletions or rearrangements in this plasmid have not been found. Thus, PCR-based approaches that detect conserved plasmid regions and genomic methods that examine pCD-METRO sequences can be used to identify pCD-METRO–containing C. difficile isolates. Of note, all isolates that carried pCD-METRO were confirmed to be metronidazole-resistant (MIC >2 mg/L) in susceptibility tests. Whereas the sequences responsible for metronidazole resistance in pCD-METRO have not yet been identified, we show that the presence of pCD-METRO in C. difficile predicts metronidazole resistance. We suggest using the invariant ORF2–6 region for PCR-based detection of pCD-METRO.
We found pCD-METRO in a metronidazole-resistant toxigenic RT005 isolate from the Netherlands and also identified pCD-METRO–associated metronidazole resistance in C. difficile isolates from North and South America. We recommend using sequence-based molecular approaches to detect pCD-METRO for plasmid-mediated metronidazole-resistant C. difficile.
Dr. Smits is an associate professor in medical microbiology at Leiden University Medical Center. His research interests focus on antimicrobial drugs, antimicrobial resistance, and plasmids of C. difficile.
We thank S. Nijssen for submitting fecal samples to the National Reference Laboratory for C. difficile, M. Delaney and K. Rainha for MIC testing, and members of the Leiden University Medical Center experimental bacteriology research group for useful discussions.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Genome sequencing and susceptibility testing at Brigham & Women’s Hospital was supported by the Hatch Family Foundation and US National Institutes of Health grant nos. R01 AI153605 and P30 DK034854.
The opinions expressed by the authors do not necessarily reflect the opinions of the institutions with which the authors are affiliated.
- Smits WK, Lyras D, Lacy DB, Wilcox MH, Kuijper EJ. Clostridium difficile infection. Nat Rev Dis Primers. 2016;2:16020.
- Murray CJL, Ikuta KS, Sharara F, Swetschinski L, Robles Aguilar G, Gray A, et al.; Antimicrobial Resistance Collaborators. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet. 2022;399:629–55.
- Boekhoud IM, Hornung BVH, Sevilla E, Harmanus C, Bos-Sanders IMJG, Terveer EM, et al. Plasmid-mediated metronidazole resistance in Clostridioides difficile. Nat Commun. 2020;11:598.
- Boekhoud IM, Sidorov I, Nooij S, Harmanus C, Bos-Sanders IMJG, Viprey V, et al.; COMBACTE-CDI Consortium. Haem is crucial for medium-dependent metronidazole resistance in clinical isolates of Clostridioides difficile. J Antimicrob Chemother. 2021;76:1731–40.
- Moura I, Spigaglia P, Barbanti F, Mastrantonio P. Analysis of metronidazole susceptibility in different Clostridium difficile PCR ribotypes. J Antimicrob Chemother. 2013;68:362–5.
- Albuquerque C, Pagnossin D, Landsgaard K, Simpson J, Brown D, Irvine J, et al. The duration of antibiotic treatment is associated with carriage of toxigenic and non-toxigenic strains of Clostridioides difficile in dogs. PLoS One. 2021;16:
- Freeman J, Vernon J, Pilling S, Morris K, Nicholson S, Shearman S, et al.; Pan-European Longitudinal Surveillance of Antibiotic Resistance among Prevalent Clostridium difficile Ribotypes Study Group. The ClosER study: results from a three-year pan-European longitudinal surveillance of antibiotic resistance among prevalent Clostridium difficile ribotypes, 2011-2014. Clin Microbiol Infect. 2018;24:724–31.
- Arkin AP, Cottingham RW, Henry CS, Harris NL, Stevens RL, Maslov S, et al. KBase: The United States Department of Energy Systems Biology Knowledgebase. Nat Biotechnol. 2018;36:566–9.
- Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018;3:124.
- Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.
- Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31:3691–3.
- Blackwell GA, Hunt M, Malone KM, Lima L, Horesh G, Alako BTF, et al. Exploring bacterial diversity via a curated and searchable snapshot of archived DNA sequences. PLoS Biol. 2021;19:
- Worley J, Delaney ML, Cummins CK, DuBois A, Klompas M, Bry L. Genomic determination of relative risks for Clostridioides difficile infection from asymptomatic carriage in intensive care unit patients. Clin Infect Dis. 2021;73:e1727–36.
- Alcock BP, Raphenya AR, Lau TTY, Tsang KK, Bouchard M, Edalatmand A, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020;48(D1):D517–25.
- Carattoli A, Hasman H. PlasmidFinder and In Silico pMLST: identification and typing of plasmid replicons in whole-genome sequencing (WGS). Methods Mol Biol. 2020;2075:285–94.
FigureCite This Article
Original Publication Date: October 13, 2022