Short-term genome stability of serial Clostridium difficile ribotype 027 isolates in an experimental gut model and recurrent human disease

Comparative Genomics of Clostridioides difficile

Clostridioides difficile, a Gram-positive spore-forming anaerobic bacterium, has rapidly emerged as the leading cause of nosocomial diarrhoea in hospitals. The availability of large numbers of genome sequences, mainly due to the use of next-generation sequencing methods, has undoubtedly shown their immense advantages in the determination of C. difficile population structure. The implementation of fine-scale comparative genomic approaches has paved the way for global transmission and recurrence studies, as well as more targeted studies, such as the PaLoc or CRISPR/Cas systems. In this chapter, we provide an overview of recent and significant findings on C. difficile using comparative genomic studies with implications for epidemiology, infection control and understanding of the evolution of C. difficile.

Transmission Patterns of Clostridioides difficile in a Non-Epidemic Setting Based on WGS Analysis

Clostridioides difficile is a major nosocomial pathogen and has a considerable burden on healthcare systems. Our objective was to determine the transmission patterns of C. difficile in a non-epidemic setting using whole-genome multi-locus sequence typing (wgMLST) and core-genome single-nucleotide polymorphism (cgSNP) analyses. A retrospective study was conducted in a 650-bed university hospital between January 2016 and February 2017. In total, 191 strains isolated from 169 symptomatic C. difficile infection (CDI) patients were analyzed by WGS. Sequences were compared using wgMLST and cgSNP analyses. Genetic data and ward movements were then combined to identify the transmission rate and the type of transmission. The transmission rate varied from 55/169 (19.5%) (wgMLST) to 33/169 (32.5%) (cgSNP). Most transmission was considered cryptic, irrespective of the genetic analysis (38/55 [69.1%] by wgMLST to 25/33 [75.8%] by cgSNP). No transmission within the same ward was observed. In a non-epidemic setting, most C. difficile transmission occurs from sources other than symptomatic CDI patients.

Plasmid-mediated metronidazole resistance in Clostridioides difficile

Metronidazole was until recently used as a first-line treatment for potentially life-threatening Clostridioides difficile (CD) infection. Although cases of metronidazole resistance have been documented, no clear mechanism for metronidazole resistance or a role for plasmids in antimicrobial resistance has been described for CD. Here, we report genome sequences of seven susceptible and sixteen resistant CD isolates from human and animal sources, including isolates from a patient with recurrent CD infection by a PCR ribotype (RT) 020 strain, which developed resistance to metronidazole over the course of treatment (minimal inhibitory concentration [MIC] = 8 mg L⁻¹). Metronidazole resistance correlates with the presence of a 7-kb plasmid, pCD-METRO. pCD-METRO is present in toxigenic and non-toxigenic resistant (n = 23), but not susceptible (n = 563), isolates from multiple countries. Introduction of a pCD-METRO-derived vector into a susceptible strain increases the MIC 25-fold. Our finding of plasmid-mediated resistance can impact diagnostics and treatment of CD infections.

Cases of C. difficile (CD) resistant to metronidazole have been reported but the mechanism remains enigmatic. Here the authors identify a plasmid, which correlates with metronidazole resistance status in a large international collection of CD isolates, and demonstrate that the plasmid can confer metronidazole resistance.

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.

Clostridium difficile trehalose metabolism variants are common and not associated with adverse patient outcomes when variably present in the same lineage

Background

Clostridium difficile ribotype-027, ribotype-078, and ribotype-017 are virulent and epidemic lineages. Trehalose metabolism variants in these ribotypes, combined with increased human trehalose consumption, have been hypothesised to have contributed to their emergence and virulence.

Methods

5232 previously whole-genome sequenced C. difficile isolates were analysed. Clinical isolates were used to investigate the impact of trehalose metabolism variants on mortality. Import data were used to estimate changes in dietary trehalose. Ribotype-027 virulence was investigated in a clinically reflective gut model.

Findings

Trehalose metabolism variants found in ribotype-027 and ribotype-017 were widely distributed throughout C. difficile clade-2 and clade-4 in 24/29 (83%) and 10/11 (91%) of sequence types (STs), respectively. The four-gene trehalose metabolism cluster described in ribotype-078 was common in genomes from all five clinically-important C. difficile clades (40/167 [24%] STs).

The four-gene cluster was variably present in 208 ribotype-015 infections (98 [47%]); 27/208 (13%) of these patients died within 30-days of diagnosis. Adjusting for age, sex, and infecting ST, there was no association between 30-day all-cause mortality and the four-gene cluster (OR 0.36 [95%CI 0.09–1.34, p = 0.13]).

Synthetic trehalose imports in the USA, UK, Germany and the EU were  < 1 g/capita/year during 2000–2006, and  < 9 g/capita/year 2007–2012, compared with dietary trehalose from natural sources of ~100 g/capita/year.

Trehalose supplementation did not increase ribotype-027 virulence in a clinically-validated gut model.

Interpretation

Trehalose metabolism variants are common in C. difficile. Increases in total dietary trehalose during the early-mid 2000s C. difficile epidemic were likely relatively minimal. Alternative explanations are required to explain why ribotype-027, ribotype-078 and ribotype-017 have been successful.

Funding

National Institute for Health Research. Gut model experiments only: Hayashibara Co. Ltd.

Increased sporulation underpins adaptation of Clostridium difficile strain 630 to a biologically-relevant faecal environment, with implications for pathogenicity

Clostridium difficile virulence is driven primarily by the processes of toxinogenesis and sporulation, however many in vitro experimental systems for studying C. difficile physiology have arguably limited relevance to the human colonic environment. We therefore created a more physiologically–relevant model of the colonic milieu to study gut pathogen biology, incorporating human faecal water (FW) into growth media and assessing the physiological effects of this on C. difficile strain 630. We identified a novel set of C. difficile–derived metabolites in culture supernatants, including hexanoyl– and pentanoyl–amino acid derivatives by LC-MSⁿ. Growth of C. difficile strain 630 in FW media resulted in increased cell length without altering growth rate and RNA sequencing identified 889 transcripts as differentially expressed (p < 0.001). Significantly, up to 300–fold increases in the expression of sporulation–associated genes were observed in FW media–grown cells, along with reductions in motility and toxin genes’ expression. Moreover, the expression of classical stress–response genes did not change, showing that C. difficile is well–adapted to this faecal milieu. Using our novel approach we have shown that interaction with FW causes fundamental changes in C. difficile biology that will lead to increased disease transmissibility.

Comparison of Control of Clostridium difficile Infection in Six English Hospitals Using Whole-Genome Sequencing

Background.

Variation in Clostridium difficile infection (CDI) rates between healthcare institutions suggests overall incidence could be reduced if the lowest rates could be achieved more widely.

Methods.

We used whole-genome sequencing (WGS) of consecutive C. difficile isolates from 6 English hospitals over 1 year (2013–14) to compare infection control performance. Fecal samples with a positive initial screen for C. difficile were sequenced. Within each hospital, we estimated the proportion of cases plausibly acquired from previous cases.

Results.

Overall, 851/971 (87.6%) sequenced samples contained toxin genes, and 451 (46.4%) were fecal-toxin-positive. Of 652 potentially toxigenic isolates >90-days after the study started, 128 (20%, 95% confidence interval [CI] 17–23%) were genetically linked (within ≤2 single nucleotide polymorphisms) to a prior patient’s isolate from the previous 90 days. Hospital 2 had the fewest linked isolates, 7/105 (7%, 3–13%), hospital 1, 9/70 (13%, 6–23%), and hospitals 3–6 had similar proportions of linked isolates (22–26%) (P ≤ .002 comparing hospital-2 vs 3–6). Results were similar adjusting for locally circulating ribotypes. Adjusting for hospital, ribotype-027 had the highest proportion of linked isolates (57%, 95% CI 29–81%). Fecal-toxin-positive and toxin-negative patients were similarly likely to be a potential transmission donor, OR = 1.01 (0.68–1.49). There was no association between the estimated proportion of linked cases and testing rates.

Conclusions.

WGS can be used as a novel surveillance tool to identify varying rates of C. difficile transmission between institutions and therefore to allow targeted efforts to reduce CDI incidence.

Comparative Genomics of Clostridium difficile

Clostridium difficile, a gram-positive spore-forming anaerobic bacterium, has rapidly emerged as the leading cause of nosocomial diarrhoea in hospitals. The availability of genome sequences in large numbers, mainly due to the use of next-generation sequencing methods, have undoubtedly shown their immense advantages in the determination of the C. difficile population structure. The implementation of fine-scale comparative genomic approaches have paved the way to global transmission and recurrence studies, but also more targeted studies such as the PaLoc or the CRISPR/Cas systems. In this chapter, we provide an overview of the recent and significant findings on C. difficile using comparative genomics studies with implication for the epidemiology, infection control and understanding of the evolution of C. difficile.

CRISPR Diversity and Microevolution in Clostridium difficile

Virulent strains of Clostridium difficile have become a global health problem associated with morbidity and mortality. Traditional typing methods do not provide ideal resolution to track outbreak strains, ascertain genetic diversity between isolates, or monitor the phylogeny of this species on a global basis. Here, we investigate the occurrence and diversity of clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated genes (cas) in C. difficile to assess the potential of CRISPR-based phylogeny and high-resolution genotyping. A single Type-IB CRISPR-Cas system was identified in 217 analyzed genomes with cas gene clusters present at conserved chromosomal locations, suggesting vertical evolution of the system, assessing a total of 1,865 CRISPR arrays. The CRISPR arrays, markedly enriched (8.5 arrays/genome) compared with other species, occur both at conserved and variable locations across strains, and thus provide a basis for typing based on locus occurrence and spacer polymorphism. Clustering of strains by array composition correlated with sequence type (ST) analysis. Spacer content and polymorphism within conserved CRISPR arrays revealed phylogenetic relationship across clades and within ST. Spacer polymorphisms of conserved arrays were instrumental for differentiating closely related strains, e.g., ST1/RT027/B1 strains and pathogenicity locus encoding ST3/RT001 strains. CRISPR spacers showed sequence similarity to phage sequences, which is consistent with the native role of CRISPR-Cas as adaptive immune systems in bacteria. Overall, CRISPR-Cas sequences constitute a valuable basis for genotyping of C. difficile isolates, provide insights into the micro-evolutionary events that occur between closely related strains, and reflect the evolutionary trajectory of these genomes.

Emergence and spread of predominantly community-onset Clostridium difficile PCR ribotype 244 infection in Australia, 2010 to 2012

We describe an Australia-wide Clostridium difficile outbreak in 2011 and 2012 involving the previously uncommon ribotype 244. In Western Australia, 14 of 25 cases were community-associated, 11 were detected in patients younger than 65 years, 14 presented to emergency/outpatient departments, and 14 to non-tertiary/community hospitals. Using whole genome sequencing, we confirm ribotype 244 is from the same C. difficile clade as the epidemic ribotype 027. Like ribotype 027, it produces toxins A, B, and binary toxin, however it is fluoroquinolone-susceptible and thousands of single nucleotide variants distinct from ribotype 027. Fifteen outbreak isolates from across Australia were sequenced. Despite their geographic separation, all were genetically highly related without evidence of geographic clustering, consistent with a point source, for example affecting the national food chain. Comparison with reference laboratory strains revealed the outbreak clone shared a common ancestor with isolates from the United States and United Kingdom (UK). A strain obtained in the UK was phylogenetically related to our outbreak. Follow-up of that case revealed the patient had recently returned from Australia. Our data demonstrate new C. difficile strains are an on-going threat, with potential for rapid spread. Active surveillance is needed to identify and control emerging lineages.

Clostridium difficile surveillance: harnessing new technologies to control transmission

Clostridium difficile surveillance allows outbreaks of cases clustered in time and space to be identified and further transmission prevented. Traditionally, manual detection of groups of cases diagnosed in the same ward or hospital, often followed by retrospective reference laboratory genotyping, has been used to identify outbreaks. However, integrated healthcare databases offer the prospect of automated real-time outbreak detection based on statistically robust methods, and accounting for contacts between cases, including those distant to the ward of diagnosis. Complementary to this, rapid benchtop whole genome sequencing, and other highly discriminatory genotyping, has the potential to distinguish which cases are part of an outbreak with high precision and in clinically relevant timescales. These new technologies are likely to shape future surveillance.

Whole-genome sequencing demonstrates that fidaxomicin is superior to vancomycin for preventing reinfection and relapse of infection with Clostridium difficile

Whole-genome sequencing was used to determine whether the reductions in recurrence of Clostridium difficile infection observed with fidaxomicin in pivotal phase 3 trials occurred by preventing relapse of the same infection, by preventing reinfection with a new strain, or by preventing both outcomes. Paired isolates of C. difficile were available from 93 of 199 participants with recurrences (28 were treated with fidaxomicin, and 65 were treated with vancomycin). Given C. difficile evolutionary rates, paired samples ≤2 single-nucleotide variants (SNVs) apart were considered relapses, paired samples >10 SNVs apart were considered reinfection, and those 3–10 SNVs apart (or without whole-genome sequences) were considered indeterminate in a competing risks survival analysis. Fidaxomicin reduced the risk of both relapse (competing risks hazard ratio [HR], 0.40 [95% confidence interval {CI}, .25–.66]; P = .0003) and reinfection (competing risks HR, 0.33 [95% CI, 0.11–1.01]; P = .05).

Comparison of multilocus variable-number tandem-repeat analysis and whole-genome sequencing for investigation of Clostridium difficile transmission

No study to date has compared multilocus variable-number tandem-repeat analysis (MLVA) and whole-genome sequencing (WGS) in an investigation of the transmission of Clostridium difficile infection. Isolates from 61 adults with ongoing and/or recurrent C. difficile infections and 17 asymptomatic carriage episodes in children (201 samples), as well as from 61 suspected outbreaks affecting 2 to 41 patients in 31 hospitals in the United Kingdom (300 samples), underwent 7-locus MLVA and WGS in parallel. When the first and last samples from the same individual taken for a median (interquartile range [IQR]) of 63 days (43 to 105 days) apart were compared, the estimated rates of the evolution of single nucleotide variants (SNVs), summed tandem-repeat differences (STRDs), and locus variants (LVs) were 0.79 (95% confidence interval [CI], 0.00 to 1.75), 1.63 (95% CI, 0.00 to 3.59), and 1.21 (95% CI, 0.00 to 2.67)/called genome/year, respectively. Differences of >2 SNVs and >10 STRDs have been used to exclude direct case-to-case transmission. With the first serial sample per individual being used to assess discriminatory power, across all pairs of samples sharing a PCR ribotype, 192/283 (68%) differed by >10 STRDs and 217/283 (77%) by >2 SNVs. Among all pairs of cases from the same suspected outbreak, 1,190/1,488 (80%) pairs had concordant results using >2 SNVs and >10 STRDs to exclude transmission. For the discordant pairs, 229 (15%) had ≥2 SNVs but ≤10 STRDs, and 69 (5%) had ≤2 SNVs but ≥10 STRDs. Discordant pairs had higher numbers of LVs than concordant pairs, supporting the more diverse measure in each type of discordant pair. Conclusions on whether the potential outbreaks were confirmed were concordant in 58/61 (95%) investigations. Overall findings using MLVA and WGS were very similar despite the fact that they analyzed different parts of the bacterial genome. With improvements in WGS technology, it is likely that MLVA locus data will be available from WGS in the near future.