- Genome report
- Open Access
Complete genome sequence of Vibrio vulnificus FORC_017 isolated from a patient with a hemorrhagic rash after consuming raw dotted gizzard shad
Gut Pathogensvolume 8, Article number: 22 (2016)
Vibrio vulnificus, a resident in the human gut, is frequently found in seafood, causing food-borne illnesses including gastroenteritis and severe septicemia. While V. vulnificus has been known to be one of the major food-borne pathogens, pathogenicity and virulence factors are not fully understood yet. To extend our understanding of the pathogenesis of V. vulnificus at the genomic level, the genome of V. vulnificus FORC_017 isolated from a female patient experiencing a hemorrhagic rash was completely sequenced and analyzed.
Three discontinuous contigs were generated from a hybrid assembly using Illumina MiSeq and PacBio platforms, revealing that the genome of the FORC_017 consists of two circular chromosomes and a plasmid. Chromosome I consists of 3,253,417-bp (GC content 46.49 %) containing 2943 predicted open reading frames (ORFs) and chromosome II of 1,905,745-bp (GC content 46.90 %) containing 1638 ORFs. The plasmid pFORC17 consists of 70,069-bp (GC content 43.77 %) containing 84 ORFs. The average nucleotide identity (ANI) value of the FORC_017 and CMCP6 strains was 98.53, suggesting that they are closely related.
Pathogenesis-associated genes including vvhA, rtx gene cluster, and various hemolysin genes were present in FORC_017. In addition, three complete secretion systems (Type I, II and VI) as well as iron uptake-related genes for virulence of the FORC_017 were detected, suggesting that this strain is pathogenic. Further comparative genome analysis revealed that FORC_017 and CMCP6 share major toxin genes including vvhA and rtx for pathogenesis activities. The genome information of the FORC_017 provides novel insights into pathogenicity and virulence factors of V. vulnificus.
Vibrio vulnificus infections in humans via either the consumption of contaminated seafood or contact of a wound with seawater generally cause mild gastroenteritis with diarrhea, vomiting, and abdominal pain . In addition, this pathogen in the human gut can cause life-threatening systemic septicemia with a mortality rate of 61 % when it crosses the intestinal mucosal barrier to enter the bloodstream, suggesting that V. vulnificus is one of many pathogenic bacteria in the human gut microbiota [2, 3]. Vibrio vulnificus has been frequently detected in sea water as well as various kinds of seafood . Vibrio vulnificus was detected in 66 of 68 oyster samples (97 %) and 27 of 30 clam samples (90 %) from Long Island Sound, California . The consumption of V. vulnificus-contaminated seafood causes clinical food-borne outbreaks. In South Korea (2000–2011), 34 individuals who consumed V. vulnificus-contaminated seafood were hospitalized and 16 died, which suggests that this pathogen is widespread and highly virulent and causes food-borne outbreaks . Therefore, virulence traits of this pathogen should be understood at the genomic level. Although outbreaks of this pathogen have been reported for many years, only five strains of V. vulnificus have had their whole genomes sequenced and analyzed to date (NCBI; http://www.ncbi.nlm.nih.gov). Genome information has revealed the presence of various virulence factors in the genomes, such as an extracellular hemolysin/cytolysin, iron uptake receptors, capsular polysaccharide related to antiphagocytosis, and secretion systems . To extend our understanding of the pathogen’s further pathogenicity and virulence factors at the genomic level, more genomes of V. vulnificus need to be sequenced and analyzed.
To elucidate the pathogenesis of V. vulnificus in South Korea, a clinical strain, V. vulnificus FORC_017, was isolated from a 63-year-old woman with the symptom of a hemorrhagic rash after consuming a raw dotted gizzard shad on Jeju Island. The genome was completely sequenced and analyzed using various bioinformatics tools. In addition, comparative genome analysis of the FORC_017 strain and other clinical isolates was performed. The results of this study would provide genomic insights into the main factors associated with food-borne outbreaks and pathogenesis, and is probably useful for further control of this highly virulent pathogen in seafood.
Growth condition and genomic DNA extraction
The FORC_017 strain was cultivated aerobically at 30 °C with modified Luria-Bertani medium (LB) supplemented with 1 % (w/v) NaCl for 12 h. The genomic DNA was extracted using DNeasy Blood & Tissue Kit (Qiagen, Valencia, CA, USA), according to the manufacturers protocol.
Genome sequencing and annotation
The whole genome sequencing and assembly were conducted at ChunLab, Inc. (Seoul, South Korea). Sequences of V. vulnificus FORC_017 were obtained using two different platforms of PacBio RS II platform (Pacific Biosciences, Menlo Park, CA, USA) and Illumina MiSeq platform (Illumina, San Diego, CA, USA). The sequencing results were summarized in Additional file 1: Table S1. The qualified sequence reads were assembled with CLC Genomic Workbench 7.5.1 (CLC Bio, Aarhus, Denmark). The prediction of open reading frames (ORFs) and their annotations were performed using the GeneMarkS program and the Rapid Annotations using Subsystems Technology (RAST) server [8, 9]. The prediction of ribosome binding sites (RBS) was conducted using the RBSfinder (J. Craig Venter Institute, Rockville, MD, USA). Subsequent predictions of the functions of ORFs and their conserved protein domains were carried out using InterProScan 5 and COG-based WebMGA programs [10, 11]. The putative virulence factors of the V. vulnificus FORC_017 strain were predicted and characterized using the basic local alignment search tool (BLAST) in the Virulence Factor Database (VFDB; http://www.mgc.ac.cn/VFs/main.htm). The circular genome maps of the FORC_017 strain were drawn using the GenVision program (DNASTAR, Madision, WI, USA).
Comparative genome analysis
The strain was identified as V. vulnificus using 16S rRNA sequencing and its identity was confirmed with comparative phylogenetic tree analysis using MEGA6  with various 16S rRNA sequences of the genus Vibrio. Phylogenetic trees were constructed from the aligned sequences using the neighbor-joining method with 1000 bootstrap replicates via MEGA6 software. Shewanella baltica OS678 was used as the out-group. ANI phylogenetic tree analysis of complete genome sequences of V. vulnificus (93U204, YJ016, CMCP6, FORC_016 and MO6-24/O) was conducted using the JSpecies program  to show the evolutionary relationship among V. vulnificus strains. The ANI values were calculated using JSpecies by comparing whole genome sequences between V. vulnificus strains, which were fragmented into 1020-bp, based on BLAST algorism. The tree was constructed using the R program. The comparative genome analysis was conducted with the Artemis Comparison Tool (ACT) .
A single colony of the FORC_017 strain was isolated and its morphology was observed using a transmission electron microscopy showing a traditional shape of V. vulnificus, including a short rod with a single flagellum (Additional file 1: Figure S1). The genomic DNA of the FORC_017 strain was extracted from the isolate and confirmed to V. vulnificus by 16S rRNA gene sequencing. The 16S rRNA genes were amplified using PCRwith a 16S rRNA universal primer set and the PCR products were sequenced using capillary sequencing technology for taxonomic identification . In addition, when the quality scores of the raw read sequences from NGS sequencers were more than 40 as cutoffs, they were selected and assembled. The assembled genome sequence was also verified with comparative ANI analysis with the published complete genome sequences of the same species.
The hybrid platforms of MiSeq and PacBio with a total fold coverage of 438.12 presented that the complete genome of V. vulnificus FORC_017 consists of two circular DNA chromosomes and a plasmid. Chromosome I consists of 3,253,417-bp with a GC content of 46.49 % containing 2943 predicted ORFs, 102 tRNA genes, and 31 rRNA genes. Among the ORFs, 2326 ORFs (79.03 %) were predicted to be functional and 617 ORFs (20.96 %) to encode hypothetical proteins. Chromosome II consists of 1,905,745-bp with a GC content of 46.90 % containing 1638 predicted ORFs, 13 tRNA genes, and 3 rRNA genes. Among the ORFs, 1256 ORFs (76.68 %) were predicted to be functional and another 382 ORFs (23.32 %) to encode hypothetical proteins. The plasmid, pFORC17, consists of 70,069-bp with a GC content of 43.77 % containing 84 predicted ORFs. 36 ORFs (42.86 %) were predicted to be functional and another 48 ORFs (57.14 %) to encode hypothetical proteins. Based on the bioinformatics analysis of these chromosomes and a plasmid, circular genome maps were drawn (Fig. 1).
Pathogenesis and virulence factors
Key virulence factors in V. vulnificus are a hemolysin/cytolysin (VvhA) exotoxin and a repeats-in-toxin (Rtx) (Additional file 1: Table S2) . The genome of the FORC_017 strain has a vvhA gene (FORC17_3853) and an rtx gene cluster (FORC17_3923, 3924, 3926 and 3927), suggesting that FORC_017 is pathogenic. The FORC_017 strain has various hemolysin genes such as FORC17_0475, 1369, 2364, 2546, 2587 and 2903 in chromosome I, and FORC17_3057, 3239 and 4284 in chromosome II. Subsequent cytotoxicity test of the FORC_017 strain supported its pathogenesis activity (Additional file 1: Figure S2).
Various secretion systems including Types I, II, IV and VI (T1SS, T2SS, T4SS and T6SS) were observed in the genome sequences of the FORC_017 strain. T1SS secretes gram-negative bacterial RtxA exotoxins in Vibrio, which play a role in pore formation in mammalian cells to cause death . T1SS was detected in the specific region of chromosome II (FORC17_4414 to FORC17_4539). The presence of rtx gene cluster and T1SS suggests the pathogenesis of the FORC_017 strain toward the human cells. T2SS secretes bacterial cholera-like toxins as well as proteases and lipases. This T2SS was found in the specific region of chromosome I (FORC17_2808 to FORC17_2819). However, the FORC_017 strain does not have this cholera-like toxin in the genome. Therefore, T2SS may be associated with protease/lipase secretions, which is probably related to the damage of the membrane in the human cells for infection of V. vulnificus . While T4SS is generally known to be associated with the translocation of bacterial virulence factors into human cells , the FORC_017 strain has only one gene encoding VirD2 (FORC17_3580), suggesting that the FORC_017 strain does not have a completely functional T4SS for translocation of virulence factors. T6SS is known to be associated with bacterial toxin protein injection machinery for virulence to human cells. This T6SS was encoded in a gene cluster of chromosome II (FORC17_3873 to FORC17_3890) . Therefore, these secretion systems may contribute to the virulence toward human cells for food-borne illnesses.
Iron uptake from host cells plays a key role in the survival of V. vulnificus since iron loss caused by V. vulnificus negatively affects the heme protein integrity, resulting in host cell death [22, 23]. The FORC_017 genome included heme receptors (FORC17_3295 in chromosome II), iron ABC transporters (FORC17_3542 to FORC17_3544 in chromosome II), and vibriobactin receptors (FORC17_4247 to FORC17_4256 in chromosome II). These iron uptake-related genes support the survival of the FORC_017 strain in the host cells.
Comparative genome analysis of the closely related FORC_017 and CMCP6 strains
Comparative phylogenetic tree analysis of the FORC_017 strain and other Vibrio strains using 16S rRNA sequences showed that the FORC_017 strain belongs to V. vulnificus (Fig. 2a). Further ANI phylogenetic tree analysis of complete genome sequences of V. vulnificus (93U204, YJ016, CMCP6, FORC_016 and MO6-24/O) revealed that the FORC_017 strain has the closest evolutionary relationship (ANI value 98.53) with a clinical isolate, V. vulnificus CMCP6, from a patient in South Korea (Fig. 2b) . Interestingly, comparative genome analysis of two closely related strains, FORC_017 and CMCP6, showed that the FORC_017 strain has additional virulence factors including thermostable hemolysin delta-VPH (FORC17_3057) ansd partial T4SS (FORC17_3580). Although this thermostable hemolysin may be involved in the pathogenesis toward humans, the presence of a single gene related to T4SS in the FORC_017 strain is not fully understood yet. Probably, this single gene encoding VirD2 may be a remnant of a complete T4SS gene cluster after its deletion event from chromosome II of the FORC_017 strain. The presence of mobile element proteins near VirD2 (FORC17_3591 and 3592) supports this hypothesis.
The genome information of the FORC_017 strain provides genomic insights into its pathogenicity and various virulence factors as well as their virulence mechanisms in V. vulnificus, which is probably involved in the seafood-borne outbreaks. In addition, the FORC_017 strain could provide further genomic information for development of rapid detection and control methods against V. vulnificus in various kinds of seafood. Therefore, genome information can be applied for rapid detection, epidemiological survey and prevention of food poisoning caused by this strain.
open reading frame
average nucleotide identity
transmission electron microscopy
Food-borne Pathogen Omics Research Center
Rapid Annotation using Subsystem Technology
ribosome binding site
Artemis Comparison Tool
type I secretion system
type II secretion system
type IV secretion system
type VI secretion system
Mark SS, Rohinee NP. Epidemiology and pathogenesis of Vibrio vulnificus. Microbes Infect. 2000;2:177–88.
Lee BC, Choi SH. Vibrio vulnificus RTX toxin plays an important role in the apoptotic death of human intestinal epithelial cells exposed to Vibrio vulnificus. Microbes Infect. 2008;10:1504–13.
Shapiro RL, Altekruse S, Hutwagner L, Bishop R, Hammond R, Wilson S, Ray B, Thompson S, Tauxe RV, Griffin PM. The role of gulf coast oysters harvested in warmer months in Vibrio vulnificus infections in the United States, 1988–1996. J Infect Dis. 1998;178:752–9.
Nadja B, Claudia J, Ralf D, Nicole B, Simone IB, Eckhard S. Virulence profiles of Vibrio vulnificus in German coastal waters, a comparison of North sea and Baltic sea isolates. Int J Environ Res Public Health. 2015;12:15943–59.
Jessica LJ, Catharina HM, John CB, Kristin DB, David HC, William H. Abundance of Vibrio cholerae, V. vulnificus, and V. parahaemolyticus in Oysters (Crassostrea virginica) and Clams (Mercenaria mercenaria) from Long Island Sound. Appl Environ Microb. 2014;80:7667–72.
Yun NR, Kim DM, Lee J, Han MA. pH level as a marker for predicting death among patients with Vibrio vulnificus infection, South Korea, 2000–2011. Emerg Infect Dis. 2015;21:259.
Linkous DA, Oliver JD. Pathogenesis of Vibrio vulnificus. FEMS Microbiol Lett. 1999;174:207–14.
Besemer J, Lomsadze A, Borodovsky M. GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001;29:2607–18.
Aziz RK, Bartels D, Best AA, DeJongh M, et al. The RAST server: rapid annotations using subsystems technology. BMC Genom. 2008;9:75.
Jones P, Binns D, Chang HY, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.
Wu S, Zhu Z, Fu L, et al. WebMGA: a customizable web server for fast metagenomic sequence analysis. BMC Genom. 2011;12:444.
Tamura K, Stecher G, Peterson D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Richter M, Rosselló MR. Shifting the genomic gold standard for the prokaryotic species definition. Proc Nat Acad Sci USA. 2009;106:19126–31.
Carver TJ, Rutherford KM, Berriman M, et al. ACT: the artemis comparison tool. Bioinformatics. 2005;21:3422–3.
Ku HJ, Lee JH. 16S rRNA universal primer, fecal DNA extraction, gastrointestinal microbiota, metagenomics, next-generation sequencing. J Microbiol Biotechnol. 2014;24:812–22.
Kim YR, Lee SE, Kook H, Yeom JA, Na HS, et al. Vibrio vulnificus RTX toxin kills host cells only after contact of the bacteria with host cells. Cell Microbiol. 2008;10:848–62.
Bethany KB, Brian MM, Karla JFS. Growth phase regulation of Vibrio cholerae RTX toxin export. J Bacteriol. 2007;189:1827–35.
Marissa EY, Konstantin VK, Jan A, Wim GJ. Structure of the minor pseudopilin EpsH from the type 2 secretion system of Vibrio cholerae. J Mol Biol. 2008;377:91–103.
Christie PJ, Krishnamohan A, Krishnamoorthy V, Jakubowski S, Cascales E. Biogenesis, architecture, and function of bacterial type IV secretion systems. Annu Rev Microbiol. 2005;59:451–85.
Coulthurst SJ. The type VI secretion system—a widespread and versatile cell targeting system. Res Microbiol. 2013;164:640–54.
Kim YR, Lee SE, Kim CM, Kim SY, Shin EK, et al. Characterization and pathogenic significance of Vibrio vulnificus antigens preferentially expressed in septicemic patients. Infect Immun. 2003;71:5461–71.
Nishina Y, Miyoshi S, Nagase A, Shinoda S. Significant role of an exocellular protease in utilization of heme by Vibrio vulnificus. Infect Immun. 1992;60:2128–32.
Helms SD, Oliver JD, Travis JC. Role of heme compounds and haptoglobin in Vibrio vulnificus pathogenicity. Infect Immun. 1984;45:345–9.
SR, HK, JHL and SHC designed the study. KHL and STH provided the strain. HYC, YTK, HJK and JHL analyzed the sequencing data. HYC, SK, EJN performed experiments. HYC, JHL and SHC discussed the results and wrote the manuscript. All authors read and approved the final manuscript.
This research was supported by a Grant (14162MFDS972) from the Ministry of Food and Drug Safety, Korea in 2015.
Availability of data and materials
The genome sequences of Vibrio vulnificus FORC_017 were deposited in the NCBI GenBank server under the accession numbers CP012739, CP012740, and CP012741 for chromosome I, II, and a plasmid, respectively.
The authors declare that they have no competing interests.
Ethics approval and consent to participate