The complete genome sequence of Cronobacter sakazakii ATCC 29544T, a food-borne pathogen, isolated from a child’s throat
© The Author(s) 2017
Received: 11 November 2016
Accepted: 7 December 2016
Published: 4 January 2017
Cronobacter sakazakii is an emerging opportunistic pathogen that is associated with rare but life-threatening cases of severe diseases: meningitis, necrotizing enterocolitis, and sepsis in premature and full-term infants. However, the pathogenesis mechanism of this pathogen remains largely unknown. To determine its pathogenesis at the genomic level, the genome of C. sakazakii ATCC 29544T was completely sequenced and analyzed.
The genomic DNA, containing a circular chromosome and three plasmids, is composed of 4,511,265 bp with a GC content of 56.71%, containing 4380 predicted open reading frames (ORFs), 22 rRNA genes, and 83 tRNA genes. The plasmids, designated pCSK29544_p1, pCSK29544_p2, and pCSK29544_p3, were 93,905-bp, 4938-bp, and 53,457-bp with GC contents of 57.02, 54.88, and 50.07%, respectively. They were also predicted to have 72, 6, and 57 ORFs without RNA genes.
The strain ATCC 29544T genome has ompA and ibeB-homologous cusC genes, probably associated with the invasion of human brain microvascular endothelial cells (BMECs). In addition, gene clusters for siderophore production (iucABCD/iutA) and the related transport system (eitCBAD) were detected in pCSK29544_p1 plasmid, indicating better iron uptake ability for survival. Furthermore, to survive under extremely dry condition like milk powder, this genome has gene clusters for biosynthesis of capsular proteins (CSK29544_00281-00284) and cellulose (CSK29544_01124-01127) for biofilm formation and a gene cluster for utilization of sialic acid in the milk (nanKTAR). The genome information of C. sakazakii ATCC 29544T would provide further understanding of its pathogenesis at the molecular level for the regulation of pathogenicity and the development of a rapid detection method using biomarkers.
KeywordsCronobacter sakazakii Complete genome sequence Pathogenesis Infant milk formula Virulence factor
Enterobacter sakazakii has been reclassified into seven species in the genus Cronobacter according to biochemical and genetic evaluations [1, 2]. Among them, Cronobacter sakazakii is a well-known opportunistic food-borne pathogen causing bacteremia, meningitis and necrotizing enterocolitis, particularly in low-birth-weight neonatal infants. This species is a Gram-negative, rod-shaped, peritrichous and yellow-pigmented facultative anaerobe belonging to the Enterobacteriaceae family . C. sakazakii has been often found in human and infant gut microbiota [4, 5]. Although C. sakazakii food-borne outbreaks are quite low, the fatality to infants generally ranges from 40 to 80% . Interestingly, C. sakazakii was reported to produce capsular or biofilm materials for its own protection from extremely dry conditions, as in formula milk powder, substantiating the high survival ability of C. sakazakii in the milk powder . After human infection, C. sakazakii can invade the intestinal epithelial cells and even the brain microvascular endothelial cells (BMECs), demonstrating its potentials to cause meningitis . Therefore, the biocontrol and regulation of C. sakazakii are urgently required. However, C. sakazakii is resistant to some antibiotics, indicating a problem with antibiotic therapies against C. sakazakii [9, 10], and its pathogenicity mechanism remains unknown. Recently, to unveil the knowledge about the Cronobacter ecology, multilocus sequence typing (MLST) analysis using seven housekeeping genes has been established to identify the diversity of the Cronobacter genus from various sources, and its application has facilitated understanding of the evolutionary relationships and environmental fitness of Cronobacter species . In this study, to understand its infection and pathogenesis at the molecular level, the genome of a representative C. sakazakii type strain, ATCC 29544T, was completely sequenced and analyzed in this study using bioinformatics. This genome information would provide the researchers with genomic insights into the virulence and pathogenicity mechanisms of this species for the further development of a rapid detection method and a novel biocontrol strategy.
Growth conditions and DNA isolation
Cronobacter sakazakii ATCC 29544T was routinely cultivated using Luria-Bertani (LB) medium at 37 °C with shaking at 220 rpm. Bacterial cells were harvested in the mid-exponential growth phase using centrifugation at 16,000×g for 1 min and its genomic DNA was isolated using G-spin™ Genomic DNA Extraction Kit for Bacteria (iNtRON Biotechnology, Seongnam, South Korea). The concentration and purity of extracted DNA were determined by NanoVue (GE healthcare, Little Chalfont, United Kingdom).
Genome sequencing and assembly
The complete genome of C. sakazakii ATCC 29544T was sequenced at Macrogen, Seoul, South Korea, using a hybrid of PacBio RS II (Pacific Biosciences, Menlo Park, CA, USA) and Illumina HiSeq 2500 (San Diego, CA, USA). The sequence reads from PacBio RS II and Illumina HiSeq 2500 platforms were assembled using HGAP (version 2.0) and ALLPATHS-LG (version r47449), respectively. The final genome coverages were on average 1321 X Illumina and 73 X PacBio, respectively.
The ORFs were predicted using Glimmer3  and GeneMark.hmm . The gene prediction results were confirmed by manual curation. The genes of rRNA and tRNA were predicted using RNAmmer 1.2  and tRNAscan-SE , respectively. The genome annotation was conducted using NCBI BLASTP  and a predicted protein analysis using InterProScan 5  for the prediction of protein functions.
The analysis of comparative genomes and phage-associated regions
Phage-associated gene clusters in the genome sequences of C. sakazakii ATCC 29544 were searched using PHAST server .
Cronobacter sakazakii ATCC 29544T was obtained from American Type Culture Collection (ATCC) and its morphological observation using a transmission electron microscopy (TEM) showed a traditional shape of C. sakazakii as a short rod (Additional file 1: Figure S1). In addition, this strain was confirmed to C. sakazakii using 16S rRNA gene sequencing. For genome sequencing, the raw read sequences were selected and assembled when their quality scores were more than 40 as cutoffs. The complete genome sequence after genome assembly was also used for confirmation using ANI analysis with previously reported complete genome sequences of C. sakazakii.
General genome properties
Pathogenesis and virulence factors
Cronobacter sakazakii infects human neonates and infants via mostly contaminated reconstituted infant formula milk, causing serious human diseases, including bacteremia, necrotizing enterocolitis, and even meningitis . The high survival rate of C. sakazakii, even under extremely dried conditions, as in milk formula powder, has not been investigated at the molecular level. To understand this extraordinary property of C. sakazakii, molecular studies have been recently performed. Interestingly, the gene clusters associated with biosynthesis of capsular polysaccharides (CSK29544_00281-00284) and cellulose (CSK29544_01124-01127) were detected. Recently, the biofilm formation and cellulose (as a component of biofilm) production of C. sakazakii were experimentally confirmed [18, 19], suggesting that these gene clusters may be involved in the biofilm formation. This biofilm formation of C. sakazakii may contribute to the survival in the infant formula conditions . However, exopolysaccharide (EPS) production was not observed in C. sakazakii ATCC 29544T , suggesting that capsular polysaccharide gene cluster may be associated with biofilm formation, not with EPS production. In particular, C. sakazakii is the only Cronobacter species that has the nanKTAR gene cluster to utilize sialic acid . Interestingly, C. sakazakii ATCC 29544T has this gene cluster (CSK29544_00587-590), indicating its ability to use it. This unique ability may be involved in its adaptation to the milk conditions because sialic acid is one of the components of milk .
Upon human infection, C. sakazakii invades brain microvascular endothelial cells (BMECs), barriers to protect the brain from infection of meningitic pathogens, including Escherichia coli K1 and Neisseria meningitides, which cause meningitis in neonates and infants [23, 24]. Outer membrane protein A (OmpA) of C. sakazakii is a critical determinant, contributing in vitro invasion of human BMECs by enhancing cell adhesion . In addition, a few genes (ibeA, ibeB, and yijP) in meningitic E. coli were also suggested to be associated with the invasion of human BMECs [25–27]. Interestingly, two genes, ompA (CSK29544_03699) for BMEC adhesion and the ibeB-homologous cusC (pCSK29544_3p0028) for the penetration of BMECs, were detected in the genome of C. sakazakii ATCC 29544T−, suggesting this strain may invade human BMECs. However, ibeA and yijP were not detected in the genome. It is noteworthy that this ibeB-homologous cusC is located in the gene cluster of the complete copper- and silver-resistance cation efflux system as encoded by cusABCF (CSK29544_3p0026–CSK29544_3p0031) in the plasmid pCSK29544_3p, suggesting this BMEC invasion ability may be transferrable to other C. sakazakii or that the strain ATCC 29544T may have been acquired from other C. sakazakii for human BMEC invasion. The presence of a conjugational transfer protein encoded by traX (CSK29544_3p0007) supports this hypothesis.
Cronobacter sakazakii ATCC 29544T harbors three plasmids, pCSK29544_1p (1p), pCSK29544_2p (2p), and pCSK29544_3p (3p). Interestingly, two large plasmids (1p and 3p) encode virulence-associated genes, related to iron uptake and BMEC invasion. Therefore, these two plasmids may contribute to host dominance and pathogenesis, respectively.
Iron is an essential element for dominant survival and colonization via bacterial competition for iron uptake because it plays an important role in the electron transport system for energy production . To accomplish this dominant survival and colonization against other bacteria, C. sakazakii has an iron acquisition system, including siderophore production (iucABCD/iutA) and an ABC-type transport system (eitCBAD) in the plasmid pESA3 . The strain ATCC 29544T also has this privileged iron acquisition system, including siderophore biosynthesis (iucABCD/iutA; CSK29544_1p0024–CSK29544_1p0028) and an ABC-type iron transport system (eitCBAD; CSK29544_1p0056–CSK29544_1p0059). However, this iron acquisition system is present in the plasmid 1p, similar to pESA3 of C. sakazakii BAA-894 . Interestingly, C. sakazakii BAA-894 has multiple copies of pESA3 plasmids in a host, which is highly homologous to the plasmid 1p, suggesting the plasmid 1p may be multi copied, too . These results indicate that the host strain may take advantage of better iron uptake ability in the given environments, probably due to the presence of multiple copies of this iron acquisition system encoded in the plasmid 1p.
Recent studies have revealed that milk formula contains an at least six-times-higher arsenic concentration than that of breast milk, suggesting bacterial survival in milk formula may require arsenic resistance . In addition to the previously mentioned BMEC invasion ability of C. sakazakii ATCC 29544T via the ibeB-homologous cusC gene in the cation efflux system, the plasmid 3p has an arsenic resistance system as encoded by arsRDABC, suggesting the strain ATCC 29544T may have arsenic resistance activity for survival in formula milk powder and that this ability may have been acquired via plasmid conjugational transfer.
Comparative genome analysis
% of Totala
Genome size (bp)
DNA coding (bp)
DNA G + C (bp)
Protein coding genes
Genes with function prediction
This genome information of C. sakazakii ATCC 29544T provides genomic insights into the presence of various virulence factors and their associated pathogenesis mechanisms at molecular level. In addition, this genome information showed that the strain ATCC 29544T has several distinct features probably contributing to dominant survival under the given habitats and even various stress conditions like extremely dry or nutrient-limited conditions. Therefore, this genome information would be useful for further development of a novel strategy for control of C. sakazakii and of rapid detection method using strain-specific DNA markers in genomic level.
brain microvascular endothelial cells
average nucleotide identity
multilocus sequence typing
transmission electron microscopy
clusters of orthologous groups
SK, JHL, and SR initiated and supervised the study. SK and JHL drafted the manuscript. SK conducted laboratory experiments, and SK and YTK worked on the genome sequencing and annotated the genome. SK, YTK, HY, JHL, and SR discussed and analyzed the data and revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The genome sequences of C. sakazakii ATCC 29544T were deposited in the NCBI GenBank server under the accession numbers CP011047, CP011048, CP011049, and CP011050 for a chromosome and three plasmids, respectively.
This research was supported by a grant (14162MFDS972) from the Ministry of Food and Drug Safety, Korea, in 2016 and by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (NRF-2015R1C1A1A01053815).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Iversen C, Lehner A, Mullane N, Bidlas E, Cleenwerck I, Marugg J, et al. The taxonomy of Enterobacter sakazakii: proposal of a new genus Cronobacter gen. nov. and descriptions of Cronobacter sakazakii comb. nov. Cronobacter sakazakii subsp. sakazakii, comb. nov., Cronobacter sakazakii subsp. malonaticus subsp. nov., Cronobacter turicensis sp. nov., Cronobacter muytjensii sp. nov., Cronobacter dublinensis sp. nov. and Cronobacter genomospecies 1. BMC Evol Biol. 2007;7:64.View ArticlePubMedPubMed CentralGoogle Scholar
- Joseph S, Forsythe SJ. Insights into the emergent bacterial pathogen Cronobacter spp., generated by multilocus sequence typing and analysis. Front Microbiol. 2012;3:397.View ArticlePubMed CentralGoogle Scholar
- Farmer JJ, Asbury MA, Hickman FW, Brenner DJ. Enterobacter sakazakii—a new species of Enterobacteriaceae isolated from clinical specimens. Int J Syst Bacteriol. 1980;30:569–84.View ArticleGoogle Scholar
- Hallstrom M, Eerola E, Vuento R, Janas M, Tammela O. Effects of mode of delivery and necrotising enterocolitis on the intestinal microflora in preterm infants. Eur J Clin Microbiol Infect Dis. 2004;23:463–70.View ArticlePubMedGoogle Scholar
- Mao B, Li D, Zhao J, Liu X, Gu Z, Chen YQ, et al. In vitro fermentation of fructooligosaccharides with human gut bacteria. Food Funct. 2015;6:947–54.View ArticlePubMedGoogle Scholar
- Healy B, Cooney S, O’Brien S, Iversen C, Whyte P, Nally J, et al. Cronobacter (Enterobacter sakazakii): an opportunistic foodborne pathogen. Foodborne Pathog Dis. 2010;7:339–50.View ArticlePubMedGoogle Scholar
- Hurrell E, Kucerova E, Loughlin M, Caubilla-Barron J, Forsythe SJ. Biofilm formation on enteral feeding tubes by Cronobacter sakazakii, Salmonella serovars and other Enterobacteriaceae. Int J Food Microbiol. 2009;136:227–31.View ArticlePubMedGoogle Scholar
- Giri CP, Shima K, Tall BD, Curtis S, Sathyamoorthy V, Hanisch B, et al. Cronobacter spp. (previously Enterobacter sakazakii) invade and translocate across both cultured human intestinal epithelial cells and human brain microvascular endothelial cells. Microb Pathog. 2012;52:140–7.View ArticlePubMedGoogle Scholar
- Kilonzo-Nthenge A, Rotich E, Godwin S, Nahashon S, Chen F. Prevalence and antimicrobial resistance of Cronobacter sakazakii isolated from domestic kitchens in middle Tennessee, United States. J Food Prot. 2012;75:1512–7.View ArticlePubMedGoogle Scholar
- Xu X, Li C, Wu Q, Zhang J, Huang J, Yang G. Prevalence, molecular characterization, and antibiotic susceptibility of Cronobacter spp. in Chinese ready-to-eat foods. Int J Food Microbiol. 2015;204:17–23.View ArticlePubMedGoogle Scholar
- Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999;27:4636–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Lukashin AV, Borodovsky M. GeneMark.hmm: new solutions for gene finding. Nucleic Acids Res. 1998;26:1107–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
- Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44:W16–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Choi Y, Kim S, Hwang H, Kim KP, Kang DH, Ryu S. Plasmid-encoded MCP is involved in virulence, motility, and biofilm formation of Cronobacter sakazakii ATCC 29544. Infect Immun. 2015;83:197–204.View ArticlePubMedGoogle Scholar
- Hu L, Grim CJ, Franco AA, Jarvis KG, Sathyamoorthy V, Kothary MH, et al. Analysis of the cellulose synthase operon genes, bcsA, bcsB, and bcsC in Cronobacter species: prevalence among species and their roles in biofilm formation and cell-cell aggregation. Food Microbiol. 2015;52:97–105.View ArticlePubMedGoogle Scholar
- Lehner A, Riedel K, Eberl L, Breeuwer P, Diep B, Stephan R. Biofilm formation, extracellular polysaccharide production, and cell-to-cell signaling in various Enterobacter sakazakii strains: aspects promoting environmental persistence. J Food Prot. 2005;68:2287–94.View ArticlePubMedGoogle Scholar
- Joseph S, Hariri S, Masood N, Forsythe S. Sialic acid utilization by Cronobacter sakazakii. Microb Inform Exp. 2013;3:3.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang B, Brand-Miller J, McVeagh P, Petocz P. Concentration and distribution of sialic acid in human milk and infant formulas. Am J Clin Nutr. 2001;74:510–5.PubMedGoogle Scholar
- Nassif X, Bourdoulous S, Eugene E, Couraud PO. How do extracellular pathogens cross the blood-brain barrier? Trends Microbiol. 2002;10:227–32.View ArticlePubMedGoogle Scholar
- Nair MK, Venkitanarayanan K, Silbart LK, Kim KS. Outer membrane protein A (OmpA) of Cronobacter sakazakii binds fibronectin and contributes to invasion of human brain microvascular endothelial cells. Foodborne Pathog Dis. 2009;6:495–501.View ArticlePubMedGoogle Scholar
- Huang SH, Wan ZS, Chen YH, Jong AY, Kim KS. Further characterization of Escherichia coli brain microvascular endothelial cell invasion gene ibeA by deletion, complementation, and protein expression. J Infect Dis. 2001;183:1071–8.View ArticlePubMedGoogle Scholar
- Huang SH, Chen YH, Fu Q, Stins M, Wang Y, Wass C, et al. Identification and characterization of an Escherichia coli invasion gene locus, ibeB, required for penetration of brain microvascular endothelial cells. Infect Immun. 1999;67:2103–9.PubMedPubMed CentralGoogle Scholar
- Wang Y, Huang SH, Wass CA, Stins MF, Kim KS. The gene locus yijP contributes to Escherichia coli K1 invasion of brain microvascular endothelial cells. Infect Immun. 1999;67:4751–6.PubMedPubMed CentralGoogle Scholar
- Mietzner TA, Morse SA. The role of iron-binding proteins in the survival of pathogenic bacteria. Annu Rev Nutr. 1994;14:471–93.View ArticlePubMedGoogle Scholar
- Franco AA, Hu L, Grim CJ, Gopinath G, Sathyamoorthy V, Jarvis KG, et al. Characterization of putative virulence genes on the related RepFIB plasmids harbored by Cronobacter spp. Appl Environ Microbiol. 2011;77:3255–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Kucerova E, Clifton SW, Xia XQ, Long F, Porwollik S, Fulton L, et al. Genome sequence of Cronobacter sakazakii BAA-894 and comparative genomic hybridization analysis with other Cronobacter species. PLoS ONE. 2010;5:e9556.View ArticlePubMedPubMed CentralGoogle Scholar
- Carignan CC, Cottingham KL, Jackson BP, Farzan SF, Gandolfi AJ, Punshon T, et al. Estimated exposure to arsenic in breastfed and formula-fed infants in a United States cohort. Environ Health Perspect. 2015;123:500–6.PubMedPubMed CentralGoogle Scholar
- Harvie DR, Ellar DJ. A ferric dicitrate uptake system is required for the full virulence of Bacillus cereus. Curr Microbiol. 2005;50:246–50.View ArticlePubMedGoogle Scholar
- Ratnaike RN. Acute and chronic arsenic toxicity. Postgrad Med J. 2003;79:391–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Petkau A, Stuart-Edwards M, Stothard P, Van Domselaar G. Interactive microbial genome visualization with GView. Bioinformatics. 2010;26:3125–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31:3691–3.View ArticlePubMedPubMed CentralGoogle Scholar