Higher genome variability within metabolism genes associates with recurrent Clostridium difficile infection
BMC Microbiology volume 21, Article number: 36 (2021)
Clostridium difficile (C. difficile) is a major source of healthcare-associated infection with a high risk of recurrence, attributable to many factors such as usage of antibiotics, older age and immunocompromised status of the patients. C. difficile has also a highly diverse genome, which may contribute to its high virulence. Herein we examined whether the genome conservation, measured as non-synonymous to synonymous mutations ratio (dN/dS) in core genes, presence of single genes, plasmids and prophages increased the risk of reinfection in a subset of 134 C. difficile isolates from our previous study in a singly hemato-oncology ward.
C. difficile isolates were subjected to whole-genome sequencing (WGS) on Ion Torrent PGM sequencer. Genomes were assembled with MIRA5 and annotated with prokka and VRprofile. Logistic regression was used to asses the relationship between single gene presence and the odds of infection recurrence. DN/dS ratios were computed with codeml. Functional annotation was conducted with eggNOG-Mapper.
We have found that the presence of certain genes, associated with carbon metabolism and oxidative phosphorylation, increased the odds of infection recurrence. More core genes were under positive selective pressure in recurrent disease isolates – they were mostly associated with the metabolism of aminoacids. Finally, prophage elements were more prevalent in single infection isolates and plasmids did not influence the odds of recurrence.
Our findings suggest higher genetic plasticity in isolates causing recurrent infection, associated mainly with metabolism. On the other hand, the presence of prophages seems to reduce the isolates’ virulence.
Clostridium difficile (reclassified in 2016 into a new Clostridiodes genus, along with Clostridium mangenotii, with which it shares a 94.7% similarity within 16 s rRNA gene ) is an anaerobic, spore-forming, Gram-positive bacterium, prevalent in the environment, as well as human gastrointestinal tract: it is mainly present in infants [2, 3], in whom the asymptomatic colonization occurs. On the other hand, in adults, C. difficile colonization is characterized by life-threatening infection, with symptoms ranging from moderate diarrhea to severe colitis and/or megacolon [4, 5]. In Europe, the majority (74.6%) of C. difficile cases are healthcare-associated (HA). The mean incidence per 10,000 patient-days is at 2.38, but some countries, such as Estonia, Lithuania and Poland present with much higher numbers . In large European hospital surveys from 10%  to 16%  Clostridium difficile infection (CDI) cases are associated with major complications, requiring admission to the intensive care unit and resulting in the death rate of 7 and 4% respectively.
Another major concern when dealing with CDI is its tendency of recurrence which, according to the European Society of Clinical Microbiology and Infectious Diseases, is defined as a relapse of CDI clinical symptoms within 2–8 weeks of successful treatment of the initial episode. The recurrent CDI may be due to a relapse of the previous CDI by the same strain or reinfection by a different strain [8,9,10]. While distinguishing of recurrence due to relapse or from recurrence due to reinfection is not feasible in daily practice, the method of choice in this distinction is bacterial genotyping . Reported recurrence rates vary between 5 and 50%, but most of them are between 10 and 20% . Various recurrence risk factors have been identified, including continued use of antibiotics not associated with CDI treatment , particularly cephalosporins [13, 14], older age [12, 15], HA diseases [15, 16], length of hospitalization  and usage of gastric acid suppressors such as proton pump inhibitors [10, 12, 17, 18]. Immunocompromised patients also typically present a higher risk of infection recurrence [19, 20].
In addition, the high plasticity of C. difficile genome may also account for its virulence and notorious recurrence. Sequenced C. difficile genomes’ size typically ranges from 4.1 to 4.3 Mb [21,22,23,24] and is much larger not only than most of the related species but also most of the Firmicutes phylum. A large genome can be a sign of exceptional adaptability to various conditions, often for prolonged periods . Indeed, a large number of differentially expressed genes during infection were found to be associated with adaptation mechanisms such as stress response and sporulation . The C. difficile 630 genome was also found to contain an unusually high proportion (11%) of mobile elements, including transposons and prophages . The horizontal gene transfer, occurring through these elements is particularly important in the acquisition of resistance to various antimicrobial agents . Finally, the C. difficile genome can be altered through point mutations and inversions: in the comparison between 3 strains: 630, R20291 and CD196 39 variations such as these were found [23, 28].
In this context, the aim of our work was twofold. Firstly, to investigate and establish the core genome and pangenome of Clostridial species recovered from patients with single and multiple CDIs hospitalized in a single hemato-oncology ward through a period of ten years, both on genetic and functional level. Secondly, to investigate the dN/dS ratio in the core genome. This data was used to compare strains, which cause recurring infections to those, which affected the patient only once.
Clostridium difficile core genome
There were 965 core genes discovered (i.e. present in more than 90% of samples) shared between isolates from ST1 and ST42. Apart from 208 proteins with poorly characterized function, the most abundant COG categories concerned metabolism (including carbohydrate and amino acid metabolism), information storage and processing (transcription and translation) and cellular processes and signaling (signal transduction and cell wall biogenesis) (Fig. 1, Supplementary Table 1). 135 KEGG Pathways (on a third level of classification) were represented in the common core genome (Supplementary Table 2). The most represented pathways (i.e. the ones with the highest ratio of present genes to the total number of KEGG orthologies) are associated with membrane transport, aminoacid, carbohydrate and lipid metabolism, bacterial cell motility, genetic information processing, energy metabolism and metabolism of cofactors and vitamins (Fig. 2). 62% pathways present within the core genome are associated with metabolism.
Specific gene presence as a predictor of recurrence
In the logistic regression model, there were 5264 genes tested; while 192 reached statistical significance at a nominal p-value of 0.05, none of them were significant after correction for multiple hypothesis testing with the FDR method (Supplementary Table 3). 7 pathways are enriched in gene set which gives higher odds of infection recurrence, however, in only one of them the enrichment score reaches the highest value at rank smaller than 192: Oxidative phosphorylation (Table 1, Supplementary Table 4).
Gene conservation differences between recurrent and single infection
515 core genes were included in this analysis. Sequences were analyzed separately for recurrent and one-time infections. For recurrent infections, 65 genes had sites under positive selection, while for single infections this number was at 17. 25 genes under positive selection only in recurrent sequences could be functionally annotated with KEGG pathways. Most of them are associated with metabolism, mainly of aminoacids and secondary metabolites. Other genes of interest include toxin B and cheC, involved in bacterial chemotaxis (Table 2, Supplementary Table 5).
Mobile genetic elements
Three different plasmids were discovered in our strains: pCD6, pCDBI1 and DSM 1296. Plasmids were present in 5% (5 out of 98) strains which cause recurrent infection and in 11% strains which did not (4 out of 36). This difference is not statistically significant (p-value 0.25, Fisher’s exact test, Table 3).
Eleven prophage elements differed in frequency between single and recurrent infections at a nominal p-value < 0.05, however, none of them remained significant after multiple hypothesis testing correction. Most of them originated from phages CD119 and C2 (Table 4, Supplementary Table 6).
The C. difficile genome is subject to constant changes, as estimated, it acquires between 1 and 2 mutations per genome per year [25, 29]. However, not all mutations are viable, and some clones may become subject to purifying selection. The selective pressure can be described with the dN/dS ratio, i.e. the ratio of non-synonymous to synonymous mutations. The dN/dS ratio significantly smaller than 1 suggests strong purifying selection, while in most C. difficile genomes the reverse situation is observed where the dN/dS is actually higher than 1 [24, 30]. This suggests less efficient purging of novel mutations, possibly contributing to C. difficile high genetic diversity and adaptability. With such plasticity and diversity, it is difficult to establish the exact size of C. difficile core genome and pangenome. Usually, the orders of magnitude of about 1000 genes are given [31, 32], but some researchers give figures as high as 3000 , which would be the most of C. difficile genome. Nevertheless, the estimated size of the core genome usually varies between 16%  and 24% , which is much lower than most bacterial species. For instance, in pathogenic Streptococcus agalactiae, the core genome constitutes about 80% of the whole genome, in Helicobacter pylori 77 and 46% in Streptococcus pneumoniae . On the other hand, the size of the core and essential genome was estimated to be composed of 404 genes , a number comparable to other bacterial species, such as Pseudomonas aeruginosa (321 genes ) and Yersinia pestis (about 500 genes ).
The C. difficile core genome is usually estimated to comprise about 1000 genes , involved mainly in pathways related to metabolism (of aminoacids and carbohydrates), genetic information processing, cell motility and signal transduction [31, 34], the unsurprising functions in the core genome. Additionally, many clostridial core genes are associated with virulence, including toxins, cell surface proteins, flagellar proteins and antibiotic resistance factors . Our study is in line with previous findings, with the core genome of 965 genes, present in the most prevalent strains. Apart from typical house-keeping pathways, we have identified several KEGG pathways associated with virulence such as beta-Lactam and vancomycin resistance, biofilm formation and flagellar assembly.
It is believed that highly adaptive metabolism is one of the key contributors to C. difficile virulence. C. difficile has two main energy sources: aminoacids and sugars. Some aminoacids (such as leucine, valine, proline) contribute to ATP formation via the so-called Stickland pathway  while other aminoacids (including cysteine, threonine, serine) and sugars contribute to energy production via central carbon metabolism and TCA cycle . Furthermore, C. difficile exhibits some autotrophic characteristics, including genes from the Wood-Ljungdahl pathway in 4 sequenced genomes that allow an autotrophic growth by generating energy from CO2 and H2 via this pathway . Production of toxin A and B is also correlated with alterations in central carbon metabolism with fluxes changing from butanoate to lactate synthesis . In our work, genes whose presence increased the odds the reinfection are mainly associated with metabolism and energy production: Oxidative phosphorylation, Carbon metabolism, 2-Oxocarboxylic acid metabolism and Carbon fixation in photosynthetic organisms. This may suggest that infection recurrence is associated with altered metabolism and alternative means of energy production rather than the presence of additional virulence factors. However, the practical significance of this discovery in the aspect of new antimicrobial targets for C. difficile remains to be uncovered. Historically, bacterial central metabolism (including carbon metabolism) has not been reported as a potential source of new antimicrobial targets, since it has been believed that homology between crucial microbial and human enzymes is simply too high . On the other hand, recent studies identified potential antimicrobial drug targets within carbon metabolism and fixation pathways in MRSA  and Mycobacterium tuberculosis . While the aforementioned works are based on in silico methods, their practical utility remains to be proven.
Virulence-conferring plasmids are common in enteropathogenic bacteria, such as Shigella spp  and Escherichia coli . They contribute to bacterial virulence by carrying genes associated with resistance against antimicrobial agents (such as plasmids R100 in Shigella flexneri 2b which contributes to resistance to sulfonamides, chloramphenicol, tetracyclines and streptomycin ) as well as host cell adhesion and invasion (plasmid pO157 in E. coli ). Two main plasmids were described in C.difficile: pCD6  and pCD630  – both are relatively small (less than 10 kb) in comparison with virulence-conferring plasmids, such as pO157, which is 93.6 kb in size . Recently, plasmids larger than 40 kb were discovered . However, while historically plasmid-typing was found to be useful in tracing and typing nosocomial CDI [53, 54], no clear associations between virulence and plasmid presence could be drawn for C. difficle [52, 55]. Only recently a plasmid pMETRO, conferring resistance to metronidazole has been discovered . In our study, plasmids were discovered in 6.7% isolates. None of the isolates were resistant to metronidazole and unsurprisingly pMETRO was not discovered in any of them. There was no statistically significant difference between plasmid fraction detected in one-time and recurrent infection isolates. Therefore, we believe our study reinforces the hypothesis that plasmids contribute little to C. difficile virulence and recurrence.
Prophages contribute to the evolution and virulence of most bacterial pathogens, including virulence and recurrence of C. difficile [57,58,59]. Prophages are abundant within C. difficile genome – up to 2018, at least 26 mobile element sequences were described . They manifest a large variety of functions in C. difficile: while CD119 represses expression of five clostridial pathogenicity locus (PaLoc) genes , other prophages may promote virulence: upon infection with CD38–2, up to two-fold rise in toxin A and B was detected in hypervirulent BI/NAP1/027 (ST1) strain . Another prophage, Semix9P1, was determined to carry a fully functional binary toxin . In our study, prophage fragments were discovered, coming from two phages: CD119 and C2. Phage CD119 contains a gene, encoding RepR protein, which downregulates the expression of toxins A and B indirectly controlling the expression of tcdR, the toxin gene regulator . RepR gene was found in almost 20% of single infection isolates and only in 5% recurrent infection isolates. This may suggest weakened virulence due to the prophage infection at least in some of the single infection isolates. On the other hand, phage C2 infection was found to affect the measured levels of toxin B in C. difficile isolates through the production of holins, proteins that disrupt the membrane and increase bacterial secretion . However, after correction for multiple hypothesis testing none of 11 prophage elements uncovered in this study, including C2 phage originated holins, differed in frequency between single and recurrent infections.
Finally, high adaptability and increased virulence may be attributed to the beneficial point mutations which do not become subject to purifying selection. While it is expected for the core genome to be highly conserved, some of the clostridial core genes were found to be under positive selection. For instance, He et al.  identified 12 such sequences, including membrane proteins and response regulators. The dN/dS ratios for core genomes were higher than 1 for both strains analyzed by Murillo et al.  Recently, Kumar et al.  have proposed that C. difficile is undergoing active speciation and characterized genes under positive selective pressure which were associated with sporulation and sugars' metabolism. In our study, more genes were found to be under positive selective pressure in recurrent infection isolates that in single ones. While this may point to higher genetic adaptability of recurrent isolates, in this case other explanations should be also considered. First of all, we had a larger number of sequences from recurrent isolates, which increases probability of detecting positive selection. In addition, in closely related lineages the dN/dS ratios were proved to be higher since the purifying selection did not have time to purge mutations . Nevertheless, genes under positive selection in recurrent infections in our study seem to share some of the characteristics with those designated by He et al. and Kumar et al. [24, 32]. In line with He et al., we have ABC transporters (cdd3, fatC and opuCC) and two-component system members (regB) with sites under positive selection. We have also observed changes in metabolism, but they concern aminoacids rather than sugars, similar to Kumar et al. study. Interestingly, toxin B also has been found to be under positive selective pressure in recurrent isolates. Toxin B was found to be crucial in C. difficile virulence , with toxinB(+) toxinA (−) mutants being fully virulent, while in the reverse situation the virulence is attenuated . Genetic variability within toxin sequences is a well-known phenomenon , with 34 toxinotypes currently defined . While our results may suggest an existing selective pressure on toxin B gene, it is also worth noting that most of the toxinotypes are known since the beginning of research on the subject and on a larger scale prevalence of alternative toxinotypes is attributable to local outbreaks .
To conclude, we have managed to thoroughly analyze how genetic mobility influences infection recurrence in CDI. We have confirmed the lack of significance of plasmid-related virulence, as well as reinforced the role of prophages in the virulence-related mechanisms. This seems to be of particular importance since phage therapy seems like a beneficial alternative due to limited antibiotics available for the treatment of CDI . We have also observed changes in metabolism-related genes, both in prevalence (shell genes), as well as in conservation (core genes).
Research involving human data was performed in accordance with the Declaration of Helsinki.
The study was approved by the Maria Skłodowska-Curie National Research Institute of Oncology Ethics Committee (number 40/2018). In line with the opinion of the Bioethics Committee at Maria Skłodowska-Curie National Research Institute of Oncology our study did not require informed consent for the following reasons: This is a retrospective study describing the genetic differences between C. difficile strains but not between patients; bacterial strains were isolated during routine diagnostics and then banked over the course of one to 10 years; most of these patients are already dead.
As reported recently , between 2008 and 2018, all patients hospitalized at the Department of Lymphoma with healthcare-associated diarrhea (defined as ≥3 stools within a 24-h period arising over the third day after hospital admission) underwent testing at the Department of Clinical Microbiology to detect pathogenic C. difficile toxins A and B. Tests were performed using the C. difficile TOX A/B kit (TechLab). Subcultured single colonies from 134 available culture-positive isolates were subjected to whole-genome sequencing (WGS) on Ion Torrent PGM sequencer. Of these, 36 isolates were recovered from patients (18 women and 18 men) with a single CDI, and 98 were recovered from 44 patients (19 women and 25 men) with multiple CDIs. Multi-locus sequence typing results are taken from this publication as well.
Genome assembly and annotation
The sequenced genomes were assembled with MIRA 5 (https://sourceforge.net/projects/mira-assembler/) genome assembler , using parameters set specific for Ion Torrent sequencing technology. The assembled sequences were annotated with prokka  version 1.13 (https://github.com/tseemann/prokka), using a minimal contig length of 1000, proteins from RT027 CD196 strain as a list of trusted proteins and a “metagenome” option to improve annotation in case of large genome fragmentation. The pan-genome calculation was conducted with roary  version 3.12 (https://sanger-pathogens.github.io/Roary/).
Functional annotation of coding sequences and mobile genetic elements
The functional annotation of identified CDS was conducted with eggNOG-Mapper  version 2.0 (https://github.com/eggnogdb/eggnog-mapper), using eggNOG  categories as well as KEGG  pathways. All visualizations were performed with R package ggplot2 . Mobile elements, specifically the prophages were annotated with VRProfile  (https://tool-mml.sjtu.edu.cn/STEP/STEP_VR.html). Plasmids were identified with PlasmidSeeker  (https://github.com/bioinfo-ut/PlasmidSeeker).
Gene presence as a predictor of disease recurrence
The influence of a single gene on odds of disease recurrence was assessed with a logistic regression model, with ST as a confounding variable. The analysis was conducted only for genes present in more than 15% and less than 90% of cases. The p-values from this model on a PHRED scale (i.e. transformed with negative logarithm with base 10) served as metric for GSEA (Gene Set Enrichment Analysis) conducted with function from Cluster Profiler  package, with 10,000 permutations and maximum gene set size of 200. The metric was negative for genes which decreased the odds of recurrence.
Gene conservation in recurrent and one-time infections
In order to compute gene conservation, the CDSs of core genes from STs 1 and 42 were translated into proteins with translate function from BioStrings R package  version 2.46 and option to resolve ambiguous codons. The sequences were then aligned with msa function from R package msa  version 1.14 (using default parameters and default aligner ClustalW ). The protein alignment was then converted to codon alignment with pal2nal script . The presence of genetic recombination was verified with PhiPack  software and analysis was not continued only if the sequences passed 2 tests present in the package. The dN/dS ratios among different sites were then assessed with codeml (part of PAML4  package - http://abacus.gene.ucl.ac.uk/software/paml.html), using a comparison of two models - nearly neutral (designed M1a in PAML manual - http://abacus.gene.ucl.ac.uk/software/pamlDOC.pdf) and positive selection (M2a). P-values were adjusted for multiple testing with Benjamini-Hochberg FDR correction .
Availability of data and materials
The datasets generated for this study can be found as raw fastq files in Sequence Read Archive with accession number PRJNA608241 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA608241/).
Lawson PA, Citron DM, Tyrrell KL, Finegold SM. Reclassification of Clostridium difficile as Clostridioides difficile (hall and O’Toole 1935) Prévot 1938. Anaerobe. 2016;40:95–9.
Leffler DA, Lamont JT. Clostridium difficile Infection. N Engl J Med. 2015;372:1539–48.
Czepiel J, et al. Clostridium difficile infection: review. Eur J Clin Microbiol Infect Dis. 2019;38:1211–21.
Ofosu A. Clostridium difficile infection: a review of current and emerging therapies. Ann Gastroenterol. 2016;29:147–54.
Bagdasarian N, Rao K, Malani PN. Diagnosis and treatment of Clostridium difficile in adults: a systematic review. JAMA. 2015;313:398–408.
Healthcare-associated infections: Clostridium difficile infections. (2018).
Bauer MP, et al. Clostridium difficile infection in Europe: a hospital-based survey. Lancet. 2011;377:63–73.
Debast SB, Bauer MP, Kuijper E. J & European Society of Clinical Microbiology and Infectious Diseases European Society of Clinical Microbiology and Infectious Diseases: update of the treatment guidance document for Clostridium difficile infection. Clin Microbiol Infect. 2014;20(Suppl 2):1–26.
Singh T, et al. Updates in treatment of recurrent Clostridium difficile infection. J Clin Med Res. 2019;11:465–71.
Song JH, Kim YS. Recurrent Clostridium difficile infection: risk factors, treatment, and prevention. Gut Liver. 2019;13:16–24.
Aslam S, Hamill RJ, Musher DM. Treatment of Clostridium difficile-associated disease: old therapies and new strategies. Lancet Infect Dis. 2005;5:549–57.
Garey KW, Sethi S, Yadav Y, DuPont HL. Meta-analysis to assess risk factors for recurrent Clostridium difficile infection. J Hosp Infect. 2008;70:298–304.
Cho SM, Lee JJ, Yoon HJ. Clinical risk factors for Clostridium difficile-associated diseases. Braz J Infect Dis. 2012;16:256–61.
Appaneal HJ, Caffrey AR, Beganovic M, Avramovic S, LaPlante KL. Predictors of Clostridioides difficile recurrence across a national cohort of veterans in outpatient, acute, and long-term care settings. Am J Health Syst Pharm. 2019;76:581–90.
Eyre DW, et al. Predictors of first recurrence of Clostridium difficile infection: implications for initial management. Clin Infect Dis. 2012;55:S77–87.
Pepin J, et al. Increasing risk of relapse after treatment of Clostridium difficile colitis in Quebec, Canada. Clin Infect Dis. 2005;40:1591–7.
Deshpande A, et al. Risk factors for recurrent Clostridium difficile infection: a systematic review and meta-analysis. Infect Control Hosp Epidemiol. 2015;36:452–60.
Abou Chakra CN, et al. Factors associated with complications of Clostridium difficile infection in a multicenter prospective cohort. Clin Infect Dis. 2015;61:1781–8.
Avni T, et al. Clostridioides difficile infection in immunocompromised hospitalized patients is associated with a high recurrence rate. Int J Infect Dis. 2020;90:237–42.
Revolinski SL, Munoz-Price LS. Clostridium difficile in Immunocompromised hosts: a review of epidemiology, risk factors, treatment, and prevention. Clin Infect Dis. 2019;68:2144–53.
Gaulton T, et al. Complete genome sequence of the Hypervirulent bacterium Clostridium difficile strain G46, Ribotype 027. Genome Announc. 2015;3:e00073–15.
Brouwer MSM, Allan E, Mullany P, Roberts AP. Draft genome sequence of the nontoxigenic Clostridium difficile strain CD37. J Bacteriol. 2012;194:2125–6.
Stabler RA, et al. Comparative genome and phenotypic analysis of Clostridium difficile 027 strains provides insight into the evolution of a hypervirulent bacterium. Genome Biol. 2009;10:R102.
He M, et al. Evolutionary dynamics of Clostridium difficile over short and long time scales. PNAS. 2010;107:7527–32.
Knight DR, Elliott B, Chang BJ, Perkins TT, Riley TV. Diversity and evolution in the genome of Clostridium difficile. Clin Microbiol Rev. 2015;28:721–41.
Kansau I, et al. Deciphering adaptation strategies of the epidemic Clostridium difficile 027 strain during infection through in vivo transcriptional analysis. PLoS One. 2016;11:e0158204.
Sebaihia M, et al. The multidrug-resistant human pathogen Clostridium difficile has a highly mobile, mosaic genome. Nat Genet. 2006;38:779–86.
Stabler RA, et al. In-depth genetic analysis of Clostridium difficile PCR-ribotype 027 strains reveals high genome fluidity including point mutations and inversions. Gut Microbes. 2010;1:269–76.
Didelot X, et al. Microevolutionary analysis of Clostridium difficile genomes to investigate transmission. Genome Biol. 2012;13:R118.
Murillo T, et al. Two groups of Cocirculating, Epidemic Clostridiodes difficile Strains Microdiversify through Different Mechanisms. Genome Biol Evol. 2018;10:982–98.
Scaria J, et al. Analysis of ultra low genome conservation in Clostridium difficile. PLoS One. 2010;5:e15147.
Kumar N, et al. Adaptation of host transmission cycle during Clostridium difficile speciation. Nat Genet. 2019;51:1315–20.
Forgetta V, et al. Fourteen-genome comparison identifies DNA markers for severe-disease-associated strains of Clostridium difficile. J Clin Microbiol. 2011;49:2230–8.
Janvilisri T, et al. Microarray identification of Clostridium difficile Core components and divergent regions associated with host origin. J Bacteriol. 2009;191:3881–91.
Hiller NL, et al. Comparative genomic analyses of seventeen Streptococcus pneumoniae strains: insights into the pneumococcal supragenome. J Bacteriol. 2007;189:8186–95.
Dembek M, et al. High-Throughput Analysis of Gene Essentiality and Sporulation in Clostridium difficile. mBio. 2015;6:02383–14.
Poulsen BE, et al. Defining the core essential genome of Pseudomonas aeruginosa. PNAS. 2019;116:10072–80.
Willcocks SJ, Stabler RA, Atkins HS, Oyston PF, Wren BW. High-throughput analysis of Yersinia pseudotuberculosis gene essentiality in optimised in vitro conditions, and implications for the speciation of Yersinia pestis. BMC Microbiol. 2018;18:46.
Stickland LH. Studies in the metabolism of the strict anaerobes (genus Clostridium). Biochem J. 1934;28:1746–59.
Neumann-Schaal M, Jahn D, Schmidt-Hohagen K. Metabolism the Difficile way: the key to the success of the pathogen Clostridioides difficile. Front Microbiol. 2019;10:219.
Köpke M, Straub M, Dürre P. Clostridium difficile is an autotrophic bacterial pathogen. PLoS One. 2013;8:e62157.
Hofmann JD, et al. Metabolic reprogramming of Clostridioides difficile during the stationary phase with the induction of toxin production. Front Microbiol. 2018;9:1970.
Murima P, McKinney JD, Pethe K. Targeting Bacterial Central Metabolism for Drug Development. Chem Biol. 2014;21:1423–32.
Haag NL, Velk KK, Wu C. Potential antibacterial targets in bacterial central metabolism. Int J Adv Life Sci. 2012;4:21–32.
Katiyar A, Singh H, Azad KK. Identification of missing carbon fixation enzymes as potential drug targets in mycobacterium tuberculosis. J Integr Bioinform. 2018;15:20170041.
Yang F, et al. Genome dynamics and diversity of Shigella species, the etiologic agents of bacillary dysentery. Nucleic Acids Res. 2005;33:6445–58.
Kaper JB, Nataro JP, Mobley HLT. Pathogenic Escherichia coli. Nat Rev Microbiol. 2004;2:123–40.
Womble DD, Rownd RH. Genetic and physical map of plasmid NR1: comparison with other IncFII antibiotic resistance plasmids. Microbiol Rev. 1988;52:433–51.
Lim JY, Yoon JW, Hovde CJ. A brief overview of Escherichia coli O157:H7 and its plasmid O157. J Microbiol Biotechnol. 2010;20:5–14.
Purdy D, et al. Conjugative transfer of clostridial shuttle vectors from Escherichia coli to Clostridium difficile through circumvention of the restriction barrier. Mol Microbiol. 2002;46:439–52.
Schmidt H, Kernbach C, Karch H. Analysis of the EHEC hly operon and its location in the physical map of the large plasmid of enterohaemorrhagic Escherichia coli O157:H7. Microbiology. 1996;142:907–14.
Amy J, et al. Identification of large cryptic plasmids in Clostridioides (Clostridium) difficile. Plasmid. 2018;96–97:25–38.
Clabots CR, Peterson LR, Gerding DN. Characterization of a nosocomial Clostridium difficile outbreak by using plasmid profile typing and clindamycin susceptibility testing. J Infect Dis. 1988;158:731–6.
Steinberg JP, Beckerdite ME, Westenfelder GO. Plasmid profiles of Clostridium difficile isolates from patients with antibiotic-associated colitis in two community hospitals. J Infect Dis. 1987;156:1036–8.
Hornung BVH, Kuijper EJ, Smits WK. An in silico survey of Clostridioides difficile extrachromosomal elements. Microb Genom. 2019;5:e000296.
Boekhoud IM, et al. Plasmid-mediated metronidazole resistance in Clostridioides difficile. Nat Commun. 2020;11:598.
Hargreaves KR, Colvin HV, Patel KV, Clokie JJP, Clokie MRJ. Genetically diverse Clostridium difficile strains harboring abundant Prophages in an estuarine environment. Appl Environ Microbiol. 2013;79:6236–43.
Sekulovic O, Meessen-Pinard M, Fortier L-C. Prophage-stimulated toxin production in Clostridium difficile NAP1/027 Lysogens. J Bacteriol. 2011;193:2726–34.
Meessen-Pinard M, Sekulovic O, Fortier L-C. Evidence of in vivo Prophage induction during Clostridium difficile infection. Appl Environ Microbiol. 2012;78:7662–70.
Fortier L-C. Bacteriophages contribute to shaping Clostridioides (Clostridium) difficile species. Front Microbiol. 2018;9:2033.
Govind R, Vediyappan G, Rolfe RD, Dupuy B, Fralick JA. Bacteriophage-mediated toxin gene regulation in Clostridium difficile. J Virol. 2009;83:12037–45.
Riedel T, et al. A Clostridioides difficile bacteriophage genome encodes functional binary toxin-associated genes. J Biotechnol. 2017;250:23–8.
Goh S, Chang BJ, Riley TV. Effect of phage infection on toxin production by Clostridium difficile. J Med Microbiol. 2005;54:129–35.
Rocha EPC, et al. Comparisons of dN/dS are time dependent for closely related bacterial genomes. J Theor Biol. 2006;239:226–35.
Carter GP, Rood JI, Lyras D. The role of toxin a and toxin B in Clostridium difficile-associated disease. Gut Microbes. 2010;1:58–64.
Carter GP, et al. Defining the Roles of TcdA and TcdB in Localized Gastrointestinal Disease, Systemic Organ Damage, and the Host Response during Clostridium difficile Infections. mBio. 2015;6:e00551.
Rupnik M. Heterogeneity of large clostridial toxins: importance of Clostridium difficile toxinotypes. FEMS Microbiol Rev. 2008;32:541–55.
Rupnik M, Janezic S. An update on Clostridium difficile Toxinotyping. J Clin Microbiol. 2016;54:13–8.
Phothichaisri W, et al. Characterization of bacteriophages infecting clinical isolates of Clostridium difficile. Front Microbiol. 2018;9:1701.
Waker E, et al. High prevalence of genetically related Clostridium Difficile strains at a single Hemato-oncology Ward over 10 years. Front Microbiol. 2020;11:1618.
Chevreux B, et al. Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 2004;14:1147–59.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.
Page AJ, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31:3691–3.
Huerta-Cepas J, et al. Fast genome-wide functional annotation through Orthology assignment by eggNOG-mapper. Mol Biol Evol. 2017;34:2115–22.
Huerta-Cepas J, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47:D309–14.
Kanehisa M, et al. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42:D199–205.
Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: springer; 2009.
Li J, et al. VRprofile: gene-cluster-detection-based profiling of virulence and antibiotic resistance traits encoded within genome sequences of pathogenic bacteria. Brief Bioinform. 2018;19:566–74.
Roosaare M, Puustusmaa M, Möls M, Vaher M, Remm M. PlasmidSeeker: identification of known plasmids from bacterial whole genome sequencing reads. PeerJ. 2018;6:e4588.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Pagès H, Aboyoun P, Gentleman R, DebRoy S. Biostrings: Efficient manipulation of biological strings. Biostrings. 2017. https://doi.org/10.18129/B9.bioc.
Bodenhofer U, Bonatesta E. Horejš-Kainrath, C. & Hochreiter, S. msa: an R package for multiple sequence alignment. Bioinformatics. 2015;31:3997–9.
Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Suyama M, Torrents D, Bork P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006;34:W609–12.
Bruen TC, Philippe H, Bryant D. A simple and robust statistical test for detecting the presence of recombination. Genetics. 2006;172:2665–81.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B (Methodological). 1995;57:289–300.
This study was supported by National Science Centre (Narodowe Centrum Nauki) [grant number the 2017/27/B/NZ5/01504, awarded to JO]. The study sponsor had no role in study design; in the collection, analysis and interpretation of data; in the writing of the report; and in the decision to submit the article for publication.
Ethics approval and consent to participate
The studies involving human participants were reviewed and approved by Bioethics Committee at Maria Skłodowska-Curie National Research Institute of Oncology. In line with the opinion of the Bioethics Committee at Maria Skłodowska-Curie National Research Institute of Oncology our study did not require informed consent for the following reasons: This is a retrospective study describing the genetic differences between C. difficile strains but not between patients; bacterial strains were isolated during routine diagnostics and then banked over the course of one to 10 years; most of these patients are already dead.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Further data are available as Supplementary Tables: Table S1. COG categories prevalence in C. difficile core genome. Table S2. KEGG pathways present in core genome. Table S3. Logistic regression results for odds of infection recurrence after adjustment for ST. Table S4. Gene set enrichment analysis for results of logistic regression. Table S5. Log-likelihood ratio test results for comparison between M2a and M1a models. Table S6. Fisher’s exact test results for prevalence difference in prophage sequence between single and recurrent infections.
About this article
Cite this article
Kulecka, M., Waker, E., Ambrozkiewicz, F. et al. Higher genome variability within metabolism genes associates with recurrent Clostridium difficile infection. BMC Microbiol 21, 36 (2021). https://doi.org/10.1186/s12866-021-02090-9
- Clostridium difficile
- Whole genome sequencing