- Research article
- Open Access
Use of single molecule sequencing for comparative genomics of an environmental and a clinical isolate of Clostridium difficile ribotype 078
BMC Genomicsvolume 17, Article number: 1020 (2016)
How the pathogen Clostridium difficile might survive, evolve and be transferred between reservoirs within the natural environment is poorly understood. Some ribotypes are found both in clinical and environmental settings. Whether these strains are distinct from each another and evolve in the specific environments is not established. The possession of a highly mobile genome has contributed to the genetic diversity and ongoing evolution of C. difficile. Interpretations of genetic diversity have been limited by fragmented assemblies resulting from short-read length sequencing approaches and by a limited understanding of epigenetic regulation of diversity. To address this, single molecule real time (SMRT) sequencing was used in this study as it produces high quality genome sequences, with resolution of repeat regions (including those found in mobile elements) and can generate data to determine methylation modifications across the sequence (the methylome).
Chromosomal rearrangements and ribosomal operon duplications were observed in both genomes. The rearrangements occurred at insertion sites within two mobile genetic elements (MGEs), Tn6164 and Tn6293, present only in the M120 and CD105HS27 genomes, respectively. The gene content of these two transposons differ considerably which could impact upon horizontal gene transfer; differences include CDSs encoding methylases and a conjugative prophage only in Tn6164. To investigate mechanisms which could affect MGE transfer, the methylome, restriction modification (RM) and the CRISPR/Cas systems were characterised for each strain. Notably, the environmental isolate, CD105HS27, does not share a consensus motif for m4C methylation, but has one additional spacer when compared to the clinical isolate M120.
These findings show key differences between the two strains in terms of their genetic capacity for MGE transfer. The carriage of horizontally transferred genes appear to have genome wide effects based on two different methylation patterns. The CRISPR/Cas system appears active although perhaps slow to evolve. Data suggests that both mechanisms are functional and impact upon horizontal gene transfer and genome evolution within C. difficile.
Clostridium difficile (reclassified as Clostrioides difficile ) is an enteric pathogenic bacterium that can cause symptomatic disease, which ranges in severity from fever and diarrhoea to the development of pseudomembranous colitis and toxic megacolon . Clostridium difficile infection (CDI) occurs following antibiotic treatment as new ecological niches become available upon disruption of the normal microbiota . CDI may arise from ingested endospores transmitted via the faecal oral route, or from vegetative cells already present in the patient, as the bacterium can be asymptomatically carried in adults and children . CDI may also be contracted outside the hospital setting , and C. difficile has been isolated from food products [4–6], on surfaces around the home  and from swimming pools . It has also been isolated from the natural environment including river water, soils, sea water and estuarine sediments [7–10]. The presence of C. difficile at these sites may be due to contamination with sewage or agricultural run-off, yet bacteria from these locations could be re-introduced to the food chain, for example via contaminated shellfish or seafood [11, 12], and they have been implicated in the infection of marine mammals .
The movement of C. difficile between reservoirs is particularly pertinent for isolates of the PCR ribotype 078 (R078). This is an epidemic strain, first identified in livestock and subsequently in clinics across Europe . Although pathogenic, it is not clear quite how much virulence versus strain fitness shapes which strains come to prominence in the hospital environment [15, 16]. R078 strains form a lineage divergent from other major ribotypes , as also determined via multilocus sequence typing (MLST) analysis [18, 19] and core genome phylogenies [20, 21]. Previously, we isolated a R078 strain, CD105HS27, from estuarine sediment  and sequenced its genome using Illumina HiSeq 2000 generating a draft assembly . The carriage of transposon Tn6293 (previously unnamed) and the absence of Tn6164 was confirmed in this study from the results of the Single Molecular Real Time (SMRT) sequencing. The accessory gene content in C. difficile as a species is high relative to the size of its core genome , and it is characterised by multiple mobile genetic elements which include transposons, integrated conjugative elements, plasmids and prophages (for recent reviews see [23–25]). The acquisition of antibiotic resistance and novel virulence factors are thought to drive C. difficile strain pathogen evolution , but its ecology outside of the human host is little understood.
Recently, SMRT technology has been applied to sequence C. difficile genomes, exploiting the long read data to determine chromosomal structure, mobile genetic content and methylation patterns [27–31]. The re-sequencing of previously analysed strain CD630 showed differences in its ribosomal operon, transposon and tRNA content [28, 31]. In this study we first determined if re-sequencing the reference strain M120 (R078) using SMRT would reveal differences in the chromosomal architecture. Next, we compared SMRT generated genome sequences of M120 with CD105HS27 in order to gain a better understanding of the differences between an environmental isolate and a clinical strain. To date, SMRT has not been applied to isolates of R078. In addition to analysing the genomic data, we compared methylation patterns across the genome. Due to the fact that the CRISPR/Cas system also can provide immunity to invading DNA elements, we assessed its potential to target MGEs for each strain. In both cases, understanding mechanisms that govern horizontal gene transfer in C. difficile provides insight into the genome evolution of this pathogen.
Results and discussion
Genome features of M120 and CD105HS27
The two genome assemblies generated using SMRT are in near-complete condition; the genome of M120 is 4,082,634 bp with an average coverage of 16.3×, an average 28.73% GC content, and is comprised of two contigs of 4,069,609 bp and 13,024 bp in length. The total sequence for CD105HS27 is 4,122,476 bp, with an overall coverage of 15.75× and an average 29.15% GC content, and consists of five contigs of 3,462,540 bp, 339,877 bp, 174,028 bp, 146,675 bp and 1156 bp, respectively.
Both assemblies were compared to the reference genome of M120, which is a single chromosome 4,047,729 bp in length with an average 28.76% GC content. The 13,024 bp size contig contains a set of 5S, 16S and 23S rRNA genes and 19 tRNA genes, and has duplicated region encoding an identical tRNA (Alanine) and 16S rRNA gene (dot plot data not shown), in addition to predicted CDSs encoding glycosyl transferases, DNA polymerase subunit and recombination protein RecR. The relative coverage of this contig is on average 1.3× (see Fig. 1). To determine whether this contig represents a sequence mobilization event and a low copy number requires experimental investigation.
Annotation of the re-sequenced M120 genome identified 3541 CDSs, 101 tRNAs and 39 rRNAs; this is consistent with the reference genome, but includes an additional 15 tRNAs and 7 rRNA genes. Similar observations were seen in a SMRT sequenced genome of CD630Δerm with additional tRNA and rRNA genes located in a novel ~5 kbp insertion . This was attributed to adaption during laboratory culture as extra ribosomal gene operon copies have been shown to affect fitness in E. coli with regards to nutrient availability . Furthermore, recombination events have been suggested as a mechanism for generating the diversity of ribotypes in C. difficile .
The genome of CD105HS27 has 3598 CDSs, 93 tRNA and 47 rRNA genes. The chromosome breaks are located in regions encoding ribosomal genes, which appear to have undergone duplication events across the genome. The application of SMRT can also improve the assembly of other regions containing repeat sequences. For example, previously, toxin gene carriage had been confirmed by PCR for CD105HS27 , but an Illumina generated draft genome assembly of its genome resulted in fragmented versions of tcdA and tcdB . Here, these genes have been resolved fully. CD105H27 has 79 CDSs that are not present in M120, most of which are encoded on Tn6293, In contrast, M120 has 103 CDSs that are not present in CD105HS27, of which 102 are encoded on Tn6164. The predicted genetic content of these two transposons suggests that they may be conjugative transposons although this has yet to be demonstrated experimentally. Therefore, these should be re-termed as putative conjugative transposons CTn6164 and CTn6293. Tn6164 is a large (~100 kbp) element that appears to be two MGEs including a prophage region which shares similarity to the Streptoccocus conjugative phage Φ1207. 3 . Φ1207. 3 has been demonstrated to transfer between strains via conjugation and was originally annotated as a conjugative transposon  but contains conserved phage genes including those predicted to encode terminases, capsid, tail and holin proteins leading to its re-designation as a conjugative prophage . Prophages transmitting via conjugation appear rarely in the literature (e.g. ). Whether these prophages also transfer via conjugation has not been established, however their discovery suggests that this mechanism may occur more widely than previously known.
The two genomes are related, sharing an average nucleotide identity of 99.98% based on the whole genome sequence (following the method described in ). Alignment of the whole genomes using MAUVE and its SNP (Single Nucleotide Polymorphism) detection tool showed that the aligned sequences differed in 85 positions by single nucleotide changes. Further comparison of the two genomes via BLASTn (Fig. 1) and within a dotplot (Fig. 2) revealed extensive sequence similarity between the two strains, with exceptions of two large indel (insertion-deletion) regions (~100 kbp) that carry the putative CTns Tn6164 and Tn6123, the movement of Tn6190, and inversion rearrangements. Use of SMRT has previously shown major chromosomal rearrangements from resequencing the genome of strain CD630 in addition to duplication of ribosomal gene operons [28, 30]. One mechanism for these rearrangements are the movement of the MGEs, as seen in the mutant CD630Δerm, where the re-mobilisation of transposon CTn5 led to the inversion of the genome sequence . What affect such chromosomal re-engineering has on the physiology of the cell in terms of gene expression is not known, but may be significant as has been described for the control of DNA elements from the chromosome in the regulation of diverse bacterial processes .
In silico typing of M120 and CD105HS27
In C. difficile, ribotyping is one of the main methods used to categorise strains. In silico ribotyping was performed to assess the outcomes from the SMRT generated genomes and to explain how the duplication events affect the ribotypes profile. As expected from the different numbers of total rRNA genes, the two profiles differ, with 11 bands predicted from M120 reference, 12 from M120 SMRT and 16 from CD105HS27 (Additional file 1: Table S1). The profiles differ by duplication of identical sized regions in addition to bands of different lengths which may affect ribotypes assigned. While ribosomal gene regions assemble poorly in Illumina datasets, the ability to generate near complete genomes using SMRT technology show how ribosome operon duplication and recombination events could be tracked.
Another method used to type C. difficile is MLST (multilocus sequence typing), a scheme that compares the sequence data for seven conserved genes . The two isolate genomes were assigned to Sequence Type (ST) 11, clade 5, which is consistent with previously typed isolates of R078 [19, 40, 41]. The C. difficile MLST tool also analysed additional key genes, such as toxins Toxin A, Toxin B and the CDT and also genes that encode for antibiotic resistance. The results confirmed both M120 and CD105HS27 have wild type toxin genes cdtAB and tcdB and a 39 bp deletion in tcdC which has been characteristic of R078 isolates from its early identification . Furthermore, tetM, predicted to encode a ribosomal protection protein (CDM120_RS02595) carried on Tn6190 in M120 , is absent in CD105HS27, which has two copies of a variant tetM, that share 67% identity at the aa level to that in M120.
Mobile genetic element content of M120 and CD105HS27
Like other isolates, those from R078 have been found to carry different sets of MGEs which encode for predicted virulence factors and antibiotic resistance [24, 25]. These include the conjugative transposons related to those in other strains of C. difficile; Tn6073 (CTn1-like), Tn6107 (CTn5-like) and CTn4 in the clinical R078 strain QCD-23 M63 , as well as those more distantly related, such as Tn6164 in the reference strain M120 . Tn6164 is a composite MGE containing a prophage and has several regions that originate from different bacterial lineages . This is considered likely to be a transposon as it can excise and circularise, and carries genes encoding products predicted to be involved in conjugation . While Tn6164 is characteristically associated with R078 strains, not all R078 isolates carry it . R078 isolates also may harbour Tn6190 (previously termed CTnCD3a ), a Tn916-related element that carries the tetracycline resistance gene tetM , as well as Tn6235 which carries aphA1, an aminoglycoside 3′-phosphotransferase suggested to confer streptomycin resistance . M120 and CD105HS27 both have Tn6190, but, as described previously, M120 has Tn6194 whereas CD105HS27 does not. However, the environmental isolate does have a different large ~104 kbp element , now assigned as Tn6293. Encoded on Tn6293 are several genes with predicted functions that could potentially enhance cell survival and growth, including homologs of aadE (which confers aminoglycoside resistance ), a LexA repressor (involved in the SOS response regulation ) and 23S rRNA methyltransferase RlmN (that could impact on cellular growth ). It has predicted transposases and conjugation transfer genes as well as homologs of plasmid maintenance and replication protein encoding genes; parA and parB, and repA, suggesting this MGE is also a composite with several origins as determined for other C. difficile transposons, Tn9194 and Tn6103 [34, 42]. Interestingly, the amino acid sequence of AadE was 100% identical to that of plasmid-carried aadE genes in Campylobacter jejuni (YP_009079621) and Pediococcus acidilactic (YP_001965484), and is present in several Firmicutes sp. sequences from WGS projects, further supporting prior observations that this resistance can transmit between bacterial genera . To determine the carriage of Tn6293 in C. difficile, its sequence was searched using BLASTn against C. difficile (taxid 1496) sequences. Homologous regions were found in the genomes of three of the seven isolates that are related to M120 (Additional file 2: Table S2); E1 and T5 (R126, human isolates) and NAP08 (R078, human isolate) . To determine its potential origin, the nt sequence was searched against the NCBI nt/nr db. It has similarity to regions in Eubacterium and Ruminococcus spp. genomes. The shared nt sequence similarity is primarily located in genes whose predicted products are involved in genetic element mobilisation and maintenance functions. These include a serine recombinase (CD105HS27_00591), DNA binding and mobilization proteins (CD105HS27_00611 and CD105HS27_00612) and plasmid recombinase (CD105HS27_00634). Both Eubacterium sp. and Ruminococcus sp. belong to the same order as C. difficile, the Clostridiales, and the shared sequence similarity observed supports previous findings of MGEs being exchanged between these genera .
Both genomes carry a predicted R-type bacteriocin. R-type bacteriocins resemble phage tail-like particles (PTLPs) and have genes predicted to encode proteins involved in structural roles for tail assembly. However, they lack predicted capsid genes and thus are not a complete virion particle. These bacteriocins, or PTLPs, have been observed in culture supernatants of diverse isolates [9, 47, 48], and been used either as typing tools or to determine their use as alternative therapeutics [49, 50]. Due to the specificity required of proteins that target the cell surface, obtaining sequence information from the genomes of clinically relevant strains could aid in using a synthetic biology approach for designer antimicrobials; this has been demonstrated for the bacteriocin carried in a R027 isolate , with subsequent genetic modification for enhancing its antimicrobial application .
It is not possible to conclude whether these strains have transferred from the environment to the patients or vice versa from the comparisons we have performed here based on a sample size of two. However, the putative origins of these CTns have been examined based on sequence homology. Tn6164 and Tn6293 are clearly distinct from one another, and to known elements in other bacterial species. For example, for Tn6164, similarity to other sequences is split over the length of the transposon into at least two major regions: the phage containing region is most closely related to a single Clostridium difficile genome Z31 (CP013196.1) based on a nt identity of 93% covering 35% of its length. In the same region, the next most closely related elements are found in the complete genome of Thermoanaerobacter spp. (CP002210 and CP000923.1) and a draft genome of Clostridium bornimense (GCA_000577895). Thermoanaeroacter strains were originally isolated from anaerobic enrichments with environmental samples from subsurface. C. bornimense is a hydrogen producing Clostridium and this species does not have an associated history with human infections, but isolated from a laboratory bioreactor . The second region of the transposon has homology to Streptococcus and Anaerococcus spp. In contrast, Tn6293 showcases sequence similarity in multiple regions across its full length to different bacterial genera including Ruminococcus, Clostridium and Eubacterium spp. It is interesting that the second region of homology in Tn6164 is to pathogenic species. However, as this is based on few sequences, it is not possible to conclusively state this has been acquired while in clinics despite its absence from CD105HS27 (and thus infer CD105HS27 has evolved outside of clinics). Whether the two isolates have evolved in isolation is one possibility. SNP analysis has been used to track the transfer of strains across the world  and in different reservoirs [19, 54], with estimated mutation rates of 1–2 sites per year, suggesting that the number of substitutions (n = 85) we observed here suggests that these two isolates have evolved from one another over some time. Increasing numbers of R078 genomes will aid in determining the movement of strains from clinics to the environment and vice versa, in addition to how these strains further evolve when in different reservoirs.
Methylome of R078 isolates
To establish genome-wide methylation patterns of the two isolates, the profiles for methylation modifications N4-methylcytosine (m4C) and N6-methyladenine (m6A) were analysed from the SMRT data . Methylation (the addition of methyl groups to bases) in bacteria may play a regulatory role in terms of gene expression , but is also one way that DNA elements can exploit to protect against their degradation by restriction modification systems in the host cell . Both strains M120 and CD105HS27 show adenine methylation of the consensus sequence CAAAAA with high efficiency of target methylation (7484/7579, or 98.75% sites in M120 and 7469/7559 or 98.8% in CD105HS27). This target specificity had been previously assigned to the N6-adenine methyltransferase named M.Cdi25 or Cdi630V (locus tag CD630_27580, protein Id YP_001089271.1) of strain CD630  and is reported in the REBASE database . The respective methyltransferases of M120 (CDM120_RS14295, WP_003422891.1) and CD105HS27 (CD105HS27_02520) are identical and show a 98% identity (565/577) to the CD630 orthologue. Strain M120 showed signatures for a N4 modified cysteine ACGGC methylation target (398/414) and a consensus sequence CGGCNTGTGNNNNNNT was identified but with unknown modified base calls (12/13). In REBASE, the ACGGC target is assigned to two tandem methyltransferases of Tn6164, M1.CdiMORFAP (CDM120_RS02255, WP_041160334.1) and M2.CdiMORFAP (CDM120_RS02260, WP_041160335.1). No further modified base was detected in strain CD105HS27. The finding that methylation pattern of m4C GCCGT/ACGGC was absent in CD105HS27 may be explained by the absence of Tn6164 and both these two methyltransferases. In contrast, both M120 and CD105HS27 encode CdiMORFEP, a homolog of M.CdiG46II (amino acid identity of 565/577 (98%)) which is predicted to recognise CAAAAA sites. Three further predicted methylases on Tn6164 are present in M120  and absent from CD105HS27, as the latter lacks this mobile element. While it was expected for the two Tn6164 m5C methyltransferases M.CdiMORFBP (CDM120_RS02360, WP_041160353.1) and M.CdiMORFCP (CDM120_RS02725, WP_041160386.1) to show no signature on the SMRT dataset, we would have expected to identify a signature for the putative m6A methyltransferase (CDM120_RS02520, WP_000662263.1). The fact that no additional adenine methylation pattern was detected could be due to one of many reasons including target identity of this enzyme and M.Cdi25/Cdi630V, lack of expression of the enzyme in CD105HS27 or inappropriate annotation of predicted CDSs.
Just as there are different sets of methylation genes functional in C. difficile, strains carry genes encoding multiple restriction enzymes . It is of interest to note that despite the fact that M120 and CD105HS27 are highly related, they share only core genome methylation systems as the adenine methylase above or the McrBC system, as they do with the strain CD630. This is due to the fact that the majority of methyltransferases are in Tn6164 which is absent from CD105HS27. In addition to methylation Restriction Modification (RM) systems, MGEs have other defence systems against super-infection . Here, Tn6164 carries three putative methylase genes on the transposon region and two on the prophage region of the element. The two sequenced strains were also found to contain defence mechanisms to combat RM systems, notably, Tn6190 carries ardA which encodes ArdA, an anti-restriction protein for type I restriction systems . Whether this system is active remains to be determined, but evidently there are multiple mechanisms employed by MGE in C. difficile to be maintained.
CRISPR/Cas system of M120 and CD105HS27
Immunity to phage infection can also be conferred via the CRISPR (Clusters of Regularly Interspaced Palindromic Repeats)/Cas system which works as an RNA based interference against invading DNA elements , but also may act as regulatory machinery for other aspects of the cell biology and genome evolution . The function of the CRISPR/Cas system depends on the action of CRISPR associated (Cas) proteins that are highly diverse in operons across prokaryotes, and ultimately involves the processing and matching of spacers to target DNA with its subsequent restriction . It comprises of arrays that have conserved direct repeat (DR) sequences that flank spacer sequences. Spacers are homologous to phage or plasmid sequences as have been incorporated into arrays following unsuccessful past invasions, and in this way they can provide information about past interactions with such elements .
In this study, six CRISPR arrays and three cassettes of Cas genes were identified in each genome. Two Cas gene operons belonged to the I-B/TNeap group and contained all gene components to be functionally complete , and the third set comprised of cas6, cas7, cas5 and cas3, but lacked cas1 and cas2. Multiple cas sets within a single genome, of both complete and incomplete operons, have been described previously in C. difficile strains CD630  and R20291, but it appears unusual that these two isolates have two complete yet distinct cassettes. The two complete sets are adjacent to CRISPR arrays CRISPR 4 and CRISPR 5.
The six CRISPR arrays are conserved between the two isolates. Five of the arrays have identical spacer contents with 17 (CRISPR_1), 44 (CRISPR 2), 13 (CRISPR 3), 32 (CRISPR 5) and 9 (CRISPR 6) spacers. The remaning array, CRISPR 4, has one additional spacer in CD105HS27 than M120, with 39 and 38 spacers, respectively (spacer number 12, indicated in by Additional file 3: Table S5. by asterisk). Previously, we showed that spacers targeted C. difficile phages . Here, we searched spacers from the six arrays against 20 C. difficile phage genomes (Fig. 3, Additional file 4: Table S3). Of the total 154 spacers present in both isolates, 19 spacers have at least one identical match to a phage sequence from 18 phages. Perfect matches were identified between spacers and phage sequences from all arrays, except CRISPR arrays 3 and 6. Spacers with matches were located throughout the arrays, but differed with regards to location and type of phage (Fig. 3). We focused on perfect matches as phages phiCDHM1, phiCDHM19, phiCDHM14 and phiCDHM13 do not produce lysis of either strain . To identify matches for the remaining spacers and to a wider range of DNA sequences, we searched the viral and plasmid databases in CRISPRTarget , the metaviromic datasets publically available on MetaVir  and C. difficile genomes (Additional file 4: Table S3 and Additional file 5: Table S4). We did not detect any perfect matches to the viromic datasets, but identified matches for spacers from all six CRISPR arrays to prophage and phage-like genes in the C. difficile bacterial genomes (Fig. 4, Additional file 3: Table S5). It has been found that CRISPR systems may also have regulatory roles in genomes . To identify if there were spacers that matched to genomic sequence, we searched the genome of CD105HS27 and identified one perfect match for a spacer in CRISPR 6. The protospacer sequence is located in CD105HS27_02420, a gene encoding a putative carboxylase. This does not have either of the previously identified CCT or CCA Protospacer Adjacent Motif (PAM) sequences  so whether this has a functional role is unknown.
We see that in C. difficile, CRISPR arrays appear to undergo horizontal exchange between strains via their presence on MGEs, including prophage, plasmids and the C. difficile sigK intervening (skinCd) element [18, 66]. In the genome of C. cellulolyticum H10, two CRISPR arrays are proximal to a transposase gene which suggests that recombination events could shift immunity profiles via the introduction of novel arrays with new spacer content . Similarly in M120 and CD105S27, two of the arrays, CRISPR 1 and CRISPR 2, are in proximity to CDSs that suggest past integration events containing either integrase or transposase domains. Whether these genes still function and these regions are mobile is not clear from annotation alone. However, these findings of arrays on MGE and signatures of past integration events nearby suggest that arrays could move following genome insertion and excision events by a variety of mechanisms.
SMRT technology has been used to generate near complete genomes for two R078 strains, allowing the comparison of clinical and environmental isolates. The two genomes differ in chromosomal structure and number of ribosomal operons. Additionally, the two genomes differ in the carriage of two transposons, Tn6164 in M120 and Tn6293 in CD105H27, which we suggest are termed as putative conjugative transposons CTn6164 and CTn6293.
The majority of unique genes are carried on the two putative CTns and include predicted methylases. The methylome analysis for each genome suggests a vastly different methylation pattern with no consensus m4C motif in CD105HS27 detected. This likely impacts the immunity of each isolate to DNA elements including phages, and to the type of HGT that may occur for each. In contrast, their CRISPR/Cas systems are highly similar with only one spacer different between the two. Our findings support previous work that the CRISPR/Cas and RM systems are not mutually exclusive , and show this indeed appears to be the case in C. difficile.
Bacterial genomic DNA extraction
Bacterial genomic DNA (gDNA) extraction was performed using 1 ml overnight culture from a single colony grown in Brain Heart Infusion (BHI) broth (Oxoid, UK). DNA was extracted using a Qiagen GenomicTip 500/G kit (Qiagen, UK) following the manufacturer’s instructions. Pulsed Field Gel Electrophoresis was performed to assess gDNA degradation, with 100 ul of each sample separated on a 1% Agarose gel (Manufacturer info) for 18 h at 6 V. Gels were stained with 10 ul of ethidium bromide and visualised using UV G Box, Syngene. Sample gDNA quantity and quality was measured using by Qubit assay on a Qubit fluorometer (Life Technologies, USA) according to the manufacturer’s instructions, and by measuring absorbance at 260 nm and 280 nm using a Nanodrop Spectrophotometer (Thermo Scientific, UK).
Genome sequencing and bioinformatics analysis
Genomic DNA sequencing using a SMRT Pacific Biosciences platform was performed at the Centre for Genomic Research, University of Liverpool. SMRTbell libraries were prepared by Margaret A. Hughes with 3 SMRT cells used per library for sequencing. High quality genome assemblies were generated using HGAP (Hierarchal Genome Assembly Processer) as part of the SMRT Portal and methylation patterns detected. Contig structure and plasmid identification was performed from dotplots generated using Gepard .
Genomes were visualised using Artemis Genome Browser . Coverage was determined from alignment of the corrected reads to the final assembly using BWA-SW , and samtools for index and conversion of file formats . Coverage was assessed using Qualimap v.1.0  and coverage plots were generated using the Artemis DNAplotter perl script . Genome annotation was performed using PROKKA v1.7 , with a custom guide database containing proteins from the reference genome of M120 (accession NC_017174.1). RNA genes were predicted using RNAmmer v1.2 . In silico ribotypes profiles were predicted using the oligonucleotide sequences from Bidet et al. . Shared gene content was identified with blast + v2.2.28 using blast-all-v-all . This publication made use of the Clostridium difficile Multi Locus Sequence Typing website (http://pubmlst.org/cdifficile/) developed by Keith Jolley and sited at the University of Oxford . The characterised C. difficile CD630 CTns were used as a reference set for the identified of similar MGEs by BLASTn. Whole genome alignment and single nucleotide differences were generated using MAUVE v.2.4.0 . Average nucleotide identity was calculated following the method described in , using the online web based tool which can be accessed at http://enve-omics.ce.gatech.edu/ani/ with parameters of min. length 700 bp, min. identity 70% and min. alignment 50. Dotplot analysis was generated using Gepard . Genome comparison maps were generated using EasyFig v.2.2.2 . Restriction modification systems were analysed using entries from REBASE (the Restriction Enzyme database) . Prophage regions were predicted using PHAST . CRISPR arrays were identified using CRISPRfinder , and the genomes CRISPR content compared using CRISPRcompar . Spacer sequences were searched against the GenBank-Phage, RefSeq-Plasmid, RefSeq-Viral and Genbank-Environmental databases (accessed 1/10/2015) using CRISPRTarget  in addition to virus metagenome datasets (Additional file 5: Table S4). Spacer protein targets were identified using a curated approach based on annotations on the NCBI genome browser at locations identified from the CRISPRTarget search. Where no annotation was available from perfect spacer-target matches on CRISPRTarget, consensus annotations from imperfect matches (up to 7 mismatches) were used.
CRISPR associated proteins
Clostridium difficile infection
Coding DNA sequence
Clostridium difficile binary toxin
Clustered regularly interspaced short palindromic repeats
Long tailed myovirus
Mobile genetic element
MultiLocus sequence typing
- NCBI nr/nt db:
National Center for Biotechnology Information non-redundant nucleotide database
Protospacer adjacent motif
Polymerase chain reaction
The Restriction Enzyme database
Ribosomal ribonucleic acid
- skincd :
Clostridium difficile sigK intervening element
Single molecule real time
Single nucleotide polymorphism
Transfer ribonucleic acid
Lawson PA, Citron DM, Tyrrell KL, Finegold SM. Reclassification of Clostridium difficile as Clostridioides difficile (Hall and O’Toole 1935) Prévot 1938. Anaerobe. 2016;40:95–9.
Leffler DA, Lamont JT. Clostridium difficile Infection. New Engl J Med. 2015;372(16):1539–48.
Gerding DN. Clindamycin, cephalosporins, fluoroquinolones, and Clostridium difficile-associated diarrhea: This is an antimicrobial resistance problem. Clin Infect Dis. 2004;38(5):646–8.
Eyre DW, Griffiths D, Vaughan A, Golubchik T, Acharya M, O’Connor L, Crook DW, Walker AS, Peto TE. Asymptomatic Clostridium difficile colonisation and onward transmission. PLoS One. 2013;8(11):e78445.
Metcalf DS, Costa MC, Dew WMV, Weese JS. Clostridium difficile in vegetables, Canada. Lett Appl Microbiol. 2010;51(5):600–2.
Weese JS, Avery BP, Rousseau J, Reid-Smith RJ. Detection and Enumeration of Clostridium difficile Spores in Retail Beef and Pork. Appl Environ Microbiol. 2009;75(15):5009–11.
Al Saif N, Brazier JS. The distribution of Clostridium difficile in the environment of South Wales. J Med Microbiol. 1996;45(2):133–7.
Zidaric V, Beigot S, Lapajne S, Rupnik M. The occurrence and high diversity of Clostridium difficile genotypes in rivers. Anaerobe. 2010;16(4):371–5.
Hargreaves KR, Colvin HV, Patel KV, Clokie JJP, Clokie MRJ. Genetically Diverse Clostridium difficile Strains Harboring Abundant Prophages in an Estuarine Environment. Appl Environ Microbiol. 2013;79(20):6236–43.
Del Mar Gamboa M, Rodriguez E, Vargas P. Diversity of mesophilic clostridia in Costa Rican soils. Anaerobe. 2005;11(6):322–6.
Pasquale V, Romano VJ, Rupnik M, Dumontet S, Ciznar I, Aliberti F, Mauri F, Saggiomo V, Krovacek K. Isolation and characterization of Clostridium difficile from shellfish and marine environments. Folia Microbiol. 2011;56(5):431–7.
Metcalf D, Avery BP, Janecko N, Matic N, Reid-Smith R, Weese JS. Clostridium difficile in seafood and fish. Anaerobe. 2011;17(2):85–6.
Miller MA, Byrne BA, Jang SS, Dodd EM, Dorfmeier E, Harris MD, Ames J, Paradies D, Worcester K, Jessup DA, et al. Enteric bacterial pathogen detection in southern sea otters (Enhydra lutris nereis) is associated with coastal urbanization and freshwater runoff. Vet Res. 2010;41(1):1.
Goorhuis A, Bakker D, Corver J, Debast S, Harmanus C, Notermans D, Bergwerff A, Dekker F, Kuijper E. Emergence of Clostridium difficile infection due to a new hypervirulent strain, polymerase chain reaction ribotype 078. Clin Infect Dis. 2008;47(9):1162–70.
Smits WK. Hype or hypervirulence. Virulence. 2013;4(7):592–6.
Barbut F, Rupnik M. 027, 078, and Others: Going Beyond the Numbers (and Away From the Hypervirulence). Clin Infect Dis. 2012;55(12):1669–72.
Cairns M, Stabler R, Shetty N, Wren B. The continually evolving Clostridium difficile species. Future Microbiol. 2012;7(8):945–57.
Sebaihia M, Wren B, Mullany P, Fairweather N, Minton N, Stabler R, Thomson N, Roberts A, Cerdeno-Tarrraga A, Wang H, et al. The multidrug-resistant human pathogen Clostridium difficile has a highly mobile, mosaic genome. Nat Genet. 2006;38(7):779–86.
Knetsch CW, Connor TR, Mutreja A, Van Dorp SM, Sanders IM, Browne HP, Harris D, Lipman L, Keessen EC, Corver J, et al. Whole genome sequencing reveals potential spread of Clostridium difficile between humans and farm animals in the Netherlands, 2002 to 2011. Eurosurveillance. 2014;19(45):30–41.
He M, Sebaihia M, Lawley T, Stabler R, Dawson L, Martin M, Holt K, Seth-Smith H, Quail M, Rance R, et al. Evolutionary dynamics of Clostridium difficile over short and long time scales. Proc Natl Acad Sci U S A. 2010;107(16):7527–32.
Kurka H, Ehrenreich A, Ludwig W, Monot M, Rupnik M, Barbut F, Indra A, Dupuy B, Liebl W. Sequence similarity of Clostridium difficile strains by analysis of conserved genes and genome content is reflected by their ribotype affiliation. PLoS One. 2014;9(1):e86535.
Hargreaves KR, Otieno JR, Thanki A, Blades MJ, Millard AD, Browne HP, Lawley TD, Clokie MRJ. As clear as mud? Determining the diverity and prevalance of prophages in the draft genomes of estuarine isolates of Clostridium difficile. Genome Biol Evol. 2015;27:evv094.
Knight DR, Elliott B, Chang BJ, Perkins TT, Riley TV. Diversity and Evolution in the Genome of Clostridium difficile. Clin Microbiol Rev. 2015;28(3):721–41.
Amy J, Johanesen P, Lyras D. Extrachromosomal and integrated genetic elements in Clostridium difficile. Plasmid. 2015;80:97–110.
Mullany P, Allan E, Roberts AP. Mobile genetic elements in Clostridium difficile and their role in genome function. Res Microbiol. 2015;166(4):361–7.
Vedantam G, Clark A, Chu M, McQuade R, Mallozzi M, Viswanathan VK. Clostridium difficile infection: toxins and non-toxin virulence factors, and their contributions to disease establishment and host response. Gut Microbes. 2012;3(2):121–34.
Luo Y, Huang C, Ye J, Fang W, Gu W, Chen Z, Li H, Wang X, Jin D. Genome Sequence and Analysis of Peptoclostridium difficile Strain ZJCDC-S82. Evol Bioinforma. 2016;12:41–9.
Van Eijk E, Anvar SY, Browne HP, Leung WY, Frank J, Schmitz AM, Roberts AP, Smits WK. Complete genome sequence of the Clostridium difficile laboratory strain 630Deltaerm reveals differences from strain 630, including translocation of the mobile element CTn5. BMC Genomics. 2015;16:31.
Gaulton T, Misra R, Rose G, Baybayan P, Hall R, Freeman J, Turton J, Picton S, Korlach J, Gharbia S, et al. Complete Genome Sequence of the Hypervirulent Bacterium Clostridium difficile Strain G46, Ribotype 027. Genome Announc. 2015;3(2):e00073–15.
Riedel T, Bunk B, Wittmann J, Thurmer A, Sproer C, Gronow S, Liesegang H, Daniel R, Overmann J. Complete Genome Sequence of the Clostridium difficile Type Strain DSM 1296T. Genome Announc. 2015;3(5):e01186–01115.
Riedel T, Bunk B, Thurmer A, Sproer C, Brzuszkiewicz E, Abt B, Gronow S, Liesegang H, Daniel R, Overmann J. Genome Resequencing of the Virulent and Multidrug-Resistant Reference Strain Clostridium difficile 630. Genome Announc. 2015;3(2):e00276–15.
Gyorfy Z, Draskovits G, Vernyik V, Blattner FF, Gaal T, Posfai G. Engineered ribosomal RNA operon copy-number variants of E-coli reveal the evolutionary trade-offs shaping rRNA operon number. Nucleic Acids Res. 2015;43(3):1783–94.
Janezic S, Ocepek M, Zidaric V, Rupnik M. Clostridium difficile genotypes other than ribotype 078 that are prevalent among human, animal and environmental isolates. BMC Microbiol. 2012;12:48.
Corver J, Bakker D, Brouwer MS, Harmanus C, Hensgens MP, Roberts AP, Lipman LJ, Kuijper EJ, Van Leeuwen HC. Analysis of a Clostridium difficile PCR ribotype 078 100 kilobase island reveals the presence of a novel transposon, Tn6164. BMC Microbiol. 2012;12:130.
Santagati M, Iannelli F, Cascone C, Campanile F, Oggioni MR, Stefani S, Pozzi G. The novel conjugative transposon tn1207.3 carries the macrolide efflux gene mef (A) in Streptococcus pyogenes. Microb Drug Resist. 2003;9(3):243–7.
Iannelli F, Santagati M, Santoro F, Oggioni MR, Stefani S, Pozzi G. Nucleotide sequence of conjugative prophage Φ1207.3 (formerly Tn1207.3) carrying the mef (A)/msr (D) genes for efflux resistance to macrolides in Streptococcus pyogenes. Front Microbiol. 2014;5:687.
Johnson SR, Romig WR. Vibrio cholerae conjugative plasmid pSJ15 contains transposable prophage dVcA1. J Bacteriol. 1981;146(2):632–8.
Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007;57:81–91.
Feiner R, Argov T, Rabinovich L, Sigal N, Borovok I, Herskovits AA. A new perspective on lysogeny: prophages as active regulatory switches of bacteria. Nat Rev Microbiol. 2015;13(10):641–50.
Griffiths D, Fawley W, Kachrimanidou M, Bowden R, Crook DW, Fung R, Golubchik T, Harding RM, Jeffery KJM, Jolley KA, et al. Multilocus Sequence Typing of Clostridium difficile. J Clin Microbiol. 2010;48(3):770–8.
Stabler RA, Dawson LF, Valiente E, Cairns MD, Martin MJ, Donahue EH, Riley TV, Songer JG, Kuijper EJ, Dingle KE, et al. Macro and micro diversity of Clostridium difficile isolates from diverse sources and geographical locations. PLoS One. 2012;7(3):e31559.
Brouwer MS, Warburton PJ, Roberts AP, Mullany P, Allan E. Genetic organisation, mobility and predicted functions of genes on integrated, mobile genetic elements in sequenced strains of Clostridium difficile. PLoS One. 2011;6(8):e23014.
Shaw KJ, Rather PN, Hare RS, Miller GH. Moleclar genetics of aminoglycoside resistance genes and familial relationships of the aminoglycoside-modifying enzymes. Microbiol Rev. 1993;57(1):138–63.
Butala M, Zgur-Bertok D, Busby SJW. The bacterial LexA transcriptional repressor. Cell Mol Life Sci. 2009;66(1):82–93.
Toh S-M, Xiong L, Bae T, Mankin AS. The methyltransferase YfgB/RlmN is responsible for modification of adenosine 2503 in 23S rRNA. RNA. 2008;14(1):98–106.
Courvalin P. Transfer of antibiotic-resistance genes between Gram-positive and Gram-negative bacteria. Antimicrob Agents Chemother. 1994;38(7):1447–51.
Nale J, Shan J, Hickenbotham P, Fawley W, Wilcox M, Clokie M. Diverse Temperate Bacteriophage Carriage in Clostridium difficile 027 Strains. PLoS One. 2012;7(5):e37263.
Fortier L, Moineau S. Morphological and genetic diversity of temperate phages in Clostridium difficile. Appl Environ Microbiol. 2007;73(22):7358–66.
Sell TL, Schaberg DR, Fekety FR. Bacteriophage and bacteriocin typing scheme for Clostridium difficile. J Clin Microbiol. 1983;17(6):1148–52.
Sangster W, Hegarty JP, Stewart DB. Phage tail-like particles kill Clostridium difficile and represent an alternative to conventional antibiotics. Surgery. 2015;157(1):96–103.
Gebhart D, Williams SR, Bishop-Lilly KA, Govoni GR, Willner KM, Butani A, Sozhamannan S, Martin D, Fortier L-C, Scholl D. Novel High-Molecular-Weight, R-Type Bacteriocins of Clostridium difficile. J Bacteriol. 2012;194(22):6240–7.
Gebhart D, Lok S, Clare S, Tomas M, Stares M, Scholl D, Donskey CJ, Lawley TD, Govoni GR. A modified R-type bacteriocin specifically targeting Clostridium difficile prevents colonization of mice without affecting gut microbiota diversity. MBio. 2015;6(2):e02368–14.
Hahnke S, Striesow J, Elvert M, Mollar XP, Klocke M. Clostridium bornimense sp. nov., isolated from a mesophilic, two-phase, laboratory-scale biogas reactor. Int J Syst Evol Microbiol. 2014;64(8):2792–7.
He M, Miyajima F, Roberts P, Ellison L, Pickard D, Martin M, Connor T, Harris S, Fairley D, Bamford K, et al. Emergence and global spread of epidemic healthcare-associated Clostridium difficile. Nat Genet. 2013;45(1):4.
Flusberg BA, Webster DR, Lee JH, Travers KJ, Olivares EC, Clark TA, Korlach J, Turner SW. Direct detection of DNA methylation during single-molecule, real-time sequencing. Nat Methods. 2010;7(6):461–U472.
Marinus MG, Casadesus J. Roles of DNA adenine methylation in host-pathogen interactions: mismatch repair, transcriptional regulation, and more. FEMS Microbiol Rev. 2009;33(3):488–503.
Oliveira PH, Touchon M, Rocha EPC. The interplay of restriction-modification systems with mobile genetic elements and their prokaryotic hosts. Nucleic Acids Res. 2014;42(16):10618–U10803.
Roberts RJ, Vincze T, Posfai J, Macelis D. REBASE-a database for DNA restriction and modification: enzymes, genes and genomes. Nucleic Acids Res. 2015;43(D1):D298–9.
Stabler RA, He M, Dawson L, Martin M, Valiente E, Corton C, Lawley TD, Sebaihia M, Quail MA, Rose G, et al. Comparative genome and phenotypic analysis of Clostridium difficile 027 strains provides insight into the evolution of a hypervirulent bacterium. Genome Biol. 2009;10(9):R102.
Kobayashi I. Behavior of restriction-modification systems as selfish mobile elements and their impact on genome evolution. Nucleic Acids Res. 2001;29(18):3742–56.
Thomas AT, Brammar WJ, Wilkins BM. Plasmid R16 ArdA protein preferentially targets restriction activity of the type I restriction-modification system EcoKI. J Bacteriol. 2003;185(6):2022–5.
Barrangou R, Fremaux C, Deveau H, Richards M, Boyaval P, Moineau S, Romero D, Horvath P. CRISPR provides acquired resistance against viruses in prokaryotes. Science. 2007;315(5819):1709–12.
Westra ER, Buckling A, Fineran PC. CRISPR-Cas systems: beyond adaptive immunity. Nat Rev Microbiol. 2014;12(5):317–26.
Makarova KS, Wolf YI, Alkhnbashi OS, Costa F, Shah SA, Saunders SJ, Barrangou R, Brouns SJJ, Charpentier E, Haft DH, et al. An updated evolutionary classification of CRISPR-Cas systems. Nat Rev Microbiol. 2015;13(11):722–36.
Marraffini LA, Sontheimer EJ. Self versus non-self discrimination during CRISPR RNA-directed immunity. Nature. 2010;463(7280):568–U194.
Hargreaves KR, Flores CO, Lawley TD, Clokie MR. Abundant and diverse clustered regularly interspaced short palindromic repeat spacers in Clostridium difficile strains and prophages target multiple phage types within this pathogen. MBio. 2014;5(5):e01045–01013.
Biswas A, Gagnon JN, Brouns SJJ, Fineran PC, Brown CM. CRISPRTarget: Bioinformatic prediction and analysis of crRNA targets. RNA Biol. 2013;10(5):817–27.
Roux S, Faubladier M, Mahul A, Paulhe N, Bernard A, Debroas D, Enault F. Metavir: a web server dedicated to virome analysis. Bioinformatics. 2011;27(21):3074–5.
Sampson TR, Weiss DS. CRISPR-Cas systems: new players in gene regulation and bacterial physiology. Front Cell Infect Microbiol. 2014;4:37.
Brown SD, Nagaraju S, Utturkar S, De Tissera S, Segovia S, Mitchell W, Land ML, Dassanayake A, Koepke M. Comparison of single-molecule sequencing and hybrid approaches for finishing the genome of Clostridium autoethanogenum and analysis of CRISPR systems in industrial relevant Clostridia. Biotechnology for Biofuels. 2014;7:40.
Dupuis M-E, Villion M, Magadan AH, Moineau S. CRISPR-Cas and restriction-modification systems are compatible and increase phage resistance. Nat Commun. 2013;4:2087.
Krumsiek J, Arnold R, Rattei T. Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics. 2007;23(8):1026–8.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream M, Barrell B. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.
Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26(5):589–95.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Garcia-Alcalde F, Okonechnikov K, Carbonell J, Cruz LM, Goetz S, Tarazona S, Dopazo J, Meyer TF, Conesa A. Qualimap: evaluating next-generation sequencing alignment data. Bioinformatics. 2012;28(20):2678–9.
Carver T, Thomson N, Bleasby A, Berriman M, Parkhill J. DNAPlotter: circular and linear interactive genome visualization. Bioinformatics. 2009;25(1):119–20.
Seamann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.
Lagesen K, Hallin P, Rodland EA, Staerfeldt H-H, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35(9):3100–8.
Bidet P, Barbut F, Lalande V, Burghoffer B, Petit JC. Development of a new PCR-ribotyping method for Clostridium difficile based on ribosomal RNA gene sequencing. FEMS Microbiol Lett. 1999;175(2):261–6.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST plus: architecture and applications. Bmc Bioinformatics. 2009;10:421.
Jolley KA, Maiden MCJ. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11:595.
Darling A, Mau B, Blattner F, Perna N. Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14(7):1394–403.
Alikhan N-F, Petty NK, Ben Zakour NL, Beatson SA. BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011;12:402.
Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: A Fast Phage Search Tool. Nucleic Acids Res. 2011;39:W347–52.
Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:W52–7.
Grissa I, Vergnaud G, Pourcel C. CRISPRcompar: a website to compare clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2008;36:W145–8.
We wish to thank Dr. Trevor Lawley, Wellcome Trust Sanger Institute, UK, for kindly providing the reference M120 strain used in this study. We would like to thank the Centre for Genomic Research, the University of Liverpool, UK, for their expertise in the sequencing and assembly of these genomes in particular Margaret A. Hughes and Xuan Liu. We also wish to thank Megan De Ste Croix, University of Leicester, UK, for her help and advice with high molecular weight DNA extraction methods and Michael Wieand, Pacific Biosciences, USA, for his generous and helpful advice on SMRT generated genomic data. This research used the ALICE High Performance Computing Facility at the University of Leicester. We also wish to thank the anonymous reviewers for their helpful comments on this manuscript.
This work has been funded by a NERC grant NBAF896.
Availability of data and materials
The two SMRT generated genomes of M120 and CD105HS27 have been deposited in the ENA in study PRJEB13565 and have accession numbers ERS1242840 and ERS1242839, respectively. The methylation data has been deposited for each isolate in REBASE as organisms #19242 and #19243, respectively.
KRH and MRJC conceived the study, AMT performed the DNA extractions, KRH, BRJ and MO performed bioinformatics analyses and interpretations. All authors contributed to the writing of the manuscript. All read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
In silico ribotype profiles for the C. difficile genomes. (DOCX 12 kb)
Novel transposon sequence similarity in C. difficile strains. (DOCX 12 kb)
CRISPR spacer matches to C. difficile genomic sequences. (DOCX 27 kb)
Genome sequences used in this study. (DOCX 13 kb)
Viral metagenome datasets used for protospacer identification. (DOCX 13 kb)