Skip to main content

Complete genome sequence of Clostridium perfringens CBA7123 isolated from a faecal sample from Korea

Abstract

Background

Clostridium perfringens is an opportunistic human pathogen that causes necrotic enteritis, mild diarrhea, clostridial myonecrosis or gas gangrene, sepsis, etc. In this study, we aim to determine the pathogenesis of this bacterium at the genomic level. The genome of strain CBA7123 was sequenced, and a comparative genomic analysis between strain CBA7123 and four other related C. perfringens strains was performed.

Results

The genome of strain CBA7123 consisted of one circular chromosome and one plasmid that were 3,088,370 and 46,640 bp long with 28.5 and 27.1 mol% G+C content, respectively. The genomic DNA was predicted to contain 2798 open reading frames (ORFs), 10 rRNA genes, and 94 tRNA genes. The genomic comparison analysis between the five strains revealed the distinctive virulence properties of strain CBA7123 by highlighting certain strain-specific genes.

Conclusions

In this study, the C. perfringens CBA7123 genome was sequenced and compared with other C. perfringens genomes. Among the various genes sequenced, the detection of antimicrobial resistance genes and those encoding various virulence factors may extend the understanding of the pathogenesis of C. perfringens strains.

Background

Clostridium perfringens is a Gram-positive, spore-forming, strict anaerobic, rod-shaped bacterium belonging to the phylum Firmicutes [1,2,3]. This bacterium inhabits diverse environments such as soil, sewage, and animal intestines [4]. Although C. perfringens does not invade healthy cells, it acts as a pathogen by producing various enzymes and toxins and is also considered a common cause of food poisoning worldwide [5]. C. perfringens toxins, mainly comprising α, β, ε, and ι extracellular toxins, are indicators for classifying its strains as A to E toxinotypes [6,7,8]. The type A of C. perfringens is the most commonly found toxinotype that possesses only α toxin-encoding gene, plc. Strains with α-toxin genes can cause necrotic enteritis, mild diarrhea, clostridial myonecrosis or gas gangrene, sepsis, and food poisoning in humans as well as various enterotoxemic diseases in livestock [7, 9, 10]. In this study, we aim to elucidate the pathogenesis of this species through its genome, specifically the virulence-related genes. The genome of C. perfringens CBA7123 was completely sequenced and analyzed using bioinformatics. Additionally, genomic data that were compared between strain CBA7123 and four other C. perfringens strains would illustrate the virulence mechanisms of these bacteria.

Methods

Isolation and DNA extraction of C. perfringens CBA7123

Clostridium perfringens CBA7123 (=KCCM 43242) was isolated from the feces of a 73-year-old man, and a pure culture was obtained using serial dilution. The colony was cultured anaerobically in ATCC medium no. 2840, modified Eggerth–Gagnon medium (10 g peptone, 4 g Na2HPO4·2H2O, 2 g porcine gastric mucin, 50 ml sheep blood, and 15 g agar per liter) at 37 °C for 24 h. The total genomic DNA of strain CBA7123 was extracted using QuickGene DNA tissue kit S (Kurabo, Japan) and the G-spin total DNA extraction kit (iNtRON Biotechnology, Korea). The purity, quality, and quantity of genomic DNA were measured using Agilent 2100 Bioanalyzer (Agilent Technologies, USA) as described the manufacturer’s instruction.

Library preparation and genome sequencing

Detailed genome sequencing was performed as described previously [11]. Briefly, the genomic DNA of strain CBA7123 was sheared according to the PacBio 20-kb Template Preparation using BluePippin Size-Selection System protocol, and SMRTbell library was prepared using P6-C4 chemistry (Pacific Biosciences, USA). The sequences were obtained using PacBio RS II system (Pacific Biosciences) as following an instruction of the manufacturer. The 150,292 reads were generated with 7090 bp of average read length by the PacBio RS II system from one SMRT cell.

Genome assembly and annotation

De novo assembly of the genome sequence was performed using Hierarchical Genome Assembly Process (HGAP) version 2 software, with default parameters supported by PacBio SMRT Analysis ver. 2.3.0 [12]. rRNA and tRNA of the assembled sequence were identified using RNAmmer 1.2 and tRNAscan-SE 1.21, respectively. Genes were predicted using Glimmer3 of Rapid Annotation using Subsystem Technology (RAST) server (http://rast.nmpdr.org), and functional gene annotations were performed using the SEED, Clusters of Orthologous Groups (COG, http://www.ncbi.nlm.nih.gov/COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) databases. PathogenFinder 1.1 [13] and ResFinder 2.1 [14] were used to estimate the pathogenicity and antimicrobial resistance genes, respectively. The virulence factors were searched using the Basic Local Alignment Search Tool (BLAST) in the virulence factors of pathogenic bacteria database [15] with default parameters and predicted with zero e-value.

Comparative genomic analysis

The C. perfringens strains for comparative genomic analysis were selected using Microbial Nucleotide BLAST in the NCBI complete genome database. The four strains with BLAST total scores over 5.5e+06 were selected: strains FORC 003, JP55, FORC 025, and JP838. To compare the genomic structures between strain CBA7123 and these four strains, the progressive alignment algorithm in MAUVE multiple genome alignment software 2.4.0 was used [16]. The OrthoANI algorithm was used to analyze genomic relatedness between strain CBA7123 and other C. perfringens species. OrthoANI percentages were calculated and a phylogenetic tree was constructed, as described by Lee et al. [17]. Orthologs between strain CBA7123 and the reference strains were predicted and mapped using the reciprocal best hit method in UBLAST [18]. Pan-genome orthologous groups (POGs) were estimated using the EzBioCloud Comparative Genomics Database (http://cg.ezbiocloud.net/) [19], and their presence was calculated using the Jaccard coefficient. UPGMA clustering was then used to show clustering between strain CBA7123 and the reference strains as a dendrogram, based on the presence or absence of gene content. The intersections and differences of POG sets of the different strains were visualized as a Venn diagram, using the jvenn program [20].

Quality assurance

The single colony of strain CBA7123 was transferred more than three times to obtain a pure culture and confirmed using a variable pressure field emission scanning electron microscope (Chuncheon Center, Korea Basic Science Institute, Korea) (Additional file 1: Figure S1) and a PCR-sequencing approach targeting 16S rRNA gene as described by Min et al. [21].

Results and discussion

General features

The genome of strain CBA7123 comprised two circular contigs that were 3,088,370 and 46,640 bp long with 28.5 and 27.1 mol% G+C content, respectively. Each contig matched to a chromosome and plasmid by the BLAST analysis with GenBank genome database. 2798 open reading frames (ORFs), 10 rRNA (10 5S rRNA, 10 16S rRNA, 10 23S rRNA) operons, and 94 tRNA genes were present. Among the various COG categories, G (carbohydrate transport and metabolism; 205 ORFs), E (amino acid transport and metabolism; 168 ORFs), J (translation, ribosomal structure, and biogenesis; 163 ORFs), K (transcription; 149 ORFs), and L (replication, recombination, and repair; 147 ORFs) made up the largest proportion (≥6% of the total COGs classifications). Detailed genome features are summarized in Table 1, and numbers of COG functional categories are shown in Table 2.

Table 1 General features of the Clostridium perfringens CBA7123 genome
Table 2 Number of general COG-associated functional genes

Comparative genomic analysis of strain CBA7123 with other C. perfringens strains

In the comparison of genomic structures between strain CBA7123 and strains FORC 003, JP55, FORC 025, and JP838, locally collinear blocks indicated high homology (Additional file 1: Figure S2). OrthoANI values between strain CBA7123 and each of strains FORC 003, JP55, FORC 025, and JP838 were 96.2, 96.2, 96.3, and 95.8, respectively. Although the phylogenetic tree constructed based on these OrthoANI values did not obviously indicate that strain CBA7123 was clustered with the other C. perfringens strains (Fig. 1a), all of the OrthoANI values were over 95%, which is the cut-off for species demarcation [17]. A POG comparison showed that the five strains shared 2392 POGs; however, 152 POGs were exclusive to strain CBA7123 (Fig. 2). Of the 152 POGs, some POGs encoded for antibiotic resistance and antiviral defense functions: two kinds of multidrug export MepA proteins, CRISPR-associated endonuclease Cas1 and Cas2, 5′ to 3′ exodeoxyribonuclease. The dendrogram created based on the presence of POGs indicated that strain CBA7123 was clustered with strains JP55, FORC 003, and FORC 025 (Fig. 1b). These results revealed that strain CBA7123 is closely related to C. perfringens, but different from other C. perfringens strains.

Fig. 1
figure 1

Phylogenetic tree depicting the relationship between Clostridium perfringens CBA7123 and four reference strains: FORC 003, FORC 025, JP55, and JP838, constructed based on a OrthoANI values using the orthologous average nucleotide identity tool and b gene content (presence or absence) using Jaccard coefficients and UPGMA clustering

Fig. 2
figure 2

Venn diagram based on POGs indicating the orthologous groups among five Clostridium perfringens strains. Number of Venn diagrams represents number of shared genes among the genomes of the five strains

Pathogenesis, virulence factors, and antimicrobial resistance genes

Strain CBA7123 was determined to be a human pathogen with 0.77 probability, using PathogenFinder 1.1. The 13 proteins except for conserved hypothetical proteins and hypothetical proteins were matched with pathogenic organisms belonging to class Clostridia (Table 3). The genome contained tetracycline-resistance genes: two tetA(P) genes, one each in the chromosome and plasmid; and one tetB(P) gene in the plasmid. Strain FORC 003 genome also contained the same tetracycline-resistance genes, whereas strain JP838 genome only contained tetA(P). No antimicrobial resistance genes were detected in the other two strains. The key virulence factors of strain CBA7123 were predicted to be exo-alpha-sialidase (nanI), sialidase (nanH), collagenase (colA, kappa-toxin), perfringolysin O (pfoA, theta-toxin), phospholipase C (plc, alpha-toxin), and alpha-clostripain (cloSI) with 96, 95, 98, 98 and 96% identity, respectively, in the chromosome and all of which were also present in the four reference strains. Sialidase hydrolyzes the α-linkage of terminal sialic acid residues from mammalian cell surface to generate free sialic acid, which is important carbon and energy source for colonization and growth of bacterial pathogen in the intestines [22]. Perfringolysin O and phospholipase C are soluble toxin and phospholipid cleaving enzyme leading to host cell lysis, respectively, and these two enzymes have synergistic effects in Clostridium perfringens-mediated gas gangrene [23, 24]. Collagenase and alpha-clostripain are not major determinant of virulence in clostridial myonecrosis, but it could be assumed that the two proteins have effects on C. perfringens-mediated disease [25, 26]. Strain CBA7123 could be classified as toxinotype A, based on the presence of the α-toxin gene (plc). Identification of tetracycline-resistance genes and various other virulence factors may help future research into reducing the pathogenesis of C. perfringens in humans and contribute to faster treatment response.

Table 3 Pathogenesis-related proteins matched to the genus Clostridium

Future directions

The genome information of strain CBA7123 can improve the understanding of C. perfringens, and the information of its antimicrobial resistance and virulence factors can contribute to the development of methods for preventing C. perfringens-related food poisoning. Future research studies should investigate the pathogenesis mechanism in detail as well as the specific roles of each virulence factor.

Abbreviations

BLAST:

Basic Local Alignment Search Tool

COGs:

Clusters of Orthologous Groups

HGAP:

Hierarchical Genome Assembly Process

KEGG:

Kyoto Encyclopedia of Genes and Genomes

ORFs:

open reading frames

POGs:

pan-genome orthologous groups

RAST:

Rapid Annotation using Subsystem Technology

References

  1. Minton NP, Ehsaan M, Humphreys CM, Little GT, Baker J, Henstra AM, et al. A roadmap for gene system development in Clostridium. Anaerobe. 2016;41:104–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Shimizu T, Ohtani K, Hirakawa H, Ohshima K, Yamashita A, Shiba T, et al. Complete genome sequence of Clostridium perfringens, an anaerobic flesh-eater. Proc Natl Acad Sci USA. 2002;99:996–1001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Wong YM, Juan JC, Gan HM, Austin CM. Draft genome sequence of Clostridium perfringens strain JJC, a highly efficient hydrogen producer isolated from landfill leachate sludge. Genome Announc. 2014;2:e00064.

    PubMed  PubMed Central  Google Scholar 

  4. Johansson A, Aspan A, Bagge E, Baverud V, Engstrom BE, Johansson KE. Genetic diversity of Clostridium perfringens type A isolates from animals, food poisoning outbreaks and sludge. BMC Microbiol. 2006;6:47.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Gross TP, Kamara LB, Hatheway CL, Powers P, Libonati JP, Harmon SM, Israel E. Clostridium perfringens food poisoning: use of serotyping in an outbreak setting. J Clin Microbiol. 1989;27:660–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Hassan KA, Elbourne LD, Tetu SG, Melville SB, Rood JI, Paulsen IT. Genomic analyses of Clostridium perfringens isolates from five toxinotypes. Res Microbiol. 2015;166:255–63.

    Article  CAS  PubMed  Google Scholar 

  7. Petit L, Gibert M, Popoff MR. Clostridium perfringens: toxinotype and genotype. Trends Microbiol. 1999;7:104–10.

    Article  CAS  PubMed  Google Scholar 

  8. Rood JI, Cole ST. Molecular genetics and pathogenesis of Clostridium perfringens. Microbiol Rev. 1991;55:621–48.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Brynestad S, Granum PE. Clostridium perfringens and foodborne infections. Int J Food Microbiol. 2002;74:195–202.

    Article  PubMed  Google Scholar 

  10. Hashiba M, Tomino A, Takenaka N, Hattori T, Kano H, Tsuda M, Takeyama N. Clostridium Perfringens Infection in a Febrile Patient with Severe Hemolytic Anemia. Am J Case Rep. 2016;17:219–23.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Kim JY, Song HS, Kim YB, Kwon J, Choi JS, Cho YJ, et al. Genome sequence of a commensal bacterium, Enterococcus faecalis CBA7120, isolated from a Korean fecal sample. Gut Pathog. 2016;8:62.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10:563–9.

    Article  CAS  PubMed  Google Scholar 

  13. Cosentino S, Voldby Larsen M, Moller Aarestrup F, Lund O. PathogenFinder–distinguishing friend from foe using bacterial whole genome sequence data. PLoS ONE. 2013;8:e77302.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67:2640–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y, Jin Q. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2005;33:D325–8.

    Article  CAS  PubMed  Google Scholar 

  16. Darling AC, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14:1394–403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Lee I, Kim YO, Park SC, Chun J. OrthoANI: an improved algorithm and software for calculating average nucleotide identity. Int J Syst Evol Microbiol. 2015;66:1100–3.

    Article  Google Scholar 

  18. Ward N, Moreno-Hagelsieb G. Quickly finding orthologs as reciprocal best hits with BLAT, LAST, and UBLAST: how much do we miss? PLoS ONE. 2014;9:e101850.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Yoon SH, Ha SM, Kwon S, Lim J, Kim Y, Seo H, Chun J. Introducing EzBioCloud: a taxonomically united database of 16S rRNA and whole genome assemblies. Int J Syst Evolut Microbiol. 2016. (In press).

  20. Bardou P, Mariette J, Escudie F, Djemiel C, Klopp C. jvenn: an interactive Venn diagram viewer. BMC Bioinform. 2014;15:293.

    Article  Google Scholar 

  21. Min UG, Kim SJ, Hong H, Kim SG, Gwak JH, Jung MY, et al. Calculibacillus koreensis gen. nov., sp. nov., an anaerobic Fe(III)-reducing bacterium isolated from sediment of mine tailings. J Microbiol. 2016;54:413–9.

    Article  PubMed  Google Scholar 

  22. Li J, Uzal FA, McClane BA. Clostridium perfringens Sialidases: potential contributors to intestinal pathogenesis and therapeutic targets. Toxins. 2016;8:341.

    Article  PubMed Central  Google Scholar 

  23. Awad MM, Ellemor DM, Boyd RL, Emmins JJ, Rood JI. Synergistic effects of alpha-toxin and perfringolysin O in Clostridium perfringens-mediated gas gangrene. Infect Immun. 2001;69:7904–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Dang TX, Hotze EM, Rouiller I, Tweten RK, Wilson-Kubalek EM. Prepore to pore transition of a cholesterol-dependent cytolysin visualized by electron microscopy. J Struct Biol. 2005;150:100–8.

    Article  CAS  PubMed  Google Scholar 

  25. Awad MM, Ellemor DM, Bryant AE, Matsushita O, Boyd RL, Stevens DL, et al. Construction and virulence testing of a collagenase mutant of Clostridium perfringens. Microb Pathog. 2000;28:107–17.

    Article  CAS  PubMed  Google Scholar 

  26. Chakravorty A, Awad MM, Hiscox TJ, Cheung JK, Carter GP, Choo JM, et al. The cysteine protease alpha-clostripain is not essential for the pathogenesis of Clostridium perfringens-mediated myonecrosis. PLoS ONE. 2011;6:e22762.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

SWR and YDN designed and coordinated all the experiments. HSS performed the microbiology experiments. YBK, JYK, and SWR performed the bioinformatic analysis and wrote the manuscript. CL, JK, JK, JKR, and MSJ checked and edited the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

This complete genome project has been deposited at DDBJ under the accession numbers AP017630–AP017631.

Ethics approval and consent to participate

The study protocol was approved by the institutional review board of the Theragen ETEX Bio Institute (700062-20160804-JR-005-02).

Funding

This research was supported by the Korea Food Research Institute (E0170602-01), the World Institute of Kimchi funded by the Ministry of Science, ICT and Future Planning, Republic of Korea (KE1702-2), the Center for Analytical Research of Disaster Science of the Korea Basic Science Institute (C37703), and the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science, and Technology (2015R1D1A1A09061039).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Seong Woon Roh.

Additional file

13099_2017_181_MOESM1_ESM.pdf

Additional file 1: Figure S1. A photomicrograph of Clostridium perfringens strain CBA7123 using Variable Pressure Field Emission Scanning Electron Microscope (VP-FE-SEM). Figure S2. Comparison of genomic structure between Clostridium perfringens CBA7123 and strains FORC 003, FORC 025, JP55, and JP838, using a progressive alignment algorithm in Mauve. The locally collinear blocks with identical colors represent highly homologous regions. The genomes were figured based on scale of the genome of strain CBA7123.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kim, Y.B., Kim, J.Y., Song, H.S. et al. Complete genome sequence of Clostridium perfringens CBA7123 isolated from a faecal sample from Korea. Gut Pathog 9, 32 (2017). https://doi.org/10.1186/s13099-017-0181-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13099-017-0181-1

Keywords