Genetic Relatedness among Environmental, Clinical, and Diseased ...

APPLIED AND ENVIRONMENTAL MICROBIOLOGY, Sept. 1998, p. 3403–3410 0099-2240/98/$04.0010 Copyright © 1998, American Society for Microbiology. All Rights Reserved.

Vol. 64, No. 9

Genetic Relatedness among Environmental, Clinical, and Diseased-Eel Vibrio vulnificus Isolates from Different Geographic Regions by Ribotyping and Randomly Amplified Polymorphic DNA PCR ´ S PUJALTE,1 ESPERANZA GARAY,1* COVADONGA R. ARIAS,1 MARIÁ JESU

AND

ROSA AZNAR1,2

Departamento de Microbiologıá, Universitat de Vale`ncia,1 and Instituto de Agroquı´mica y Tecnologıá de Alimentos, Consejo Superior de Investigaciones Cientı´ficas,2 Burjassot, E-46100 Valencia, Spain Received 9 February 1998/Accepted 23 June 1998

Genetic relationships among 132 strains of Vibrio vulnificus (clinical, environmental, and diseased-eel isolates from different geographic origins, as well as seawater and shellfish isolates from the western Mediterranean coast, including reference strains) were analyzed by random amplified polymorphic DNA (RAPD) PCR. Results were validated by ribotyping. For ribotyping, DNAs were digested with KpnI and hybridized with an oligonucleotide probe complementary to a highly conserved sequence in the 23S rRNA gene. Random amplification of DNA was performed with M13 and T3 universal primers. The comparison between ribotyping and RAPD PCR revealed an overall agreement regarding the high level of homogeneity of diseased-eel isolates in contrast to the genetic heterogeneity of Mediterranean isolates. The latter suggests the existence of autochthonous clones present in Mediterranean coastal waters. Both techniques have revealed a genetic proximity among Spanish fish farm isolates and a close relationship between four Spanish eel farm isolates and some Mediterranean isolates. Whereas the differentiation within diseased-eel isolates was only possible by ribotyping, RAPD PCR was able to differentiate phenotypically atypical isolates of V. vulnificus. On the basis of our results, RAPD PCR is proposed as a better technique than ribotyping for rapid typing in the routine analysis of new V. vulnificus isolates. Vibrio vulnificus is an autochthonous marine and estuarine bacterium which is able to colonize surfaces and internal organs of invertebrate and vertebrate marine animals (21). It may cause a rapid, invasive, and highly lethal disease associated with oyster consumption (17, 18), especially in individuals who have a preexisting chronic illness or who are immunocompromised, and may cause fatal wound infections through exposure to seawater even in people without preexisting disease (13, 22). Its isolation and its role in human infections were first reported in the United States, Japan, and Taiwan (17, 20, 32), but since 1992 reports from other geographical areas including Australia, Brazil, and Europe are increasing (4, 10, 15, 22, 25, 31, 40, 41). Infection occurs either after consumption of contaminated seafood, mainly oysters, or through wounds exposed to seawater (13, 22). V. vulnificus, being a human pathogen, includes some strains, originally defined by Tison et al. (39) as biotype 2, that constitute a serologically homogeneous group recently defined as serovar E by Biosca et al. (9) and that are responsible for severe epizootics in eel farms (3, 7). The strains pathogenic to humans are serologically and genetically diverse and have been isolated from seawater and shellfish, whereas serovar E strains have been isolated only from diseased eels or from people who handle eels (3, 40). Epidemiological reports have revealed that, in most cases, the incidence of infections by V. vulnificus is related to the consumption of raw oysters (17, 18, 27). In Spain, clinical reports of infections caused by V. vulnificus are very scarce (16, 31), in accordance with the low incidence of this species in the Mediterranean Sea (4, 5). Temperature and salinity have proven to be the two most relevant physicochemical parameters correlated with the distribution of this species (23, 32). In the Med-

iterranean, summer temperatures are always above 20°C, clearly favorable for this species. Winter temperatures are never below 10°C, and therefore the presence of viable but nonculturable cells in response to low temperatures is not expected. In contrast, salinity values (around 35‰) are much above the optimal values described for this species (23). Until very recently, the only V. vulnificus isolates at the Spanish Mediterranean coast corresponded to serovar E strains recovered from diseased eels at an eel farm with an intensive culture system placed close to the sea, although not directly connected to it, and to a few non-serovar E strains recovered from healthy eels or tank water at the same farm (3, 7, 8). This farm presently uses recirculated freshwater, but at the time of the first outbreaks (1989 to 1990) it used well water of 1.7‰ salinity. In a very recent study, and after a specific search for this species at several sea sites, a search combining culture methods for isolation and identification by PCR using specific primers (4), we were able to recover non-serovar E strains from seawater and shellfish for the first time in our coastal waters. The Mediterranean isolates together with fish farm, clinical, and environmental isolates from different geographic origins (132 strains) were subjected to ribotyping and randomly amplified polymorphic DNA (RAPD) PCR. Both techniques have been widely used for the typing of bacteria in epidemiological and ecological studies, and their advantages and disadvantages are well known. Ribotyping offers highly reproducible patterns and has already been used for intraspecific differentiation of V. vulnificus from different origins (2, 6, 9, 19, 38). We have used this technique successfully in a previous study of the fish farm isolates, together with strains from other origins (2). RAPD PCR is less laborious and time-consuming than other DNA-based techniques and seems to be the fastest genetic typing method that could be employed for a rapid identification, although it is less reproducible among different laboratories. It has also been widely used as a typing technique for both gram-positive and gram-negative bacteria (11, 26, 27, 30),

* Corresponding author. Mailing address: Departamento de Microbiologıá y Ecologıá, Universitat de Vale`ncia, Burjassot, E-46100 Valencia, Spain. Phone: 34-6-398 31 43. Fax: 34-6-386 43 72. E-mail: [email protected]. 3403

3404

ARIAS ET AL.

APPL. ENVIRON. MICROBIOL.

and more specifically to differentiate pathogenic Vibrio species, including V. vulnificus (6, 19, 42). In the present study we have used both techniques on all our serovar E and non-serovar E isolates in order to (i) analyze the genetic diversity of the Mediterranean isolates by ribotyping and RAPD PCR, by comparing the obtained profiles with the ones shown by V. vulnificus strains from different origins, (ii) determine the genetic proximity of the V. vulnificus strains isolated at the Spanish fish farm, including serovar E and non-serovar E strains, over several years, and (iii) evaluate the usefulness of RAPD as a typing technique for our isolates, by comparing the results with the ones obtained by ribotyping. MATERIALS AND METHODS Bacterial strains and growth conditions. A total of 132 strains of V. vulnificus, including diseased-eel, clinical, and environmental isolates from different origins, were used in this study. These and their sources, when known, are listed in Table 1. Strains were grown in nutrient agar (Oxoid), supplemented with 0.5% NaCl (wt/vol). Clinical strains were incubated at 37°C, and environmental isolates were incubated at 25°C. DNA isolation. Chromosomal DNA was extracted by the guanidinium thiocyanate method of Pitcher et al. (34) and further purified by RNase and proteinase K treatments (only for ribotyping) as described by Sambrook et al. (35). Ribotyping. Ribotypes were obtained as described in a preliminary study by Aznar et al. (6), where different enzymes and probes were tested to determine the combination of KpnI digests and 23S rRNA directed probe (1038) yielding the best band discrimination. In accordance with the methods of Aznar et al. (6), 5 mg of chromosomal DNA was digested with endonuclease KpnI (GIBCO BRL) as recommended by the manufacturer and the DNA restriction fragments were separated by electrophoresis in 0.8% (wt/vol) agarose gels with TAE (Trisacetate-EDTA) buffer. DNA was transferred to a noncharged nylon membrane (Qiabrane; Qiagen) under vacuum (Vacu-AID System; Hybaid). After transfer, DNA was hybridized with the 18-mer universal 1038 probe, complementary to a highly conserved region of eubacterial 23S rRNA genes that were digoxigenin labeled. Hybrid detection was performed by chemiluminescence (DIG Luminescent Detection Kit; Boehringer Mannheim). The hybridization temperature was 48°C. Band patterns displayed on X-OMAT 5 (Kodak) films were recorded with a video camera (Gel Station; Technologia para Diagno ´stico e Investigacio ń, Madrid, Spain) and stored as TIFF files. RAPD PCR. A new database, based on the RAPD PCR fingerprintings of the 131 strains included in this study, was created. Universal primers M13 (59GAA ACAGCTATGACCATG39) and T3 (59ATTAACCCTCACTAAAGG39) used in this study, as well as the PCR conditions, were described previously by Aznar et al. (6). PCR was conducted in a total volume of 50 ml containing 5 mM universal primer, 0.03 U of Taq polymerase (SuperTherm), 5 ml of Taq reaction buffer, 10 ml of 25 mM MgCl2, a 0.2 mM concentration of each deoxynucleoside triphosphate, and 100 ng of template DNA. Reaction mixtures were overlaid with 30 ml of mineral oil (Sigma, St. Louis, Mo.) and subjected to one cycle (OMNIGENE, St. Louis, Mo.) of 94°C for 5 min, 40°C for 5 min, and 72°C for 5 min. This was followed by 33 cycles of 94°C for 20 s, 48°C for 30 s, and 72°C for 45 s. The amplification products were electrophoresed on a 1.6% agarose gel, stained with ethidium bromide, and photographed under UV light. Gel images were recorded with a video camera (Gel Station; TDI) and stored as TIFF files. Banding pattern analysis. Digitized images were converted, normalized, analyzed, and combined with the software package Gel Compar, version 4.0 (Applied Maths, Kortrijk, Belgium). In order to normalize the banding patterns, molecular weight markers were included every four to six tracks. The levels of similarity between pairs of traces were computed by using the Pearson productmoment correlation coefficient for RAPD PCR and the Dice similarity coefficient (SD) for ribotyping (37). Data were clustered by using the unweighted pair group method by arithmetic averaging (UPGMA) algorithm (37).

RESULTS Ribotyping. A total of 132 V. vulnificus strains, including 80 isolates previously analyzed (2), 47 Mediterranean non-serovar-E isolates (4), and 5 new isolates recovered at the Spanish eel farm during a recent outbreak, were analyzed in the present work. For the eel farm isolates, one (A1) was recovered from a diseased eel and belonged to serovar E and four (P3 to P6) were non-serovar E isolates. In all, 34 ribotypes were distinguished among the strains analyzed. They have been des-

ignated RT 1 to RT 34 (Table 1). The Mediterranean isolates exhibited 14 ribotypes different from the ones already recorded in our database (2). The new isolates from the fish farm displayed two different ribotypes: RT 33, shown by strain A1, which corresponded to an already-described serovar E ribotype, and RT 10 (shared by the rest of strains), which was not previously recorded. Figure 1 shows the dendrogram obtained with the 34 different profiles, after the UPGMA clustering. Despite the high level of diversity observed, the inclusion of new ribotypes in the cluster analysis did not perturb the previous groupings (2). Indeed, RT 32 to RT 34 comprised all serovar E strains, irrespective of geographic origin or outbreak. None of the 19 seawater or 28 bivalve isolates belonged either to serovar E or to RT 32 to RT 34. Non-serovar E strains isolated at the Spanish eel farm between 1990 and 1997 exhibited RT 6 to RT 10. RT 10, which was displayed by strains isolated in 1997, clustered with RT 7 and RT 9 (isolated in 1990) at 82% similarity. Eight seawater isolates from a place close to the fish farm, recovered in three different samplings over one year, showed a new ribotype (RT 18), which clustered with that of four eel farm isolates (RT 8). The shellfish isolates showed five new ribotypes, RT 20, RT 26, and RT 29 to RT 31, and 22 strains shared the same ribotype (RT 20), which was the closest neighbor to serovar E ribotypes. The rest of new ribotypes were quite diverse as demonstrated by the low levels of similarity among them and could not be related to geographical origins or sources. RAPD PCR. RAPD PCR with primers M13 and T3 rendered reproducible profiles consisting of 7 to 14 bands ranging from 162 to 3,800 bp. Results obtained with both primers allowed the differentiation of the isolates at the intraspecific level. Figure 2 shows the 13 clusters defined at the 77% similarity level by RAPD patterns obtained with primer M13 and the 10 strains which remained ungrouped. Cluster 1 includes mainly Spanish eel farm and shellfish isolates and a U.S. clinical isolate (strain UMH1). Cluster 3 includes all serovar E isolates, with the exception of strains 171 and VIB 523 (both positive for indole production), which grouped apart. Cluster 4 includes Belgium eel isolates together with the type strain (ATCC 27562). Cluster 5 includes clinical and shellfish strains. Clusters 6 and 7 include Spanish seawater and eel farm isolates. Clusters 8 and 9 include eel and shellfish isolates from different origins. Cluster 10 includes Spanish shellfish and seawater isolates. Clusters 11, 12, and 13 include Spanish seawater isolates which exhibited atypical responses for some phenotypic traits. A cluster analysis of RAPD profiles obtained with primer T3 yielded 12 clusters at the 85% similarity level, and seven strains remained ungrouped. As with primer M13, all serovar E isolates except strains 171 and VIB 523 grouped together (cluster 8). The rest of the clusters included clinical and environmental isolates, and the phenotypically atypical strains formed several small clusters (Table 1). Comparison of ribotyping and RAPD PCR results. The comparison of the two fingerprinting techniques was performed by determining the location of all the strains with the same ribotype in the RAPD PCR dendrogram. For serovar E, a genetically homogeneous group, results with both techniques were in agreement. Ribotyping produced three very closely related profiles (RT 32 to RT 34), which were included in clearly defined clusters by RAPD analysis (cluster 3 with primer M13 or cluster 8 with primer T3). The only two serovar E strains that were not included in any of the two clusters corresponded to atypical biogroup 2 isolates, according to the criteria of Tison et al. (39), now serovar E according to Biosca et al. (9). These isolates shared the main serovar E ribotype (RT 32) but

VIBRIO VULNIFICUS DIVERSITY BY RIBOTYPING AND RAPD PCR

VOL. 64, 1998

3405

TABLE 1. Strains used in this study Source (date [mo-yr])

Straina T

ATCC 27562 C7184b UMH1b Vv1b 374b L180b TW1 LMG 12095 M631 D10 M76 E109 E110 M40 M47 M63 M79 M89 M90 M324 M328 M335 M626 M647 E112 E113 E114 P3 P4 P5 P6 143c LMG 12090c 147c 149c LMG 12092 155c 161c 169c 160c 167c 248 249c 170c LMG 13637 VIB 537d VIB 521d VIB 535d VIB 536d C1 C2 C32 C33 C36 C37 C38 C39 C3 C4–C25 C26 C27 C28 C29 C30 C31 C34 C35

Human wound infection Human blood Fatal wound infection Fatal wound infection Septicemia case Septicemia case Tank water fish farm (10-90) Eel Eel (9-90) Sea bream fish farm (7-90) Eel (7-90) Eel (7-90) Eel (7-90) Eel (3-90) Eel (3-90) Eel (3-90) Eel (3-90) Eel (3-90) Eel (3-90) Eel (3-90) Eel (3-90) Eel (3-90) Eel (9-90) Eel (9-90) Eel (9-90) Eel (9-90) Eel (9-90) Tank water fish farm (6-97) Tank water fish farm (6-97) Tank water fish farm (6-97) Tank water fish farm (6-97) Eel Eel Eel Eel Eel Eel Eel Eel Eel Glass eel Eel Eel Human blood Shrimp Shrimp Unknown Eel Eel Seawater (4-96) Seawater (4-96) Seawater (12-96) Seawater (12-96) Seawater (5-97) Seawater (5-97) Seawater (5-97) Seawater (5-97) Seawater (4-96) Wedge shell (5-96) Seawater (5-96) Seawater (5-96) Seawater (5-96) Seawater (5-96) Seawater (5-96) Seawater (5-96) Seawater (5-97) Seawater (5-97)

Country i

USA USA USA USA USA USA Spain Belgium Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Belgium Belgium Belgium Belgium Belgium Belgium Belgium Belgium Belgium Belgium Belgium Belgium Belgium Thailand Thailand Unknown Sweden Sweden Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain

Serovar E e

2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2

Ribotype f

1 2 3 4 5 NDj 6 6 7 8 8 8 8 9 9 9 9 9 9 9 9 9 9 9 9 9 9 10 10 10 10 11 11 11 11 11 11 11 11 12 13 13 13 14 15 15 16 17 17 18 18 18 18 18 18 18 18 19 20 21 22 22 23 23 24 25 25

RAPD PCR result g with primer: M13

T3

4 5 1 5 5 — 1 8 2 6 6 6 6 1 1 1 1 1 1 1 1 1 1 2 2 2 1 1 1 1 1 4 4 4 4 4 4 4 — 8 9 9 9 — 8 — — — — 7 7 7 7 6 6 6 6 1 1 11 11 — 13 13 7 12 12

—h 5 — 1 5 7 5 7 2 7 7 7 7 3 3 3 3 3 3 3 3 3 1 1 1 4 4 4 4 7 7 7 7 7 7 7 — 6 1 2 2 — 6 6 7 — 7 7 7 7 7 7 7 7 7 4 7 12 — — 12 12 7 9 9

Continued on following page

3406

ARIAS ET AL.

APPL. ENVIRON. MICROBIOL. TABLE 1—Continued

Straina

C40 C41 C42 C43 C44 C45 C46 C47 ATCC 33149 NCIMB 2136 NCIMB 2137 E4 E12 E22 E24 E32 E37 E39 E40 E52 E56 E58 E64 E80 E86 E92 M206 M338 E103 E105 E106 ER1 R1 R2 VIB 525d VIB 526d VIB 527d VIB 524d 171c VIB 523d NCIMB 2138 LMG 13638 A1 VIB 522d

Source (date [mo-yr])

Wedge shell (5-97) Seawater (6-97) Seawater (6-97) Common littleneck (6-97) Common littleneck (8-97) Wedge shell (8-97) Wedge shell (8-97) Wedge shell (8-97) Diseased eel Diseased eel Diseased eel Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (89) Diseased eel (90) Diseased eel (90) Diseased eel (90) Diseased eel (90) Diseased eel (90) Diseased eel (90) Diseased eel (90) Diseased eel (90) Diseased eel (90) Diseased eel (92) Diseased eel (94) Diseased eel (94) Diseased eel Diseased eel Diseased eel Diseased eel Diseased eel Diseased eel Diseased eel Shrimp Diseased eel (97) Human blood

Country

Spain Spain Spain Spain Spain Spain Spain Spain Japan Japan Japan Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Spain Sweden Sweden Sweden Norway Belgium Sweden Japan Unknown Spain Unknown

Serovar E e

2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

Ribotype f

26 27 28 ND 29 30 31 31 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 33 33 33 34

RAPD PCR result g with primer: M13

T3

— 10 10 10 10 10 5 5 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 — — 3 3 3 3

— 11 10 11 — 10 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 8 8 8 8

a ATCC, American Type Culture Collection, Manassas, Va.; NCIMB, National Collection of Industrial and Marine Bacteria, Aberdeen, Scotland; LMG, Laboratorium voor Microbiologie, Rijksuniversiteit, Ghent, Belgium. Strains in boldface correspond to the new V. vulnificus isolates (C1 to C47 are Mediterranean isolates; P3 to P6 and A1 are new eel farm isolates). b Supplied by J. D. Oliver, University of North Carolina, Charlotte. c Supplied by L. Grisez, Katholieke Universiteit, Leuven, Belgium. d Supplied by L. Verdonck, Laboratorium voor Microbiologie, Rijksuniversiteit, Ghent, Belgium. e Ascription to serovar E was performed by slide agglutination with serum anti-E39 (diseased-eel isolate). 1, serovar E; 2, not serovar E. f Ribotype designations are each preceded by RT (e.g., RT 1). g Clusters obtained by RAPD PCR after the UPGMA analysis. h —, ungrouped strain. i USA, United States. j ND, not determined.

clustered quite far from the rest of the serovar E isolates with both primers by RAPD PCR. Isolates from healthy eels or tank water at the Spanish eel farm that showed closely related ribotypes (RT 7, RT 9, and RT 10) belonged to the same cluster (cluster 1), as determined with primer M13, or to closely related ones (clusters 1 to 4), as determined with T3. RT 8, exhibited by four isolates from healthy eels at the Spanish eel farm (strains D10, M76, E109, and E110), clustered with RT 18 (main ribotype exhibited by Spanish seawater iso-

lates) (Table 1); by RAPD PCR with primer M13 these strains were included in cluster 6 together with Spanish seawater isolates belonging to RT 18. When primer T3 was used, these strains were placed in cluster 7 together with strains from other origins, including the same Spanish seawater isolates (Table 1). Shellfish strains showing RT 20, which clustered close to serovar E and non-serovar E Spanish eel farm ribotypes (Fig. 1), were also determined to be related to these strains by RAPD PCR. In general, the correspondence between results obtained


VOL. 64, 1998

FIG. 1. Dendrogram based on UPGMA cluster analysis of the 34 different ribotypes obtained in this study. The scale measures the percentage of similarity. Numbers next to the dendrogram are ribotype designations, and numbers of strains sharing the same ribotype are shown in parentheses.

by ribotyping and RAPD analysis was better when primer M13 was employed to generate the profiles. DISCUSSION In recent years several studies have been performed in order to analyze the intraspecific diversity of V. vulnificus. This characterization has been carried out by different techniques including molecular typing methods such as lipopolysaccharide typing, total protein profiles (sodium dodecyl sulfate-polyacrylamide gel electrophoresis), and nucleic acid-based methods: DNA sequencing, plasmid profiling, ribotyping, pulsedfield gel electrophoresis, RAPD PCR, and amplified-fragment length polymorphisms (AFLP). Several serotypes could be differentiated based on lipopolysaccharide or capsular antigens (9, 28, 36), although many isolates remained untypeable. Total protein profiles revealed a high level of homogeneity within the species. This is not useful for intraspecific differentiation (3). Plasmid profiles appeared to be good markers for serovar E isolates (8), but their value for the rest of the species is limited (14). Isolates related by common geographic origin or by serotype could be grouped by ribotyping (2, 9, 38) or RAPD profiles (6). Pulsed-field gel electrophoresis and AFLP allowed a finer differentiation of isolates from the same origin, revealing a high level of intraspecific diversity (2, 12). Recently, some of these typing techniques have been combined in order to obtain a better picture of the species diversity (3, 19, 38), although different protocols employed for the same technique by different authors have led to different results (6, 19). In this work we have searched for the genetic relationships among clinical, environmental, and diseased-eel isolates of V. vulnificus, including a high number of strains from an eel farm that has suffered several infections by this bacterium and seawater strains from the Spanish Mediterranean coast. We had previously used ribotyping for 80 V. vulnificus strains (2) and RAPD analysis for 21 V. vulnificus strains and 5 strains of

3407

other Vibrio species (6). Both were found to be good techniques for intraspecific differentiation, with RAPD analysis being more rapid and simpler to perform. The 52 new isolates yielded 15 new ribotypes in addition to the 19 already obtained in our previous study (2). Most of the new isolates used in the present work represent the first seawater and shellfish isolates of V. vulnificus recovered so far on the west coast of the Mediterranean (4). The rest (five strains) are recent isolates from the same eel farm where the Spanish diseased-eel isolates analyzed were recovered (2). Among the new eel farm isolates, strain A1 was recovered from a diseased eel during an outbreak that occurred in July 1997 and was assigned to serovar E by slide agglutination (unpublished data). It displays RT 33, which was previously assigned to two non-Spanish isolates and which is one of the ribotypes displayed by serovar E isolates. All strains isolated from diseased eels during epizootics that occurred from 1989 to 1994 in the Spanish eel farm shared RT 32. This ribotype is also exhibited by three of four Japanese reference strains (ATCC 33149, NCIMB 2136, and NCIMB 2137), originally isolated from diseased eels, and by diseased-eel isolates from other European countries as well. The fact that the strain isolated during the last outbreak at the Spanish farm exhibits the same ribotype as the fourth Japanese reference strain (NCIMB 2138) could be explained by the existence of a few original clones that have become widely distributed and that are responsible for outbreaks in several countries. This finding, together with the lack of reports on the isolation of serovar E strains from seawater samples, even in the vicinity of the fish farm that has suffered several outbreaks caused by this serovar, supports the hypothesis of animal-to-animal transmission (eel to eel or eel to man) and casts doubts on the hypothesis of water transmission of the disease, a hypothesis based only on a laboratory study (1). The existence of one serovar E strain isolated from human blood would support the animal-to-animal transmission hypothesis, since human infections related to the handling of eels have been reported (13, 40, 41). In addition to strain A1, some non-serovar E strains (P3 to P6) were isolated from tank water during the last outbreak at the Spanish fish farm. All of them shared RT 10, which is very similar to other eel farm ribotypes (RT 7 and RT 9) described for strains isolated 7 years ago, revealing the close relationship among fish farm isolates. As expected, diversity among seawater isolates was higher and the ribotypes obtained did not cluster, but eight strains isolated in three different samplings over a year showed the same ribotype (RT 18). Interestingly, they were recovered from a sampling site located at the seaside close to the fish farm and clustered at 90% similarity with an eel farm ribotype (RT 8). The diversity in ribotypes found in the seawater isolates suggests the existence of autochthonous clones of V. vulnificus present on the west coast of the Mediterranean. A V. vulnificus RAPD profile database has been created with the 131 strains, and the results of this technique have been compared with the ones obtained by ribotyping. In all, a good correspondence between the two typing techniques was observed. Of the two primers used, M13 allowed a more accurate definition of clusters and a better correspondence with the results obtained by ribotyping. Both primers differentiated serovar E strains as a tight group, clearly separated from the rest of strains. Furthermore, with both primers the serovar E strain recovered from the last outbreak (A1) clustered with the rest of the strains of serovar E isolated several years before with a very high level of similarity, in accordance with the results obtained by ribotyping. Nevertheless, its differentiation from the rest of the strains of serovar E was only possible through

3408

ARIAS ET AL.


FIG. 2. RAPD PCR patterns of V. vulnificus. The dendrogram was derived from an UPGMA cluster analysis of the RAPD PCR profiles of 131 V. vulnificus strains. The tracks show the processed band patterns after conversion, normalization, and subtraction of the background.

ribotyping. The close relationship between some eel farm isolates (displaying RT 8) and seawater isolates obtained at the sampling site located nearest to the fish farm (displaying RT 18) was identified as well. The genetic proximity among the Spanish fish farm isolates, irrespective of the source (eel or tank water), detected by ribotyping analysis was also detected by RAPD PCR with both primers. Phenotypically atypical strains clustered separately on the basis of RAPD analysis, whereas they were not differentiated by ribotyping. In this sense, it is noteworthy that the two atypical (indole-positive) serovar E strains (171 and VIB 523) did not cluster with the rest of strains belonging to this serovar when RAPD PCR was

used, whereas they displayed the main ribotype (RT 32). This result is in accordance with the ones obtained by AFLP fingerprinting in a previous work (2) and reveals the usefulness of RAPD analysis for intraspecific differentiation of V. vulnificus, as found in other studies (6, 42). Nevertheless, other authors did not find a good correspondence between the two techniques: RAPD PCR was found to be less discriminatory (19). In any case, reports on comparisons among techniques have to be taken with caution, and only when exactly the same protocols are employed are results comparable. The existence of large databases for the different techniques that allow the comparison of new isolates is also highly desirable.

VOL. 64, 1998


3409

FIG. 2—Continued.

The results obtained in this study have demonstrated that both ribotyping and RAPD PCR are good typing techniques for V. vulnificus at the intraspecific level; they have allowed us to differentiate among V. vulnificus strains from different types of samples and geographic origins. They have revealed a high level of genomic diversity among seawater isolates from the western Mediterranean coast; these isolates seem to constitute an autochthonous population, some members of which are closely related to environmental strains recovered from a fish farm in the same area. By contrast, the high level of homogeneity among diseased-eel isolates, irrespective of geographic origin, was confirmed by both techniques. In a previous work we had reported the usefulness of ribotyping for intraspecific differentiation of V. vulnificus biotypes, which has allowed in this study the differentiation of ribotypes for diseased-eel isolates from the same eel farm. The RAPD PCR technique has been compared with and validated by ribotyping; it yielded similar overall results and was especially useful for the differentiation of phenotypically atypical strains. On the basis of our experience, and taking into account the results of previous studies, we propose the use of RAPD PCR for rapid typing in the routine analysis of new isolates of V. vulnificus and ribotyping or AFLP for a finer discrimination within diseased-eel isolates.

ACKNOWLEDGMENTS This work has been supported by “Comisio ń Interministerial de Ciencia y Tecnologıá” grant AGF95-0264. C.R.A. is the recipient of a Ph.D. fellowship from the Ministerio de Educacio ń y Ciencia. REFERENCES 1. Amaro, C., E. G. Biosca, B. Fouz, E. Alcaide, and C. Esteve. 1995. Evidence that water transmits Vibrio vulnificus biotype 2 infections to eels. Appl. Environ. Microbiol. 61:1133–1137. 2. Arias, C. R., L. Verdonck, J. Swings, E. Garay, and R. Aznar. 1997. Intraspecific differentiation of Vibrio vulnificus biotypes by amplified fragment length polymorphism and ribotyping. Appl. Environ. Microbiol. 63:2600–2606. 3. Arias, C. R., L. Verdonck, J. Swings, R. Aznar, and E. Garay. 1997. A polyphasic approach to study the intraspecific diversity amongst Vibrio vulnificus isolates. Syst. Appl. Microbiol. 20:622–633. 4. Arias, C. R., R. Aznar, M. J. Pujalte, and E. Garay. 1998. A comparison of strategies for the detection and recovery of Vibrio vulnificus from marine samples of the Western Mediterranean coast. Syst. Appl. Microbiol. 21:128– 134. 5. Arias, C. R., M. C. Macia ń, R. Aznar, E. Garay, and M. J. Pujalte. Low incidence of V. vulnificus among Vibrio isolates from seawater and shellfish of the Western Mediterranean coast. J. Appl. Microbiol., in press. 6. Aznar, R., W. Ludwig, and K.-H. Schleifer. 1993. Ribotyping and randomly amplified polymorphic DNA analysis of Vibrio vulnificus biotypes. Syst. Appl. Microbiol. 16:303–309. 7. Biosca, E., C. Amaro, C. Esteve, E. Alcaide, and E. Garay. 1991. First record of Vibrio vulnificus biotype 2 from diseased European eel, Anguilla anguilla L. J. Fish Dis. 14:103–109. 8. Biosca, E., J. Oliver, and C. Amaro. 1996. Phenotypic characterization of

3410

9.

10. 11.

12.

13. 14. 15. 16. 17. 18. 19.

20. 21. 22.

23. 24.

25.

ARIAS ET AL.

Vibrio vulnificus biotype 2, a lipopolysaccharide-based homogeneous O serogroup within Vibrio vulnificus. Appl. Environ. Microbiol. 62:918–927. Biosca, E. G., C. Amaro, J. L. Larsen, and K. Pedersen. 1997. Phenotypic and genotypic characterization of Vibrio vulnificus: proposal for the substitution of the subspecific taxon biotype for serovar. Appl. Environ. Microbiol. 63:1460–1466. Bock, T., N. Christensen, N. H. R. Eriksen, S. Winter, H. Rygaard, and F. Jørgensen. 1994. The first fatal case of Vibrio vulnificus in Denmark. APMIS 102:874–876. Brousseau, R., A. Saint-Onge, G. Pre´fontaine, L. Masson, and J. Cabana. 1993. Arbitrary primer polymerase chain reaction, a powerful method to identify Bacillus thuringiensis serovars and strains. Appl. Environ. Microbiol. 59:114–119. Buchrieser, C., V. V. Gangar, R. L. Murpree, M. L. Tamplin, and C. W. Kaspar. 1995. Multiple Vibrio vulnificus strains in oysters as demonstrated by clamped homogeneous electric field gel electrophoresis. Appl. Environ. Microbiol. 61:1163–1168. Dalsgaard, A., N. Frimodt-Møller, B. Brunn, L. Høi, and J. L. Larsen. 1996. Clinical manifestations and molecular epidemiology of Vibrio vulnificus infections in Denmark. Eur. J. Microbiol. Infect. Dis. 15:227–232. Davidson, L. S., and J. D. Oliver. 1986. Plasmid carriage in Vibrio vulnificus and other lactose-fermenting marine vibrios. Appl. Environ. Microbiol. 51: 211–213. Ghosh, H. K., and T. E. Bowen. 1980. Halophilic vibrios from human tissue infections on the Pacific coast of Australia. Pathology 12:397–402. Herna éz, J., F. Gonza ´lez, M. Provencio, M. A. Romera, M. F. Portero, R. Pe´rez-Maestu, C. Masa-Va ´zquez, and J. Martıńez de Letona. 1991. Septicemia por Vibrio vulnificus. Rev. Esp. Microbiol. Clin. 6:144–145. Hlady, W. G., R. C. Mullen, and R. S. Hopkins. 1993. Vibrio vulnificus from raw oysters leading cause of reported deaths from foodborne illness in Florida. J. Florida. Med. Assoc. 80:2–4. Hlady, W. G., and K. C. Klontz. 1996. The epidemiology of Vibrio infections in Florida, 1991–1993. J. Infect. Dis. 173:1176–1183. Høi, L., A. Dalsgaard, J. L. Larsen, J. M. Warner, and J. D. Oliver. 1997. Comparison of ribotyping and randomly amplified polymorphic DNA PCR for characterization of Vibrio vulnificus. Appl. Environ. Microbiol. 63:1674– 1678. Hor, L. I., C. T. Goo, and L. Wan. 1995. Isolation and characterization of Vibrio vulnificus inhabiting the marine environment of the Southwestern area of Taiwan. J. Biomed. Sci. 2:384–389. Horre´, R., G. Marklein, and K. P. Schaal. 1996. Vibrio vulnificus, an emerging human pathogen. Zentbl. Bakteriol. 284:273–284. Hoyer, J., E. Engelmann, R.-M. Liehr, A. Distler, H. Hahn, and T. Shimada. 1995. Septic shock due to Vibrio vulnificus serogroup O4 wound infection acquired from the Baltic Sea. Eur. J. Clin. Microbiol. Infect. Dis. 14:1016– 1018. Kaspar, C. W., and M. L. Tamplin. 1993. Effects of temperature and salinity on the survival of Vibrio vulnificus in seawater and shellfish. Appl. Environ. Microbiol. 59:2425–2429. Klontz, K. C., S. Lieb, M. Schreiber, H. T. Janowski, L. M. Baldy, and R. A. Gunn. 1988. Syndromes of Vibrio vulnificus infections: clinical and epidemiological features in Florida cases, 1981–1987. Ann. Intern. Med. 109:318– 323. Landgraf, M., K. B. P. Leme, and M. L. Garcia-Moreno. 1996. Occurrence


26. 27.

28. 29. 30. 31. 32. 33.

34. 35. 36. 37. 38.

39. 40. 41. 42.

of emerging pathogenic Vibrio spp. in seafood consumed in Saõ Paulo city, Brazil. Rev Microbiol. 27:126–130. Lawrence, L. M., J. Harvey, and A. Gilmour. 1993. Development of a random amplification of polymorphic DNA typing method for Lysteria monocytogenes. Appl. Environ. Microbiol. 59:3117–3119. Linton, C. J., A. D. Smart, J. P. Leeming, H. Jalal, A. Telenti, T. Bodmer, and M. R. Millar. 1995. Comparison of random amplified polymorphic DNA with restriction fragment length polymorphism as epidemiological typing methods for Mycobacterium tuberculosis. J. Clin. Pathol. 48:M133–M135. Martin, S. J., and R. J. Siebeling. 1991. Identification of Vibrio vulnificus O serovars with antilipopolysaccharide monoclonal antibody. J. Clin. Microbiol. 29:1684–1688. Mascola, L., M. Torney, and D. Dassey. 1996. Vibrio vulnificus infections associated with eating raw oysters. JAMA 276:937–938. Mazurier, S., A. van de Giessen, K. Heuvelman, and K. Wearnars. 1992. RAPD analysis of Campylobacter isolates: DNA fingerprinting without the need to purify DNA. Lett. Appl. Microbiol. 14:260–262. Melhus, Å., T. Holmdahl, and I. Tjernberg. 1995. First documented case of bacteremia with Vibrio vulnificus in Sweden. Scand. J. Infect. Dis. 27:81–82. Oliver, J. D., R. A. Warner, and D. R. Cleland. 1983. Distribution of Vibrio vulnificus and other lactose-fermenting vibrios in the marine environment. Appl. Environ. Microbiol. 45:985–998. Pe´rez-Moreno, M. M. O., G. Romera, G. Pous, A. M. Jardı´, J. Zaragoza, J. I. Buj, and F. J. Pe´rez. 1996. Bacteremia por Vibrio vulnificus en paciente con u ´lcera cutańea expuesta a agua de mar. Enferm. Infecc. Microbiol. Clin. 14: 66–67. Pitcher, D. G., N. A. Saunders, and R. J. Owen. 1989. Rapid extraction of bacterial genomic DNA with guanidinium thiocyanate. Lett. Appl. Microbiol. 8:151–156. Sambrook, J., E. F. Fritsch, and T. Maniatis. 1989. Molecular cloning: a laboratory manual, 2nd ed. Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y. Simonson, J., and R. J. Siebeling. 1993. Immunogenicity of Vibrio vulnificus capsular polysaccharides and polysaccharide-protein conjugates. Infect. Immun. 61:2053–2058. Sneath, P. H. A., and R. R. Sokal. 1973. Numerical taxonomy: the principles and practice of numerical classification. W. H. Freeman, San Francisco, Calif. Tamplin, M. L., J. K. Jackson, C. Buchrieser, R. L. Murphfee, K. M. Portier, V. Gangar, L. G. Miller, and C. W. Kaspar. 1996. Pulsed-field gel electrophoresis and ribotype profiles of clinical and environmental Vibrio vulnificus isolates. Appl. Environ. Microbiol. 62:3572–3580. Tison, D. L., M. Nishibushi, J. D. Greenwood, and R. J. Seidler. 1982. Vibrio vulnificus biotype 2: new biogroup pathogenic for eels. Appl. Environ. Microbiol. 44:640–646. Veenstra, J., P. J. Rietra, C. P. Stoutenbeek, J. M. Coster, H. H. de Gier, and S. Dirks-Go. 1992. Infection by an indole-negative variant of Vibrio vulnificus transmitted by eels. J. Infect. Dis. 166:209–210. Veenstra, J., P. J. G. M. Rietra, J. M. Coster, E. Slaats, and S. Dirks-Go. 1994. Seasonal variations in the occurrence of Vibrio vulnificus along the Dutch coast. Epidemiol. Infect. 112:285–290. Wu, J.-J., L.-I. Hor, and S.-L. Shiau. 1995. Differentiation of Vibrio vulnificus strains by an arbitrarily primed polymerase chain reaction. Chin. J. Microbiol. Immunol. 28:70–78.