•Multi-omics reveals markers of CDI in pediatric IBD patients
•Identification of metabolites reveals distinctive features for IBD and CDI
•Isocaproyltaurine is made by C. difficile and associates with active IBD
•Identifies biomarkers potentially useful for distinguishing disease processes
Children with inflammatory bowel diseases (IBD) are particularly vulnerable to infection with Clostridioides difficile (CDI). IBD and IBD + CDI have overlapping symptoms but respond to distinctive treatments, highlighting the need for diagnostic biomarkers. Here, we studied pediatric patients with IBD and IBD + CDI, comparing longitudinal data on the gut microbiome, metabolome, and other measures. The microbiome is dysbiotic and heterogeneous in both disease states, but the metabolome reveals disease-specific patterns. The IBD group shows increased concentrations of markers of inflammation and tissue damage compared with healthy controls, and metabolic changes associate with susceptibility to CDI. In IBD + CDI, we detect both metabolites associated with inflammation/tissue damage and fermentation products produced by C. difficile. The most discriminating metabolite found is isocaproyltaurine, a covalent conjugate of a distinctive C. difficile fermentation product (isocaproate) and an amino acid associated with tissue damage (taurine), which may be useful as a joint marker of the two disease processes.
To be redirected to the article in its entirety, please click on the following link:
World Patient Safety Day, is an initiative by the World Health Organisation to raise awareness of the need to continuously improve patient safety, building on the foundations of a safer culture and safer systems.
World Patient Safety Day 2020 will bring together health workers, patients, families, caregivers, communities, health care leaders and policy-makers, all showing commitment to health worker safety and patient safety. Health workers include all people engaged in actions whose primary intent is to enhance health. The year 2020 has been designated as the International Year of the Nurse and Midwife
‘Safe health workers, Safe patients’
Members of the C Diff Foundation recognize the importance of this annual WHO “World Patient Safety Day” campaign and we ask fellow patients, health workers, health care leaders, patients, and healthcare professionals to actively participate and share this information with others.
Share using any of the following hashtags: #PatientSafety#WorldPatientSafetyDay#HealthWorker
Thursday, September 17th is World Patient Safety Day. No one should be harmed in health care. And yet, every day, thousands of patients suffer avoidable harm while receiving care. The ongoing COVID-19 pandemic has exerted unprecedented pressure on health systems and health workers. An alarming number of health workers have also been harmed while responding to the COVID-19 pandemic. Safety of health workers has a significant impact on the safety of patients In view of this, the theme of World Patient Safety Day this year is Health Worker Safety: A Priority for Patient Safety. If you or a loved one are a patient, be actively involved in your own care. Enhance your knowledge about safety in health care and raise safety concerns with your health providers. Ask questions to be aware of your health condition and treatment. Speak up for the safety of your care as well as the safety of health workers!
Sprinkle one side with salt, pepper, as desired. Spices can be omitted or replaced with spices of choice.
Place the breast in the air fryer basket seasoned side down. Then season the other side.
Cook at 360° for 9 minutes (for 8 oz chicken breast). Then flip the chicken breast over and cook for another 9 minutes. Larger chicken breasts need more cook time, and smaller ones need less time. Always use a meat thermometer to check the internal temperature.
Open the air fryer immediately so it doesn’t continue cooking in the heat.
* Food Safety: Check the internal temperature by placing the instant-read meat thermometer in the thickest part of each piece of chicken. The temperature should read 158°-160° F. If it is below that, close the lid of the air fryer to let the chicken cook in the residual heat for a few minutes (2 to 4 additional minutes). Then check it again with the meat thermometer.
Place chicken pieces on a plate and loosely cover with foil and let it rest for 5 minutes. It will continue to cook in the residual heat (and the juices will go back into the meat). Then after 5 minutes take the temperature again. The internal temperature should not be below 165° F.
Serve immediately, or prepare left-over pieces of chicken to be placed in the refrigerator immediately to be served within two days. Left-over chicken pieces can also be prepared for other recipes (e.g. chicken soup, chicken salad, chicken and pasta).
Additional Cooking Notes:
Air Fryer Chicken Breast Cook Temp & Times Turn chicken pieces over halfway through cook time (Your model and size of the air fryer may vary from this slightly)
Please check your Air Fryer Manual for cooking time and instructions.
Average Cooking Temperature and Times: 6 oz chicken breast 360° for 15 minutes 8 oz chicken breast 360° for 18 minutes 10 oz chicken breast 360° for 20 minutes
When cooking multiple chicken breasts at once, use the average weight of the chicken pieces added together.
Martinique Frentrup1, Zhemin Zhou2, Matthias Steglich1,3, Jan P. Meier-Kolthoff1, Markus Göker1, Thomas Riedel1,3, Boyke Bunk1, Cathrin Spröer1, Jörg Overmann1,3,4, Marion Blaschitz5, Alexander Indra5, Lutz von Müller6, Thomas A. Kohl7,8, Stefan Niemann7,8, Christian Seyboldt9, Frank Klawonn10,11, Nitin Kumar12, Trevor D. Lawley12, Sergio García-Fernández13,14, Rafael Cantón13,14, Rosa del Campo13,14, Ortrud Zimmermann15, Uwe Groß15, Mark Achtman2, Ulrich Nübel1,3,4
Clostridioides difficile is the primary infectious cause of antibiotic-associated diarrhea. Local transmissions and international outbreaks of this pathogen have been previously elucidated by bacterial whole-genome sequencing, but comparative genomic analyses at the global scale were hampered by the lack of specific bioinformatic tools. Here we introduce a publicly accessible database within EnteroBase (http://enterobase.warwick.ac.uk) that automatically retrieves and assembles
C. difficile short-reads from the public domain, and calls alleles for core-genome multilocus sequence typing (cgMLST). We demonstrate that comparable levels of resolution and precision are attained by EnteroBase cgMLST and single-nucleotide polymorphism analysis. EnteroBase currently contains 18 254 quality-controlled
C. difficile genomes, which have been assigned to hierarchical sets of single-linkage clusters by cgMLST distances. This hierarchical clustering is used to identify and name populations of C. difficile at all epidemiological levels, from recent transmission chains through to epidemic and endemic strains. Moreover, it puts newly collected isolates into phylogenetic and epidemiological context by identifying related strains among all previously published genome data. For example, HC2 clusters (i.e. chains of genomes with pairwise distances of up to two cgMLST alleles) were statistically associated with specific hospitals (P<10−4) or single wards (P=0.01) within hospitals, indicating they represented local transmission clusters. We also detected several HC2 clusters spanning more than one hospital that by retrospective epidemiological analysis were confirmed to be associated with inter-hospital patient transfers. In contrast, clustering at level HC150 correlated with k-mer-based classification and was largely compatible with PCR ribotyping, thus enabling comparisons to earlier surveillance data. EnteroBase enables contextual interpretation of a growing collection of assembled, quality-controlled C. difficile genome sequences and their associated metadata. Hierarchical clustering rapidly identifies database entries that are related at multiple levels of genetic distance, facilitating communication among researchers, clinicians, and public-health officials who are combatting disease caused by C. difficile.
Clostridioides (Clostridium) difficile,
genomic population structure,
Author NotesAll supporting data, code and protocols have been provided within the article or through supplementary data files. Three supplementary tables and eight supplementary figures are available with the online version of this article.Sequence accession numbers: PRJEB33768, PRJEB33779, PRJEB33780.†These authors contributed equally to this work‡These authors also contributed equally to this work
Clostridioides difficile is a major cause of healthcare-associated diarrhea and causes large infection outbreaks. Whole-genome sequencing is increasingly applied for genotyping C. difficile, with the objectives to monitor and curb the pathogen’s spread. We present a publicly accessible database for quality-controlled genome sequences from C. difficile that enables contextual interpretation of newly collected isolates by identifying related strains among published data. It also provides a nomenclature for genomic types to facilitate communication about transmission chains, epidemics, and phylogenetic lineages. Finally, we demonstrate that genome-based hierarchical clustering is largely compatible with previously used molecular typing techniques, thus enabling comparisons to earlier surveillance data.
All genome sequencing data were submitted to the European Nucleotide Archive (www.ebi.ac.uk/ena) under study numbers PRJEB33768, PRJEB33779 and PRJEB33780. The Clostridioides database within EnteroBase is publicly accessible at http://enterobase.warwick.ac.uk. In addition, stand-alone versions of all EnteroBase tools are available at https://github.com/zheminzhou/EToKi.
The anaerobic gut bacterium Clostridioides difficile (formerly Clostridium difficile)  is the primary cause of antibiotic-associated diarrhea in Europe and North America . Molecular genotyping of C. difficile isolates has demonstrated international dissemination of diverse strains through healthcare systems [3–5], the community  and livestock production facilities [7, 8]. Previously, genotyping was commonly performed by PCR ribotyping or DNA macrorestriction. More recent publications have documented that genome-wide single-nucleotide polymorphisms (SNPs) from whole-genome sequences provide improved discrimination, and such analyses have enabled dramatic progress in our understanding of the emergence and spread of epidemic strains [9–12] and the epidemiology of local transmission [13, 14]. Eyre and colleagues have argued that transmission of C. difficile isolates within a hospital environment can be recognized with high probability as chains of genomes, which differ by up to two SNPs whereas genomes, which differ by at least ten genomic SNPs represent unrelated bacteria [13, 15]. However, SNP analyses require sophisticated bioinformatic tools and are difficult to standardize [16, 17]. A convenient alternative to SNP-based genotyping is offered by the commercial software SeqSphere, which implements a core-genome multilocus sequence typing scheme (cgMLST) for the analysis of genomic diversity in
C. difficile  and other organisms. Indeed, cgMLST  confirmed the prior conclusion from genomic SNP analyses  that a common clone of C. difficile had been isolated over two successive years at a hospital in China . However, a recent quantitative comparison of the two methods showed that SeqSphere’s cgMLST achieved a low predictive value (41 %) for identifying isolate pairs that were closely related by the ≤2 SNPs’ criterion . cgMLST of genomic sequences of a variety of bacterial pathogens can also be performed with EnteroBase (http://enterobase.warwick.ac.uk/), which has been developed over the last few years with the goal of facilitating genomic analyses by microbiologists . EnteroBase automatically retrieves Illumina short-read sequences from public short-read archives. It uses a consistent assembly pipeline to automatically assemble these short-reads into draft genomes consisting of multiple contigs, and presents the assembled genomes together with their metadata for public access . It also performs the same procedures on sequencing data uploaded by its registered users. Assembled genomes that pass quality control are genotyped by MLST at the levels of seven-gene MLST, ribosomal MLST (rMLST), cgMLST and whole-genome MLST (wgMLST) [21, 22]. EnteroBase supports subsequent analyses based on either SNPs or cgMLST alleles using the GrapeTree or Dendrogram visualization tools . EnteroBase also assigns these genotypes to populations by hierarchical clustering (HierCC), which supports the identification of close relatives at the global level . Originally, EnteroBase was restricted to the bacterial genera Salmonella, Escherichia, Yersinia and Moraxella but since January 2018, EnteroBase has included a database for genomes and their metadata for the genus Clostridioides. In June 2020, EnteroBase contained 18 254 draft genomes of C. difficile plus one genome of C. mangenotii. These included over 900 unpublished draft genomes that were sequenced at the Leibniz Institute DSMZ, as well as 80 complete genome sequences based on Pacific Biosciences plus Illumina sequencing technologies. It also included 862 unpublished draft genomes that were sequenced at the Wellcome Sanger Institute.
Here we show that comparable levels of resolution and precision are attained by EnteroBase cgMLST as by SNP analyses. We also summarize the genomic diversity that accumulated during recurring infections within single patients as well as transmission chains within individual hospitals and between neighbouring hospitals in Germany, and show that it can be detected by HierCC. We also demonstrate that HierCC can be used to identify bacterial populations at various epidemiological levels ranging from recent transmission chains through to epidemic and endemic spread, and relate these HierCC clusters to genotypes that were identified by PCR ribotyping and k-mer-based diversity analysis. These observations indicate that cgMLST and HierCC within EnteroBase can provide a common language for communications and interactions by the global community who is combatting disease caused by C. difficile.Results
Implementation of MLST schemes in EnteroBase
cgMLST in EnteroBase consists of a defined subset of genes within a whole-genome MLST scheme that represents all single-copy orthologues within the pan-genome of a representative set of bacterial isolates. To this end, we assembled the draft genomes of 5232 isolates of C. difficile from public short-read archives, and assigned them to ribosomal sequence types (rSTs) according to rMLST, which indexes diversity at 53 loci encoding ribosomal protein subunits on the basis of existing exemplar alleles at PubMLST . We then created a reference set of 442 genomes consisting of one genome of C. mangenotii , 18 complete genomes from GenBank, 81 closed genomes from our work and the draft genome with the smallest number of contigs from each of the 343 rSTs (https://tinyurl.com/Cdiff-ref). The Clostridioides pan-genome was calculated with PEPPA  and used to define a wgMLST scheme consisting of 13 763 genetic loci (http://enterobase.warwick.ac.uk/species/clostridium/download_data). EnteroBase uses the wgMLST scheme to call loci and alleles from each assembly and extracts the allelic assignments for the subsets corresponding to cgMLST, rMLST and seven-gene MLST from those allelic calls . The cgMLST subset consists of 2556 core genes, which were present in ≥98% of the reference set, intact in ≥94%, and were not excessively divergent (Fig. 1).Fig. 1.Criteria for inclusion in a cgMLST scheme of a subset of wgMLST genes based on their properties in a reference set of 442 genomes (https://tinyurl.com/Cdiff-ref). (a) Numbers of genes versus frequency (% presence) within the reference set. In total, 2634 genes satisfied the cut-off criterion of ≥98 % presence (dashed line). (b) Numbers of genes versus intact ORF (% intact ORF) within the 2634 genes from (a). Overall, 2560 genes satisfied the cut-off criterion of ≥94 % intact ORF (dashed line). (c) Frequency of allelic variants versus gene size among the 2560 genes from (b). The genetic diversity was calculated using the GaussianProcessRegressor function in the sklearn module in Python. This function calculates the Gaussian process regression of the frequency of genetic variants on gene sizes, using a linear combination of a radial basis function kernel (RBF) and a white kernel . The shadowed region shows a single-tailed 99.9% confidence interval (≤3 sigma) of the prediction. Altogether, 2556 loci fell within this area and were retained for the cgMLST scheme, while four were excluded due to excessive numbers of alleles.Click to view
Comparison of cgMLST and SNPs for analyses of transmission chains
We compared the numbers of cgMLST allelic differences and the numbers of non-recombinant SNPs in isolates from multiple epidemiological chains. These included 176 isolates from four patients with recurring CDI (C. difficile infection), 63 isolates from four transmission chains in multiple hospitals [14, 19, 26], and a comprehensive sample of 1158 isolates collected over several years in four hospitals in Oxfordshire, UK . A strong linear relationship (R2, 0.71–0.93) was found in all three analyses between the pairwise differences in cgMLST alleles and non-recombinant SNPs (Fig. S1, available in the online version of this article). The slope of the regression lines was close to 1.0, indicating a 1 : 1 increase in cgMLST allelic differences with numbers of SNPs. The same data were also investigated with cgMLST calculated with the commercial program SeqSphere , with similar correlation coefficients but a lower slope due to lesser discriminatory power of the SeqSphere cgMLST scheme (lower panels in Fig. S1).
Eyre et al.  concluded that direct transmission between two hospital patients can be detected because their bacterial genomes differ by two SNPs or less. Our analysis indicated that these transmission chains in the Oxfordshire dataset would also have been recognized by cgMLST in EnteroBase. Genomes that differed by two cgMLST alleles usually also differed by ≤2 SNPs according to a binary logistic regression model (probability=89%; 95% confidence interval, 88–89%) (Fig. 2). Of 3807 pairs of genomes with ≤2 allelic differences, 3474 also differed by ≤2 SNPs, yielding a positive predictive value of 91 % for identifying isolate pairs with ≤2 SNPs by EnteroBase cgMLST and a sensitivity of 62 % (≤2 cgMLST allelic differences were found in 3474 of 5707 pairs with ≤2 SNPs). The comparable values for SeqSphere were 78 % positive predictive value and 99% sensitivity.Fig. 2.Binary logistic regression model to determine the probability that two genomes are related at ≤2 SNPs, given a certain difference in their cgMLST allelic profiles, based on the Oxfordshire dataset . The number of SNPs was encoded as a binary dependent variable (1 if ≤2 SNPs, 0 if otherwise) and the number of allelic differences was used as a predictor variable.Click to view
We also compared the genetic distances between 242 genomes from Oxfordshire, which had been isolated during the initial 6 months and 916 genomes from the actual testing period (April 2008 to March 2011) . Overall, 35% (318/916) of the latter genomes matched at least one genome collected earlier by two or less EnteroBase cgMLST alleles and 34% (316/916) matched an earlier genome by ≤2 SNPs. The two sets of genomes were 89% concordant. Thus, cgMLST is equivalent to SNP analysis for detecting inter-patient transmission chains.
Hierarchical clustering for tracing local and regional spread
SNP analyses are computer intensive, and are only feasible with limited numbers of genomes . cgMLST-based relationships can be analysed for up to 100 000 genomes with GrapeTree, but analyses involving more than 10 000 genomes remain computer intensive . EnteroBase implements single-linkage hierarchical clustering (HierCC V1) of cgMLST data in pairwise comparisons at multiple levels of relationship after excluding missing data . These are designated as HC0 for hierarchical clusters of indistinguishable core-genome sequence types (cgSTs), HC2 for clusters with pairwise distances of up to two cgMLST alleles, etc. EnteroBase presents cluster assignments for C. difficile at the levels of HC0, HC2, HC5, HC10, HC20, HC50, HC100, HC150, HC200, HC500, HC950, HC200 and HC2500. Here we address the nature of the genetic relationships that are associated with these multiple levels of HierCC among 13 515 publicly available C. difficile genomes, and examine which levels of pairwise allelic distances correspond to epidemic outbreaks and to endemic populations.
In our analyses of 176 C. difficile isolates from four patients with two recurrent episodes of CDI, multiple genomes were assigned to patient-specific HC2 clusters, some of which were isolated from the initial episode as well as the recurrence 80–153 days later (Fig. 3, patients D, F and G; 4 to 36 isolates had been collected per episode; Table S1). For these patients, relapsing disease likely reflected continued colonization after initially successful therapy. However, some isolates from patient F differed by 12–21 cgMLST alleles from the bulk population (Fig. 3), which indicates that the patient was co-infected simultaneously with multiple related strains. In patient E, the two genomes from the two CDI episodes differed by >2000 allelic differences (Fig. 3), which indicates that the second incident of CDI represented an independent infection with an unrelated strain. Hence, discrimination between relapse and reinfection based on cgMLST appears to be straightforward except that two episodes of CDI might arise by reinfection with identical strains from an environment that is heavily contaminated with C. difficile spores . We note that the time intervals (16–22 weeks) investigated here exceeded the currently recommended threshold of 8 weeks for surveillance-based detection of CDI relapses [29, 30] but still yielded almost identical strains in three of four patients.Fig. 3.Minimum-spanning trees indicating the population structure of C. difficile in four patients with recurrent CDI episodes. Red, first episode; blue, second episode.Click to view
Our examinations of multiple local outbreaks have revealed individual, outbreak-specific HC2 clusters. However, it is also conceivable that multiple HC2 clusters might be isolated from a single epidemiological outbreak due to the accumulation of genetic diversity over time. Alternatively, multiple HC2 clusters within a single outbreak may represent the absence of crucial links due to incomplete sampling. Incomplete sampling of outbreaks is not unlikely because asymptomatic patients are only rarely examined for colonization with C. difficile [31–33] even though they may constitute an important reservoir for transmission. Indeed, some of the outbreaks investigated here did consist of more than one HC2 cluster (Fig. 4). For example, nine isolates from a recently reported ribotype 018 (RT018) outbreak in Germany  encompassed four related HC2 clusters, and outbreaks with RT027 and RT106 in a hospital in Spain  were each affiliated with two or three HC2 clusters (Fig. 4).Fig. 4.Neighbour-joining trees based on cgMLST showing the phylogenetic relationships among C. difficile isolates from previously published CDI outbreaks as indicated [14, 19, 26]. Nodes are coloured by HC2. CC, cgST complex, i.e. related at level HC150; RT, PCR ribotype. The scale, indicating one allelic difference, applies to all trees.Click to view
We identified 23 HC2 clusters encompassing 133 genome sequences in a dataset of 309 C. difficile genome sequences collected from CDI patients in six neighbouring hospitals in Germany. These HC2 clusters were associated with individual hospitals (Χ2, P=8.6×10−5; Shannon entropy, P=4.2×10−5) and even with single wards in these hospitals (Χ2, P=0.01; Shannon entropy, P=6.2×10−3). We investigated whether these HC2 clusters reflected the local spread of C. difficile within institutions by retrospective analyses of patient location data. Sixty-six patients (50 %) were found to have had ward contacts with another patient with the same HC2 cluster (median time interval between ward occupancy: 63 days; range, 0 to 521). These results are consistent with the direct transmission on wards of C. difficile isolates of the same HC2 cluster (Fig. 5). For patients such as P1 and P2 where the shared ward contacts were separated in time (Fig. 5), transmission may have occurred indirectly through asymptomatically colonized patients or from a common reservoir, such as environmental spore contamination [14, 31, 32]. We also detected 15 HC2 clusters that included isolates from two or more hospitals in the region. Subsequent analyses of patient location data confirmed that some of these HC2 clusters were associated with patient transferrals between the hospitals (Fig. 5). Hence, hierarchical clustering of C. difficile genome sequences in conjunction with retrospective analysis of patient movements revealed multiple likely nosocomial transmission events, none of which had been detected previously by routine surveillance.Fig. 5. Timelines of two transmission chains discovered retrospectively through inspection of files from CDI patients with closely related C. difficile isolates (HC2). Colours indicate hospital wards, ‘X’ indicates a diagnosis of CDI, and arrows indicate presumed transmission pathways. Minimum-spanning trees indicating genomic distances among C. difficile isolates are shown on the right. Upper panel: patient P1 was diagnosed with CDI in hospital H2 and transferred to hospital H3 15 days later. Another five and 6 days later, respectively, patients P2 in hospital H2 and P3 in hospital H3 got diagnosed with CDI with closely related strains. Both these patients were on the same wards as the initial patient, who probably had been the source for the pathogen. Since there was no temporal overlap between patient P2 and the other patients in hospital H2, the transmission may have occurred indirectly, possibly through environmental contamination. Lower panel: another putative transmission chain involved three patients that had shared time in hospital H2. Patients P4 and P5 got diagnosed with CDI on the same day after they had shared 7 days in this hospital, albeit on different medical wards. The third patient developed CDI with the same C. difficile cgST 4 days after being transferred to another hospital (h5), but had previously stayed at hospital H2 during the time when CDI got diagnosed in the first two patients. Since the three patients stayed on different wards in hospital H2, transmission presumably occurred indirectly.Click to view
Hierarchical clustering for identification of epidemic strains and endemic populations
The international epidemic spread of C. difficile over up to 25 years has been inferred previously on the basis of molecular epidemiology with lower resolution techniques . For multiple representatives of those epidemic strains in EnteroBase, the majority of these epidemic groups corresponded to HC10 clusters, including epidemic RT017  (HC10_17), the two fluoroquinolone-resistant lineages of RT027  (HC10_4, HC10_9), or livestock-associated RT078/126  (HC10_1) (Fig. 6).Fig. 6.Phylogenetic structure of three international C. difficile epidemics, each of which has spread for about 25 years [9, 11]. Within each epidemic, the majority of isolates are related at level HC10, as indicated by the colours. CC, cgST complex, i.e. related at level HC150; RT, PCR ribotype.Click to view
Endemic populations have also been described by ribotyping and phylogenetic analyses, some of which have acted as sources for the emergence of epidemic strains [2, 9]. Many endemic populations seem to be represented by HC150 clusters. Clustering at HC150 was well supported statistically (Fig. S2), and the frequency distribution of pairwise genomic distances indicated that multiple database entries clustered at <150 cgMLST allelic differences (Fig. S3). HC150 clusters also correlated well with k-mer-based classification . When applied to the dataset of 309 C. difficile genomes from six hospitals in Germany, the two methods implemented in EnteroBase and PopPUNK found 51 and 48 clusters, respectively, the majority of which coincided (adjusted Rand coefficient, 0.97).
A cgMLST-based phylogenetic tree of 13515 C. difficile genomes showed 201 well-separated HC150 clusters, each encompassing a set of related isolates, plus 209 singletons (Fig. 7). Because these HC150 clusters are based on cgMLST genetic distances, we refer to them as ‘cgST complexes’, abbreviated as CCs. Genomes from each of the major CCs have been collected over many years in multiple countries, indicating their long-term persistence over wide geographic ranges (Table 1).Fig. 7.Rapid-neighbour-joining phylogenetic tree based on cgMLST variation from 13 515 C. difficile genomes. Colours and numerals indicate CCs (HC150 clusters) with ≥10 entries, and information on predominant PCR ribotypes is provided in brackets.Click to view
We compared HC150 clustering with PCR ribotyping for 2263 genomes spanning 84 PCR ribotypes for which PCR ribotyping data were available in EnteroBase. These included 905 genomes, which we ribotyped (Table S2), as well as several hundred other genomes for which ribotype information was manually retrieved from published data. The correlation between HC150 clustering and ribotyping was high (adjusted Rand coefficient, 0.92; 95% confidence interval, 0.90–0.93). However, our analysis also revealed that PCR ribotypes did not always correspond to phylogenetically coherent groupings. PCR ribotypes 002, 015 and 018 were each distributed across multiple phylogenetic branches (Fig. 8). Furthermore, some genomes with indistinguishable cgMLST alleles were assigned to multiple ribotypes, including RT001/RT241, RT106/RT500 and RT126/RT078 (Fig. 8, Table 1). In these cases, both ribotypes occurred in several, closely related clades (Fig. 8), indicating that similar ribotype banding patterns had evolved multiple times. In contrast, HC150 clusters corresponded to clear-cut phylogenetic groupings within a phylogenetic tree of core genes (Fig. 8b).Fig. 8.Rapid-neighbour-joining phylogenetic tree based on cgMLST variation from 2263 C. difficile genomes, for which PCR ribotyping information is available. Upper panel: nodes are coloured by PCR ribotype as indicated. Lower panel: nodes are coloured by CC (HC150 clusters).Click to view
Higher population levels
HierCC can also identify clusters at still higher taxonomic levels, up to the levels of species and sub-species . In C. difficile, HC950 clusters seem to correspond to deep evolutionary branches (Fig. S4) and HC2000 clusters were congruent with the major clades reported previously , except that cluster HC2000_2 encompassed clade 1 plus clade 2 (Fig. S5). Finally, HC2500 may correspond to the subspecies level, because it distinguished between C. difficile and distantly related ‘cryptic clades’ (Fig. S6).Discussion
Infectious disease epidemiologists frequently seek to know if new isolates of bacterial pathogens are closely related to others from the different geographical origin, i.e. if they are part of a widespread outbreak. Unlike a previous cgMLST implementation , EnteroBase supports this goal by taking full advantage of rapidly growing, public repositories of short-read genome sequences . In contrast to short-read archives, however, where stored sequence data are not readily interpretable without specialized bioinformatic tools , EnteroBase enables contextual interpretation of a growing collection (18 254 entries as of June 2020) of assembled, quality-controlled C. difficile genome sequences and their associated metadata. At least the collection date (year), the geographic origin (country), and the source (host species) are available for the majority of database entries. Importantly, phylogenetic trees based on cgMLST allelic profiles from many thousand bacterial genomes can be reconstructed within a few minutes, whereas such calculations are currently prohibitively slow based on SNP alignments . Genome-sequencing reads from newly sampled C. difficile isolates can be uploaded to EnteroBase and compared to all publicly available genome data within hours, without requiring any command-line skills.
We demonstrate that the application of cgMLST to investigations of local C. difficile epidemiology yields results that are quantitatively equivalent to those from SNP analyses. This is a major advance because SNP analyses require specific bioinformatic skills and infrastructure, are time-consuming and not easily standardized . A web platform for centralized, automated SNP analyses on bacterial genomes is limited to food pathogens currently and does not offer any analyses on C. difficile genomes . Even though a cgMLST scheme for C. difficile had been published recently , its ability to identify closely related isolates and the inferred genomic distances was shown to be inferior to SNP analyses due to an excess of errors introduced by the de novo assembly of sequencing reads and a lack of per-base quality control . In EnteroBase, cgMLST is also based on de novo assembly, but EnteroBase uses Pilon  to polish the assembled scaffolds and evaluate the reliability of consensus bases of the scaffolds, thereby achieving comparable accuracy to mapping-based SNP analyses. When applied to a large dataset of C. difficile genomes from hospital patients in the Oxfordshire region (UK), cgMLST and SNP analysis were largely consistent (89% match) at discriminating between isolates that were sufficiently closely related to have arisen during transmissions chains from others that were epidemiologically unrelated.
After assembly, draft genomes contain missing data and many cgSTs have unique cgST numbers but are identical to other cgSTs, except for missing data. Hence, individual cgST numbers are only rarely informative. However, indistinguishable cgSTs are clustered in common HierCC HC0 clusters, which ignore missing data. In June 2020, the Clostridioides database contained >12 000 HC0 clusters, indicating that the majority of genomes was unique. Similarly, EnteroBase provides cluster designations at multiple levels of HierCC, enabling rapid identification of all cgSTs that are related at multiple levels of genetic distance. The data presented here shows that HierCC designations can facilitate communications between researchers, clinicians and public-health officials about transmission chains, epidemic outbreaks, endemic populations and higher phylogenetic lineages up to the level of subspecies.
EnteroBase cgMLST identified numerous HC2 clusters of strains in C. difficile isolates that seem to have arisen during transmission chains in six neighboring hospitals in Germany. These assignments were in part consistent with retrospective investigation of patient location data, although none of the nosocomial outbreaks (defined by German law as two or more infections with likely epidemiological connections [http://www.gesetze-im-internet.de/ifsg/]) had been detected previously by standard epidemiological surveillance by skilled clinical microbiologists. Recent publications propose that prospective genome sequencing of nosocomial pathogens should be applied routinely at the hospital level to guide epidemiological surveillance . Our data indicates that the combination of genome sequencing with cgMLST and HierCC may identify nosocomial transmission routes of C. difficile more effectively than presently common practice, and hence could help to reduce pathogen spread and the burden of disease. Reliable identification of transmission chains requires interpretation of pathogen genome sequence data in its epidemiological context, however .
HierCC will also enable comparisons to previously published data because we have provided a correspondence table between HC150 clusters and PCR ribotypes (Table 1). Rarefaction analysis indicated that the currently available genome sequences represent about two-thirds of extant HC150 (CC) diversity, which extrapolated to about 600 CCs (Fig. S7). At least some of this enormous diversity may be due to the occupation of multiple, distinct ecological niches, as exemplified by differential propensities for colonizing non-human host species (Table 1) [43, 44]. Individual CCs may also differ in their aptitudes for epidemic spread, as indicated by drastically different proportions of genomes assigned to HC2 chains: only 7% of CC141 were assigned to HC2 clusters versus 35% of CC34 and 77% of CC4 (Table 1). A full understanding of the population structure of C. difficile and its relationship to epidemiological patterns will require additional study because many of the clusters described here have not yet been studied or described. However, this task can be addressed by the global community due to the free public access to such an unprecedented amount of genomic data from this important pathogen.Methods
In total, 309 C. difficile isolates were collected at a diagnostic laboratory providing clinical microbiology services to several hospitals in central Germany. To assemble a representative sample, we included the first 20 isolates from each of six hospitals from each of three consecutive calendar years (Table S2). For the investigation of recurrent CDI, a set of 176 C. difficile isolates were collected in a diagnostic laboratory in Saarland, Germany. Here, primary stool culture agar plates were stored at 4 °C for 5 months to eventually enable the analysis of multiple plates representing episodes of recurrent C. difficile infection from individual patients, who had developed the recurrent disease by then and could be chosen with hindsight. It was attempted to pick and cultivate as many bacterial colonies from each selected plate as possible, resulting in 6 to 36 isolates per CDI episode (Table S1). In addition, we sequenced the genomes from 383 isolates that had been characterized by PCR ribotyping previously, including 184 isolates sampled from piglets , 71 isolates from various hospitals in Germany , and 108 isolates from stool samples collected from nursing home residents (unpublished; Table S2).
PCR ribotyping was performed as described previously , applying an ABI Prism 3100 apparatus for capillary electrophoresis and comparing banding patterns to the Webribo database (https://webribo.ages.at/).Whole-genome sequencing
For Illumina sequencing, genomic DNA was extracted from bacterial isolates by using the DNeasy Blood and Tissue kit (Qiagen), and libraries were prepared as described previously  and sequenced on an Illumina NextSeq 500 machine using a Mid-Output kit (Illumina) with 300 cycles. For generating complete genome sequences, we applied SMRT long-read sequencing on an RSII instrument (Pacific Biosciences) in combination with Illumina sequencing as reported previously . All genome sequencing data were submitted to the European Nucleotide Archive (www.ebi.ac.uk/ena) under study numbers PRJEB33768, PRJEB33779 and PRJEB33780
.SNP detection and phylogenetic analysis
Sequencing reads were mapped to the reference genome sequence from C. difficile strain R20291 (sequence accession number FN545816) by using BWA-MEM and sequence variation was detected by applying VarScan2 as reported previously . Sequence variation likely generated by recombination was detected through analysis with ClonalFrameML  and removed prior to determination of pairwise sequence distances  and to construction of maximum-likelihood phylogenetic trees with RAxML (version 8.2.9) .Genome assembly, quality control and wgMLST allele calling
Genomic data was processed by automated pipelines within EnteroBase, which were described in detail previously . Briefly, Illumina sequencing reads were assembled by using Spades v3.10  and assemblies were improved by applying Pilon . To pass quality control, assemblies were required to comply with the following thresholds: total length, 3.6 to 4.8 Mbp; N50, ≥20 000; the number of contigs, ≤600; number of unresolved nucleotides, ≤3%; the proportion of Clostridioides sequences, >65 % (as determined by Kraken with MiniKraken database ). Assemblies were aligned to exemplar alleles by using blastn  and the usearch module UblastP , and allele numbers, STs, and HC numbers assigned by using the EnteroBase module MLSType . All EnteroBase tools are available at https://github.com/zheminzhou/EToKi.
To determine the probability that two genomes are related at ≤2 SNPs, given a certain difference in their cgMLST allelic profiles, we inferred a logistic regression model using R (, pp. 593–609). Genomic relatedness was encoded as a binary response variable (1 if ≤2 SNPs, 0 if otherwise) and the number of core-genome allelic differences was used as a predictor variable. We applied this model to a dataset of 1158 genome sequences from a previous study, representing almost all symptomatic CDI patients in Oxfordshire, UK, from 2007 through 2011 . While that original study had encompassed a slightly larger number of sequences, we restricted our analysis to the data (95 %) that had passed quality control as implemented in EnteroBase . We used the SNP data from Eyre’s report .
The hierarchical single-linkage clustering of cgMLST sequence types was carried out as described  for all levels of allelic distances between 0 and 2556. We searched for stable levels of differentiation by HierCC according to the Silhouette index , a measure of uniformity of the divergence within clusters. The Silhouette index was calculated based on d^’, a normalized genetic distance between pairs of STs, which was calculated from their allelic distance d as follows: d^’=1-(1-d)^(1/l), where l is the average length (937 bp) of the genes in the cgMLST scheme.
We further evaluated the ‘stability’ of hierarchical clustering using two other criteria. The Shannon index is a measure of diversity in a given population. The Shannon index drops from nearly 1 in HC0, because most cgSTs are assigned to a unique HC0 cluster, to 0 in HC2500, which assigns all sequence types to one cluster. The gradient of the Shannon index between the two extremes reflects the frequencies of coalescence of multiple clusters at a lower HC level. Thus, the plateaus in the curve correspond to stable hierarchical levels, where the Shannon index does not change dramatically with HC level. We also evaluated the stability of hierarchical clustering by pairwise comparison of the results from different levels based on the normalized mutual information score  (Fig. S3).
For clustering C. difficile diversity with PopPUNK , we used a sketch size of 105 and a K value (maximum number of mixture components) of 15. Of note, the resulting number of clusters for the tested dataset was identical for all K between 15 and 30.
To estimate concordance between cgMLST-based hierarchical clustering and PCR ribotyping or PopPUNK clustering, respectively, we calculated the adjusted Rand coefficient  by using the online tool available at http://www.comparingpartitions.info/. To test statistical associations of HC2 clusters with specific hospitals and hospital wards, respectively, we compared Χ2 values and normalized Shannon entropy values (R package ‘entropy’ v.1.2.1) from contingency tables containing real isolate distributions (Table S3) and randomly permuted distributions (n=1000), by using the non-parametric, two-sided Mann–Whitney U test (R package ‘stats’ v.3.5.0).
This work was partially funded by the German Center for Infection Research (DZIF), by the Federal State of Lower Saxony (Niedersächsisches Vorab VWZN2889/3215/3266), by the EU Horizon 2020 programme (grant agreement number 643476), the Wellcome Trust (098051), and the UK Medical Research Council (PF451). EnteroBase development was funded by the BBSRC (BB/L020319/1) and the Wellcome Trust (202792/Z/16/Z), and the salary of Z.Z. was also provided by The Wellcome Trust. The funders had no role in the study design, preparation of the article, or decision to publish.
We thank Vera Junker, Simone Severitt, Nicole Heyer and Carola Berg for excellent technical assistance, Johannes Sikorski for help with R, and David Eyre for supplying SNP data from his 2013 paper in tabular format.Author contributions
M.A. and U.N. designed the study. M.F., Z.Z., M.S., J.P.M., M.G., T.R., B.B., C.S., J.O., M.B., A.I., L.v.M., F.K., S.G.F., O.Z., U.G., M.A., and U.N. analyzed the data. L.v.M., C.S., O.Z., U.G., R.C., R.d.C., N.K., T.D.L., T.A.K. and S.N. acquired bacterial isolates and sequence data. M.F., Z.Z., M.A., and U.N. wrote the manuscript and all authors reviewed and edited the manuscript and approved the final version.
– SER-109 met Phase 3 primary endpoint, showing a highly statistically significant 30.2% absolute reduction in the rate of C. difficile infection recurrence compared to placebo –
– SER-109 was well tolerated, with a safety profile comparable to placebo –
– Efficacy results substantially exceeded FDA regulatory guidance to support BLA filing as a single pivotal trial; Company to meet with the agency to discuss filing for product approval as soon as possible –
–Positive SER-109 Phase 3 data provide validation for Seres’ microbiome therapeutics platform and further development of its pipeline of product candidates –
Seres Therapeutics, Inc. reported on August 10, 2020, positive topline results from the pivotal Phase 3 ECOSPOR III study evaluating its investigational oral microbiome therapeutic SER-109 for recurrent C. difficile infection (CDI). The study showed that SER-109 administration resulted in a highly statistically significant absolute decrease of 30.2% in the proportion of patients who experienced a recurrence in CDI within eight weeks of administration versus placebo, the study’s primary endpoint. 11.1% of patients administered SER-109 experienced a CDI recurrence, versus 41.3% of placebo patients. The study results were equally compelling when characterized by the alternative metric of sustained clinical response, where 88.9% of patients in the SER-109 arm achieved this objective.
The study’s efficacy results exceeded the statistical threshold previously provided in consultation with the U.S. Food and Drug Administration (FDA) that could allow this single clinical study to fulfill efficacy requirements for a Biologics License Application (BLA). The SER-109 safety results were favorable, with an adverse event profile comparable to placebo.
“We are extremely pleased with these highly clinically meaningful SER-109 Phase 3 study results, greatly exceeding the statistical threshold provided by the FDA. Based on our prior discussions with the FDA, we believe this trial should provide the efficacy basis for submitting an application for product approval. We look forward to meeting with the FDA as soon as possible to discuss the regulatory path forward with the goal of bringing SER-109 to patients as a first-in-class microbiome therapeutic,” said Eric D. Shaff, President and Chief Executive Officer of Seres. “Our results represent the first-ever positive pivotal clinical study results for a targeted microbiome drug candidate. We believe these Phase 3 data provide strong validation for our underlying microbiome therapeutics platform, which has been the scientific basis for the Company, as well as persuasive clinical evidence supporting our other active pipeline programs.”
“We would like to thank all those who participated in this landmark study. Based on these highly positive SER-109 ECOSPOR III results, we believe that this novel microbiome therapeutic candidate could potentially provide a much-needed effective oral treatment option for the approximately 170,000 patients in the U.S. that suffer from recurrent CDI annually,” said Lisa von Moltke, M.D., FCP, Chief Medical Officer of Seres. “Seres applied a data-driven and scientifically rigorous approach to develop SER-109. The proprietary scientific learnings we have obtained continue to drive our overall R&D efforts and the advancement of our other ongoing microbiome therapeutic programs.”
“Recurrent C. difficile infection is a serious disease that devastates patients’ quality of life, and in many severe cases may result in a patient’s death. Today’s treatment options have important shortcomings related to efficacy, safety and route of administration, and novel approaches that target the root causes of the disease are urgently needed. The SER-109 Phase 3 results are highly impressive and represent an exceptional advance in the fight against this disease. I believe that SER-109 has the potential to fundamentally transform the treatment of recurrent C. difficile infection,” said Mark Wilcox, M.D., Professor of Medical Microbiology, University of Leeds.
ECOSPOR III Study Design and Results
The ECOSPOR III study (ClinicalTrials.gov identifier: NCT03183128) is a multicenter, randomized, placebo-controlled study that enrolled 182 patients with multiply recurrent CDI. Patients were randomized 1:1 to receive either SER-109 or placebo, after standard of care antibiotic treatment. SER-109, or placebo, was administered orally for three consecutive days. All patients were required to have a positive C. difficile toxin diagnostic test both at study entry and in the case of suspected recurrence to ensure the selection of individuals with active disease and to confirm the accuracy of the primary endpoint.
The primary efficacy endpoint of ECOSPOR III was the proportion of patients with recurrent CDI at up to eight weeks following administration of SER-109 or placebo. As a secondary endpoint, patients are evaluated for CDI recurrence through 24 weeks post-treatment, and the Company plans to present those results at a future date.
SER-109 met the study’s primary endpoint with a significantly lower recurrence rate of 11.1% in SER-109 patients versus 41.3% in placebo patients at eight weeks; p<0.001 tested at the one-sided 0.25 level. Patients administered SER-109 experienced a 30.2% lower rate of recurrence, on an absolute basis, compared to placebo. The SER-109 treatment arm relative risk was 0.27 (95% CI=0.15 to 0.51) versus placebo. The ECOSPOR III recurrence rates translate into a sustained clinical response rate of 88.9% versus 58.7% with SER-109 and placebo, respectively. The SER-109 Number Needed to Treat (NNT) was approximately 3.
In prior discussions, the FDA communicated that demonstration of a statistically very persuasive efficacy finding in the ECOSPOR III primary endpoint, defined as demonstrating a 95% upper confidence level of relative risk lower than 0.833, could support a BLA submission on the basis of this single study. The results of ECOSPOR III demonstrated a SER-109 relative risk of 0.27 (95% CI=0.15 to 0.51) compared to placebo. As a result, Seres believes that this study should support the efficacy basis for BLA submission. SER-109 has obtained FDA Breakthrough Therapy and Orphan Drug designations.
SER-109 was well tolerated, with no treatment-related serious adverse events (SAEs) observed in the active arm, and an adverse event profile similar to placebo. The overall incidence of patients who experienced AEs during the eight-week study period was similar between SER-109 and placebo arms. The most commonly observed treatment-related AEs were flatulence, abdominal distention and abdominal pain, which were generally mild to moderate in nature, and these were observed at a similar rate in both the SER-109 and placebo arms.
A SER-109 open-label study is ongoing ( clinicaltrials.gov identifier: NCT03183141) at selected clinical sites that participated in the ECOSPOR III study, and the Company may initiate the program at additional clinical sites. The FDA has previously indicated that SER-109 administration to at least 300 patients, consistent with standard FDA guidance, would be required to support BLA submission. The ongoing SER-109 open-label study is continuing to contribute to the SER-109 safety database.
The Company plans to immediately request a Breakthrough Therapy designation meeting with the FDA to discuss the requirements to submit a BLA seeking regulatory approval of SER-109. Given the favorable efficacy and safety results seen in ECOSPOR III, the safety results observed in prior SER-109 clinical studies, and the critical unmet need for a therapeutic option for recurrent CDI patients, the Company plans to discuss with the FDA the safety data requirements for a BLA filing.
Seres continues to advance its commercial readiness for the potential launch of SER-109. In June 2020, Seres appointed Terri Young, Ph.D., R.Ph., as Chief Commercial and Strategy Officer. The Company has been conducting activities to support successful future potential commercialization. Seres believes that the commercial opportunity for SER-109 could be substantial, given the dire need for an effective, safe, oral therapeutic, and the strength of the SER-109 Phase 3 study results.
Conference Call Information
Seres’ management will host a conference call today, August 10, 2020, at 8:30 a.m. ET. To access the conference call, please dial 844-277-9450 (domestic) or 336-525-7139 (international) and reference the conference ID number 3216859. Accompanying slides will be posted on the Seres website ahead of the conference call. To join the live webcast, and to view the accompanying slides, please visit the “Investors and Media” section of the Seres website at www.serestherapeutics.com.
A webcast replay will be available on the Seres website beginning approximately two hours after the event and will be archived for approximately 21 days.
SER-109 is an investigational, oral, biologically-derived microbiome therapeutic that is designed to reduce recurrence of C. difficile infection (CDI), enabling patients to achieve a sustained clinical response by breaking the vicious cycle of CDI recurrence and restoring the diversity of the gastrointestinal microbiome. SER-109 is a consortium of purified bacterial spores of multiple Firmicute species, manufactured by fractionating targeted bacteria from the stool of healthy human donors with further steps to inactivate potential pathogens. The FDA has granted SER-109 Breakthrough Therapy designation and Orphan Drug designation for the treatment of CDI.
SER-109 is fundamentally distinct from fecal microbiota transplantation (FMT). SER-109 is comprised of a highly-purified consortia of spore-based commensal bacteria and designed to be manufactured in accordance with Good Manufacturing Practice conditions using stringent standards to ensure product quality and consistency. To support product safety, Seres utilizes a unique manufacturing process that inactivates numerous potential pathogens, including species of non-spore bacteria, such as Escherichia coli, and viruses such as SARS-CoV-2.
About C. difficile Infection (CDI) and Current Treatments
C. difficile infection (CDI) is one of the top three most urgent antibiotic-resistant bacterial threats in the U.S., according to the Centers for Disease Control, and is a leading cause of hospital-acquired infection in the U.S. It is responsible for the deaths of approximately 20,000 Americans each year. CDI is associated with debilitating diarrhea, which significantly impacts quality of life in every functional domain. Since the discovery of C. difficile more than four decades ago, vancomycin has been the most commonly used drug for patient management. Current approaches provide only modest improvements in sustained clinical response rates, leaving behind a significant pool of patients with recurrent disease. Unapproved FMT, used in cases that are not responsive to approved drugs, remains poorly characterized clinically and has been associated with serious safety concerns, including the transmission of bacterial pathogens and the potential transmission of viruses such as SARS-CoV-2, the virus that causes COVID-19. The recent quarantine and shipping hold of FMT from a major stool bank highlights the urgent need for an approved effective and safe treatment for recurrent CDI.
About Seres Therapeutics
Seres Therapeutics, Inc., (Nasdaq: MCRB) is a leading microbiome therapeutics platform company developing a novel class of multifunctional bacterial consortia that are designed to functionally interact with host cells and tissues to treat disease. Seres’ SER-109 program achieved the first-ever positive pivotal clinical results for a targeted microbiome drug candidate and has obtained Breakthrough Therapy and Orphan Drug designations from the FDA. The SER-109 program is being advanced for the treatment of recurrent C. difficile infection and has the potential to become a first-in-class FDA-approved microbiome therapeutic. Seres’ SER-287 program has obtained Fast Track and Orphan Drug designations from the FDA and is being evaluated in a Phase 2b study in patients with active mild-to-moderate ulcerative colitis. Seres is developing SER-401 in a Phase 1b study in patients with metastatic melanoma, SER-301 for ulcerative colitis and SER-155 to prevent mortality due to gastrointestinal infections, bacteremia, and graft versus host disease. For more information, please visit www.serestherapeutics.com
SOURCE: Seres Therapeutics