A global to local genomics analysis of Clostridioides difficile ST1/RT027 identifies cryptic transmission events in a northern Arizona healthcare network

Whole-genome sequencing of toxigenic Clostridioides difficile reveals multidrug resistance and virulence genes in strains of environmental and animal origin

BACKGROUND

Clostridioides difficile has been recognized as an emerging pathogen in both humans and animals. In this context, antimicrobial resistance plays a major role in driving the spread of this disease, often leading to therapeutic failure. Moreover, recent increases in community-acquired C. difficile infections have led to greater numbers of investigations into the animal origin of the disease. The aim of this study was to evaluate the genetic similarities between 23 environmental and animal isolates by using whole-genome sequencing and to determine antimicrobial resistance and virulence factor genes in toxigenic C. difficile strains to provide important data for the development of diagnostic methods or treatment guidelines.

RESULTS

The most common sequence type was ST11 (87%), followed by ST2 (9%) and ST19 (4%). In addition, 86.95% of the strains exhibited multidrug resistance, with antimicrobial resistance to mainly aminoglycosides, fluoroquinolones, tetracycline and B-lactams; nevertheless, one strain also carried other resistance genes that conferred resistance to lincosamide, macrolides, streptogramin a, streptogramin b, pleuromutilin, oxazolidinone and amphenicol. In addition, a wide range of virulence factor genes, such as those encoding adherence factors, exoenzymes and toxins, were found. However, we observed variations between toxinotypes, ribotypes and sequence types.

CONCLUSIONS

The results of this study demonstrated significant genetic similarity between ST11 strains isolated from environmental sampling and from animal origin; these strains may represent a reservoir for community-acquired C. difficile infection, which is becoming a growing public health threat due to the development of multridug resistant (MDR) bacteria and the number of virulence factors detected.

A local-scale One Health genomic surveillance of Clostridioides difficile demonstrates highly related strains from humans, canines, and the environment

Although infections caused by Clostridioides difficile have historically been attributed to hospital acquisition, growing evidence supports the role of community acquisition in C. difficile infection (CDI). Symptoms of CDI can range from mild, self-resolving diarrhoea to toxic megacolon, pseudomembranous colitis, and death. In this study, we sampled C. difficile from clinical, environmental, and canine reservoirs in Flagstaff, Arizona, USA, to understand the distribution and transmission of the pathogen in a One Health framework; Flagstaff is a medium-sized, geographically isolated city with a single hospital system, making it an ideal site to characterize genomic overlap between sequenced C. difficile isolates across reservoirs. An analysis of 562 genomes from Flagstaff isolates identified 65 sequence types (STs), with eight STs being found across all three reservoirs and another nine found across two reservoirs. A screen of toxin genes in the pathogenicity locus identified nine STs where all isolates lost the toxin genes needed for CDI manifestation (tcdB, tcdA), demonstrating the widespread distribution of non-toxigenic C. difficile (NTCD) isolates in all three reservoirs; 15 NTCD genomes were sequenced from symptomatic, clinical samples, including two from mixed infections that contained both tcdB+ and tcdB- isolates. A comparative single nucleotide polymorphism (SNP) analysis of clinically derived isolates identified 78 genomes falling within clusters separated by ≤2 SNPs, indicating that ~19 % of clinical isolates are associated with potential healthcare-associated transmission clusters; only symptomatic cases were sampled in this study, and we did not sample asymptomatic transmission. Using this same SNP threshold, we identified genomic overlap between canine and soil isolates, as well as putative transmission between environmental and human reservoirs. The core genome of isolates sequenced in this study plus a representative set of public C. difficile genomes (n=136), was 2690 coding region sequences, which constitutes ~70 % of an individual C. difficile genome; this number is significantly higher than has been published in some other studies, suggesting that genome data quality is important in understanding the minimal number of genes needed by C. difficile. This study demonstrates the close genomic overlap among isolates sampled across reservoirs, which was facilitated by maximizing the genomic search space used for comprehensive identification of potential transmission events. Understanding the distribution of toxigenic and non-toxigenic C. difficile across reservoirs has implications for surveillance sampling strategies, characterizing routes of infections, and implementing mitigation measures to limit human infection.

K-mer based prediction of Clostridioides difficile relatedness and ribotypes

Comparative analysis of Clostridioides difficile whole-genome sequencing (WGS) data enables fine scaled investigation of transmission and is increasingly becoming part of routine surveillance. However, these analyses are constrained by the computational requirements of the large volumes of data involved. By decomposing WGS reads or assemblies into k-mers and using the dimensionality reduction technique MinHash, it is possible to rapidly approximate genomic distances without alignment. Here we assessed the performance of MinHash, as implemented by sourmash, in predicting single nucleotide differences between genomes (SNPs) and C. difficile ribotypes (RTs). For a set of 1905 diverse C. difficile genomes (differing by 0–168 519 SNPs), using sourmash to screen for closely related genomes, at a sensitivity of 100 % for pairs ≤10 SNPs, sourmash reduced the number of pairs from 1 813 560 overall to 161 934, i.e. by 91 %, with a positive predictive value of 32 % to correctly identify pairs ≤10 SNPs (maximum SNP distance 4144). At a sensitivity of 95 %, pairs were reduced by 94 % to 108 266 and PPV increased to 45 % (maximum SNP distance 1009). Increasing the MinHash sketch size above 2000 produced minimal performance improvement. We also explored a MinHash similarity-based ribotype prediction method. Genomes with known ribotypes (n=3937) were split into a training set (2937) and test set (1000) randomly. The training set was used to construct a sourmash index against which genomes from the test set were compared. If the closest five genomes in the index had the same ribotype this was taken to predict the searched genome’s ribotype. Using our MinHash ribotype index, predicted ribotypes were correct in 780/1000 (78 %) genomes, incorrect in 20 (2 %), and indeterminant in 200 (20 %). Relaxing the classifier to 4/5 closest matches with the same RT improved the correct predictions to 87 %. Using MinHash it is possible to subsample C. difficile genome k-mer hashes and use them to approximate small genomic differences within minutes, significantly reducing the search space for further analysis.

Clostridioides difficile phage biology and application

Clostridium difficile, now reclassified as Clostridioides difficile, is the causative agent of C. difficile infections (CDI). CDI is particularly challenging in healthcare settings because highly resistant spores of the bacterium can persist in the environment, making it difficult to curb outbreaks. Dysbiosis of the microbiota caused by the use of antibiotics is the primary factor that allows C. difficile to colonize the gut and cause diarrhea and colitis. For this reason, antibiotics targeting C. difficile can be ineffective at preventing recurrent episodes because they exacerbate and prolong dysbiosis. The emergence of antibiotic resistance in C. difficile also presents a significant threat. The diverse array of bacteriophages (phages) that infect C. difficile could offer new treatment strategies and greater insight into the biology of the pathogen. In this review, we summarize the current knowledge regarding C. difficile phages and discuss what is understood about their lifestyles and genomics. Then, we examine how phage infection modifies bacterial gene expression and pathogenicity. Finally, we discuss the potential clinical applications of C. difficile phages such as whole phage therapy and phage-derived products, and we highlight the most promising strategies for further development.

Development and Implementation of Whole Genome Sequencing-Based Typing Schemes for Clostridioides difficile

Clostridioides difficile is an important nosocomial pathogen increasingly observed in the community and in different non-human reservoirs. The epidemiology and transmissibility of C. difficile has been studied using a variety of typing methods, including more recently developed whole-genome sequence (WGS) analysis that is becoming used routinely for bacterial typing worldwide. Here we review the schemes for WGS-based typing methods available for C. difficile and their applications in the field of human C. difficile infection (CDI). The two main approaches to discover genomic variations are single nucleotide variant (SNV) analysis and methods based on gene-by-gene comparisons (frequently called core genome or whole genome MLST, cgMLST, or wgMLST). SNV analysis currently provides the ultimate resolution, however, typing nomenclature and standardized methodology are missing. On the other hand, gene-by-gene approaches allow portability and standardized nomenclature, and are therefore becoming increasingly popular in bacterial epidemiology and outbreak investigation. Two commercial software packages (BioNumerics and Ridom SeqSphere+) and an open source database (EnteroBase) for allele and sequence type determination for C. difficile are currently available. Proof-of-concept WGS studies have already enabled advances in the investigation of the population structure of C. difficile species, microevolution within the epidemic strains, intercontinental transmission over time and in tracking of transmission events. WGS of clinical C. difficile isolates demonstrated a considerable genetic diversity suggesting diverse reservoirs for CDI. WGS was also shown to aid in resolving relapses and reinfections in recurrent CDI and has potential for use as a tool for assessing hospital infection prevention and control performance.

Integrated genomic epidemiology and phenotypic profiling of Clostridium difficile across intra-hospital and community populations in Colombia

Clostridium difficile, the causal agent of antibiotic-associated diarrhea, has a complex epidemiology poorly studied in Latin America. We performed a robust genomic and phenotypic profiling of 53 C. difficile clinical isolates established from diarrheal samples from either intrahospital (IH) or community (CO) populations in central Colombia. In vitro tests were conducted to evaluate the cytopathic effect, the minimum inhibitory concentration of ten antimicrobial agents, the sporulation efficiency and the colony forming ability. Eleven different sequence types (STs) were found, the majority present individually in each sample, however in three samples two different STs were isolated. Interestingly, CO patients were infected with STs associated with hypervirulent strains (ST-1 in Clade-2). Three coexistence events (two STs simultaneously detected in the same sample) were observed always involving ST-8 from Clade-1. A total of 2,502 genes were present in 99% of the isolates with 95% of identity or more, it represents a core genome of 28.6% of the 8,735 total genes identified in the set of genomes. A high cytopathic effect was observed for the isolates positive for the two main toxins but negative for binary toxin (TcdA+/TcdB+/CDT- toxin production type), found only in Clade-1. Molecular markers conferring resistance to fluoroquinolones (cdeA and gyrA) and to sulfonamides (folP) were the most frequent in the analyzed genomes. In addition, 15 other markers were found mostly in Clade-2 isolates. These results highlight the regional differences that C. difficile isolates display, being in this case the CO isolates the ones having a greater number of accessory genes and virulence-associated factors.