# Genome-driven evaluation and redesign of PCR tools for improving the detection of virulence-associated genes in aeromonads

Emilie Talagrand-Reboul, Fadua Latif-Eugenín, Roxana Beaz-Hidalgo, Sophie Colston, Maria-Jose Figueras, Joerg Graf, Estelle Jumas-Bilak, Brigitte Lamy

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

## To cite this version:

Emilie Talagrand-Reboul, Fadua Latif-Eugenín, Roxana Beaz-Hidalgo, Sophie Colston, Maria-Jose Figueras, et al.. Genome-driven evaluation and redesign of PCR tools for improving the detection of virulence-associated genes in aeromonads. PLoS ONE, 2018, 13 (8), pp.e0201428. 10.1371/journal.pone.0201428. hal-02094315

HAL Id: hal-02094315

https://hal.umontpellier.fr/hal-02094315

Submitted on 15 Feb 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Genome-driven evaluation and redesign of PCR tools for improving the detection of virulence-associated genes in aeromonads

Emilie Talagrand-Reboul<sup>1,2*</sup>, Fadua Latif-Eugenín<sup>3</sup>, Roxana Beaz-Hidalgo<sup>3</sup>, Sophie Colston<sup>4</sup>, Maria-Jose Figueras<sup>3</sup>, Joerg Graf<sup>4</sup>, Estelle Jumas-Bilak<sup>1,5</sup>, Brigitte Lamy<sup>1,6*</sup>

1 Équipe Pathogènes Hydriques Santé Environnements, UMR 5569 HSM, Université de Montpellier, Montpellier, France, 2 Laboratoire de Bactériologie, Hôpitaux universitaires de Strasbourg, Strasbourg, France, 3 Unidad de Microbiología, Departamento de Ciencias Médicas Básicas, Facultad de Medicina y Ciencias de la Salud, IISPV, Universidad Rovira i Virgili, Reus, Spain, 4 Department of Molecular and Cell Biology, University of Connecticut, Storrs, Connecticut, United States of America, 5 Département d'Hygiène Hospitalière, CHRU de Montpellier, Montpellier, France, 6 Département de Bactériologie, CHU de Nice, Nice, France

* brigitte_lamy@yahoo.fr (BL); emilie.talagrand@wanadoo.fr (ETR)

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

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

OPEN ACCESS

**Citation:** Talagrand-Reboul E, Latif-Eugenín F, Beaz-Hidalgo R, Colston S, Figueras M-J, Graf J, et al. (2018) Genome-driven evaluation and redesign of PCR tools for improving the detection of virulence-associated genes in aeromonads. PLoS ONE 13(8): e0201428. https://doi.org/10.1371/journal.pone.0201428

**Editor:** Daniela Flavia Hozbor, Universidad Nacional de la Plata, ARGENTINA

**Received:** February 26, 2018

**Accepted:** July 16, 2018

**Published:** August 15, 2018

**Copyright:** © 2018 Talagrand-Reboul et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by the Association des Biologistes de l'Ouest (ABO) and the Association pour la recherche et le développement en microbiologie et pharmacie (ADEREMPHA) to ETR. Part of this research was supported by USDA and ARS agreement 58-1930-4-002 to JG.

## Abstract

Many virulence factors have been described for opportunistic pathogens within the genus *Aeromonas*. Polymerase Chain Reactions (PCRs) are commonly used in population studies of aeromonads to detect virulence-associated genes in order to better understand the epidemiology and emergence of *Aeromonas* from the environment to host, but their performances have never been thoroughly evaluated. We aimed to determine diagnostic sensitivity and specificity of PCR assays for the detection of virulence-associated genes in a collection of *Aeromonas* isolates representative for the genetic diversity in the genus. Thirty-nine *Aeromonas* strains belonging to 27 recognized species were screened by published PCR assays for virulence-associated genes (*act*, *aerA*, *aexT*, *alt*, *ascFG*, *ascV*, *ast*, *lafA*, *lip*, *ser*, *stx1*, *stx2A*). In parallel, homologues of the 12 putative virulence genes were searched from the genomes of the 39 strains. Of the 12 published PCR assays for virulence factors, the comparison of PCR results and genome analysis estimated diagnostic sensitivities ranging from 34% to 100% and diagnostic specificities ranged from 71% to 100% depending upon the gene. To improve the detection of virulence-associated genes in aeromonads, we have designed new primer pairs for *aerA/act*, *ser*, *lafA*, *ascFG* and *ascV*, which showed excellent diagnostic sensitivity and specificity. Altogether, the analysis of high quality genomic data, which are more and more easy to obtain, provides significant improvements in the genetic detection of virulence factors in bacterial strains.

## Introduction

*Aeromonas* are common inhabitants of aquatic environments and can be involved in fish and human diseases. They are frequently found in drinking water and food, including meat, fish,
**Competing interests:** The authors have declared that no competing interests exist.

vegetables, and, more recently in ready-to-eat foods [1–6]. Among other illnesses, aeromonads are responsible for gastrointestinal syndromes following ingestion of contaminated food, and wound infection following water exposure in human [1–6], and for furunculosis and septicemia in fish causing major losses in the aquaculture sector [7].

The pathogenesis of *Aeromonas* infection is only partially elucidated [1,5,6,8] although a wide repertoire of virulence factors contributing to biofilm formation, cell adherence, invasion and cytotoxicity have been described [8]. It has been suggested that only certain subsets of strains within a species, referred as “pathotypes”, produce disease in certain individuals [1,9,10]; presumably due to differences in the content of virulence-associated genes. From this assumption, many studies leading to an abundant literature searched for virulent strains from water, fish, food or clinical samples using Polymerase Chain Reaction (PCR) for detecting virulence factors [11–14], notably for epidemiologic surveys (e. g., [15–17]). The accurate inference of virulence in such an approach depends on how accurate the PCR-based methods are for detecting genetic markers of virulence in aeromonads. But these methods have often been developed from a small number of strains, few species and for a specific purpose, and their performances have never been evaluated for the whole genus. Before raising conclusions on pathotype patterns in aeromonads, PCR performances at the genus level deserve to be questioned considering the high level of genetic polymorphism in the genus *Aeromonas* [18,19].

The recent availability of whole genome sequences (WGS) data from many *Aeromonas* isolates from various species [20,21] allows the search for virulence-associated genes by sequence comparison and the evaluation of performance of molecular methods. In the present study, we aimed to evaluate the accuracy of several widely used virulence PCRs by matching PCR results against available genomes. When necessary, accuracy was improved by the genome-driven design of new PCR tools.

# Materials and methods

## Bacterial strains and DNA extraction

Thirty-nine *Aeromonas* spp. strains were chosen to cover the whole genus. They belonged to 27 species (Table 1) among the 30 validated species in the *Aeromonas* genus at the time of January 2018. The current taxonomic affiliations are indicated for every strain with previously published names of strains indicated inside braces, where appropriate. Strains were grown on Trypticase Soy Agar at 35 °C for 16–24 h, and genomic DNA was extracted using the Master-Pure™ DNA Purification Kit (Epicentre, USA).

## Genome sequences

Thirty-seven draft genomes of *Aeromonas* spp., sequenced during previous works [20–22], were provided at the Microbial Analysis, Resources and Services (MARS) facility at the University of Connecticut (Storrs, USA) (Table 1) with a high quality level of draft genome on the basis of the quality of the sequencing and assembling (e.g., number of average genomes coverage and number of scaffolds) and of the verification of the automated annotation (16 housekeeping genes and 47 ribosomal protein-coding genes checked in every genome). Two other genomes were downloaded from GenBank (Table 1). To check the taxonomical assignment of genomes employed, we performed a ML phylogenetic tree of all the genomes (n = 39) calculated by the program kSNP3.1 [25], a kmer based alignment free method (see S1 Fig).
<table><caption>Table 1. General features of the strains and genomes used in this study.</caption><thead><tr><th>Strain (n = 39)</th><th>Source of isolation</th><th>Geno-me size (Mbp)</th><th>No of scaffolds</th><th>Average genome coverage</th><th>N50 (nt)</th><th>G+C content (%)</th><th>No of predicted CDSs</th><th>Level of assem-bly</th><th>Genome accession number</th><th>Reference</th></tr></thead><tbody><tr><td>A. allosacharophila CECT 4199<sup>T</sup></td><td>Infected eel</td><td>4.66</td><td>120</td><td>87</td><td>114,541</td><td>58.4</td><td>4,173</td><td>IHQ</td><td>PRJEB7019<sup>a</sup></td><td>[20]</td></tr><tr><td>A. australiensis CECT 8023<sup>T</sup></td><td>Irrigation water system</td><td>4.11</td><td>113</td><td>128</td><td>95,095</td><td>58.1</td><td>3,733</td><td>IHQ</td><td>PRJEB7021<sup>a</sup></td><td></td></tr><tr><td>A. bestiarum CECT 4227<sup>T</sup></td><td>Fish</td><td>4.68</td><td>41</td><td>53</td><td>237,067</td><td>60.5</td><td>4,223</td><td>IHQ</td><td>PRJEB7022<sup>a</sup></td><td></td></tr><tr><td>A. bivalvium CECT 7113<sup>T</sup></td><td>Cockles</td><td>4.28</td><td>69</td><td>30</td><td>149,050</td><td>62.3</td><td>3,909</td><td>IHQ</td><td>PRJEB7023<sup>a</sup></td><td></td></tr><tr><td>A. caviae CECT 838<sup>T</sup></td><td>Guinea pig</td><td>4.47</td><td>111</td><td>95</td><td>101,663</td><td>61.6</td><td>4,081</td><td>IHQ</td><td>PRJEB7024<sup>a</sup></td><td></td></tr><tr><td>A. dhakensis CECT 7289<sup>T</sup> {A. aquariorum}<sup>†</sup></td><td>Aquaria of ornamental fish</td><td>4.69</td><td>78</td><td>117</td><td>163,504</td><td>61.7</td><td>4,266</td><td>IHQ</td><td>PRJEB7020<sup>a</sup></td><td></td></tr><tr><td>A. dhakensis CIP 107500 {A. hydrophila subsp. dhakensis}<sup>†</sup></td><td>Human diarrheic stool</td><td>4.71</td><td>73</td><td>84</td><td>165,885</td><td>61.8</td><td>4,284</td><td>IHQ</td><td>PRJEB7048<sup>a</sup></td><td></td></tr><tr><td>A. dhakensis BVH28b</td><td>Human wound</td><td>4.89</td><td>68</td><td>130</td><td>150,860</td><td>61.7</td><td>4,466</td><td>IHQ</td><td>PRJEB9016<sup>a</sup></td><td>[22]</td></tr><tr><td>A. diversa CECT 4254<sup>T</sup></td><td>Human wound</td><td>4.06</td><td>37</td><td>116</td><td>203,531</td><td>61.5</td><td>3,711</td><td>IHQ</td><td>PRJEB7026<sup>a</sup></td><td>[20]</td></tr><tr><td>A. encheleia CECT 4342<sup>T</sup></td><td>Fish</td><td>4.47</td><td>35</td><td>112</td><td>380,984</td><td>61.9</td><td>4,076</td><td>IHQ</td><td>PRJEB7027<sup>a</sup></td><td></td></tr><tr><td>A. enteropelogenes CECT 4255<sup>T</sup> {A. trota}<sup>†</sup></td><td>Human stool</td><td>4.34</td><td>27</td><td>66</td><td>640,249</td><td>60.0</td><td>3,917</td><td>IHQ</td><td>PRJEB7043<sup>a</sup></td><td></td></tr><tr><td>A. eucrenophila CECT 4224<sup>T</sup></td><td>Fresh water fish</td><td>4.54</td><td>22</td><td>50</td><td>441,212</td><td>61.1</td><td>4,113</td><td>IHQ</td><td>PRJEB7029<sup>a</sup></td><td></td></tr><tr><td>A. fluvialis LMG 24681<sup>T</sup></td><td>River water</td><td>3.90</td><td>76</td><td>48</td><td>108,949</td><td>58.2</td><td>3,609</td><td>IHQ</td><td>PRJEB7030<sup>a</sup></td><td></td></tr><tr><td>A. hydrophila subsp. hydrophila CECT 839<sup>T</sup></td><td>Tin of milk with fishy odor</td><td>4.74</td><td>1</td><td>-</td><td>4,744,448</td><td>61.5</td><td>4,119</td><td>C</td><td>CP000462<sup>b</sup></td><td>[23]</td></tr><tr><td>A. hydrophila BVH25a</td><td>Human respiratory tract</td><td>5.10</td><td>130</td><td>44</td><td>84,371</td><td>60.9</td><td>4,598</td><td>IHQ</td><td>PRJEB9013<sup>a</sup></td><td>[22]</td></tr><tr><td>A. jandaei CECT 4228<sup>T</sup></td><td>Human stool</td><td>4.50</td><td>58</td><td>55</td><td>161,393</td><td>58.7</td><td>4,065</td><td>IHQ</td><td>PRJEB7031<sup>a</sup></td><td>[20]</td></tr><tr><td>A. media CECT 4232<sup>T</sup></td><td>River water</td><td>4.48</td><td>233</td><td>60</td><td>37,608</td><td>60.9</td><td>4,075</td><td>IHQ</td><td>PRJEB7032<sup>a</sup></td><td></td></tr><tr><td>A. media LMG 13464 {A. caviae}<sup>†</sup></td><td>Infected fish</td><td>4.45</td><td>99</td><td>87</td><td>103,746</td><td>61.3</td><td>4,014</td><td>IHQ</td><td>PRJEB12347<sup>a</sup></td><td>[21]</td></tr><tr><td>A. media CECT 7111</td><td>Oyster</td><td>4.41</td><td>92</td><td>70</td><td>108,504</td><td>61.6</td><td>3,998</td><td>IHQ</td><td>PRJEB12345<sup>a</sup></td><td></td></tr><tr><td>A. molluscorum CIP 108876<sup>T</sup></td><td>Wedge-shells</td><td>4.23</td><td>309</td><td>9</td><td>21,565</td><td>59.2</td><td>3,946</td><td>IHQ</td><td>AQGQ01<sup>b</sup></td><td>[24]</td></tr><tr><td>Aeromonas sp. genomospecies paramedia CECT 8838</td><td>Human diarrheic stool</td><td>4.46</td><td>128</td><td>99</td><td>78,349</td><td>62.2</td><td>4,086</td><td>IHQ</td><td>PRJEB12349<sup>a</sup></td><td>[21]</td></tr><tr><td>A. piscicola LMG 24783<sup>T</sup></td><td>Salmon</td><td>5.18</td><td>91</td><td>99</td><td>150,424</td><td>59.0</td><td>4,713</td><td>IHQ</td><td>PRJEB7033<sup>a</sup></td><td>[20]</td></tr><tr><td>A. popoffii CIP 105493<sup>T</sup></td><td>Drinking water production plant</td><td>4.76</td><td>105</td><td>67</td><td>113,495</td><td>58.4</td><td>4,331</td><td>IHQ</td><td>PRJEB7034<sup>a</sup></td><td></td></tr><tr><td>A. rivipollensis LMG 13459 {A. caviae}<sup>†</sup></td><td>Infected fish</td><td>4.49</td><td>111</td><td>76</td><td>107,760</td><td>61.7</td><td>4,091</td><td>IHQ</td><td>PRJEB12346<sup>a</sup></td><td>[21]</td></tr><tr><td>A. rivipollensis 76c {A. media}<sup>†</sup></td><td>Human diarrheic stool</td><td>4.69</td><td>137</td><td>79</td><td>93,768</td><td>61.3</td><td>4,255</td><td>IHQ</td><td>PRJEB8966<sup>a</sup></td><td>[22]</td></tr><tr><td>A. rivipollensis BVH40 {A. media}</td><td>Human stool</td><td>4.70</td><td>123</td><td>79</td><td>105,841</td><td>61.4</td><td>4,204</td><td>IHQ</td><td>PRJEB9017<sup>a</sup></td><td></td></tr></tbody></table>

(Continued)
<table><caption>**Table 1.** (Continued)</caption><thead><tr><th>Strain (n = 39)</th><th>Source of isolation</th><th>Geno-me size (Mbp)</th><th>No of scaffolds</th><th>Average genome coverage</th><th>N50 (nt)</th><th>G+C content (%)</th><th>No of predicted CDSs</th><th>Level of assem-bly</th><th>Genome accession number</th><th>Reference</th></tr></thead><tbody><tr><td>A. rivuli DSM 22539<sup>T</sup></td><td>Freshwater</td><td>4.53</td><td>102</td><td>99</td><td>155,151</td><td>60.0</td><td>4,149</td><td>IHQ</td><td>PRJEB7035<sup>a</sup></td><td rowspan="2">[20]</td></tr><tr><td>A. salmonicida subsp. salmonicida CIP 103209<sup>T</sup></td><td>Salmon</td><td>4.74</td><td>128</td><td>117</td><td>89,543</td><td>58.5</td><td>4,442</td><td>IHQ</td><td>PRJEB7036<sup>a</sup></td></tr><tr><td>A. sanarellii LMG 24682<sup>T</sup></td><td>Human wound</td><td>4.19</td><td>98</td><td>121</td><td>82,664</td><td>63.1</td><td>3,828</td><td>IHQ</td><td>PRJEB7037<sup>a</sup></td><td></td></tr><tr><td>A. schubertii CECT 4240<sup>T</sup></td><td>Human wound</td><td>4.13</td><td>111</td><td>260</td><td>108,810</td><td>61.7</td><td>3,808</td><td>IHQ</td><td>PRJEB7038<sup>a</sup></td><td></td></tr><tr><td>A. simiae CIP 107798<sup>T</sup></td><td>Healthy monkey</td><td>3.99</td><td>100</td><td>86</td><td>73,112</td><td>61.1</td><td>3,654</td><td>IHQ</td><td>PRJEB7039<sup>a</sup></td><td></td></tr><tr><td>A. sobria CECT 4245<sup>T</sup></td><td>Fish</td><td>4.68</td><td>52</td><td>34</td><td>188,072</td><td>58.6</td><td>4,160</td><td>IHQ</td><td>PRJEB7040<sup>a</sup></td><td></td></tr><tr><td>A. taiwanensis LMG 24683<sup>T</sup></td><td>Human wound</td><td>4.24</td><td>106</td><td>66</td><td>85,294</td><td>62.8</td><td>3,884</td><td>IHQ</td><td>PRJEB7041<sup>a</sup></td><td></td></tr><tr><td>A. tecta CECT 7082<sup>T</sup></td><td>Human diarrheic stool</td><td>4.76</td><td>51</td><td>89</td><td>238,229</td><td>60.1</td><td>4,278</td><td>IHQ</td><td>PRJEB7042<sup>a</sup></td><td></td></tr><tr><td>A. veronii bv. veronii CECT 4257<sup>T</sup></td><td>Human respiratory tract</td><td>4.52</td><td>52</td><td>59</td><td>181,171</td><td>58.8</td><td>4,070</td><td>IHQ</td><td>PRJEB7044<sup>a</sup></td><td></td></tr><tr><td>A. veronii bv. sobria LMG 13067</td><td>Environment</td><td>4.74</td><td>72</td><td>46</td><td>147,470</td><td>58.3</td><td>4,171</td><td>IHQ</td><td>PRJEB7051<sup>a</sup></td><td></td></tr><tr><td>A. veronii BVH25b</td><td>Human respiratory tract</td><td>4.66</td><td>35</td><td>63</td><td>241,725</td><td>58.7</td><td>4,185</td><td>IHQ</td><td>PRJEB9014<sup>a</sup></td><td rowspan="3">[22]</td></tr><tr><td>A. veronii BVH26b</td><td>Human wound</td><td>4.58</td><td>48</td><td>73</td><td>180,501</td><td>58.7</td><td>4,107</td><td>IHQ</td><td>PRJEB9015<sup>a</sup></td></tr><tr><td>A. veronii 77c</td><td>Human diarrheic stool</td><td>4.61</td><td>42</td><td>78</td><td>230,104</td><td>58.6</td><td>4,124</td><td>IHQ</td><td>PRJEB9012<sup>a</sup></td></tr></tbody></table>

## Genome analysis

Genomic contigs and circular genomes were annotated by using the RAST annotation server to identify RNAs and protein-coding genes [26]. The genomes were queried for genes encoding for a set of virulence factors described in *Aeromonas* spp. by using reference protein sequences (Table 2), either using translated sequences of the validated subset of UniProtKB-SwissProt or annotated genes in UniprotKB-TrEMBL database. Sequence comparisons with reference protein sequences were performed with SEEDviewer that uses bidirectional protein-protein BLAST (BlastP) sequence comparison of translated open reading frames. Proteins with amino acid sequence similarities ≥65% and E-value ≤10⁻¹⁰ were considered to be homologous [27].

## Polymerase Chain Reaction (PCR)

Both selected and new primers were used to amplify DNA from aerolysin/enterotoxin cytotoxic (*aerA/act*), heat-stable cytotonic enterotoxin (*ast*), heat-labile cytotonic enterotoxin (*alt*), lipase (*lip*), serine protease (*ser*), ADP-ribosylating toxins (*aexT* and *aexU*), shiga toxins 1 (*stx1*) and 2 (*stx2a*), T3SS needle proteins (*ascFG*), T3SS inner membrane channel protein (*ascV*) and lateral flagellin A (*lafA*) genes. PCR assays were carried out in a 25 µL reaction mixture containing 0.2 µM of each primer (Sigma-Aldrich, Saint-Louis, US), 0.2 mM of each deoxynucleoside triphosphate (dNTP) (Euromedex, Souffelweyersheim, France), 2.5 mM MgCl2

<footnote><sup>a</sup>: Performed at the Microbial Analysis, Resources and Services (MARS) facility at the University of Connecticut (Storrs, USA)</footnote>

<footnote><sup>b</sup>: Obtained from GenBank, National Center for Biotechnology Information</footnote>

<footnote><sup>T</sup>: Type strain</footnote>

<footnote><sup>†</sup> Previously published names are indicated inside braces.</footnote>

<footnote>Abbreviations: IHQ, Improved high quality draft genome; C, complete genome.</footnote>
<table><caption>Table 2. Virulence-associated genes detected by genome analysis.</caption><thead><tr><th rowspan="2">Strain (n = 39)</th><th>aerA/act</th><th>ser</th><th>ast</th><th>alt/pla</th><th>stx1a, stx2a</th><th>aexT</th><th>aexU</th><th>ascFG‡</th><th>ascV</th><th>lafA</th></tr><tr><th>P09167<sup>S</sup><br>(AerA)<br>Q44063<sup>E</sup><br>(Act)</th><th>A4SNU7<sup>E</sup><br>(Ahe2) P31339<sup>S</sup><br>(AspA)<br>Q9RG23<sup>E</sup><br>(Ahe2)</th><th>Q8VRN3<sup>E</sup></th><th>Q44061<sup>E</sup><br>(Alt)<br>O87651<sup>E</sup><br>(Pla)</th><th>E2DQN2<sup>E</sup><br>(Stx1a)<br>E2DQN6<sup>E</sup><br>(Stx2a)</th><th>Q93Q17<sup>S</sup></th><th>D5LUP3<sup>E</sup></th><th>Q6WG33<sup>E</sup><br>(AscF)<br>Q6WG32<sup>E</sup><br>(AscG)</th><th>A4SUH2<sup>E</sup></th><th>Q93TL9<sup>E</sup></th></tr></thead><tbody><tr><td>A. allosacharophila<br>CECT 4199<sup>T</sup></td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. australiensis<br>CECT 8023<sup>T</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. bestiarum<br>CECT 4227<sup>T</sup></td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. bivalvium<br>CECT 7113<sup>T</sup></td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. caviae<br>CECT 838<sup>T</sup></td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. dhakensis<br>CECT 7289<sup>T</sup> {A.<br>aquariorum}<sup>†</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. dhakensis<br>CIP 107500 {A. hydrophila<br>subsp. dhakensis}<sup>†</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. dhakensis<br>BVH28b</td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. diversa<br>CECT 4254<sup>T</sup></td><td>+</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. encheleia<br>CECT 4342<sup>T</sup></td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. enteropelogenes<br>CECT 4255<sup>T</sup> {A. trota}</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. eucrenophila<br>CECT 4224<sup>T</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. fluvialis<br>LMG 24681<sup>T</sup></td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. hydrophila<br>subsp.<br>hydrophila<br>CECT 839<sup>T</sup></td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. hydrophila<br>BVH25a</td><td>-</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. jandaei<br>CECT 4228<sup>T</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. media<br>CECT 4232<sup>T</sup></td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. media<br>LMG 13464<br>{A. caviae}<sup>†</sup></td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. media<br>CECT 7111</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. molluscorum<br>CIP 108876<sup>T</sup></td><td>+</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Aeromonas sp.<br>genomospecies<br>paramedia<br>CECT 8838</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. piscicola<br>LMG 24783<sup>T</sup></td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. popoffii<br>CIP 105493<sup>T</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. rivipollensis<br>LMG 13459 {A. caviae}<sup>†</sup></td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. rivipollensis<br>76c {A.<br>media}<sup>†</sup></td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

(Continued)
<table><caption>**Table 2. (Continued)**</caption><thead><tr><th rowspan="2">Strain (n = 39)</th><th>aerA/act</th><th>ser</th><th>ast</th><th>alt/pla</th><th>stx1a, stx2a</th><th>aexT</th><th>aexU</th><th>ascFG‡</th><th>ascV</th><th>lafA</th></tr><tr><th>P09167<sup>S</sup><br>(AerA)<br>Q44063<sup>E</sup><br>(Act)</th><th>A4SNU7<sup>E</sup><br>(Ahe2) P31339<sup>S</sup><br>(AspA)<br>Q9RG23<sup>E</sup><br>(Ahe2)</th><th>Q8VRN3<sup>E</sup></th><th>Q44061<sup>E</sup><br>(Alt)<br>O87651<sup>E</sup><br>(Pla)</th><th>E2DQN2<sup>E</sup><br>(Stx1a)<br>E2DQN6<sup>E</sup><br>(Stx2a)</th><th>Q93Q17<sup>S</sup></th><th>D5LUP3<sup>E</sup></th><th>Q6WG33<sup>E</sup><br>(AscF)<br>Q6WG32<sup>E</sup><br>(AscG)</th><th>A4SUH2<sup>E</sup></th><th>Q93TL9<sup>E</sup></th></tr></thead><tbody><tr><td>A. rivipollensis BVH40<br>{A. media}†</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. rivuli DSM 22539<sup>T</sup></td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. salmonicida subsp.<br>salmonicida CIP<br>103209<sup>T</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. sanarellii LMG<br>24682<sup>T</sup></td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. schubertii CECT<br>4240<sup>T</sup></td><td>+</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. simiae CIP 107798<sup>T</sup></td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. sobria CECT 4245<sup>T</sup></td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. taiwanensis LMG<br>24683<sup>T</sup></td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td></tr><tr><td>A. tecta CECT 7082<sup>T</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. veronii bv. veronii<br>CECT 4257<sup>T</sup></td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. veronii bv. sobria<br>LMG 13067</td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. veronii BVH25b</td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A. veronii BVH26b</td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>A. veronii 77c</td><td>+</td><td>+</td><td>-</td><td>+</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr></tbody></table>

S/F: Accession numbers correspond to protein sequences in Swiss-prot (S) or TrEMBL databases (E).

T: Type strain

† Previously published names are indicated inside braces.

‡ ascF and ascG are flanking genes.

https://doi.org/10.1371/journal.pone.0201428.t002

(Promega, Madison, US), and 1.25 U of GoTaq DNA polymerase (Promega, Madison, WI) in the appropriate reaction buffer containing 1.5 mM MgCl2 and 25 ng of genomic DNA as the template. The amplification conditions performed on a GeneAmp PCR System 9700 (Applied Biosystems, Foster city, US) were as follows: initial denaturation for 3 min at 95°C, followed by (1) main amplification program of the reference source specified in the Table 3 for already published PCR, with length of extension adapted to the efficacy of the DNA polymerase used (1 min/kb) when necessary, and a final extension step at 72°C for 10 min or (2) 32 amplification cycles as indicated in Table 4 for newly developed PCR. The PCR products and the 50 bp DNA ladder (New England BioLabs, Ipswich, US) were separated in 1.5% agarose gels in 0.5X TBE buffer-ethidium bromide 500 µg/mL and revealed under UV. Positive and negative controls were added in each PCR to assess the validity and specificity of amplification reaction. The positive controls of amplification corresponded to genomic DNA of one of the strains included in the study and for which the presence of gene of interest was checked on the basis of genomic analysis as follows: A. veronii 77C (aer/act, aexT and ser), A. dhakensis CIP 107500<sup>T</sup> (alt, ascFG, ascV and lafA), A. hydrophila BVH 25a (ast). Specificity of the assays was confirmed using PCR-grade water and DNA from A. fluvialis LMG 24681<sup>T</sup> which harbors in its genome none of the genes of interest. All strains were tested for all primer sets. For all the
<table><caption>Table 3. Performance of several PCRs used for the detection of virulence-associated genes in the genus Aeromonas.</caption><thead><tr><th>Gene</th><th>Virulence factor</th><th>Prevalence estimated from WGS analysis (n = 39)</th><th>PCRRef.</th><th>Primers (nt.)</th><th>No of sequences in align-ment</th><th>No of primer mismatches† (min-max)</th><th>DiagnosticSe % (CI95%)</th><th>Comments on primers / diagnostic sensitivity relationship</th><th>DiagnosticSp (%) (CI95%)</th><th>Comments on primers / specificity relationship</th></tr></thead><tbody><tr><td rowspan="3">aerA/ act</td><td rowspan="3">Aerolysin AerA/ Cytotoxic enterotoxin Act</td><td rowspan="3">56% (n = 22)</td><td>[15]</td><td>AHCF1 (22) AHCR1 (22)</td><td>32 32</td><td>0-6 0-6</td><td>64 (41;83)</td><td>High variability in hybridization sequences of forward and reverse primers</td><td>100 (80;100)</td><td>-</td></tr><tr><td>[30]</td><td>aer-f (20) aer-r (20)</td><td>32 32</td><td>0-4 0-2</td><td>91 (71;99)</td><td>-</td><td>82 (57;96)</td><td>aer-r: high GC content and self-end dimers</td></tr><tr><td>This study</td><td>aer-1FX (20) aer-2R (19)</td><td>32 32</td><td>0-1 0-2</td><td>100 (85;100)</td><td>aer-1FX: degenerate nucleotides in the forward primer aer-f</td><td>94 (71;100)</td><td>Degenerate bases did not decrease specificity aer-2R: identical to aer-r but one 5' end base deleted</td></tr><tr><td rowspan="2">ser</td><td rowspan="2">Serine protease</td><td rowspan="2">69% (n = 27)</td><td>[13]</td><td>Serine-f (20) Serine-r (20)</td><td>39 39</td><td>0-7 0-11</td><td>59 (39;78)</td><td>High variability in hybridization sequences of forward and reverse primers One nucleotide deleted in the reference sequences used to design the reverse primer (ENA X67043 and ENA AAF22245)</td><td>92 (62;100)</td><td>-</td></tr><tr><td>This study</td><td>ser-1FX (21) ser-2RX (20)</td><td>39 39</td><td>0-2 0-0</td><td>96 (81;100)</td><td>Selection of less variable sequences Degenerate bases</td><td>92 (62;100)</td><td>Degenerate bases did not decrease specificity</td></tr><tr><td>ast</td><td>Cytotonic heat-stable enterotoxin</td><td>15% (n = 6)</td><td>[31]</td><td>ast-F (20) ast-R (20)</td><td>7 7</td><td>0-4 0-3</td><td>100 (54;100)</td><td>-</td><td>91 (76;98)</td><td>-</td></tr><tr><td rowspan="2">alt</td><td rowspan="2">cytotonic heat-labile enterotoxin</td><td rowspan="2">97% (n = 38)</td><td>[31]</td><td>alt-F (20) alt-R (19)</td><td>39 39</td><td>0-10 0-5</td><td>34 (20;51)</td><td>High variability in hybridization sequences of forward and reverse primers</td><td>-</td><td>Unevaluable specificity: insufficient number of strains for which alt was not detected in WGS (1/39)</td></tr><tr><td>[13]</td><td>lip-f (18) lip-r (20)</td><td>39 39</td><td>0-1 0-2</td><td>68 (51;83)</td><td>-</td><td>-</td><td></td></tr><tr><td>stx1a</td><td>Shiga toxin 1 subunit A</td><td>0% (n = 0)</td><td>[32]</td><td>Stx1-a (20) Stx1-b (20)</td><td>-</td><td>-</td><td>-</td><td>Not evaluable sensitivity: no strain for which stx1a or stx2a was detected in WGS</td><td>100 (82;100)</td><td>-</td></tr><tr><td>stx2a</td><td>Shiga toxin 2 subunit A</td><td>0% (n = 0)</td><td>[33]</td><td>S2Aup (19) S2Alp (19)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>100 (82;100)</td><td>-</td></tr><tr><td>aexT</td><td>ADP-ribosylating toxin AexT</td><td>5% (n = 2)</td><td>[34]</td><td>RASEXOS-L (18) RASEXOS-R (18)</td><td>9 9</td><td>0-1 0-1</td><td>2/2</td><td>-</td><td>100 (91;100)</td><td>-</td></tr><tr><td>aexT and/or aexU</td><td>ADP-ribosylating toxins</td><td>23% (n = 9)</td><td>This study</td><td>aexTU-1FX (19) aexTU-2RX (20)</td><td>21 21</td><td>0-2 0-2</td><td>100 (66;100)</td><td>-</td><td>100 (88;100)</td><td>-</td></tr><tr><td rowspan="2">asc-FG</td><td rowspan="2">T3SS needle proteins</td><td rowspan="2">28% (n = 11)</td><td>[16]</td><td>ascF-G-fwd (20) ascF-G-rev (20)</td><td>38 23</td><td>0-5 na</td><td>45 (17;77)</td><td>High variability in hybridization sequences of forward and reverse primers</td><td>74 (54;89)</td><td>ascF-G-fwd: high GC content, self-end dimers and 3'-end unstability</td></tr><tr><td>This study</td><td>ascFG-1F (21) ascFG-2RX (20)</td><td>38 23</td><td>0-4 0-1</td><td>91 (59;100)</td><td>Selection of less variable sequences Degenerate bases</td><td>100 (88;100)</td><td>Degenerate bases did not decrease specificity</td></tr></tbody></table>

(Continued)
<table><caption>**Table 3.** (Continued)</caption><thead><tr><th>Gene</th><th>Virulence factor</th><th>Prevalence estimated from WGS analysis (n = 39)</th><th>PCRRef.</th><th>Primers (nt.)</th><th>No of sequences in align-ment</th><th>No of primer mis-matches† (min-max)</th><th>DiagnosticSe % (CI95%)</th><th>Comments on primers / diagnostic sensitivity relationship</th><th>DiagnosticSp (%) (CI95%)</th><th>Comments on primers / specificity relationship</th></tr></thead><tbody><tr><td rowspan="2">asc-V</td><td rowspan="2">T3SS inner membrane channel protein</td><td rowspan="2">28% (n = 11)</td><td>[16]</td><td>ascV-fwd (20)</td><td>21</td><td>1-2</td><td rowspan="2">55 (23;83)</td><td rowspan="2">High variability in hybridization sequences of forward and reverse primers</td><td rowspan="2">71 (51;87)</td><td rowspan="2">ascV-fwd: high GC content</td></tr><tr><td>ascV-rev (20)</td><td>21</td><td>4-7</td></tr><tr><td rowspan="2">lafA</td><td rowspan="2">Lateral flagellin A</td><td rowspan="2">54% (n = 21)</td><td>This study</td><td>ascV-1F (19)</td><td>21</td><td>0-1</td><td rowspan="2">100 (72;100)</td><td rowspan="2">Selection of less variable sequences</td><td rowspan="2">100 (88;100)</td><td rowspan="2">-</td></tr><tr><td>ascV-2R (20)</td><td>21</td><td>0-2</td></tr><tr><td rowspan="2"></td><td rowspan="2"></td><td rowspan="2"></td><td>[35]</td><td>Laf1 (18)</td><td>58</td><td>4-6</td><td rowspan="2">55 (34;80)</td><td rowspan="2">High variability in hybridization sequences of forward and reverse primers</td><td rowspan="2">100 (82;100)</td><td rowspan="2">-</td></tr><tr><td>Laf2 (17)</td><td>58</td><td>0-3</td></tr><tr><td rowspan="2"></td><td rowspan="2"></td><td rowspan="2"></td><td>This study</td><td>lafA-1FX (19)</td><td>58</td><td>0-5</td><td rowspan="2">86 (64;97)</td><td rowspan="2">Selection of less variable sequences</td><td rowspan="2">100 (81;100)</td><td rowspan="2">Degenerate bases did not decrease specificity</td></tr><tr><td>lafA-2RX (18)</td><td>58</td><td>0-3</td></tr></tbody></table>

Abbreviations: CI95%, confidence interval 95%; na, not available; PCR, Polymerase Chain Reaction; Ref., References; WGS, Whole Genome Sequences.

https://doi.org/10.1371/journal.pone.0201428.t003

PCR assays and for each strain, the results of absence/presence of amplification was determined by two different readers and in two independent experiments on the basis of the expected size of amplified products and of the positive control.

## Primer design strategy

Nucleic sequences were aligned using the Clustal ω2 program within Seaview 4 package [28] (see S1 Files). Primers displaying a length of 18–21 nucleotides were designed from the conserved regions, and AmplifX 1.7.0 (CNRS, Aix-Marseille Université, France, http://crn2m.univ-mrs.fr/pub/amplifx-dist) was used to assess their intrinsic properties such as complexity (polyX and triplet repetitions), 3' stability and self-dimer formation. The specificity of the new primers was checked with NCBI nucleotide-nucleotide BLAST.

## Challenging PCR accuracy characteristics

The accuracy of the PCR assays were evaluated by comparing the results of the PCR with those obtained by genome analysis, with WGS results considered as gold standard. Diagnostic sensitivity was defined as the percentage of 'positive PCR results among all the strains for which the gene was detected in WGS', while specificity was defined as the percentage of 'negative PCR results among all the strains for which the gene was not detected in WGS'. Diagnostic sensitivity and specificity values below 90% were considered as insufficient. Exact confidence interval of a frequency, i.e., binomial confidence interval, was determined for each percentage, as described elsewhere [29] by using the web server http://statpages.info/confint.html.

<footnote>†Evaluated from the virulence-associated-gene/primers multiple alignments</footnote>

<footnote>‡ Evaluated from *lafA* multiple alignments</footnote>
<table><caption>Table 4. Primers designed in this study and amplification conditions of PCRs.</caption><thead><tr><th rowspan="2">Virulence-associated gene</th><th rowspan="2">Primer sequences (5'-3')</th><th rowspan="2">Amplified sequence length (bp)</th><th colspan="3">Amplification conditions†</th></tr><tr><th>°C</th><th>Time (sec)</th><th>Number of cycles</th></tr></thead><tbody><tr><td rowspan="3">aerA<br>(syn: act)</td><td>aer-1FX: 5'-CCTAYRGYCTSAGCGAGAAG-3'</td><td rowspan="3">430</td><td>94</td><td>60</td><td rowspan="3">32</td></tr><tr><td>aer-2R: 5'-CAGTTCCAGTCCCACCACT-3'</td><td>56</td><td>60</td></tr><tr><td></td><td>72</td><td>60</td></tr><tr><td rowspan="3">ser</td><td>ser-1FX: 5'-GACAAYCGVGTSTTCAAAGAG-3'</td><td rowspan="3">262</td><td>94</td><td>60</td><td rowspan="3">32</td></tr><tr><td>ser-2RX: 5'-ACCACCARGTTCAGAAGTT-3'</td><td>59</td><td>60</td></tr><tr><td></td><td>72</td><td>30</td></tr><tr><td rowspan="3">aexT/aexU</td><td>aexTU-1FX: 5'-TGGCVMTSAAAGAGTGGAT-3'</td><td rowspan="3">225</td><td>94</td><td>60</td><td rowspan="3">32</td></tr><tr><td>aexTU-2RX: 5'-GCARDGSRCCRTTGCCRGTC-3'</td><td>63</td><td>60</td></tr><tr><td></td><td>72</td><td>45</td></tr><tr><td rowspan="3">asc-FG</td><td>ascFG-1F: 5'-CAAGATCAACAAATGGTCGGT-3'</td><td rowspan="3">262</td><td>94</td><td>60</td><td rowspan="3">32</td></tr><tr><td>ascFG-2RX: 5'-TTCAAYCARRGADGASAGGCG-3'</td><td>60</td><td>60</td></tr><tr><td></td><td>72</td><td>30</td></tr><tr><td rowspan="3">asc-V</td><td>ascV-1F: 5'-CGCAAGGACATCATGCTGG-3'</td><td rowspan="3">578</td><td>94</td><td>60</td><td rowspan="3">32</td></tr><tr><td>ascV-2R: 5'-ATGATGATGAGGCCCGCGAT-3'</td><td>65</td><td>60</td></tr><tr><td></td><td>72</td><td>45</td></tr><tr><td rowspan="3">lafA</td><td>lafA-1FX: 5'-GATGYTGRGYACYGCCATG-3'</td><td rowspan="3">619</td><td>94</td><td>60</td><td rowspan="3">32</td></tr><tr><td>lafA-2RX: 5'-CATRTTGGARAGGTTRTTGAC-3'</td><td>64</td><td>60</td></tr><tr><td></td><td>72</td><td>60</td></tr></tbody></table>

Degenerate bases are underlined.

https://doi.org/10.1371/journal.pone.0201428.t004

# Results

## Virulence-associated gene selection

The literature analysis was used to inventory virulence genes commonly searched in *Aeromonas* and to identify the reference of the PCRs used in the different studies, as summarized in Fig 1. We selected 61 original research studies that performed virulence-associated gene PCRs for populations of at least 15 *Aeromonas* strains (Fig 1). From these 61 population studies, the origin of strains and distribution of *Aeromonas* species, the virulence-associated genes and the references of the PCRs method used are presented in the S1 Table. Among these 61 articles (see S1 Table for references), the most studied genes were *aerA* (n = 51), *alt* (n = 35), *act* (n = 30), *ast* (n = 30), *lip* (n = 19), *hlyA* (n = 17), *fla* (n = 15), *aexT* (n = 11), *aexU* (n = 2), *ascV* (n = 11), *ser* (n = 10), *gcat* (n = 9), *lipH3* (n = 8), *pla* (n = 8), *lafA* (n = 8), *ascFG* (n = 6), *exu* (n = 6), *ahyB* (n = 5), *stx1* (n = 5), *stx2* (n = 5) and *alp-1* (n = 5). These genes corresponded to different categories of virulence factors that are included in the pathotypes definition: hemolysins and related genes (*aerA*, *act*, *hlyA*), cytotonic enterotoxins genes (*alt*, *ast*), lipases genes (*lip*, *gcat*, *lipH3*, *pla*, *alp-1*), flagellar appendages genes (*fla*, *lafA*), type III secretion system effectors *AexT* and *AexU* genes (*aexT*, *aexU*), type III secretion system components genes (*ascV*, *ascFG*), serine protease gene (*ser*), nuclease gene (*exu*), elastase gene (*ahyB*) and shiga-toxins genes (*stx1*, *stx2*). We selected one or several genes from these virulence categories: *aerA* (aerolysin), *act* (cytolytic enterotoxin), *ast* (heat-stable cytotonic enterotoxin), *alt* (heat-labile cytotonic enterotoxin), *pla* (phospholipase A), *ser* (serine protease), *aexT* and *aexU* (ADP-ribosylating toxins), *stx1* and *stx2* (shiga-toxins), *ascFG* (type III secretion system needle

<footnote>†The amplification conditions of PCR are preceded by an initial denaturation step at 95°C during 3 minutes and followed by a final elongation step at 72°C during 10 minutes.</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Flowchart for the literature analysis. The automatic Pubmed query using «*Aeromonas* », « PCR » and « virulence » as keywords generated 126 articles (October 27, 2016), among which 46 corresponded to the entry criteria. The reading of the relevant references quoted in these 46 selected articles led to the manual inclusion of 15 additional articles.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201428.g001

proteins), *ascV* (type III secretion system inner membrane channel protein) and *lafA* (lateral flagellin A).

## Genome content in virulence-associated genes

We detected homologues for all the studied genes except for *stx1* and *stx2* (Table 2). The reference sequences of aerolysin gene *aerA* from *A. bestiarum Ah65* (considered in the original source as *A. hydrophila Ah65*) (SwissProt P09167) and Cytolytic enterotoxin gene *act* from *A. dhakensis* SSU (considered in the original source as *A. hydrophila* SSU) (UniprotKB-TrEMBL Q44063) shared a very high level of identity (95.1%) and similarity (97.2%) in protein sequence (493 amino-acids both). Similarly, the phospholipase A gene *pla* from *A. piscicola* AH-3 (considered in the original source as *A. hydrophila* AH-3) (UniprotKB-TrEMBL O87651) and the heat-labile cytotonic enterotoxin gene *alt* from *A. dhakensis* SSU (UniprotKB-TrEMBL Q44061) shared a very high level of identity (96.9%) and similarity (98.5%) in protein sequence, from an alignment of 89% query-cover of *alt* coding-sequence (368/385 amino-acids), *pla* coding-sequence spanning 805 amino-acids. These homology data and the presence
of common flanking genes (data not shown), led to consider *aerA/act* and *alt/pla* as two single virulence-associated genes for all the strains included in the study. The *ser* gene that is putatively involved in the aerolysin activation by proteolytic cleavage was detected in 14 out of 17 *aerA/act* positive strains (Table 2). The genes coding for structural components of type III secretion system (T3SS) were detected simultaneously in the 11 T3SS gene positive strains (Table 2). The 248 residues from the N-terminal domain of the two ADP-ribosyltransferase toxins, AexT from *A. salmonicida* subsp. *salmonicida* ATCC 33658<sup>T</sup> (SwissProt Q93Q17; 475 amino-acids) and AexU from *A. veronii* bv. *sobria* AeG1 (UniprotKB-TrEMBL D5LUP3, 512 amino-acids) presented a high level of identity (89.5%) and similarity (94.8%) in protein sequence, but they differed significantly in their C-terminal domain (&lt; 20% similarity). All the strains that harbor the T3SS effector genes *aexT/aexU*, except *A. salmonicida* subsp. *salmonicida* CIP 103209<sup>T</sup>, also possessed the T3SS structural components. For the lateral flagellin gene, from 1 to 4 copies were detected in the genomes of the 21-*lafA* positive strains.

## PCR performance assessment

Results of the evaluation of PCR accuracy are presented in Table 3. Ability to amplify gene of interest was assessed for all genes but *stx1a* or *stx2a* because there were no strains for which the genes were detected in WGS (Table 2). The ability of the PCRs to amplify the gene of interest (i.e., diagnostic sensitivity, hereafter sensitivity) ranged from 34 to 100%. The PCR screening *aerA/act* that uses primers “aer”, and the PCRs screening the *ast* and *aexT* genes displayed high sensitivities with the tested set of strains, with values of 91%, 100% and 100%, respectively. However, the corresponding confidence intervals for *ast* and *aexT* PCR were large (range &gt;40%) because the number of positive genomes was low, 6 and 2 out of the 39 genomes studied, respectively. On the opposite, with a sensitivity value of 64%, the PCR usually employed for screening *aerA/act* genes using primers “AHC” was disappointing and was associated with frequent false negative results. Similarly, low values of sensitivity were observed for the following genes: *ser* (59%), *alt/pla* (primers “alt”, 34%; primers “lip”, 68%), *ascFG* (45%), *ascV* (55%) and *lafA* (55%). Depending on the gene, the lack of sensitivity resulted from either an excessive variability in annealing of forward and/or reverse primers, or even errors in primer sequence transcript (e.g., Serine-r and Laf1; Table 3).

The specificities of the PCRs that could be evaluated ranged from 71% to 100% depending on the gene considered (Table 3). Specificity was not determined for *alt* because only one out of the 39 WGS lacked *alt* gene (*A. fluvialis* LMG 24681<sup>T</sup>). PCRs used for screening the following genes were highly specific: *aerA/act*-using primers “AHC” (100%), *ser* (92%), *ast* (91%), *aexT* (100%) and *lafA* (100%). The test method was insufficiently specific for *aerA/act* genes using primers “aer” (82%), *ascF-G* (74%) and *ascV* (71%) genes. Therefore, an overestimation of the prevalence of these virulence-associated factors by generation of false positive results was expected. Depending on the gene, the default of specificity likely resulted from PCR primers properties such as high GC content, self-end dimer formation and/or 3'-end instability (Table 3).

## Genome driven design and evaluation of new PCR primers

Considering the low performance of some PCRs, we designed new primer pairs on the basis of multiple alignments (Table 4) with partial and full-length sequences retrieved from Genbank or EMBL databases (S2 Table), and from the WGS provided in this study (Table 1). To design the new primers, we selected low variable areas of sequence alignments and used degenerate bases in case of polymorphism. Newly designed primers were associated with an increased sensitivity of the virulence gene detection (Table 3). Meanwhile, the use of degenerate primers did
not impair the specificity (Table 3). For screening of ADP-ribosylating toxins, we have
designed primers matching in the homologous part shared by aexT and aexU genes that were
associated with a sensitivity of 100%.

Discussion

Many virulence factors are described in the genus *Aeromonas*, but to date, the pathogenicity of specific strains still cannot be predicted from the genome content. The genus *Aeromonas* is characterized by intraspecific and interspecific genetic diversities of both ribosomal and housekeeping genes that impair the definition of species and lead to the population structure of the genus in several complexes of species [19,36]. However, using the whole genome sequences and the average nucleotide identity index to compare them, all the species could be clearly delineated [20,21]. The huge variability observed in some housekeeping genes also occurs in virulence-associated genes, as already shown for the “S-layer protein” gene (*vapA*) that harbored variability up to 15% at the amino-acid level among the strains of different sub-species of *A. salmonicida* [18].

This study focuses on the performances of molecular techniques used for the detection of virulence-associated genes in *Aeromonas* strains. Some previous work aiming to detect virulence factors in the genus *Aeromonas* reported discrepancies between studies due to molecular methods used (e.g., [37]). The choice of strains used for the primer design was particularly questioned. The present study does not question the validity of previously published PCRs, especially when the original study aimed to screen the presence of virulence-associated genes in particular groups or species (e.g., in *A. caviae*, [38]), but aims to evaluate PCR performances at the genus level. The comparison of PCR results and high-quality WGS with automatic annotation and careful manual checking allows for the first time the proper evaluation of the PCR performances for detecting virulence-associated genes in aeromonads. Most new primers and PCR conditions proposed herein improve PCR performances although the design of degenerate primers was required for some virulence genes, i.e. when sequence polymorphisms are scattered all along the gene. Neither diagnostic sensitivity nor specificity was hampered in these cases.

Similar to most other population studies, our study suffers from some limitations. Studying whether strains are virulent may be challenging for opportunistic pathogens but at first, tools need to be accurate in establishing whether a gene is present or absent from genomes. Our study may have some potential bias in the strain sampling. But the collection of strains included in the study aimed to be as representative of the diversity found in the genus as possible, covering the interspecies diversity in the whole genus, given the studied genes are sought in a wide range of species [39], and the intra-species diversity of species that are the most frequently isolated in human clinical contexts [1,40]. Illustrating the unavoidable sampling bias, the shiga-toxin coding genes stx were not detected in this study by PCR nor by genome analysis. However, the stx1 gene seems to be quite rare in the genus [41]. Another explanation would be stx instability due to the mobility of the stx phage that can lead to the gene loss after subculture [41] or may occur during storage before genome sequencing and PCR. The prevalence of some virulence markers was either too low (e.g., stx) or too high (e.g., *lip*, *alt*) so that PCR assay diagnostic sensitivity or specificity could not be evaluated, respectively. Despite the limitation in strain and sequence collections, this study provides a meaningful highlight in accuracy of widely used tools together with virulence gene content in a large collection of strains.

Besides performance assessment and method improvements, this study provides some
interesting data on virulence in aeromonads, like the possible matching between several
virulence markers identified in different *Aeromonas* species that could be established or confirmed. For instance, the homology between aerolysin gene *aerA* from *A. hydrophila* and cytotoxic enterotoxin gene *act* from *A. dhakensis* has previously been considered [42]. In fact, both toxins exhibit hemolytic, cytotoxic and enterotoxic activities and both cause lethality in mice [43–47]. Moreover, a similar mechanism of action, i.e., oligomerization followed by pore-formation, has been reported for these 2 toxins [46,48]. These two toxins are closely related but it is still possible that the aerolysin/cytotonic enterotoxin from *A. dhakensis* SSU possesses some structural and functional originality [46], and further characterization of the purified homologous proteins is required to clarify this point. The *aerA/act* PCR based on newly designed primers improved the diagnostic sensitivity for the panel of aeromonads covering the whole genus from 64–91% to 100%, and dramatically improved sensitivity for the species *A. dhakensis* (from 0–67% to 100%). Given the clinical importance of this species [49], an optimized performance for detecting the well-studied *aerA* virulence gene is critical, and should provide advances in the patho-epidemiological knowledge of the species. The putative correspondence between phospholipase A gene (*pla*) from *A. piscicola* AH-3 and heat-labile cytotonic enterotoxin gene (*alt*) from *A. dhakensis* SSU, was also observed by Balsalobre et al. [50]. In this case, the nomenclature of the virulence factor depends on the study [51,52]. The biological features are possibly divergent as mentioned by Merino et al. [52] and could depend on the allele of the gene *alt/pla* but this requires further studies.

Biological interpretation of the presence of virulence genes needs to be cautious because the presence of a virulence gene does not imply its expression in the host. In addition, we have to consider the pathogenic and non-pathogenic interactions depending on the context. For example, the aeromonads T3SS system is important for both pathogenicity and mutualism as demonstrated inside the leech microbiota [53]. In parallel, post-translational modifications for activation or effector translocation should be examined. The joint presence of the *ser* gene and *aer/act* should be considered because aerolysin is activated by a serine protease [54]. The *ser* gene was absent from the genomic data of three *aerA/act* positive strains included in this study: *A. diversa* CECT 4254T, *A. molluscorum* CIP 108876T and *A. schubertii* CECT 4240T. Similarly, the presence of T3SS effectors AexT/U should be considered with respect to the presence of T3SS components, e.g., *ascFG* and *ascV*. Their presence is required for the delivery of AexT/U [55,56]. The presence of both genes *aexT* and *aexU* that we observed with the strains *A. veronii* 77C and BVH26b has already been mentioned for other *A. veronii* isolates and appears to be a frequent situation for strains carrying T3SS in the *A. veronii* group [37]. Therefore, for screening *aexT* and *aexU* genes in aeromonads, we propose an *aexT/U* PCR, followed by a molecular search of the *aexT* gene alone when the screening assay is positive.

In epidemiologic and pathogenesis studies of *Aeromonas*, large collections of strains were screened for their content in virulence-associated genes. Given the polymorphisms in virulence-associated genes highlighted in this study and because PCR assays to detect virulence associated genes were widely used without prior evaluation of their performance (diagnostic sensitivity and specificity) at the genus level (S1 Table), the published prevalence of virulence genes should be interpreted with caution. When the PCR performance is poor, the estimated prevalence is blurred with assay errors and does not reflect the true prevalence of the virulence factors. Importantly, assumptions on the pathogenic behavior of aeromonads have been made from virulence patterns obtained from defective/inaccurate tool results, and this may have led to wrong conclusions. This may contribute to the current complexity of virulence patterns in aeromonads, to the poor understanding of aeromonad virulence and to the frequent lack of links between the virulence factor pattern and the pathogenic behavior. Optimization of PCR conditions and primer is crucial to avoid biased data in applied and clinical microbiology, and in epidemiology.
The bacterial genomes, today widespread and easy to obtain, have a very high potential in numerous research areas, including the design of consolidated microbiological diagnostic tools, although they are maybe currently underused for this purpose. In our study, genomes allowed an accuracy evaluation of virulence-associated gene PCRs commonly used in the genus *Aeromonas* and their improvement by designing several improved oligonucleotide primers. WGS could obviously be used as primary data for searching virulence associated genes but WGS are still expensive when a large panel of strains is studied. Therefore, we propose here an accurate approach based on these genomic data to improve PCR methods aimed to detect virulence-associated genes in environmental bacteria.

Further advances in the knowledge of the genetic diversity and of the evolution of virulence-associated genes should improve the understanding of aeromonad adaptation to pathogenic behavior or patho-adaptation, and thereby the revision of pathotype definition.

**Supporting information**

**S1 Fig. Genome-based ML phylogenetic tree.**
(DOCX)

**S1 Files. Genomic DNA and primer alignments.**
(ZIP)

**S1 Table. Literature survey and analysis.**
(DOCX)

**S2 Table. Additional nucleic sequence accession numbers used to design new primers.**
(DOCX)

**Author Contributions**

**Conceptualization:** Emilie Talagrand-Reboul, Brigitte Lamy.

**Data curation:** Emilie Talagrand-Reboul, Fadua Latif-Eugenín, Roxana Beaz-Hidalgo, Sophie Colston.

**Formal analysis:** Emilie Talagrand-Reboul.

**Investigation:** Emilie Talagrand-Reboul, Fadua Latif-Eugenín.

**Methodology:** Emilie Talagrand-Reboul, Brigitte Lamy.

**Resources:** Maria-Jose Figueras, Joerg Graf, Estelle Jumas-Bilak, Brigitte Lamy.

**Visualization:** Emilie Talagrand-Reboul, Estelle Jumas-Bilak, Brigitte Lamy.

**Writing – original draft:** Emilie Talagrand-Reboul, Estelle Jumas-Bilak, Brigitte Lamy.

**Writing – review & editing:** Emilie Talagrand-Reboul, Fadua Latif-Eugenín, Roxana Beaz-Hidalgo, Sophie Colston, Maria-Jose Figueras, Joerg Graf, Estelle Jumas-Bilak, Brigitte Lamy.

**References**

1. Janda JM, Abbott SL. The genus Aeromonas: taxonomy, pathogenicity, and infection. Clin Microbiol Rev. 2010; 23: 35–73. https://doi.org/10.1128/CMR.00039-09 PMID: 20065325

2. Khajanchi BK, Fadl AA, Borchardt MA, Berg RL, Horneman AJ, Stemper ME, et al. Distribution of Virulence Factors and Molecular Fingerprinting of Aeromonas Species Isolates from Water and Clinical Samples: Suggestive Evidence of Water-to-Human Transmission. Appl Environ Microbiol. 2010; 76: 2313–2325. https://doi.org/10.1128/AEM.02535-09 PMID: 20154106
3. Parker JL, Shaw JG. Aeromonas spp. clinical microbiology and disease. J Infect. 2011; 62: 109–118.
https://doi.org/10.1016/j.jinf.2010.12.003 PMID: 21163298

4. Nyenje ME, Odjadjare CE, Tanih NF, Green E, Ndip RN. Foodborne pathogens recovered from ready-to-eat foods from roadside cafeterias and retail outlets in Alice, Eastern Cape Province, South Africa: public health implications. Int J Environ Res Public Health. 2012; 9: 2608–2619. https://doi.org/10.3390/ijerph9082608 PMID: 23066386

5. Figueras MJ, Beaz-Hidalgo R. Aeromonas infections in humans. Aeromonas. Caister Academic Press, Norfolk, UK, ed Graf J; 2015. pp. 65–68.

6. Teunis P, Figueras MJ. Reassessment of the Enteropathogenicity of Mesophilic Aeromonas Species. Front Microbiol. 2016; 7: 1395. https://doi.org/10.3389/fmicb.2016.01395 PMID: 27708621

7. Hossain MJ, Sun D, McGarey DJ, Wrenn S, Alexander LM, Martino ME, et al. An Asian Origin of Virulent Aeromonas hydrophila Responsible for Disease Epidemics in United States-Farmed Catfish. mBio. 2014; 5: e00848–14. https://doi.org/10.1128/mBio.00848-14 PMID: 24895303

8. Tomás JM. The main Aeromonas pathogenic factors. ISRN Microbiol. 2012; 2012: 256261. https://doi.org/10.5402/2012/256261 PMID: 23724321

9. Janda JM, Abbott SL. Evolving concepts regarding the genus Aeromonas: an expanding Panorama of species, disease presentations, and unanswered questions. Clin Infect Dis Off Publ Infect Dis Soc Am. 1998; 27: 332–344.

10. Joseph S.W., Carnahan A. M. Update on the genus Aeromonas. ASM News. 2000; 218–223.

11. Albert MJ, Ansaruzzaman M, Talukder KA, Chopra AK, Kuhn I, Rahman M, et al. Prevalence of entero-toxin genes in Aeromonas spp. isolated from children with diarrhea, healthy controls, and the environment. J Clin Microbiol. 2000; 38: 3785–3790. PMID: 11015403

12. Martins LM, Marquez RF, Yano T. Incidence of toxic Aeromonas isolated from food and human infection. FEMS Immunol Med Microbiol. 2002; 32: 237–242. PMID: 11934569

13. Chacón MR, Figueras MJ, Castro-Escarpulli G, Soler L, Guarro J. Distribution of virulence genes in clinical and environmental isolates of Aeromonas spp. Antonie Van Leeuwenhoek. 2003; 84: 269–278.
PMID: 14574104

14. Senderovich Y, Ken-Dror S, Vainblat I, Blau D, Izhaki I, Halpern M. A molecular study on the prevalence and virulence potential of Aeromonas spp. recovered from patients suffering from diarrhea in Israel. PloS One. 2012; 7: e30070. https://doi.org/10.1371/journal.pone.0030070 PMID: 22355306

15. Kingombe CI, Huys G, Tonolla M, Albert MJ, Swings J, Peduzzi R, et al. PCR detection, characterization, and distribution of virulence genes in Aeromonas spp. Appl Environ Microbiol. 1999; 65: 5293–5302. PMID: 10583979

16. Chacón MR, Soler L, Groisman EA, Guarro J, Figueras MJ. Type III secretion system genes in clinical Aeromonas isolates. J Clin Microbiol. 2004; 42: 1285–1287. https://doi.org/10.1128/JCM.42.3.1285-1287.2004 PMID: 15004096

17. Nawaz M, Khan SA, Khan AA, Sung K, Tran Q, Kerdahi K, et al. Detection and characterization of virulence genes and integrons in Aeromonas veronii isolated from catfish. Food Microbiol. 2010; 27: 327–331. https://doi.org/10.1016/j.fm.2009.11.007 PMID: 20227596

18. Lund V, Mikkelsen H. Genetic diversity among A-proteins of atypical strains of Aeromonas salmonicida. Dis Aquat Organ. 2004; 61: 257–262. https://doi.org/10.3354/dao061257 PMID: 15609880

19. Roger F, Marchandin H, Jumas-Bilak E, Kodjo A, Lamy B. Multilocus genetics to reconstruct aeromonad evolution. BMC Microbiol. 2012; 12: 62. https://doi.org/10.1186/1471-2180-12-62 PMID: 22545815

20. Colston SM, Fullmer MS, Beka L, Lamy B, Gogarten JP, Graf J. Bioinformatic genome comparisons for taxonomic and phylogenetic assignments using Aeromonas as a test case. mBio. 2014; 5. https://doi.org/10.1128/mBio.02136-14 PMID: 25406383

21. Talagrand-Reboul E, Roger F, Kimper J-L, Colston SM, Graf J, Latif-Eugenín F, et al. Delineation of Taxonomic Species within Complex of Species: Aeromonas media and Related Species as a Test Case. Front Microbiol. 2017; 8: 621. https://doi.org/10.3389/fmicb.2017.00621 PMID: 28458658

22. Mosser T, Talagrand-Reboul E, Colston SM, Graf J, Figueras MJ, Jumas-Bilak E, et al. Exposure to pairs of Aeromonas strains enhances virulence in the Caenorhabditis elegans infection model. Front Microbiol. 2015; 6: 1218. https://doi.org/10.3389/fmicb.2015.01218 PMID: 26583012

23. Seshadri R, Joseph SW, Chopra AK, Sha J, Shaw J, Graf J, et al. Genome sequence of Aeromonas hydrophila ATCC 7966T: jack of all trades. J Bacteriol. 2006; 188: 8272–8282. https://doi.org/10.1128/JB.00621-06 PMID: 16980456

24. Spataro N, Farfán M, Albarral V, Sanglas A, Lorén JG, Fusté MC, et al. Draft Genome Sequence of Aeromonas molluscorum Strain 848TT, Isolated from Bivalve Molluscs. Genome Announc. 2013; 1. https://doi.org/10.1128/genomeA.00382-13 PMID: 23788549
25. Gardner SN, Slezak T, Hall BG. kSNP3.0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome. *Bioinforma Oxf Engl.* 2015; 31: 2877–2878. https://doi.org/10.1093/bioinformatics/btv271 PMID: 25913206

26. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). *Nucleic Acids Res.* 2014; 42: D206–214. https://doi.org/10.1093/nar/gkt1226 PMID: 24293654

27. Altschul SF, Lipman DJ. Protein database searches for multiple alignments. *Proc Natl Acad Sci U S A.* 1990; 87: 5509–5513. PMID: 2196570

28. Gouy M, Guindon S, Gascuel O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. *Mol Biol Evol.* 2010; 27: 221–224. https://doi.org/10.1093/molbev/msp259 PMID: 19854763

29. Lamy B, Delignette-Muller ML, Baty F, Carret G. Simple table for estimating confidence interval of discrepancy frequencies in microbiological safety evaluation. *J Microbiol Methods.* 2004; 56: 137–139. https://doi.org/10.1016/j.mimet.2003.10.001 PMID: 14706759

30. Soler L, Figueras MJ, Chacón MR, Vila J, Marco F, Martinez-Murcia AJ, et al. Potential virulence and antimicrobial susceptibility of Aeromonas popoffii recovered from freshwater and seawater. *FEMS Immunol Med Microbiol. 2002*; 32: 243–247. PMID: 11934570

31. Aguilera-Areola MG, Hernández-Rodríguez C, Zúñiga G, Figueras MJ, Castro-Escarpulli G. Aeromonas hydrophila clinical and environmental ecotypes as revealed by genetic diversity and virulence genes. *FEMS Microbiol Lett.* 2005; 242: 231–240. https://doi.org/10.1016/j.femsle.2004.11.011 PMID: 15621443

32. Wang G, Clark CG, Rodgers FG. Detection in Escherichia coli of the genes encoding the major virulence factors, the genes defining the O157:H7 serotype, and components of the type 2 Shiga toxin family by multiplex PCR. *J Clin Microbiol.* 2002; 40: 3613–3619. https://doi.org/10.1128/JCM.40.10.3613-3619.2002 PMID: 12354854

33. Muniesa M, de Simon M, Prats G, Ferrer D, Pañella H, Jofre J. Shiga toxin 2-converting bacteriophages associated with clonal variability in Escherichia coli O157:H7 strains of human origin isolated from a single outbreak. *Infect Immun.* 2003; 71: 4554–4562. https://doi.org/10.1128/IAI.71.8.4554-4562.2003 PMID: 12874335

34. Braun M, Stuber K, Schlatter Y, Wahli T, Kuhnert P, Frey J. Characterization of an ADP-ribosyltransferase toxin (AexT) from Aeromonas salmonicida subsp. salmonicida. *J Bacteriol.* 2002; 184: 1851–1858. https://doi.org/10.1128/JB.184.7.1851-1858.2002 PMID: 11889090

35. Merino S, Gavin R, Vilches S, Shaw JG, Tomás JM. A colonization factor (production of lateral flagella) of mesophilic Aeromonas spp. is inactive in Aeromonas salmonicida strains. *Appl Environ Microbiol.* 2003; 69: 663–667. https://doi.org/10.1128/AEM.69.1.663-667.2003 PMID: 12514057

36. Roger F, Lamy B, Jumas-Bilak E, Kodjo A, colBVH study group, Marchandin H. Ribosomal multi-operon diversity: an original perspective on the genus Aeromonas. *PloS One.* 2012; 7: e46268. https://doi.org/10.1371/journal.pone.0046268 PMID: 23032081

37. Silver AC, Graf J. Prevalence of genes encoding the type three secretion system and the effectors AexT and AexU in the Aeromonas veronii group. *DNA Cell Biol.* 2009; 28: 383–388. https://doi.org/10.1089/dna.2009.0867 PMID: 19534604

38. Wang G, Tyler KD, Munro CK, Johnson WM. Characterization of cytotoxic, hemolytic Aeromonas caviae clinical isolates and their identification by determining presence of a unique hemolysin gene. *J Clin Microbiol.* 1996; 34: 3203–3205. PMID: 8940472

39. Li F, Wang W, Zhu Z, Chen A, Du P, Wang R, et al. Distribution, virulence-associated genes and antimicrobial resistance of Aeromonas isolates from diarrheal patients and water, China. *J Infect.* 2015; 70: 600–608. https://doi.org/10.1016/j.jinf.2014.11.004 PMID: 25447712

40. Chen P-L, Wu C-J, Chen C-S, Tsai P-J, Tang H-J, Ko W-C. A comparative study of clinical Aeromonas dhakensis and Aeromonas hydrophila isolates in southern Taiwan: A. dhakensis is more predominant and virulent. *Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis.* 2013; https://doi.org/10.1111/1469-0691.12456 PMID: 24237662

41. Alperi A, Figueras MJ. Human isolates of Aeromonas possess Shiga toxin genes (stx1 and stx2) highly similar to the most virulent gene variants of Escherichia coli. *Clin Microbiol Infect.* 2010; 16: 1563–1567. https://doi.org/10.1111/j.1469-0691.2010.03203.x PMID: 20219084

42. Buckley JT, Howard SP. The Cytotoxic Enterotoxin of Aeromonas hydrophila Is Aerolysin. *Infect Immun.* 1999; 67: 466–467. PMID: 9925450

43. Bernheimer AW, Avigad LS. Partial characterization of aerolysin, a lytic exotoxin from Aeromonas hydrophila. *Infect Immun.* 1974; 9: 1016–1021. PMID: 4208526
44. Chakraborty T, Huhle B, Bergbauer H, Goebel W. Cloning, expression, and mapping of the Aeromonas hydrophila aerolysin gene determinant in Escherichia coli K-12. J Bacteriol. 1986; 167: 368–374. PMID: 3522552

45. Rose JM, Houston CW, Coppenhaver DH, Dixon JD, Kurosky A. Purification and chemical characterization of a cholera toxin-cross-reactive cytolytic enterotoxin produced by a human isolate of Aeromonas hydrophila. Infect Immun. 1989; 57: 1165–1169. PMID: 2925244

46. Ferguson MR, Xu XJ, Houston CW, Peterson JW, Coppenhaver DH, Popov VL, et al. Hyperproduction, purification, and mechanism of action of the cytotoxic enterotoxin produced by Aeromonas hydrophila. Infect Immun. 1997; 65: 4299–4308. PMID: 9317040

47. Bücker R, Krug SM, Rosenthal R, Günzel D, Fromm A, Zeitz M, et al. Aerolysin from Aeromonas hydrophila perturbs tight junction integrity and cell lesion repair in intestinal epithelial HT-29/B6 cells. J Infect Dis. 2011; 204: 1283–1292. https://doi.org/10.1093/infdis/jir504 PMID: 21917902

48. Wilmsen HU, Pattus F, Buckley JT. Aerolysin, a hemolysin from Aeromonas hydrophila, forms voltage-gated channels in planar lipid bilayers. J Membr Biol. 1990; 115: 71–81. PMID: 1692344

49. Chen P-L, Lamy B, Ko W-C. Aeromonas dhakensis, an Increasingly Recognized Human Pathogen. Front Microbiol. 2016; 7: 793. https://doi.org/10.3389/fmicb.2016.00793 PMID: 27303382

50. Balsalobre LC, Dropa M, Matté GR, Matté MH. Molecular detection of enterotoxins in environmental strains of Aeromonas hydrophila and Aeromonas jandaei. J Water Health. 2009; 7: 685–691. https://doi.org/10.2166/wh.2009.082 PMID: 19590136

51. Chopra AK, Peterson JW, Xu XJ, Coppenhaver DH, Houston CW. Molecular and biochemical characterization of a heat-labile cytotonic enterotoxin from Aeromonas hydrophila. Microb Pathog. 1996; 21: 357–377. https://doi.org/10.1006/mpat.1996.0068 PMID: 8938643

52. Merino S, Aguilar A, Nogueras MM, Regue M, Swift S, Tomás JM. Cloning, sequencing, and role in virulence of two phospholipases (A1 and C) from mesophilic Aeromonas sp. serogroup O:34. Infect Immun. 1999; 67: 4008–4013. PMID: 10417167

53. Silver AC, Kikuchi Y, Fadl AA, Sha J, Chopra AK, Graf J. Interaction between innate immune cells and a bacterial type III secretion system in mutualistic and pathogenic associations. Proc Natl Acad Sci U S A. 2007; 104: 9481–9486. https://doi.org/10.1073/pnas.0700286104 PMID: 17517651

54. Abrami L, Fivaz M, Decroly E, Seidah NG, Jean F, Thomas G, et al. The pore-forming toxin proaerolysin is activated by furin. J Biol Chem. 1998; 273: 32656–32661. PMID: 9830006

55. Burr SE, Stuber K, Frey J. The ADP-ribosylating toxin, AexT, from Aeromonas salmonicida subsp. salmonicida is translocated via a type III secretion pathway. J Bacteriol. 2003; 185: 6583–6591. https://doi.org/10.1128/JB.185.22.6583-6591.2003 PMID: 14594831

56. Sha J, Wang SF, Suarez G, Sierra JC, Fadl AA, Erova TE, et al. Further characterization of a type III secretion system (T3SS) and of a new effector protein from a clinical isolate of Aeromonas hydrophila— part I. Microb Pathog. 2007; 43: 127–146. https://doi.org/10.1016/j.micpath.2007.05.002 PMID: 17644303