Comparative genomic analysis of toxin-negative strains of Clostridium difficile from humans and animals with symptoms of gastrointestinal disease
BMC Microbiology volume 16, Article number: 41 (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.
Strains that resided within clade C-I showed greater than 5 % sequence divergence from the reference genome at homologous sites. Clades that have diverged more than 5 % from the reference genome can be poorly resolved in workflows based on Illumina read mapping to the reference . Given this and the uncertain ancestry of the isolates included in clade C-I in the Dingle et al.  study, members of clade C-I were excluded from the subsequent analysis based on a cohort of 86 strains depicted in Fig. 1. The branching order of the five clades in our phylogenetic tree was congruent with that reported earlier  with identical clustering of strains in the different sub-clades within each of the five clades. The five toxin-negative isolates from our study clustered with strains in Clades 1 (C. difficile 5.3, and H3) and 4 (C. difficile 19.3, 22.1 and P29) that are known to contain toxin-negative strains . Our toxin-negative isolates (sourced both from animal and human sources) also clustered with strains of the same sequence type that included both toxin-positive and toxin-negative strains isolated from human clinical specimens in an earlier study .
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
In toxin-negative strains, a 115 bp sequence replaces the PaLoc locus in phylogenetic clades 1 and 4 and the genetic context in which the 115 bp sequence resides varies within the different clades that represent the C. difficile population structure . An analysis of the genetic context of the 115 bp sequence in our toxin-negative isolates compared to GenBank sequence entries representing chromosomal regions adjacent to the PaLoc insertion site in toxin-negative strains belonging to Clades 1 (GenBank Accession no HG002393) and 4 (GenBank Accession no HG002391) is depicted in Fig. 2. We identified several SNPs within the 115 bp-conserved fragment and an 80 nt long insertion in strains 5.3 and H3 downstream of the gene designated CD06642 in the reference genomes (Fig. 2). In addition, we identified a 68 nt tandem repeat of the sequence adjacent to the 80 nt insertion site in strain H3 (see Additional file 2: Figure S7). Within members of clade 4, porcine isolate P29 had a deletion of the hypothetical gene seen in the reference region while human isolates 19.3 and 22.1 had a significant decrease in the nucleotide identity in hypothetical genes in the reference genome. Consistent with observations reported earlier, these differences indicate ongoing micro-evolutionary events within the locus that flanks the PaLoc insertion site .
Comparative BLASTp analyses of toxin-negative strains to identify functional similarities between groups of isolates
An all versus all BLASTp based protein-clustering analysis was used to identify the putative core proteome of a subset of Clade 1 strains comprising five genomes and Clade 4 strains comprising 10 genomes (see boxed region in Fig. 1). To maintain uniformity in the input data used in the analysis, raw reads representing each of the genomes from an earlier study  were reassembled and annotated using the protocols that were used to assemble and annotate the five toxin-negative Australian isolates of C. difficile reported in this study, as described above. The assembly statistics and a comparison of the assembly outputs are presented in Additional file 3: Table S8. On an average, RAST predicted 3700 proteins per genome and these were included in the protein clustering analysis. A bit score cut off of 50 was used to cluster homologous protein sequences. An overview of the data is presented in Table 1.
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).
We queried the viral and prophage database in GenBank with the genome sequences of all five toxin-negative isolates using PHAST . The database confirmed the presence of several regions contain phage DNA in each of the five genomes in Clade 1 (see Additional file 2: Figure S10 and Additional file 4: Table S11). Table 2 lists a summary of the complete prophage sequences (PHAST scores > 110 and over) identified in the five Australian toxin-negative isolates. Notably, isolates 19.3 and 22.1 returned identical phage profiles (Additional file 2: Figure S10). Both genomes carried an identical and complete 56.8 kb phage that is a close relative of phiC2, a 56.5 kb phage first identified in C. difficile strain CD242 . phiC2 is one of the first completely sequenced temperate phages in C. difficile and regions of phiC2 are present in the majority of C. difficile genomes of clinical origin . We also identified a phage in isolates P29 and H3 that carries sequence identity with phiC2. The prophage in P29 was larger (97.4 kb) than the original phiC2 sequence (Additional file 2: Figure S10). Isolate P29 also carries two other complete prophage sequences. One of these, prophage region 9 has significant sequence identity to the C. difficile temperate bacteriophage phiCD6356 that belongs to the Siphoviridae family . The Siphoviridae family prophage identified in P29 is 52.8 kb (Additional file 4: Table S11) and is considerably larger than the first report of this phage sequence at 37.6 kb . Isolate H3 also carries an identical variant of phage phiCD6356 (Additional file 2: Figure S12) comprising 41.1 kb. Evidence of other complete prophage genomes is listed in Table 2. Comparative BLASTp analysis of the four genomes (19.3, 5.3, H3 and P29) also aligned to the phage search protocol and confirmed the data generated by PHAST analyses, reinforcing the observation that the major differences in our Australian toxin-negative C. difficile isolates have a prophage origin.
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
A homology-based functional prediction analysis of proteins that have been experimentally verified, or predicted to play a role in the colonization of C. difficile strain CD630 with homologous proteins in the genomes of the toxin-negative isolates in our study is shown in Table 3 [70, 71]. Most of the proteins i had near perfect (98–100 %) protein sequence identity including proteins encoded by spo0A, which serves as a positive regulator for genes required in spore germination and groEL, a chaperone that also functions as an adhesin . Several serine-proteases and other metalloproteases which may contribute to the disease development process  were also highly conserved. However, some membrane associated proteins that have been shown to play significant roles in the disease development process had lower alignment scores including SlpA, a surface layer protein that is proposed to facilitate host cell adherence  and FliC, an adhesin necessary for the colonization of gut epithelium .
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.
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.
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.
Peterson LR. Bad bugs, no drugs: no ESCAPE revisited. Clin Infect Dis. 2009;49(6):992–3.
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.
Rupnik M, Wilcox MH, Gerding DN. Clostridium difficile infection: new developments in epidemiology and pathogenesis. Nat Rev Microbiol. 2009;7(7):526–36.
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.
Bauer MP, Kuijper EJ. Potential sources of Clostridium difficile in human infection. Infect Dis Clin North Am. 2015;29(1):29–35.
Songer JG, Anderson MA. Clostridium difficile: an important pathogen of food animals. Anaerobe. 2006;12(1):1–4.
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.
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.
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.
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.
Cairns MD, Stabler RA, Shetty N, Wren BW. The continually evolving Clostridium difficile species. Future Microbiol. 2012;7(8):945–57.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Rupnik M. How to detect Clostridium difficile variant strains in a routine laboratory. Clin Microbiol Infect. 2001;7(8):417–20.
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.
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.
Matamouros S, England P, Dupuy B. Clostridium difficile toxin expression is inhibited by the novel regulator TcdC. Mol Microbiol. 2007;64(5):1274–88.
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.
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.
Collins DA, Elliott B, Riley TV. Molecular methods for detecting and typing of Clostridium difficile. Pathology. 2015;47(3):211–8.
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.
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.
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.
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.
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.
Wilson KH, Sheagren JN. Antagonism of toxigenic Clostridium difficile by nontoxigenic C. difficile. J Infect Dis. 1983;147(4):733–6.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.
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.
Stamatakis A, editor. Phylogenetic models of rate heteroginity: a high performance computing perspective. Parallel and Distributed Processing Symposium, 2006. Rhodes Island: IEEE; 2006.
Stamatakis A, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol. 2008;57(5):758–71.
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.
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.
Salichos L, Rokas A. Evaluating ortholog prediction algorithms in a yeast model clade. PLoS One. 2011;6(4):e18755. Pubmed Central PMCID: 3076445.
Frith MC, Hamada M, Horton P. Parameters for accurate genome alignment. BMC Bioinformatics. 2010;11:80. Pubmed Central PMCID: 2829014.
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.
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.
Blodel VD, Guillaume JL, Lambiotte R, Lefebvre E. Fast unfolding of communities in large networks. J Stat Mech. 2008;P10008:P10008.
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.
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.
Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10. Pubmed Central PMCID: 3065679.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The authors declare that they have no competing interests, including financial competing interests.
PRC assembled the genome sequences, conceived, analysed and designed the presentation of data in the manuscript, performed all comparative genomic analysis and benchwork to confirm features identified from the comparative genomic analysis and drafted the manuscript. AED and MD wrote data analysis software packages and scripts and assisted with data interpretation. TC cultured the isolates, performed PCR analyses for toxin genes strains, ribotyping and fingerprinting assays and designed aspects the project. PW constructed the libraries for genome sequencing. IGC assisted with experimental design, data analysis and writing the manuscript. SPD conceived and designed aspects of the project, managed the project and drafted the final manuscript. All authors read and approved the final manuscript.
Isolation history, genomic assembly and initial PCR screening results of the five isolates analysed in this manuscript. (DOCX 36 kb)
Preliminary phylogeny with reference genomes to identify the most closely related reference genome for tiling of genomic scaffolds. Figure S3. Mauve alignments of genomes 19.3 and 22.1. Figure S7. DNA sequence alignment of repetitive regions in PaLoc of isolates in clade 1. Figure S10. Phage profiles in the genomes of the five Australian toxin-negative C. difficile isolates included in this study. Figure S12. Preliminary Phylogenetic tree with all isolates included in Dingle et al’s publication in 2011 . Isolates that form clade C-I is highlighted in blue. (DOCX 2963 kb)
SNPs identified by Mauve alignments of genomes of C. difficile strains 19.3 and 22.1. Table S5. Bidirectional BLASTp analysis of peptide sequences predicted from genomes of C. difficile strains 19.3 and 22.1 to identify amino acid changes caused by SNPs in the genomes. Table S8. Assembly Statistics downloaded from short read archives. Table S9. Unique proteins identified from the All Vs All BLASTp analysis of C. difficile genomes. (XLSX 632 kb)
Regions identified in each C. difficile genome by PHAST. Sheet 1. C. difficile strain P29 genome, Sheet 2. C. difficile strain H3 genome, Sheet 3. C. difficile strain 5.3 genome, Sheet 4. C. difficile strain 19.3 genome. Sheet 5. C. difficile strain 22.1 genome. (XLSX 134 kb)
Sheet 1 = Comparative BLASTp analysis of 4 genomes, Sheet 2 = Alignment of scaffold 11.1 of genome P29 against C. difficile CD630 genome. (XLSX 682 kb)
About this article
Cite this article
Roy Chowdhury, P., DeMaere, M., Chapman, T. et al. Comparative genomic analysis of toxin-negative strains of Clostridium difficile from humans and animals with symptoms of gastrointestinal disease. BMC Microbiol 16, 41 (2016). https://doi.org/10.1186/s12866-016-0653-3