Prevalence, Molecular Characterization, and Antibiotic Susceptibility of Vibrio parahaemolyticus from Ready-to-Eat Foods in China

Abstract

Vibrio parahaemolyticus is the leading cause of foodborne outbreaks, particularly outbreaks associated with consumption of fish and shellfish, and represents a major threat to human health worldwide. This bacterium harbors two main virulence factors: the thermostable direct hemolysin (TDH) and TDH-related hemolysin (TRH). Additionally, various serotypes have been identified. The extensive use of antibiotics is a contributing factor to the increasing incidence of antimicrobial-resistant V. parahaemolyticus. In the current study, we aimed to determine the incidence and features of V. parahaemolyticus in ready-to-eat (RTE) foods in China. We found 39 V. parahaemolyticus strains on Chinese RTE foods through investigation of 511 RTE foods samples from 24 cities in China. All isolates were analyzed for the presence of tdh and trh gene by PCR, serotyping was performed using multiplex PCR, antibiotic susceptibility analysis was carried out using the disk diffusion method, and molecular typing was performed using enterobacterial repetitive intergenic consensus sequence PCR (ERIC-PCR) typing and multilocus sequence typing (MLST). The results showed that none of the isolates were positive for tdh and trh. Most of the isolates (33.3%) were serotype O2. Antimicrobial susceptibility results indicated that most strains were resistant to streptomycin (89.7%), cefazolin (51.3%), and ampicillin (51.3%). The isolates were grouped into five clusters by ERIC-PCR and four clusters by MLST. We updated 10 novel loci and 33 sequence types (STs) in the MLST database. Thus, our findings demonstrated the presence of V. parahaemolyticus in Chinese RTE foods, provided insights into the dissemination of antibiotic-resistant strains, and improved our knowledge of methods of microbiological risk assessment in RTE foods.

Introduction

Vibrio parahaemolyticus is a gram-negative halophilic bacterium that naturally occurs worldwide in estuarine environments. This microorganism is recognized as one of the leading causes of foodborne illness worldwide and has been shown to cause acute gastroenteritis in humans. Complications such as septicemia can sometimes lead to death in patients with V. parahaemolyticus infection (Qadri et al., 2005; Lopatek et al., 2015). Previous studies have focused on the prevalence of V. parahaemolyticus in shellfish, oysters, water, seafood, and shrimp (Khouadja et al., 2013; Lopatek et al., 2015). However, no reports have identified isolates of V. parahaemolyticus on ready-to-eat (RTE) foods such as cooked meat, roasted poultry, and cold vegetable dishes in sauce, that are highly popular in China. Various foodborne pathogens may be present in RTE foods and may cause illnesses in consumers because RTE foods do not require heat treatment or other forms of curing before eating (Shi et al., 2015). Moreover, the high genetic diversity on RTE foods can facilitate identification of strain relatedness and epidemiological investigations.

Traditionally, V. parahaemolyticus is considered susceptible to antimicrobials. However, during the past few decades, antimicrobial resistance has emerged and evolved in many bacterial genera owing the excessive use of antimicrobials in human, agriculture, and aquaculture systems (Cabello, 2006; Kang et al., 2016). For example, tetracyclines are recommended as antibiotics in the treatment of severe Vibrio infections, and third-generation cephalosporins with doxycycline or fluoroquinolone alone are also used on occasion. Moreover, antibiotics are commonly used to treat fish (Han et al., 2007; Devi et al., 2009). Some V. parahaemolyticus isolates from seafood and other environments are commonly resistant to ampicillin, aminoglycosides (streptomycin and gentamicin), ciprofloxacin, chloramphenicol, and other antibiotics (Oh et al., 2011; Raissy et al., 2012; Shaw et al., 2014; Lopatek et al., 2015). Therefore, the potential presence of antibiotic-resistant V. parahaemolyticus on RTE foods may be an important public health problem related to disease management and control.

The virulence of V. parahaemolyticus is mainly attributed to the production of two major factors: thermo-stable direct hemolysin (TDH) encoded by the tdh gene and TDH-related hemolysin encoded by the trh gene (Honda and Iida, 1993). Clinical strains commonly contain either these genes, and the presence of these genes is associated with the pathogenicity of the strain in humans (Su and Liu, 2007; Jones et al., 2012; Pazhani et al., 2014). Detection of V. parahaemolyticus virulence factors is typically based on molecular biological analysis and amplification of V. parahaemolyticus-specific sequences (Shirai et al., 1990; West et al., 2013). To date, there are 13 O-serogroups and over 70 K-serogroups, differentiated on the basis of somatic (O) and capsular (K) antigens in V. parahaemolyticus (Ishibashi et al., 2000; Jones et al., 2012). The emergence of the first pandemic strain belonging to serovar O3:K6 (Okuda et al., 1997) supported the view that the serotype of V. parahaemolyticus is correlated with virulence. A multiplex polymerase chain reaction (PCR)-based O-antigen serotyping method has been adopted because other antiserum-based approaches are time consuming, expensive and can be associated with the risk of cross-reactions (Chen M. et al., 2012; Xu et al., 2014).

Molecular typing of V. parahaemolyticus, including pulsed field gel electrophoresis (PFGE) (Marshall et al., 1999), repetitive extragenic palindromic sequence PCR (REP-PCR) (Wong and Lin, 2001), and ribotyping (Bag et al., 1999), has been shown to be a useful tool for providing genetic relatedness information (Olive and Bean, 1999). Enterobacterial repetitive intergenic consensus sequence PCR (ERIC-PCR) has been proven useful for subtyping of V. parahaemolyticus strains with highly conserved repetitive intergenic consensus sequences (De Bruijn, 1992; Chen W. et al., 2012). Additionally, multilocus sequence typing (MLST), which is based on sequence analysis of selected housekeeping genes (e.g., recA, gyrB, dnaE, dtdS, pntA, pyrC, and tnaA), is becoming an important method for investigation of the evolution and epidemiology of V. parahaemolyticus owing to its high repeatability (González-Escalona et al., 2008; Banerjee et al., 2014).

V. parahaemolyticus is not frequently found in RTE foods, but is still considered hazardous in humans. Therefore, the objectives of the current study were as follows: (i) to determine the prevalence and contamination level of V. parahaemolyticus in RTE foods in China and (ii) to determine the genetic variation and phenotypic characteristics of V. parahaemolyticus isolates from RTE foods. The information generated in this study will provide insights into the distribution and population of V. parahaemolyticus across Chinese RTE foods and differentiation of V. parahaemolyticus strains on different RTE foods.

Materials and methods

Sample collection

From November 2011 to May 2014, a total of 511 RTE food samples were collected from retail markets in 24 cities, covering most provincial capitals of China (Figure 1). The samples consisted of 371 deli meat samples, 97 cold vegetable dishes or noodles in sauce, and 43 fried rice or noodle samples. The samples were placed in sterile sealed plastic bags, transported to the laboratory in a cold box below 4°C, and analyzed immediately.

Figure 1

Qualitative detection

The bacteriological media used herein, unless indicated, were purchased from Guangdong Huankai Co. Ltd. (Guangzhou, China). For qualitative detection, the V. parahaemolyticus was isolated according to GB 4789.7-2013 for food microbiological examination of V. parahaemolyticus (National Food Safety Standards of China) with minor modifications. In brief, 25 g of each sample was homogenized for 60 s in a stomacher bag (Huankai Co. Ltd., Guangzhou, China) with 225 mL of alkaline peptone water (APW) containing 3% NaCl. Homogenates were incubated at 37°C for 16−18 h. After incubation, a loopful from the top 1 cm was streaked onto thiosulfate-citrate-bile salts-sucrose (TCBS) agar plates and incubated at 37°C for 18−24 h. Presumptive colonies (green or blue green colonies, 2−3 mm in diameter) were streaked onto Chromogenic Vibrio Medium and incubated at 37°C for 24 h. The mauve colonies were selected for further V. parahaemolyticus identification by analysis of oxidase activity, Gram staining, 3.5% NaCl triple sugar iron (TSI) tests, halophilism tests, and API 20E diagnostic strips (Biomerieux Company, France).

MPN enumeration (MPN) of V. parahaemolyticus

Enumeration of V. parahaemolyticus in RTE food samples was performed using the MPN method according to the Bacteriological Analytical Manual standard method (Kaysner and Depaola, 2004) and a previous study (Xu et al., 2014).

Detection of tdh and trh genes

Detection of the V. parahaemolyticus tdh and trh genes was performed by PCR as described previously (Gutierrez West et al., 2013). The oligonucleotide primers were synthesized by Sangon Biotech (Shanghai, China) (Tdh-F: CTGTCCCTTTTCCTGCCCCCG, Tdh-R: AGC CAGACACCGCTGCCATTG; Trh-F: ACCTTTTCCTTCTCCWGGKTCSG, Trh-F: CCGCTC TCATATGCYTCGACAKT). Each reaction mixture included the following (total volume, 25 μL): 2 × PCR Mix (Qiagen), 12.5 μL; 0.5 μM each primer, dd H₂O, 9.5 μL; and DNA template, 1 μL. Both genes were amplified using the following thermal-cycling program: denaturation at 95°C for 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 of 72°C for 2 min. PCR was conducted in a Bio-Rad PTC-200 Thermal Cycler (Bio-Rad, Hercules, CA, USA). The amplified products were then analyzed electrophoretically on a 2% agarose gel containing Gold View. The images were captured digitally and analyzed using a Gel Image system (Bio-Rad). V. parahaemolyticus strains ATCC 33847 (tdh+) and ATCC 17802 (trh+) were used as positive controls, and distilled water was used as the negative control.

Multiplex serotyping PCR

The serotypes of V. parahaemolyticus isolates were identified using the PCR-based O-antigen serotyping technique. The primer concentrations and amplification conditions used were as described previously (Chen M. et al., 2012). The primers used for this assay were synthesized by Sangon Biotech (Shanghai, China). The 12 primer pairs were divided into two groups to amplify target DNA; PCR group 1 was used to detect serogroups O1, O2, O4, O5, O6, and O10, whereas PCR group 2 was used to detect serogroups O3/O13, O7, O8, O9, O11, and O12. The PCR was performed in a 25-μL reaction mixture containing the following: 2 × PCR mix (Dongshen, Guangzhou, China), 12.5 μL; 0.5 μM each primer, dd H₂O, 9.5 μL; and DNA template, 1 μL. All amplifications were carried out with the following protocol: 30 cycles of 95°C for 30 s, 60°C for 45 s, and 72°C for 1 min. The thermal cycler was prewarmed to 80°C before all the reaction tubes were added in order to reduce nonspecific amplification. PCR was conducted in a Bio-Rad PTC-200 Thermal Cycler (Bio-Rad, California, USA). The amplified products were then analyzed electrophoretically on a 2% agarose gel containing Gold View. The images were captured digitally and analyzed using the Gel Image system (Bio-Rad, California, USA). V. parahaemolyticus ATCC 17802 and ATCC 33847 were used as control strains.

Antimicrobial susceptibility

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 (CLSI, 2012) and a previous study (Xie et al., 2015). Briefly, Muller-Hinton agar and a panel of 12 antibiotics disks were selected for resistance tests. These 12 antibiotic disks (Oxoid, Hampshire, UK) were ampicillin (10 μg), azithromycin (15 μg), cefazolin (30 μg), cephalothin (30 μg), chloramphenicol (30 μg), ciprofloxacin (5 μg), gentamicin (10 μg), kanamycin (30 μg), nalidixic acid (30 μg), streptomycin (10 μg), trimethoprim-sulfamethoxazole (25 μg), and tetracycline (30 μg). The results were expressed as sensitive (S), intermediate (I), or resistant (R) following the methods of the CLSI. Escherichia coli ATCC 25922 and V. parahaemolyticus ATCC 17802 were used as quality control organisms.

ERIC-PCR analysis

ERIC-PCR analysis was performed on the V. parahaemolyticus isolates using a previously described protocol with some modifications (Chen W. et al., 2012; Xie et al., 2015). The reaction mixture (25 μL per reaction) consisted of 12.5 μL 2 × Long Taq Mix (Dongsheng Biotech, Guangzhou, China), 0.6 μM of each primer (5′-ATGTAAGCTCCTGGGGATTCAC-3′ and 5′-AAGTAAGTGACTGGGGTGAGCG-3′), and 100 ng template DNA. PCR was performed in a DNA thermocycler (Applied Biosystems, CA, USA) using the following procedure: one cycle of denaturation at 95°C for 5 min; 35 cycles of 94°C for 45 s, 52°C for 1 min, and 72°C for 3 min; and a final extension at 72°C for 10 min. The PCR products were separated by electrophoresis on 2.0% agarose gels followed by Goldview staining (0.005%, v/v; SBS Genetech, Beijing) and imaging with a UV Imaging System (GE Healthcare, WI, USA). The images were captured in TIFF file format for further analysis.

The size of each band in the ERIC patterns was determined, and the data were coded as 0 (absence) or 1 (presence). Cluster analysis was performed with NTSYS-pc (Version 2.10), a numerical taxonomy and multivariate analysis software package (Rohlf, 2000), based on Dice's similarity coefficient (SD), with a 1% position tolerance and the unweighted pair group method using arithmetic averages (UPGMA).

MLST

MLST analysis was conducted via the V. parahaemolyticus MLST website and database (http://pubmlst.org/vparahaemolyticus/) (Jolley et al., 2004). PCR conditions were denaturation at 96°C for 1 min; primer (Table S1; synthesized by Sangon Biotech, Shanghai, China) annealing at 58°C for 1 min; and extension at 72°C for 1 min, for 30 cycles; with a final extension step at 72°C for 10 min. PCR was performed with a Bio- Rad PTC-200 Thermal Cycler (Bio-Rad, California, USA) according to the manufacturer's directions. PCR products were sequenced on a BGI instrument (Shenzhen, China). The alignments of these sequences were determined using BioEdit. The sequences were analyzed online (http://pubmlst.org/vparahaemolyticus/) to assign allele numbers and define STs. New sequences for alleles and new ST profiles were submitted to the V. parahaemolyticus MLST database.

The evolution tree of the concatenated sequences of the seven loci was built based on the method of the Kimura-2-parameter in Mega 6.0 (Tamura et al., 2013). The ratio between the number of synonymous and nonsynonymous substitutions, showing the type of selection at each locus, was calculated using the method of Nei and Gojobori in Mega 6. The hypotheses of neutrality (dS = dN), purifying selection (dS/dN >1), and positive selection (dS/dN < 1) were tested using DNAsp 5.10 (Lüdeke et al., 2015).

Results

V. parahaemolyticus in RTE food samples

Of the 511 samples tested, eight were positive by both qualitative and MPN methods; 12 showed positive results by the qualitative method only, while 10 were positive with the MPN method only. Thirty (5.9%) samples positive for V. parahaemolyticus were detected among 511 samples after qualitative and MPN analyses, including 22 (5.9%) of the 371 deli meat samples, seven (7.2%) of the 97 samples of cold vegetable dishes or noodles in sauce, and one (2.3%) of the 43 fried rice or noodle samples. In 18 positive samples detected by the MPN method, V. parahaemolyticus densities ranged between 3 and 100 MPN/g in all of the samples. The densities of rest positive samples were below 3.0 MPN/g. According to the National Food Safety Standards of China (GB 2727-2005, GB2726-2005), the pathogenic bacteria should not be detected in deli meat such as cooked meat or roasted meat.

Detection of tdh and trh genes in V. parahaemolyticus isolates

A total of 39 V. parahaemolyticus isolates were confirmed and tested for the presence of the trh and tdh genes. None of the isolates possessed the tdh and trh genes.

Serotyping by multiplex PCR

With the exception of serotypes O9, O10, and O11, all other serotypes were detected among the isolates. Serotype O2 was the most prevalent (13 isolates), followed by serotype O4 (eight isolates). The results of the O-antigen serotyping for all 39 isolates are shown in Table 1. The serotypes of V. parahaemolyticus ATCC 17802 and ATCC 33847 were O1 and O4, respectively (Figure S2).

Antimicrobial susceptibility

The susceptible, intermediate, and resistance rates of the 39 examined V. parahaemolyticus isolates with respect to 12 antibiotics are shown in Table 2. The resistance to streptomycin (S), ampicillin (AMP), cefazolin (KZ), cephalothin (KF), kanamycin (K), gentamicin (CN), and trimethoprim-sulfamethoxazole (SXT) reached 89.7, 51.3, 51.3, 41.0, 41.0, 2.6, and 2.6%, respectively. None of the strains showed resistance to azitromycin (AZM), chloramphenicol (C), ciprofloxacin (CIP), nalidixic acid (NA), or tetracycline (TE). Most of the isolates were resistant to streptomycin, with resistance and intermediate rates of 89.3 and 10.3%, respectively, which was consistent with a previous study (Xie et al., 2015). The next-highest susceptible rates were observed for trimethoprim-sulfamethoxazole (87.1%) and gentamicin (51.3%). In addition, there were four multidrug-resistant isolates showing resistance to five antibiotics, of which one isolate was collected from deli meat in Xiamen, two isolates were obtained from deli meat in Chengdu, and one isolate was obtained from deli meat in Jinan (Table S2).

ERIC-PCR

The results of ERIC-PCR analysis of the 39 isolates are shown in Figure 2. ERIC-PCR resulted in 3–8 amplification bands, with sizes ranging from 130 to about 6000 bp. The ERIC image shows that bands with molecular sizes of 500, 1500, and 2500 bp were common (Figure S1). All the isolates were classified into five clusters (designated as A, B, C, D, and E). Most isolates were distributed on cluster A. Only one strain (no. 34) from Zhanjiang divided into cluster D. The standard strain ATCC 17802 was in cluster E alone. The other standard strain (ATCC 33847) was belonged on cluster A and on the same sub-cluster along with isolates nos. 46 and 47.

Figure 2

MLST

All V. parahaemolyticus isolates were analyzed by MLST using the sequences generated from internal fragments of seven HK genes. Numbers for alleles and sequence types (STs) were assigned according to the database created for V. parahaemolyticus on submitting the sequence results (Table 3). The potential novel loci would appear some mismatching base. Seven gene locus mismatching the ST types on database maybe the novel one. The administrator of the website was notified when we obtained novel loci and ST types; after confirmation, new numbers will be given. The numbers of alleles observed for each MLST locus in our study were distributed as follows: 22 (dna E), 22 (gyr B), 27 (rec A), 20 (dtd S), 18 (pnt A), 24 (pyr C), and 18 (tna A). Thirty-three STs were observed among the 39 isolates. There were ten novel loci: dna E 302, 303; gyr B 428; rec A 306, 307; dtd S 371; pnt A 213; pyr C 339, 340; and tna A 238. Most of the isolates (26) were novel STs, namely 1228–1251, 1301, and 1302. Among all 3682 sites, the number of variable sites was 230; the total number of mutations was 245; the haplotype number was 33; and the haplotype diversity was 0.989. All loci showed ratios of synonymous and nonsynonymous substitutions (dN/dS) below 1. A minimum evolution tree was constructed using the concatenated sequences of each allele (Figure 3).

Figure 3

Discussion

V. parahaemolyticus is a major seafood-borne gastroenteritis-causing bacterium that is frequently isolated from aquatic products (Letchumanan et al., 2015b). In our study, we analyzed 511 samples from Chinese RTE foods and isolated 39 V. parahaemolyticus strains. In our study, some samples only detected positive in the qualitative method. It may be related to the variation in culture concentration. Enrichment culture was streaked onto TCBS agar plates with undiluted liquid in qualitative method, while culture represent 1/10, 1/100, 1/1000 dilution were streaked onto TCBS agar plates in MPN method. For some samples were only positive by the MPN method, the reason maybe attribute to the incomplete homogeneous of the sample culture solution. When a loopful (10 μL) of bacterial suspension was streaked onto the TCBS agar plate could not certain to form colonies on the agar. A combination of qualitative and MPN methods avoid missing detection of V. parahaemolyticus effectively. Unlike Listeria monocytogenes (Shi et al., 2015) and Salmonella (Yang et al., 2016), V. parahaemolyticus was not frequently detected on RTE foods. Previous study reported that no V. parahaemolyticus was positive in 145 samples of RTE food in Korea (Chung et al., 2010), similar result was found in Iran (Zarei et al., 2012). However, a report indicated that the prevalence of V. parahaemolyticus in RTE foods was an important cause of food poisoning in Shanghai, China (Tian et al., 2008). As we know, RTE foods do not need further processing before consumption. Thus, identification of contamination in RTE foods is critical for assuring food safety. Our research can provide insights into food safety in RTE foods.

TDH and TRH are considered major virulence factors in V. parahaemolyticus (Ceccarelli et al., 2013; Letchumanan et al., 2014; Raghunath, 2014). Moreover, the presence of tdh- and/or trh-positive V. parahaemolyticus strains represents a major public health risk. In our study, strains identified in Chinese RTE foods strains were negative for both trh and tdh genes. These results are consistent with the findings of a previous study reported in India (Raghunath et al., 2008), but this was contradictory to most findings from other previous studies (Zhao et al., 2011; Letchumanan et al., 2015a). The overall mechanism of V. parahaemolyticus pathogenesis remains unclear (Ceccarelli et al., 2013); although TDH and TRH have been shown to be correlated with pathogenic strains, they do not fully explain the pathogenicity of V. parahaemolyticus (Lynch et al., 2005). Several studies have reported the presence of clinical strains without tdh and trh (Shirai et al., 1990). Thus, even in the absence of tdh or trh, V. parahaemolyticus still remains pathogenic, and some environmental isolates lacking tdh and/or trh can produce putative virulence factors. For example, some oyster isolates will contain T3SS1 genes without tdh and/or trh (Mahoney et al., 2010; Jones et al., 2012).

In 1996, pandemic O3:K6 serovars were shown to be responsible for V. parahaemolyticus outbreaks. The serotype O3:K6 V. parahaemolyticus emerged from India and spread throughout the world, including to China and the USA (Okuda et al., 1997; Honda et al., 2008; Tan et al., 2010). Our study indicated that serovar O2 was the predominant serotype among the strains isolated from RTE foods. These findings are in contrast with the results of a previous study, in which the O3 serotype was identified as the predominant serotype in China (Zhang et al., 2006; Zhao et al., 2011). O3:K6, O1:Kut, O4:K8, and O2:K3 were also the dominant serovars identified in outbreaks of V. parahaemolyticus in China (Zhang et al., 2013; Ma et al., 2014). The O3 isolate was the same as the most frequent serotype among our clinical isolates (source: Shenzhen Centres for Disease Control). We also showed the presence of other O-type serovars of V. parahaemolyticus from RTE foods; the results demonstrated the diverse distributions in different RTE food types and locations in China. The changes in pandemic serogroups of V. parahaemolyticus have been reported to occur over time, an increasing number of nonpandemic serogroups have been shown to carry pandemic marker genes (Matsumoto et al., 2000). Therefore, these O2 strains from Chinese RTE foods may have the potential for pathogenesis in humans.

With the steady expansion of the Asian aquaculture industry, in order to increase production, aquaculture farmers are using different antibiotics to prevent (prophylactic use) and treat (therapeutic use) pathogenic bacterial infections in aquatic products (Cabello et al., 2013; Huang et al., 2015; Tan et al., 2016). Furthermore, the continuous and extensive use of antibiotics in humans, has led to the emergence of antimicrobial-resistant V. parahaemolyticus strains worldwide (Sani et al., 2013; Yano et al., 2014). In our study, the highest resistance rate was observed for streptomycin. Similarly, previous study have demonstrated the occurrence of streptomycin- and ampicillin-resistant V. parahaemolyticus isolates (Pazhani et al., 2014). We also found a small number of isolates showing resistance to gentamicin and trimethoprim-sulfamethoxazole, which are first-line drugs used in clinic treatment. Moreover, we found that most of the strains (21/39 isolates) were multidrug resistant. As RTE foods are eaten without cooking, the presence of these strains will increase the health risks of consuming such foods in humans. Thus, it may be important to evaluate variations in antimicrobial susceptibility profiles in V. parahaemolyticus strains.

Molecular subtyping is widely used for the analysis of genetic diversity. ERIC-PCR provides discriminatory values and can be used for rapid V. parahaemolyticus typing (Khan et al., 2002). Compared to V. parahaemolyticus isolates from the 1997 Canadian outbreak using three method, they found that ERIC-PCR and ribotyping were the most informative typing methods (Marshall et al., 1999). Using this approach, the isolates were classified into five clusters at 0.65 similarity. Most of the strains were on cluster A; which showed they were may be genetically related. The reference strain ATCC 17804 harbors trh grouped into a single cluster E, exhibiting differences with other strains. Strain (no. 34), from Zhanjiang, was grouped into cluster D alone and was genetically diverse from other isolates. Clustering based on ERIC-PCR did not coincide with the isolation sources or patterns of antibiotic resistance. The result was similar to other studies, showing the high genetic diversity in V. parahaemolyticus strains. MLST is a good method for typing owing to its reproducibility, as shown in the sequencing of seven housekeeping genes; this method has been widely used for analysis of the V. parahaemolyticus sequence (Lüdeke et al., 2015). In our study, all the isolates could be grouped into four main clusters (I, II, III, and IV). ST 291, ST 396, and ST 300 were identified in a public database as environmental isolates from China; ST1013 was separated from environmental samples from USA. Thus, some V. parahaemolyticus isolates from RTE foods were similar to environment strains. Additionally, one ST162 strain was reported as a clinical isolate from the USA. With the identification of 10 novel loci and 22 novel STs, this study substantially contributed to the diversity in the MLST database. As most strains identified in RTE foods V. parahaemolyticus stains were novel STs, these results suggested that V. parahaemolyticus strains found in RTE foods may be distinct from other stains. In our study, both ERIC-PCR and MLST confirmed the genetic diversity within strains.

In summary, this is the first comprehensive study that described the prevalence, serotype, virulence, antibiotic resistance phenotype, and molecular subtyping of V. parahaemolyticus from RTE foods in China. This study showed that none of the isolates possessed tdh and trh, and serotype O2 was found to be prevalent. The antimicrobial-resistance patterns revealed that the streptomycin-resistant were widespread (89.7%) and the isolates resistance to some clinical antibiotics such as cephalothin, gentamicin. ERIC-PCR and MLST typing showed genetic diversity. The novel loci and STs indicated genetically diverse on RTE foods isolates. As RTE foods are common and popular food choices in China, therefore the continuous monitoring of food-borne pathogens such as V. parahaemolyticus are vital to ensure the safety of these food products.

References

• 1

  BagP. K.NandiS.BhadraR. K.RamamurthyT.BhattacharyaS. K.NishibuchiM.et al. (1999). Clonal diversity among recently emerged strains of Vibrio parahaemolyticus O3: K6 associated with pandemic spread. J. Clin. Microbiol.37, 2354–2357.

  • Pubmed Abstract
  • Google Scholar

• 2

  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.52, 1081–1088. 10.1128/JCM.03047-13

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  CabelloF. C. (2006). Heavy use of prophylactic antibiotics in aquaculture: a growing problem for human and animal health and for the environment. Environ. Microbiol.8, 1137–1144. 10.1111/j.1462-2920.2006.01054.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  CabelloF. C.GodfreyH. P.TomovaA.IvanovaL.DölzH.MillanaoA.et al. (2013). Antimicrobial use in aquaculture re-examined: its relevance to antimicrobial resistance and to animal and human health. Environ. Microbiol.15, 1917–1942. 10.1111/1462-2920.12134

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  CeccarelliD.HasanN. A.HuqA.ColwellR. R. (2013). Distribution and dynamics of epidemic and pandemic Vibrio parahaemolyticus virulence factors. Front. Cell. Infect. Microbiol.3:97. 10.3389/fcimb.2013.00097

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  ChenM.GuoD.WongH.-C.ZhangX.LiuF.ChenH.et al. (2012). Development of O-serogroup specific PCR assay for detection and identification of Vibrio parahaemolyticus. Int. J. Food Microbiol.159, 122–129. 10.1016/j.ijfoodmicro.2012.08.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  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.153, 402–412. 10.1016/j.ijfoodmicro.2011.12.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  ChungM. S.KimC. M.HaS. (2010). Detection and enumeration of microorganisms in ready-to-eat foods, ready-to-cook foods and fresh-cut produce in Korea. J. Food Saf.30, 480–489. 10.1111/j.1745-4565.2010.00221.x

  • CrossRef
  • Google Scholar

• 9

  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.

  • Pubmed Abstract
  • Google Scholar

• 10

  De BruijnF. J. (1992). Use of repetitive (repetitive extragenic palindromic and enterobacterial repetitive intergeneric consensus) sequences and the polymerase chain reaction to fingerprint the genomes of Rhizobium meliloti isolates and other soil bacteria. Appl. Environ. Microbiol.58, 2180–2187.

  • Pubmed Abstract
  • Google Scholar

• 11

  DeviR.SurendranP.ChakrabortyK. (2009). Antibiotic resistance and plasmid profiling of Vibrio parahaemolyticus isolated from shrimp farms along the southwest coast of India. World J. Microbiol. Biotechnol.25, 2005–2012. 10.1007/s11274-009-0101-8

  • CrossRef
  • Google Scholar

• 12

  González-EscalonaN.Martinez-UrtazaJ.RomeroJ.EspejoR. T.JaykusL.-A.DepaolaA. (2008). Determination of molecular phylogenetics of Vibrio parahaemolyticus strains by multilocus sequence typing. J. Bacteriol.190, 2831–2840. 10.1128/JB.01808-07

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 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.79, 2247–2252. 10.1128/AEM.03792-12

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  HanF.WalkerR. D.JanesM. E.PrinyawiwatkulW.GeB. (2007). Antimicrobial susceptibilities of Vibrio parahaemolyticus and Vibrio vulnificus isolates from Louisiana Gulf and retail raw oysters. Appl. Environ. Microbiol.73, 7096–7098. 10.1128/AEM.01116-07

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  HondaT.IidaT. (1993). The pathogenicity of Vibrio parahaemolyticus and the role of the thermostable direct haemolysin and related haemolysins. Rev. Med. Microbiol.4, 106–113. 10.1097/00013542-199304000-00006

  • CrossRef
  • Google Scholar

• 16

  HondaT.IidaT.AkedaY.KadamaT. (2008). Sixty Years of Vibrio parahaemolyticus Research-Despite decades of research following a severe foodborne outbreak in Japan, this pathogen remains mysterious at the molecular level. Microbe3, 462.

  • Google Scholar

• 17

  HuangY.ZhangL.TiuL.WangH. H. (2015). Characterization of antibiotic resistance in commensal bacteria from an aquaculture ecosystem. Front. Microbiol.6:00914. 10.3389/fmicb.2015.00914

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  IshibashiM.OhtaK.ShimadaT.HondaT.SugiyamaJ.MiwataniT.et al. (2000). Current status of OK serotype combinations of Vibrio parahaemolyticus. Nippon Saikingaku Zasshi.55, 539–541. 10.3412/jsb.55.539

  • CrossRef
  • Google Scholar

• 19

  JolleyK. A.ChanM.-S.MaidenM. C. (2004). mlstdbNet–distributed multi-locus sequence typing (MLST) databases. BMC Bioinformatics5:86. 10.1186/1471-2105-5-86

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  JonesJ. L.LüdekeC. H. M.BowersJ. C.GarrettN.FischerM.ParsonsM. B.et al. (2012). Biochemical, serological, and virulence characterization of clinical and oyster Vibrio parahaemolyticus isolates. J. Clinl. Microbiol.50, 2343–2352. 10.1128/JCM.00196-12

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  KangC.-H.ShinY.KimW.KimY.SongK.OhE.-G.et al. (2016). Prevalence and antimicrobial susceptibility of Vibrio parahaemolyticus isolated from oysters in Korea. Environ.Sci. Pollut. Res.23, 918–926. 10.1007/s11356-015-5650-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  KaysnerC. A.DepaolaA. J. (2004). U. S. Food and Drug Administration; Bacteriological Analytical Manual; Methods for Specific Pathogens; Chapter 9 Vibrio. Available online at http://www.fda.gov/food/foodscienceresearch/laboratorymethods/ucm070830.htm

  • Google Scholar

• 23

  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.206, 209–214. 10.1111/j.1574-6968.2002.tb11011.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  KhouadjaS.SuffrediniE.SpagnolettiM.CrociL.ColomboM. M.AminaB. (2013). Presence of pathogenic Vibrio Parahaemolyticus in waters and seafood from the Tunisian Sea. World J. Microbiol. Biotechnol.29, 1341–1348. 10.1007/s11274-013-1297-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  LetchumananV.ChanK.-G.LeeL.-H. (2014). Vibrio parahaemolyticus: a review on the pathogenesis, prevalence, and advance molecular identification techniques. Front. Microbiol.5:705. 10.3389/fmicb.2014.00705

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  LetchumananV.PusparajahP.HernT.LohT.YinW. F.LeeL.-H.et al. (2015a). Occurrence and antibiotic resistance of Vibrio parahaemolyticus from shellfish in Selangor, Malaysia. Front. Microbiol.6:1417. 10.3389/fmicb.2015.01417

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  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

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  LopatekM.WieczorekK.OsekJ. (2015). Prevalence and antimicrobial resistance of Vibrio parahaemolyticus Isolated from Raw Shellfish in Poland. J. Food Protect.78, 1029–1033. 10.4315/0362-028X.JFP-14-437

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  LüdekeC.González-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

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  LynchT.LivingstoneS.BuenaventuraE.LutterE.FedwickJ.BuretA. G.et al. (2005). Vibrio parahaemolyticus disruption of epithelial cell tight junctions occurs independently of toxin production. Infect. Immun.73, 1275–1283. 10.1128/IAI.73.3.1275-1283.2005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  MaC.DengX.KeC.HeD.LiangZ.LiW.et al. (2014). Epidemiology and etiology characteristics of foodborne outbreaks caused by Vibrio parahaemolyticus during 2008–2010 in Guangdong province, China. Foodborne Pathog. Dis.11, 21–29. 10.1089/fpd.2013.1522

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  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.76, 7459–7465. 10.1128/AEM.01450-10

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  MarshallS.ClarkC. G.WangG.MulveyM.KellyM. T.JohnsonW. M. (1999). Comparison of Molecular Methods for Typing Vibrio parahaemolyticus. J. Clin. Microbiol.37, 2473–2478.

  • Pubmed Abstract
  • Google Scholar

• 34

  MatsumotoC.OkudaJ.IshibashiM.IwanagaM.GargP.RammamurthyT.et al. (2000). Pandemic spread of an O3: K6 clone of Vibrio parahaemolyticus and emergence of related strains evidenced by arbitrarily primed PCR and toxRS sequence analyses. J. Clin. Microbiol.38, 578–585.

  • Pubmed Abstract
  • Google Scholar

• 35

  OhE.-G.SonK.-T.YuH.LeeT.-S.LeeH.-J.ShinS.et al. (2011). Antimicrobial resistance of Vibrio parahaemolyticus and Vibrio alginolyticus strains isolated from farmed fish in Korea from 2005 through 2007. J. Food Prot.74, 380–386. 10.4315/0362-028X.JFP-10-307

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  OkudaJ.IshibashiM.HayakawaE.NishinoT.TakedaY.MukhopadhyayA. K.et al. (1997). Emergence of a unique O3: K6 clone of Vibrio parahaemolyticus in Calcutta, India, and isolation of strains from the same clonal group from Southeast Asian travelers arriving in Japan. J. Clin. Microbiol.35, 3150–3155.

  • Pubmed Abstract
  • Google Scholar

• 37

  OliveD. M.BeanP. (1999). Principles and applications of methods for DNA-based typing of microbial organisms. J. Clin. Microbiol.37, 1661–1669.

  • Pubmed Abstract
  • Google Scholar

• 38

  PazhaniG. P.BhowmikS. K.GhoshS.GuinS.DuttaS.RajendranK.et al. (2014). Trends in the epidemiology of pandemic and non-pandemic strains of Vibrio parahaemolyticus isolated from diarrheal patients in Kolkata, India. PLoS. Negl. Trop. Dis.8:2815. 10.1371/journal.pntd.0002815

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  QadriF.ChowdhuryN. R.TakedaY.NairG. B. (2005). Vibrio parahaemolyticus—seafood safety and associations with higher organisms, in Oceans and Health Pathogens in the Marine Environment, ed BelkinS. S.ColwellR. R. (New York, NY: Springer), 277–295.

  • Google Scholar

• 40

  RaghunathP. (2014). Roles of thermostable direct hemolysin (TDH) and TDH-related hemolysin (TRH) in Vibrio parahaemolyticus. Front. Microbiol.5:805. 10.3389/fmicb.2014.00805

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  RaghunathP.AcharyaS.BhanumathiA.KarunasagarI.KarunasagarI. (2008). Detection and molecular characterization of Vibrio parahaemolyticus isolated from seafood harvested along the southwest coast of India. Food Microbiol.25, 824–830. 10.1016/j.fm.2008.04.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  RaissyM.MoumeniM.AnsariM.RahimiE. (2012). Antibiotic resistance pattern of some Vibrio strains isolated. Iran. J. Fish. Sci.11, 618–626.

  • Google Scholar

• 43

  RohlfF. J. (2000). NTSYS-pc. Numerical Taxonomy and Multi-Variate Analysis System, Version 2.1. Setauket, NY: Exeter Software.

  • Google Scholar

• 44

  SaniN. A.AriyawansaS.BabjiA. S.HashimJ. K. (2013). The risk assessment of Vibrio parahaemolyticus in cooked black tiger shrimps (Penaeus monodon) in Malaysia. Food Control31, 546–552. 10.1016/j.foodcont.2012.10.018

  • CrossRef
  • Google Scholar

• 45

  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

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  ShiW.QingpingW.JumeiZ.MoutongC.ZéanY. (2015). Prevalence, antibiotic resistance and genetic diversity of Listeria monocytogenes isolated from retail ready-to-eat foods in China. Food Control47, 340–347. 10.1016/j.foodcont.2014.07.028

  • CrossRef
  • Google Scholar

• 47

  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.58, 3568–3573.

  • Pubmed Abstract
  • Google Scholar

• 48

  SuY.-C.LiuC. (2007). Vibrio parahaemolyticus: a concern of seafood safety. Food Microbiol.24, 549–558. 10.1016/j.fm.2007.01.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol.30, 2725–2729. 10.1093/molbev/mst197

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  TanH.LiB.DengX.FangW.KeB.KeC. (2010). Molecular characteristics of Vibrio parahaemolyticus isolated from food-poisoning events in Guangdong Province, 2008. South China J. Preven. Med.36, 23–27.

  • Google Scholar

• 51

  TanL. T.-H.ChanK.-G.LeeL.-H.GohB.-H. (2016). Streptomyces bacteria as potential probiotics in aquaculture. Front. Microbiol.7:79. 10.3389/fmicb.2016.00079

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  TianM.-S.ZhengL.-J.PengS.-J.LiJ.ZhangL.QiL.-B.et al. (2008). Analysis and measures on collective food poisoning caused by Vibrio parahaemolyticus in Shanghai Municipality During 2000-2007 [J]. China J. Food Hyg.6, 017.

  • Google Scholar

• 53

  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.79, 2247–2252. 10.1128/AEM.03792-12

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  WongH.-C.LinC.-H. (2001). Evaluation of typing of Vibrio parahaemolyticus by three PCR methods using specific primers. J. Clin. Microbiol.39, 4233–4240. 10.1128/JCM.39.12.4233-4240.2001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  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

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  XuX.WuQ.ZhangJ.ChengJ.ZhangS.WuK. (2014). Prevalence, pathogenicity, and serotypes of Vibrio parahaemolyticus in shrimp from Chinese retail markets. Food Control46, 81–85. 10.1016/j.foodcont.2014.04.042

  • CrossRef
  • Google Scholar

• 57

  YangX.HuangJ.WuQ.ZhangJ.LiuS.GuoW.et al. (2016). Prevalence, antimicrobial resistance and genetic diversity of Salmonella isolated from retail ready-to-eat foods in China. Food Control60, 50–56. 10.1016/j.foodcont.2015.07.019

  • CrossRef
  • Google Scholar

• 58

  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 Control38, 30–36. 10.1016/j.foodcont.2013.09.019

  • CrossRef
  • Google Scholar

• 59

  ZareiM.MaktabiS.GhorbanpourM. (2012). Prevalence of Listeria monocytogenes, Vibrio parahaemolyticus, Staphylococcus aureus, and Salmonella spp. in seafood products using multiplex polymerase chain reaction. Foodborne Pathog. Dis.9, 108–112. 10.1089/fpd.2011.0989

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  ZhangH.SunS.ShiW.CuiL.GuQ. (2013). Serotype, virulence, and genetic traits of foodborne and clinical Vibrio parahaemolyticus isolates in Shanghai, China. Foodborne Pathog. Dis.10, 796–804. 10.1089/fpd.2012.1378

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  ZhangW.PanJ.MengD.ChenK. (2006). [Molecular typing on Vibrio parahaemolyticus isolates from Hangzhou, China, during 2000-2002]. Zhonghua Liu Xing Bing Xue Za Zhi27, 343–346.

  • Pubmed Abstract
  • Google Scholar

• 62

  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 Control22, 1095–1100. 10.1016/j.foodcont.2010.12.017

  • CrossRef
  • Google Scholar