Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2008, p. 265-274, Vol. 74, No. 1
0099-2240/08/$08.00+0 doi:10.1128/AEM.01307-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Instituto de Acuicultura,1 Instituto de Investigaciones Tecnológicas, Universidad de Santiago de Compostela, Campus Universitario Sur, 15782 Santiago de Compostela, Spain,2 Instituto Tecnolóxico para o Control do Medio Mariño de Galicia, Peirao de Vilaxoán, 36611 Vilagarcia de Arousa, Spain3
Received 12 June 2007/ Accepted 22 October 2007
|
|
|---|
|
|
|---|
V. parahaemolyticus infections were rarely reported in Europe before 1998 (2). The presence of V. parahaemolyticus was reduced to sporadic cases reported in different countries without any evidence of epidemiological connection. The epidemiology of the organism in Europe changed significantly when a large number of illnesses associated with V. parahaemolyticus were reported in Galicia (in northwestern Spain) in 1999 and 2000. Isolates obtained from the outbreaks in Galicia and from hospitals in other regions of Spain were characteristically tdh positive, serotype O4:K11, and belonged to a distinctive clone of V. parahaemolyticus not related to other clinical clones that are dominant in other parts of the world (22). Nevertheless, pathogenic isolates of this clone have never been directly recovered from environmental sources or seafood from coastal areas close to the places of emergence of the outbreaks, which makes the environmental origin of the Spanish clinical isolates uncertain.
Despite the importance of the emergence of V. parahaemolyticus infections in Galicia, little is known about the presence of the organism in the marine environments of the region. Galicia is one of the most important shellfish-producing regions in Europe, with extensive harvesting of mollusks in the estuarine portions of the coasts, called rias. In the present study, we examined the presence and abundance of V. parahaemolyticus in the four main rias in Galicia and evaluated the ecological aspects related to the dynamics of the organism in the coastal areas of this temperate Atlantic region. Additionally, V. parahaemolyticus isolates were subjected to intensive analysis to assess the presence of pathogenic specimens in environmental sources.
|
|
|---|
![]() View larger version (30K): [in a new window] |
FIG. 1. Area of study and locations of the sampling stations in the coastal areas of Galicia.
|
Bacteriological analyses.
One kilogram of mussels was collected at each sampling site. The samples were placed in sterile bags and transported to the laboratory for analysis. In the laboratory, the mussels were immediately removed from the bags and washed in running potable water. Dead mussels or those with broken shells were discarded.
Enumeration of V. parahaemolyticus was performed according to the American Public Health Association three-tube most probable number (MPN) procedure (10). Fifty grams of shell liquid and meat was collected in a sterile jar and added to 450 ml of phosphate-buffered saline to make a 1:10 dilution, and the mixture was blended for 120 s with a stomacher. The shellfish homogenate was added to alkaline peptone water in a three-tube MPN dilution series. Due to the low densities of V. parahaemolyticus expected in temperate and cold waters like those of Galicia, 10 ml of one 1:10 dilution (containing 1 g of shellfish) and 1 ml of the 1:10 and 1:100 dilutions were selected to inoculate the tubes of three MPN series. The tubes were incubated at 35°C for 24 h. After incubation, a 3-mm loopful from the top 1 cm of each broth tube showing growth was streaked onto two plates of thiosulfate-citrate-bile salt-sucrose (TCBS) agar (Oxoid, Hampshire, United Kingdom). The TCBS plates were incubated at 35°C for 24 h. At least three typical sucrose-negative colonies from each plate were isolated and subjected to identification by biochemical tests on API 20E strips (bioMérieux, Marcy-l'Etoile, France).
Investigation of tlh, tdh, and trh genes.
Isolates were confirmed by the presence of the species-specific gene tlh. The presence of tlh and the virulence genes tdh and trh were determined by multiplex PCR according to the procedure described by Bej et al. (4). For DNA extraction, the isolate was cultured overnight on a tryptone soy agar plate containing 3% NaCl at 37°C. Several well-grown colonies were chosen and resuspended in 300 µl sterile distilled water and boiled for 15 min to lyse the cells. The lysate was centrifuged, and the supernatant containing DNA was used directly as a template in the PCR. PCRs were carried out in a PTC200 thermocycler (MJ Research, Waltham, MA) with the following reaction conditions: denaturation at 94°C for 3 min, followed by 30 cycles of denaturation at 94° for 1 min, primer annealing at 58°C for 1 min, and primer extension at 72°C for 1 min. A final extension was performed at 75°C for 5 min. The amplicons were analyzed in a 1.8% agarose gel.
PFGE.
Environmental isolates of V. parahaemolyticus were compared with human isolates obtained from clinical sources in Spain (22). Pulsed-field gel electrophoresis (PFGE) was performed according to the 1-day (24- to 28-h) standardized laboratory protocol for molecular subtyping of nontyphoidal salmonella by PFGE (7) following a previously described method (22). Chromosomal DNA was digested with 30 U of NotI (Promega, Southampton, United Kingdom) at 37°C for 4 h. DNA macrorestriction fragments were resolved on 1% SeaKem Gold Agarose (Cambrex, Baltimore, MD) in 0.5x Tris-borate-EDTA buffer. DNA from Salmonella sp. strain Braenderup H9812 restricted with 50 U of XbaI (Promega, Madison, WI) at 37°C for 2 h was used as a size marker. Pulse times were ramped from 2 to 40 s during an 18-h run at 6.0 V/cm. Restriction patterns were compared by the use of BioNumerics software (Applied Maths, Sint-Martens-Latem, Belgium).
Serotyping.
Lipopolysaccharide (O) and capsular (K) serotypes of the tdh-positive isolates were determined by agglutination tests with specific antisera according to the manufacturer's instructions (Denka-Seiken Ltd., Tokyo, Japan).
Environmental parameters.
The environmental parameters included in the study were air temperature, rain, wind, hours of sunshine, rainfall, river flow, upwelling, chlorophyll, phytoplankton counts, seawater temperature, and salinity. The daily ambient temperature was taken as the average of temperatures registered in a day. Wind direction was measured as the time in hours that the wind blew in each of the four prevailing quadrants (northwest, northeast, southwest, and southeast) or no wind (calm). Wind speed was measured as kilometers per day. Rainfall was measured as mm of daily precipitation, and river flow was calculated as the daily average volume of water in m3/s from the river Ulla. Oceanic parameters (chlorophyll, phytoplankton counts, seawater temperature, and salinity) for each sampling site were obtained from the data provided by the marine environmental-monitoring program that exists in Galicia and is carried out by the Instituto Tecnolóxico para o Control do Medio Mariño de Galicia (Xunta de Galicia Vilagarcia de Arousa, Spain). Daily values of upwelling were estimated for the point with the coordinates 10°30'W, 42°30'N from the 6-h upwelling values of the Pacific Fisheries Environmental Laboratory (National Oceanic and Atmospheric Administration, Pacific Grove, CA) according to the criteria described by the Pacific Fisheries Environmental Laboratory (http://www.pfeg.noaa.gov/products/las/docs/wind_calc_details.html#interp).
Ambient temperature, rain, wind, hours of sunshine, and rainfall data were provided by the National Weather Institute of the Ministry of Environment and were collected from weather station number 1844 located in the ria of Pontevedra (coordinates, 8°36'59"W, 42°26'24"N; altitude, 107 m). River flow data were obtained from station 544 on the river Ulla (coordinates, 4°20'15"W, 42°50'48"N) and were provided by the Galicia-Coast Network of the Department of Hydraulic Public Domain Management of the Galician Water Department of the Xunta de Galicia.
Spatial analysis.
The results of the analyses were processed with the Geographical Information System (GIS) software ArcGIS version 9.1 and the extension Spatial Analyst by ESRI (Redlands, CA). Data were analyzed by the inverse-distance weighted-interpolation algorithm with the following settings: power, 2 (Euclidean distance); 12 points; and variable search radius. The formats of the data were Shapefile (vector data) and GRID (raster data) by ESRI, whereas the vector data source was BCN200 (Base Cartografica Numerica scale, 1:200,000) from the Instituto Geografico Nacional of Spain.
Statistical analysis.
The differences in the frequencies of V. parahaemolyticus present at different sites and periods were evaluated by chi-square and Fisher's exact tests. Significant differences in the geographical distributions of environmental parameters between rias and inside the rias (north-south and east-west) were evaluated by comparison of pairs of means by analysis of variance and Bonferroni-adjusted multiple t tests. When dependent variables were not normally distributed, differences among groups were also assessed by the Kruskall-Wallis analysis of variance of ranks test.
The associations between environmental factors and the presence of V. parahaemolyticus were initially analyzed by Pearson correlation coefficients. Environmental conditions that affected the appearance of V. parahaemolyticus in the different areas and those related to the increase in V. parahaemolyticus counts were investigated separately by the use of independent regression models. Relationships between the presence of V. parahaemolyticus and each of the environmental parameters included in the study were initially surveyed by simple logistic regression analysis. Once the significant variables at an individual level were selected, a multiple logistic regression model was conducted with them. In addition, a multiple linear regression was performed with the environmental parameter values and V. parahaemolyticus counts to identify the environmental conditions that affected the abundance of V. parahaemolyticus in the area. Predicted probabilities and odds ratios (OR) were estimated by logistic regression analysis. The OR was defined as the predicted change in the odds for a unit increase in the corresponding independent variable.
V. parahaemolyticus counts were log10 transformed for statistical purposes, and nondetectable MPN values were replaced by half of the limit of detection. All statistical analyses were carried out with SPSS version 14.0.1 (SPSS Inc.), and the level of significance was set at a P value of <0.05.
|
|
|---|
![]() View larger version (43K): [in a new window] |
FIG. 2. Variations in salinity and seawater temperature throughout the period of study (A) and distribution of the presence and abundance of V. parahaemolyticus during the same period (B).
|
Although V. parahaemolyticus was detected at all the sampling sites studied, its spatial distribution showed a distinctive pattern in each of the four rias investigated (Table 1). A significantly high presence was observed in the rias of Vigo and Muros-Noia, with 26.3% and 18.7% positive samples, respectively, whereas the rias of Arousa and Pontevedra showed significantly lower incidences, with values of 9.8% and 9.2%, respectively. The highest incidences of V. parahaemolyticus were detected at sampling stations 3 (20%), 4 (24.2%), 6 (25%), 7 (45%), and 8 (26.3%). Sampling site 7 showed the significantly highest incidence and abundance among all the areas investigated, with V. parahaemolyticus present in up to 45% of the samples and a mean density of 3,232 MPN/100 g. All of these sampling sites with maximum occurrence of V. parahaemolyticus were in the rias of Vigo and Muros-Noia, and all were located on the northern side of the inner areas of the rias, close to the mouths of the rivers, where the lowest values of salinity were characteristically recorded over the period of the study (Fig. 3). The strong correspondence between the presence of V. parahaemolyticus and reduced salinity was also observed in the pattern of distribution of V. parahaemolyticus among and within the different rias. The ria of Vigo was the least saline and the coldest of the rias and displayed the highest values of incidence and abundance of V. parahaemolyticus (26.3%). By contrast, the ria of Pontevedra was significantly more saline and displayed the highest temperatures, with a significantly lower incidence of V. parahaemolyticus (9.2%) than in the other rias. Similarly, the sampling sites with the highest incidence of V. parahaemolyticus (sites 4, 6, 7, and 8) displayed significantly lower values of salinity and temperature than the other areas. Station 7, the site of maximum presence of V. parahaemolyticus, was also the least saline of all of the coastal areas investigated.
|
View this table: [in a new window] |
TABLE 1. Average incidence and abundance of V. parahaemolyticus, salinity, and seawater temperature in the different sampling stations over the period of the study
|
![]() View larger version (63K): [in a new window] |
FIG. 3. GIS maps showing the spatial distributions of V. parahaemolyticus, salinity, and seawater temperatures in the rias of Muros-Noia, Arousa, Pontevedra, and Vigo. Note the different values associated with the color scales in the maps.
|
Simple logistic regression analysis of the presence of V. parahaemolyticus and environmental variables showed significant relationships (P < 0.05) with northwest and southwest winds, salinity at deeper levels, diatoms, and upwelling, with a lag of 1 day. Salinity was the dominant factor affecting the presence of V. parahaemolyticus and showed a strong influence on the presence of V. parahaemolyticus in the zone. An increment of 1 unit (ppt) in salinity accounted for a reduction of more than 0.5-fold in the probability of detection of V. parahaemolyticus (P < 0.0001).
To identify those environmental variables that distinctively influenced the abundance of V. parahaemolyticus, associations between the explanatory variables and V. parahaemolyticus counts were evaluated statistically, first for all of the coastal areas studied and subsequently for the ria of Vigo, where the occurrence of V. parahaemolyticus was highest. Apart from the northwest and southwest winds, salinity at deeper levels, and upwelling with a lag of 1 day, the results from multiple linear regression analysis identified seawater temperature as a distinctive significant factor governing the increase in V. parahaemolyticus levels (r = 0.26; r2 = 0.06). Decrease in salinity was the dominant factor in the model and had a highly significant negative effect on the abundance of V. parahaemolyticus (P < 0.0001). Results obtained from the analysis of only data from the ria of Vigo offered a better explanatory model (r = 0.50; r2 = 0.25) and identified northwest winds, calm periods, salinity, and seawater temperature as the only significant factors that affected the abundance of V. parahaemolyticus in the ria. Salinity was the main factor in the model (P < 0.0001).
The influences of seawater temperature and salinity on the V. parahaemolyticus dynamics in the area were further evaluated by calculation of the probability of detection of V. parahaemolyticus in association with variations in salinity and seawater temperature and with both variables simultaneously (Fig. 4). The presence of V. parahaemolyticus was mainly influenced by salinity (OR = 0.561; P < 10–10). A reduction in salinity of 2 units doubled the probability of the presence of V. parahaemolyticus. By contrast, seawater temperature showed a weaker effect and contributed to lower predictability of the presence of V. parahaemolyticus (OR = 1.169; P < 0.01). The combined effects of both of these environmental variables corroborated previous observations about their influence on the temporal and spatial distribution of V. parahaemolyticus. The probability of detection increased greatly as values of salinity dropped below 34 ppt. However, seawater temperature alone exhibited a secondary effect in the model and influenced the predicted detection only if the increase in temperature was combined with reduced values of salinity.
![]() View larger version (26K): [in a new window] |
FIG. 4. Predicted probabilities for the detection of V. parahaemolyticus in association with the values of salinity (A) and seawater temperature (B) and both variables simultaneously (C).
|
Molecular characterization by PFGE of 23 selected tdh-negative environmental isolates and the two tdh-positive strains clearly discriminated the isolates according to their genetic traits (Fig. 5). The two tdh-positive isolates shared identical restriction patterns and were differentiated from the rest of the environmental isolates. They showed identical serotypes and a restriction pattern that was indistinguishable from that of the clinical strains of V. parahaemolyticus isolated in different locations in Spain over the last 10 years and belonging to the Spanish clinical clone. By contrast, all the tdh-negative environmental isolates presented high genetic diversity and showed unrelated PFGE profiles. No relationships among isolates from the same ria or among isolates obtained during the same period of study were observed.
![]() View larger version (110K): [in a new window] |
FIG. 5. Dendrogram generated by Bionumeric software showing the relationship of PFGE patterns for V. parahaemolyticus isolates included in the comparison. The numbers at the top of the figure indicate molecular sizes in kbp. <<, pandemic O3:K6 isolates.
|
|
|
|---|
The influence of seawater temperature in the seasonal pattern of V. parahaemolyticus was primarily controlled by the effect of salinity. Temperature could modulate the levels of V. parahaemolyticus only when warm waters were present during periods of reduced salinity. This pattern of the combined effects of low salinity and warm temperature has also been observed to influence the abundance of V. parahaemolyticus in formal investigations carried out in different regions of North America (11, 21). The specific dependence on salinity minimized the presence of V. parahaemolyticus during the summer, when the waters were warmest. The low occurrence and abundance of V. parahaemolyticus during the summer contrasts with previously reported patterns of presence in other areas of the world, where the highest densities were obtained during the warmest months of summer and low or undetectable levels were observed in autumn and winter (13, 18, 19, 21). The areas of study in previous investigations were mostly confined to estuaries with highly variable values of salinity, ranging from 5 ppt to 30 ppt (9, 12, 20, 21). This range of salinity is probably more suitable for growth of V. parahaemolyticus, and under such circumstances, V. parahaemolyticus populations are exposed to the influence of temperature over the entire year. Results inferred from the data of this study estimated the maximum probability of detection of V. parahaemolyticus at around a salinity of 25 ppt, similar to the levels of salinity reported to enhance the survival, growth, and cellular activity of Vibrio cholerae in microcosms (28) and also close to the optimal salinity for V. parahaemolyticus abundance in oysters of 23 ppt (13).
The dependence of V. parahaemolyticus on the lowest values of salinity was also observed in the pattern of distribution of the organism in the coastal areas of the rias. Results from GIS showed that the gradients of salinity operated as natural boundaries for the presence of V. parahaemolyticus. Each ria characteristically has a main river at its inner point that is the primary source of freshwater. The importance of this river, together with the dimensions of each ria, control the circulation of oceanic water into the ria and the levels of salinity in the estuary. Sampling station 7 in the ria of Vigo is located in the interior of the ria, an estuarine area separated from the rest of the ria by a narrow strait that promotes the retention of freshwater from the river and enables relatively high temperatures to be reached due to the reduced inflow of cold oceanic water. These characteristics convert this area into the point of maximum presence of V. parahaemolyticus among the rias.
The rias are subjected to seasonal upwelling-downwelling sequences that determine the patterns of circulation and exchange of water in the inner areas of the rias. Upwelling periods predominate in March-April and September-October, while downwelling conditions prevail the rest of the year (1). Results from this study associated the presence of V. parahaemolyticus with downwelling periods. During the upwelling season, the continental runoff has a reduced influence on the circulation and the rias show characteristics of coastal upwelling systems rather than of estuaries. By contrast, continental runoff into the rias is important during the downwelling periods, and the rias then resemble estuaries with reduced salinity. Another factor characteristic of downwelling periods is the generation of downwelling fronts through the convergence of salty oceanic water and runoff water with reduced salinity (8). The arrival of oceanic water may promote the influx of zooplankton from the ocean into the interior of the rias and its accumulation in these areas of density transition (26), as has been reported for the accumulation of toxic dinoflagellates in the interior of the rias (29). The presence of patches of zooplankton in the density fronts may lead to the accumulation of V. parahaemolyticus, due to the ability of the organism to survive in the marine environment in association with plankton (3, 17).
Molecular characterization of the isolates obtained throughout this study enabled the identification of only two pathogenic isolates of V. parahaemolyticus. The extremely low presence of pathogenic populations of V. parahaemolyticus in environmental samples is a constant characteristic in most of the investigations carried out in different regions of the world (6, 9, 11, 12, 14, 25). The two isolates identified here are, to our knowledge, the first tdh-positive strains of V. parahaemolyticus isolated directly from the marine environment of the Atlantic coast of Europe. The presence of pathogenic isolates of V. parahaemolyticus has been reported in previous investigations carried out on the coasts of France (16, 27) that identified trh-positive strains from marine sources, although no tdh-positive isolates were found during these studies. The environmental tdh-positive strains obtained from the ria of Arousa showed serotypes and genetic characteristics identical to those of isolates obtained from clinical sources during the outbreaks of infection by V. parahaemolyticus recorded in Spain in the same period (22). Most of these clinical strains were obtained from the largest outbreaks of V. parahaemolyticus illness detected in Galicia in 1999 and 2000 (22). This finding is the first evidence that identifies the marine environment of Galicia as the source of the pathogenic V. parahaemolyticus involved in the large outbreaks of 1999 and 2000.
One factor that may restrict the detection of pathogenic strains is the limitation of the analytical procedures. Analysis of shellfish samples in this study showed an overwhelming dominance of V. alginolyticus colonies over V. parahaemolyticus populations on TCBS plates. Yellow colonies of V. alginolyticus spread rapidly onto TCBS agar, almost completely covering the plates. The dominance of V. alginolyticus may have limited the detection of V. parahaemolyticus colonies, which were less frequent and grew more slowly. The consecutive streaking of two TCBS plates enhanced the isolation of Vibrio species other than V. alginolyticus on the plate with the most dilute inoculum and minimized the possibility of false-negative assays during the investigation.
The results of this study reveal the extraordinary importance of salinity in governing the seasonal pattern and the spatial distribution of V. parahaemolyticus in the marine environment of the Atlantic coast of Europe. This information provides a novel perspective regarding the ecological conditions regulating the dynamic of V. parahaemolyticus in temperate regions and is of practical use in managing shellfish-harvesting areas.
Published ahead of print on 2 November 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»