Clostridioides difficile is a pathogen often associated with hospital-acquired infection or antimicrobial-induced disease; however, increasing evidence indicates infections can result from community or environmental sources. Most genomic sequencing of C. difficile has focused on clinical strains, although evidence is growing that C. difficile spores are widespread in soil and water in the environment. In this study, we sequenced 38 genomes collected from soil and water isolates in Flagstaff (AZ, USA) and Slovenia in an effort targeted towards environmental surveillance of C. difficile . At the average nucleotide identity (ANI) level, the genomes were divergent to C. difficile at a threshold consistent with different species. A phylogenetic analysis of these divergent genomes together with Clostridioides genomes available in public repositories confirmed the presence of three previously described, cryptic Clostridioide s species and added two additional clades. One of the cryptic species (C-III) was almost entirely composed of Arizona and Slovenia genomes, and contained distinct sub-groups from each region (evidenced by SNP and gene-content differences). A comparative genomics analysis identified multiple unique coding sequences per clade, which can serve as markers for subsequent environmental surveys of these cryptic species. Homologues to the C. difficile toxin genes, tcdA and tcdB, were found in cryptic species genomes, although they were not part of the typical pathogenicity locus observed in C. difficile , and in silico PCR suggested that some would not amplify with widely used PCR diagnostic tests. We also identified gene homologues in the binary toxin cluster, including some present on phage and, for what is believed to be the first time, on a plasmid. All isolates were obtained from environmental samples, so the function and disease potential of these toxin homologues is currently unknown. Enzymatic profiles of a subset of cryptic isolates (n=5) demonstrated differences, suggesting that these isolates contain substantial metabolic diversity. Antimicrobial resistance (AMR) was observed across a subset of isolates (n=4), suggesting that AMR mechanisms are intrinsic to the genus, perhaps originating from a shared environmental origin. This study greatly expands our understanding of the genomic diversity of Clostridioides . These results have implications for C. difficile One Health research, for more sensitive C. difficile diagnostics, as well as for understanding the evolutionary history of C. difficile and the development of pathogenesis.

No abstract

All supporting data, code and protocols have been provided within the article or through supplementary data files. Two supplementary figures and eighteen supplementary tables are available with the online version of this article and deposited in figshare:

A single-chromosome closed genome of Peptacetobacter (Clostridium) hiranonis strain DGF055142 was generated using Illumina MiSeq short reads paired with Oxford Nanopore MinION long reads. This isolate was obtained from a canine in Flagstaff, Arizona, in 2019. Peptacetobacter (C.) hiranonis was hypothesized to contribute to canine Clostridium difficile infection resistance.

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.