- Research article
- Open Access
Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus
BMC Genomicsvolume 15, Article number: 1135 (2014)
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 (http://kwanlab.bio.cuhk.edu.hk/vp) 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 . 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 .
V. parahaemolyticus is well known as the causative agent of the most prevalent food poisoning in Asia since the mackerel-borne outbreak in 1959 . V. parahaemolyticus infections arise from the consumption of raw or undercooked seafood, typically causing gastroenteritis . Infections can become systemic and can be life-threatening to immunocompromised individuals .
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 . 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 . 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 . 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 . 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 . Recent multilocus sequence typing (MLST) data have further confirmed that multiple serotypes occur in a single genetic lineage .
Comparative genomics analysis of merely clinical strains of V. parahaemolyticus was described before . 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 , and the existing regulatory mechanisms of environmental V. parahaemolyticus could favor regulations of foreign virulence genes . 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 [17–20], 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 .
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  and further confirm clonality of new O3:K6 and its serovariants.
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 .
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.
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 , 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.
Previously, microarray-based approach failed to detect genes responsible for pathogenesis of V. parahaemolyticus. 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 . 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 . 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.
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 . 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. ArcC family transcriptional regulators could regulate T3SS genes to modulate bacterial virulence , 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 , Effective T4 , and T3SS effector prediction , we found that the third protein in this genomic island showed positive results in all of the prediction tools.
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 . We annotated CRISPR elements in V. parahaemolyticus using CRISPR finder  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.
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  (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.
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 .
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  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 . 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.
toxR gene, located next to toxS gene, involves in the regulation of many virulence-associated genes of V. parahaemolyticus. 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 . 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.
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 (http://kwanlab.bio.cuhk.edu.hk/vp).
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.
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 . BRIG was used for genome alignment visualization . Prophage elements of V. parahaemolyticus were predicted using PHAST .
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  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 . 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 .
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 .
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: http://purl.org/phylo/treebase/phylows/study/TB2:S16780).
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. http://www.chp.gov.hk/files/pdf/foodborne_illness-intersection_between_clinical_and_public_health_approaches_r.pdf,
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Liu M, Chen S: Draft genome sequence of Vibrio parahaemolyticus V110, isolated from shrimp in Hong Kong. Genome Announc. 2013, 1: e00300-e00313.
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
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
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.
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.
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.
Li L: OrthoMCL: identification of Ortholog groups for Eukaryotic Genomes. Genome Res. 2003, 13: 2178-2189. 10.1101/gr.1224503.
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.
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.
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.
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.
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.
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.
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.
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-
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Abramowicz DA: Aerobic and anaerobic PCB biodegradation in the environment. Environ Health Perspect. 1995, 103: 97-99. 10.1289/ehp.95103s297.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The authors declare that they have no competing interests.
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.