- Research article
- Open Access
Comparative genomic analysis of toxin-negative strains of Clostridium difficile from humans and animals with symptoms of gastrointestinal disease
© Roy Chowdhury et al. 2016
Received: 29 October 2015
Accepted: 2 March 2016
Published: 12 March 2016
Clostridium difficile infections (CDI) are a significant health problem to humans and food animals. Clostridial toxins ToxA and ToxB encoded by genes tcdA and tcdB are located on a pathogenicity locus known as the PaLoc and are the major virulence factors of C. difficile. While toxin-negative strains of C. difficile are often isolated from faeces of animals and patients suffering from CDI, they are not considered to play a role in disease. Toxin-negative strains of C. difficile have been used successfully to treat recurring CDI but their propensity to acquire the PaLoc via lateral gene transfer and express clinically relevant levels of toxins has reinforced the need to characterise them genetically. In addition, further studies that examine the pathogenic potential of toxin-negative strains of C. difficile and the frequency by which toxin-negative strains may acquire the PaLoc are needed.
We undertook a comparative genomic analysis of five Australian toxin-negative isolates of C. difficile that lack tcdA, tcdB and both binary toxin genes cdtA and cdtB that were recovered from humans and farm animals with symptoms of gastrointestinal disease. Our analyses show that the five C. difficile isolates cluster closely with virulent toxigenic strains of C. difficile belonging to the same sequence type (ST) and have virulence gene profiles akin to those in toxigenic strains. Furthermore, phage acquisition appears to have played a key role in the evolution of C. difficile.
Our results are consistent with the C. difficile global population structure comprising six clades each containing both toxin-positive and toxin-negative strains. Our data also suggests that toxin-negative strains of C. difficile encode a repertoire of putative virulence factors that are similar to those found in toxigenic strains of C. difficile, raising the possibility that acquisition of PaLoc by toxin-negative strains poses a threat to human health. Studies in appropriate animal models are needed to examine the pathogenic potential of toxin-negative strains of C. difficile and to determine the frequency by which toxin-negative strains may acquire the PaLoc.
Clostridium difficile is a Gram-positive pathogen that has emerged to become one of the leading causes of infectious diarrhoea in adult humans, securing its inclusion in the ESCAPE group of pathogens [1–4]. C. difficile infections range from being asymptomatic to causing mild or severe diarrhoea and occasionally life-threatening conditions such as pseudomembranous colitis and toxic megacolon [1, 5]. However, community-acquired C. difficile infection is being reported with increasing frequency  and C. difficile is also emerging as a pathogen in animals particularly cattle, pigs and horses [5, 7–10]. Molecular epidemiological studies show that infections in humans and animals can share the same ribotype or multilocus sequence type (ST)  suggesting that pathogenic C. difficile may traffic between humans and animals, although further studies are needed to confirm these linkages.
C. difficile is a genetically diverse and globally dispersed species [11–16] having a clonal structure comprising six major clades (clades 1, 2, 3, 4, 5 and C-I). Clade C-I is the most phylogenetically divergent clade and may represent of a new subspecies of C. difficile . Clade C-I typically comprise toxin-negative strains of C. difficile  but toxigenic variants that reside in Clade C-I have recently been described . Representatives from most clades have been associated with CDI in humans and comprise toxigenic strains with A+/B+, A−/B+ toxin types [11, 14, 17, 19–22]. Non-toxigenic strains of C. difficile are represented in all six clades .
Toxin expression is considered mandatory for the development of C. difficile disease [23, 24]. Two large clostridial toxins known as toxins A (308 kDa) and B (260 kDa) encoded by tcdA and tcdB and the genes implicated in regulating their expression (tcdC, tcdE and tcdR) reside on a 19.6-kb pathogenicity locus known as the PaLoc [25, 26]. The PaLoc is replaced by 115/75 base pair non-coding region in toxin negative strains of C. difficile . Approximately 20 % of C. difficile strains express a third toxin, known as the binary toxin (CDT) . Genes encoding binary toxin (cdtA and cdtB) and a regulator gene (cdtR) are usually located on a locus (CdtLoc) that is physically separated from the PaLoc. A recent study described six toxin-negative (A-/B-) isolates of C. difficile that were positive for CDT from patient with symptoms of CDI .
Assays that detect toxin genes or the products of their expression dominate laboratory-based tests used to diagnose infections caused by C. difficile [29, 30]. Diagnostic tests that target tcd genes underestimate the frequency of detection of toxin-negative strains (including those that express binary toxin) in C. difficile disease and as such, their role in disease is poorly understood. Phylogenetic studies show that toxin-negative strains of C. difficile cluster tightly with toxin-positive isolates belonging to the same ST  suggesting that presence and absence of the PaLoc may be one of the major defining features that differentiate toxin-negative strains from toxin producing strains of C. difficile. Notably, oral bacteriotherapy with toxin-negative strains or their spores has been used successfully to treat patients undergoing long-term antibiotic regimes and prevent colonisation by toxigenic strains of C. difficile [31–33]. The utility of this therapeutic strategy is supported by previous studies in hamsters which showed that exposure of the gastrointestinal tract to toxin-negative C. difficile strains prevented colonisation by toxin-positive strains [34, 35]. Interestingly, challenge studies in hamsters have shown that toxin-negative strains can effectively colonise the gut [36, 37] suggesting that toxin production may be of little consequence in determining the success of colonisation of the gastrointestinal tract. Notably, the toxin-negative strain CD1342 (tcdA − , tcdB − , cdtA − and cdtB − ) was reported to elicit an innate immune response in the caecum resulting in neutrophil infiltration, damage to epithelial mucosa and localised haemorrhagic congestion . These findings suggest that virulence factors are carried by C. difficile in addition to the known toxins that can induce host pathology.
Studies of toxin-negative C. difficile strains have focused on the characterisation of functional binary toxins and their roles in pathogenesis [28, 33, 38]. The binary toxins cdtA and cdtB have adenosine diphosphate ribosyltransferase activity but their capacity to induce symptoms of C. difficile infection remains unclear [39–42]. Several adhesins, ECM-binding proteins, proteases, motility proteins, hydrolytic enzymes and other surface-associated proteins have been described in C. difficile and these factors are likely to contribute significantly to the establishment, progression and severity of C. difficile disease [11, 43]. Therefore, further studies are needed to examine the pathogenic potential of toxin-negative strains of C. difficile and to determine the frequency at which toxin-negative strains may acquire the PaLoc and express toxins.
Studies that seek to understand the evolutionary history of the PaLoc highlight the complex nature of the multiple clade-specific acquisitions that have occurred after clonal expansion of each clade in populations of C. difficile . Those studies report homologous and site-specific recombination events as having played an important role in the loss and gain of the PaLoc . The PaLoc is proposed to be a mobile element that can transfer to toxin-negative strains rendering the recipient with the ability to produce clinically relevant concentrations of ToxA and ToxB . Toxin-negative strains are purported to be ancestral to modern C. difficile but lateral genetic events complicate phylogenetic interpretation and alternate hypotheses have been proposed . Genomic studies incorporating a greater diversity of toxin-negative strains of C. difficile are needed to shed light on their potential to cause disease.
Isolation and culture of Clostridium difficile
All C. difficile isolates analysed in this study (P29, 5.3, 19.3, 22.1, H3) were obtained from watery diarrhoea stool samples from their respective hosts (Additional file 1: Table S1). The porcine and equine C. difficile isolates analysed in this study were sourced in 2008 from different geographical locations in New South Wales, Australia. The porcine isolate P29 was isolated from a stool sample submitted by the veterinarian attending a piglet with severe but non-fatal diarrhoea. The equine isolate H3 was isolated from a live neonatal foal suffering from non-fatal watery diarrhoea. Stool samples were tested with PCR targeting major ETEC virulence genes  and common viruses known to cause diarrhoea in neonatal animals and were plated on blood agar plates to select for enteric pathogens. The stool specimens were initially tested for Escherichia coli, Clostridium perfringens and C. difficile using species-specific PCR primers . Briefly, DNA was extracted from 500 μl of stool sample using a FastDNA spin kit (QBiogene, California, USA) and used as a template for PCR using primers specific for C. difficile and C. perfringens 16S rDNA , tcdA and tcdB genes (see below) and for E. coli . To enrich for C. difficile 100 μl of each faecal sample was added to 10 ml cooked meat medium (TM0102 Oxoid Australia) and incubated anaerobically at 37 °C for 24 h using the anoxomat system (MART Microbiology B.B., The Netherlands).
Two hundred μl of culture samples that tested positive for C. difficile by PCR were transferred (from cooked meat media enrichment broth) into an Eppendorf tube and centrifuged (10,000 rpm, 5 min). The pellet was resuspended in 1 ml of absolute ethanol (room temperature, 2 h with periodic inversions), harvested by centrifugation (10,000 rpm, 5 min), resuspended in brain heart infusion broth (100 μl) and plated onto C. difficile selective agar (CC-BHIA + Taurocholate, PP2362 Oxoid Australia). Plates were incubated under anaerobic conditions at 37 °C for 24 h. Colonies morphologically representing C. difficile from each plate were selected and sub-cultured onto CC-BHIA + Taurocholate until pure cultures were achieved.
For routine PCR, template DNA was extracted with Chelex (BIO-RAD) from 2 ml brain heart infusion broth cultures grown under anaerobic conditions at 37 °C for 48 h. Briefly, cell pellets were obtained by centrifuging (10,000 rpm for 5 min) 200 μl aliquots of liquid culture, washed 2 × with 500 μl of sterile water and resuspended in 200 μl of 6 % Chelex solution made in Tris-EDTA buffer (pH 7.5). The samples were incubated at 56 °C for 20 min, vortexed for 10 s and incubated at 100 °C for 8 min. After incubation, the sample was immediately transferred to ice. One aliquot was stored at 4 °C for routine PCR tests while the other aliquots were archived at −20 °C.
Sequencing-quality genomic DNA was prepared from 2 ml brain heart infusion broth culture of isolates grown under anaerobic conditions at 37 °C for 48 h. The overnight culture was harvested by centrifugation (10,000 rpm for 10 min), washed in sterile PBS and resuspended in 180 μl of lysis buffer comprising 20 mM Tris-HCl, pH 8.0, 2 mM EDTA, 1.2 % Triton X-100 and lysozyme (20 mg ml−1) and incubated for 45 min at 37 °C. DNA was isolated using a DNeasy® Blood and Tissue Kit (Qiagen) by adhering to the manufacturer’s instructions for the extraction of DNA from Gram-positive bacteria.
C. difficile specific 16S rDNA primers, C.diff-F: 5′-TTGAGCGATTTACTTCGGTAAAGA-3′ and C.diff-R: 5′-CCATCCTGTACTGGCTCACCT-3′ were used for identification and confirmation of C. difficile in enrichment as well as pure cultures. The presence of the tpi gene (encoding Triose Phosphate Isomerase), tcdA gene (encoding Toxin A) and tcdB gene (encoding Toxin B) were tested using previously published primer pairs. Conditions for PCR were as described previously  with minor modifications. Briefly, PCR was carried out in 25 μl volumes containing 2 μl of Chelex extracted DNA, 2.5 μl of 10 × PCR buffer, 1.5 mM of MgCl2, 1 mM of each dATP, dGTP, dCTP and dDTP (Bioline, Australia), 0.5 μM of each primer and 1 U of BioRad Taq polymerase (Bioline, Australia). PCR cycling conditions consisted of an initial denaturation cycle (2 min, 95 °C) followed by 30 cycles of denaturation (94 °C, 1 min), annealing (55 °C, 1 min) and extension (72 °C, 2 min). The cycling process was completed with a final extension of 72 °C for 5 min.
Whole genome sequencing, data assembly and phylogenetic analysis
Sequencing was performed at the Next Generation Sequencing facility within the ithree institute at the University of Technology Sydney using a bench top Illumina MiSeq® sequencer and MiSeq V3 chemistry. Genomic DNA stocks shipped to the sequencing facility at concentrations between 1.8 and of 3.7 ng μl−1 were used as template for the preparation of sequencing libraries. The genomes were sequenced and assembled de novo using published protocols . Raw data and assembled genome sequences were submitted in GenBank under the following Bio-project numbers, 5.3: PRJNA232267, 19.3: PRJNA239262, 22.1: PRJNA239264, P29: PRJNA239265 and H3: PRJNA238844.
PhyloSift was used to conduct a phylogenetic analysis of the five C. difficile genomes (P29, 5.3, 19.3, 22.1, H3) with nine closed C. difficile genomes including strains M120, CF5, M68, 2007855, BI1, CD196, R20291, ATCC43255 and CD630 available in the NCBI genome database on the 18th of December 2014 . FigTree version 1.4.0 (http://tree.bio.ed.ac.uk/software/figtree/) was used to draw phylogenetic trees. Genome sequences of C. perfringens (ATCC13124), Clostridium botulinum (ATCC19397) and Clostridium tetani (E88) were included as outgroups in the analysis. To improve visual resolution of the evolutionary distances between test and reference strains of C. difficile the final figure was generated without the out-groups.
For reference-genome based phylogenetic inference, raw Illumina reads from all taxa were mapped to a single reference (strain CD630) using BWA-MEM (ver0.7.9a) (Li unpublished, github commit: 3efc33160c) and consensus sequences generated using the samtools/bcftools (ver0.1.19-96b5f2294a) tool-chain . The complete set of consensus sequences were combined into a multiple sequence alignment. 1,216,986 alignment columns containing unresolved nucleotides (N) were removed using Mothur (ver1.33.3) . A total of 3,073,266 (72 %) polymorphic and non-polymorphic sites were retained for further analysis. The inclusion of invariant sites has been demonstrated to improve accuracy of whole genome phylogeny . Maximum likelihood phylogenetic inference was employed using RAxML (ver8.0.20)  with the following options: raxmlHPC-PTHREADS-SSE3 -T 40 -f a -x 2136841 -p 1486312 -N autoMRE -m GTRCAT. Inference was carried out under a general time reversible (GTR) substitution model with an infinite mixture model for substitutional heterogeneity (CAT), following the suggestion of the RAxML user guide for datasets of this size. The CAT approximation has been previously demonstrated to be an accurate and highly efficient alternative to Gamma-distributed rate heterogeneity on data sets with many taxa (73 – 1663) [53, 54]. Confidence in each clade of the Maximum Likelihood tree was estimated using the rapid bootstrap procedure  with automatic extended majority-rule criterion (100 bootstraps) and the resulting tree and bootstrap confidence estimates were visualized with FigTree version 1.4.0 (http://tree.bio.ed.ac.uk/software/figtree/).
Multi locus sequence typing and comparative genomic analysis
The online version of C. difficile PubMLST database (http://pubmlst.org/cdifficile/) was used to sequence type the isolates from the assembled genome sequences. The database was also exploited to locate certain genes of interest.
The online version of the RAST annotation server (http://rast.nmpdr.org/)  was used to annotate the genomes. The Classic RAST annotation scheme and FigFAM release 70 were used to predict genes (5.3 = RAST-ID 6666666.71923, 19.3 = RAST-ID 6666666.71924, H3 = RAST-ID 6666666.72094, P29 = RAST-ID 1440056.4 and 22.1 = RAST-ID 6666666.72093). Amino acid sequences corresponding to translated peptide products of all open reading frames predicted by RAST  from each of the five genomes were used in the ‘all Vs all’ homology search protocol deposited in the github repository as cRBLH (https://github.com/cerebis/crblh/tree/v0.1). The protocol included clustering of predicted peptide sequences using a modified reciprocal best hit method, where simplicity was favoured for the apparent advantage in identifying orthogroups . The all vs. all homology search was carried out with LAST [59–61] using runtime parameters (−T 1 -f 0 -e 100). The best hits were used to generate a directed graph with genes as vertices and best hits as edges. Unidirectional links between any two nodes were then pruned. Sets of disconnected subgraphs were then analysed for weak intra-cluster linkages, which likely represented overlap between partially homologous protein clusters. Each subgraph was subjected to modularity optimisation  and further decomposed until modularity scores of constituent elements fell below a given threshold (0.2). The nodes of the resulting subgraphs were then written out as protein clusters. Singletons defined as nodes without a single edge to any other were deemed unique/isolated genes.
Whole genome comparisons were performed using Mauve version 2.3.1 [63, 64] and iterative BLASTn analysis (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Inter-isolate regions of interest identified from genome-wide comparisons using Mauve, BLASTp and protein clustering analyses were analysed further using iterative BLASTn and BLASTp searches. Figures of comparative genomic analysis, including comparisons of the PaLoc, were compiled using locally downloaded version of EasyFig version 2.1 .
The genome of the epidemic C. difficile CD630 strain was used in whole genome BLASTp analysis (in RAST) with our test C. difficile genomes to identify genes that have been correlated with pathogenicity. All genes deemed as candidate alternative virulence genes or genes for which the products could potentially confer pathogenic traits were individually interrogated using BLASTp and setting amino acid alignment cut off set to 100 % of input query sequence to avoid any data extrapolation.
Toxin-negative C. difficile from animals and humans with clinical disease
The original stool samples and primary enrichment cultures of the stool samples tested negative for C. difficile toxins A and B. PCR assays using DNA from enrichment broths tested negative for enteric (other than C. difficile) and viral pathogens associated with neonatal diarrhoea. The porcine faecal sample was negative for the enterotoxigenic E. coli genes STa, STb and LT and C. perfringens and the disease symptomology did not correlate with viral disease as diagnosed by the attending veterinarian. Similarly, the foal sample was tested for E. coli, Salmonella enterica and rotavirus and none were detected. While gastrointestinal disease was most likely associated with the presence of toxin-negative C. difficile we cannot rule out the possibility that disease was caused by unculturable/unknown pathogens present in the gastrointestinal tract of these animals. Toxin-negative human C. difficile isolates 5.3, 19.3 and 22.1 were collected in the course of routine diagnostic tests for C. difficile-associated diarrhoea in patients presenting typical symptoms of the disease at a gastrointestinal clinic in Sydney, Australia in 2008.
Interrogation of the C. difficile PubMLST database confirmed that none of the toxin-negative isolates in our cohort (P29, H3, 5.3, 19.3, 22.1) had homologs of the known C. difficile PubMLST toxin genes (tcdA, tcdB, cdtA and cdtB) confirming our initial diagnostic PCR data for toxin A and B genes (Additional file 1: Table S1). Isolates 5.3 (ST15), P29 (ST109) and H3 (ST29) were distinct from each other and from ST types of Australian isolates included in a recent phylogenetic study of C. difficile (Additional file 1: Table S1) .
Phylogenetic analysis of toxin-negative isolates of C. difficile
A study of the evolution of the C difficile pathogenicity locus (PaLoc) identified an extremely divergent clade C-I that exclusively comprised toxin-negative isolates predominantly of Australian origin . A maximum-likelihood phylogenetic tree using a reference-based, whole genome alignment protocol (see methods section for protocol details) that incorporates both variant and invariant sites of the C. difficile genome sequences was used to verify the ancestry of our toxin-negative isolates. Our approach uses approximately 72 % of the C. difficile genome for the analysis, considerably more than what was used in the original study . All 73 genomes and the reference genome CD630 used in the previous study  as well as additional closed genomes of C. difficile (strains 2007855, ATCC43255, BI1, CF5, M68, M120 and R20291 from the GenBank database) were used in our initial phylogenetic analysis. A preliminary phylogenetic tree (Additional file 2: Figure S12) revealed that our genome-based phylogeny was largely congruent with that described in an earlier study  and clearly indicated that none of the five genomes that were the subject of our study clustered within the divergent clade C-I.
Our genome sequences were assembled with a de novo assembler using A5 . Prior to conducting a detailed analysis of the toxin-negative C. difficile isolates we identified the closest reference genome for tiling genomic scaffolds. A preliminary phylogeny generated using PhyloSift and FastTree (Additional file 2: Figure S2) indicated that C. difficile strain CF5 (toxin-positive ST86) was the most appropriate reference to order genomic scaffolds of isolates 19.3 (ST39), 22.1 (ST39) and P29 (ST109) while C. difficile strain 630 (toxin-positive, ST54 (PCR ribotype 012) was appropriate to order genomic scaffolds of isolates 5.3 (ST15) and H3 (ST29). Strain CF5 was isolated from a patient in Belgium in 1995 while CD630 is a highly virulent, multiple antibiotic resistant strain of C. difficile that caused pseudomembranous colitis in a human patient and later caused an epidemic of C. difficile infection in a Swiss hospital ward in 1982. All down-stream analyses of the genomes presented in this study were performed on genomic assemblies with scaffolds ordered to match the reference genomes.
Homology based functional similarity in the toxin-negative isolates
Initially a Progressive Mauve alignment performed (Additional file 2: Figure S3) on genome sequences of human ST39 isolates 19.3 and 22.1 revealed a high level of nucleotide identity across the genomes with 324 SNP differences. Most of the SNPs were clustered into 7 groups (see Additional file 3: Table S4) suggesting that lateral gene transfer or homologous recombination-mediated genomic rearrangements may be responsible for the differences and were not considered further. Only 27 SNPs were identified that could generate changes in the amino acid sequence of the predicted proteins in the table Additional file 3: Table S4.
To identify genes affected by the SNP changes, a bi-directional BLASTp comparison of 3772 proteins comprising the predicted proteome of strain 19.3 was performed with 3764 predicted protein sequences from strain 22.1 in RAST. The analysis identified 3750 protein sequences that were identical in both the genomes. Nine protein sequences had greater than 99 % sequence identity, six others had greater than 97 % sequence identity and one ORF encoding a hypothetical protein showed 49 % sequence identity (Additional file 3: Table S5). Proteins sharing 97 and 99 % sequence identity predominantly encoded components of the bacterial cell surface including N-acetylmuramoyl-L-alanine amidase, flagellar assembly protein FliH, lipoprotein signal peptidase, putative ABC transporters and permeases. Eight ORFs in isolate 19.3 predominantly encoding hypothetical proteins were missing in isolate 22.1. The genomes of isolates 19.3 and 22.1 comprise 4,181,809 and 4,180,898 bp respectively.
Since isolates 22.1 and 19.3 had high levels of homology, predicted proteins only from isolate 19.3 were included in a pairwise bi-directional BLASTp analysis that seeks to identify conserved genes and major differences among toxin-negative isolates of Australian origin. Isolate P29 had the largest predicted proteome in our collection and was used as the reference. Comparisons of predicted proteomes of human isolates 5.3 & 19.3 and equine isolate H3 with P29 identified major differences in regions harbouring prophage-associated proteins, hypothetical proteins (Additional file 3: Table S5) and putative transposases associated with mobile genetic elements.
Chromosomal context of the PaLoc insertion site
Comparative BLASTp analyses of toxin-negative strains to identify functional similarities between groups of isolates
Summary of protein clustering results within the different sub-clades containing the five toxin negative isolates included in this study. Summary of protein clusters within the different sub-clades of C. difficile
No of predicted proteins
Total number of unique peptides in:
C. difficile C0000509
C. difficile C0000541
C. difficile C0000562
C. difficile C00005945
C. difficile H3
Total number of unique peptides in:
C. difficile 19.3
C. difficile 22.1
C. difficile C00000089
C. difficile C000011286
C. difficile C00006473
C. difficile C00007671
C. difficile C00011267
C. difficile CF5
C. difficile M68
C. difficile P29
Ten clade 4 genomes boxed in Fig. 1 shared 3357 proteins (Table 1). Isolate P29 had the highest number (299) of unshared/unique proteins within the clade 4 cohort and most of the 299 proteins were phage-related (Additional file 3: Table S9). Some of the unique proteins clustered together in the same scaffold indicating lateral movement of phage-associated genomic DNA. C. difficile isolates 19.3 and 22.1 carried eight and seven unique proteins respectively. The handful of unshared proteins in 19.3 and 22.1 were attributed to mobile genetic elements or were designated to encode proteins of unknown function. The five C. difficile strains within clade 1 shared 3323 proteins. Equine isolate H3 had 111 unique proteins, most of which were phage related or hypothetical with some clustered in single scaffolds (Additional file 3: Table S9).
Summary of Phage related regions identified by PHAST in the 5 genomes. Phage sequences identified in this study
PHAST region identifier
Length of prophage
No of predicted CDS
Relative position on genome
Location on Genomic scaffolds
C difficile P29 genome
27.1, 30.1, 40.1, 34.1, 36,1, 5.1, 47.1 and 16.1
5.1, 19.1 and 22.1
8.1, 26.1, 41.1 and 42.1
In over 20 very small scaffolds
C. difficile H3 genome
22.1 and 31.1
18.1 and 5.1
26.1, 27.1, 28.1, 32.1, 33.1
C. difficile 5.3 genome
C.difficle 19.3 genome
C. difficile 22.1 genome:
Comparative BLASTp analysis of isolates 19.3, 5.3, H3 and P29 identified a 119.3 kb region on contig 11 in P29 (Additional file 5: Table S6). An all versus all protein clustering analysis also identified a subset of unique proteins on contig 11 of the P29 genome but not in the 10 human C. difficile genomes within clade 4 (Fig. 1). Equine isolate H3 was not included in this analysis as it was on a different clade. BLASTn analysis of the 119.3 kb region against the C. difficile genome database in GenBank identified similarity at the DNA level to parts of C. difficile strain CD630 indicating a phage-mediated lateral movement of parts of the genome of CD630 into the genomes of isolates H3 and P29.
Homology based functional prediction of Putative C. difficile virulence factors implicated in host colonization
Proteins derived from C. difficile CD630 that are predicted to play a role in pathogenesis
Selected gene and product
C. difficile 630 locus tag
RAST Annotation identifiers
In C. difficile 5.3 (% identity)
In C. difficile 19.3 (% identity)
In C. difficile H3 (% identity)
In Cc difficile P29 (% identity)
Flagellin C gene fliC
Flagellin D gene fliD
Precursor S-layer protein gene slpA
fig|6666666.71924.peg.3725 (59) *
Stage 0 Sporulation gene spoA
Fibrinectin binding proten encoding fbpA gene
GroEL encoding gene groL
Cell surface protein cwp66
Cell wall binding protein encoding cwp2
Cell wall binding protein encoding cwp12
Cell wall binding protein encoding cwp11
Cell wall binding protein encoding cwp9
Cell wall hydrolase (LPXTG)
Cell wall binding protein encoding cwp25 gene
N-acetylmuramoyl-L-analini amidase encoding cwp16
Cell wall hydrolase encoding gene (invasin)
yes, RNAseq and proteome
LmbE-like deacetylase encoding gene
Invasin/Sh3 domain containing surface protein
Cell wall hydrolase/Invasin associated protein
Autolysin acd gene homolog/mannosyl-glycoprotein endo neta N acetylglucosamine
Protease/Serine protease, HrtA family
Intracellular serine protease
Ser-type protease/subtilisin-like serine germination related protease
yes, Mass spectrometry
Serine protease precursor/Subtilinase subfamily
Membrane-associated zinc metalloprotease/M50 family peptidase
Zinc Protease/M16 family peptidase
C. difficile colonisation in humans is age dependent. While asymptomatic carriage is common in infants less than three years of age it is rare in adults . As such, infants can be a major reservoir of both pathogenic and toxin-negative strains in a community setting . We isolated five toxin-negative isolates of C. difficile including three from humans (22.1, 19.3, 5.3), one from a pig (P29) and one (H3) from a horse all showing symptoms of gastrointestinal disease. Despite efforts to identify C. difficile toxin genes or toxin gene products in the stool samples during the course of the isolation of these strains, none were detected. Phylogenetic studies showed that human ST39 isolates 19.3 and 22.1 and porcine ST109 isolate P29 grouped with clinical human toxigenic strains of ST39 and ST109 respectively in Clade 4. Furthermore, human ST15 isolate 5.3 and equine ST29 isolate H3 grouped with human clinical toxigenic strains with ST15 and ST29 respectively in Clade 1. Comparative genome analyses showed that our toxin-negative isolates displayed virulence gene profiles akin to those identified in toxigenic strains. The animals from which samples were collected in this study exhibited gastrointestinal disease and we were unable to attribute these symptoms to the presence of toxin-positive strains of C. difficile. Given the mobility of the PaLoc  and evidence that the acquisition or loss of the PaLoc via recombination  has occurred multiple times during the evolution of the five major clades of C. difficile [11, 17], our data reinforces calls to include toxin-negative strains in genomic epidemiological studies of C. difficile [17, 36] and to better characterise asymptomatic carriage of closely related Clostridia in gut microbiome surveys such as the Human Microbiome Project and MetaHIT, both in humans and close animal contacts.
Our reference-based whole genome alignment and phylogeny analyses support the global population structure of C. difficile as described by Dingle et al. in 2014 [14, 17, 19]. Each clade has been shown previously to have representatives of both toxin-positive and toxin-negative strains . Our toxin-negative isolates (19.3, 22.1, 5.3, P29, H3) belonged to STs that are distinct from those reported in an earlier study . The role of toxins in C. difficile infection has been extensively studied but factors that enable C. difficile to efficiently colonise the human gastrointestinal tract are relatively poorly understood and are not associated with genes encoded on the PaLoc. It is not known why some toxigenic strains evolve into dominant hypervirulent clones. Thus, considering the genetic diversity inherent within the phylogenetic structure of C. difficile  a sub-population of toxin-negative strains of C. difficile that are efficient colonisers of the host gastrointestinal tract may readily acquire the PaLoc and evolve to become future hypervirulent strains. Several proteins have been suggested to play crucial roles in the colonization of gastrointestinal epithelium and disease progression [43, 71, 73, 74, 76–78]. A recent global proteome study of C. difficile strains CD630 and R20291 has identified numerous extracellular proteins from culture supernatants that may contribute to the virulence attributes of these strains .
Our study reinforced the important role played by phage in the evolution of C. difficile. While PHAST analysis was useful for identifying phage sequences, the analysis may not have identified the full extent of lysogenic phage because our draft genomes remain in multiple scaffolds. Although the complete sequence of phage phiC2 was identified in isolates 19.3 and 22.1 the regions that had significant homology with phiC2 in isolates P29 and H3 were located on different scaffolds. We used a scaffold tiling approach against the closed genome of a reference strain to create the input file for PHAST analysis (PHAST converts the scaffolded genomes into a concatenated artificial chromosome prior to predicting the phage content) and as such it remains a possibility that the partial matches are a consequence of the data handling process. Phage phiC2 is present in the majority of human isolates of C. difficile . However, we detected regions of phiC2 in strains P29 and H3 suggesting that further studies are needed to address issues surrounding the association of phiC2 in C. difficile of animal origin. We also identified the C. difficile temperate bacteriophage phiCD6356 from the Siphoviridae family in isolates P29 and H3 but not in our human isolates of C. difficile. Genomes of bacteriophages belonging to the Siphoviridae family range in size from 14 to 50 kb [79, 80] and this broad range may be a reflection of the stringency governing the amount of DNA that can be packaged by phiCD6356. In addition to the acquisition of phage-associated genes, a 119.3-kb region on contig 11 in isolate P29 was also identified in the course of this analysis. This region is unique to the P29 genome and displayed significant DNA sequence identity to portions of the CD630 genome. It remains unknown if the 119.3-kb region exists in C. difficile strains of porcine origin. Further analyses with greater numbers of genomes from both human and animal sources are required to conclusively address these questions.
Our studies reinforce calls to improve our understanding of the physiological conditions that promote lateral transfer of the PaLoc in the gastrointestinal tract . This is important because the conditions that facilitate movement of fragments of DNA carrying the PaLoc and their recombination into the chromosome are also conducive to the movement of conjugative transposons that carry antibiotic resistance genes and putative virulence factors as independent genetic events .
Genome sequences reported in this analysis were submitted to GenBank and are available via the accession numbers provided. The bioinformatics softwares are made available through the GitHub repository links.
This work is a product of the ausgem partnership. The authors wish to acknowledge Prof Thomas Borody for kindly providing the human isolate included in this analysis.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Karadsheh Z, Sule S. Fecal transplantation for the treatment of recurrent clostridium difficile infection. N Am J Med Sci. 2013;5(6):339–43. Pubmed Central PMCID: 3731863.View ArticlePubMedPubMed CentralGoogle Scholar
- Lessa FC, Mu Y, Bamberg WM, Beldavs ZG, Dumyati GK, Dunn JR, et al. Burden of Clostridium difficile infection in the United States. N Engl J Med. 2015;372(9):825–34.View ArticlePubMedGoogle Scholar
- Peterson LR. Bad bugs, no drugs: no ESCAPE revisited. Clin Infect Dis. 2009;49(6):992–3.View ArticlePubMedGoogle Scholar
- Redelings MD, Sorvillo F, Mascola L. Increase in Clostridium difficile-related mortality rates, United States, 1999–2004. Emerg Infect Dis. 2007;13(9):1417–9. Pubmed Central PMCID: 2857309.View ArticlePubMedPubMed CentralGoogle Scholar
- Rupnik M, Wilcox MH, Gerding DN. Clostridium difficile infection: new developments in epidemiology and pathogenesis. Nat Rev Microbiol. 2009;7(7):526–36.View ArticlePubMedGoogle Scholar
- Khanna S, Pardi DS, Aronson SL, Kammer PP, Baddour LM. Outcomes in community-acquired Clostridium difficile infection. Aliment Pharmacol Ther. 2012;35(5):613–8. Pubmed Central PMCID: 3293482.View ArticlePubMedPubMed CentralGoogle Scholar
- Bauer MP, Kuijper EJ. Potential sources of Clostridium difficile in human infection. Infect Dis Clin North Am. 2015;29(1):29–35.View ArticlePubMedGoogle Scholar
- Songer JG, Anderson MA. Clostridium difficile: an important pathogen of food animals. Anaerobe. 2006;12(1):1–4.View ArticlePubMedGoogle Scholar
- Hensgens MP, Keessen EC, Squire MM, Riley TV, Koene MG, de Boer E, et al. Clostridium difficile infection in the community: a zoonotic disease? Clin Microbiol Infect. 2012;18(7):635–45.View ArticlePubMedGoogle Scholar
- Goorhuis A, Debast SB, van Leengoed LA, Harmanus C, Notermans DW, Bergwerff AA, et al. Clostridium difficile PCR ribotype 078: an emerging strain in humans and in pigs? J Clin Microbiol. 2008;46(3):1157. Pubmed Central PMCID: 2268365, author reply 8.View ArticlePubMedPubMed CentralGoogle Scholar
- Dingle KE, Griffiths D, Didelot X, Evans J, Vaughan A, Kachrimanidou M, et al. Clinical Clostridium difficile: clonality and pathogenicity locus diversity. PLoS One. 2011;6(5):e19993. Pubmed Central PMCID: 3098275.View ArticlePubMedPubMed CentralGoogle Scholar
- Walk ST, Micic D, Jain R, Lo ES, Trivedi I, Liu EW, et al. Clostridium difficile ribotype does not predict severe infection. Clin Infect Dis. 2012;55(12):1661–8. Pubmed Central PMCID: 3501335.View ArticlePubMedPubMed CentralGoogle Scholar
- Cairns MD, Stabler RA, Shetty N, Wren BW. The continually evolving Clostridium difficile species. Future Microbiol. 2012;7(8):945–57.View ArticlePubMedGoogle Scholar
- Stabler RA, Dawson LF, Valiente E, Cairns MD, Martin MJ, Donahue EH, et al. Macro and micro diversity of Clostridium difficile isolates from diverse sources and geographical locations. PLoS One. 2012;7(3):e31559. Pubmed Central PMCID: 3292544.View ArticlePubMedPubMed CentralGoogle Scholar
- Behroozian AA, Chludzinski JP, Lo ES, Ewing SA, Waslawski S, Newton DW, et al. Detection of mixed populations of Clostridium difficile from symptomatic patients using capillary-based polymerase chain reaction ribotyping. Infect Control Hosp Epidemiol. 2013;34(9):961–6. Pubmed Central PMCID: 4016961.View ArticlePubMedPubMed CentralGoogle Scholar
- Waslawski S, Lo ES, Ewing SA, Young VB, Aronoff DM, Sharp SE, et al. Clostridium difficile ribotype diversity at six health care institutions in the United States. J Clin Microbiol. 2013;51(6):1938–41. Pubmed Central PMCID: 3716112.View ArticlePubMedPubMed CentralGoogle Scholar
- Dingle KE, Elliott B, Robinson E, Griffiths D, Eyre DW, Stoesser N, et al. Evolutionary history of the Clostridium difficile pathogenicity locus. Genome Biol Evol. 2014;6(1):36–52. Pubmed Central PMCID: 3914685.View ArticlePubMedPubMed CentralGoogle Scholar
- Monot M, Eckert C, Lemire A, Hamiot A, Dubois T, Tessier C, et al. Clostridium difficile: New insights into the evolution of the pathogenicity locus. Sci Rep. 2015;5:15023. Pubmed Central PMCID: 4597214.View ArticlePubMedPubMed CentralGoogle Scholar
- Griffiths D, Fawley W, Kachrimanidou M, Bowden R, Crook DW, Fung R, et al. Multilocus sequence typing of Clostridium difficile. J Clin Microbiol. 2010;48(3):770–8. Pubmed Central PMCID: 2832416.View ArticlePubMedPubMed CentralGoogle Scholar
- Lemee L, Bourgeois I, Ruffin E, Collignon A, Lemeland JF, Pons JL. Multilocus sequence analysis and comparative evolution of virulence-associated genes and housekeeping genes of Clostridium difficile. Microbiology. 2005;151(Pt 10):3171–80.View ArticlePubMedGoogle Scholar
- Lemee L, Dhalluin A, Pestel-Caron M, Lemeland JF, Pons JL. Multilocus sequence typing analysis of human and animal Clostridium difficile isolates of various toxigenic types. J Clin Microbiol. 2004;42(6):2609–17. Pubmed Central PMCID: 427854.View ArticlePubMedPubMed CentralGoogle Scholar
- Lemee L, Dhalluin A, Testelin S, Mattrat MA, Maillard K, Lemeland JF, et al. Multiplex PCR targeting tpi (triose phosphate isomerase), tcdA (Toxin A), and tcdB (Toxin B) genes for toxigenic culture of Clostridium difficile. J Clin Microbiol. 2004;42(12):5710–4. Pubmed Central PMCID: 535266.View ArticlePubMedPubMed CentralGoogle Scholar
- Rupnik M. How to detect Clostridium difficile variant strains in a routine laboratory. Clin Microbiol Infect. 2001;7(8):417–20.View ArticlePubMedGoogle Scholar
- Voth DE, Ballard JD. Clostridium difficile toxins: mechanism of action and role in disease. Clin Microbiol Rev. 2005;18(2):247–63. Pubmed Central PMCID: 1082799.View ArticlePubMedPubMed CentralGoogle Scholar
- Hundsberger T, Braun V, Weidmann M, Leukel P, Sauerborn M, von Eichel-Streiber C. Transcription analysis of the genes tcdA-E of the pathogenicity locus of Clostridium difficile. Eur J Biochem/FEBS. 1997;244(3):735–42.View ArticleGoogle Scholar
- Matamouros S, England P, Dupuy B. Clostridium difficile toxin expression is inhibited by the novel regulator TcdC. Mol Microbiol. 2007;64(5):1274–88.View ArticlePubMedGoogle Scholar
- Braun V, Hundsberger T, Leukel P, Sauerborn M, von Eichel-Streiber C. Definition of the single integration site of the pathogenicity locus in Clostridium difficile. Gene. 1996;181(1–2):29–38.View ArticlePubMedGoogle Scholar
- Eckert C, Emirian A, Le Monnier A, Cathala L, De Montclos H, Goret J, et al. Prevalence and pathogenicity of binary toxin-positive Clostridium difficile strains that do not produce toxins A and B. New Microbes New Infect. 2015;3:12–7. Pubmed Central PMCID: 4337936.View ArticlePubMedPubMed CentralGoogle Scholar
- Collins DA, Elliott B, Riley TV. Molecular methods for detecting and typing of Clostridium difficile. Pathology. 2015;47(3):211–8.View ArticlePubMedGoogle Scholar
- Rupnik M, Brazier JS, Duerden BI, Grabnar M, Stubbs SL. Comparison of toxinotyping and PCR ribotyping of Clostridium difficile strains and description of novel toxinotypes. Microbiology. 2001;147(Pt 2):439–47.View ArticlePubMedGoogle Scholar
- Villano SA, Seiberling M, Tatarowicz W, Monnot-Chase E, Gerding DN. Evaluation of an oral suspension of VP20621, spores of nontoxigenic Clostridium difficile strain M3, in healthy subjects. Antimicrob Agents Chemother. 2012;56(10):5224–9. Pubmed Central PMCID: 3457387.View ArticlePubMedPubMed CentralGoogle Scholar
- Nagaro KJ, Phillips ST, Cheknis AK, Sambol SP, Zukowski WE, Johnson S, et al. Nontoxigenic Clostridium difficile protects hamsters against challenge with historic and epidemic strains of toxigenic BI/NAP1/027 C. difficile. Antimicrob Agents Chemother. 2013;57(11):5266–70. Pubmed Central PMCID: 3811292.View ArticlePubMedPubMed CentralGoogle Scholar
- Natarajan M, Walk ST, Young VB, Aronoff DM. A clinical and epidemiological review of non-toxigenic Clostridium difficile. Anaerobe. 2013;22:1–5. Pubmed Central PMCID: 3729612.View ArticlePubMedPubMed CentralGoogle Scholar
- Seal D, Borriello SP, Barclay F, Welch A, Piper M, Bonnycastle M. Treatment of relapsing Clostridium difficile diarrhoea by administration of a non-toxigenic strain. Eur J Clin Microbiol. 1987;6(1):51–3.View ArticlePubMedGoogle Scholar
- Wilson KH, Sheagren JN. Antagonism of toxigenic Clostridium difficile by nontoxigenic C. difficile. J Infect Dis. 1983;147(4):733–6.View ArticlePubMedGoogle Scholar
- Buckley AM, Spencer J, Maclellan LM, Candlish D, Irvine JJ, Douce GR. Susceptibility of hamsters to Clostridium difficile isolates of differing toxinotype. PLoS One. 2013;8(5):e64121. Pubmed Central PMCID: 3660315.View ArticlePubMedPubMed CentralGoogle Scholar
- Sambol SP, Merrigan MM, Tang JK, Johnson S, Gerding DN. Colonization for the prevention of Clostridium difficile disease in hamsters. J Infect Dis. 2002;186(12):1781–9.View ArticlePubMedGoogle Scholar
- Hung YP, Lin HJ, Wu TC, Liu HC, Lee JC, Lee CI, et al. Risk factors of fecal toxigenic or non-toxigenic Clostridium difficile colonization: impact of Toll-like receptor polymorphisms and prior antibiotic exposure. PLoS One. 2013;8(7):e69577. Pubmed Central PMCID: 3723847.View ArticlePubMedPubMed CentralGoogle Scholar
- Gerding DN, Johnson S, Rupnik M, Aktories K. Clostridium difficile binary toxin CDT: mechanism, epidemiology, and potential clinical importance. Gut Microbes. 2014;5(1):15–27. Pubmed Central PMCID: 4049931.View ArticlePubMedPubMed CentralGoogle Scholar
- Geric B, Carman RJ, Rupnik M, Genheimer CW, Sambol SP, Lyerly DM, et al. Binary toxin-producing, large clostridial toxin-negative Clostridium difficile strains are enterotoxic but do not cause disease in hamsters. J Infect Dis. 2006;193(8):1143–50.View ArticlePubMedGoogle Scholar
- Bacci S, Molbak K, Kjeldsen MK, Olsen KE. Binary toxin and death after Clostridium difficile infection. Emerg Infect Dis. 2011;17(6):976–82. Pubmed Central PMCID: 3358205.View ArticlePubMedPubMed CentralGoogle Scholar
- Barbut F, Decre D, Lalande V, Burghoffer B, Noussair L, Gigandon A, et al. Clinical features of Clostridium difficile-associated diarrhoea due to binary toxin (actin-specific ADP-ribosyltransferase)-producing strains. J Med Microbiol. 2005;54(Pt 2):181–5.View ArticlePubMedGoogle Scholar
- Barketi-Klai A, Monot M, Hoys S, Lambert-Bordes S, Kuehne SA, Minton N, et al. The flagellin FliC of Clostridium difficile is responsible for pleiotropic gene regulation during in vivo infection. PLoS One. 2014;9(5):e96876. Pubmed Central PMCID: 4026244.View ArticlePubMedPubMed CentralGoogle Scholar
- Brouwer MS, Roberts AP, Hussain H, Williams RJ, Allan E, Mullany P. Horizontal gene transfer converts non-toxigenic Clostridium difficile strains into toxin producers. Nat Commun. 2013;4:2601. Pubmed Central PMCID: 3826655.View ArticlePubMedPubMed CentralGoogle Scholar
- Casey TA, Bosworth BT. Design and evaluation of a multiplex polymerase chain reaction assay for the simultaneous identification of genes for nine different virulence factors associated with Escherichia coli that cause diarrhea and edema disease in swine. J Vet Diagn Investig. 2009;21(1):25–30.View ArticleGoogle Scholar
- Rinttila T, Kassinen A, Malinen E, Krogius L, Palva A. Development of an extensive set of 16S rDNA-targeted primers for quantification of pathogenic and indigenous bacteria in faecal samples by real-time PCR. J Appl Microbiol. 2004;97(6):1166–77.View ArticlePubMedGoogle Scholar
- Darling AE, Worden P, Chapman TA, Roy Chowdhury P, Charles IG, Djordjevic SP. The genome of Clostridium difficile 5.3. Gut pathogens. 2014;6(1):4. Pubmed Central PMCID: 4234979.View ArticlePubMedPubMed CentralGoogle Scholar
- Darling AE, Jospin G, Lowe E, Matsen FA, Bik HM, Eisen JA. PhyloSift: phylogenetic analysis of genomes and metagenomes. PeerJ. 2014;2:e243. Pubmed Central PMCID: 3897386.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. Pubmed Central PMCID: 2723002.View ArticlePubMedPubMed CentralGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41. Pubmed Central PMCID: 2786419.View ArticlePubMedPubMed CentralGoogle Scholar
- Bertels F, Silander OK, Pachkov M, Rainey PB, van Nimwegen E. Automated reconstruction of whole-genome phylogenies from short-sequence reads. Mol Biol Evol. 2014;31(5):1077–88. Pubmed Central PMCID: 3995342.View ArticlePubMedPubMed CentralGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.View ArticlePubMedGoogle Scholar
- Lartillot N, Philippe H. A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process. Mol Biol Evol. 2004;21(6):1095–109.View ArticlePubMedGoogle Scholar
- Stamatakis A, editor. Phylogenetic models of rate heteroginity: a high performance computing perspective. Parallel and Distributed Processing Symposium, 2006. Rhodes Island: IEEE; 2006.Google Scholar
- Stamatakis A, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol. 2008;57(5):758–71.View ArticlePubMedGoogle Scholar
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75. Pubmed Central PMCID: 2265698.View ArticlePubMedPubMed CentralGoogle Scholar
- Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014;42(Database issue):D206–14. Pubmed Central PMCID: 3965101.View ArticlePubMedPubMed CentralGoogle Scholar
- Salichos L, Rokas A. Evaluating ortholog prediction algorithms in a yeast model clade. PLoS One. 2011;6(4):e18755. Pubmed Central PMCID: 3076445.View ArticlePubMedPubMed CentralGoogle Scholar
- Frith MC, Hamada M, Horton P. Parameters for accurate genome alignment. BMC Bioinformatics. 2010;11:80. Pubmed Central PMCID: 2829014.View ArticlePubMedPubMed CentralGoogle Scholar
- Frith MC, Wan R, Horton P. Incorporating sequence quality data into alignment improves DNA read mapping. Nucleic Acids Res. 2010;38(7):e100. Pubmed Central PMCID: 2853142.View ArticlePubMedPubMed CentralGoogle Scholar
- Kielbasa SM, Wan R, Sato K, Horton P, Frith MC. Adaptive seeds tame genomic sequence comparison. Genome Res. 2011;21(3):487–93. Pubmed Central PMCID: 3044862.View ArticlePubMedPubMed CentralGoogle Scholar
- Blodel VD, Guillaume JL, Lambiotte R, Lefebvre E. Fast unfolding of communities in large networks. J Stat Mech. 2008;P10008:P10008.View ArticleGoogle Scholar
- Darling AE, Treangen TJ, Messeguer X, Perna NT. Analyzing patterns of microbial evolution using the mauve genome alignment system. Methods Mol Biol. 2007;396:135–52.View ArticlePubMedGoogle Scholar
- Rissman AI, Mau B, Biehl BS, Darling AE, Glasner JD, Perna NT. Reordering contigs of draft genomes using the Mauve aligner. Bioinformatics. 2009;25(16):2071–3. Pubmed Central PMCID: 2723005.View ArticlePubMedPubMed CentralGoogle Scholar
- Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10. Pubmed Central PMCID: 3065679.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: a fast phage search tool. Nucleic Acids Res. 2011;39(Web Server issue):W347–52. Pubmed Central PMCID: 3125810.View ArticlePubMedPubMed CentralGoogle Scholar
- Goh S, Chang BJ, Riley TV. Effect of phage infection on toxin production by Clostridium difficile. J Med Microbiol. 2005;54(Pt 2):129–35.View ArticlePubMedGoogle Scholar
- Goh S, Ong PF, Song KP, Riley TV, Chang BJ. The complete genome sequence of Clostridium difficile phage phiC2 and comparisons to phiCD119 and inducible prophages of CD630. Microbiology. 2007;153(Pt 3):676–85.View ArticlePubMedGoogle Scholar
- Horgan M, O'Sullivan O, Coffey A, Fitzgerald GF, van Sinderen D, McAuliffe O, et al. Genome analysis of the Clostridium difficile phage PhiCD6356, a temperate phage of the Siphoviridae family. Gene. 2010;462(1–2):34–43.View ArticlePubMedGoogle Scholar
- Cafardi V, Biagini M, Martinelli M, Leuzzi R, Rubino JT, Cantini F, et al. Identification of a novel zinc metalloprotease through a global analysis of Clostridium difficile extracellular proteins. PLoS One. 2013;8(11):e81306. Pubmed Central PMCID: 3841139.View ArticlePubMedPubMed CentralGoogle Scholar
- Pettit LJ, Browne HP, Yu L, Smits WK, Fagan RP, Barquist L, et al. Functional genomics reveals that Clostridium difficile Spo0A coordinates sporulation, virulence and metabolism. BMC Genomics. 2014;15:160. Pubmed Central PMCID: 4028888.View ArticlePubMedPubMed CentralGoogle Scholar
- Hennequin C, Porcheray F, Waligora-Dupriet A, Collignon A, Barc M, Bourlioux P, et al. GroEL (Hsp60) of Clostridium difficile is involved in cell adherence. Microbiology. 2001;147(Pt 1):87–96.View ArticlePubMedGoogle Scholar
- Merrigan MM, Venugopal A, Roxas JL, Anwar F, Mallozzi MJ, Roxas BA, et al. Surface-layer protein A (SlpA) is a major contributor to host-cell adherence of Clostridium difficile. PLoS One. 2013;8(11):e78404. Pubmed Central PMCID: 3827033.View ArticlePubMedPubMed CentralGoogle Scholar
- Baban ST, Kuehne SA, Barketi-Klai A, Cartman ST, Kelly ML, Hardie KR, et al. The role of flagella in Clostridium difficile pathogenesis: comparison between a non-epidemic and an epidemic strain. PLoS One. 2013;8(9):e73026. Pubmed Central PMCID: 3781105.View ArticlePubMedPubMed CentralGoogle Scholar
- Rousseau C, Poilane I, De Pontual L, Maherault AC, Le Monnier A, Collignon A. Clostridium difficile carriage in healthy infants in the community: a potential reservoir for pathogenic strains. Clin Infect Dis. 2012;55(9):1209–15.View ArticlePubMedGoogle Scholar
- Dawson LF, Valiente E, Faulds-Pain A, Donahue EH, Wren BW. Characterisation of Clostridium difficile biofilm formation, a role for Spo0A. PLoS One. 2012;7(12):e50527. Pubmed Central PMCID: 3517584.View ArticlePubMedPubMed CentralGoogle Scholar
- Deakin LJ, Clare S, Fagan RP, Dawson LF, Pickard DJ, West MR, et al. The Clostridium difficile spo0A gene is a persistence and transmission factor. Infect Immun. 2012;80(8):2704–11. Pubmed Central PMCID: 3434595.View ArticlePubMedPubMed CentralGoogle Scholar
- Ethapa T, Leuzzi R, Ng YK, Baban ST, Adamo R, Kuehne SA, et al. Multiple factors modulate biofilm formation by the anaerobic pathogen Clostridium difficile. J Bacteriol. 2013;195(3):545–55. Pubmed Central PMCID: 3554014.View ArticlePubMedPubMed CentralGoogle Scholar
- Petrovski S, Dyson ZA, Seviour RJ, Tillett D. Small but sufficient: the Rhodococcus phage RRH1 has the smallest known Siphoviridae genome at 14.2 kilobases. J Virol. 2012;86(1):358–63. Pubmed Central PMCID: 3255915.View ArticlePubMedPubMed CentralGoogle Scholar
- Sekulovic O, Garneau JR, Neron A, Fortier LC. Characterization of temperate phages infecting Clostridium difficile isolates of human and animal origins. Appl Environ Microbiol. 2014;80(8):2555–63. Pubmed Central PMCID: 3993186.View ArticlePubMedPubMed CentralGoogle Scholar