Open Access

Complete genome sequence of Vibrio vulnificus FORC_017 isolated from a patient with a hemorrhagic rash after consuming raw dotted gizzard shad

  • Han Young Chung1,
  • You-Tae Kim2,
  • Suyeon Kim1,
  • Eun Jung Na1,
  • Hye-Jin Ku2,
  • Keun Hwa Lee3,
  • Sang Taek Heo3,
  • Sangryeol Ryu1,
  • Heebal Kim4,
  • Sang Ho Choi1 and
  • Ju-Hoon Lee2Email author
Gut Pathogens20168:22

https://doi.org/10.1186/s13099-016-0104-6

Received: 10 April 2016

Accepted: 10 May 2016

Published: 20 June 2016

Abstract

Background

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.

Results

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.

Conclusions

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.

Keywords

Vibrio vulnificus Genome sequencingHybrid assemblyPathogenesisVirulence factorHemolysinSecretion system

Background

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 [1]. 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 [4]. Vibrio vulnificus was detected in 66 of 68 oyster samples (97 %) and 27 of 30 clam samples (90 %) from Long Island Sound, California [5]. 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 [6]. 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 [7]. 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.

Methods

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 [12] 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 [13] 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) [14].

Quality assurance

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 [15]. 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.

Initial findings

Genome properties

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).
Fig. 1

The circular genome maps of the strain FORC_017 using the GenVision program (DNASTAR, Madison, WI, USA). a Chromosome I. b Chromosome II. c Plasmid FORC_017. The outer circle indicates the locations of all annotated ORFs and the inner circle with the red peaks indicates GC content. Between these circles, sky blue arrows indicate the rRNA operons and orange arrows indicate the tRNAs. All annotated ORFs were colored differently according to the COG assignments. Genes with specialized functions labeled with different colors as follows; virulence-related genes in red and other functional genes in black

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) [16]. 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 [17]. 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 [18]. While T4SS is generally known to be associated with the translocation of bacterial virulence factors into human cells [19], 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) [20]. 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) [21]. 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.
Fig. 2

a The phylogenetic tree of the FORC_017 strain with closely related Vibrio species aligned by using ClustalW. Shewanella baltica OS678 was used as the out-group. b The ANI genome tree of various V. vulnificus strains were constructed using the ANI values. The 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

Future directions

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.

Abbreviations

ORF: 

open reading frame

ANI: 

average nucleotide identity

TEM: 

transmission electron microscopy

FORC: 

Food-borne Pathogen Omics Research Center

VvhA: 

hemolysin/cytolysin exotoxin

Rtx: 

a repeats-toxin

LB: 

Luria-Bertani medium

RAST: 

Rapid Annotation using Subsystem Technology

RBS: 

ribosome binding site

ACT: 

Artemis Comparison Tool

T1SS: 

type I secretion system

T2SS: 

type II secretion system

T4SS: 

type IV secretion system

T6SS: 

type VI secretion system

Declarations

Authors’ contributions

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.

Acknowledgements

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.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

Not applicable.

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.

Authors’ Affiliations

(1)
Department of Agricultural Biotechnology, Center for Food Safety and Toxicology, Seoul National University
(2)
Department of Food Science and Biotechnology, Kyung Hee University
(3)
Department of Internal Medicine, Jeju National University School of Medicine
(4)
Department of Animal Science and Biotechnology, Seoul National University

References

  1. Mark SS, Rohinee NP. Epidemiology and pathogenesis of Vibrio vulnificus. Microbes Infect. 2000;2:177–88.View ArticleGoogle Scholar
  2. 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.View ArticlePubMedGoogle Scholar
  3. 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.View ArticlePubMedGoogle Scholar
  4. 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.View ArticleGoogle Scholar
  5. 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.View ArticleGoogle Scholar
  6. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Linkous DA, Oliver JD. Pathogenesis of Vibrio vulnificus. FEMS Microbiol Lett. 1999;174:207–14.View ArticlePubMedGoogle Scholar
  8. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Aziz RK, Bartels D, Best AA, DeJongh M, et al. The RAST server: rapid annotations using subsystems technology. BMC Genom. 2008;9:75.View ArticleGoogle Scholar
  10. Jones P, Binns D, Chang HY, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.View ArticlePubMedPubMed CentralGoogle Scholar
  11. Wu S, Zhu Z, Fu L, et al. WebMGA: a customizable web server for fast metagenomic sequence analysis. BMC Genom. 2011;12:444.View ArticleGoogle Scholar
  12. Tamura K, Stecher G, Peterson D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Richter M, Rosselló MR. Shifting the genomic gold standard for the prokaryotic species definition. Proc Nat Acad Sci USA. 2009;106:19126–31.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Carver TJ, Rutherford KM, Berriman M, et al. ACT: the artemis comparison tool. Bioinformatics. 2005;21:3422–3.View ArticlePubMedGoogle Scholar
  15. Ku HJ, Lee JH. 16S rRNA universal primer, fecal DNA extraction, gastrointestinal microbiota, metagenomics, next-generation sequencing. J Microbiol Biotechnol. 2014;24:812–22.View ArticlePubMedGoogle Scholar
  16. 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.View ArticlePubMedGoogle Scholar
  17. Bethany KB, Brian MM, Karla JFS. Growth phase regulation of Vibrio cholerae RTX toxin export. J Bacteriol. 2007;189:1827–35.View ArticleGoogle Scholar
  18. 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.View ArticleGoogle Scholar
  19. 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.View ArticlePubMedGoogle Scholar
  20. Coulthurst SJ. The type VI secretion system—a widespread and versatile cell targeting system. Res Microbiol. 2013;164:640–54.View ArticlePubMedGoogle Scholar
  21. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  22. 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.PubMedPubMed CentralGoogle Scholar
  23. Helms SD, Oliver JD, Travis JC. Role of heme compounds and haptoglobin in Vibrio vulnificus pathogenicity. Infect Immun. 1984;45:345–9.PubMedPubMed CentralGoogle Scholar

Copyright

© The Author(s) 2016