# Highly diverse recombining populations of Vibrio cholerae and Vibrio parahaemolyticus in French Mediterranean coastal lagoons

E. Esteves, T. Mosser, F. Aujoulat, D. Hervio-Heath, P. Monfort, E. Jumas-Bilak

► To cite this version:

E. Esteves, T. Mosser, F. Aujoulat, D. Hervio-Heath, P. Monfort, et al.. Highly diverse recombining populations of Vibrio cholerae and Vibrio parahaemolyticus in French Mediterranean coastal lagoons. Frontiers in Microbiology, 2015, 6, pp.708. 10.3389/fmicb.2015.00708 . hal-01790230

HAL Id: hal-01790230

https://hal.science/hal-01790230

Submitted on 16 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Highly diverse recombining populations of *Vibrio cholerae* and *Vibrio parahaemolyticus* in French Mediterranean coastal lagoons

Kévin Esteves<fnref n="1" />, Thomas Mosser<fnref n="1" />, Fabien Aujoulat<fnref n="1" />, Dominique Hervio-Heath², Patrick Monfort<fnref n="1" /> and Estelle Jumas-Bilak<fnref n="1" />,³*

## OPEN ACCESS

**Edited by:**
Télesphore Sime-Ngando,
Centre National de la Recherche,
France

**Reviewed by:**
Martin W. Hahn,
University of Innsbruck, Austria
Emma Jane Rochelle-Newall,
Institut de Recherche Pour le
développement, France

**\*Correspondence:**
Estelle Jumas-Bilak,
"HydroSciences", UMR 5569, Equipe
"Pathogènes Hydriques Santé
Environnement", Faculté de
Pharmacie, 15 Avenue Charles
Flahault, BP 14491,
34093 Montpellier, France
ebilak@univ-montp1.fr

**Specialty section:**
This article was submitted to
Aquatic Microbiology,
a section of the journal
Frontiers in Microbiology

**Received:** 11 February 2015
**Accepted:** 28 June 2015
**Published:** 16 July 2015

**Citation:**
Esteves K, Mosser T, Aujoulat F, Hervio-Heath D, Monfort P and Jumas-Bilak E (2015) Highly diverse recombining populations of Vibrio cholerae and Vibrio parahaemolyticus in French Mediterranean coastal lagoons. Front. Microbiol. 6:708. doi: 10.3389/fmicb.2015.00708

Vibrio parahaemolyticus and Vibrio cholerae are ubiquitous to estuarine and marine environments. These two species found in Mediterranean coastal systems can induce infections in humans. Environmental isolates of V. cholerae (n = 109) and V. parahaemolyticus (n = 89) sampled at different dates, stations and water salinities were investigated for virulence genes and by a multilocus sequence-based analysis (MLSA). V. cholerae isolates were all ctxA negative and only one isolate of V. parahaemolyticus displayed trh2 gene. Most Sequence Types (ST) corresponded to unique ST isolated at one date or one station. Frequent recombination events were detected among different pathogenic species, V. parahaemolyticus, V. cholerae, Vibrio mimicus, and Vibrio metoecus. Recombination had a major impact on the diversification of lineages. The genetic diversity assessed by the number of ST/strain was higher in low salinity condition for V. parahaemolyticus and V. cholerae whereas the frequency of recombination events in V. cholerae was lower in low salinity condition. Mediterranean coastal lagoon systems housed V. cholerae and V. parahaemolyticus with genetic diversities equivalent to the worldwide diversity described so far. The presence of STs found in human infections as well as the frequency of recombination events in environmental vibrios populations could predict a potential epidemiological risk.

**Keywords:** Vibrio, human pathogens, French coastal lagoons, phylogeny, recombination, multi-locus sequence analysis, virulence factor

## Introduction

The genus Vibrio groups gram-negative halophilic bacteria that are found in marine and estuarine environments, among which some species are responsible for human infections (Colwell et al., 1977). Particularly, V. cholerae pandemic serogroups O1 and O139 cause the choleric form of the disease. Many other O serogroups have been reported for cases of V. cholerae associated with diarrhea (Octavia et al., 2013) but the pandemic serogroups O1 and O139 are distinguished by the presence of two genetic elements, the lysogenic bacteriophage designated CTXΦ containing the genes encoding cholera toxin and the Vibrio Pathogenicity Island (VPI) (Sack et al., 2004). Most of V. cholerae isolated from the environment belong to non-toxigenic non-O1/non-O139 serotypes,

<footnote n="1">¹ Team "Pathogènes Hydriques Santé, Environnements", HydroSciences Montpellier, UMR 5569, Centre National de la Recherche, IRD, Université de Montpellier, Montpellier, France, ² Ifremer, RBE, SG2M, Laboratoire Santé, Environnement et Microbiologie, Plouzané, France, ³ Département d'Hygiène Hospitalière, Centre Hospitalier Universitaire, Montpellier, France</footnote>
which have been involved in cases of diarrhea (Dutta et al., 2013), some of them carrying various repertories of virulence genes (Octavia et al., 2013). The majority of environmental strains of *Vibrio parahaemolyticus* are innocuous but subpopulations are pathogens for humans. Potentially virulent strains are characterized by the presence of the thermostable direct (tdh) and/or tdh-related (trh1 and trh2) hemolysin genes. Some pathogenic clones encountered an epidemic success that led to the dissemination of the pandemic O3:K6 clone in many countries (Nair et al., 2007) and to the worldwide emergence of other clones such as O4:K8, O4:K11, O4:K12, and O4:K13 as shown in Peru (Gavilan et al., 2013). These vibrios, especially *V. parahaemolyticus*, are recognized worldwide as agents of gastroenteritis resulting from consumption of raw or undercooked seafood (Tantillo et al., 2004; Pruzzo et al., 2005).

In most temperate coastal areas, dynamics of vibrios are seasonal. Temperature, salinity and phytoplankton, and zooplankton populations are major factors determining the abundance of vibrios (Johnson et al., 2012). The increase of *Vibrio* densities and the spread of pathogenic species in coastal marine systems have been linked to climate anomalies and in particular, to the increase of sea surface temperature (Vezzulli et al., 2013). *V. cholerae* non-O1/non-O139 and *V. parahaemolyticus* are found in coastal lagoons of Southern France (Mediterranean) (Cantet et al., 2013). Interfacing between watersheds and the sea, these transitional aquatic ecosystems are exposed to important seasonal variations of temperature and to sudden and intense Mediterranean rainfalls that lead to flash floods bringing large volumes of freshwater into lagoons, thus reducing their salinity (Pecqueur et al., 2011).

The population structure of vibrios has been increasingly studied by the means of multilocus genetics. MLST schemes and databases are available for *V. cholerae* (Teh et al., 2011; Islam et al., 2013; Octavia et al., 2013) and *V. parahaemolyticus* (Chowdhury et al., 2004a,b; Gonzalez-Escalona et al., 2008; Yan et al., 2011). Most studies on *V. cholerae* and *V. parahaemolyticus* populations are centered on human pathogenicity and/or shellfish contamination. Consequently, research efforts are concentrated on geographic regions of cholera endemicity or vibriosis outbreaks, and on pandemic serotypes or strains harboring virulence traits. It is now obvious that the ecology of potentially pathogenic vibrios extends beyond the human body (Haley et al., 2012, 2014; Islam et al., 2013) and that these vibrios do not require human transmission for persistence (Schuster et al., 2011; Haley et al., 2012). Beside investigation of vibriosis outbreaks (Chowdhury et al., 2004a; Marin et al., 2013), several MLST studies on *Vibrio* spp. considered strain collections of very diverse geographical origins and thus gave a snapshot of the global structure of *Vibrio* populations (Octavia et al., 2013) or compared clinical and environmental populations (Teh et al., 2011; Yan et al., 2011; Theethakaew et al., 2013; Turner et al., 2013; Zhang et al., 2014).

Studies involving environmental *V. cholerae* non-O1/non-O139 in cholera-free regions (Zo et al., 2009; Keymer and Boehm, 2011; Schuster et al., 2011; Octavia et al., 2013) and *V. parahaemolyticus* (Gamble and Lovell, 2011; Ellis et al., 2012; Urmersbach et al., 2014) are more recent. Findings support

that some pathogenicity factors in vibrios may have an adaptive role in the natural environment (Haley et al., 2012; Islam et al., 2013). Furthermore, more and more human clinical infections are caused by non-toxigenic non-O1/non-O139 *V. cholerae* (Chatterjee et al., 2009; Hasan et al., 2012; Luo et al., 2013) and by *V. parahaemolyticus* isolates lacking tdh, trh genes or the type three secretion system (Chao et al., 2009; Garcia et al., 2009; Harth et al., 2009; Marin et al., 2013). From a genetic point of view, relatedness in housekeeping gene sequences between some toxigenic and non-toxigenic strains provides clues to the emergence of toxigenic strains from non-toxigenic progenitors by the acquisition of virulence genes as suggested for *V. cholerae* (Mohapatra et al., 2009; Zhou et al., 2014) and for *V. parahaemolyticus* (Chao et al., 2011; Ellingsen et al., 2013; Gavilan et al., 2013). The emergence of successful pathogenic clones from environmental populations was suspected for both species (Moore et al., 2014; Velazquez-Roman et al., 2014). All these findings suggest that, beside successful pandemic and harmful strains, environmental *V. cholerae* and *V. parahaemolyticus* species behave as non-specialist pathogens with great capacities for genetic exchange. It is therefore interesting to study populations of *V. cholerae* and *V. parahaemolyticus* sampled independently from their pathogenic behavior and/or from their virulence genes repertory.

In the present work, we focused on the genetic structure and dynamics of *V. parahaemolyticus* and *V. cholerae* populations collected in coastal lagoons in the South of France (Cantet et al., 2013) in a context of local environmental changes. The genetic diversity and the recombination rate in these two populations were considered independently from any recent history of cholera or *Vibrio* outbreaks. The effects of salinity and sampling site on genetic diversity and on intra- and interspecific recombinations were particularly studied.

# Materials and Methods

## Area of Study and Sampling

The studied lagoon system is presented in **Figure 1**. The lagoons are supplied with fresh water from four main watersheds. To the east, the water from the large watershed of the Vidourle River (826 km²) and its tributary the Vistre River (598 km²) goes into the Ponant and Vidourle lagoons before flowing to the sea via the Ponant grau and the Grau du Roi channel, respectively. Numerous streams flow directly into the Mauguio or the Prévost lagoons, forming the watersheds of Mauguio and Palavas lagoons, respectively. To the west, the Lez watershed (535 km²) is formed by the Lez River and by its main tributary the Mosson River. Fresh water from Lez watershed doesn't flow directly into the Prévost lagoons but cross them via the Lez channel. In the latter system, salinity variations during heavy rainfalls are lower than those observed for the Ponant and Vidourle lagoons that directly receive large rivers. Sampling stations were situated in lagoons (stations 1, 5, and 12), in channels connecting lagoons (stations 2, 3, 6, 8, and 11) or rivers (stations 4, and 10) to the sea, and on beaches (stations 7, and 9). Surface water (5 L) was collected in July and September 2011 at 12 stations, and in November 2011 at six stations. Temperature and salinity
<figure><img src="image_2.png" /></figure>

FIGURE 1 | Location of the sampling stations along the French Mediterranean coast (Languedoc area) in lagoons (diamond) (1, 5, and 12), in channels connecting lagoons (square) (2, 3, 6, 8, and 11) or rivers (triangle) (4 and 10) to the sea, and on beaches (circle) (7 and 9).
were recorded *in situ* using a WTW LF 196 conductimeter at the time of sampling. Samples were transported in coolers (15–18°C) to the laboratory and processed within 4 h after collection.

## Vibrio Isolation

A Most Probable Number (MPN) method was applied to detect *V. parahaemolyticus*, and *V. cholerae* in water samples according Cantet et al. (2013). Enrichment of water samples (1 L, 100 and 10 mL filtered on 0.45 µm pore size nitrocellulose membranes and volumes of 1 mL, 0.1 and 0.01 mL) was performed in alkaline peptone water (APW). After incubation at 41.5°C for 18 h, 10 µl of each APW broth were spread onto CHROMagar™ Vibrio medium (Humeau, La Chapelle-sur-Erdre, France) and incubated at 37°C for 24 h. Typical colonies of *V. cholerae* and *V. parahaemolyticus* were selected for further identification.

## Molecular Identification of *V. cholerae* and *V. parahaemolyticus* and MLST Analysis

Prior to DNA analysis strains were grown overnight in APW at 37°C with shaking. Bacterial DNA was extracted from 1 mL of culture by heat lysis. Molecular identification of *V. cholerae* and *V. parahaemolyticus* were performed as previously described by PCR using specific primers for 16S–23S rRNA Intergenic Spacer Regions (ISR) (Chun et al., 1999) and *toxR* (Cantet et al., 2013), respectively.

The *ctxA* gene in *V. cholerae* and *tdh*, *trh1*, *trh2* genes in *V. parahaemolyticus* were detected as previously described by Hoshino et al. (1998) and Cantet et al. (2013), respectively. The genes included in the MLST schemes of *V. cholerae* (*adk*, *gyrB*, *mdh*, *metE*, *pntA*, *purM*, *pyrC*) and *V. parahaemolyticus* (*dnaE*, *gyrB*, *tnaA*, *recA*, *pyrC*, *pntA*, *dtdS*) were amplified as previously described by Octavia et al. (2013) and Gonzalez-Escalona et al. (2008), respectively.

PCR products and molecular weight marker (BenchTop 100 bp DNA Ladder, Promega) were separated in 1.5% (w/v) agarose gel in 0.5X Tris Borate EDTA (TBE) buffer. Amplification products were sequenced using forward primers on an ABI 3730xl automatic sequencer (LGC Genomics, Germany). Sequencing using reverse primer was performed when the sequences obtained with forward primers were not found in databases. The ssDNA sequences were deposited in the GenBank database with accession numbers from KJ656184 to KJ657569.

## Data Analysis
### Phylogeny

Gene sequences were codon-aligned using ClustalX after translation with TRANSLATE (http://www.expasy.org). The length of sequences used for further analyses is indicated in Table 1. Phylogenetic trees were reconstructed for each individual locus and for the 7 loci concatenated nucleotide sequences. Phylogenetic reconstructions were performed on IQ Tree server (Minh et al., 2013). The best-fit model was chosen for each dataset according to Akaike information criterion and used to reconstruct Maximum likelihood (ML) trees.

Ultrafast bootstrap approximation approach (UFBoot) was used to compute the support of ML phylogenetic groups (Minh et al., 2013).

## Multi Locus Sequence Typing (MLST) Analysis

The alignment obtained for phylogenetic treeing was used for assignation of each isolate to a Sequence Type (ST) according to their allelic profiles and compared to accessible sequences on pubmlst databases (http://www.pubmlst.org). Allele profiles were also analyzed using goeBURST v3 software (http://goeburst.phyloviz.net/) to determine clonal complexes (CC) defined as sets of related strains that share at least five identical alleles on the 7 loci. The singleton (S) STs corresponded to STs differing from all the others by 3 or more of the 7 loci. The Phyloviz program employing eBURST algorithm was used to determine Minimum Spanning Tree that showed the relationships between STs (http://goeburst.phyloviz.net/). Dataset for each species are constituted of strains of this study and pubMLST isolates with complete MLST information: i.e., 228 isolates for *V. cholerae*, and 1513 isolates for *V. parahaemolyticus*.

## Population Genetics and Recombination Analyses

The number of polymorphic sites was calculated using DnaSP version 5.10 (http://www.ub.edu/dnasp/). The number of synonymous (dS) and non-synonymous (dN) substitutions per site was determined on codon-aligned sequences using SNAP software (http://www.hiv.lanl.gov/content/sequence/SNAP/SNAP.html). Program LIAN 3.5 (http://pubmlst.org/) was used to determine the *h* score of genetic diversity at each locus. The mean genetic diversity (H) was then calculated as the arithmetic mean of the *h* values for all loci. Rarefaction analysis was carried out using the freeware program aRarefactWin (http://www.uga.edu/strata/software/Software.html). Program LIAN 3.5 was also used to evaluate the recombination by calculation of the standardized index of association $I_A(sI_A)$ and to test the null hypothesis of linkage disequilibrium. The homoplasy index $\phi w$ test was done on pub MLST web site (http://www.pubmlst.org) or split tree. Recombination events were detected from aligned concatenated sequences using seven different methods with default parameters of the RDP v3.44 (Martin et al., 2010) software package: general, RDP, GENECONV, Bootscan, MAxChi, CHIMAERA, and Siscan. The r/m ratio, i.e., the relative impact of recombination and mutation in the diversification of the lineages, was calculated by LDhat (McVean et al., 2004) implemented in RDP v3.44. Recombination events were also visualized by decomposition analyses with SplitsTree 4.0 software using a distance matrix in nexus format generated from the concatenated sequence (http://www.splittree.org/).

## Statistics

Data were analyzed using R software available at http://www.R-project.org and R package FactoMineR (Lê et al., 2008). Distribution of data was determined by Shapiro-Wilk test. Wilcoxon test was used to determine differences between groups. ANOVA test was realized using R software and R package FactoMineR.
<table><caption>TABLE 1 | Sequence analysis of the seven loci studies for Vibrio parahaemolyticus and Vibrio cholerae.</caption><thead><tr><th></th><th>Locus (length in pb)</th><th>Number of alleles</th><th>GC content</th><th>Number of polymorphic sites (%)</th><th>Number of non-synonymous codons</th><th>Genetic diversity (h)</th><th>dN<sup>a</sup></th><th>dS<sup>b</sup></th><th>dN/dS</th></tr></thead><tbody><tr><td rowspan="6">Vibrio cholerae</td><td>adk (416 bp)</td><td>27</td><td>0.485</td><td>27 (6.49%)</td><td>4</td><td>0.8862</td><td>0.0037</td><td>0.0267</td><td>0.1386</td></tr><tr><td>gyrB (431 bp)</td><td>36</td><td>0.497</td><td>40 (9.28%)</td><td>1</td><td>0.9344</td><td>0.0031</td><td>0.0556</td><td>0.0558</td></tr><tr><td>mdh (421 bp)</td><td>30</td><td>0.517</td><td>33 (7.84%)</td><td>0</td><td>0.9269</td><td>-</td><td>-</td><td>-</td></tr><tr><td>metE (591 bp)</td><td>53</td><td>0.49</td><td>109 (18.44%)</td><td>16</td><td>0.9659</td><td>0.0041</td><td>0.1675</td><td>0.0245</td></tr><tr><td>pntA (431 bp)</td><td>42</td><td>0.513</td><td>38 (8.82%)</td><td>4</td><td>0.9709</td><td>0.0034</td><td>0.0425</td><td>0.0800</td></tr><tr><td>purM (476 bp)</td><td>18</td><td>0.531</td><td>62 (13.03%)</td><td>0</td><td>0.6984</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="5">Vibrio parahaemolyticus</td><td>pyrC (449 bp)</td><td>50</td><td>0.504</td><td>105 (23.39%)</td><td>18</td><td>0.9743</td><td>0.0108</td><td>0.1854</td><td>0.0583</td></tr><tr><td>dnaE (557 bp)</td><td>53</td><td>0.486</td><td>127 (22.80%)</td><td>16</td><td>0.9837</td><td>0.0045</td><td>0.1173</td><td>0.0384</td></tr><tr><td>gyrB (592 bp)</td><td>54</td><td>0.477</td><td>62 (10.47%)</td><td>4</td><td>0.9860</td><td>0.0088</td><td>0.1045</td><td>0.0842</td></tr><tr><td>tnaA (423 bp)</td><td>34</td><td>0.488</td><td>27 (6.38%)</td><td>1</td><td>0.9308</td><td>0.0030</td><td>0.0430</td><td>0.0698</td></tr><tr><td>recA (729 bp)</td><td>51</td><td>0.452</td><td>150 (20.58%)</td><td>16</td><td>0.9847</td><td>0.0034</td><td>0.1590</td><td>0.0214</td></tr><tr><td rowspan="4">Vibrio parahaemolyticus</td><td>pyrC (493 bp)</td><td>49</td><td>0.483</td><td>41 (8.32%)</td><td>7</td><td>0.9747</td><td>0.0036</td><td>0.0511</td><td>0.0705</td></tr><tr><td>pntA (430 bp)</td><td>35</td><td>0.438</td><td>29 (6.74%)</td><td>4</td><td>0.9512</td><td>0.0031</td><td>0.0389</td><td>0.0797</td></tr><tr><td>dtdS (458 bp)</td><td>49</td><td>0.501</td><td>62 (13.54%)</td><td>3</td><td>0.9831</td><td>0.0101</td><td>0.4641</td><td>0.0218</td></tr></tbody></table>

# Results

## Population Structure and Genetics of V. cholerae and V. Parahaemolyticus in Mediterranean Coastal Lagoons

During the three sampling campaigns (12 stations in July and September 2011 and 6 stations in November 2011) (Figure 1), 109 strains of non-toxigenic V. cholerae, i.e., negative in ctxA-specific PCR and 89 isolates of V. parahaemolyticus were isolated. All V. parahaemolyticus isolates were negative for tdh, trh1 and trh2 genes but one, strain 8.15, was positive for trh2 gene only (Supplementary Tables S1, S2). The two Vibrio species were detected in all sites and approximately half of the strains were isolated from water collected in channels.

The 109 V. cholerae isolates displayed high genetic variability with 78 sequence types (STs) detected by MLST, the majority being represented by unique STs found only once and in one isolate (Supplementary Table S1). Only 14 STs (18%) were represented by multiple isolates. Three of them (ST255, 259, and 242) grouped strains isolated at different dates and sites, and in different salinity conditions. Seven other STs were each represented by isolates collected at the same date but in different sites (ST200, 201, 223, 246, 265, 270, 273). The 4 STs formed by strains isolated in a same site at the same date could probably be considered as replicates of a unique strain (ST 218, 219, 221, 226) (Supplementary Table S1). None of the STs described in this study except ST2 were documented in databases.

A similarly high genetic diversity was observed for V. parahaemolyticus isolates with 68 STs, 53 STs being represented by a single isolate and 15 STs (22.1%) by 2 to 5 isolates. All STs but two, represented by multiple isolates (n = 13), grouped strains originating from different sites, the most represented being ST481, ST2000 and ST2005 (Supplementary

Table S2). Eleven of the 68 STs were previously described in databases.

The STs clustering in clonal complexes (CCs) by goeBURST confirmed the high level of genetic diversity within the two species. The population of V. cholerae displayed seven small CCs containing 2 to 3 STs (Figure 2) (Supplementary Table S1). Seven CCs, each formed by only 2 STs, were also detected in the population of V. parahaemolyticus (Figure 3) (Supplementary Table S2).

The genetic diversity (H) was 0.9082 ± 0.0369 and 0.9706 ± 0.0081 for V. cholerae and V. parahaemolyticus, respectively. The genetic diversity at individual loci (h) and the number of single nucleotide polymorphism (SNP) are given in Table 1. For the two species, the rarefaction curve for each locus did not reach a plateau suggesting that we only partially explored the genetic diversity present in the ecosystem, thus confirming the high genetic diversity of the populations studied (Supplementary Figure S1).

In the two species, the level of polymorphic sites varied according to the genes from about 6 to 23% (Table 1). The multi-locus schemes used for V. cholerae and V. parahaemolyticus shared 3 identical loci, gyrB, pntA, and pyrC. The levels of polymorphism of gyrB and pntA were similar in the two species whereas pyrC displayed a much higher polymorphism in V. cholerae than in V. parahaemolyticus. In the two species, the numbers of SNPs at the 7 loci were close, 414 (12.9%) for V. cholerae and 498 (13.5%) for V. parahaemolyticus. None of these SNPs produced premature stop codon. The frequency of non-synonymous mutations was low for each locus in both species, the highest dN/dS ratio being 0.1386 for the gene adk in V. cholerae (Table 1).

Altogether, genetic data from V. cholerae and V. parahaemolyticus multi-locus schemes showed high genetic diversity in genes subjected to low selective pressure.

<footnote><sup>a</sup>dN, non-synonymous substitutions per non-synonymous site.</footnote>

<footnote><sup>b</sup>dS, synonymous substitutions per synonymous site.</footnote>
<figure><img src="image_3.png" /></figure>

**FIGURE 2 |** Population structure of *Vibrio cholerae* based on MLST data after clustering by goeBURST software. Each Sequence Type (ST) is represented by a circle with a size proportional to the number of strains by ST. STs that differed by 2 alleles (connected by dotted lines) or one

(connected by full lines) are considered to belong to the same clonal complex (CC), which is surrounded by a dotted line. Gray ellipses group ST collected at the same date in July, September or November. The gray background of the circles representing the STs is related to the salinity at the time of sampling.

## Genetic Diversity of *V. cholerae* and *V. parahaemolyticus* from Mediterranean Lagoons Compared to Worldwide Population

We compared the multi-locus genetic data of the *Vibrio* population studied herein with those available in pubMLST databases for clinical and environmental isolates originating from various countries and collected at different times. Reconstruction of minimum-spanning trees from MLST profiles obtained in pubMLST databases and in this study showed no specificity in the distribution of the strains from Mediterranean coastal lagoons. They were equally distributed in the different clusters and globally reflected the overall diversity of *V. cholerae* and *V. parahaemolyticus* strains available in pubMLST databases (data not shown). H and ST/strain scores were similar in a collection sampled during 6 months in Mediterranean lagoons and in *Vibrio* collections from databases. The H score was 0.8836 in database vs. 0.9082 in this study for *V. cholerae* and H was 0.9536 in database vs. 0.9705 in this study for

*V. parahaemolyticus*. Therefore, the very high genetic diversity observed locally in the coastal lagoon systems for the two *Vibrio* species was similar to that observed in strains collected from many countries (five continents) over the past 60 years (data not shown).

Among the 78 STs detected in this study for *V. cholerae* population, ST2 was the only ST previously deposited in the database. It corresponded to *ctx*-strains of serotype O49 isolated in Bangladesh in 2000 in a context of sporadic digestive infections (Octavia et al., 2013). It belonged to the same CC as ST1 (serotype O49) previously isolated in the same epidemiological context. ST223 and ST238 belonged to CCs that grouped clinical strains of *V. cholerae*. ST223 was related to ST65 corresponding to O37 strains isolated in the environment in Germany in 1993 and in sporadic digestive infections in Bangladesh in 2000. ST238 was related to ST30 isolated in human in Peru in 1994.

For *V. parahaemolyticus* population, 11 STs among the 68 detected herein were previously deposited in pubMLST
<figure><img src="image_4.png" /><figcaption>FIGURE 3 | Population structure of Vibrio parahaemolyticus based on MLST data after clustering by goeBURST software.</figcaption></figure>

Each Sequence Type (ST) is represented by a circle with a size proportional to the number of strains by ST. STs that differed by one (connected by full lines) or two (connected by dotted lines)

alleles are considered to belong to the same clonal complex (CC), which is surrounded by a dotted line. Gray ellipses group ST collected at the same date in July, September or November. The gray background of the circles representing the STs is related to the salinity at the time of sampling.

database, which is larger than the *V. cholerae* database. Most of them corresponded to strains isolated from water or shellfish, worldwide (USA, China, Thailand, Vietnam and Northern Europe) from 1997 to 2012. They presented the virulence genotype *tdh-/trh-,* except for one environmental strain isolated in USA (Pacific coast) in 2006 that was *tdh+/trh+.* The ST363 that contained 2 strains isolated in sites 5 and 10 in July 2011 was composed of clinical strains (O3, *tdh+/trh+) isolated in a patient with gastro-enteritis in Thailand in 1991. The ST3, which represented 15% of the strains deposited in MLST database and corresponded to the pandemic virulent strain O3:K6, was not detected in this study.

## Genetic Diversity of *V. cholerae* and *V. parahaemolyticus* according to Stations, Times of Sampling and Salinity

The 109 isolates of *V. cholerae* were equally distributed whatever the month of sampling, the station or the salinity (Supplementary Table S3) while *V. parahaemolyticus* population was not, particularly according to salinity conditions. Indeed, *V. parahaemolyticus* was distributed in two groups, a high salinity group composed of 73 isolates and a low salinity group of 16 isolates (Supplementary Table S4).

Genetic data of strains of the two species were compared according to the month of sampling (July, September, November
2011), the sampling stations (lagoon, beach, channel between river and sea, and channel between lagoon and sea) and water salinity (high, i.e., > 20‰ and low, i.e., < 20‰) (Supplementary Table S5). For *V. cholerae*, the diversity assessed by the ST/strain score showed significant differences according to salinity, with the number of ST/strain significantly higher (*p* < 0.005) when the salinity was low. For *V. parahaemolyticus*, the ST/strain score was higher in September 2011 (0.95) than in July 2011 (0.795) but the difference was not significant (*p* = 0.41) (Supplementary Table S4). However, as observed for *V. cholerae*, the ST/strain score for *V. parahaemolyticus* was significantly higher (*p* < 0.005) when the salinity was low. Both species' genetic diversity was also significantly higher in seawater from beach stations than in samples from other sampling stations (*p* < 0.005) but this could be due to a sampling bias because of the low number of strains collected from seawater.

**Figures 2, 3** show the qualitative distribution of STs by goeBURST according to sampling dates and salinity. For the two species, clustering by reconstruction of minimum spanning tree gave the same results as goeBURST clustering (data not shown). A large majority of *V. cholerae* STs and strains was represented at only one date of sampling. Indeed, only one ST (ST259), containing five isolates, was found during the three campaigns in four different sites while two others (ST242 and ST255) were found in July and September 2011 (**Figure 2**). As expected, high salinities were recorded in July and low salinities in November whereas the high and low salinities observed in September were dependent on the sampling stations. Most STs (all but ST273) found only in September were detected either at high or low salinity whereas only a few STs found in July and September (ST242, ST255, and ST259) were detected in both salinity conditions. The STs forming CCs were mostly distributed among two or three sampling campaigns, except for CC1 that contained only isolates harvested in July. Most CCs contained STs detected in high salinity samples and closely related STs detected in low salinity samples (**Figure 2**). They also contained the same or related STs detected in different stations (Supplementary Figure S2). This suggested that slight variations in the genotypes could be related to changes in environmental conditions.

When matched with the stations, the types of station and the watersheds feeding the station (Supplementary Figure S2), the population structure and distribution of STs appeared totally independent from the sampling stations. Most STs that grouped several isolates were composed of strains isolated from different stations, different types of station and different watersheds. For instance, isolates from channels and lagoons could share the same ST and/or CC and so could isolates from channels and seawater.

We observed roughly the same results for *V. parahaemolyticus* with July and September populations seldom overlapping (**Figure 3**). Among the eight STs shared between July and September, four were detected in both high and low salinities. As observed for *V. cholerae*, all CCs but one (CC7) grouped STs detected in different sampling campaigns but mostly at similar low salinities (**Figure 3**). Again, the distribution of STs and CCs appeared totally independent from the sampling stations (Supplementary Figure S3).

## Phylogeny of *V. cholerae* and *V. parahaemolyticus*

Maximum-Likelihood phylogenies based on concatenated sequences of the seven housekeeping gene fragments of the MSLT scheme (3215 bp for *V. cholerae* and 3682 bp for *V. parahaemolyticus*) are shown in **Figures 4, 5**. Major phylogenetic clusters were congruent with the CCs determined by goeBURST with the exception of *V. cholerae* CC6 and *V. parahaemolyticus* CC4, which split into distinct parts in ML phylogeny. As observed above with goeBURST structuration, phylogeny showed only a few sporadic clusters of *V. cholerae* and *V. parahaemolyticus* composed by a few strains from a same date or a same station (**Figure 4** and Supplementary Figure S2 part B).

The *V. cholerae* ML tree (**Figure 4**) displayed 3 main clades. The main one (clade A) included the majority of isolates and STs and inside this clade, a sub-clade named *A'* displayed a higher level of genetic diversity with emerging long branches. Clade B was atypical due to an external position and long branches. In monogenic-based trees (data not shown), clade B isolates were either in the *V. cholerae* clade (*adk*, *mdh*, *pntA*, and *pyrC*-based trees) or in *V. metoecus* and *Vibrio parilis* clades (*gyrB* and *purM*-based trees) remote from *V. cholerae* or on an independent lineage (*metE*-based tree). Comparison of *metE* and *purM* gene sequences of strains N1-12 and N5-13 with sequences of databases showed a best match of only 93% with other published sequences of *Vibrio* spp. The *gyrB* gene of N1-12 and N5-13 gave the best match with strains *Vibrio* sp. W0706-81 (HM009674), L8M (HM009590), and HB0308-2A5 (HM009568) non-affiliated to the species level and detected in coastal marine water in California and Hawaii. The third member of clade B, N6-9, was generally related to N1-12 and N5-13, except for *purM* phylogeny that placed it on a long independent branch remote from all other *Vibrio* spp. *purM* sequences, and for *metE* phylogeny that placed it in the *V. cholerae* clade.

Multi-locus and mono-locus phylogenies showed that clade B could not be affiliated with certainty to *V. cholerae* in spite of the positive *V. cholerae* 16S-23S rRNA ISR PCR for isolates N1-12, N5-13, and N6-9 (**Figure 6**). These conflictual phylogenies suggested that Horizontal Gene Transfer (HGT) was involved in clade B emergence. HGT of *pyrC* gene among *V. cholerae*, *V. metoecus*, *V. mimicus*, and *V. parilis* was also suspected for most strains of the sub-clade *A'*.

Two main clades were shown in the *V. parahaemolyticus* tree (**Figure 5**). The structure of the tree appeared independent from dates and sampling stations. The long branch observed for isolate S5.61 was caused by a very atypical *dnaE* sequence, which displayed less than 84% similarity with other members of the species but 98.7% identity with *V. mimicus* ATCC 33655 *dnaE* sequence. Other best matches were *V. cholerae* strain 7449 (98.5%), and *Vibrio* sp. HB0308 and L8M (96.3%), which were previously found related to *V. cholerae* clade B strains (N5-13, N1-12) in *gyrB* phylogeny. This result indicated genetic exchanges between *Vibrio* sp. related to *V. cholerae* and *V. parahaemolyticus*. It is noteworthy that *Vibrio* sp. N5-13 and *V. parahaemolyticus* S5.61 were isolated in the same station in November and September, respectively. *V. parahaemolyticus* J8.24 and S11.38 were also placed on long branches due to very polymorphic *dtdS* and *gyrB* alleles, respectively. Comparison
<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>FIGURE 6 | Maximum-Likelihood tree based on concatenated sequences of the seven housekeeping gene fragments of the MSLT scheme indicating the relative placement of *V. cholerae*</figcaption></figure>

clade B isolates regarding *V. metoecus*, *V. mimicus*, and *V. parilis*. The scale bar indicates the number of substitutions per nucleotide position.

to databases showed that the best matches of isolate J8.24 for *dtdS* was 94% with *Vibrio campbelli* and *Vibrio harveyi* type strains. The *recA* gene of strain S11.38 displayed 98.8% identity with *V. mimicus* strain LMG 7896T and only 86% with the type strain of *V. parahaemolyticus* LMG2850T. Therefore, we described a genetic exchange occurring between two species of pathogenic vibrios, *V. parahaemolyticus* and *V. mimicus* but also between *V. parahaemolyticus* and particular genotypes related to *V. cholerae*.

## Genetic Recombination in Vibrio Populations of Mediterranean Coastal Lagoon Systems

Events of recombination by HGT showed by phylogenetic incongruences were also supported by decomposition analysis using Neighbor-Net reconstruction, which revealed interconnected networks for both *V. cholerae* (Figure 7) and *V. parahaemolyticus* (Figure 8). Recombination events occurred in the CCs, among CCs, among singletons but also between CCs and singletons (Figures 7, 8). The two main clades detected in *V. parahaemolyticus* ML phylogeny were also seen by decomposition analysis showing that recombination events occurred more frequently inside each clade than between clades

(Figure 8). Most CCs observed in the two species involved branches with recombining signals. In *V. cholerae*, major events of recombination led to the emergence of CC5, CC6 and clade B (Figure 7). These emerging branches were similar to the branches corresponding to the speciation of *V. mimicus*, *V. metoecus*, and *V. parilis* (Figure 7). These results suggested that recombination was a major force driving diversification and speciation of vibrios in Mediterranean coastal lagoon systems.

Recombination was then evaluated in *V. cholerae* and *V. parahaemolyticus* by the determination of standardized index of association (sI<sub>A</sub>), the homoplasy index φw test and by using RDP software. The sI<sub>A</sub>, expected to be null when a population is at linkage equilibrium, was significantly different from zero for both *V. cholerae* (sI<sub>A</sub> = 0.1697; *p* = 0.0002) and *V. parahaemolyticus* (sI<sub>A</sub> = 0.3307; *p* < 10<sup>3</sup>). No significant sI<sub>A</sub> difference was observed according to dates, stations, and salinities for both species (data not shown). The sI<sub>A</sub> values indicated significant linkage disequilibrium suggesting that *V. cholerae* and *V. parahaemolyticus* population conserved a clonal backbone despite frequent events of recombination. In contrast, the homoplasy index φw test, which
<figure><img src="image_8.png" /></figure>

**FIGURE 7** | Neighbor-net graph reconstructed from the concatenated sequences of the 109 strains of *V. cholerae* using Splits Tree 4.0. A network-like graph indicates recombination events. ST numbering and type of sampling site (as indicated in the legend) are shown at the branch tips. The

position of *V. metoecus*, *V. parilis*, and *V. mimicus* regarding the emerging sub-clade B is indicated in the frame B. Dotted circles indicated clonal complexes as determined by goeBURST. The scale bars indicate the number of substitutions per nucleotide position.

discriminates between recurrent mutation and recombination, found statistically significant evidence for recombination (*p* = 0.0) in both *V. cholerae* and *V. parahaemolyticus* populations.

Recombination events within these two *Vibrio* populations were further analyzed by RDP software. In *V. cholerae* population detected by at least 4 out of the 7 methods performed by RDP, 34 strains (31.2% of isolates) underwent
<figure><img src="image_9.png" /><figcaption>**FIGURE 8** | Neighbor-net graph reconstructed from the concatenated sequences of the 89 strains of *V. parahaemolyticus* using Splits Tree 4.0. A network-like graph indicates recombination events. ST numbering and</figcaption></figure>

type of sampling site (as indicated in the legend) are shown at the branch tips. Dotted circles indicated clonal complexes as determined by goeBURST. The scale bars indicate the number of substitutions per nucleotide position.

at least one recombination event (Supplementary Table S6). Recombination events involved *metE*, *pyrC*, and *pntA* genes. For *V. cholerae*, the sampling month and the salinity influenced slightly the frequency of recombining isolates, a low salinity being associated with a lower frequency of recombining isolates. In *V. parahaemolyticus* population, 52 strains (58.4%) underwent at least one recombination event (Supplementary Table S7). The gene *recA* was involved in recombination events for 56.1% of the strains. No particular environmental condition influenced the frequency of recombination.

Finally, the relative impact of recombination and mutation in the diversification of the lineages was calculated. The r/m ratio was 37 and 12 for *V. cholerae* and *V. parahaemolyticus*, respectively. These values indicate a very high impact of recombination as defined by Vos and Didelot (2009).

## Discussion

The outstanding clinical and epidemic manifestations of cholera have led to consider *V. cholerae* separately from other potentially pathogen vibrios. However, *V. cholerae* and *V. parahaemolyticus* share environmental niches and display similar lifestyles, as halophilic bacteria that occasionally infect humans. Previous data showed that both species presented a high genetic diversity with a clonal structure on a background of recombination, i.e., a classical epidemic structure (Luo et al., 2013). The ecological and population structure similarities between *V. cholerae* and *V. parahaemolyticus* prompted us to compare the two species living in the same ecosystem. Herein, we explored the genetic structure of *V. cholerae* and *V. parahaemolyticus* populations collected together in a cholera-free region without any recent history of major outbreaks of
human vibriosis. The presence of non-O1/non-O139 *V. cholerae* and *V. parahaemolyticus* had been previously demonstrated in a French coastal region (Deter et al., 2010; Cantet et al., 2013). In this study, the sampling and analysis of waters from 12 different sites in a small region of about 30 km² showed the very high degree of genetic diversity of *V. cholerae* and *V. parahaemolyticus* in Mediterranean coastal lagoon systems.

The isolates harvested in Mediterranean lagoons did not carry any major virulence factors, except for one isolate of *V. parahaemolyticus* that carried only the *trh2* gene. The frequency of virulence factors in environmental populations of *V. parahaemolyticus* (1.1%) and *V. cholerae* (0%) was lower, especially for *V. parahaemolyticus*, than that observed in waters collected from other French coastal areas (English Channel, Atlantic and Mediterranean coasts)(Hervio-Heath et al., 2002; Deter et al., 2010). Studies in cholera-free temperate or cold regions are scarce but all have underlined the absence of *ctx*+ isolates in environmental strains (Zo et al., 2009; Schuster et al., 2011; Haley et al., 2012) whereas toxigenic strains were mainly detected in outbreak contexts (Tobin-D'Angelo et al., 2008; Haley et al., 2014).

Databases comparison showed that most STs described in this study were detected for the first time, demonstrating that the *V. cholerae* and *V. parahaemolyticus* populations are extremely diverse and their genetic diversity under-reported. In this study, the distribution of STs, as well as the H and ST/strain scores, are in agreement with the overall diversity of *V. cholerae* and *V. parahaemolyticus*. Indeed, we observed locally a genetic diversity similar to or higher than that observed worldwide in the last fifty years. This has been previously shown for *V. parahaemolyticus* populations (Gonzalez-Escalona et al., 2008; Yan et al., 2011) and for *V. cholerae* populations. Zhang et al. (2014) demonstrated that non-toxigenic isolates exhibited a greater diversity than did their toxigenic counterparts, suggesting that the acquisition of virulence factors and pathogenic behavior changes the population genetics of vibrios, probably by clonal selection and success during human infectious outbreaks

The vibrios diversity appeared independent from the type of sampling site and the type of lagoon system. Moreover, we did not detect any clones specialized locally according to a particular site in the lagoon system, but this could be due to the hydrodynamic conditions in these systems where water is extensively mixed. Some STs detected in Mediterranean lagoon systems were previously detected elsewhere in the world, in very different environmental conditions or in clinical samples. This suggests the absence of geographic or lifestyle specialization of *V. cholerae* and *V. parahaemolyticus* populations collected in cholera-free or vibriosis outbreak-free regions. Moreover, some STs without virulence factors identified in this study had been previously detected as virulent strains, suggesting that the population described here could provide the genetic background for the emergence of pathogenic clones (Li et al., 2014).

To our knowledge, the study of Ellis et al. (2012) and the present study are the only ones considering the genetic diversity of *V. parahaemolyticus* according to local variation of environmental conditions. The level of genetic diversity of *V. cholerae* and *V. parahaemolyticus* was higher when the salinity was low. Some of *V. parahaemolyticus* CCs detected in both low and high salinity conditions are composed of closely related STs, each being detected in only one condition of salinity. Similarly, *V. cholerae* CCs detected in different sites are formed by closely related STs. These slight differences in nucleotide sequences suggest microevolution events spurred by adaptation to changing environmental conditions. In New Hampshire shellfish waters, the diversity in *V. parahaemolyticus* populations increased as water temperature increased (Ellis et al., 2012). The increased diversity associated with warming waters and low salinity associated with precipitations suggests again that new pathogenic lineages may emerge with climate change.

It is well-known that the emergence of pathogenic clones in vibrios is linked to HGT of virulence genes or genetic elements among strains as for phage CTX in *V. cholerae* (Hazen et al., 2010). Several factors have been shown to influence the recombination rate in vibrios. In *V. cholerae*, high chitin concentrations induce an increase of recombination rate (Meibom et al., 2005). Moreover, SOS system activation by environmental factors stimulates integrase activity that in turn increases the frequency of transformation (Guerin et al., 2009). We studied the recombination rate according to the site and salinity, and showed that recombination events were more frequent in high-salinity conditions for *V. cholerae*. This could be related to the positive association between salinity and the rate of HGT detected by genomics (Jain et al., 2003). More generally, we describe a high rate of recombination in *V. parahaemolyticus* and *V. cholerae* populations. For some of the recombining STs, the donor STs were detected in the Mediterranean lagoon ecosystems. This demonstrates that, beside SNPs described above as related to changes in environmental conditions, recombination is a diversity-generating mechanism ongoing in Mediterranean lagoons. This has been previously described for environmental strains but not for clinical strains of *V. parahaemolyticus* (Chowdhury et al., 2004a; Ellis et al., 2012), which is consistent with an epidemic spread of a subset of virulent clones. Similarly, previous studies on disease-causing strains have shown *V. cholerae* to be a primarily clonal bacterial species, but the vast diversity in environmental strains showed a recombination rate of 6.5 times the mutation rate (Keymer and Boehm, 2011). One could hypothesize that environmental *Vibrio* spp. populations form a reservoir where the diversity is efficiently enhanced by recombination in order to produce a generalist population from which specialized clones can emerge. Moreover, we showed that HGT of housekeeping genes occurred among species of potentially pathogenic vibrios in Mediterranean coastal systems. In multilocus phylogeny, such HGT leads to emerging branches similar to those supporting newly described vibrios species, *V. mimicus* (Davis et al., 1981), and *V. metoecus*
(Kirchberger et al., 2014). This suggests that recombination events such as the ones we describe in Mediterranean coastal systems may drive speciation among potentially pathogenic vibrios.

## Conclusive Remarks

Mediterranean coastal lagoons host *V. cholerae* and *V. parahaemolyticus* showing a genetic diversity equivalent to the worldwide diversity described so far. The current very low frequency of virulence genes in these populations indicates a low risk for human infection. But the presence of STs involved in human infections elsewhere, as well as the frequent recombinations observed in these transitional aquatic ecosystems, indicates that autochthonous *V. cholerae* and *V. parahaemolyticus* populations could be involved in epidemiological cycles and infectious hazard.

## References

Cantet, F., Hervio-Heath, D., Caro, A., Le Mennec, C., Monteil, C., Quemere, C., et al. (2013). Quantification of *Vibrio parahaemolyticus*, *Vibrio vulnificus* and *Vibrio cholerae* in French Mediterranean coastal lagoons. *Res. Microbiol.* **164**, 867-874. doi: 10.1016/j.resmic.2013.06.005

Chao, G., Jiao, X., Zhou, X., Yang, Z., Huang, J., Pan, Z., et al. (2009). Serodiversity, pandemic O3:K6 clone, molecular typing, and antibiotic susceptibility of foodborne and clinical *Vibrio parahaemolyticus* isolates in Jiangsu, China. *Foodborne Pathog. Dis.* **6**, 1021-1028. doi: 10.1089/fpd. 2009.0295

Chao, G. X., Wang, F., Zhou, X. H., Jiao, X. N., Huang, J. L., Pan, Z. M., et al. (2011). Origin of *Vibrio parahaemolyticus* O3:K6 pandemic clone. *Int. J. Food Microbiol.* **145**, 459-463. doi: 10.1016/j.ijfoodmicro.2011. 01.022

Chatterjee, S., Patra, T., Ghosh, K., Raychoudhuri, A., Pazhani, G. P., Das, M., et al. (2009). *Vibrio cholerae* O1 clinical strains isolated in 1992 in Kolkata with progenitor traits of the 2004 Mozambique variant. *J. Med. Microbiol.* **58**, 239-247. doi: 10.1099/jmm.0.003780-0

Chowdhury, A., Ishibashi, M., Thiem, V., Tuyet, D., Van Tung, T., Chien, B., et al. (2004a). Emergence and serovar transition of *Vibrio parahaemolyticus* pandemic strains isolated during a diarrhea outbreak in Vietnam between 1997 and 1999. *Microbiol. Immunol.* **48**, 319-327. doi: 10.1111/j.1348- 0421.2004.tb03513.x

Chowdhury, N. R., Stine, O. C., Morris, J. G., and Nair, G. B. (2004b). Assessment of evolution of pandemic *Vibrio parahaemolyticus* by multilocus sequence typing. *J. Clin. Microbiol.* **42**, 1280-1282. doi: 10.1128/JCM.42.3.1280-1282.2004

Chun, J., Huq, A., and Colwell, R. R. (1999). Analysis of 16S-23S rRNA intergenic spacer regions of *Vibrio cholerae* and *Vibrio mimicus*. *Appl. Environ. Microbiol.* **65**, 2202-2208.

Colwell, R. R., Kaper, J., and Joseph, S. W. (1977). *Vibrio cholerae*, *Vibrio parahaemolyticus* and other vibrios: occurrence and distribution in Chesapeake Bay. *Science* **198**, 394-396. doi: 10.1126/science.198.4315.394-a

Davis, B. R., Fanning, G. R., Madden, J. M., Steigerwalt, A. G., Bradford, H. B., Smith, H. L., et al. (1981). Characterization of biochemically atypical *Vibrio cholerae* strains and designation of a new pathogenic species, *Vibrio mimicus*. *J. Clin. Microbiol.* **14**, 631-639.

Deter, J., Lozach, S., Derrien, A., Véron, A., Chollet, J., and Hervio-Heath, D. (2010). Chlorophyll a might structure a community of potentially pathogenic culturable *Vibrioaceae*. Insights from a one-year study of water and mussels surveyed on the French Atlantic coast. *Environ. Microbiol. Rep.* **2**, 185-191. doi: 10.1111/j.1758-2229.2010.00133.x

Dutta, D., Chowdhury, G., Pazhani, G. P., Guin, S., Dutta, S., Ghosh, S., et al. (2013). *Vibrio cholerae* non-O1, non-O139 serogroups and

## Acknowledgments

This work was supported by funding provided by the GIS "Climat Environnement Société", by the Observatoire Homme Milieux (OHM) "Littoral méditerranéen", by the Programme National EC2CO "Ecosphère Continentale et Côteière", and by the association ADEREMPHA (Sauzet, France), and a doctoral fellowship (KE.) from CNRS and Université Montpellier 2. The authors are thankful to Mylène Toubiana and Anne Lotthé for their involvement in the critical reading of the manuscript.

## Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb. 2015.00708

cholera-like diarrhea, Kolkata, India. *Emerg. Infect. Dis.* **19**, 464-467. doi: 10.3201/eid1903.121156

Ellingsen, A. B., Olsen, J. S., Granum, P. E., Rorvik, L. M., and Gonzalez-Escalona, N. (2013). Genetic characterization of trh positive *Vibrio* spp. isolated from Norway. *Front. Cell. Infect. Microbiol.* **3**:107. doi: 10.3389/fcimb.2013. 00107

Ellis, C. N., Schuster, B. M., Striplin, M. J., Jones, S. H., Whistler, C. A., and Cooper, V. S. (2012). Influence of seasonality on the genetic diversity of *Vibrio parahaemolyticus* in New Hampshire shellfish waters as determined by Multilocus sequence analysis. *Appl. Environ. Microbiol.* **78**, 3778-3782. doi: 10.1128/AEM.07794-11

Gamble, M. D., and Lovell, C. R. (2011). Infaunal burrows are enrichment zones for *Vibrio parahaemolyticus*. *Appl. Environ. Microbiol.* **77**, 3703-3714. doi: 10.1128/AEM.02897-10

Garcia, K., Torres, R., Uribe, P., Hernandez, C., Rioseco, M. L., Romero, J., et al. (2009). Dynamics of clinical and environmental *Vibrio parahaemolyticus* strains during seafood-related summer diarrhea outbreaks in southern Chile. *Appl. Environ. Microbiol.* **75**, 7482-7487. doi: 10.1128/AEM. 01662-09

Gavilan, R. G., Zamudio, M. L., and Martinez-Urtaza, J. (2013). Molecular epidemiology and genetic variation of pathogenic *Vibrio parahaemolyticus* in Peru. *PLoS Negl. Trop. Dis.* **7:e2210**. doi: 10.1371/journal.pntd. 0002210

Gonzalez-Escalona, N., Martinez-Urtaza, J., Romero, J., Espejo, R. T., Jaykus, L. A., and DePaola, A. (2008). Determination of molecular phylogenetics of *Vibrio parahaemolyticus* strains by multilocus sequence typing. *J. Bacteriol.* **190**, 2831-2840. doi: 10.1128/JB.01808-07

Guerin, E., Cambray, G., Sanchez-Alberola, N., Campoy, S., Erill, I., Da Re, S., et al. (2009). The SOS response controls integron recombination. *Science* **324**, 1034. doi: 10.1126/science.1172914

Haley, B. J., Chen, A., Grim, C. J., Clark, P., Diaz, C. M., Taviani, E., et al. (2012). *Vibrio cholerae* in a historically cholera-free country. *Environ. Microbiol. Rep.* **4**, 381-389. doi: 10.1111/j.1758-2229.2012. 00332.x

Haley, B. J., Kokashvili, T., Tskshvediani, A., Janelidze, N., Mitaishvili, N., Grim, C. J., et al. (2014). Molecular diversity and predictability of *Vibrio parahaemolyticus* along the Georgian coastal zone of the Black Sea. *Front. Microbiol.* **5**:45. doi: 10.3389/fmicb.2014.00045

Harth, E., Matsuda, L., Hernandez, C., Rioseco, M. L., Romero, J., Gonzalez-Escalona, N., et al. (2009). Epidemiology of *Vibrio parahaemolyticus* outbreaks, Southern Chile. *Emerg. Infect. Dis.* **15**, 163-168. doi: 10.3201/eid1502. 071269

Hasan, N. A., Choi, S. Y., Eppinger, M., Clark, P. W., Chen, A., Alam, M., et al. (2012). Genomic diversity of 2010 Haitian cholera outbreak strains.
Proc. Nat. Acad. Sci. U.S.A. 109, E2010–E2017. doi: 10.1073/pnas.12073
59109

Hazen, T. H., Pan, L., Gu, J. D., and Sobecky, P. A. (2010). The contribution of mobile genetic elements to the evolution and ecology of vibrios. *FEMS Microbiol. Ecol.* 74, 485–499. doi: 10.1111/j.1574-6941.2010.00937.x

Hervio-Heath, D., Colwell, R. R., Derrien, A., Robert-Pillot, A., Fournier, J. M., and Pommepuy, M. (2002). Occurrence of pathogenic vibrios in coastal areas of France. *J. Appl. Microbiol.* 92, 1123–1135. doi: 10.1046/j.1365-2672.2002.01663.x

Hoshino, K., Yamasaki, S., Mukhopadhyay, A. K., Chakraborty, S., Basu, A., Bhattacharya, S. K., et al. (1998). Development and evaluation of a multiplex PCR assay for rapid detection of toxigenic Vibrio cholerae O1 and O139. *FEMS Immunol. Med. Microbiol.* 20, 201–207. doi: 10.1111/j.1574-695X.1998.tb01128.x

Islam, A., Labbate, M., Djordjevic, S. P., Alam, M., Darling, A., Melvold, J., et al. (2013). Indigenous Vibrio cholerae strains from a non-endemic region are pathogenic. *Open Biol.* 3:120181. doi: 10.1098/rsob.120181

Jain, R., Rivera, M. C., Moore, J. E., and Lake, J. A. (2003). Horizontal gene transfer accelerates genome innovation and evolution. *Mol. Biol. Evol.* 20, 1598–1602. doi: 10.1093/molbev/mmsg154

Johnson, C. N., Bowers, J. C., Griffitt, K. J., Molina, V., Clostio, R. W., Pei, S., et al. (2012). Ecology of Vibrio parahaemolyticus and Vibrio vulnificus in the coastal and estuarine waters of Louisiana, Maryland, Mississippi, and Washington (United States). *Appl. Environ. Microbiol.* 78, 7249–7257. doi: 10.1128/AEM.01296-12

Keymer, D. P., and Boehm, A. B. (2011). Recombination shapes the structure of an environmental Vibrio cholerae population. *Appl. Environ. Microbiol.* 77, 537–544. doi: 10.1128/AEM.02062-10

Kirchberger, P. C., Turnsek, M., Hunt, D. E., Haley, B. J., Colwell, R. R., Polz, M. F., et al. (2014). Vibrio metoeocsp. nov., a close relative of Vibrio cholerae isolated from coastal brackish ponds and clinical specimens. *Int. J. Syst. Evol. Microbiol.* 64, 3208–3214. doi: 10.1099/ijs.0.060145-0

Lê, S., Josse, J., and Husson, F. (2008). FactoMineR: an R package for multivariate analysis. *J. Stat. Softw.* 25, 1–18.

Li, F., Du, P., Li, B., Ke, C., Chen, A., Chen, J., et al. (2014). Distribution of virulence-associated genes and genetic relationships in non-O1/O139 Vibrio cholerae aquatic isolates from China. *Appl. Environ. Microbiol.* 80, 4987–4992. doi: 10.1128/AEM.01021-14

Luo, Y., Ye, J., Jin, D., Ding, G., Zhang, Z., Mei, L., et al. (2013). Molecular analysis of non-O1/non-O139 Vibrio cholerae isolated from hospitalised patients in China. *BMC Microbiol.* 13:52. doi: 10.1186/1471-2180-13-52

Marin, M. A., Thompson, C. C., Freitas, F. S., Fonseca, E. L., Aboderin, A. O., Zailani, S. B., et al. (2013). Cholera outbreaks in Nigeria are associated with multidrug resistant atypical El Tor and non-O1/non-O139 Vibrio cholerae. *PLoS Negl. Trop. Dis.* 7:e2049. doi: 10.1371/journal.pntd.0002049

Martin, D. P., Lemey, P., Lott, M., Moulton, V., Posada, D., and Lefeuvre, P. (2010). RDP3: a flexible and fast computer program for analyzing recombination. *Bioinformatics* 26, 2462–2463. doi: 10.1093/bioinformatics/btq467

McVean, G. A., Myers, S. R., Hunt, S., Deloukas, P., Bentley, D. R., and Donnelly, P. (2004). The fine-scale structure of recombination rate variation in the human genome. *Science* 304, 581–584. doi: 10.1126/science.1092500

Meibom, K. L., Blokesch, M., Dolganov, N. A., Wu, C. Y., and Schoolnik, G. K. (2005). Chitin induces natural competence in Vibrio cholerae. *Science* 310, 1824–1827. doi: 10.1126/science.1120096

Minh, B. Q., Nguyen, M. A., and von Haeseler, A. (2013). Ultrafast approximation for phylogenetic bootstrap. *Mol. Biol. Evol.* 30, 1188–1195. doi: 10.1093/molbev/mst024

Mohapatra, S. S., Ramachandran, D., Mantri, C. K., Colwell, R. R., and Singh, D. V. (2009). Determination of relationships among non-toxigenic Vibrio cholerae O1 biotype El Tor strains from housekeeping gene sequences and ribotype patterns. *Res. Microbiol.* 160, 57–62. doi: 10.1016/j.resmic.2008.10.008

Moore, S., Thomson, N., Mutreja, A., and Piarroux, R. (2014). Widespread epidemic cholera caused by a restricted subset of Vibrio cholerae clones. *Clin. Microbiol. Infect.* 20, 373–379. doi: 10.1111/1469-0691.12610

Nair, G. B., Ramamurthy, T., Bhattacharya, S. K., Dutta, B., Takeda, Y., and Sack, D. A. (2007). Global dissemination of Vibrio parahaemolyticus serotype O3:K6 and its serovariants. *Clin. Microbiol. Rev.* 20, 39–48. doi: 10.1128/CMR.00025-06

Octavia, S., Salim, A., Kurniawan, J., Lam, C., Leung, Q., Ahsan, S., et al. (2013). Population structure and evolution of non-O1/non-O139 Vibrio cholerae by multilocus sequence typing. *PLoS ONE* 8:e65342. doi: 10.1371/journal.pone.0065342

Pecqueur, D., Vidussi, F., Fouilland, E., Floc'h, E. L., Mas, S., Roques, C., et al. (2011). Dynamics of microbial planktonic food web components during a river flash flood in a Mediterranean coastal lagoon. *Hydrobiologia* 673, 13–27. doi: 10.1007/s10750-011-0745-x

Pruzzo, C., Huq, A., Colwell, R. R., and Donelli, G. (2005). "Pathogenic Vibrio species in the marine and estuarine environment," in *Oceans and Health: Pathogens in the Marine Environment*, eds S. Belkin and R. R. Colwell R. R. (New York, NY: Springer), 217–252.

Sack, D. A., Sack, R. B., Nair, G. B., and Siddique, A. K. (2004). Cholera. *Lancet* 363, 223–233. doi: 10.1016/S0140-6736(03)15328-7

Schuster, B. M., Tyzik, A. L., Donner, R. A., Striplin, M. J., Almagro-Moreno, S., Jones, S. H., et al. (2011). Ecology and genetic structure of a northern temperate Vibrio cholerae population related to toxigenic isolates. *Appl. Environ. Microbiol.* 77, 7568–7575. doi: 10.1128/AEM.00378-11

Tantillo, G. M., Fontanarosa, M., Di Pinto, A., and Musti, M. (2004). Updated perspectives on emerging vibrios associated with human infections. *Lett. Appl. Microbiol.* 39, 117–126 doi: 10.1111/j.1472-765x.2004.01568.x

Teh, C. S., Chua, K. H., and Thong, K. L. (2011). Genetic variation analysis of Vibrio cholerae using multilocus sequencing typing and multivirulence locus sequencing typing. *Infect. Genet. Evol.* 11, 1121–1128. doi: 10.1016/j.meegid.2011.04.005

Theethakaew, C., Feil, E. J., Castillo-Ramirez, S., Aanensen, D. M., Suthienkul, O., Neil, D. M., et al. (2013). Genetic relationships of Vibrio parahaemolyticus isolates from clinical, human carrier, and environmental sources in Thailand, determined by multilocus sequence analysis. *Appl. Environ. Microbiol.* 79, 2358–2370. doi: 10.1128/AEM.03067-12

Tobin-D'Angelo, M., Smith, A. R., Bulens, S. N., Thomas, S., Hodel, M., Izumiya, H., et al. (2008). Severe diarrhea caused by cholera toxin-producing Vibrio cholerae serogroup O75 infections acquired in the southeastern United States. *Clin. Infect. Dis.* 47, 1035–1040. doi: 10.1086/591973

Turner, J. W., Paranjpye, R. N., Landis, E. D., Biryukov, S. V., Gonzalez-Escalona, N., Nilsson, W. B., et al. (2013). Population structure of clinical and environmental Vibrio parahaemolyticus from the Pacific Northwest Coast of the United States. *PLoS ONE* 8:e55726. doi: 10.1371/journal.pone.0055726

Urmersbach, S., Alter, T., Koralage, M. S., Sperling, L., Gerdts, G., and Messelhausser, et al. (2014). Population analysis of Vibrio parahaemolyticus originating from different geographical regions demonstrates a high genetic diversity. *BMC Microbiol.* 14:59. doi: 10.1186/1471-21-80-14-59

Velazquez-Roman, J., León-Sicairos, N., de Jesus Hernández-Diaz, L., and Canizalez-Roman, A. (2014). Pandemic Vibrio parahaemolyticus O3:K6 from the American continent. *Front. Cell. Infect. Microbiol.* 2:110. doi: 10.3389/fcimb.2013.00110

Vezzulli, L., Colwell, R. R., and Pruzzo, C. (2013). Ocean warming and spread of pathogenic vibrios in the aquatic environment. *Microbial. Ecol.* 65, 817–825. doi: 10.1007/s00248-012-0163-2

Vos, M., and Didelot, X. (2009). A comparison of homologous recombination rates in bacteria and archaea. *ISME J.* 3, 199–208. doi: 10.1038/ismej.2008.93

Yan, Y., Cui, Y., Han, H., Xiao, X., Wong, H. C., Tan, Y., et al. (2011). Extended MLST-based population genetics and phylogeny of Vibrio parahaemolyticus with high levels of recombination. *Int. J. Food Microbiol.* 145, 106–112. doi: 10.1016/j.ijfoodmicro.2010.11.038

Zhang, P., Zhou, H., Diao, B., Li, F., Du, P., Li, J., et al. (2014). A molecular surveillance reveals the prevalence of Vibrio cholerae O139 isolates in China from 1993 to 2012. *J. Clin. Microbiol.* 52, 1146–1152. doi: 10.1128/JCM.03354-13
Zhou, H., Zhao, X., Wu, R., Cui, Z., Diao, B., Li, J., et al. (2014). Population structural analysis of O1 El Tor Vibrio cholerae isolated in China among the seventh cholera pandemic on the basis of multilocus sequence typing and virulence gene profiles. *Infect. Genet. Evol.* 22, 72–80. doi: 10.1016/j.meegid.2013.12.016

Zo, Y. G., Chokesajjawatee, N., Grim, C., Arakawa, E., Watanabe, H., and Colwell, R. R. (2009). Diversity and seasonality of bioluminescent Vibrio cholerae populations in Chesapeake Bay. *Appl. Environ. Microbiol.* 75, 135–146. doi: 10.1128/AEM.02894-07

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2015 Esteves, Mosser, Aujoulat, Hervio-Heath, Monfort and Jumas-Bilak. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.