Abstract
Vibrio parahaemolyticus is a leading cause of foodborne infections in China and a threat to human health worldwide. The main objective of this study is to determine the prevalence and characteristic of V. parahaemolyticus isolates in fish, oyster and shrimp samples from the South China domestic consumer market. To accomplish this, we examined 504 seafood samples from 11 provinces of China. The prevalence rates were 9.38, 30.36, and 25.60%, respectively. In summer (33.33%), the prevalence of V. parahaemolyticus was more common than that detected in the winter (14.01%). In addition, we identified 98 V. parahaemolyticus strains. The antimicrobial resistance trends of our seafood isolates to 15 antimicrobial agents revealed that major isolates were resistant to ampicillin (79.59%). Furthermore, 68.38% of the isolates were identified as being multidrug resistance. The prevalence of tdh or trh genes among the isolates was 8.16 and 12.24%, respectively. ERIC-PCR and multilocus sequence typing (MLST) results enabled classification of the isolates (n = 98) into different clusters, revealing genetic variation and relatedness among the isolates. Thus, our findings demonstrate the prevalence of V. parahaemolyticus in a variety of common seafood consumed domestically in China and provides insights into the dissemination of antibiotic-resistant strains, which should improve our microbiological risk assessment knowledge associated with V. parahaemolyticus in seafoods.
Introduction
Foodborne disease remains a threat to public health worldwide. Vibrio parahaemolyticus is a halophilic Gram-negative bacterium that is commonly associated with food-borne gastroenteritis. This bacterium has been found in seawater, marine fish, shellfish, etc. (; ; ) and causes outbreaks worldwide. In China, V. parahaemolyticus is also an important cause of food poisoning associated with consumption of seafood (; ). It was also reported that the population levels of V. parahaemolyticus tends to show strong seasonal trends (). Consequently, foodborne outbreaks caused by V. parahaemolyticus typically show a seasonal difference, peaking in the warmer months of the United States (). However, to date, the presence and contamination levels of V. parahaemolyticus in seafood from different seasons has received less attention, and little information is available among various types of samples.
As a result of the excessive use of antibiotics in human and aquaculture systems, some antibiotics have no longer been effective in controlling pathogens infections during the past several decades (). Variation in antibiotic resistance patterns among V. parahaemolyticus strains isolated from different countries and sources has been observed (). Recent studies regarding a V. parahaemolyticus isolate demonstrated resistance to certain common antibiotics such as ampicillin and kanamycin (; ). The frequent occurrence of multi-drug resistance in V. parahaemolyticus is among the most important of public health concerns (). In order to improve yield and treatment of infectious diseases in farm products, increase levels of antibiotics in animal feed and water is a common behavior (). Therefore, it is very important to set up a monitoring system to research on antimicrobial-resistance trends.
An important virulence genes, toxR is involved in the regulation of many genes in V. parahaemolyticus (). The pathogenicity of V. parahaemolyticus is strongly correlated with the production of either thermostable direct hemolysin (TDH), TDH related hemolysin (TRH), or both (; ). V. parahaemolyticus strains isolated from diarrheal patients produce thermostable TDH which is a toxin with several biological properties, including hemolytic activity, enterotoxicity and cytotoxicity, TDH is also present in most of the Kanagawa phenomenon (KP)-positive clinical strains. TRH is believed to act similarly to TDH. (; ; ). Presently, PCR-based methods is useful in identified of tdh and trh gene in V. parahaemolyticus isolates ().
The utility of molecular techniques has already been established in epidemiological studies of V. parahaemolyticus infections, such as pulsed-field gel electrophoresis (PFGE) () and random amplified polymorphic DNA (RAPD) analysis (). Enterobacterial repetitive intergenic consensus sequence PCR (ERIC-PCR) is an easy, fast, and relatively cheap method that is comparable to PFGE in terms of its discriminatory power and reproducibility, It is also superior to other PCR-based fingerprinting techniques (; ). Multilocus sequence typing (MLST) is based on sequence analysis of housekeeping (HK) genes and is a useful method to trace the global epidemiology of bacterial pathogens. MLST, which is sequence based, provides a clear determination of the genetic features of strains that is consistent from one laboratory to another. Previous MLST studies of V. parahaemolyticus isolates have provided a better understanding of the genetic relatedness within these bacteria. These studies also revealed the relative evolutionary importance of mutations and lateral gene transfer events among these strains (; ).
Overall, V. parahaemolyticus strains are often isolated from fish, shrimp and other types of seafoods. Recently, V. parahaemolyticus was also found in ready to eat food (; ). In the Zhejiang province of China alone, there were 71 outbreaks caused by V. parahaemolyticus resulting in 933 illnesses and 117 hospitalizations from 2010 to 2014. In our study, we analyzed 224 fish, 112 oyster, and 168 shrimp samples. In China, seafood is very popular and higher consumption is correlated with an increase in the overall standard of living in the population. Our study aimed to find the differences in the seasonal prevalence of V. parahaemolyticus in all provinces of South China. We next characterized the prevalence all of the identified isolates to determine their genetic relatedness by phenotyping and genotyping methods. We attempted to confirm the virulence and antibiotic resistance trends. We also performed ERIC and MLST typing of the isolated strains, which revealed the molecular diversity of the V. parahaemolyticus isolates. These findings serve as baseline values that support the establishment of a national seafood surveillance system to ensure food safety in China.
Materials and Methods
Sample Collection
From July 2015 to July 2017, a total of 504 aquatic products samples, including 224 fish samples, 112 oyster samples and 168 shrimp samples were collected from retail markets from 14 different cities in China, belonging to 11 provinces (Figure 1). In these regions, the climate is cold from September to March (winter), and it is hot from March to September (summer). From each city, we collected eight fish samples, four oyster samples, and six shrimp samples (from four different retail markets) during both summer and winter months. We put the samples in sterile sealed plastic bags and sent them to the lab in an ice box (4°C). Samples were analyzed immediately thereafter.
FIGURE 1
Detection and Enumeration of V. parahaemolyticus
The prevalence and bacterial load of V. parahaemolyticus in the samples were tested using the most probable number (MPN) method according to the National Food Safety Standards of China (GB4789.7-2013); this method was reported previously by . In brief, we first weighed and homogenized 25 g of each sample. Then, 225 mL of alkaline peptone water (APW) with 3% NaCl was added; 3 × 1 mL of each sample dilution was inoculated into 9 mL of APW (3% NaCl) when prepared up to a 1:1000 dilution. Samples were incubated at 37°C for 18 h. After that, we used an inoculation loop to streaked the solution onto thiosulfate-citrate-bile salts-sucrose (TCBS) agar plates and incubated them at 37°C for 24 h. The r green or blue green colonies of 2–3 mm in diameter were presumed to be V. parahaemolyticus isolates, three to five colonies (if present) were chose from each plate, streaked onto Chromogenic Vibrio Medium and culture at 37°C for 24 h. One (if present) colony (mauve) from each Chromogenic Vibrio Medium plate was selected for further identification tests, such as Gram staining, NaCl tolerance tests, oxidase activity and analysis with API 20E diagnostic strips (BioMerieux, Marcy-l’Étoile, France). The pollution result in the samples were determined using an MPN table. The levels of contamination in the samples were determined using an MPN table.
Antimicrobial Susceptibility Testing
The susceptibility of the V. parahaemolyticus isolates to antibiotics was examined by the disk-diffusion method, according to the guidelines of the Clinical and Laboratory Standards Institute (; ). Muller – Hinton agar and a panel of 15 antibiotics disks were selected for the resistance analysis. The 15 common antimicrobials used in this study belonged to six classes and were as follows: β-lactam [ampicillin: AMP (10 μg); piperacillin: PRL (20 μg); cefotaxime: CTX (30 μg); cefoxitin: FOX (30 μg); cephazolin: KZ (30 μg); imipenem: IPM (20 μg); meropenem: MEM (20 μg)], aminoglycoside [gentamicin: CN (10 μg); kanamycin: K (30 μg); streptomycin: S (10 μg)], tetracycline [tetracycline: TET (30 μg)], quinolone [ciprofloxacin:CIP (5 μg); levofloxacin: LEV (20 μg)], sulfonamides [trimethoprim-sulfamethoxazole: SXT (25 μg)], chloramphenicol [chloramphenicol:C (30 μg)]. Following the methods of the CLSI, sensitive (S), intermediate (I), or resistant (R) classification system was used.
Detection of Virulence Genes: toxR, tdh, and trh
According to the manufacturer’s instructions, genomic DNA was extracted from each V. parahaemolyticus isolate using a bacterial DNA extraction kit (Sangon, Shanghai, China).
The toxR gene is important and appears to be well conserved among V. parahaemolyticus isolates. The primers were as follows F: GTCTTCTGACGCAATCGTTG, R: ATACGAGTGGTTGCTGTCATG (). The tdh and trh genes detection was executed as previously reported (). The following primers were used: Tdh-F:CTGTCCCTTTTCCTGCCCCCG, Tdh-R:AGCCAGACACCGCTGCCATTG;Trh-F:ACCTTTTCCTTCTCCWGGKTCSG,Trh-F:CCGCTCTCATATGCYTCGACAKT). All the oligonucleotide primers were synthesized by Sangon Biotech (Shanghai, China). The PCR reaction system contained DNA template (1 μL), 0.5 μM of each primer, 12.5 μL of :2 × PCR mix (Qiagen), and ddH2O (9.5 μL). The amplified thermal-cycling program was set with the following conditions: denaturation at 95°C (5 min); 40 cycles of 95°C for 1 min, 62°C for 1 min, and 72°C for 1 min; and a final extension step of 72°C for 5 min. PCR amplicons were electrophoresed on 2.0% (wt/vol) agarose gels containing GoldView. The images were captured digitally and analyzed using a gel imaging system. V. parahaemolyticus strains ATCC33847 (tdh+) and ATCC17802 (trh+) were used as positive controls, and distilled water was used as a negative control.
ERIC-PCR Analysis
ERIC-PCR analysis was performed on the V. parahaemolyticus isolates using a previously reported method (; ). A pair of primers, ERIC - F (5-ATGTAAGCTCCTGGGGATTCAC-3) and ERIC -R (5-AAGTAAGTGACTGGGGTGAGCG-3) were used. ERIC-PCR typing was performed as follows: 1 μl of genomic DNA (100 ng) was added to a master mixture that contained 0.6 μmol/L of each primer, 12.5 μL of 2 × Long Taq mix (Dongsheng Biotech, Guangzhou, China), and ddH2O to a total volume of 25 μL per reaction. After an initial denaturation at 95°C for 5 min, 35 cycles of amplification were performed under the following conditions: denaturation at 94°C for 45 s, annealing 52°C for 1 min and extension at 72°C for 3 min followed by a final extension step at 72°C for 10 min. The amplicons were electrophoresed on 2.0% (wt/vol) agarose gels containing GoldView. The images were captured digitally in TIFF file format for further analysis.
MLST Analysis
The V. parahaemolyticus MLST website and database1 shows the MLST analysis procedure (). Seven genes were selected (recA, gyrB, dnaE, dtdS, pntA, pyrC, and tnaA). The amplified thermal-cycling program was designed with the following conditions: denaturation at 95°C (5 min); 30 cycles of 95°C for 1 min, 58°C for 1 min, and 72°C for 1 min; and a final extension step of 72°C for 10 min. A BGI instrument (Shenzhen, China) was used to sequence the PCR amplicons. The assigning allele numbers and defined sequence types (STs)1. BioEdit was used to determine the alignments of these sequences.
Statistical Analysis
The size of each band in the ERIC patterns was determined and the data were coded as 0 (absent) or 1 (present). Cluster analysis was performed with NTSYS-pc (version 2.10), a numerical taxonomy and multivariate analysis software package (), based on Dice’s similarity coefficient (SD), with a 1% position tolerance and the unweighted-pair group method with arithmetic averages (UPGMA). The evolution tree of the concatenated sequences of the seven loci (MLST) was built based on the method involving the Kimura-2-parameter in Mega 6.0 ().
Results
Prevalence and Level of V. parahaemolyticus
Of the 504 seafood samples that we tested, 98 (19.44%) were identified as V. parahaemolyticus positive. This included 21 (9.38%) of the 224 fish samples, 34 (30.36%) of the 112 oyster samples and 43 (25.60%) of the 168 shrimp samples. These results are shown in Table 1. Overall, the degree of V. parahaemolyticus contamination in these samples varied from 1.50 to 1000 MPN/g. Most of the positive samples had a level of 3 to 10 MPN/g (62/98, 63.27%). Moreover, 19.39% (19/98) of the positive samples exceeded 100 MPN/g, and 17 samples were below 3 MPN/g. In total, 98 V. parahaemolyticus isolates were confirmed (Supplementary Table S1).
Table 1
| Samples | Number of samples analyzed | Number of samples positive (%) | Number of samples containing the pathogen | ||
|---|---|---|---|---|---|
| <3 (MPN/g) | 3–100 (MPN/g) | >100–103 (MPN/g) | |||
| Fish | 224 | 21 (9.38) | 3 | 13 | 5 |
| Oyster | 112 | 34 (30.36) | 8 | 16 | 10 |
| Shrimp | 168 | 43 (25.60) | 6 | 33 | 4 |
| Total | 504 | 98 (19.44) | 17 | 62 | 19 |
Prevalence of Vibrio parahaemolyticus in food samples from South China.
Considering the influence of the season, the number of positive summer samples was much higher than the number of positive winter samples (Table 2). The isolation rate of V. parahaemolyticus in summer reached 33.33%, whereas it was 14.01% in the winter. Moreover, the V. parahaemolyticus population differed significantly for samples collected during the winter and summer. Densities among samples collected in the summer varied less, and the contamination level was higher than that observed in the winter samples. In summer, the mean levels of V. parahaemolyticus in samples were 99.08 MPN/g, which is higher than that observed in the winter (22.13 MPN/g). A total of 14 samples with densities above 100 MPN/g were detected.
Table 2
| Samples | Number of samples analyzed | Number of samples positive (%) | Number of samples containing the pathogen | ||
|---|---|---|---|---|---|
| <3 (MPN/g) | 3–100 (MPN/g) | >100–103 (MPN/g) | |||
| Summer | 207 | 69 (33.33) | 9 | 46 | 14 |
| Winter | 207 | 29 (14.01) | 8 | 16 | 5 |
| Total | 905 | 98 (19.44) | 17 | 62 | 19 |
Prevalence of Vibrio parahaemolyticus during different seasons.
Antimicrobial Susceptibility of the V. parahaemolyticus Isolates
Isolates of V. parahaemolyticus were tested for antibiotic susceptibility, the resistance patterns to 15 antibiotics are shown in Table 3. The isolates were most resistant to ampicillin, with resistance and intermediate rates of 79.59 and 12.24%, respectively. In addition, the isolates exhibited relatively high resistance rates, of 68.37, 39.80, and 39.80%, for streptomycin, cefazolin, and kanamycin, respectively. Fortunately, all the examined isolates were susceptible to cefotaxime, cefoxitin, imipenem and meropenem. Among the remaining tested antibiotics, the isolates were partially susceptible to piperacillin, ciprofloxacin, levofloxacin, or chloramphenicol. However, among all the isolates, there were multidrug-resistant isolates (Vps15) showing resistance to seven antibiotics and some strains showing resistance to five antibiotics. In addition, 68.38% (67/98) of the isolates were resistance to more than three antibiotics (Supplementary Table S1).
Table 3
| Antimicrobial agents | Vibrio parahaemolyticus (n = 98) | ||
|---|---|---|---|
| Number | Number | Number | |
| (%) of R | (%) of I | (%) of S | |
| Ampicillin (AMP) | 78 (79.59) | 12 (12.24) | 8 (8.16) |
| Piperacillin (PRL) | 8 (8.16) | 1 (1.02) | 89 (90.82) |
| Cefotaxime (CTX) | 0 (0.00) | 2 (2.04) | 96 (97.96) |
| Cefoxitin (FOX) | 0 (0.00) | 1 (1.02) | 97 (98.98) |
| Cefazolin (KZ) | 39 (39.80) | 18 (18.37) | 41 (41.84) |
| Imipenem (IPM) | 0 (0.00) | 2 (2.04) | 96 (97.96) |
| Meropenem (MEM) | 0 (0.00) | 0 (0.00) | 98 (100.00) |
| Gentamicin (CN) | 17 (17.35) | 20 (20.41) | 61 (62.24) |
| Kanamycin (K) | 39 (39.80) | 24 (24.49) | 35 (35.71) |
| Streptomycin (S) | 67 (68.37) | 26 (26.53) | 5 (5.10) |
| Tetracycline (TE) | 13 (13.27) | 8 (8.16) | 77 (78.57) |
| Ciprofloxacin (CIP) | 5 (5.10) | 2 (2.04) | 91 (92.86) |
| Levofloxacin (LEV) | 1 (1.02) | 5 (5.10) | 92 (93.88) |
| Trimethoprim-sulfamethoxazole (SXT) | 9 (9.18) | 12 (12.24) | 77 (78.57) |
| Chloramphenicol (C) | 2 (2.04) | 1 (1.02) | 95 (96.94) |
Antimicrobial resistance profiles of Vibrio parahaemolyticus isolates.
The toxR, tdh, and trh Genes Test
To detect pathogenic isolates, all of the 98 V. parahaemolyticus isolates were examined for the presence of the toxR, trh, and tdh genes. This analysis is summarized in Table 4. All of the isolates were positive for tox R. Among them, 8.16% (8/98) and 12.24% (12/98) of the V. parahaemolyticus strains carried the tdh, or trh gene, respectively; in contrast, none of the isolates harbored both the tdh and trh genes. For the tdh -positive strains, two were found in a fish source, two were identified from oysters, and three were found from shrimp. Of the trh-positive isolates, five samples were from fish, four samples were from oyster, and three samples were from shrimp.
Table 4
| Sample category | Number of isolates | tox R-positive (%) | tdh- positive (%) | trh- positive (%) |
|---|---|---|---|---|
| Fish | 21 | 100 (21/21) | 2 | 5 |
| Oyster | 34 | 100 (34/34) | 2 | 4 |
| Shrimp | 43 | 100 (43/43) | 4 | 3 |
| Total | 98 | 100 (98/98) | 8.16 (8/98) | 12.24 (12/98) |
Distribution of the tox R, tdh, and trh genes in V. parahaemolyticus isolates.
ERIC-PCR and MLST
Figure 2 shows the results from ERIC-PCR analysis of the 98 isolates. There were 5 to10 amplification bands, with sizes ranging between 100 bp to about 6000 bp. Analysis of the ERIC-PCR patterns revealed that those strains could be divided into six clusters, A, B, C, D, E, and F at a relative similarity coefficient of 0.67. The isolates proved to be very diverse genetically. Most isolates were distributed among the A and D clusters. However, we did not find a relationship between the different sample types and the ERIC-PCR clusters
FIGURE 2
The isolates were analyzed by MLST using the sequences generated from the internal fragments of the seven HK genes. The alleles numbers and STs were assigned according to a database created for V. parahaemolyticus isolates upon submitting the sequence results. A total of 86 STs were observed among the 98 isolates. The numbers of alleles observed for each MLST locus in this study showed statistics as follows: dna E: 61; gyr B: 70; rec A: 62; dtd S: 57; pnt A: 43; pyr C: 60; and tna A: 42. The haplotype diversity was 0.987. A minimum evolutionary tree was constructed using the concatenated sequences of each allele. The MLST results divided the V. parahaemolyticus isolates into five clusters (designated as I, II, III, IV, and V). Most isolates were distributed in cluster I. Interestingly, Figure 3 shows that all the isolates in cluster III were from Fujian, whereas, cluster IV contained the V. parahaemolyticus isolates from Hainan. In cluster V, most of the isolates were from three coastal cities. For strain Vps18 and Vps06 (from Guangdong), the ST type (ST212) showed a greater distance on the evolutionary tree.
FIGURE 3
Discussion
According to epidemiologic reports, V. parahaemolyticus is a major cause of bacterial infections associated with the consumption of raw or undercooked aquatic products including different kinds of shellfish and fish (). Our results show that the prevalence of V. parahaemolyticus was higher in oyster samples (30.36%) than in the other seafood samples. This prevalence is also lower than that found in Germany (41.8%) (). The prevalence of V. parahaemolyticus isolates detected off the eastern coast of China was also reported to be higher than that observed in our study (). The reason for this discrepancy maybe that most of the sampled cities were inland, with much fewer kinds of seafood. Consistent with this hypothesis, the prevalence of this bacteria in coastal cities (Guangdong, Hainan, and Fujian) was higher (Supplementary Table S1).
Notably, the prevalence of V. parahaemolyticus in summer (33.33%) was higher than that in winter (14.01%). Our results are in accordance with previous reports that V. parahaemolyticus outbreaks typically show a seasonal pattern, peaking in the warmer months (; ), and that V. parahaemolyticus isolates are cold sensitive, thus, lower temperature in winter may limit or reduce their growth (). As our results were not only obtained from samples collected from a majority of regions in South China, but also from different seasons and kinds of seafood, the data can be useful for determining optimal temperature values for control of V. parahaemolyticus contamination ().
As there is an increase in the number of resistance genes and the spread of antimicrobial-resistant V. parahaemolyticus isolates worldwide, the misuse and overuse of antibiotics are considered the most important factors (; ). In the Hunan province of China, the average utilization rate of antibiotics reached 64.56% and more than 50 antibiotics were available for sale (). Susceptibility tests show that isolates V. parahaemolyticus in South China appear to a high level of resistance to ampicillin. This result is similar to a report by in which 82% of the isolates from shrimp samples were also resistant to ampicillin. In addition, many isolates were resistant to streptomycin (68.37%) and kanamycin (39.80%). Similarly, previous studies have suggested that resistance to these antibiotics is common in V. parahaemolyticus isolates (; ). Susceptibility profiles to antibiotic classes such as cefotaxime, cefoxitin, imipenem, and meropenem were examined. Notably, we found some isolates showing resistance to chloramphenicol, tetracycline or ciprofloxacin, which are first-line drugs used in clinical treatment (; ). In this work, 68.38% of the strains were multi-drug resistant and some even showed resistance to seven antibiotics. This ratio is higher than that detected in previous reports (; ). In general, evaluating variations in the antimicrobial susceptibility trends of V. parahaemolyticus strains is important give that infection emergence of microbial resistance to multiple drugs is a serious clinical problem.
The toxR gene, which is involved in the regulation of many other genes, was ubiquitous in V. parahaemolyticus isolates. As we know, not all V. parahaemolyticus strains are pathogenic in humans; however, the products of the hemolysin genes (tdh and trh) are believed to rapidly induce inflammatory gastroenteritis and are often detected in clinical strains (; ). Thus, detection of the hemolysin genes could be an important way to infer the virulence potential of food isolates. In our study, up to 8.16 and 12.24% of the strains were tdh or trh positive. Our findings are similar to those of previous studies on environmental V. parahaemolyticus isolates (; ). The distributions and evolution of tdh+ and/or trh+ strains may differ depending on sample source, geographical region, and other environmental factors (; ). However, the relatively high percentages of tdh and trh positive V. parahaemolyticus isolates identified in food samples represent a potential public health risk.
Recently, the utility of molecular techniques has been established in epidemiological studies of V. parahaemolyticus infections and used for the analysis of genetic diversity. As a relatively simple method, ERIC-PCR was easier to perform than PFGE. It also was appropriate for the analysis of the large population of strains (; ). Using this approach, the isolates were classified into five clusters when a 0.67 similarity cutoff value was used. However, we did not identify significant association between ERIC-PCR clusters and collection sites (source and provinces). Nevertheless, these result revealed high level of genetic diversity among V. parahaemolyticus isolates. showed a link between environment and likelihood of an outbreak (). We also obtained similar results. MLST is also a useful method for typing strains because of its reproducibility and comparability. The approach has been commonly used for analysis of V. parahaemolyticus isolates (). In this study, our isolates were grouped into five main clusters (I, II, III, IV, and V). MLST analysis revealed that the same source strains were in the same clusters (III, IV, and V) and suggested that these isolates have similar genotypes. As previously reported, the same ST types generally have the same serotype (). Multiple ST types were identified in our research. For example, ST291, ST300, and ST396 were identified in a public database as environmental isolates from China; ST1228, ST1251, and ST1231 were isolated from RTE foods; ST40 and ST162 were separated from environmental and clinical samples from the United States. There were also some clinical strains identified such as ST112 and ST739. ST212, which was isolated from clinical samples in Peru, showed greater evolutionary distance compared with other strains in our study. MLST analysis confirmed the genetic relatedness and genetic diversity occurring within those strains. Thus, ERIC-PCR and MLST are useful as phylogenetic research tools and as important methods for investigating outbreaks of V. parahaemolyticus.
Diarrhea due to the seafood borne pathogen V. parahaemolyticus has been a longstanding problem. In general, we provided the first comprehensive research on V. parahaemolyticus isolates from seafood in South China in both the summer and the winter by describing the prevalence, antibiotic resistance phenotypes, virulence, and molecular diversity. Our study showed that the prevalence of V. parahaemolyticus is 19.44%. The percentages of the isolates possessing tdh or trh were 8.16 and 12.24%, respectively. The antimicrobial-resistance patterns showed that ampicillin-resistance was widespread (79.59%). ERIC-PCR and MLST typing showed genetic relatedness and genetic diversity within these isolates. As seafood is a common and popular food choice in China, our findings should serve as a reference for understanding antibiotic-resistant trends in V. parahaemolyticus strains, providing information for evaluating this bacteria at the level of consumption, and establishing the appropriate monitoring programs, which are important in ensuring the safety of seafood.
Statements
Author contributions
YY finished the experiments and wrote the article. HY gave the idea and experimental support. JX and HL provided assistance and guidance throughout the research. ST and YC helped to finish the experiments. All authors assisted in writing the article and manuscript editing.
Acknowledgments
Our work was supported by research grants from the National Natural Science Foundation of Guangdong, China (2017A030310642) and Science and Technology Planning Project of Guangdong Province, China (2016A020210141).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/articles/10.3389/fmicb.2017.02566/full#supplementary-material
Footnotes
References
1
BanerjeeS. K.KearneyA. K.NadonC. A.PetersonC. L.TylerK.BakoucheL.et al (2014). Phenotypic and genotypic characterization of Canadian clinical isolates of Vibrio parahaemolyticus collected from 2000 to 2009.J. Clin. Microbiol.521081–1088. 10.1128/JCM.03047-13
2
ChenJ.ZhangR.QiX.ZhouB.WangJ.ChenY.et al (2017). Epidemiology of foodborne disease outbreaks caused by Vibrio parahaemolyticus during 2010–2014 in Zhejiang Province, China.Food Control77110–115. 10.1016/j.foodcont.2017.02.004
3
ChenW.XieY.XuJ.WangQ.GuM.YangJ.et al (2012). Molecular typing of Vibrio parahaemolyticus isolates from the middle-east coastline of China.Int. J. Food Microbiol.153402–412. 10.1016/j.ijfoodmicro.2011.12.001
4
ChoT.KimN.KimS.SongJ.RheeM.-S. (2016). Survival of foodborne pathogens (Escherichia coli O157: H7, Salmonella Typhimurium, Staphylococcus aureus, Listeria monocytogenes, and Vibrio parahaemolyticus) in raw ready-to-eat crab marinated in soy sauce.Int. J. Food Microbiol.23850–55. 10.1016/j.ijfoodmicro.2016.08.041
5
CLSI (2012). Methods for Antimicrobial Dilution and Disk Susceptibility Testing of Infrequently. Isolated or Fastidious Bacteria. Approved Standard M45-A.Wayne, PA: Clinical and Laboratory Standards Institute.
6
DaimeiC.DaiC. (2016). The present situation and countermeasures of antibiotic use in rural aquaculture industry - taking Zhangjiajie as an example [J].Mod. Agric. Technol.18246–250.
7
DanielsN. A.MackinnonL.BishopR.AltekruseS.RayB.HammondR. M.et al (2000). Vibrio parahaemolyticus Infections in the United States, 1973–1998.J. Infect. Dis.1811661–1666. 10.1086/315459
8
ElexsonN.Afsah-HejriL.RukayadiY.SoopnaP.LeeH. Y.ZainazorT. C. T.et al (2014). Effect of detergents as antibacterial agents on biofilm of antibiotics-resistant Vibrio parahaemolyticus isolates.Food Control35378–385. 10.1016/j.foodcont.2013.07.020
9
ElmahdiS.DasilvaL. V.ParveenS. (2016). Antibiotic resistance of Vibrio parahaemolyticus and Vibrio vulnificus in various countries: a review.Food Microbiol.57128–134. 10.1016/j.fm.2016.02.008
10
FujinoT.OkunoY.NakadaD.AoyamaA.FukaiK.MukaiT.et al (1953). On the bacteriological examination of shirasu-food poisoning.Med. J. Osaka Univ.4299–304.
11
GonzálezescalonaN.MartinezurtazaJ.RomeroJ.EspejoR. T.JaykusL. A.DepaolaA. (2008). Determination of molecular phylogenetics of Vibrio parahaemolyticus strains by multilocus sequence typing.J. Bacteriol.1902831–2840. 10.1128/JB.01808-07
12
GuerreroA.LiceanavarroA. F.WongchangI.GonzálezsánchezR. (2017). Genetic analysis of Vibrio parahaemolyticus O3:K6 strains that have been isolated in Mexico since 1998.PLOS ONE12:e0169722. 10.1371/journal.pone.0169722
13
Gutierrez WestC. K.KleinS. L.LovellC. R. (2013). High frequency of virulence factor genes tdh, trh, and tlh in Vibrio parahaemolyticus strains isolated from a pristine estuary.Appl. Environ. Microbiol.792247–2252. 10.1128/AEM.03792-12
14
HondaT.IidaT. (1993). The pathogenicity of Vibrio parahaemolyticus and the role of the thermostable direct haemolysin and related haemolysins.Rev. Med. Microbiol.4106–113. 10.1097/00013542-199304000-00006
15
JolleyK. A.ChanM.-S.MaidenM. C. (2004). mlstdbNet–distributed multi-locus sequence typing (MLST) databases.BMC Bioinformatics5:86. 10.1186/1471-2105-5-86
16
JonesJ. L.LüdekeC. H.BowersJ. C.GarrettN.FischerM.ParsonsM. B.et al (2012). Biochemical, serological, and virulence characterization of clinical and oyster Vibrio parahaemolyticus isolates.J. Clin. Microbiol.502343–2352. 10.1128/JCM.00196-12
17
JosephS. W.ColwellR. R.KaperJ. B. (1982). Vibrio parahaemolyticus and related halophilic Vibrios.Crit. Rev. Microbiol1077–124. 10.3109/10408418209113506
18
JunJ. W.KimJ. H.ChorescaC. H.Jr.ShinS. P.HanJ. E.HanS. Y.et al (2012). Isolation, molecular characterization, and antibiotic susceptibility of Vibrio parahaemolyticus in Korean seafood.Foodborne Pathog. Dis.9224–231. 10.1089/fpd.2011.1018
19
KangC. H.ShinY. J.KimW. R.KimY. G.SongK. C.OhE. G.et al (2016). Prevalence and antimicrobial susceptibility of Vibrio parahaemolyticus isolated from oysters in Korea.Environ. Sci. Pollut. Res.231–9. 10.1007/s11356-015-5650-9
20
KhanA. A.McCarthyS.WangR. F.CernigliaC. E. (2002). Characterization of United States outbreak isolates of Vibrio parahaemolyticus using enterobacterial repetitive intergenic consensus (ERIC) PCR and development of a rapid PCR method for detection of O3:K6 isolates.FEMS Microbiol. Lett.206209–214. 10.1111/j.1574-6968.2002.tb11011.x
21
KimY. B.OkudaJ.MatsumotoC.TakahashiN.HashimotoS.NishibuchiM. (1999). Identification of Vibrio parahaemolyticus strains at the species level by PCR targeted to the toxR gene.J. Clin. Microbiol.371173–1177.
22
LeoniF.TaleviG.MasiniL.OttavianiD.RocchegianiE. (2016). Trh (tdh-/trh+) gene analysis of clinical, environmental and food isolates of Vibrio parahaemolyticus as a tool for investigating pathogenicity.Int. J. Food Microbiol.22543–53. 10.1016/j.ijfoodmicro.2016.02.016
23
LetchumananV.PusparajahP.TanT. H.YinW. F.LeeL. H.ChanK. G. (2015a). Occurrence and antibiotic resistance of Vibrio parahaemolyticus from shellfish in Selangor, Malaysia.Front. Microbiol.6:1417. 10.3389/fmicb.2015.01417
24
LetchumananV.YinW. F.LeeL. H.ChanK. G. (2015b). Prevalence and antimicrobial susceptibility of Vibrio parahaemolyticus isolated from retail shrimps in Malaysia.Front. Microbiol.6:33. 10.3389/fmicb.2015.00033
25
LinZ.KumagaiK.BabaK.MekalanosJ.NishibuchiM. (1993). Vibrio parahaemolyticus has a homolog of the Vibrio cholerae toxRS operon that mediates environmentally induced regulation of the thermostable direct hemolysin gene.J. Bacteriol.1753844–3855. 10.1128/jb.175.12.3844-3855.1993
26
LopatekM.WieczorekK.OsekJ. (2015). Prevalence and antimicrobial resistance of Vibrio parahaemolyticus isolated from raw shellfish in Poland.J. Food Prot.781029–1033. 10.4315/0362-028X.JFP-14-437
27
LüdekeC. H.Gonzalez-EscalonaN.FischerM.JonesJ. L. (2015). Examination of clinical and environmental Vibrio parahaemolyticus isolates by multi-locus sequence typing (MLST) and multiple-locus variable-number tandem-repeat analysis (MLVA).Front. Microbiol.6:564. 10.3389/fmicb.2015.00564
28
MahoneyJ. C.GerdingM. J.JonesS. H.WhistlerC. A. (2010). Comparison of the pathogenic potentials of environmental and clinical Vibrio parahaemolyticus strains indicates a role for temperature regulation in virulence.Appl. Environ. Microbiol.767459–7465. 10.1128/AEM.01450-10
29
MalaW.AlamM.AngkititrakulS.WongwajanaS.LulitanondV.HuttayananontS.et al (2016). Serogroup, virulence, and molecular traits of Vibrio parahaemolyticus isolated from clinical and cockle sources in northeastern Thailand.Infect. Genet. Evol.39212–218. 10.1016/j.meegid.2016.01.006
30
MarshallS.ClarkC. G.WangG.MulveyM.KellyM. T.JohnsonW. M. (1999). Comparison of molecular methods for typing Vibrio parahaemolyticus.J. Clin. Microbiol.372473–2478.
31
MatsudaS.KodamaT.OkadaN.OkayamaK.HondaT.IidaT. (2010). Association of Vibrio parahaemolyticus thermostable direct hemolysin with lipid rafts is essential for cytotoxicity but not hemolytic activity.Infect. Immun.78603–610. 10.1128/IAI.00946-09
32
MillerR.CarsonJ.DalsgaardI.GauntP.GiesekerC.HawkeJ.et al (2014). CLSI Performance Standards for Antimicrobial Susceptibility Testing of Bacteria Isolated from Aquatic Animals; Second Information Supplement.Wayne, PA: CLSI.
33
OdeyemiO. A.StratevD. (2016). Occurrence of antimicrobial resistant or pathogenic Vibrio parahaemolyticus in seafood. A mini review.Rev. Méd. Vét.16793–98.
34
OttavianiD.LeoniF.TaleviG.MasiniL.SantarelliS.RocchegianiE.et al (2013). Extensive investigation of antimicrobial resistance in Vibrio parahaemolyticus from shellfish and clinical sources, Italy.Int. J. Antimicrob. Agents42191–193. 10.1016/j.ijantimicag.2013.05.003
35
RaghunathP. (2015). Roles of thermostable direct hemolysin (TDH) and TDH-related hemolysin (TRH) in Vibrio parahaemolyticus.Front. Microbiol.5:805. 10.3389/fmicb.2014.00805
36
RohlfF. J. (2000). Ntsys-Pc: Numerical Taxonomy and Multivariate Analysis System, Version 2.2.Setauket, NY: Exeter Publishing.
37
ShawK. S.Rosenberg GoldsteinR. E.HeX.JacobsJ. M.CrumpB. C.SapkotaA. R. (2014). Antimicrobial susceptibility of Vibrio vulnificus and Vibrio parahaemolyticus recovered from recreational and commercial areas of Chesapeake Bay and Maryland Coastal Bays.PLOS ONE9:e89616. 10.1371/journal.pone.0089616
38
ShenX.LiuW.LiuC. (2010). “Effects of cold storage and thermal treatment on growth and survival of pathogenicVibrio parahaemolyticus,” in Proceedings of the International Conference on Bioinformatics and Biomedical Technology, Chengdu, 371–373.
39
ShiraiH.ItoH.HirayamaT.NakamotoY.NakabayashiN.KumagaiK.et al (1990). Molecular epidemiologic evidence for association of thermostable direct hemolysin (TDH) and TDH-related hemolysin of Vibrio parahaemolyticus with gastroenteritis.Infect. Immun.583568–3573.
40
SuY. C.LiuC. (2007). Vibrio parahaemolyticus: a concern of seafood safety.Food Microbiol.24549–558. 10.1016/j.fm.2007.01.005
41
TakahashiA.KenjyoN.ImuraK.MyonsunY.HondaT. (2000). Cl- secretion in colonic epithelial cells induced by the Vibrio parahaemolyticus hemolytic toxin related to thermostable direct hemolysin.Infect. Immun.685435–5438. 10.1128/IAI.68.9.5435-5438.2000
42
TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0.Mol. Biol. Evol.302725–2729. 10.1093/molbev/mst197
43
TanC. W.MalcolmT. T. H.KuanC. H.ThungT. Y.ChangW. S.LooY. Y.et al (2017). Prevalence and antimicrobial susceptibility of Vibrio parahaemolyticus isolated from short mackerels (Rastrelliger brachysoma) in Malaysia.Front. Microbiol.8:1087. 10.3389/fmicb.2017.01087
44
TanL. T.ChanK. G.LeeL. H.GohB. H. (2016). Streptomyces bacteria as potential probiotics in aquaculture.Front. Microbiol.7:79. 10.3389/fmicb.2016.00079
45
WestC. K. G.KleinS. L.LovellC. R. (2013). High frequency of virulence factor genes tdh, trh, and tlh in Vibrio parahaemolyticus strains isolated from a pristine estuary.Appl. Environ. Microbiol.792247–2252. 10.1128/AEM.03792-12
46
WilsonB. A.SalyersA. A. (2003). Is the evolution of bacterial pathogens an out-of-body experience?Trends Microbiol.11347–350.
47
WongH. C.ChenM. C.LiuS. H.LiuD. P. (1999). Incidence of highly genetically diversified Vibrio parahaemolyticus in seafood imported from Asian countries.Int. J. Food Microbiol.52181–188. 10.1016/S0168-1605(99)00143-9
48
WongH. C.LinC. H. (2001). Evaluation of typing of Vibrio parahaemolyticus by three PCR methods using specific primers.J. Clin. Microbiol.394233–4240.
49
World Health Organization (2014). EN JEMRA Joint FAO/WHO Expert Meetings on Microbiological Risk Assessment.Geneva: WHO.
50
XieT.WuQ.XuX.ZhangJ.GuoW. (2015). Prevalence and population analysis of Vibrio parahaemolyticus in aquatic products from South China markets.FEMS Microbiol. Lett.362:fnv178. 10.1093/femsle/fnv178
51
XieT.WuQ.ZhangJ.XuX.ChengJ. (2017). Comparison of Vibrio parahaemolyticus isolates from aquatic products and clinical by antibiotic susceptibility, virulence, and molecular characterisation.Food Control71315–321. 10.1016/j.foodcont.2016.06.046
52
XieT.XuX.WuQ.ZhangJ.ChengJ. (2016). Prevalence, molecular characterization, and antibiotic susceptibility of Vibrio parahaemolyticus from ready-to-eat foods in China.Front. Microbiol.7:549. 10.3389/fmicb.2016.00549
53
XuX.ChengJ.WuQ.ZhangJ.XieT. (2016). Prevalence, characterization, and antibiotic susceptibility of Vibrio parahaemolyticus isolated from retail aquatic products in North China.BMC Microbiol.16:32. 10.1186/s12866-016-0650-6
54
YangZ. Q.JiaoX. A.ZhouX. H.CaoG. X.FangW. M.GuR. X. (2008). Isolation and molecular characterization of Vibrio parahaemolyticus from fresh, low-temperature preserved, dried, and salted seafood products in two coastal areas of eastern China.Int. J. Food Microbiol.125279–285. 10.1016/j.ijfoodmicro.2008.04.007
55
YanoY.HamanoK.SatomiM.TsutsuiI.BanM.Aue-UmneoyD. (2014). Prevalence and antimicrobial susceptibility of Vibrio species related to food safety isolated from shrimp cultured at inland ponds in Thailand.Food Control3830–36. 10.1016/j.foodcont.2013.09.019
56
YuQ.NiuM.YuM.LiuY.WangD.ShiX. (2015). Prevalence and antimicrobial susceptibility of Vibrio parahaemolyticus isolated from retail shellfish in Shanghai.Food Control60263–268. 10.1016/j.foodcont.2015.08.005
57
ZhaoF.ZhouD. Q.CaoH. H.MaL. P.JiangY. H. (2011). Distribution, serological and molecular characterization of Vibrio parahaemolyticus from shellfish in the eastern coast of China.Food Control221095–1100. 10.1016/j.foodcont.2010.12.017
Summary
Keywords
Vibrio parahaemolyticus, prevalence, antibiotic resistance, virulence gene, MLST
Citation
Yang Y, Xie J, Li H, Tan S, Chen Y and Yu H (2017) Prevalence, Antibiotic Susceptibility and Diversity of Vibrio parahaemolyticus Isolates in Seafood from South China. Front. Microbiol. 8:2566. doi: 10.3389/fmicb.2017.02566
Received
27 September 2017
Accepted
11 December 2017
Published
20 December 2017
Volume
8 - 2017
Edited by
Aldo Corsetti, Università di Teramo, Italy
Reviewed by
Anwar Huq, University of Maryland, College Park, United States; Ben D. Tall, U.S. Food and Drug Administration, United States
Updates
Copyright
© 2017 Yang, Xie, Li, Tan, Chen and Yu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hui Yu, yu88hui@qq.com
This article was submitted to Food Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.