Complete genome sequence of the Vibrio vulnificus strain VV2014DJH, a human-pathogenic bacterium isolated from a death case in China
© The Author(s) 2017
Received: 24 July 2017
Accepted: 8 November 2017
Published: 21 November 2017
Vibrio vulnificus, an opportunistic pathogen, is the causative agent of life-threatening septicemia and severe wound infections. However, the pathogenicity and virulence factors of V. vulnificus are not fully understood. Here we report the complete genome sequence of V. vulnificus VV2014DJH, which was isolated from a death case.
The genome of the V. vulnificus VV2014DJH contains two circular chromosomes with a mean G+C content of 46.8%, but does not consists of any plasmids. The chromosome I and chromosome II consist of 3,303,590 and 1,770,972 bp, respectively. In addition, the genome consists of 4617 protein coding genes, 172 RNA genes and type I, II and III secretion systems were predicted.
In this study, the genomic information of the V. vulnificus VV2014DJH has been described. The information would contribute to the increasing scope and depth of Vibrio genome database, and provide insights into the pathogenicity and virulence factors of V. vulnificus.
Vibrio vulnificus is an opportunistic Gram-negative pathogen, which is widely distributed in marine environments around the world and has been isolated form water, sediments and seafood . V. vulnificus is the leading cause of deaths reported from seafood in the United States with approximately 40 cases reported per year, and the fatality rate exceeds 50% . However, the mortality rate ranges from 18 to 56% in China [3, 4]. In the recent years, more cases of amputation and death caused by V. vulnificus infections were reported in China, and most of these cases occurred in the coastal area of Zhejiang province . In this paper, we report the complete genome sequence of V. vulnificus strain VV2014DJH, which was isolated from the blood culture of the death case of a shellfish aquaculture worker in Taizhou, Zhejiang province . Determination of the genomic information of V. vulnificus strain VV2014DJH is necessary to understand the pathogenesis of the strain. Our study was approved by the ethics committee of the Zhejiang Provincial Center for Disease Control and Prevention (ZJCDC), China.
The genomic DNA of VV2014DJH was extracted using the DNeasy Blood and Tissue Kit (QIAGEN, Germany) for Gram-negative bacteria, according to the manufacturer’s instructions. Total DNA was subjected to quality control by 2% agarose gel electrophoresis and quantified by a NanoDrop™ spectrophotometer. The sequencing of VV2014DJH was performed on Pacbio RS with Single Molecule, Real-Time (SMRT) technology. SMRT Analysis 2.3.0 was used to filter low quality reads, whose quality scores were under 40, and the filtered reads were assembled by SOAPdenovo2 to generate one contig without gaps . According to Clusters of Orthologous Group (COG) category, circular representations with the annotated genes were obtained using Circos software and displayed the diagram of the VV2014DJH genome.
Transfer RNA (tRNA) genes were predicted by tRNAscan-SE . Ribosome RNA (rRNA) genes were predicted with rRNAmmer  and sRNAs were obtained by BLAST against Rfam database . RepeatMasker  (http://www.repeatmasker.org/) was used to detect repetitive sequences, and Tandem Repeat Finder  (http://tandem.bu.edu/trf/trf.html) was used to find Tandem Repeats. Gene prediction was performed on the VV2014DJH genome assembly by GeneMarkS  (http://topaz.gatech.edu/) with integrated model which combine the GeneMarkS generated (native) and Heuristic model parameters.
A whole genome BLAST  search (E-value less than 1E−5, minimal alignment length percentage larger than 40%) was performed against five databases. They are KEGG [15–17] (Kyoto Encyclopedia of Genes and Genomes), COG [18, 19], NR (Non-Redundant Protein Database), Swiss-Prot , and GO  (Gene Ontology) . Secretory proteins were detected in the genome assembly by SignalP . Type I–VII secretion system related proteins were extracted from all the annotation results. Type III secretion system effector proteins were detected by EffectiveT3 also . The MUMmer system (version 3.0) was used applied for colinearity analysis of to compare the genome sequences of VV2014DJH, with CMCP6 and YJ016 . The virulence factors annotation in VV2014DJH was performed by blasting the genome sequence to the Virulence factors Pathogenic Bacteria database (VFDB, http://www.mgc.ac.cn/cgi-bin/VFs/jsif/main.cgi) .
The genomic DNA used for sequencing was isolated from a single colony of the VV2014DJH. The 16S rRNA gene was sequenced and BLAST was performed against the NCBI database. In addition, the raw read sequences were selected and assembled only when their quality scores were more than 40 as cutoffs.
A PacBio RS sequencing run resulted in 56,579 reads with a mean read length of 8107 bp and a N50 read length of 11,532 bp. PacBio RS reads were assembled into one polished contigs with a N50 contig length of 3,307,273. The genome of VV2014DJH has a size of 5,074,562 bp with a mean G+C content of 46.8% and is composed of two circular chromosomes (chromosome I: 3,303,590 bp, 46.6% G+C content; chromosome II: 1,770,972 bp, 47.3% G+C content). It does not contain any plasmids.
General features of strain VV2014DJH genome
16 S rRNA
23 S rRNA
Protein coding genes
Interspersed repetitive sequences
Total length (bp)
In genome (%)
Average length (bp)
Tandem repetitive sequences
Repeat size (bp)
Total length (bp)
In genome (%)
Three hundred and fifty-seven secretory proteins were detected in the 4617 proteins, and 45 proteins related to type I, II, III, IV and VI secretion systems were predicted, no protein related to type V and VII secretion systems was found. Type III secretion systems are complex bacterial structures that provide Gram-negative pathogens with unique virulence mechanisms. In the strain VV2014DJH, 184 effector proteins were detected by EffectiveT3 software.
The results of colinearity analysis for these three genome sequences of VV2014DJH, CMCP6 and YJ016 showed that the majority of the genomes were conserved in all three strains, although these were some transpositions between VV2014DJH and CMCP6, or VV2014DJH and YJ016. The Blast results showed that the pathogenic factors of VV2014DJH including the aforementioned type II and type VI secretion systems, flagellin C/D/E, and other pathogenic factors were also found in V. vulnificus VV2014DJH (Additional file 1: Fig S1, Table S1).
Our report provides an extended understanding on Vibrio vulnificus at genomic level and may provide useful information for epidemiological investigation, pathogenicity and virulence factors. In addition, the genome information is probably helpful for rapid detection and prevention of food poisoning by the strain.
JHP was responsible for the most of the experiments; YS, WWY and HYM collected and analyzed the sequencing data; YJZ and MYZ reviewed the article. All authors read and approved the final manuscript.
This work was supported by the Major Science and Technology Special Projects of Zhejiang, China (No. 2013C03045-1) and the Medical Scientific Research Foundation of Zhejiang Province, China (No. 2012RCA013).
The authors declare that they have no competing interests.
Availability of data and materials
This whole genome sequences of V. vulnificus VV2014DJH has been deposited in the DDBJ/EMBL/GenBank database under the accession numbers CP019320 (ChI) and CP019321 (ChII).
Consent for publication
Ethics approval and consent to participate
Our study was approved by the ethics committee of the Zhejiang Provincial Center for Disease Control and Prevention (ZJCDC), China.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Gulig PA, Bourdage KL, Starks AM. Molecular pathogenesis of Vibrio vulnificus. J Microbiol. 2005;43 Spec No:118–31.PubMedGoogle Scholar
- Jones MK, Oliver JD. Vibrio vulnificus: disease and pathogenesis. Infect Immun. 2009;77:1723–33. https://doi.org/10.1128/iai.01046-08.View ArticlePubMedPubMed CentralGoogle Scholar
- Hong G, et al. Emergency treatment of 16 patients with necrotizing fasciitis caused by Vibrio vulnificus infection complicated with septic shock. Chin Med J. 2014;127:1984–6.PubMedGoogle Scholar
- Chao WN, et al. Impact of timing of surgery on outcome of Vibrio vulnificus-related necrotizing fasciitis. Am J Surg. 2013;206:32–9. https://doi.org/10.1016/j.amjsurg.2012.08.008.View ArticlePubMedGoogle Scholar
- Zhao H, et al. Correlations between clinical features and mortality in patients with Vibrio vulnificus infection. PLoS ONE. 2015;10:e0136019. https://doi.org/10.1371/journal.pone.0136019.View ArticlePubMedPubMed CentralGoogle Scholar
- Pan JH, et al. Drug resistance phenotype and molecular characteristics of a virulent Vibrio vulnificus strain isolated fom a fatal case. Dis Surveill. 2016;31:229–32.Google Scholar
- Luo R, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 2012;1:18. https://doi.org/10.1186/2047-217x-1-18.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
- Lagesen K, et al. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8. https://doi.org/10.1093/nar/gkm160.View ArticlePubMedPubMed CentralGoogle Scholar
- Gardner PP, et al. Rfam: updates to the RNA families database. Nucleic Acids Res. 2009;37:D136–40. https://doi.org/10.1093/nar/gkn766.View ArticlePubMedGoogle Scholar
- Saha S, Bridges S, Magbanua ZV, Peterson DG. Empirical comparison of ab initio repeat finding programs. Nucleic Acids Res. 2008;36:2284–94. https://doi.org/10.1093/nar/gkn064.View ArticlePubMedPubMed CentralGoogle Scholar
- Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27:573–80.View ArticlePubMedPubMed CentralGoogle Scholar
- 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
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10. https://doi.org/10.1016/S0022-2836(05)80360-2.View ArticlePubMedGoogle Scholar
- Kanehisa M. A database for post-genome analysis. Trends Genet. 1997;13:375–6.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:D277–80. https://doi.org/10.1093/nar/gkh063.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanehisa M, et al. From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 2006;34:D354–7. https://doi.org/10.1093/nar/gkj102.View ArticlePubMedGoogle Scholar
- Tatusov RL, Koonin EV, Lipman DJ. A genomic perspective on protein families. Science. 1997;278:631–7.View ArticlePubMedGoogle Scholar
- Tatusov RL, et al. The COG database: an updated version includes eukaryotes. BMC Bioinform. 2003;4:41. https://doi.org/10.1186/1471-2105-4-41.View ArticleGoogle Scholar
- Magrane M. UniProt Knowledgebase: a hub of integrated protein data. Database. 2011;2011:bar009. https://doi.org/10.1093/database/bar009.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashburner M, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9. https://doi.org/10.1038/75556.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu K, Zhang T. Construction of customized sub-databases from NCBI-nr database for rapid annotation of huge metagenomic datasets using a combined BLAST and MEGAN approach. PLoS ONE. 2013;8:e59831. https://doi.org/10.1371/journal.pone.0059831.View ArticlePubMedPubMed CentralGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6. https://doi.org/10.1038/nmeth.1701.View ArticlePubMedGoogle Scholar
- Jehl MA, Arnold R, Rattei T. Effective—a database of predicted secreted bacterial proteins. Nucleic Acids Res. 2011;39:D591–5. https://doi.org/10.1093/nar/gkq1154.View ArticlePubMedGoogle Scholar
- Kurtz S, et al. Versatile and open software for comparing large genomes. Genome Biol. 2004;5:R12. https://doi.org/10.1186/gb-2004-5-2-r12.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen L, Zheng D, Liu B, Yang J, Jin Q. VFDB 2016: hierarchical and refined dataset for big data analysis—10 years on. Nucleic Acids Res. 2016;44:D694–7. https://doi.org/10.1093/nar/gkv1239.View ArticlePubMedGoogle Scholar
- Krzywinski M, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45. https://doi.org/10.1101/gr.092759.109.View ArticlePubMedPubMed CentralGoogle Scholar