Skip to main content

Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus



Vibrio parahaemolyticus is a Gram-negative halophilic bacterium. Infections with the bacterium could become systemic and can be life-threatening to immunocompromised individuals. Genome sequences of a few clinical isolates of V. parahaemolyticus are currently available, but the genome dynamics across the species and virulence potential of environmental strains on a genome-scale have not been described before.


Here we present genome sequences of four V. parahaemolyticus clinical strains from stool samples of patients and five environmental strains in Hong Kong. Phylogenomics analysis based on single nucleotide polymorphisms revealed a clear distinction between the clinical and environmental isolates. A new gene cluster belonging to the biofilm associated proteins of V. parahaemolyticus was found in clincial strains. In addition, a novel small genomic island frequently found among clinical isolates was reported. A few environmental strains were found harboring virulence genes and prophage elements, indicating their virulence potential. A unique biphenyl degradation pathway was also reported. A database for V. parahaemolyticus ( was constructed here as a platform to access and analyze genome sequences and annotations of the bacterium.


We have performed a comparative genomics analysis of clinical and environmental strains of V. parahaemolyticus. Our analyses could facilitate understanding of the phylogenetic diversity and niche adaptation of this bacterium.


Foodborne diseases remain a serious public health problem worldwide. In Hong Kong, infections due to foodborne pathogens are also a common and important public health issue. Among all causative agents resulting in foodborne outbreaks, bacteria have caused more than 80% of the cases. In 2006, the Centre for Health Protection of Hong Kong revealed more than 800 local foodborne outbreaks due to bacterial causative agents, inflicting more than 3000 people [1]. Vibrio parahaemolyticus, a Gram-negative pathogenic halophilic bacterium, is the most prevalent in certain Asian areas and causes food poisoning, occasionally outbreaks, especially in the hot season [2].

V. parahaemolyticus is well known as the causative agent of the most prevalent food poisoning in Asia since the mackerel-borne outbreak in 1959 [3]. V. parahaemolyticus infections arise from the consumption of raw or undercooked seafood, typically causing gastroenteritis [4]. Infections can become systemic and can be life-threatening to immunocompromised individuals [5].

V. parahaemolyticus strains isolated from diarrheal patients produce thermostable direct hemolysin (TDH), TDH-related hemolysin (TRH), or both, while isolates from the environment rarely contain genes encoding these proteins. The TDH is the major virulence gene of V. parahaemolyticus and present in most of the clinical Kanagawa phenonemon (KP)-positive strains. In addition to its hemolytic activity, enterotoxicity, cytotoxicity as well as cardiotoxicity have been confirmed in TDH [6]. TDH-related hemolysin (TRH) is produced by KP-negative strains, encoded by trh gene, with about 67% amino acid sequence homology and immunologically related with TDH [7]. TRH is hemolytic, heat labile and is suspected to play an unclarified role in the diarrhea caused by these KP-negative strains. After the bacteria produce TDH and lipopolysaccharides released from a dead cell, human host would show strong systemic responses and produce mucosal B-cells against these antigens. However, recent study has also revealed that V. parahaemolyticus could profoundly disturb epithelial barrier function in Caco-2 cells with other virulence factors, in the absence of TDH [8]. V. parahaemolyticus also encodes two type III secretion systems (T3SS), which are located in different chromosomes. The contribution of these two T3SS has previously been characterized in animal models [9, 10]. These two T3SSs are involved in distinct pathogenic mechanisms during infection. The type III secretion system 1 (T3SS1) is required for cytotoxicity, whereas the type III secretion system 2 (T3SS2) is often responsible for enterotoxicity and intestinal fluid accumulation. The overall mechanism of pathogenesis in V. parahaemolyticus remains unclear.

V. parahaemolyticus infections are associated with multiple serotypes. Analysis of an outbreak in Calcutta in 1996 has identified a unique serotype clone, O3:K6, which is not previously isolated in that area [11]. After this report, O3:K6 isolates and its serovariants, defined to have identical genotype and molecular characteristics to those isolated in Calcutta, have been reported from foodborne outbreaks across the world. These serotypes were thought to be clonal derivatives of the O3:K6 serotype [12]. Recent multilocus sequence typing (MLST) data have further confirmed that multiple serotypes occur in a single genetic lineage [13].

Comparative genomics analysis of merely clinical strains of V. parahaemolyticus was described before [14]. However, the phylogenetic diversity and niche adaptation of V. parahaemolyticus at the species level is currently unknown. Strains isolated from environmental reservoirs could be as virulent as clinical strains in animal models [15], and the existing regulatory mechanisms of environmental V. parahaemolyticus could favor regulations of foreign virulence genes [16]. Nonetheless, the virulence potential of V. parahaemolyticus at the genome-wide level has not been described before. We sequenced the genomes of four O3:K6 and its serovariant clinical strains and five environmental strains of V. parahaemolyticus in Hong Kong. Genomes of a few environmental V. parahaemolyticus isolates have been announced recently [1720], however, comparative genomics study has not been reported until now. Phylogenomic analysis of our local strains and other available genomes revealed a clear distinction between clinical and environmental strains.

Results and discussion

Summary of genome sequencing data

As summarized in Table 1, we have completed 454 genome sequencing of nine V. parahaemolyticus isolates obtained in Hong Kong, including four clinical and five environmental strains. The average sequencing coverage was between 44× to 66×. Genome sequences of these local V. parahaemolyticus strains were de novo assembled using Newbler 2.7 (Roche Diagnostics). Genome assembly of V. parahaemolyticus yielded 130 to 1205 contigs of N50 contig length ranging from 7 kb to 305 kb. The clinical strains we sequenced were new O3:K6 and its serovariant (O3:K59), and they were positive for the tdh gene but negative for the trh gene. Our environmental strains were all tdh and trh negative. We have also included in our analyses six publicly available V. parahaemolyticus genomes, which are all clinical strains (Table 1) and two of them (AQ3810, AQ4037) pre-pandemic strains [14].

Table 1 List of V. parahaemolyticus clinical and environmental isolates analyzed in this study

Phylogenetic relationships

Phylogenetic relationships among our V. parahaemolyticus isolates and the reference strains (Table 1) were examined using a genome-wide approach based on 169,998 single nucleotide polymorphisms (SNPs). Eighty nine percent of these SNPs were in coding regions, of which 25.8% were missense, 0.6% was nonsense, and 73.6% were synonymous. Genome-wide comparison readily resolved the relationship among the clinical and environmental isolates (Figure 1). Clinical strains were found to be more related compared to the environmental strains. Moreover, our local clinical isolates were highly similar to RIMD 2210633 isolated from Thailand in 1996, consistent with previous results using various molecular methods [21][22] and further confirm clonality of new O3:K6 and its serovariants.

Figure 1
figure 1

Neighbor-joining phylogenetic tree showing relationships among V. parahaemolyticus isolates inferred using SNP sites. Nodal supports were calculated from 500 bootstrap pseudoreplicates. The tree was rooted using Vibrio harveyi ATCC BAA-1116. The scale bar represents 0.002 substitutions per nucleotide position. Red color indicates clinical strains, and grey color indicates environmental strains.

Pan-genome analysis

Comparative analysis of nine newly sequenced genomes and six publicly available genomes of V. parahaemolyticus could help us determine the global gene repertoire of the species. This can be described by its ‘pan-genome’ that includes a core genome containing genes present in all strains and a dispensable genome composed of genes absent from one or more strains and genes that are unique to some strains. Here, the core and pan-genomes of V. parahaemolyticus were identified using OrthoMCL [23].

Our result shows that the number of pan-genome gene families increased with the number of genomes analysed, indicating that V. parahaemolyticus harbors an open pan-genome (Figure 2). We have also analyzed the trend of new gene families, and we found that new genes will continue to be found with increasing of much more genomes.

Figure 2
figure 2

Gene repertoire analysis of V. parahaemolyticus. (A) Pan-genome and core genome size accumulation. (B) New gene family accumulation. The green bars indicate the number of expected new gene family detected for a particular number of genomes analysed, and the orange line indicates the trend of expected new gene family with an increasing number of genomes.

Gene-content analysis of V. parahaemolyticus

By investigating the presence and absence of genes in V. parahaemolyticus, we found that most regions within their genomes were conserved (Figure 3). We then further analyzed 24 genomic regions that are unique to V. parahaemolyticus RIMD 2210633 [24], which mainly include two T3SS, seven pathogenic islands, and f237 prophages. Comparative genomics revealed that these 24 regions accounted for most of the varying regions. Five pathogenic islands (VPal-1, VPal-3, VPal-4, VPal-5, and VPal-6) were found absent in not only the environmental strains, but also two pre-pandemic clinical strains. VPal-7 was found absent in AQ4037 and other environmental strains.

Figure 3
figure 3

BRIG visualization of V. parahaemolyticus genomes. The innermost circles represent the reference sequence of V. parahaemolyticus 2210633. Outer rings illustrate shared identity with other V. parahaemolyticus isolates.

Previously, microarray-based approach failed to detect genes responsible for pathogenesis of V. parahaemolyticus[25]. Traditional features also failed to differentiate clinical strains from environmental strains [26, 27]. Here, analysis of these varying regions revealed the presence of only 11 regions in the clinical group, which were absent in the environmental group (Table 2). As expected, genes in T3SS2 were not found in all of the clinical strains, which further confirm that it is not reliable to simply use tdh, trh or genes in T3SS as virulence features [25]. However, a region located in chromosome II (984909–999901) contained genes from VPA0950 to VPA0957, which encode putative biofilm associated proteins and outer membrane proteins. A few genes including two type IV pili were previously described in biofilm formation of V. parahaemolyticus, however details of the process are still unexplored in this bacterium [28, 29]. In other Vibrio species such as V. fischeri, biofilm formation was found to play an important role in host colonization [30]. It is suggested that gene clusters found in clinical strains of V. parahaemolyticus could also provide this bacterium with advantages in host adaptation. These gene clusters have not been found previously, and we suggest that they could act as new candidate virulence markers.

Table 2 Genomic regions that found in clinical group only

A novel genomic island commonly found in clinical strains

O3:K59 is a new serovariant of serotype O3:K6. It is recently reported that it could replace local serovariants and has caused pandemics in Chile and Peru [31, 32]. Here we analyzed the genome of our O3:K59 strain VIP4-0407 and found a novel genomic island that is absent in V. parahaemolyticus RIMD 2210633 [33]. This genomic island was found in most of our clinical strains including AQ4037, K5030, AN-5034, Peru-466, VIP4-0439, VIP4-0395, VIP4-0407 and AQ3810, but was absent in our environmental strains. We defined this genomic island as Vibrio parahaemolyticus island-8 (VPal-8). It contained six CDSs (Table 3), one of which showing high similarity with an AraC family transcriptional regulator that caused virulence in a mouse infection model in Mycobacterium tuberculosis[34]. ArcC family transcriptional regulators could regulate T3SS genes to modulate bacterial virulence [35], We analyzed sequences in this genomic island to determine whether there is any T3SS secreted protein. Using a combination of available T3SS prediction tools including T3_MM [36], Effective T4 [37], and T3SS effector prediction [38], we found that the third protein in this genomic island showed positive results in all of the prediction tools.

Table 3 Gene clusters in small genomic island found in V. parahaemolyticus

We further investigated whether VPal-8 is widespread across the population of clinical strains. Using this genomic island as a probe to perform southern hybridization on an additional population of 40 clinical strains, we found that this novel genomic island can be found in 29 of the clinical strains (Additional file 1: Table S1). Therefore, we suggest that the genomic island found in this study may be related to the virulence of V. parahaemolyticus.

sRNAs and CRISPR element analysis

CRISPRs (Clustered regularly interspaced short palindromic repeats) could play important roles in the interaction of bacteria and mobile genetic elements [39]. We annotated CRISPR elements in V. parahaemolyticus using CRISPR finder [40] and found a total of six CRISPR elements (Table 4). Our environmental strains were found to harbor fewer CRISPR types, whereas at least two CRISPRs were found in each of the clinical strains.

Table 4 Distribution of CRISPRs in V. parahaemolyticus strains

sRNAs are a class of non-coding RNAs in bacteria and are important post-transcriptional regulators in multiple crucial biological processes such as biofilm formation, quorum sensing and virulence [41, 42] . We determined the sRNA sequences in our V. parahaemolyticus strains and found that most sRNAs were highly conserved (Additional file 2: Table S2). However, a few sRNAs such as GcvB and STnc1460 were only present in a few strains, indicating that these sRNAs may mediate some strain-specific regulations.

Biphenyl degradation pathway identified in environmental strains

We annotated the subsystems of V. parahaemolyticus genomes using the RAST Server [43] (Figure 4). Five subsystems were found enriched in the environmental strains, including Metabolism of Aromatic Compounds (p = 1.12E06), Protein metabolism (p = 0.01), Iron acquisition and metabolism (p = 0.027), Phages and Prophages (p = 0.041) and Sulfur Metabolism (p = 0.046). The most enriched subsystem “Metabolism of Aromatic Compounds” contained an average of six genes in clinical strains, but an average of 10 genes in the environmental strains.

Figure 4
figure 4

Subsystem annotation of V. parahaemolyticus. The y-axis represents subsystems annotated in the RAST server whereas the x-axis shows the respective -log(P) values. P-values were calculated to examine if there is significant difference in the numbers between the subsystems of environmental strains and clinical strains. The subsystem “metabolism of aromatic compounds” was shown to have a most significant P-value (P-value <0.05, −logP >1.3).

We further examined this category to investigate which pathway has contributed most of the differences. Genes involved in biphenyl degradation pathway (BphE1, BphC, Bphj2) were found in all of our environmental strains but not in the clinical strains, except the pre-pandemic strain AQ3810 (Table 5). The presence of the unique biphenyl degradation pathway in our environmental strains is possibly caused by the widespread of polychlorinated biphenyls in the environment [44].

Table 5 Detection of biphenyl degradation pathway genes

Virulence potential of environmental strains

V. parahaemolyticus harbors many virulence factors including TDH, TRH and two T3SSs, which are generally pathogenic. Here we report the distribution of virulence factors and prophages in the genomes of our V. parahaemolyticus strains to reveal their pathogenic potential.

We predicted prophage regions in V. parahaemolyticus using the PHAST webserver [45] and found that clinical strains were restricted to a few prophage elements. The pandemic group of clinical strains possessed a filamentous prophage f237, which could increase virulence by adhering to intestinal cells [46]. A similar prophage was also found in the AQ3047 strain. In contrast, many diverse prophage elements were found in our environmental strains (Table 6). Six new prophages were found in four of our environmental strains, which harbor genes that could enhance the fitness for survival in specific environments. We further examined the relationship between the number of prophages and CRISPR in our strains. We analyzed the spacer sequences of CRISPRs and eight unique spacers sequence were found. However, we could not find any homologous prophage regions that could match with these spacer sequences in the respective strains. This could be explained by the possibility that the CRISPRs of Vibrio parahaemolyticus could inhibit the prophage insertion to the host genome. Indeed, an inverse correlation of CRISPRs with the number of prophage was also described in Streptococcus pyogenes[47].

Table 6 Distribution of prophages in V. parahaemolyticus strains

toxR gene, located next to toxS gene, involves in the regulation of many virulence-associated genes of V. parahaemolyticus[48]. Variation in the toxR sequence can be used to differentiate phylogenetically distinct clusters in V. parahaemolyticus. Seven base positions were previously identified to distinguish the O3:K6 isolates before 1995 from the new O3:K6 clones within a 1346-bp region [49]. We further investigated the variations between these base positions within our V. parahaemolyticus isolates. In all of the seven new O3:K6 clones, the bases were perfectly mapped to the reference ones (Table 7). However, another environmental strain (VIP4-0443) also had the same sequence bases with the new O3:K6 clones, suggesting that it may also harbor virulence potential. We then further examined whether there are any other virulence proteins in this strain. A total of 264 candidate virulence genes were found (Additional file 3: Table S3), which further confirm that this environmental strain could have virulence potential.

Table 7 Base variations in the toxRS sequence of selected V. parahaemolyticus strains

Construction of a genome sequence database for V. parahaemolyticus

A V. parahaemolyticus genome database was constructed here based on the Ensembl genome annotation system using Perl scripts (Figure 5). Genome DNA sequences can be input into the platform via a user-friendly web-based interface. Our web-based database allows comparative analysis and data mining using available V. parahaemolyticus genome data. A suite of useful computational tools is now available for data analysis such as MLST and detection of genetic variations. The web-based database is now launched as a publicly accessible domain (

Figure 5
figure 5

Architecture of the Ensembl-based genome database for V. parahaemolyticus.


This is the first study on the genome dynamics of V. parahaemolyticus at the species level. We have three main conclusions. First, by comparing the genome sequences, we found that our clinical strains were phylogenetically distinct from the environmental strains. We discovered a new gene cluster belonging to the biofilm associated proteins of V. parahaemolyticus. We suggest that this novel gene cluster is a new virulence marker differentiating clinical from environmental strains. We also found a novel genomic island (VPal-8) that frequently distributed in the clincial strains. Second, by analyzing the virulence features and prophage elements of our environmental strains on a genome-wide scale, we discovered a new type of toxRS in one of the environmental strains and a diverse array of prophage elements in the environmental strains. These results indicate that the environmental strains studied could have virulence potential. Third, we set up an Ensembl-based genome database for V. parahaemolyticus to provide a unified and user-friendly platform to access all currently available genome sequences and annotations of this bacterium.


Ethics statement and isolate collection

Nine isolates of V. parahaemolyticus were selected from the bacterial archive maintained at the microbiology laboratory of the Department of Health of Hong Kong and were characterized for genome sequencing. The study was approved by the Joint CUHK-NTEC Clinical Research Ethical Committee at the Prince of Wales Hospital in Hong Kong. Written informed consent was obtained from all the studied subjects for sample collection and subsequent analysis.

DNA sequencing

Whole-genome shotgun sequencing was performed using a 454 Genome Sequencer (GS) FLX-Titanium (Roche Diagnostics, US). DNA of each selected isolate was extracted from culture and subjected to sequencing library preparation according to the manufacturer’s recommended protocols. Libraries of processed genomic DNA fragments immobilized on DNA capture beads were individually sequenced on a PicoTiterPlate device. Bases sequenced and the corresponding quality values were called and delivered in a standard format by GS-FLX-Titanium system for downstream bioinformatic analyses. The remaining genome gaps were filled by the GS-FLX-provided paired-end approach and/or primer walking. The primer walking approach involved design of primers flanking the gaps and PCR to amplify the DNA fragments covering the gaps. PCR products were then directly sequenced using the ABI DNA Sequencer (Applied Biosystems, US).

Sequence assembly and annotation

Raw sequence reads were first filtered to remove low-quality reads. Sequence assembly was performed on the remaining clean reads de novo using the GS Newbler 2.7. The genome sequences were annotated using the RAST webserver. Contigs were reordered and genome comparisons were carried out using mauve program [50]. BRIG was used for genome alignment visualization [51]. Prophage elements of V. parahaemolyticus were predicted using PHAST [45].

Phylogenetic tree construction

SNPs among our isolates and several foreign strains were identified using GS Reference Mapper (Roche Diagnostics). SNPs detected among all strains were concatenated together and a neighbor-joining phylogenetic tree was constructed using MEGA 6 [52] with 500 bootstrap pseudoreplicates. Vibrio harveyi ATCC BAA-1116 was used as outgroup for defining the root.

Small RNAs and CRISPRs

A list of small RNA sequences was retrieved from the BSRD database [42]. sRNA homologs in our V. parahaemolyticus isolates were then identified using BLASTn. The e-value of BLAST was set to 1e-5. Annotation of CRISPR elements was done using CRISPRfinder [40].

Experimental validation of genomic island

Two sets of primers were designed for PCR-based detection of genomic island in a population of clinical V. parahaemolyticus. These clinical strains were previously used for microarray-based analysis [53].

Availability of supporting data

Data of this Whole Genome Shotgun project have been deposited at GenBank under the accessions AXNJ00000000, AXNK00000000, AXNL00000000, AXNM00000000, AXNO00000000, AXNP00000000, AXNQ00000000, AXNR00000000, and AXNS00000000. The phylogenetic tree and associated data matrix are available in TreeBASE database (Accession URL:


  1. Scientific Committee on Enteric Infections and Foodborne Diseases. 2008. Foodborne illness - Intersection between clinical and public health approaches, Centre for Health Protection, Hong Kong SAR.,

  2. Nair GB, Ramamurthy T, Bhattacharya SK, Dutta B, Takeda Y, Sack DA: Global dissemination of Vibrio parahaemolyticus serotype O3:K6 and its serovariants. Clin Microbiol Rev. 2007, 20: 39-48. 10.1128/CMR.00025-06.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  3. Miyamoto Y, Nakamura K, Takizawa K: Seasonal distribution of oceanomonas spp., halophilic bacteria, in the coastal sea. its significance in epidemiology and marine industry. Jpn J Microbiol. 1962, 6: 141-158. 10.1111/j.1348-0421.1962.tb00231.x.

    Article  Google Scholar 

  4. Miyamoto Y, Kato T, Obara Y, Akiyama S, Takizawa K, Yamai S: In vitro hemolytic characteristic of Vibrio parahaemolyticus: its close correlation with human pathogenicity. J Bacteriol. 1969, 100: 1147-1149.

    CAS  PubMed Central  PubMed  Google Scholar 

  5. Baker-Austin C, Stockley L, Rangdale R, Martínez-Urtaza J: Environmental occurrence and clinical impact of Vibrio vulnificus and Vibrio parahaemolyticus: a European perspective. Environ Microbiol Rep. 2010, 2: 7-18. 10.1111/j.1758-2229.2009.00096.x.

    Article  PubMed  Google Scholar 

  6. Raimondi F, Kao JP, Fiorentini C, Fabbri A, Donelli G, Gasparini N, Rubino A, Fasano A: Enterotoxicity and cytotoxicity of Vibrio parahaemolyticus thermostable direct hemolysin in in vitro systems. Infect Immun. 2000, 68: 3180-3185. 10.1128/IAI.68.6.3180-3185.2000.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Honda T, Ni YX, Miwatani T: Purification and characterization of a hemolysin produced by a clinical isolate of Kanagawa phenomenon-negative Vibrio parahaemolyticus and related to the thermostable direct hemolysin. Infect Immun. 1988, 56: 961-965.

    CAS  PubMed Central  PubMed  Google Scholar 

  8. Lynch T, Livingstone S, Buenaventura E, Lutter E, Fedwick J, Buret AG, Graham D, DeVinney R: Vibrio parahaemolyticus disruption of epithelial cell tight junctions occurs independently of toxin production. Infect Immun. 2005, 73: 1275-1283. 10.1128/IAI.73.3.1275-1283.2005.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  9. Piñeyro P, Zhou X, Orfe LH, Friel PJ, Lahmers K, Call DR: Development of two animal models to study the function of Vibrio parahaemolyticus type III secretion systems. Infect Immun. 2010, 78: 4551-4559. 10.1128/IAI.00461-10.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Hiyoshi H, Kodama T, Iida T, Honda T: Contribution of Vibrio parahaemolyticus virulence factors to cytotoxicity, enterotoxicity, and lethality in mice. Infect Immun. 2010, 78: 1772-1780. 10.1128/IAI.01051-09.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Okuda J, Ishibashi M, Hayakawa E, Nishino T, Takeda Y, Mukhopadhyay AK, Garg S, Bhattacharya SK, Nair GB, Nishibuchi M: Emergence of a unique O3:K6 clone of Vibrio parahaemolyticus in Calcutta, India, and isolation of strains from the same clonal group from southeast asian travelers arriving in Japan. J Clin Microbiol. 1997, 35: 3150-3155.

    CAS  PubMed Central  PubMed  Google Scholar 

  12. Chowdhury NR, Chakraborty S, Ramamurthy T, Nishibuchi M, Yamasaki S, Takeda Y, Nair GB: Molecular evidence of clonal Vibrio parahaemolyticus pandemic strains. Emerg Infect Dis. 2000, 6: 631-636. 10.3201/eid0606.000612.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Chowdhury NR, Stine OC, Morris JG, Nair GB: Assessment of evolution of pandemic Vibrio parahaemolyticus by multilocus sequence typing. J Clin Microbiol. 2004, 42: 1280-1282. 10.1128/JCM.42.3.1280-1282.2004.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Chen Y, Stine OC, Badger JH, Gil AI, Nair GB, Nishibuchi M, Fouts DE: Comparative genomic analysis of Vibrio parahaemolyticus: serotype conversion and virulence. BMC Genomics. 2011, 12: 294-10.1186/1471-2164-12-294.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. Chatzidaki-Livanis M, Hubbard MA, Gordon K, Harwood VJ, Wright AC: Genetic distinctions among clinical and environmental strains of Vibrio vulnificus. Appl Environ Microbiol. 2006, 72: 6136-6141. 10.1128/AEM.00341-06.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Mahoney JC, Gerding MJ, Jones SH, Whistler CA: Comparison of the pathogenic potentials of environmental and clinical Vibrio parahaemolyticus strains indicates a role for temperature regulation in virulence. Appl Environ Microbiol. 2010, 76: 7459-7465. 10.1128/AEM.01450-10.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Liu M, Chen S: Draft genome sequence of Vibrio parahaemolyticus V110, isolated from shrimp in Hong Kong. Genome Announc. 2013, 1: e00300-e00313.

    PubMed Central  PubMed  Google Scholar 

  18. Jun JW, Kim JH, Choresca CH, Shin SP, Han JE, Park SC: Draft genome sequence of Vibrio parahaemolyticus SNUVpS-1 isolated from Korean Seafood. Genome Announc. 2013, 1: doi:10.1128/genomeA.00132-12. Epub 2013 Feb 7

    Google Scholar 

  19. Kalburge SS, Polson SW, Boyd Crotty K, Katz L, Turnsek M, Tarr CL, Martinez-Urtaza J, Boyd EF: Complete genome sequence of Vibrio parahaemolyticus environmental strain UCM-V493. Genome Announc. 2014, 2: doi:10.1128/genomeA.00159-14

    Google Scholar 

  20. Kumar BK, Deekshit VK, Rai P, Gurtler V, Karunasagar I, Karunasagar I: Draft genome sequence of trh + Vibrio parahaemolyticus VP-49, isolated from seafood harvested along the Mangalore Coast, India. Genome Announc. 2014, 2: e00607-e00614.

    PubMed Central  PubMed  Google Scholar 

  21. Wong HC, Liu SH, Wang TK, Lee CL, Chiou CS, Liu DP, Nishibuchi M, Lee BK: Characteristics of Vibrio parahaemolyticus O3:K6 from Asia. Appl Environ Microbiol. 2000, 66: 3981-3986. 10.1128/AEM.66.9.3981-3986.2000.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Matsumoto C, Okuda J, Ishibashi M, Iwanaga M, Garg P, Rammamurthy T, Wong HC, DePaola A, Kim YB, Albert MJ, Nishibuchi M: Pandemic spread of an O3:K6 clone of Vibrio parahaemolyticus and emergence of related strains evidenced by arbitrarily primed PCR and toxRS sequence analyses. J Clin Microbiol. 2000, 38: 578-585.

    CAS  PubMed Central  PubMed  Google Scholar 

  23. Li L: OrthoMCL: identification of Ortholog groups for Eukaryotic Genomes. Genome Res. 2003, 13: 2178-2189. 10.1101/gr.1224503.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Boyd EF, Cohen A, Naughton LM, Ussery DW, Binnewies TT, Stine OC, Parent MA: Molecular analysis of the emergence of pandemic Vibrio parahaemolyticus. BMC Microbiol. 2008, 8: 110-10.1186/1471-2180-8-110.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Izutsu K, Kurokawa K, Tashiro K, Kuhara S, Hayashi T, Honda T, Iida T: Comparative genomic analysis using microarray demonstrates a strong correlation between the presence of the 80-kilobase pathogenicity island and pathogenicity in Kanagawa phenomenon-positive Vibrio parahaemolyticus strains. Infect Immun. 2008, 76: 1016-1023. 10.1128/IAI.01535-07.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Tsai S-E, Jong K-J, Tey YH, Yu W-T, Chiou C-S, Lee Y-S, Wong H-C: Molecular characterization of clinical and environmental Vibrio parahaemolyticus isolates in Taiwan. Int J Food Microbiol. 2013, 165: 18-26. 10.1016/j.ijfoodmicro.2013.04.017.

    Article  CAS  PubMed  Google Scholar 

  27. Yu W-T, Jong K-J, Lin Y-R, Tsai S-E, Tey YH, Wong H-C: Prevalence of Vibrio parahaemolyticus in oyster and clam culturing environments in Taiwan. Int J Food Microbiol. 2013, 160: 185-192. 10.1016/j.ijfoodmicro.2012.11.002.

    Article  PubMed  Google Scholar 

  28. Enos-Berlage JL, Guvener ZT, Keenan CE, McCarter LL: Genetic determinants of biofilm development of opaque and translucent Vibrio parahaemolyticus. Mol Microbiol. 2005, 55: 1160-1182.

    Article  CAS  PubMed  Google Scholar 

  29. Shime-Hattori A, Iida T, Arita M, Park K-S, Kodama T, Honda T: Two type IV pili of Vibrio parahaemolyticus play different roles in biofilm formation. FEMS Microbiol Lett. 2006, 264: 89-97. 10.1111/j.1574-6968.2006.00438.x.

    Article  CAS  PubMed  Google Scholar 

  30. Yildiz FH, Visick KL: Vibrio biofilms: so much the same yet so different. Trends Microbiol. 2009, 17: 109-118. 10.1016/j.tim.2008.12.004.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Velazquez-Roman J, León-Sicairos N, de Jesus Hernández-Díaz L, Canizalez-Roman A: Pandemic Vibrio parahaemolyticus O3:K6 on the American continent. Front Cell Infect Microbiol. 2013, 3: 110-

    PubMed Central  Google Scholar 

  32. García K, Torres R, Uribe P, Hernández C, Rioseco ML, Romero J, Espejo RT: Dynamics of clinical and environmental Vibrio parahaemolyticus strains during seafood-related summer diarrhea outbreaks in southern Chile. Appl Environ Microbiol. 2009, 75: 7482-7487. 10.1128/AEM.01662-09.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Makino K, Oshima K, Kurokawa K, Yokoyama K, Uda T, Tagomori K, Iijima Y, Najima M, Nakano M, Yamashita A, Kubota Y, Kimura S, Yasunaga T, Honda T, Shinagawa H, Hattori M, Iida T: Genome sequence of Vibrio parahaemolyticus: a pathogenic mechanism distinct from that of V cholerae. Lancet. 2003, 361: 743-749. 10.1016/S0140-6736(03)12659-1.

    Article  CAS  PubMed  Google Scholar 

  34. Frota CC, Papavinasasundaram KG, Davis EO, Colston MJ: The AraC family transcriptional regulator Rv1931c plays a role in the virulence of Mycobacterium tuberculosis. Infect Immun. 2004, 72: 5483-5486. 10.1128/IAI.72.9.5483-5486.2004.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Plano GV: Modulation of AraC family member activity by protein ligands. Mol Microbiol. 2004, 54: 287-290. 10.1111/j.1365-2958.2004.04306.x.

    Article  CAS  PubMed  Google Scholar 

  36. Wang Y, Sun M, Bao H, White AP: T3_MM: a Markov model effectively classifies bacterial type III secretion signals. PLoS One. 2013, 8: e58173-10.1371/journal.pone.0058173.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Arnold R, Brandmaier S, Kleine F, Tischler P, Heinz E, Behrens S, Niinikoski A, Mewes H-W, Horn M, Rattei T: Sequence-based prediction of type III secreted proteins. PLoS Pathog. 2009, 5: e1000376-10.1371/journal.ppat.1000376.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Löwer M, Schneider G: Prediction of type III secretion signals in genomes of gram-negative bacteria. PLoS One. 2009, 4: e5917-10.1371/journal.pone.0005917.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Westra ER, Swarts DC, Staals RHJ, Jore MM, Brouns SJJ, van der Oost J: The CRISPRs, they are a-changin’: how prokaryotes generate adaptive immunity. Annu Rev Genet. 2012, 46: 311-339. 10.1146/annurev-genet-110711-155447.

    Article  CAS  PubMed  Google Scholar 

  40. Grissa I, Vergnaud G, Pourcel C: CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007, 35: W52-W57. 10.1093/nar/gkm360.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Storz G, Vogel J, Wassarman KM: Regulation by small RNAs in bacteria: expanding frontiers. Mol Cell. 2011, 43: 880-891. 10.1016/j.molcel.2011.08.022.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Li L, Huang D, Cheung MK, Nong W, Huang Q, Kwan H-S: BSRD: a repository for bacterial small regulatory RNA. Nucleic Acids Res. 2013, 41: D233-D238. 10.1093/nar/gks1264.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, Edwards RA, Gerdes S, Parrello B, Shukla M, Vonstein V, Wattam AR, Xia F, Stevens R: The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014, 42: D206-D214. 10.1093/nar/gkt1226.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. Abramowicz DA: Aerobic and anaerobic PCB biodegradation in the environment. Environ Health Perspect. 1995, 103: 97-99. 10.1289/ehp.95103s297.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS: PHAST: a fast phage search tool. Nucleic Acids Res. 2011, 39: W347-W352. 10.1093/nar/gkr485.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Nasu H, Iida T, Sugahara T, Yamaichi Y, Park KS, Yokoyama K, Makino K, Shinagawa H, Honda T: A filamentous phage associated with recent pandemic Vibrio parahaemolyticus O3:K6 strains. J Clin Microbiol. 2000, 38: 2156-2161.

    CAS  PubMed Central  PubMed  Google Scholar 

  47. Nozawa T, Furukawa N, Aikawa C, Watanabe T, Haobam B, Kurokawa K, Maruyama F, Nakagawa I: CRISPR inhibition of prophage acquisition in Streptococcus pyogenes. PLoS One. 2011, 6: e19543-10.1371/journal.pone.0019543.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Whitaker WB, Parent MA, Boyd A, Richards GP, Boyd EF: The Vibrio parahaemolyticus ToxRS regulator is required for stress tolerance and colonization in a novel orogastric streptomycin-induced adult murine model. Infect Immun. 2012, 80: 1834-1845. 10.1128/IAI.06284-11.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Okura M, Osawa R, Iguchi A, Arakawa E, Terajima J: Genotypic analyses of Vibrio parahaemolyticus and development of a pandemic group-specific multiplex PCR assay. J Clin Microbiol. 2003, 41: 4676-4682. 10.1128/JCM.41.10.4676-4682.2003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  50. Darling AE, Mau B, Perna NT: progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010, 5: e11147-10.1371/journal.pone.0011147.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Alikhan N-F, Petty NK, Ben Zakour NL, Beatson SA: BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011, 12: 402-10.1186/1471-2164-12-402.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  52. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013, 30: 2725-2729. 10.1093/molbev/mst197.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  53. Han H, Wong H-C, Kan B, Guo Z, Zeng X, Yin S, Liu X, Yang R, Zhou D: Genome plasticity of Vibrio parahaemolyticus: microevolution of the “pandemic group.”. BMC Genomics. 2008, 9: 570-10.1186/1471-2164-9-570.

    Article  PubMed Central  PubMed  Google Scholar 

Download references


This work was supported by the Research Fund for the Control of Infectious Diseases (CHP-PH-06) from the Food and Health Bureau of the Hong Kong SAR, ROC.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Hoi Shan Kwan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the project: HSK. Analyzed and interpreted the data: LL, HCW. Performed the genomic island experiments: HCW. Constructed the genome database: WN. Performed 454 sequencing: PTWL. Provided the bacterial strains: KMK. Drafted the manuscript: LL. Revised the manuscript: LL, HSK, HCW, MKC. All authors have read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, L., Wong, Hc., Nong, W. et al. Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus. BMC Genomics 15, 1135 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: