Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Applied and Environmental Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Genetics and Molecular Biology

Genomic Features of Environmental and Clinical Vibrio parahaemolyticus Isolates Lacking Recognized Virulence Factors Are Dissimilar

J. Ronholm, N. Petronella, C. Chew Leung, A. W. Pightling, S. K. Banerjee
H. L. Drake, Editor
J. Ronholm
aMicrobiology Research Division, Bureau of Microbial Hazards, Food Directorate, Health Products and Food Branch, Health Canada, Ottawa, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
N. Petronella
bBiostatistics and Modelling Division, Bureau of Food Surveillance and Science Integration, Food Directorate, Health Products and Food Branch, Health Canada, Ottawa, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
C. Chew Leung
aMicrobiology Research Division, Bureau of Microbial Hazards, Food Directorate, Health Products and Food Branch, Health Canada, Ottawa, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
A. W. Pightling
aMicrobiology Research Division, Bureau of Microbial Hazards, Food Directorate, Health Products and Food Branch, Health Canada, Ottawa, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
S. K. Banerjee
aMicrobiology Research Division, Bureau of Microbial Hazards, Food Directorate, Health Products and Food Branch, Health Canada, Ottawa, Ontario, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
H. L. Drake
University of Bayreuth
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AEM.03465-15
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Vibrio parahaemolyticus is a bacterial pathogen that can cause illness after the consumption or handling of contaminated seafood. The primary virulence factors associated with V. parahaemolyticus illness are thermostable direct hemolysin (TDH) and Tdh-related hemolysin (TRH). However, clinical strains lacking tdh and trh have recently been isolated, and these clinical isolates are poorly understood. To help understand the emergence of clinical tdh- and trh-negative isolates, a genomic approach was used to comprehensively compare 4 clinical tdh- and trh-negative isolates with 16 environmental tdh- and trh-negative isolates and 34 clinical isolates positive for tdh or trh, or both, with the objective of identifying genomic features that are unique to clinical tdh- and trh-negative isolates. The prevalence of pathogenicity islands (PAIs) common to clinical isolates was thoroughly examined in each of the clinical tdh- and trh-negative isolates. The tdh PAI was not present in any clinical or environmental tdh- and trh-negative isolates. The trh PAI was not present in any environmental isolates; however, in clinical tdh- and trh-negative isolate 10-4238, the majority of the trh PAI including a partial trh1 gene was present, which resulted in reclassification of this isolate as a tdh-negative and trh-positive isolate. In the other clinical tdh- and trh-negative isolates, neither the trh gene nor the trh PAI was present. We identified 862 genes in clinical tdh- and trh-negative isolates but not in environmental tdh- and trh-negative isolates. Many of these genes are highly homologous to genes found in common enteric bacteria and included genes encoding a number of chemotaxis proteins and a novel putative type VI secretion system (T6SS) effector and immunity protein (T6SS1). The availability of genome sequences from clinical V. parahaemolyticus tdh- and trh-negative isolates and the comparative analysis may help provide an understanding of how this pathotype is able to survive in vivo during clinical illness.

INTRODUCTION

Vibrio parahaemolyticus is a Gram-negative, halophilic bacterium that is ubiquitous in marine and estuarine environments and is often found colonizing shellfish or shrimp. While most strains are nonpathogenic, many have acquired virulence factors that result in illness when individuals are exposed to V. parahaemolyticus strains carrying these virulence factors (1). V. parahaemolyticus is recognized to be a leading cause of foodborne illness worldwide and is transmitted via the handling and consumption of raw or undercooked contaminated seafood (1). Infections occur both sporadically and in very large outbreaks. The most common manifestation of V. parahaemolyticus infection is acute, watery diarrhea accompanied by abdominal pain and nausea, although symptoms can also be severe and include a dysentery-like illness or septicemia (2). Since most cases of illness caused by V. parahaemolyticus are self-limiting, rates of infection are probably underestimated due to underreporting.

Clinical V. parahaemolyticus isolates generally have at least one of two major toxigenic virulence factors, thermostable direct hemolysin (TDH) (3) and TDH-related hemolysin (TRH) (4). TDH has hemolytic activity on a blood-containing medium, Wagatsuma agar, and the process is referred to as the Kanagawa phenomenon (KP) (3). During infection, TDH is involved in cytotoxicity and hemolytic activity, and on the basis of the sequence similarity between TDH and TRH, TRH is believed to act similarly (5–7). The presence of tdh and/or trh is common in pathogenic isolates but relatively rare in environmental strains; therefore, the presence of these genes is used to assess the virulence potential of V. parahaemolyticus isolates (8, 9). In addition to tdh and trh, whole-genome sequencing (WGS) of V. parahaemolyticus led to the identification of two nonredundant type III secretion system (T3SS) gene clusters, dubbed T3SS1 and T3SS2, on chromosome 1 and chromosome 2, respectively, which are also involved in virulence (5, 10). T3SSs are needle-like apparatuses that inject bacterial effector proteins, such as toxins or hemolysins, directly through the membrane and into the cytoplasm of eukaryotic cells (11). On the basis of its G+C content and its high degree of sequence identity with the T3SSs of other Vibrio species, T3SS1 appears to have been ancestrally acquired and is present in all V. parahaemolyticus isolates, even nonpathogenic strains (10). Once it is in vivo, T3SS1 appears to inject effectors, such as VopQ, VopR, VopS, and VPA0450, directly into eukaryotic cells, resulting in cytotoxicity (5, 10). T3SS2 derives from two separate lineages: T3SS2α is typically found on a pathogenicity island (PAI) with the tdh gene, and T3SS2β is found with the trh gene. T3SS2 appears to inject VopA, VopC, VopL, and VopT into eukaryotic cells, resulting in both cytotoxicity and enterotoxicity (5, 12–14). T3SS2 has a lower G+C content than the genomic average, which is indicative of a recent acquisition of this region through lateral gene transfer (15).

Some V. parahaemolyticus strains also have two type VI secretion systems (T6SSs) (16). The T6SS was recently defined functionally in V. cholerae and is structurally similar to a contractile phage tail, but in the reverse orientation, that is fully assembled inside the bacterial cell and injects effectors directly into the recipient cell (17, 18). It is composed of 13 essential genes and a variable number of nonessential genes, including various effectors (19). In V. parahaemolyticus, T6SS2 is found in all strains, while T6SS1 is mostly associated with clinical isolates (20) and may play a role in virulence (21), although this has not yet been demonstrated conclusively. The T6SS1 also appears to have antibacterial activity in the environment, which may give isolates containing this system a competitive advantage (16).

Recent studies have reported the isolation of strains that lack both tdh and trh (tdh- and trh-negative strains) from clinical samples (22–27). However, the ability of tdh- and trh-negative strains to independently cause clinical illness is still controversial. For example, coinfection with multiple V. parahaemolyticus strains is known to occur, and if at least one infecting strain carries tdh or trh, it is possible that nonpathogenic tdh- and trh-negative strains could be isolated from a sick individual without having a direct involvement in illness (24). The coinfection model is supported by the finding that a single seafood sample is often contaminated by several different V. parahaemolyticus strains, some of which appear to be nonpathogenic (24). Alternatively, there are several lines of evidence that support the opinion that some tdh- and trh-negative isolates are able to induce clinical infection. During a coinfection study, three sick patients produced 30 tdh- and trh-negative isolates, and despite multiple culturing attempts, no other enteric pathogen or tdh- and trh-positive V. parahaemolyticus strain could be isolated from these patients (24). However, regardless of their independent pathogenicity, clinical tdh- and trh-negative isolates have a demonstrated ability to survive in vivo during illness, and we do not understand if the role of these isolates in human illness is as a causative agent, an innocent bystander, or an active participant in a multistrain infection. Therefore, this investigation was undertaken to better understand clinical tdh- and trh-negative isolates and how they compare on a genomic level to traditional pathogenic V. parahaemolyticus isolates and environmental tdh- and trh-negative isolates. Here we present a thorough comparative genomic analysis of multiple clinical tdh- and trh-negative isolates. This comprehensive approach has provided several insights into the pangenomics of clinical V. parahaemolyticus isolates and led to the identification of a novel putative T6SS effector.

MATERIALS AND METHODS

Genome sequencing, assembly, and annotation.Each of the clinical isolates sequenced in this study were Canadian clinical strains originating from provincial public health laboratories and submitted to the National Microbiology Laboratory (Public Health Agency of Canada), the British Columbia Centre for Disease Control (BCCDC), or the Bureau of Microbial Hazards (BMH) (Health Canada). All isolates were routinely propagated on TSA-2N agar (Difco BD, NJ, USA). The environmental isolates were each isolated from seafood by the Vibrio Reference Laboratory in Canada. DNA for whole-genome sequencing was extracted using a Maxwell 16 SEV cell DNA purification kit (Promega, Madison, WI). The short-read sequence data were generated by preparing a paired-end library with a Nextera XT DNA sample preparation kit (Illumina, San Diego, CA) and sequencing the library on a MiSeq benchtop sequencer (Illumina) for 500 or 600 cycles. Previously, Banerjee et al. (26) performed PCR amplification of the tdh and trh genes, and this allowed us to select isolates for WGS and group isolates according to clinical pathotype (tdh and trh positive, tdh negative and trh positive, tdh positive and trh negative, or tdh and trh negative). This grouping scheme was used during the analysis in our study.

Genomes were assembled de novo using the SPAdes (version 3.1.1) program, and quality was assessed using QUAST, as outlined by Ronholm et al. (Table 1) (47). De novo assembly resulted in the PAIs of interest, including the tdh PAI and the trh PAI, being split between contigs. For analysis of pathogenicity islands, raw MiSeq reads were reassembled by reference-guided assembly using the Burrows-Wheeler aligner (BWA; version 07.05a) (29) and then visualized and inspected using the Tablet (version 1.14.10.20) program (30).

View this table:
  • View inline
  • View popup
TABLE 1

Sequencing and annotation results for V. parahaemolyticus strains analyzed during this studya

Annotation of the function of the protein-coding sequences was performed using the Basic Local Alignment Search Tool (BLAST), which was used to compare the sequences obtained against the sequences in the COG (Clusters of Orthologous Groups) of proteins databases and the NCBI nr protein database.

Comparative genomics.An in-house Perl script was used to identify all orthologous, accessory, and unique genes between Vibrio genomes (the source code can be found at https://github.com/bfssi-nicholas-petronella/HARDCore.git.) Using the BLAST program, genes having 60% sequence identity over 80% of their length were considered orthologs (31, 32). The set of orthologous genes shared by all genomes was defined as the core genome, total genes identified within all genomes were defined as the pangenome, accessory genomes were defined as the set of genes possessed by a subset of genomes (pathotypes), and unique genomes were the subset of genes possessed only by a single strain's genome. To construct whole-genome rarefaction curves, an in-house script was also used in an iterative manner to obtain the number of genes contained in the pangenome relative to the number of genomes analyzed. For example, for the data points for the 16 genomes, the size of the pangenome was calculated for 16 randomly selected genomes (from the 38 possibilities). This was repeated 1,000 times using 16 genomes randomly selected from 38 genomes, and the average was taken. The purpose was to determine if V. parahaemolyticus had an open or closed pangenome and verify that the genomes from the number of strains selected for the experiment could adequately represent the core genome.

Phylogenetics.To construct phylogenetic trees for the clinical isolates, the largest representative sequence of each core gene was retrieved from the pangenome. The homologue for each representative sequence was then retrieved from the whole genome of each strain. The Prokka software tool was used to generate the corresponding amino acid sequence (33). Trees were constructed from the core genome amino acid sequence. Each gene was aligned using the Muscle (version 3.8.31) program with default parameters (34), and then RAxML software was used to construct a maximum likelihood phylogeny from concatenated alignments with 100 bootstrap values (35).

To construct phylogenetic trees for each isolate used in the study, the multilocus sequence typing (MLST) sequences were retrieved from the WGS. Trees were constructed from the MLST alleles. Each gene was aligned using the Muscle (version 3.8.31) program with default parameters (34), and then RAxML was used to construct a maximum likelihood phylogeny from concatenated alignments with 1,000 bootstrap values (35).

Genetic islands and secretion systems.The presence and absence of the tdh and trh pathogenicity islands and various secretion systems were investigated using BLASTn. The EasyFig (version 2.1) program was used for visualization of genomic islands (36). For larger genomic islands (tdh PAI and trh PAI), a Burrows-Wheeler transform reference-guided assembly was used on raw fastq data with the canonical sequence of the genomic island (29). The assembly was inspected for accuracy and validity using Tablet, prior to visualization (30). To possibly identify novel pathogenicity islands in the tdh- and trh-negative pathotype, the completed results of an IslandViewer analysis on strain CDC_K4557 were downloaded from the IslandViewer3 website (37). BLASTn was used to search each genome included in this study (Table 1) for each island that was identified.

Nucleotide sequence accession numbers.All nucleotide sequence data referred to in this article have been deposited in the DDBJ/EMBL/GenBank database under the various BioProject accession numbers provided in Table 1. Additional important data are included in the supplemental material.

RESULTS AND DISCUSSION

Type III secretion systems and tdh and trh pathogenicity islands.Traditionally, a pathogenic V. parahaemolyticus strain has been defined by the presence of tdh or trh, or both (5). These two virulence factors generally occur near a T3SS. Two T3SSs have been reported in V. parahaemolyticus, and two variants of T3SS2 (T3SS2α and T3SS2β) have been described (12). T3SS2α is associated with tdh, while T3SS2β is typically found with trh (14), though exceptions exist (27). T3SS1 is composed of 42 genes (VPA1656 to VAP1696 and VPA0450 in V. parahaemolyticus RIMD2210633 [15]). T3SS1 was present in each of the clinical and environmental isolates in this investigation. The tdh PAI (also known as VPAI-7) is composed of 87 coding sequences (in V. parahaemolyticus RIMD2210633) and includes tdhA, tdhS, and the T3SS2α genes (38). In our current work, due to multiple homologous areas in the tdh PAI and trh PAI areas, the use of de novo assembly led to this large PAI being split between multiple contigs. Therefore, to analyze this region properly, a reference-guided assembly based on the tdh PAI sequence (VPA1310 to VPA1396) from V. parahaemolyticus RIMD2210663 was used (39). Using this method, a complete tdh PAI including T3SS2α was consistently observed in all clinical tdh-positive and trh-negative isolates (Fig. 1A). Clinical tdh- and trh-positive isolates contained homologues to some coding sequences typically found in this island (VPA1310, VPA1311, VPA1313 to VPA1318, VPA1320, VPA1321, VPA1329, VPA1342, VPA1347, and VPA1382 to VPA1393), and their presence was consistent among all clinical tdh- and trh-positive isolates (see Fig. 3A). Isolates categorized as tdh negative and trh positive or tdh and trh negative did not contain a PAI with homology to tdh PAI genes (Fig. 1A).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Pathogenicity islands and T3SS2. tdh PAI (A) and trh PAI (B) are the two canonical genomic islands present in clinical isolates of V. parahaemolyticus. The presence or absence of these islands is strongly correlated with the presence of the tdh and trh genes. In three of the four tdh- and trh-negative isolates, these genes are entirely absent. (C) In the fourth tdh- and trh-negative isolate, 10-4328, a trh PAI including a partial trh1 gene is mostly present. Abbreviations: neg, negative; pos, positive.

The trh PAI, composed of 81 coding sequences, was also independently assembled using reference-guided assembly and V. parahaemolyticus VIPARAQ4037 residues 1748 to 1830 as a reference and was consistently observed in isolates of both the tdh-negative and trh-positive pathotype and the tdh- and trh-positive pathotype and in all instances included trh1 and the T3SS2β genes (Fig. 1B). The association between the trh PAI and T3SS2β with the tdh-negative and trh-positive pathotype and the tdh- and trh-positive pathotype agrees with the findings of previous studies (12, 27). Three clinical tdh- and trh-negative isolates, isolates 04-2548, 09-5357, and 10-4239, did not contain the tdh PAI, the trh PAI, or a T3SS2; however, an almost complete trh PAI including a partial trh gene, urease gene cluster, and the T3SS2β genes was identified in isolate 10-4238, which was categorized as tdh and trh negative by PCR analysis by Banerjee et al. (26) (Fig. 1C). In addition to the TRH hemolysins, the T3SSβ present in this PAI also contains effectors thought to be involved in enterotoxicity and cytotoxicity (5). T3SS2β appears to be a recent acquisition by V. parahaemolyticus and is sometimes found in non-O1, non-O139, and CTX V. cholerae strains (12). The presence of the trh PAI in 10-4238 made us question whether this strain should indeed be classified as a clinical tdh- and trh-negative isolate or if it would be more accurate, on the basis of its WGS, to classify it as a tdh-negative and trh-positive strain. Therefore, this strain was removed from the remaining analysis specific to clinical tdh- and trh-negative isolates.

General genomic features of V. parahaemolyticus clinical isolates.To depict the genetic diversity of pathogenic V. parahaemolyticus strains, 38 clinical isolates (Table 1) representing each of the four previously described genotypes (tdh and trh positive, tdh negative and trh positive, tdh positive and trh negative, and tdh and trh negative) (40, 41) were extensively compared. The pangenome of clinical V. parahaemolyticus isolates was calculated and consisted of 8,399 protein-coding genes (Fig. 2A). To assess the accuracy of computing of a pangenome using draft genomes, three clinical isolates with closed genomes (see Table S1 in the supplemental material) were added to our data set and the pangenome size was recalculated. The addition of three closed genomes caused the size of the pangenome of the clinical isolates to increase to 8,609 genes. This increase would also have been expected after the addition of three draft genomes, and therefore, we concluded that our calculations based on the draft genomes were accurate.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

General genomic features of clinical V. parahaemolyticus isolates. (A) The pangenome of clinical V. parahaemolyticus isolates was constructed using the de novo assembly of 38 clinical isolates and contained 8,399 genes. Pangenomes were also assembled for each pathotype and compared. (B) A phylogenetic tree, constructed from concatenated core genes, shows the phylogeny of clinical isolates and demonstrates that each of the pathotypes is polyphyletic.

A gene intersection analysis of the accessory genomes was performed. From this analysis, we defined a core genome of orthologous genes that were shared by all clinical V. parahaemolyticus isolates. The core genome contained 3,807 protein-coding genes, which represented between 76 to 81% of each isolate's genome (Fig. 2A). The size of the core genome remained relatively stable for each additional genome added after the first four. When the closed genomes of three additional clinical isolates (see Table S1 in the supplemental material) were added to the data set and the core genome was recalculated, it decreased to a size of 3,803 protein-coding genes. This indicated that the 38 draft genomes sequenced here provided an excellent estimation of the true core genome of clinical V. parahaemolyticus isolates.

Accessory genes that were unique to each pathotype were also identified, and 654, 520, 918, and 1,097 genes were specific to the tdh- and trh-positive, tdh-negative and trh-positive, tdh-positive and trh-negative, and tdh- and trh-negative pathotypes of clinical isolates, respectively (Fig. 2A). Strain-specific unique genes were also identified, and the sizes of the unique genome varied between strains. For example, seven strains (09-3216, 10-4303, 10-7197, 09-3218, 09-4434, 09-4660, and 09-4664) possessed no unique genes, while one strain, 04-2548 (tdh and trh negative), had 405 unique genes. Clinical isolates of the tdh- and trh-negative pathotype consistently had large unique genomes.

To determine if clinical tdh- and trh-negative isolates are monophyletic and possibly the result of a single loss of a pathogenicity island, an unrooted phylogenetic tree was constructed from the concatenated amino acid sequences of the 3,899 core genes (Fig. 2B). Use of the core genes provides a high-resolution view of phylogeny. This revealed that isolates of each of the pathotypes, tdh- and trh-negative isolates, are polyphyletic. This may indicate a high degree of mobility of pathogenicity elements between V. parahaemolyticus isolates, rather than a single PAI deletion event.

The distribution of clusters of orthologous groups (COGs) of proteins was determined for the core genome, the accessory genome of each pathotype, and the unique genome of each of the tdh- and trh-negative isolates, to determine if there were differences in the proportion of the genome attributable to particular cellular processes in clinical isolates (see Table S2 in the supplemental material). Almost 20% of the core genome was classified as having an unknown function; this proportion dropped to less than 10% in each of the pathotype-specific unique genomes (Fig. 3A and B). The tdh- and trh-negative isolates had a large proportion of genes involved in cell motility relative to the proportion for the other pathotypes (Fig. 3B). Individual tdh- and trh-negative isolates also had large and functionally consistent unique genomes (Fig. 3C). Within this functional category of cell motility, each of the three clinical tdh- and trh-negative isolates contained a methyl-accepting chemotaxis protein not observed in the other pathotypes or in environmental isolates. Cell motility is generally considered to be a factor associated with the ability of V. parahaemolyticus to survive in vivo (5); therefore, the finding of an increase in the number of genes involved in motility in clinical isolates is logical.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

COG profiles of the core (A), accessory (B), and unique (C) genomes of clinical tdh- and trh-negative isolates. The numbers at the top of each column denote the number of genes in the unique genome of the corresponding strain listed at the bottom of the column. Abbreviations: neg, negative; pos, positive.

General genomic features of V. parahaemolyticus environmental isolates.To provide a basis for comparison of the clinical tdh- and trh-negative isolates, we sequenced the genomes of 5 environmental tdh- and trh-negative isolates and collected the genomes of an additional 11 from the NCBI database (Table 1). To demonstrate the wide diversity and the relationships of the strains used in this study, a rooted phylogenetic tree was constructed from the concatenated nucleotide sequences of seven housekeeping genes (recA, gyrB, dnaE, dtdS, pntA, pyrC, and tnaA) traditionally used in V. parahaemolyticus MLST analysis (Fig. 4A). Clinical and environmental isolates shared several common lineages, again demonstrating the dynamic nature of virulence factors in this species.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

General genomic features of environmental V. parahaemolyticus isolates. (A) A phylogenetic tree, constructed from concatenated MLST sequences, demonstrates the diversity of the strains used in this study, as well as the relationships between each clinical and environmental strain included in this study. (B) A rarefaction curve of the genetic diversity of clinical and environmental V. parahaemolyticus strains was created. Environmental isolates have a much greater genetic diversity than clinical isolates. In addition, these curves demonstrate that the pangenome of V. parahaemolyticus is open.

The sizes of the core genome and pangenome were calculated for the 16 environmental isolates. The core genome of the environmental tdh- and trh-negative isolates was composed of 2,773 protein-coding genes, and though the pangenome was constructed from fewer genomes, it was much larger than that of the clinical isolates at 11,669 protein-coding genes (Fig. 4B). For both the clinical and environmental isolates, the pangenome size increased, in terms of protein-coding gene number, after the addition of each genome, which indicates an open pangenome (42), and this is in agreement with the findings of earlier studies (43) (Fig. 4B). The increased size of the pangenome in environmental isolates is likely due to their ability to survive in diverse niches in situ. For example, environmental isolates were collected from colonized shrimp, muscle, oyster, and clams as well as from marine water and sediments, while clinical isolates have been preselected by use of a very narrow criterion, which is the ability to colonize humans and cause illness. This likely leads to the low level of genetic diversity observed in clinical isolates.

Whole-genome comparison of clinical and environmental tdh- and trh-negative isolates.The genomes of the 3 clinical tdh- and trh-negative isolates were compared with the genomes of 16 environmental tdh- and trh-negative isolates, resulting in the identification of 862 protein-coding genes unique to clinical isolates (Fig. 5; see also Table S3 in the supplemental material). From these genes, 529 were annotated as hypothetical proteins. A large portion of both hypothetical proteins and annotated proteins had high sequence similarity with other genes from other enteric bacteria, such as V. cholerae, Listeria monocytogenes, Campylobacter jejuni, Salmonella enterica, Escherichia coli, and Enterobacter. While this finding is not conclusive, it may indicate that the acquisition of genes from other enteric bacteria may contribute to the ability of clinical tdh- and trh-negative isolates to colonize humans during a gastrointestinal illness.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Genomic comparison of environmental and clinical tdh- and trh-negative isolates. On the basis of a comparison of protein-coding genes, 862 protein-coding genes that are unique to clinical tdh- and trh-negative isolates were identified. neg, negative.

Pathogenicity islands.By examining the complete genome of V. parahaemolyticus RIMD2210633, Hurley et al. (2006) identified seven genomic islands which occur in pathogenic V. parahaemolyticus isolates (38). VPAI-1, VPAI-4, VPAI-5, and VPAI-6 appeared to represent DNA acquired by pandemic V. parahaemolyticus isolates (38). The presence of these islands was variable in our clinical and environmental isolates (Tables 2 and 3). VPAI-1 was present in most clinical tdh-positive and trh-negative isolates but was largely absent from isolates of the other pathotypes. The entire VPAI-2 was present in most tdh-positive and trh-negative isolates and tdh-negative and trh-positive isolates but was only partially present in tdh- and trh-positive and tdh- and trh-negative strains. VPAI-6 was consistently present in each of the clinical isolates. Strains 07-2965 and 08-0278 did not typically have the same PAI profile as the other tdh-positive and trh-negative strains, and this was also reflected in the phylogenetic tree (Fig. 2B). The 10-4238 tdh- and trh-negative strain carried VPAIs more similar to those carried by the tdh-negative and trh-positive isolates, agreeing with findings presented earlier in this paper that this strain is likely a tdh-negative and trh-positive strain and was misidentified by PCR analysis. The environmental tdh- and trh-negative isolates carried pathogenicity island profiles similar to those of the clinical isolates (Table 3), this observation raises several questions about the true roles of these islands and whether their inclusion in pandemic strains is associated more with fitness in the environment than with pathogenicity.

View this table:
  • View inline
  • View popup
TABLE 2

Distribution of pathogenicity islands in clinical isolates

View this table:
  • View inline
  • View popup
TABLE 3

Distribution of pathogenicity islands in environmental isolates

IslandViewer3 was used to search the closed genome of CDC_K4557, which is also a clinical V. parahaemolyticus tdh- and trh-negative isolate, for genomic islands that are common to clinical isolates but that are not found in environmental tdh- and trh-negative isolates (37). We reasoned that if an island was present in at least some of the clinical tdh- and trh-negative isolates but absent from all of the environmental tdh- and trh-negative isolates, it would be a reasonable candidate for evaluation as a novel pathogenicity island. IslandViewer3 was used to identify 29 genomic islands on chromosome 1 and 8 genomic islands on chromosome 2 (Fig. 6). These genomic islands were assessed for their presence or absence across our 43 genomes (see Table S4 in the supplemental material). We found that islands 1, 3, 33, 35, and 36 (as denoted in Fig. 6) were present in almost every strain of V. parahaemolyticus. An island that was present in clinical tdh- and trh-negative strains but not in environmental tdh- and trh-negative strains was not identified.

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Genomic islands identified in CDC_K4557. A closed genome is required to identify novel genomic islands. To determine if a novel pathogenicity island is responsible for virulence in the tdh- and trh-negative pathotype, the closed genome of V. parahaemolyticus CDC_K4557, another clinical tdh- and trh-negative isolate, was searched for genomic islands using IslandViewer3, and the images shown above were modified from the IslandViewer3 (37) website (http://www.pathogenomics.sfu.ca/islandviewer/accession/NC_021848.1/ [top panel] and http://www.pathogenomics.sfu.ca/islandviewer/accession/NC_021822.1/ [bottom panel]). BLASTn was used to search for each of the islands in the clinical and environmental isolates from Table 1. IslandViewer3 uses multiple algorithms (IslandPick, SIGI-HMM, and IslandPath-DIMOB) to predict the presence of genomic islands. Red, islands identified by an integrative algorithm incorporating multiple mechanisms of island prediction; blue, islands predicted by the IslandPath algorithm, yellow, islands predicted by the SIGI-HMM algorithm; green, islands predicted by the IslandPick algorithm. M, millions; k, thousands.

Type VI secretion systems.Two T6SSs have previously been found in the pangenome of V. parahaemolyticus. T6SS1 (VP1387 to VP1414) is found on chromosome 1 and is commonly associated with clinical isolates, while T6SS2 (VPA1025 to VPA1046) is found on chromosome 2 and has been found in all tested strains (16, 20). In the current study, six clinical isolates (10-4238, 08-0278, 09-5357, 07-2965, 09-3217, and 10-4255) and eight environmental isolates (ISF-54-12, S171, VPCR-2009, VPTS-2009, SG176, 22720, J-C2-34, and VH3) did not have a T6SS1. Of the isolates that had a T6SS1, variance in this gene cluster was observed between isolates, although variance occurred only in the VP1388, VP1389, and VP1390 genes (Fig. 7). There were five different alleles of these genes, and a phylogenetic tree is shown to demonstrate which isolates contained which alleles (Fig. 7). The 10-4239 isolate had unique alleles for each of these genes which were not observed in any of the other isolates. The N terminus of VP1388 was conserved in all isolates that had a T6SS1, while variance was observed in the C terminus. There is biological significance underlying the variation observed in VP1388 to VP1390. VP1388 has previously been identified to be an antibacterial effector, and VP1389 is its associated immunity protein (21). Changes in the effector must be accompanied by changes in the immunity protein to maintain self-protection. The putative functionality of the third larger gene downstream (VP1390) is still unknown, but on the basis of its association with VP1388 and VP1390, it may have a role in antimicrobial activity. The finding of a novel putative effector in a V. parahaemolyticus clinical tdh- and trh-negative isolate indicates that this protein should be further investigated for roles during infection.

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

T6SS1. The gene cluster homologous to T6SS1 is present in 32 of our clinical isolates; however, variance in the VP1388 to VP1390 genes was observed between isolates.

Conclusion.The ability of tdh- and trh-negative strains to cause clinical illness is still controversial, and several theories have been proposed to explain why tdh- and trh-negative strains are sometimes isolated from clinical cases, including coinfection with pathogenic V. parahaemolyticus strains, the loss of virulence genes during infection, the presence of novel and uncharacterized virulence factors, or the fact that they play a role in a multistrain infection. However, in this investigation we have identified 862 genes that are present in clinical tdh- and trh-negative isolates but that are not present in environmental isolates. Several of these genes are highly homologous to genes from other enteric bacteria, indicating that horizontal gene transfer may play an important role in the ability of tdh- and trh-negative isolates to survive in the human gastrointestinal tract. In addition, tdh- and trh-negative isolate 10-4239 contains a unique T6SS1 effector/immunity gene combination that should be investigated further.

ACKNOWLEDGMENTS

We are very grateful to Robyn Kenwell for technical assistance with the sequencing of the genomes analyzed in this paper. We also thank John Austin and Franco Pagotto of the BMH Research Division of Health Canada for reviewing the manuscript and offering helpful comments.

This work was funded (A-base) by Health Canada in support of Canada's Food Safety Programs. J.R. is supported by the Visiting Fellow in a Government Laboratory program.

J.R., N.P., C.C.L., and S.K.B. conceived of and designed the experiments. J.R. and C.C.L. performed all wet laboratory experimental, molecular, and next-generation sequencing work. J.R. and N.P. conducted the in silico analysis, including the use of novel tools contributed by A.W.P.; J.R. wrote the paper with input from all other authors.

FOOTNOTES

    • Received 22 October 2015.
    • Accepted 25 November 2015.
    • Accepted manuscript posted online 4 December 2015.
  • Supplemental material for this article may be found at http://dx.doi.org/10.1128/AEM.03465-15.

  • Copyright © 2016, American Society for Microbiology. All Rights Reserved.

REFERENCES

  1. 1.↵
    1. Ceccarelli D,
    2. Hasan NA,
    3. Huq A,
    4. Colwell RR
    . 2013. Distribution and dynamics of epidemic and pandemic Vibrio parahaemolyticus virulence factors. Front Cell Infect Microbiol 3:97. doi:10.3389/fcimb.2013.00097.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Joseph SW,
    2. Colwell RR,
    3. Kaper JB
    . 1982. Vibrio parahaemolyticus and related halophilic vibrios. Crit Rev Microbiol 10:77–124. doi:10.3109/10408418209113506.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Nishibuchi M,
    2. Fasano A,
    3. Russell RG,
    4. Kaper JB
    . 1992. Enterotoxigenicity of Vibrio parahaemolyticus with and without genes encoding thermostable direct hemolysin. Infect Immun 60:3539–3545.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Honda T,
    2. Ni Y,
    3. Miwatani T,
    4. Adachi T,
    5. Kim J
    . 1992. The thermostable direct hemolysin of Vibrio parahaemolyticus is a pore-forming toxin. Can J Microbiol 38:1175–1180. doi:10.1139/m92-192.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Broberg CA,
    2. Calder TJ,
    3. Orth K
    . 2011. Vibrio parahaemolyticus cell biology and pathogenicity determinants. Microbes Infect 13:992–1001. doi:10.1016/j.micinf.2011.06.013.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Hiyoshi H,
    2. Kodama T,
    3. Iida T,
    4. Honda T
    . 2010. Contribution of Vibrio parahaemolyticus virulence factors to cytotoxicity, enterotoxicity, and lethality in mice. Infect Immun 78:1772–1780. doi:10.1128/IAI.01051-09.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Matsuda S,
    2. Kodama T,
    3. Okada N,
    4. Okayama K,
    5. Honda T,
    6. Iida T
    . 2010. Association of Vibrio parahaemolyticus thermostable direct hemolysin with lipid rafts is essential for cytotoxicity but not hemolytic activity. Infect Immun 78:603–610. doi:10.1128/IAI.00946-09.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Letchumanan V,
    2. Chan K-G,
    3. Lee L-H
    . 2014. Vibrio parahaemolyticus: a review on the pathogenesis, prevalence, and advance molecular identification techniques. Front Microbiol 5:705. doi:10.3389/fmicb.2014.00705.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Miyamoto Y,
    2. Kato T,
    3. Obara Y,
    4. Akiyama S,
    5. Takizawa K,
    6. Yamai S
    . 1969. In vitro hemolytic characteristic of Vibrio parahaemolyticus: its close correlation with human pathogenicity. J Bacteriol 100:1147–1149.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Makino K,
    2. Oshima K,
    3. Kurokawa K,
    4. Yokoyama K,
    5. Uda T,
    6. Tagomori K,
    7. Iijima Y,
    8. Najima M,
    9. Nakano M,
    10. Yamashita A,
    11. Kubota Y,
    12. Kimura S,
    13. Yasunaga T,
    14. Honda T,
    15. Shinagawa H,
    16. Hattori M,
    17. Iida T
    . 2003. Genome sequence of Vibrio parahaemolyticus: a pathogenic mechanism distinct from that of V. cholerae. Lancet 361:743–749. doi:10.1016/S0140-6736(03)12659-1.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Izoré T,
    2. Perdu C,
    3. Job V,
    4. Attree I,
    5. Faudry E,
    6. Dessen A
    . 2011. Structural characterization and membrane localization of ExsB from the type III secretion system (T3SS) of Pseudomonas aeruginosa. J Mol Biol 413:236–246. doi:10.1016/j.jmb.2011.07.043.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Okada N,
    2. Iida T,
    3. Park K-S,
    4. Goto N,
    5. Yasunaga T,
    6. Hiyoshi H,
    7. Matsuda S,
    8. Kodama T,
    9. Honda T
    . 2009. Identification and characterization of a novel type III secretion system in trh-positive Vibrio parahaemolyticus strain TH3996 reveal genetic lineage and diversity of pathogenic machinery beyond the species level. Infect Immun 77:904–913. doi:10.1128/IAI.01184-08.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Park K-S,
    2. Ono T,
    3. Rokuda M,
    4. Jang M-H,
    5. Okada K,
    6. Iida T,
    7. Honda T
    . 2004. Functional characterization of two type III secretion systems of Vibrio parahaemolyticus. Infect Immun 72:6659–6665. doi:10.1128/IAI.72.11.6659-6665.2004.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Noriea NF, III,
    2. Johnson CN,
    3. Griffitt KJ,
    4. Grimes DJ
    . 2010. Distribution of type III secretion systems in Vibrio parahaemolyticus from the northern Gulf of Mexico. J Appl Microbiol 109:953–962. doi:10.1111/j.1365-2672.2010.04722.x.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Ham H,
    2. Orth K
    . 2012. The role of type III secretion system 2 in Vibrio parahaemolyticus pathogenicity. J Microbiol 50:719–725. doi:10.1007/s12275-012-2550-2.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Salomon D,
    2. Gonzalez H,
    3. Updegraff BL,
    4. Orth K
    . 2013. Vibrio parahaemolyticus type VI secretion system 1 is activated in marine conditions to target bacteria, and is differentially regulated from system 2. PLoS One 8:e61086. doi:10.1371/journal.pone.0061086.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Salomon D,
    2. Orth K
    . 2015. Type VI secretion system. Curr Biol 25:R265–R266. doi:10.1016/j.cub.2015.02.031.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Pukatzki S,
    2. McAuley SB,
    3. Miyata ST
    . 2009. The type VI secretion system: translocation of effectors and effector-domains. Curr Opin Microbiol 12:11–17. doi:10.1016/j.mib.2008.11.010.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    1. Boyer F,
    2. Fichant G,
    3. Berthod J,
    4. Vandenbrouck Y,
    5. Attree I
    . 2009. Dissecting the bacterial type VI secretion system by a genome wide in silico analysis: what can be learned from available microbial genomic resources? BMC Genomics 10:104. doi:10.1186/1471-2164-10-104.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Yu Y,
    2. Yang H,
    3. Li J,
    4. Zhang P,
    5. Wu B,
    6. Zhu B,
    7. Zhang Y,
    8. Fang W
    . 2012. Putative type VI secretion systems of Vibrio parahaemolyticus contribute to adhesion to cultured cell monolayers. Arch Microbiol 194:827–835. doi:10.1007/s00203-012-0816-z.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Salomon D,
    2. Kinch LN,
    3. Trudgian DC,
    4. Guo X,
    5. Klimko JA,
    6. Grishin NV,
    7. Mirzaei H,
    8. Orth K
    . 2014. Marker for type VI secretion system effectors. Proc Natl Acad Sci U S A 111:9271–9276. doi:10.1073/pnas.1406110111.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Lüdeke CHM,
    2. Kong N,
    3. Weimer BC,
    4. Fischer M,
    5. Jones JL
    . 2015. Complete genome sequences of a clinical isolate and an environmental isolate of Vibrio parahaemolyticus. Genome Announc 3(2):e00216-15. doi:10.1128/genomeA.00216-15.
    OpenUrlCrossRef
  23. 23.↵
    1. Meador CE,
    2. Parsons MM,
    3. Bopp CA,
    4. Gerner-Smidt P,
    5. Painter JA,
    6. Vora GJ
    . 2007. Virulence gene- and pandemic group-specific marker profiling of clinical Vibrio parahaemolyticus isolates. J Clin Microbiol 45:1133–1139. doi:10.1128/JCM.00042-07.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Bhoopong P,
    2. Palittapongarnpim P,
    3. Pomwised R,
    4. Kiatkittipong A,
    5. Kamruzzaman M,
    6. Nakaguchi Y,
    7. Nishibuchi M,
    8. Ishibashi M,
    9. Vuddhakul V
    . 2007. Variability of properties of Vibrio parahaemolyticus strains isolated from individual patients. J Clin Microbiol 45:1544–1550. doi:10.1128/JCM.02371-06.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Jones JL,
    2. Lüdeke CHM,
    3. Bowers JC,
    4. Garrett N,
    5. Fischer M,
    6. Parsons MB,
    7. Bopp CA,
    8. DePaola A
    . 2012. Biochemical, serological, and virulence characterization of clinical and oyster Vibrio parahaemolyticus isolates. J Clin Microbiol 50:2343–2352. doi:10.1128/JCM.00196-12.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Banerjee SK,
    2. Kearney AK,
    3. Nadon CA,
    4. Peterson C-L,
    5. Tyler K,
    6. Bakouche L,
    7. Clark CG,
    8. Hoang L,
    9. Gilmour MW,
    10. Farber JM
    . 2014. Phenotypic and genotypic characterization of Canadian clinical isolates of Vibrio parahaemolyticus collected from 2000 to 2009. J Clin Microbiol 52:1081–1088. doi:10.1128/JCM.03047-13.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Hazen TH,
    2. Lafon PC,
    3. Garrett NM,
    4. Lowe TM,
    5. Silberger DJ,
    6. Rowe LA,
    7. Frace M,
    8. Parsons MB,
    9. Bopp CA,
    10. Rasko DA,
    11. Sobecky PA
    . 2015. Insights into the environmental reservoir of pathogenic Vibrio parahaemolyticus using comparative genomics. Front Microbiol 6:204. doi:10.3389/fmicb.2015.00204.
    OpenUrlCrossRefPubMed
  28. 28.
    1. Bankevich A,
    2. Nurk S,
    3. Antipov D,
    4. Gurevich AA,
    5. Dvorkin M,
    6. Kulikov AS,
    7. Lesin VM,
    8. Nikolenko SI,
    9. Pham S,
    10. Prjibelski AD,
    11. Pyshkin AV,
    12. Sirotkin AV,
    13. Vyahhi N,
    14. Tesler G,
    15. Alekseyev MA,
    16. Pevzner PA
    . 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455–477. doi:10.1089/cmb.2012.0021.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Li H,
    2. Durbin R
    . 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25:1754–1760. doi:10.1093/bioinformatics/btp324.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    1. Milne I,
    2. Stephen G,
    3. Bayer M,
    4. Cock PJA,
    5. Pritchard L,
    6. Cardle L,
    7. Shaw PD,
    8. Marshall D
    . 2013. Using Tablet for visual exploration of second-generation sequencing data. Brief Bioinform 14:193–202. doi:10.1093/bib/bbs012.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Marchesi JR
    . 2012. Metagenomics: current innovations and future trends. Future Microbiol 7:813–814. doi:10.2217/fmb.12.41.
    OpenUrlCrossRef
  32. 32.↵
    1. Chen JC,
    2. Sun S,
    3. Li W,
    4. Wooley JC
    . 2011. A community cyberinfrastructure for metagenomics: CAMERA 2.0. In Metagenomics: current innovations and future trends. Caister Academic Press, Norfolk, United Kingdom.
  33. 33.↵
    1. Seemann T
    . 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069. doi:10.1093/bioinformatics/btu153.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    1. Edgar RC
    . 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32:1792–1797. doi:10.1093/nar/gkh340.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Stamatakis A
    . 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30:1312–1313. doi:10.1093/bioinformatics/btu033.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Sullivan MJ,
    2. Petty NK,
    3. Beatson SA
    . 2011. EasyFig: a genome comparison visualizer. Bioinformatics 27:1009–1010. doi:10.1093/bioinformatics/btr039.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Dhillon BK,
    2. Laird MR,
    3. Shay JA,
    4. Winsor GL,
    5. Lo R,
    6. Nizam F,
    7. Pereira SK,
    8. Waglechner N,
    9. McArthur AG,
    10. Langille MGI,
    11. Brinkman FSL
    . 2015. IslandViewer 3: more flexible, interactive genomic island discovery, visualization and analysis. Nucleic Acids Res 43(W1):W104–W108. doi:10.1093/nar/gkv401.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Hurley CC,
    2. Quirke A,
    3. Reen FJ,
    4. Boyd EF
    . 2006. Four genomic islands that mark post-1995 pandemic Vibrio parahaemolyticus isolates. BMC Genomics 7:104. doi:10.1186/1471-2164-7-104.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Boyd EF,
    2. Cohen ALV,
    3. Naughton LM,
    4. Ussery DW,
    5. Binnewies TT,
    6. Stine OC,
    7. Parent MA
    . 2008. Molecular analysis of the emergence of pandemic Vibrio parahaemolyticus. BMC Microbiol 8:110. doi:10.1186/1471-2180-8-110.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Nishibuchi M,
    2. Kaper JB
    . 1995. Thermostable direct hemolysin gene of Vibrio parahaemolyticus: a virulence gene acquired by a marine bacterium. Infect Immun 63:2093–2099.
    OpenUrlFREE Full Text
  41. 41.↵
    1. Honda T,
    2. Ni Y,
    3. Miwatani T
    . 1988. 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 56:961–965.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Tettelin H,
    2. Riley D,
    3. Cattuto C,
    4. Medini D
    . 2008. Comparative genomics: the bacterial pan-genome. Curr Opin Microbiol 11:472–477. doi:10.1016/j.mib.2008.09.006.
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.↵
    1. Li L,
    2. Wong H,
    3. Nong W,
    4. Cheung MK,
    5. Law PTW,
    6. Kam KM,
    7. Kwan HS
    . 2014. Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus. BMC Genomics 15:1135. doi:10.1186/1471-2164-15-1135.
    OpenUrlCrossRefPubMed
  44. 44.
    1. Banerjee S,
    2. Petronella N,
    3. Chew Leung C,
    4. Farber J
    . 2015. Draft genome sequences of four Vibrio parahaemolyticus isolates from clinical cases in Canada. Genome Announc 3:e01482-14. doi:10.1128/genomeA.01482-14.
    OpenUrlAbstract/FREE Full Text
  45. 45.
    1. Liu M,
    2. Chen S
    . 2013. Draft genome sequence of Vibrio parahaemolyticus V110, isolated from shrimp in Hong Kong. Genome Announc 1:e00300-13. doi:10.1128/genomeA.00300-13.
    OpenUrlAbstract/FREE Full Text
  46. 46.
    1. Castillo D,
    2. Jun JW,
    3. D'Alvise P,
    4. Middelboe M,
    5. Gram L,
    6. Liu S,
    7. Katharios P
    . 2015. Draft genome sequence of Vibrio parahaemolyticus VH3, isolated from an aquaculture environment in Greece. Genome Announc 3:e00731-15. doi:10.1128/genomeA.00731-15.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Ronholm J,
    2. Petronella N,
    3. Kenwell R,
    4. Banerjee S
    . 2015. Draft whole-genome sequences of 14 Vibrio parahaemolyticus clinical isolates with an ambiguous K serogroup. Genome Announc 3:e00217-15. doi:10.1128/genomeA.00217-15.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Download PDF
Citation Tools
Genomic Features of Environmental and Clinical Vibrio parahaemolyticus Isolates Lacking Recognized Virulence Factors Are Dissimilar
J. Ronholm, N. Petronella, C. Chew Leung, A. W. Pightling, S. K. Banerjee
Applied and Environmental Microbiology Feb 2016, 82 (4) 1102-1113; DOI: 10.1128/AEM.03465-15

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Applied and Environmental Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Genomic Features of Environmental and Clinical Vibrio parahaemolyticus Isolates Lacking Recognized Virulence Factors Are Dissimilar
(Your Name) has forwarded a page to you from Applied and Environmental Microbiology
(Your Name) thought you would be interested in this article in Applied and Environmental Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Genomic Features of Environmental and Clinical Vibrio parahaemolyticus Isolates Lacking Recognized Virulence Factors Are Dissimilar
J. Ronholm, N. Petronella, C. Chew Leung, A. W. Pightling, S. K. Banerjee
Applied and Environmental Microbiology Feb 2016, 82 (4) 1102-1113; DOI: 10.1128/AEM.03465-15
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • MATERIALS AND METHODS
    • RESULTS AND DISCUSSION
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

Cited By...

About

  • About AEM
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AppEnvMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

 

Print ISSN: 0099-2240; Online ISSN: 1098-5336