- Genome Announcement
- Open Access
The genome of Clostridium difficile 5.3
Gut Pathogens volume 6, Article number: 4 (2014)
Clostridium difficile is the leading cause of infectious diarrhea in humans and responsible for large outbreaks of enteritis in neonatal pigs in both North America and Europe. Disease caused by C. difficile typically occurs during antibiotic therapy and its emergence over the past 40 years is linked with the widespread use of broad-spectrum antibiotics in both human and veterinary medicine.
We sequenced the genome of Clostridium difficile 5.3 using the Illumina Nextera XT and MiSeq technologies. Assembly of the sequence data reconstructed a 4,009,318 bp genome in 27 scaffolds with an N50 of 786 kbp. The genome has extensive similarity to other sequenced C. difficile genomes, but also has several genes that are potentially related to virulence and pathogenicity that are not present in the reference C. difficile strain.
Genome sequencing of human and animal isolates is needed to understand the molecular events driving the emergence of C. difficile as a gastrointestinal pathogen of humans and food animals and to better define its zoonotic potential.
Clostridium difficile is a Gram-positive, spore-forming bacterium that over the last 40 years has emerged from obscurity to become a leading gastrointestinal pathogen in hospital environments . The clinical features of disease range from non-haemorrhagic watery diarrhea to fatal episodes of pseudomembranous colitis and toxic megacolon. C. difficile disease is also community-acquired with incidence rates increasing . In the past decade C. difficile has caused widespread outbreaks of neonatal diarrhea in piglets raising concerns that it may have a food-borne etiology. One porcine strain of C. difficile is frequently isolated from human cases of C. difficile disease in human in Europe, underscoring its zoonotic potential .
Non-toxigenic strains of Clostridium difficile have been isolated from patients with diarrhoeal disease . Recent studies showing that nontoxigenic strains are capable of blocking colonisation and prevent disease caused by toxigenic strains indicates that much remains to be learned about virulence factors that play important roles in colonization of the gastrointestinal tract . Here we report the sequence of nontoxigenic C. difficile strain 5.3 isolated from a patient with diarrhoea. The clinical presentation of this patient was of suspected Clostridium difficile infection. This was confirmed at the pathology laboratory with the isolation of only non-toxigenic C. difficile.
Isolation and DNA preparation
Clostridium difficile 5.3 originated from a faecal sample submitted from a female patient with the clinical symptomology of C. difficile infection. From the faecal sample 100 μ l of sample was added to cooked meat medium (TM0102 Oxoid Australia) and incubated anaerobically using the anoxomat system (MART Microbiology B.B., The Netherlands) at 37C for 24 hours. From here 200 μ l of the culture was transferred aseptically into a centrifuge tube and centrifuged at 10,000 rpm for 5 mins. The supernatant was discarded and the pellet resuspended in 1 ml of absolute ethanol and incubated at room temperature with periodic inversion for two hours.
A pellet was recovered by centrifuging at 10,000 rpm for 5 mins, the ethanol was discarded and the pellet resuspended in 100 μ l of brain heart infusion broth and spread plated onto Clostridium difficile selective agar (CC-BHIA + Taurocholate, PP2362 Oxoid Australia) and incubated anaerobically at 37 C for 24 hours. All plates were examined for colonies that morphologically represented C. difficile, from each plate colonies were selected and subcultured onto CC-BHIA + Taurocholate until pure cultures were achieved.
DNA was extracted from pure cultures by taking a single colony into 2 ml of brain heart infusion agar and incubating anaerobically at 37 C for upto 48 hours. The bacterial suspension was then vortexed and the DNA extracted using the DNeasy Blood and Tissue Kit (Qiagen, 69581) using the manufacturers instructions for the extraction of gram positive bacteria and the use of lysozyme.
DNA was quantified using qubit flourimetry and 0.5 ng gDNA was used as input to the Illumina Nextera XT library preparation protocol. Tagmentation of Genomic DNA, and PCR amplification of tagged DNA were performed as per manufacturer’s instructions. However the “PCR Clean-Up” and “Library Normalization” steps were omitted and size selection was instead performed by running balanced and pooled samples in a one percent agarose gel and excising the 600 bp to 1200 bp region of interest. The DNA was then purified from the agarose using Promega’s Wizard SV Gel and PCR Clean-Up System. Finally, an Agilent 2100 Bioanalyzer, with a High Sensitivity DNA Kit, was used to quantitate the pooled DNA library before loading onto the MiSeq with 23 other multiplexed samples. Paired-end 300 nt reads were generated using MiSeq V3 chemistry.
Assembly and annotation
The genome was assembled using a version of the A5 pipeline  called A5-miseq  that has been revised to process reads up to 500 nt long. Briefly, A5-miseq consists of five stages: (1) read quality filtering and error correction, (2) contig assembly, (3) permissive draft scaffolding, (4) misassembly detection, and (5) conservative scaffolding. The revised A5 pipeline uses a new version of idba_ud that uses read pairing information, and that has been modified to accept reads up to 500 nt long and to construct de Bruijn graphs with k-mers up to 500 nt. These modifications provide substantial improvements in assembly contiguity.
The genome was annotated with the RAST annotation system using FigFAM release 70 .
A5-miseq includes a quality checking step that detects putative misassemblies by identifying clusters of read pairs that map to disjoint locations in the assembled genome. This method did not detect any putative misassemblies.
Sequencing generated 925,170 read pairs for a total of 555,102,000 nt that were assembled to reconstruct the 4,009,318 bp genome of C. difficile 5.3 in 27 scaffolds, with a scaffold N50 of 786 kbp and an N90 of 135 kbp. The raw (unfiltered) coverage is 138x, and after read filtering the assembly has a median depth of coverage of 60. The annotation of this assembly identified 3647 predicted CDS and 125 predicted RNA genes. 13 genes were identified as possibly missing from the assembly by the RAST system. The overall functional profile of the genome is shown in Figure 1. We conducted a phylogenetic analysis of C. difficile 5.3 using the PhyloSift software  to identify the most closely related organism with an available reference genome, using a phylogenetic placement  on a tree constructed from a set of 37 universally conserved genes. The resulting analysis assigned a posterior probability of 100% for C. difficile 5.3 diverging on the same lineage as Clostridium difficile 6503. Because Clostridium difficile 6503 has a draft-quality genome we used the closely related genome of Clostridium difficile 630 as a reference for further comparative analysis.
The scaffolds of C. difficile 5.3 were reordered to match the order in the finished genome of the closely related strain C. difficile 630 using the Mauve Contig Mover . After reordering we find the genomes are free from large-scale rearrangement. Rearrangement breakpoints for some small translocated regions occur at sites annotated with phage-related gene functions, suggesting that bacteriophage may have lysogenized at different locations in each strain. The genome of C. difficile 5.3 was aligned to those of C. difficile BI1, C. difficile 2007855, and C. difficile 630 to characterize the extent of gene content sharing and genomic rearrangement. A visualization of the genome comparison produced by the Mauve software is shown in Figure 2. From this figure we can see that the majority of the genome is conserved among all four isolates, although some strain-specific regions exist.
Comparison of the gene content between C. difficile 5.3 and the finished C. difficile 630 reference genome identified 105 annotated gene functions predicted to be present only in C. difficile 5.3. In particular, genes related to metal resistance, bacitracin resistance, colicin resistance, drug efflux, and several phage-related genes are present only in C. difficile 5.3. These genes may contribute to the pathogenic phenotype of C. difficile 5.3. A full list of these genes has been provided in Additional file 1.
Improved efficiency of the clinical genomics pipeline will enable fine-scale epidemiological monitoring of Clostridium difficile outbreaks.
Availability of supporting data
The draft genome assembly has been submitted to NCBI and assigned accession PRJNA232267. Genome annotations are available from the RAST web server under accession 6666666.54620 when logged in with username guest, password guest. The Illumina sequence reads have been deposited to the Short Read Archive and under accession SRX396630.
Rapid annotation using subsystem technology
Andrew and Aaron’s Awesome Assembly
Squire MM, Riley TV: One Health: The human-animal-environment interfaces in emerging infectious diseases. Curr Top Microbiol Immunol. 2013, 365: 299-314.http://link.springer.com/chapter/10.1007/82_2012_237,
Khanna S, Pardi DS, Aronson SL, Kammer PP, Orenstein R, St Sauver JL, Harmsen WS, Zinsmeister AR: The epidemiology of community-acquired clostridium difficile infection: a population-based study. Am J Gastroenterol. 2011, doi:10.1038/ajg.2011.398
Goorhuis A, Bakker D, Corver J, Debast SB, Harmanus C, Notermans DW, Bergwerff AA, Dekker FW, Kuijper EJ: Emergence of clostridium difficile infection due to a new hypervirulent strain, polymerase chain reaction ribotype 078. Clin Infect Dis. 2008, 47 (9): 1162-1170. 10.1086/592257.
Martirosian G, Szczesny A, Cohen SH, Silva Jr J: Isolation of non-toxigenic strains of clostridium difficile from cases of diarrhea among patients hospitalized in hematology/oncology ward. Polish J Microbiol. 2004, 53 (3): 197-200.
Nagaro KJ, Phillips ST, Cheknis AK, Sambol SP, Zukowski WE, Johnson S, Gerding DN: 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-5270. 10.1128/AAC.00580-13.
Tritt A, Eisen JA, Facciotti MT, Darling AE: An integrated pipeline for de novo assembly of microbial genomes. PLoS ONE. 2012, 7 (9): 42304-10.1371/journal.pone.0042304. doi:10.1371/journal.pone.0042304
Coil D, Jospin G, Darling AE: A5-miseq: an updated pipeline to assemble microbial genomes from illumina miseq data. arXiv preprint arXiv:1401.5130. 2014,http://arxiv.org/abs/1401.5130,
Aziz R, Bartels D, Best A, DeJongh M, Disz T, Edwards R, Formsma K, Gerdes S, Glass E, Kubal M, Meyer F, Olsen G, Olson R, Osterman A, Overbeek R, McNeil L, Paarmann D, Paczian T, Parrello B, Pusch G, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST server: Rapid Annotations using Subsystems Technology. BMC Genomics. 2008, 9 (1): 75-10.1186/1471-2164-9-75. doi:10.1186/1471-2164-9-75
Darling AE, Jospin G, Lowe E, Bik HM, Eisen JA, Matsen IV FA: Phylosift: phylogenetic analysis of genomes and metagenomes. PeerJ. 2014, 2: 243-
Matsen F, Kodner R, Armbrust EV: pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree. BMC Bioinformatics. 2010, 11 (1): 538-10.1186/1471-2105-11-538. doi:10.1186/1471-2105-11-538
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-2073. 10.1093/bioinformatics/btp356.
This work was supported by a collaboration between the NSW Department of Primary Industries and the ithree institute.
The authors declare that they have no competing interests.
TAC provided the DNA. PW constructed Illumina libraries and sequenced them. AED conducted the assembly, analysis, and data deposition. AED, TAC, SPD, and PW wrote the paper. IGC and PRC provided general direction. All authors read and approved the final manuscript.