# A Phylogenetic Re-Analysis of Groupers with Applications for Ciguatera Fish Poisoning

Charlotte Schoelinck, Damien D Hinsinger, Agnès Dettaï, Corinne Cruaud, Jean-Lou Justine

► To cite this version:

Charlotte Schoelinck, Damien D Hinsinger, Agnès Dettaï, Corinne Cruaud, Jean-Lou Justine. A Phylogenetic Re-Analysis of Groupers with Applications for Ciguatera Fish Poisoning. PLoS ONE, 2014, 9 (8), pp.e98198. 10.1371/journal.pone.0098198. cea-04646651

HAL Id: cea-04646651

https://cea.hal.science/cea-04646651v1

Submitted on 12 Jul 2024

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
# A Phylogenetic Re-Analysis of Groupers with Applications for Ciguatera Fish Poisoning

Charlotte Schoelinck¹,²*, Damien D. Hinsinger¹, Agnès Dettai¹, Corinne Cruaud³, Jean-Lou Justine¹

1 UMR 7138 "Systématique, Adaptation, Évolution", Muséum National d'Histoire Naturelle, Département Systématique et Évolution, Paris, France, 2 Fisheries and Oceans Canada, Molecular biology, Aquatic animal health, Moncton, Canada, 3 Génoscope, Centre National de Séquençage, Évry, France

## Abstract

**Background:** Ciguatera fish poisoning (CFP) is a significant public health problem due to dinoflagellates. It is responsible for one of the highest reported incidence of seafood-borne illness and Groupers are commonly reported as a source of CFP due to their position in the food chain. With the role of recent climate change on harmful algal blooms, CFP cases might become more frequent and more geographically widespread. Since there is no appropriate treatment for CFP, the most efficient solution is to regulate fish consumption. Such a strategy can only work if the fish sold are correctly identified, and it has been repeatedly shown that misidentifications and species substitutions occur in fish markets.

**Methods:** We provide here both a DNA-barcoding reference for groupers, and a new phylogenetic reconstruction based on five genes and a comprehensive taxonomical sampling. We analyse the correlation between geographic range of species and their susceptibility to ciguatera accumulation, and the co-occurrence of ciguatoxins in closely related species, using both character mapping and statistical methods.

**Results:** Misidentifications were encountered in public databases, precluding accurate species identifications. Epinephelinae now includes only twelve genera (vs. 15 previously). Comparisons with the ciguatera incidences show that in some genera most species are ciguateric, but statistical tests display only a moderate correlation with the phylogeny. Atlantic species were rarely contaminated, with ciguatera occurrences being restricted to the South Pacific.

**Conclusions:** The recent changes in classification based on the reanalyses of the relationships within Epinephelidae have an impact on the interpretation of the ciguatera distribution in the genera. In this context and to improve the monitoring of fish trade and safety, we need to obtain extensive data on contamination at the species level. Accurate species identifications through DNA barcoding are thus an essential tool in controlling CFP since meal remnants in CFP cases can be easily identified with molecular tools.

**Citation:** Schoelinck C, Hinsinger DD, Dettai A, Cruaud C, Justine J-L (2014) A Phylogenetic Re-Analysis of Groupers with Applications for Ciguatera Fish Poisoning. PLoS ONE 9(8): e98198. doi:10.1371/journal.pone.0098198

**Editor:** James P. Meador, Northwest Fisheries Science Center, NOAA Fisheries, United States of America

**Received** September 18, 2013; **Accepted** April 29, 2014; **Published** August 5, 2014

**Copyright:** © 2014 Schoelinck et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This study was funded by MNHN ATM Barcode (2010, 2011), MNHN ATM Biodiversité Actuelle et Fossile (2010) and MNHN BQR (2011), awarded to CS and JLJ. This work was supported by the "Consortium National de Recherche en Génomique", and the "Service de Systématique Moléculaire" of the Muséum National d'Histoire Naturelle (CNRS UMS 2700). It is part of the agreement n°2005/67 between the Genoscope and the Muséum National d'Histoire Naturelle on the project "Macrophylogeny of life" directed by Guillaume Lecointre. This work is part of the project @ SPEED-ID "Accurate SPECiEs Delimitation and iDentification of eukaryotic biodiversity using DNA markers" proposed by F-BoL, the French Barcode of life initiative. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* Email: schoelinck@mnhn.fr

## Introduction

Large carnivorous fishes associated with coral reefs are frequently contaminated by toxins responsible for ciguatera fish poisoning (CFP) in tropical and subtropical waters [1,2]. CFP is a food-borne disease contracted by the consumption of finfish that have accumulated lipid-soluble toxins produced by microalgae (dinoflagellates) of the genus *Gambierdiscus* in their flesh and viscera. Dinoflagellates produce gambiertoxins which are first accumulated in the viscera of herbivorous fish and are further accumulated and converted to ciguatoxins in the flesh of larger carnivorous species. For the purposes of this report, we define ciguateric as possessing the ability to accumulate ciguatoxins and cause ciguatera fish poisoning. At least three groups of ciguatoxins have been identified: Pacific (P-CTX), Indian Ocean (I-CTX) and

Caribbean (C-CTX) [3,4]. While the gambiertoxin precursors for P-CTX have been identified, the corresponding precursors for I-CTX and C-CTX have yet to be identified, let alone a thorough examination of which dinoflagellates produce them. This disease produces several gastrointestinal, neurological and cardiac symptoms a few minutes to a few hours after ingestion of contaminated seafood [5]. Although there are reports of symptom amelioration with some interventions (e.g. IV mannitol), no efficient treatment exists so far [6]. It is a significant public health problem, especially in the South Pacific but also in the United States, where it is responsible for one of the highest reported incidence of seafood-borne illness [1]. Although CFP was historically restricted to tropical and sub-tropical regions, case reports are increasingly seen in higher latitudes with escalating global trade and movement of seafood products [1,7]. The incidence of ciguatera, as well as the
species of fish that are potentially poisonous, vary from region to region [2]. Precise information about the distribution of ciguatera-carrying species can be obtained from epidemiological data collected by research and health organisations in each country or region, but this depends heavily on correct species identification and is highly dependent on the intensity of data collection. An additional problem pointed out by several authors is the role of recent climate change on harmful algal blooms (HAB), including *Gambierdiscus* ssp. [1,8,9]. The abundance of *G.* spp. correlates positively with elevated sea surface temperature [8]. CFP cases might therefore become more frequent and more geographically widespread as an indirect consequence of climate change (review in [1]). Moreover, coral reefs perturbations, such as hurricanes or bleaching events, also free up space for microalgae to colonize. Even human activities altering the environment such as petroleum production platform building can contribute to the HAB [10]. Therefore, populations from developing countries, already facing these disturbances, appeared to be particularly exposed to the intensification of CFP.

Since there is no appropriate treatment for CFP (for a review, see [6]), the most efficient solution is to regulate fish consumption [11,12]. Lewis [13], and more recently Clua et al. [11] recommended banning some specific species and sizes from fish markets. However, such a strategy can only work if the fish sold are correctly identified and labelled, and it has been repeatedly shown that misidentifications and species substitutions commonly occur in fish markets [14–16].

Groupers (Epinephelidae: rockcods, coralgroupers, hinds, and lyretails) are one of the families most commonly reported as a source of ciguatera poisoning [11]. Some grouper species, like *Plectropomus laevis* and *Cephalopholis argus*, are known to be especially contaminated by ciguatera toxins [17–19]. Large individuals are generally more toxic than small ones since ciguatoxins accumulate in fish via the food chain [11,20]. For instance, specimens belonging to the potentially ciguatoxic fish species *Epinephelus fuscoguttatus* and *Variola louti* are considered dangerous only if they weigh more than 13 and 1.7 kg respectively [20]. Because they are widely distributed in warm and temperate shore waters, from surface to deep-sea, and adults of some species reach 3 m in length and 400 kg [18], groupers represent a considerable economic value in tropical and subtropical regions and most particularly in south-east Asia [21–24]. They are a major component of the artisanal fisheries resource especially in the south Pacific [18]. Global capture fisheries production has increased from approximately 214,000 tons in 1999 to more than 275,000 tons in 2009 [25]. Grouper aquaculture was first introduced in the early 1970s and is now widely practised throughout Southeast Asia [21]. Global grouper aquaculture production has increased tremendously due to increasing demand, from 60,000 tons in 1990 to 200,000 tons in 2007. The premium price of groupers can reach US$ 100/kg in the Chinese live fish markets [21].

Although groupers are large fish and supposedly easily identifiable, comprehensive and reliable species identification tools are rare and a good taxonomic framework is also necessary. Even when intact adult specimens are available (which is generally not the case for food-borne poisoning cases) the morphological characters used to discriminate species can be subtle, making identification difficult even for trained taxonomists. Moreover, accessing the historical literature and assessing the validity of species with a controversial taxonomic history are challenging tasks, even for experts [26].

Some rapid and reliable species identification tools such as DNA barcoding have been developed to facilitate species identification

[26–30]. Given the estimated $US200 billion annual value of fisheries worldwide, the Fish Barcode of Life campaign (FISH-BOL) initiative, as a part of the International Barcode of Life Project (iBOL; http://www.ibolproject.org), is addressing socially relevant questions concerning market substitution and quota management of commercial fisheries (http://www.fishbol.org), with a special focus on developing countries [31]. However, species identification tools require complete and reliable databases. Indeed, DNA databases play a key role for the species identification of groupers, and more generally, for seafood, as non-specialists use essentially those databases to identify species for which they often have access to tissue samples only. The Epinephelidae comprise about 163 species [32] among which 106 are recorded in BOLD (942 public sequences in February 2014). The incompleteness of the reference datasets is a well-identified problem for species identification [33], which can be slightly alleviated if the marker used for identification is also relevant for phylogeny. In such a case, and if the taxonomic framework is accurate, species not represented in the database might still be assigned to clades or higher rank groups, like genera. Completing the largest molecular identification dataset (the cytochrome oxidase 1 of the Barcode of Life project), combined with an accurate study of the relationships of groupers, will help the management of grouper diversity through easier and more accurate identifications. Much remains to be done on both of these aspects, as the phylogenetic framework of the group has undergone many changes recently, and is yet incomplete.

The relationships of the Epinephelidae, recently raised to family rank by Smith and Craig [34], are indeed not yet totally resolved. Epinephelidae were previously a subfamily (Epinephelinae) included with Serraninae and Anthiinae among Serranidae [18,35]. The relationships within the former Serranidae, as well as the composition of the family, have been the object of much discussion. Two molecular studies including the Serranidae showed the non-monophyly of the family [34,36]. Smith and Craig [34] grouped Serraninae and Anthiinae in the Serranidae and raised the subfamily Epinephelinae to the family rank Epinephelidae. On the other hand, Lautrédou et al. [36] showed the polyphyly of Serranidae (with the Serraninae – Anthiinae composition), while recovering an Anthiinae and Epinephelidae clade. Craig et al. [32] defined four subfamilies in Epinephelidae: Diploprioninae, Epinephelinae, Grammistinae and Liopropominae, corresponding to the four previous tribes Diploprionini, Epinephelini, Grammistini, Liopropomini. In their molecular phylogeny, Craig and Hastings [37] attempted to resolve the phylogeny of the Epinephelidae using an almost complete species sampling within the genus *Epinephelus*, and several specimens of other subfamilies, using two mitochondrial and two nuclear markers. They proposed taxonomic changes for species of the subfamily Epinephelinae to reflect their phylogenetic position. For instance, they included *Cromileptes altivelis* and *Anyperodon leucogrammicus* in *Epinephelus* and they moved *Epinephelus septemfasciatus* and *E. ergastularius* to *Hyporthodus*. However, many nodes of their phylogeny lacked robustness. Because of the absence of morphological differences between the genera *Anyperodon*, *Cromileptes* and *Epinephelus*, Craig et al. [32] retained the monotypic genera *Anyperodon* and *Cromileptes*.

To further the study of the relationships between the genera, we sequenced five markers, two mitochondrial, Cytochrome Oxidase Subunit I (COI) and 16S ribosomal RNA (16S) and three nuclear, Rhodopsin (Rh), Titin-like protein (TMO-4C4) and Polycystic kidney disease 1 protein (Pkd1). We choose to include the reference barcoding marker COI to provide a simple and reliable tool for species identification to the non-specialist community
<table><caption>Table 1. List of fish species, specimen vouchers localities and sequence accession numbers.</caption><thead><tr><th rowspan="2">Family</th><th rowspan="2">Group: subfamily /tribe</th><th rowspan="2">Species</th><th colspan="4">Accession Number</th><th colspan="2">Specimen</th><th rowspan="2">Locality</th></tr><tr><th colspan="2">Mitochondrial</th><th colspan="2">Nuclear</th><th colspan="2">Voucher</th></tr><tr><th></th><th></th><th></th><th>COI</th><th>16S</th><th>TMO-4C</th><th>Rhodopsin</th><th>Pkd1</th><th></th><th></th><th></th></tr></thead><tbody><tr><td rowspan="12">Epinephelidae</td><td rowspan="12">Diplopriorninae</td><td>Belonoprica chabanaudi</td><td>JQ431484,<sup>*,v</sup></td><td>JX094024,<sup>*,v</sup></td><td>JX093971,<sup>*,v</sup></td><td></td><td></td><td>MNHN 2008-1159</td><td>Moorea, French Polynesia</td></tr><tr><td>Diploprion bifasciatum</td><td>KM077912,<sup>*,v</sup></td><td>KM077970,<sup>*,v</sup></td><td>KM078001,<sup>*,v</sup></td><td></td><td></td><td>MNHN-icti-2815</td><td>Queensland, Australia</td></tr><tr><td rowspan="5">Liopropominae</td><td>JX093903*</td><td>JX093999*</td><td>JX093972*</td><td>JX093952*</td><td>JX093928*</td><td></td><td>Ecuador</td></tr><tr><td>JX0931888,<sup>v</sup></td><td>JX094023,<sup>*,v</sup></td><td>JX093974,<sup>*,v</sup></td><td>JX093953,<sup>*,v</sup></td><td>JX093929,<sup>*,v</sup></td><td></td><td>MNHN 2008-1023</td><td>Moorea, French Polynesia</td></tr><tr><td>Liopropoma lunulatum</td><td>JQ431890</td><td>JX094020,<sup>*,v</sup></td><td>JX093973,<sup>*,v</sup></td><td></td><td></td><td>MNHN 2008-0793</td><td>Moorea, French Polynesia</td></tr><tr><td>Liopropoma pallidum</td><td>JQ431457</td><td>JX094016,<sup>*,v</sup></td><td>AY949271</td><td></td><td></td><td>MNHN 2008-0307</td><td>Moorea, French Polynesia</td></tr><tr><td>Aporops bilinearis</td><td>JQ431457,<sup>v</sup></td><td>AY539050</td><td>AY539458</td><td></td><td></td><td>MNHN 2008-1105</td><td>Moorea, French Polynesia</td></tr><tr><td rowspan="3">Grammistinae</td><td>Grammistes sexlineatus</td><td>JQ431776,<sup>v</sup></td><td>AY539050</td><td>AY539458</td><td></td><td></td><td>USNM 391102</td><td>Moorea, French Polynesia</td></tr><tr><td>Grammistops ocellatus</td><td>JQ431778,<sup>v</sup></td><td>JX094021,<sup>*,v</sup></td><td>X093975,<sup>*,v</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>Pogonoprica punctata</td><td>JX093904*</td><td>JX093998*</td><td>JX093976*</td><td>JX093951*</td><td>JX093927*</td><td></td><td></td></tr><tr><td rowspan="2">Pseudogramma gregoryi</td><td rowspan="2">GU225013</td><td rowspan="2">AY947571</td><td rowspan="2">AY949213</td><td></td><td></td><td>MBIO509.4</td><td>Moorea, French Polynesia</td></tr><tr><td></td><td></td><td></td><td></td><td>Ghana</td></tr><tr><td rowspan="2">Pseudogramma polyanantha</td><td>JQ432063,<sup>v</sup></td><td>JX094018,<sup>*,v</sup></td><td>JX093977,<sup>*,v</sup></td><td></td><td></td><td>MNHN2002-158</td><td></td></tr><tr><td>JX093905,<sup>*,v</sup></td><td>JX094000,<sup>*,v</sup></td><td>JX093978,<sup>*,v</sup></td><td>JX093956,<sup>*,v</sup></td><td>JX093932,<sup>*,v</sup></td><td>MNHN2002-158</td><td></td></tr><tr><td rowspan="2">Rypticus saponaceus</td><td rowspan="2">JQ432178</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="6">Epinephelinae</td><td>Cephalopholis argus</td><td>JQ431565,<sup>v</sup></td><td>JX094015,<sup>*,v</sup></td><td>JX093979,<sup>*,v</sup></td><td>JX093958,<sup>*,v</sup></td><td>JX093934,<sup>*,v</sup></td><td>MNHN 2008-0229</td><td>Moorea, French Polynesia</td></tr><tr><td>Cephalopholis boenak</td><td>KM077907,<sup>*,v</sup></td><td>KM077965,<sup>*,v</sup></td><td>KM077996,<sup>*,v</sup></td><td>KM077936,<sup>*,v</sup></td><td>KM077879,<sup>*,v</sup></td><td>MNHN-icti-2875</td><td>New Caledonia</td></tr><tr><td>Cephalopholis colonus</td><td>GU440449</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cephalopholis cruentata</td><td>GU225172</td><td>AF297323</td><td>AY949266</td><td></td><td></td><td>MNHN-icti-2821</td><td>Queensland, Australia</td></tr><tr><td>Cephalopholis cyanostigma</td><td>KM077908,<sup>*,v</sup></td><td>KM077966,<sup>*,v</sup></td><td>KM077997,<sup>*,v</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>Cephalopholis formosa</td><td>FJ583004</td><td>AY947603</td><td>EF517741</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Cephalopholis fulva</td><td rowspan="2">FJ583007</td><td rowspan="2">AF297292</td><td rowspan="2">AY949282</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Cephalopholis hemistiktos</td><td>HQ149822</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cephalopholis igarashiensis</td><td>EU871685</td><td>AY947599</td><td>AY949292</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Cephalopholis leopardus</td><td>FJ583010</td><td>AY947560</td><td>AY949323</td><td></td><td></td><td></td><td></td></tr><tr><td>Cephalopholis micropion</td><td>FJ237608</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Cephalopholis miniata</td><td>KM077909*</td><td>KM077967*</td><td>KM077998,<sup>*,v</sup></td><td>KM077937,<sup>*,v</sup></td><td>KM077880,<sup>*,v</sup></td><td>MNHN-icti-3008</td><td>New Caledonia</td></tr><tr><td>Cephalopholis rogaa</td><td>JQ349677</td><td>EF503626</td><td>EF517737</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Cephalopholis sexmaculata</td><td>JQ431572,<sup>v</sup></td><td>JX094019,<sup>*,v</sup></td><td>JX093980,<sup>*,v</sup></td><td>JX093959,<sup>*,v</sup></td><td>JX093935,<sup>*,v</sup></td><td>MNHN 2008-0754</td><td>Moorea, French Polynesia</td></tr><tr><td>Cephalopholis sonnerati</td><td>JX093918*</td><td>JX094007*</td><td>JX093981*</td><td>JX093960*</td><td>JX093936*</td><td></td><td>New Caledonia</td></tr><tr><td rowspan="2">Cephalopholis spilopararea</td><td>KM077910,<sup>*,v</sup></td><td>KM077968,<sup>*,v</sup></td><td>KM077999,<sup>*,v</sup></td><td>KM077938,<sup>*,v</sup></td><td>KM077881,<sup>*,v</sup></td><td>MNHN-icti-2961</td><td>New Caledonia</td></tr><tr><td>Cephalopholis urodeta</td><td>KM077911,<sup>*,v</sup></td><td>KM077969,<sup>*,v</sup></td><td>KM078000,<sup>*,v</sup></td><td>KM077939,<sup>*,v</sup></td><td>KM077882,<sup>*,v</sup></td><td>MNHN-icti-2886</td><td>New Caledonia</td></tr><tr><td rowspan="2">Dermatolepis dermatolepis</td><td>JX093917*</td><td>JX094006*</td><td>JX093982*</td><td>JX093955*</td><td>JX093931*</td><td></td><td>Ecuador</td></tr><tr><td>Epinephelus adscensionis</td><td>FJ583396</td><td>AY539049</td><td>AY949284</td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 1. Cont.</caption><thead><tr><th rowspan="2">Family</th><th rowspan="2">Group: subfamily /tribe</th><th rowspan="2">Species</th><th colspan="2">Accession Number</th><th colspan="2">Specimen</th><th rowspan="2">Locality</th></tr><tr><th>COI</th><th>16S</th><th>Nuclear</th><th>Voucher</th></tr></thead><tbody><tr><td></td><td></td><td>Epinephelus aeneus</td><td>KM077913<sub>v</sub>*</td><td>KM077971<sub>v</sub>*</td><td>KM078002<sub>v</sub>*</td><td>KM077940<sub>v</sub>*</td><td>KM077883<sub>v</sub>*</td><td>MNHN-icti-2842</td><td>Senegal<sup>1</sup></td></tr><tr><td></td><td></td><td>Epinephelus akaara</td><td>GBGC7768</td><td>DQ154107</td><td>EF517707</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus altivelis</td><td>JX093906*</td><td>JX094001*</td><td>JX093983*</td><td>JX093957*</td><td>JX093933*</td><td></td><td>Origin unknown, Aquarium</td></tr><tr><td></td><td></td><td>Epinephelus amblycephalus</td><td>JX093910*</td><td>JX094009*</td><td>JX093984*</td><td>JX093961*</td><td>JX093937*</td><td></td><td>South China Sea<sup>1</sup></td></tr><tr><td></td><td></td><td>Epinephelus analogus</td><td>JX093915*</td><td>JX094003*</td><td>AY949220</td><td></td><td></td><td></td><td>Ecuador</td></tr><tr><td></td><td></td><td>Epinephelus areolatus</td><td>KM077914<sub>v</sub>*</td><td>KM077972<sub>v</sub>*</td><td>KM078003<sub>v</sub>*</td><td>KM077941<sub>v</sub>*</td><td>KM077884<sub>v</sub>*</td><td>MNHN-icti-2996</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus awoara</td><td>JX093913*</td><td>JX094013*</td><td>JX093985*</td><td>JX093969*</td><td>JX093945*</td><td></td><td>South China Sea<sup>1</sup></td></tr><tr><td></td><td></td><td>Epinephelus bleekeri</td><td>JX093911*</td><td>JX094011*</td><td>JX093986*</td><td>JX093962*</td><td>JX093938*</td><td></td><td>South China Sea<sup>1</sup></td></tr><tr><td></td><td></td><td>Epinephelus brunus</td><td>JX093909*</td><td>JX094010*</td><td>JX093992*</td><td>JX093968*</td><td>JX093944*</td><td></td><td>South China Sea<sup>1</sup></td></tr><tr><td></td><td></td><td>Epinephelus chlorostigma</td><td>JQ412501<sub>v</sub></td><td>KM077973<sub>v</sub>*</td><td>KM078004<sub>v</sub>*</td><td>KM077942<sub>v</sub>*</td><td>KM077885<sub>v</sub>*</td><td>MNHN-icti-2145</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus clippertonensis</td><td>JX093914*</td><td>JX094002*</td><td>JX093987*</td><td>JX093964*</td><td>JX093940*</td><td></td><td>Clipperton Island</td></tr><tr><td></td><td></td><td>Epinephelus coeruleopunctatus</td><td>KM077915*,<sub>v</sub></td><td>KM077974*,<sub>v</sub></td><td>KM078005*,<sub>v</sub></td><td>KM077943*,<sub>v</sub></td><td>KM077886*,<sub>v</sub></td><td>MNHN-2006-1706</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus coioides</td><td>KM077916*,<sub>v</sub></td><td>KM077975*,<sub>v</sub></td><td>KM078006*,<sub>v</sub></td><td>KM077944*,<sub>v</sub></td><td>KM077887*,<sub>v</sub></td><td>MNHN-icti-2994</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus corallicola</td><td>JX093908*</td><td>JX094008*</td><td>JX093988*</td><td>JX093967*</td><td>JX093943*</td><td></td><td>South China Sea<sup>1</sup></td></tr><tr><td></td><td></td><td>Epinephelus cyanopodus</td><td>JQ412502<sub>v</sub></td><td>KM077976*,<sub>v</sub></td><td>KM078007*,<sub>v</sub></td><td>KM077945*,<sub>v</sub></td><td>KM077888*,<sub>v</sub></td><td>MNHN-icti-2148</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus diacanthus</td><td>EF609517</td><td>AY947619</td><td>AY949274</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus fasciatomaculosus</td><td>EF607565</td><td>AY947622</td><td>EF517717</td><td></td><td></td><td></td><td>Gulf of Aqaba</td></tr><tr><td></td><td></td><td>Epinephelus fasciatatus</td><td>JX093907*</td><td>JX094005*</td><td>JX093990*</td><td>JX093970*</td><td>JX093946*</td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus flavocaeruleus</td><td>SAIAB330-06<sub>b</sub></td><td>AY947607</td><td>EF517731</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus fuscoguttatus</td><td>EU600139</td><td>JF750752</td><td>EF517713</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus hexagonatus</td><td>JQ431719<sub>v</sub></td><td>JX094017*,<sub>v</sub></td><td>JX093991*,<sub>v</sub></td><td>JX093966*,<sub>v</sub></td><td>JX093942*,<sub>v</sub></td><td>MNHN 2008-0338</td><td>Moorea, French Polynesia</td></tr><tr><td></td><td></td><td>Epinephelus howlandi</td><td>KM077917*,<sub>v</sub></td><td>KM077977*,<sub>v</sub></td><td>KM078008*,<sub>v</sub></td><td>KM077946*,<sub>v</sub></td><td>KM077889*,<sub>v</sub></td><td>MNHN-icti-2981</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus lanceolatus</td><td>NC_011715</td><td>AY947588</td><td>EF517736</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus latifasciatus</td><td>EU014219</td><td>DQ088044</td><td>EF517724</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus leucogrammicus</td><td>KM077918*,<sub>v</sub></td><td>KM077978*,<sub>v</sub></td><td>KM078009*,<sub>v</sub></td><td>KM077947*,<sub>v</sub></td><td>KM077890*,<sub>v</sub></td><td>MNHN-icti-2993</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus longispinis</td><td>EF609522</td><td>EF213704</td><td>EF517706</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus macrospilos</td><td>SAIAB547-07<sub>b</sub></td><td>AY731072</td><td>AY949238</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus maculatus</td><td>KM077919*,<sub>v</sub></td><td>KM077979*,<sub>v</sub></td><td>KM078010*,<sub>v</sub></td><td>KM077948*,<sub>v</sub></td><td>KM077891*,<sub>v</sub></td><td>MNHN-icti-2907</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus malabaricus</td><td>KM077920*,<sub>v</sub></td><td>KM077980*,<sub>v</sub></td><td>KM078011*,<sub>v</sub></td><td>KM077949*,<sub>v</sub></td><td>KM077892*,<sub>v</sub></td><td>MNHN-icti-3018</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus melanostigma</td><td>JQ349966</td><td>EF503624</td><td>EF517730</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus merra</td><td>KM077921*,<sub>v</sub></td><td>KM077981*,<sub>v</sub></td><td>KM078012*,<sub>v</sub></td><td>KM077950*,<sub>v</sub></td><td>KM077893*,<sub>v</sub></td><td>MNHN-icti-2963</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus multinotatus</td><td>SAIAB149-06<sub>b</sub></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 1. Cont.</caption><thead><tr><th rowspan="2">Family</th><th rowspan="2">Group: subfamily /tribe</th><th rowspan="2">Species</th><th colspan="2">Accession Number</th><th colspan="2">Specimen</th><th rowspan="2">Locality</th></tr><tr><th>Mitochondrial</th><th>16S</th><th>Nuclear</th><th>Pkd1</th></tr></thead><tbody><tr><td></td><td></td><td>Epinephelus ongus</td><td>KM077922<sub>v</sub>*</td><td>KM077982<sub>v</sub>*</td><td>KM078013<sub>v</sub>*</td><td>KM077951<sub>v</sub>*</td><td>KM077894<sub>v</sub>*</td><td>MNHN-icti-3030</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus polyphekadion</td><td>KM077923<sub>v</sub></td><td>KM077983<sub>v</sub></td><td>KM078014<sub>v</sub></td><td>KM077952<sub>v</sub>*</td><td>KM077895<sub>v</sub></td><td>MNHN-icti-2979</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus polyphekadion</td><td></td><td></td><td></td><td></td><td></td><td>MNHN-icti-3015</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus quoyanus</td><td>DQ107861</td><td>JX094014*</td><td>JX093993*</td><td></td><td></td><td>MNHN-icti-2819</td><td>Queensland, Australia</td></tr><tr><td></td><td></td><td>Epinephelus retouti</td><td>KM077924*</td><td>KM077984*</td><td>KM078015*</td><td>KM077954*</td><td>KM077896*</td><td>MNHN-icti-2967</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus rivulatus</td><td>KM077925<sub>v</sub>*</td><td>KM077985<sub>v</sub></td><td>KM078016<sub>v</sub></td><td>KM077955<sub>v</sub></td><td>KM077897<sub>v</sub></td><td>MNHN-icti-2890</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Epinephelus sexfasciatus</td><td>EF607564</td><td>DQ067310</td><td>EF517716</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Epinephelus spintoceps</td><td>KM077926<sub>v</sub>*</td><td>KM077986<sub>v</sub></td><td>KM078017<sub>v</sub>*</td><td>KM077956<sub>v</sub>*</td><td>KM077898<sub>v</sub>*</td><td>MNHN-icti-3050</td><td>Maldives<sup>1</sup></td></tr><tr><td></td><td></td><td>Epinephelus tauvina</td><td>JX093916*</td><td>JX094004*</td><td>JX093994*</td><td>JX093965*</td><td>JX093941*</td><td></td><td>Gulf of Aqaba</td></tr><tr><td></td><td></td><td>Epinephelus undulosus</td><td>EF609352</td><td>DQ088041</td><td>EF517704</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Hyporthodus haifensis</td><td>CSFOM034-10<sub>b</sub></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Hyporthodus septemfasciatus</td><td>DQ107851</td><td>AY947559</td><td>AY949247</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca bonaci</td><td>GU225646</td><td>DQ267150</td><td>AY949270</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca canina</td><td>CSFOM032-10<sub>b</sub></td><td>AY947585</td><td>AY949294</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca costae</td><td>KM077928<sub>v</sub>*</td><td>KM077988<sub>v</sub></td><td>KM078019<sub>v</sub></td><td>KM077957<sub>v</sub>*</td><td></td><td>MNHN-icti-2868</td><td>Tunisia<sup>1</sup></td></tr><tr><td></td><td></td><td>Mycateroperca costae</td><td></td><td></td><td></td><td></td><td>KM077899<sub>v</sub>*</td><td>MNHN-icti-2859</td><td>Tunisia<sup>1</sup></td></tr><tr><td></td><td></td><td>Mycateroperca interstitialis</td><td>FJ583668</td><td>AY947632</td><td>AY949221</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca jordani</td><td>GU440412</td><td>AF297329</td><td>AY949303</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca marginata</td><td>KM077929<sub>v</sub>*</td><td>KM077989<sub>v</sub></td><td>KM078020<sub>v</sub></td><td>KM077958<sub>v</sub></td><td>KM077900<sub>v</sub>*</td><td>MNHN-icti-2854</td><td>Senegal</td></tr><tr><td></td><td></td><td>Mycateroperca microlepis</td><td>JN021310</td><td>AF297312</td><td>AY949253</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca morrhua</td><td>KM077930<sub>v</sub>*</td><td>KM077990<sub>v</sub>*</td><td>KM078021<sub>v</sub>*</td><td>KM077959<sub>v</sub>*</td><td>KM077901<sub>v</sub>*</td><td>MNHN-icti-2947</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Mycateroperca poecilonota</td><td>SAJAB576-07<sub>b</sub></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca rubra</td><td>CSFOM051-10<sub>b</sub></td><td>AY947587</td><td>AY949255</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Mycateroperca xenarcha</td><td>GU440413</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Plectropomus areolatus</td><td>JN242591</td><td>EF213706</td><td>EF517750</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Plectropomus laevis</td><td>KM077932<sub>v</sub>*</td><td>KM077992<sub>v</sub>*</td><td>KM078023<sub>v</sub>*</td><td>KM077961<sub>v</sub>*</td><td>KM077903<sub>v</sub>*</td><td>MNHN-icti-3012</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Plectropomus leopardus</td><td>KM077933<sub>v</sub>*</td><td>KM077993<sub>v</sub></td><td>KM078024<sub>v</sub></td><td>KM077962<sub>v</sub></td><td>KM077904<sub>v</sub></td><td>MNHN-icti-2883</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Plectropomus maculatus</td><td>DQ107911</td><td>JF750755</td><td>EF517751</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Saloptia powelli</td><td>JQ432090<sub>v</sub></td><td>JX094022<sub>v</sub>*</td><td>JX093995<sub>v</sub></td><td>JX093954<sub>v</sub>*</td><td>JX093930<sub>v</sub>*</td><td>MNHN 2008-1014</td><td>Moorea, French Polynesia</td></tr><tr><td></td><td></td><td>Triso dermopterus</td><td>DQ107934</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Variola albimarginata</td><td>KM077934<sub>v</sub>*</td><td>KM077994<sub>v</sub></td><td>KM078025<sub>v</sub></td><td>KM077963<sub>v</sub></td><td>KM077905<sub>v</sub></td><td>MNHN-icti-2999</td><td>New Caledonia</td></tr><tr><td></td><td></td><td>Variola louti</td><td>KM077935<sub>v</sub>*</td><td>KM077995<sub>v</sub>*</td><td>KM078026<sub>v</sub></td><td>KM077964<sub>v</sub></td><td>KM077906<sub>v</sub></td><td>MNHN-icti-2964</td><td>New Caledonia</td></tr></tbody></table>
<table><caption>Table 1. Cont.</caption><thead><tr><th rowspan="2">Family</th><th rowspan="2">Group: subfamily /tribe</th><th rowspan="2">Species</th><th colspan="2">Accession Number</th><th colspan="2">Specimen</th><th rowspan="2">Locality</th></tr><tr><th>Mitochondrial</th><th>COI</th><th>Nuclear</th><th>Pkd1</th></tr></thead><tbody><tr><td rowspan="2">Anthiinae</td><td rowspan="2"></td><td>Pseudanthias hyspelosoma</td><td>JX093919*</td><td>JX094027*</td><td>JX093996*</td><td>JX093950*</td><td>JX093926*</td><td rowspan="2">New Caledonia¹</td></tr><tr><td>Pseudanthias pleuroataenia</td><td>JX093920*</td><td>JX094026*</td><td>AY949308*</td><td>JX093949*</td><td>JX093925*</td></tr><tr><td rowspan="4">Cirrhitidae</td><td rowspan="4"></td><td>Pseudanthias tuka</td><td>JX093921*</td><td>JX094025*</td><td>JX093997*</td><td>JX093948*</td><td>JX093924*</td><td></td></tr><tr><td>Cirrhitus pinnulatus</td><td>JQ431641</td><td></td><td></td><td>KC222240b</td><td>JX628396</td><td></td></tr><tr><td>Cirrhitus rivulatus</td><td></td><td>AY539059</td><td>AY539467</td><td></td><td></td><td></td></tr><tr><td>Harpagifer kerguelenensis</td><td>EATF605b</td><td>AY539063</td><td>AY539471</td><td>EATF605b</td><td>JQ688766</td><td></td></tr><tr><td>Harpigiferidae</td><td></td><td>Niphon spinosus</td><td>EF143386</td><td>AY947575</td><td>AY949210</td><td></td><td></td></tr><tr><td>Niphonidae</td><td></td><td>Perca fluviatilis</td><td>AP040-12b</td><td>GU018097</td><td>AP040-12b</td><td>JX628360</td><td></td></tr><tr><td rowspan="10">Scorpaenidae</td><td rowspan="10"></td><td>Perca flavescens</td><td></td><td></td><td>AY539463</td><td></td><td></td></tr><tr><td>Zanclorhynchus spinifer</td><td>AP139-12b</td><td>AY538999</td><td>AY539413</td><td>EU638021</td><td>JX628429</td><td></td></tr><tr><td>Cyclopterus lumpus</td><td>AP041-12b</td><td>AY539043</td><td>AY539451</td><td>AP041-12b</td><td>JX628359</td><td></td></tr><tr><td>Scorpaeninae</td><td>Pontinus longispinis</td><td>AY538982</td><td>AY539398</td><td></td><td></td><td></td></tr><tr><td>Pontinus macrocephalus</td><td>JX093922*</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Scarpaenopsis macrochir</td><td></td><td>AY538987</td><td>AY539402</td><td></td><td></td><td>Philippines</td></tr><tr><td>Scarpaenopsis possi</td><td>JQ432137</td><td></td><td></td><td>KC222244b</td><td>JX628403</td><td></td></tr><tr><td>Helicolenus dactylopterus</td><td>AP121-12b</td><td>AY538975</td><td>AY539391</td><td>AP121-12b</td><td>JX628383</td><td></td></tr><tr><td>Trachyscorpia cristulata</td><td>AP111-12b</td><td>AY538980</td><td>AY539396</td><td>AP111-12b</td><td>JX628374</td><td></td></tr><tr><td>Synanceinae</td><td>Centropistis striata</td><td>JQ432179</td><td>AY538995</td><td>AY539410</td><td>KC222242b</td><td>JX628405</td><td></td></tr><tr><td rowspan="2">Serraninae</td><td>Serranus tigrinus</td><td>HQ024935</td><td>AY072667</td><td>AY949216</td><td></td><td></td><td></td></tr><tr><td>Trachinus draco</td><td>FJ584106</td><td>AY072688</td><td>AY949259</td><td>AP104-12b</td><td>JX628367</td><td></td></tr></tbody></table>

<footnote>(1) Fish specimen collected at Nouméa fish market or obtained through colleagues, (*) new sequences, (‡) BOLD accession number, (‡) sequence corresponds to voucher indicated in Table.
doi:10.1371/journal.pone.0098198.t001</footnote>
(fisheries, governmental organisations, etc.). As ciguatera occurrence has not yet been studied with regard to the evolutionary relationships of the ciguateric fish replaced in their evolutionary context, the second aim of this paper is to map the high risk species for ciguatera fish poisoning into the phylogeny using published information about ciguatera-prone species. We statistically test whether high risk species are closely related and could therefore have inherited their susceptibility to ciguatera from their common ancestor.

## Methods

### Taxonomic sampling

**Sequencing fresh material collection.** Fishes collected from different localities (Table 1) in 2009–2011 were dead at the time we acquired them for study, having been commercially caught, and available for purchase at the Nouméa fish market. Each individual was morphologically identified according to Heemstra & Randall [18], measured, weighed, and photographed, and a tissue sample was collected and preserved in absolute ethanol until DNA extraction. Several specimens per species were sequenced to evaluate intraspecific variation and to corroborate identification (data not shown). Additional fish tissues were obtained from colleagues (see acknowledgements) or bought at fish markets, and also preserved in absolute ethanol. For these tissues, no photograph was available; the identification of these tissues was checked by a BLAST search in BOLD [38] followed by a thorough evaluation of the results. Samples and results not corresponding to a higher sequence and identification quality standard were discarded.

**Publicly available sequences.** All available COI, 16S and TMO-4C4 sequences of Epinephelidae and some outgroup sequences were downloaded from public sequence databases (Barcode of life Database, GenBank Nucleotide). All sequences were controlled for contamination, indels, and stop-codons indicating possible pseudogenes [39]. We followed the classification of Craig and Hastings [37]. When necessary, we amended the species name to agree with the genus gender (*Cephalopholis rogaa*, *Mycteroperca canina* and *M. marginata*) as per the recommendations of the International Code of Zoological Nomenclature (Fourth Edition).

Problems were identified within both GenBank and BOLD sequences, and some sequences available in those databases were therefore not integrated in our dataset. Database sequences presented either: (i) taxonomic problems, such as high genetic

divergence within species or erroneous species identification and (ii) nomenclatural problems like the use of invalid species names.

The specimen identification for some sequences was problematic. While there were often several specimens attributed to a species in the databases, high intraspecific variability for COI between specimens was observed for some of them. For instance, the species *E. macrospilos* and *E. tauvina* have a COI genetic divergence within species of 7.4 and 6.2% respectively, values very largely above what is known of fish intraspecific diversity for this marker [40]. These species, and some others, were represented by multiple, disjoint clusters. For example, for *Epinephelus tauvina* one cluster of database sequences included our seven *E. tauvina* sequences, while another cluster matched with our eleven sequences of *E. coioides*. The species *E. akaara*, *E. amblycephalus*, *E. diacanthus*, *E. longispinis*, *E. macrospilos*, *E. sexfasciatus* and *Variola louti* were especially subject to such high intraspecific divergence, a well-known indication of misidentification or unresolved taxonomic issues. These reliability problems are all the more important when the molecular identification has medical applications, such as determining whether a sample can come from a ciguateric species. In order to select the more reliably identified specimens, we checked the supplementary information (voucher, geographic information) to corroborate the identification of the specimens. Moreover, in case of relatively high COI intraspecific variability within species (>1%), we selected the specimen collected closest to the type-locality to minimize taxonomic misidentifications and errors linked to possible cryptic species.

Second, the databases use some invalid names. While spelling errors are relatively straightforward to identify, like the use of *Cephalopholis miniatus* instead of *C. miniata* or *Pseudogramma polyacanthum* instead of *P. polyacantha*, nomenclaturally invalid names were also present in the databases. For instance, *E. 'fario'* (Thunberg, 1793), represented by eight specimens in GenBank including two COI sequences, has been pointed out by Heemstra and Randall [18] as a synonym of *E. longispinis* (Kner, 1864) and as a *nomen dubium*. Randall and Heemstra [41] also regarded this species as unidentifiable. We did not include these sequences in our dataset.

### Ethics statement

Fish were dead at the time we acquired them for study, having been commercially caught, and available for purchase at the Nouméa fish market; no permits were required for the described study, which complied with all relevant regulations.

<table><caption>**Table 2.** List of the primers used in this study; T° of hyb: temperature of hybridisation used to amplify the marker.</caption><thead><tr><th>Gene</th><th>Fragment size</th><th>Name</th><th>Primers</th><th>T° of hyb</th><th>Sources</th></tr></thead><tbody><tr><td rowspan="2">COI</td><td rowspan="2">≈650 bp</td><td>FishF1</td><td>5'-TCAACCAACCACAAAGACATTGGCAC-3'</td><td rowspan="2">48°C</td><td rowspan="2">[30]</td></tr><tr><td>FishR1</td><td>5'-TAGACTTCTGGGTGGCCAAAGAATCA3'</td></tr><tr><td rowspan="2">16S</td><td rowspan="2">≈410 bp</td><td>16SarL</td><td>5'-CGCCTGTTTATCAAAAACAT-3'</td><td rowspan="2">54°C</td><td rowspan="2">[75]</td></tr><tr><td>16SbrH</td><td>5'-CCGGTCTGAACTCAGATCACGT-3'</td></tr><tr><td rowspan="2">Rhodo</td><td rowspan="2">≈720 bp</td><td>Rh193</td><td>5'-CNTATGAATAYCCTCAGTACTACC-3'</td><td rowspan="2">52°C</td><td rowspan="2">[76]</td></tr><tr><td>Rh1039r</td><td>5'-TGCTTGTTCATGCAGATGTAGA-3'</td></tr><tr><td rowspan="2">Pkd1</td><td rowspan="2">≈850 bp</td><td>Pkd1F62</td><td>5'-CATGAGYGTCTACAGCATCCT-3'</td><td rowspan="2">50°C</td><td rowspan="2">[77]</td></tr><tr><td>Pkd1R952</td><td>5'-YCCTCTNCCAAAGTCCCACT-3'</td></tr><tr><td rowspan="2">TMO-4C4</td><td rowspan="2">≈540 bp</td><td>TMOF1</td><td>5'-CCTCCGGCCTTCCTAAAACCTCTC-3'</td><td rowspan="2">55°C</td><td rowspan="2">[78]</td></tr><tr><td>TMOR1</td><td>5'-CATCGTGCTCCTGGGTGACAAAGT-3'</td></tr></tbody></table>

doi:10.1371/journal.pone.0098198.t002
## DNA extraction, amplification, and sequencing

DNA was extracted from tissue samples using NucleoSpin 96 tissue kit (Macherey-Nagel, Düren, Germany) and five genes were amplified (Table 2).

Each PCR reaction was performed in a 20 µl final volume, containing 2 ng of DNA, 1X reaction buffer, 0.26 mM dNTP, 0.8 µM of each primer, 5% DMSO and 1.5 units of Taq polymerase (Qiagen). Thermocycles consisted of an initial denaturation step at 94°C for 2', followed by 37-55 cycles of denaturation at 94°C for 30", annealing at 48-55°C for 40" (Table 2) and extension at 72°C for 1'. The final extension was conducted at 72°C for 10'. Purification and cycle-sequencing reactions were performed at the Génoscope (Évry, France), using the BigDye Terminator version 3 sequencing kit, the GeneAmp PCR System 9700 and a capillary ABI3730 DNA Analyser, all from Applied Biosystems. Sequences were edited and assembled using Sequencher 4.9 (Gene Codes Corporation, Ann Arbor, MI, USA). Sequences for this study were deposited in GenBank (Table 1).

## Phylogenetic analyses and mapping

All markers were sequenced for both directions to confirm accuracy for each individual specimen. However, a few sequences could not be obtained (Table 1). Sequences were aligned using ClustalW as implemented in BioEdit version 7.0.5.3 [42] or by eye. The accuracy of automatic alignments was assessed by eye. 16S sequences were aligned manually and two portions (between positions 240-300 and 347-450) were removed from the alignment due to hypervariable regions that could not be aligned reliably.

94 species of Epinephelidae were included in our analyses, 2 species for Diplopriinae, 81 species for Epinephelinae, 8 species for Grammistinae and 3 species for Liopropominae.

To study (i) the monophyly of the Epinephelidae (*i.e* according to traditional classification Epinephelinae minus *Niphon*) and (ii) the relationships within this family, outgroups were chosen in the sub-families of Serranidae, the Serraninae and the Anthiinae [34]. We also included multiple, non-monophyletic outgroups: specimens from Perciformes (Cirrhitidae, Harpagiferidae, Niphonidae, Percidae, Trachinidae) and from Scorpaeniformes (Congiopodidae, Cyclopteridae, Scorpaeninae, Sebastinae, Synanceinae) placed as close relatives of Epinephelidae in the study of Lautrédou et al. [36].

## Phylogenetic analyses of the Epinephelidae

All markers were first analysed separately, tested for incongruence and, since none was detected, concatenated in two different datasets. The combined and separate analyses have three different aims (i) obtaining the most robust phylogenetic reconstruction while including a maximal species representation, *i.e.* including downloaded sequences (dataset 1), (ii) maximising the number of markers, even if it includes less species (dataset 2), and (iii) best representing the diversity of the Epinephelidae with the COI gene alone (dataset 3). Dataset 1 included the concatenated COI, 16S and TMO-4C4 sequences. Dataset 2 included all five concatenated markers, *i.e.* COI, 16S, TMO-4C4, Rhodopsin and Pkd1 sequences.

The best-fitting models of nucleotide evolution for each gene and for the concatenation were determined based on the Akaike Information Criterion (AIC) implemented in ModelTest 2.3 [43] in conjunction with PAUP 4.0b10 [44]. The GTR + I + Γ model was selected for each marker and for the two datasets.

Trees were inferred using two probabilistic approaches: maximum likelihood with a non-parametric bootstrap (BP) using

RAxML 7.2.8 [45,46] and Bayesian Inference [47]. Maximum likelihood (ML) analyses were carried out online on the CIPRES Science Gateway (The CIPRES Portals. URL: http://www.phylo.org/sub sections/portal) with RAxML-HPC BlackBox (7.2.7) [45]. Datasets were partitioned by codon position for each marker (except 16S) and by marker and by codon position for datasets 2 and 3. BI analyses were performed with MRBAYES version 3.2.1 [47] using 75,000,000, 75,000,000, and 30,000,000 generations for datasets 1, 2 and 3 respectively, with sampling every 1,000 generations and four Metropolis-coupled Markov chains Monte Carlo (MCMCMC). The parameter estimates and convergence were checked using Tracer version 1.4 [48]. The first 25% of sampled trees were considered burn-in trees, and discarded prior to constructing a 50% majority rule consensus trees. Posterior probabilities (PP) (Bayesian analysis) and Bootstrap probabilities (BP) were used as indicators of node credibility; P95% was considered significant [49]. Two independent analyses were conducted to check for convergence of the results.

## Mapping of potential ciguatera poisonous species

Potentially ciguatera affected species were mapped onto the phylogenetic tree inferred from all the available reliable sequences for the COI fragment. Although several relationships were not resolved compared to the combined analyses, the COI topology is congruent with the combined analyses topologies, and moreover best represents the diversity of the Epinephelidae.

A large number of studies on ciguatera were reviewed to establish a list of ciguatera affected species [7, 11, 12, 19, 20, 50-59]. However, there were very few publications where species names were precisely indicated. We ended up by using five [11, 20, 53, 54, 58], as well as the available ciguatera data in Fishbase database [17]. To detect a potential pattern in the evolution of the occurrence of ciguatera in Epinephelinae, we quantified the strength of phylogeny-trait association. The MrBayes phylogeny from the dataset 3 was used as an input in the BaTS software [60], with the occurrence of the ciguatera coded as absent/present. 10.000 trees from the posterior set of trees from MrBayes analysis were randomly selected after removing the first 7.5 million generations as a burnin according to Tracer [48], and re-rooted using the outgroups and a custom-made R script [61] (script available on request). BaTS estimates several statistics: the parsimony score (PS), the association index (AI) and the monophyletic clade (MC) and tests their significance against a null distribution (obtained by reshuffling 100 times the ciguatera states on the tips). A strong phylogeny-trait association is identified by low PS and AI scores and a high MC score.

## Evolution of the ciguatera fish poisoning (CFP)

Ancestral character state reconstructions for ciguatera fish poisoning were conducted using the maximum-likelihood method implemented in Mesquite 2.75 [62]. Recognized species were assigned different states (following a bibliographic survey, see above): absence (0) and presence (1) of CFP. Ancestral states were reconstructed for all Bayesian trees retained from the analysis of the combined data set and their mean likelihood was then plotted on the maximum clade credibility tree.

## Results

### Phylogenetic relationships within Epinephelidae

Dataset length and number of variable sites are reported in Table 3. 47 and 58 sequences were obtained for the mitochondrial markers COI and 16S, respectively. 54 sequences were obtained for TMO-4C4, 50 for Rhodopsin and 50 for Pkd1.
<table><caption>**Table 3.** Dataset composition, marker size and information.</caption><thead><tr><th>Dataset</th><th>Number of species</th><th>Markers</th><th>Length (bp)</th><th>Conserved sites</th><th>Variable sites</th></tr></thead><tbody><tr><td rowspan="3">1</td><td rowspan="3">102</td><td>COI</td><td>651</td><td>372</td><td>279</td></tr><tr><td>16S</td><td>410</td><td>254</td><td>156</td></tr><tr><td>TMO-4C4</td><td>511</td><td>271</td><td>240</td></tr><tr><td rowspan="5">2</td><td rowspan="5">60</td><td>COI</td><td>651</td><td>382</td><td>269</td></tr><tr><td>16S</td><td>410</td><td>278</td><td>132</td></tr><tr><td>TMO-4C4</td><td>536</td><td>330</td><td>206</td></tr><tr><td>Pkd1</td><td>861</td><td>422</td><td>439</td></tr><tr><td>Rhodopsin</td><td>723</td><td>470</td><td>253</td></tr><tr><td>3</td><td>112</td><td>COI</td><td>651</td><td>371</td><td>280</td></tr></tbody></table>

doi:10.1371/journal.pone.0098198.t003

The concatenated phylogeny (Fig. 1A) based on the three markers with the largest sampling (COI, 16S and TMO-4C4) recovered the monophyly of Epinephelidae comprising the four subfamilies, Diploprioninae, Epinephelinae, Grammistinae, and Liopropominae with a weak support. The Grammistinae was monophyletic. The three species of the genus *Liopropoma* of the Lioprominae were grouped. The Diploprioninae and the Epinephelinae were together not monophyletic because of the inclusion of *Belonoperca chabanaudi* within Epinephelinae. However, to the exception of the insertion of *B. chabanaudi*, all Epinephelinae specimens were grouped. In the Epinephelinae, the genera *Variola* and *Plectropomus* were together monophyletic with a strong support. *Saloptia* was sister-group to *Plectropomus. Cephalopholis* was monophyletic but divided in two robust groups with the inclusion of *C. colonus* and *C. rogaa* (earlier combinations: *Paranthias colonus* and *Aethaloperca rogaa*). *Epinephelus leuco-grammicus* and *E. altivelis* (earlier combinations: *Anyperodon leucogrammicus* and *Cromileptes altivelis*, both monotypic genera) were included within the clade *Epinephelus. Mycteroperca* was sister-group to *Epinephelus* and monophyletic with a strong support. *D. dermatolepis* was sister-group to *Hyporthodus, Mycteroperca* and *Epinephelus*. The subfamily Grammistinae was sister-group to the Liopropominae and *Diploprion bifasciatum* from the Diploprioninae.

The concatenated phylogeny (Fig. 1B) based on all five markers (COI, 16S, TMO-4C4, Rhodopsin and Pkd1) confirmed the monophyly of Epinephelinae as well as the three subfamilies Epinephelinae, Grammistinae and Liopropominae. No specimen from the subfamily Diploprioninae was available for the Rhodopsin and Pkd1 markers; the monophyly of this subfamily cannot be evaluated with this second dataset.

Within Epinephelinae and as suggested by dataset 1, the monophyly of *Mycteroperca, Plectropomus* and *Variola* were well supported. *Cephalopholis* also constituted a robust monophylum with two clades, one containing *C. argus* and *C. boenak* and the other the rest of the species (as observed with the previous tree). *Epinephelus* was monophyletic with a weak support. *D. derma-tolepis* was sister-group to a clade including *Epinephelus* and *Mycteroperca*.

Although several relationships were not resolved compared to the combined analyses (Fig. 1), the COI topology (Figure S1) is congruent with them, and has a much larger species sampling. The subfamilies Grammistinae and Liopropominae are monophyletic.

With the exception of the insertion of the species *Epinephelus poecilonotus* and *E. haifensis* included in the genera *Mycteroperca* and *Hyporthodus* respectively, all *Mycteroperca* and *Hyporthodus* species were grouped.

Like in the concatenated analyses (Figs. 1A-B), the monotypic genera *Anyperodon* and *Cromileptes* were together included in *Epinephelus*. The subfamily Diploprioninae including the three genera *Beloperca, Diploprion and Aulacocephalus* was not monophyletic, with the exclusion of *Beloperca chabanaudi*. *Diploprion bifasciatum* and *Aulacocephalus temminckii* constituted a robust clade.

## Mapping of potentially ciguatera affected species on the COI tree

Twenty nine ciguateric species were found in a review of the literature. Most of the species affected by ciguatera belong to the subfamily Epinephelinae (Fig. 2). Only two species, outside the Epinephelinae, *Grammistes sexlineatus* and *Rypticus saponaceus* are ciguateric. Within Epinephelinae, three genera contain multiple ciguateric species. All species included in the genera *Plectropomus* (4) and *Variola* (2) are ciguateric. In the genus *Epinephelus*, four species out of seven are ciguateric in one clade (*E. fasciatus*, *E. hexagonatus*, *E. melanostigma*, *E. merra*, *E. retouti*, *E. spilotoceps* and *E. tauvina*) (Fig. 2). In another, three species out of four are ciguateric (*E. quoyanus*, *E. macrospilos*, *E. howlandi* and *E. rivulatus*). The other ciguateric species are dispersed in clades where most species are not known to be affected. Within the other genera of Epinephelinae, *Mycteroperca* and *Hyporthodus*, only one species is ciguateric, *M. bonaci*.

The analysis of the presently available data on presence or absence of ciguatera fish poisoning (CFP) for each species mapped the absence of CFP as the ancestral state of Epinephelidae. The ancestral state analysis underlines the multiple appearances of CFP in this group with at least 10 events in the subfamily Epinephelinae (Fig. 2).

The Bayesian analysis of the phylogenetic signal (Table 4) revealed the occurrences of ciguatera in groupers species have a moderate non-random association with phylogeny; only the MC statistic for the occurrence of ciguatera shown significant P-value (p = 0.04), whereas the PS statistic shown a marginally significant P-value (p = 0.080). Both the AI statistic and the MC value for the ancestral state were not significant (p = 0.24 and p = 0.55, respectively).
<figure><img src="image_2.png" /><figcaption>**Figure 1. Phylogenetic relationships within the Epinephelidae.** Bayesian inference phylogram obtained from phylogenetic analyses of the dataset 1 (A) based on the concatenation of three genes, COI, 16S and TMO-4C4, and the dataset 2 (B) obtained with the concatenation of five genes, COI, 16S, TMO-4C4, Rhodopsin and Pkd1, both under the GTR + I + Γ model. Epinephelidae are highlighted in colour (pink and blue). Each sub families are shown in alternate blue and pink colours. Values at nodes indicate Bayesian posterior probabilities (PP) and maximum likelihood bootstrap percentages (BP). Black circles indicate nodes supported by posterior probability ≥95% and ML bootstrap probability ≥75%. doi:10.1371/journal.pone.0098198.g001</figcaption></figure>

Discussion

Systematics and taxonomy within Epinephelidae

As in Craig and Hastings [37] (but contradicted by Craig et al.
[32], because of the absence of morphological synapomorphy), our

result show the inclusion of *Epinephelus leucogrammicus* and *E. altivelis* (previous combinations *Anyperodon leucogrammicus* and *Cromileptes altivelis*) within *Epinephelus* but also the inclusion of *Cephalopholis rogaa* and *C. colonus* within *Cephalopholis* (previous combinations *Aethaloperca rogaa* and *Paranthias colonus*).
<figure><img src="image_3.png" /><figcaption>**Figure 2. Ancestral ciguatera fish poisoning (CFP) reconstruction of Epinephelidae.** Bayesian cladogram of the COI dataset with maximum likelihood estimates of ancestral CFP states. Pie charts correspond to average likelihoods for each state. Percentage values are given for nodes of interest.</figcaption></figure>

doi:10.1371/journal.pone.0098198.g002

Craig and Hastings [37] included in the genus *Mycteroperca* several species usually considered as members of *Epinephelus*. Based on our results, we confirm these new combinations. In addition, we proposed to transfer *E. poecilonotus* (Temminck and Schlegel, 1842) in *Mycteroperca* as *M. poecilonota*.

They also observed a monophyletic lineage distinct from the remaining species of *Epinephelus* and from *Mycteroperca*, for which they resurrected the oldest available generic name *Hyporthodus*. Our topologies corroborate this lineage, which includes at least two species, *H. septemfasciatus* and *H. haifensis* (previously *Epinephelus haifensis*).
<table><caption>**Table 4.** Results of the Bayesian phylogeny-trait association.</caption><thead><tr><th rowspan="2">Statistic</th><th colspan="2">Observed Distribution</th><th colspan="2">Null Distribution</th><th rowspan="2">P-value<br>significance</th></tr><tr><th>observed mean</th><th>95% CI</th><th>null mean</th><th>95% CI</th></tr></thead><tbody><tr><td>AI</td><td>4.2</td><td>3.7-4.6</td><td>4.6</td><td>3.6-5.5</td><td>0.240</td></tr><tr><td>PS</td><td>23.7</td><td>23.0-24.0</td><td>26.3</td><td>23.3-28.7</td><td>0.080*</td></tr><tr><td>MC (absence of ciguatera)</td><td>6.6</td><td>6.0-9.0</td><td>6.1</td><td>4.1-7.9</td><td>0.550</td></tr><tr><td>MC (occurrence of ciguatera)</td><td>4.0</td><td>4.0-4.0</td><td>2.3</td><td>1.6-3.1</td><td>0.040**</td></tr></tbody></table>

Association index (AI), parsimony score (PS), and monophyletic clade (MC) and their significance. While AI and PS indices test for the overall phylogeny and all the characters at once, MC is drawn to specifically quantify the phylogenetic signal for each specific character (occurrence of ciguatera contamination). Asterisk indicates significant values (\*: p≤0.1; \*: p≤0.05).

doi:10.1371/journal.pone.0098198.t004

The relationships between *Belonoperca chabanaudi*, *Diploprion bifasciatum* and *Aulacocephalus temminckii* in Figure S1 question the monophyly of the Diploprionini. In Figure 1A where *A. temminckii* was not represented, *B. chabanaudi* is included in the subfamily Epinephelinae (PP = 0.96, BP<sub>ML</sub> = 54). *D. bifasciatum* constituted, with a weak support, the sister-group of the Grammistinae with the Liopropominae. Since the support for several deeper clades is low at least for maximum likelihood analyses, we prefer not to discuss further the phylogeny of these groups.

Within Epinephelidae, there is no modification of the subfam-
ilies Grammistinae and Liopropominae. The Epinephelinae now
include only twelve genera (minus *Anyperodon*, *Cromileptes* and
*Paranthias* vs. 15 previously). In our results (Fig. 1), *Variola* and
*Plectropomus-Saloptia* are in basal position relative to the other
Epinephelinae genera.

Ciguatera fish poisoning

In the studies on ciguatera reviewed, most were fairly unhelpful
to establish a list of ciguatera-affected species with scientific names.
Most publications use vernacular names only, and many of these
designate several fish species (see [2,6]). Vernacular names might
be important for local communication and consumer warnings.
However, scientific names should be systematically associated to
them to enhance precision and communication between localities,
as vernacular name use varies with geography. Feedback
questionnaires [63,64] were provided to public health and fisheries
department staff to collect ciguatera poisoning data and then to
put more monitoring and research into place. Traceback
investigations of fish associated with outbreaks provide valuable
information regarding fishing areas associated with CFP. Howev-
er, in a weekly report, the Centers for Disease Control and
Prevention pointed to limitations in the personal feedbacks [65];
where physician reports were unavailable, the symptoms were
based entirely on self-report or second hand reports from family
members and may be wrong. Moreover, additional cases might
have occurred but were unrecognized because involved physicians
were not aware of the need to make an appropriate diagnosis and
to report, especially in countries where the public health network is
weakly organized.

With the limitations of our current knowledge in mind, 29
among 163 grouper species are considered ciguateric in the world.
While there is no correlation between the phylogeny and the
currently known ciguatera status of the species, some clades
include a higher number of possibly ciguateric species, and might
be interesting to investigate in order to determine the ciguateric
status of species where it has not been described yet. The recent

changes in classification based on the reanalyses of the relation-
ships within Epinephelidae ([37], this study) have an impact on the
interpretation of the ciguatera distribution in the genera. Multiple
species from the genus *Epinephelus* that were not reported to be
ciguateric belong in fact to other genera that have very few
ciguateric species (*Hyporthodus* and *Mycteroperca*). As for the
genera *Plectropomus* and *Variola*, the species included in our
datasets are all considered potentially ciguateric (4 and 2 species
respectively). In the three other species contained in the genus
*Plectropomus*, only one, *P. oligacanthus*, is ciguateric according to
FishBase.

In our study, even if the COI gene gives little support for deeper
nodes in phylogenetic analyses, it performs well in Epinephelidae
for interspecific relationships, as already observed for other
Teleosts and taxonomic groups [66,67]. Better phylogenetic
performance means that even if an identical sequence is not
available for identification in the database, the position of the
unknown sequence in a phylogenetic analysis can provide
information about the genus it belongs to, and which species it
is most closely related to. Consequently, samples from species that
cannot be readily identified, but fall within these groups in an
analysis might best be considered potentially ciguateric when
tested for ciguatera case suspicions. This, added to the large
dataset already available, and the more stringent guidelines of the
Barcode of Life project compared to other sequence databases,
makes it a very good choice for identification in the group.

Many questions remain to be answered pertaining to the
production and accumulation of ciguatoxins and the subsequent
occurrence of ciguatera fish poisoning. For example, why does
Ciguatera affect only some species in a given locality, why does
Ciguatera not affect species over their whole range? Explanations
might be found on the life history traits of groupers. The diet, the
location on the reef and the size of grouper species could explain
some of these ciguateric patterns. Yet no clear tendency based on
these data appears. On the other hand, some species like
*Epinephelus macrospilos* or *Cephalopholis miniata* are ciguateric,
for example, in New Caledonia (Fig. 2) but not in their whole
range. Richlen et al. [68] were some of the first authors to
correlate the geographic patterns of fish toxicity with the
preponderance of highly toxic strains of *Gambierdiscus* spp. They
showed that *G. toxicus* was not a single cosmopolitan species, but
instead was a species complex comprised of several distantly
related groups co-occurring across geography. Thus the range of
*G. toxicus* appears to be much smaller than the range of the fish
species, but this also depends on the correctness of our knowledge
about the systematics and distribution of the fish species. The
presence and relative abundance of the members of this species
complex among geographic regions may help explaining patterns of ciguatera toxicity, particularly if differences in physiology and/or toxin-producing capabilities also exist among these groups. In addition, a better understanding of the three groups of ciguatoxins produced by *Gambierdiscus* species will be very useful.

As suggested by several authors, due to recent climate change or human activities, harmful algal blooms (HAB), including *Gambierdiscus* spp. [1,8,9] might become more frequent and more geographically widespread. A consequence would be that ciguatera starts affecting fish species in other localities in their range where there was no previously recorded problem, affecting populations that will be not aware of the ciguateric risk of fish consumption. The first reports of consumer illness and detection of ciguatoxic fish from the Canary Islands [55] seems to be consistent with such a geographic expansion of ciguatera. Another recent expansion of dinoflagellates, and subsequent incidences of CFP illness, was recorded in the northern Gulf of Mexico, USA [10]. Of the 29 ciguateric grouper species, 16 have a very large distribution in the Indo-Pacific Ocean. The ranges extend from the East of Africa to the West of America, including the Red Sea. Current distribution areas might help us determine future changes in the prevalence area of ciguatera.

## Barcoding, a useful tool for various applications

As all our knowledge on ciguateric species hinges on correct identification of the species involved in CFP, DNA barcoding is an essential tool in controlling CFP, but also investigating mislabelling of seafood or endangered species monitoring. Fishermen and restaurant owners do not hesitate to sell fish of lower quality under erroneous labelling [16], a problem compounded in groupers by overfishing. Because of the very high price of live groupers, the mislabelling of fishes and the dangers of ciguatera, accurate species identification is part of what makes the quality of the fish meat [16]. The misidentification of several species was largely discussed by Stewart et al. [12]. They encountered nine misidentifications in their molecular analyses, and remarked that eight of the nine misnamed fish contain ciguatera toxins. While the authors did not suggest deliberate substitutions, they wondered about this high rate. Some confusions were not surprising in regard of the similar morphology of the species (Spanish mackerel - *Scomberomorus commerson* and *Scomberomorus queenslandicus* for instance). However, they pointed out that selling the Spanish mackerel – a ciguatera-prone species banned in Platypus Bay, Queensland – as red snapper and swordfish, neither known to be problematic for ciguatera, represents a potential risk of bypassing ciguatera prevention strategies [12].

In such a case, DNA barcoding has proved to be a very useful tool to quickly and easily identify seafood species sold on fish markets or on restaurants, thus helping to avoid at least some of the risky species. The mislabelling studies have also had a positive effect on at least some parts of the retail sector [69]. Moreover, meal remnants in CFP cases often cannot be identified without molecular tools, where proper identification is critical for improving the ciguatera affected species list and warnings. DNA barcoding-based traceability procedures were implemented in several U.S. state and federal laboratories [65]. As the symptoms appear in minutes to hours after the ingestion of contaminated seafood [13], recovery of samples for sequencing is possible, at least in some cases. For most types of cooked seafood, the full barcode can be obtained. For severely degraded or heavily processed products (e.g. canned), the sequencing of shorter sequences (i.e. Mini barcodes, see [70]), also works for identification.

## The key role of the public databases

For DNA barcoding to be a useful tool for identifying species, it needs complete (or at least large) and reliable reference sequences available in public databases. For an example of both the usefulness of barcoding methods and the limits of the available databases, a sample of fish tissue collected at the fish market in Nouméa (New Caledonia) showed 100% identity in the databases with those of two distinct species. Both are deep-sea species, as was the grouper from the market according to the fishmongers. *H. ergastularius* is not recorded from New Caledonia, but *H. octofasciatus* was recorded, sometimes under other names [71]. Unfortunately, identical sequences for *H. ergastularius* and *H. octofasciatus* in databases preclude an identification of this sample. Currently, available grouper sequences contain several cases of misidentification both in GenBank and in BOLD. The numerous grey-brown species with dots (*E. tauvina*, *E. akaara*, *E. diacanthus*, *E. amblycephalus*, *E. longispinis*, *E. sexfasciatus* and *E. macrospilos*) are particularly affected, with the additional problem that young specimens of some species resemble adults of smaller species. This is not an isolated case. According to Vilgalys [72] up to 20% of the named sequences in public databases may be misidentified. Hassanin et al. [73] suggested to annotate database sequences through an additional "external expertise" field, and there is indeed a possibility to add comments to data in BOLD. However, a lot remains to be done before all available identifications can be trusted, even with the additional geographical and voucher information required by BOLD. There is an ongoing effort by the BOLD crew to flag dubious sequences. They developed a system to grade the level of reliability of the identification in BOLD [74], but it has not been applied to all sequences yet. In the end, a morphological study of the voucher specimen remains necessary, but this is not practical when fast identification is needed such as for ciguateric sample identifications. While currently neither reliability nor comprehensivity are at hand for Epinephelinae sequences, our study has added 47 COI sequences for carefully identified and vouchered specimens.

## Supporting Information

**Figure S1** Phylogenetic relationships within Epinephelidae. Bayesian inference phylogram obtained from phylogenetic analyses of the COI under the GTR + I + Γ model. Values at nodes indicate Bayesian posterior probabilities (PP) and maximum likelihood bootstrap percentages (BP). Black circles indicate nodes supported by posterior probability ≥95% and ML bootstrap probability ≥75%.

(TIF)

## Acknowledgments

Anne-Claire Lautrédou (MNHN, Paris), Paolo Galli (University of Milano, Italy: material from Red Sea), Philippe Béarez (MNHN, Paris: Ecuador, Philippines), Serges Planes (Moorea Biocode Project, French Polynesia), Monica Mwale (SAIAB, South Africa), Dominique Ponton (IRD, Nouméa, New Caledonia), Jean-Dominique Durand (IRD, Dakar, Senegal), Rishen Liang (Guangzhou University, China), Lassâd Neifar (Sfax University, Tunisia) provided tissues samples and/or sequences. The Australian Node of the CReefs Global Research Initiative, Thomas Cribb, Terrence Miller, and Holly Heiniger helped in obtaining tissue samples during a mission in Heron Island. Martine Desoutter (MNHN, Paris) and Philippe Bouchet (MNHN, Paris) discussed nomenclatural questions and Gaël Denys, Frédéric Busson and Gaël Lancelot contributed to the analyses, the data collection and text corrections respectively. Sincere thanks to the staff of the "Service de Systématique Moleéculaire" (UMS2700 CNRS-MNHN). This work is part of the project @ SPEED-ID "Accurate SPEciEs Delimitation and IDentification of eukaryotic biodiversity using DNA
markers" proposed by F-BoL, the French Barcode of life initiative, and the network "Bibliothèque du Vivant" funded by the CNRS, the Muséum National d'Histoire Naturelle, the INRA and the CEA (Genoscope). We thank the anonymous reviewers for their very helpful comments on the manuscript.

## References

1. Dickey RW, Plakas SM (2010) Ciguatera: A public health perspective. Toxicon 56: 123-136.

2. Laurent D, Yeeting B, Labrosse P, Gaudechoux JP (2005) Ciguatera: a Field Reference Guide.(Ciguatera: un guide pratique). Secretariat of the Pacific Community. 91 p.

3. De Fouw JC, Van Egmond HP, Speijers GJA (2001) Ciguatera fish poisoning: a review. Available: http://wwwrivmnl.bibliothek/rapporten/388802021pdf. Accessed 25 June 2014.

4. Lewis RJ, Jones A (1997) Characterization of ciguatoxins and ciguatoxin congeners present in ciguateric fish by gradient reverse-phase high-performance liquid chromatography mass spectrometry. Toxicon 35: 159-168.

5. Lewis RJ (2001) The changing face of ciguatera. Toxicon 39: 97-106.

6. Friedman MA, Fleming LE, Fernandez M, Bienfang P, Schrank K, et al. (2008) Ciguatera fish poisoning: Treatment, prevention and management. Marine Drugs 6: 456-479.

7. Kipping R, Eastcott H, Sarangi J (2006) Tropical fish poisoning in temperate climates: food poisoning from ciguatera toxin presenting in Avonmouth. Journal of Public Health 28: 343-346.

8. Chateau-Degat ML, Chinain M, Cerf N, Gingras S, Hubert B, et al. (2005) Seawater temperature, *Gambierdiscus* spp. variability and incidence of ciguatera poisoning in French Polynesia. Harmful Algae 4: 1053-1062.

9. Moore SK, Trainer VL, Mantua NJ, Parker MS, Laws EA, et al. (2008) Impacts of climate variability and future climate change on harmful algal blooms and human health. Environmental Health 7: (Suppl 2):S4.

10. Villarreal TA, Hanson S, Qualia S, Jester ELE, Granade HR, et al. (2007) Petroleum production platforms as sites for the expansion of ciguatera in the northwestern Gulf of Mexico. Harmful Algae 6: 253-259.

11. Clua E, Brena PF, Lecasble C, Ghnassia R, Chauvet C (2011) Prevalence and proposal for cost-effective management of the ciguatera risk in the Noumea fish market, New Caledonia (South Pacific). Toxicon 58: 591-601.

12. Stewart I, Eaglesham GK, Poole S, Graham G, Paulo C, et al. (2010) Establishing a public health analytical service based on chemical methods for detecting and quantifying Pacific ciguatoxin in fish samples. Toxicon 56: 804-812.

13. Lewis RJ (2006) Ciguatera: Australian perspectives on a global problem. Toxicon 48: 799-809.

14. Carvalho DC, Neto DAP, Brasil BSAF, Oliveira DAA (2011) DNA barcoding unveils a high rate of mislabeling in a commercial freshwater catfish from Brazil. Mitochondrial DNA 22: 97-105.

15. Cline E (2012) Marketplace substitution of Atlantic salmon for Pacific salmon in Washington State detected by DNA barcoding. Food Research International 45: 388-393.

16. Wong EHK, Hanner RH (2008) DNA barcoding detects market substitution in North American seafood. Food Research International 41: 828-837.

17. Froese R, Pauly D (2012) FishBase. World Wide Web electronic publication. Available: www.fishbase.org.

18. Heemstra PC, Randall JE (1993) FAO Species Catalogue. Vol. 16. Groupers of the world (Family Serranidae, Subfamily Epinephelinae). An annotated and illustrated catalogue of the grouper, rockcod, hind, coral grouper and lyretail species known to date. Rome: FAO. 382 p.

19. Laboute P, Grandperrin R (2000) Poissons de Nouvelle-Calédonie. Nouméa, New Caledonia: Éditions Catherine Ledru. 520 p.

20. Oshiro N, Yogi K, Asato S, Sasaki T, Tamanaha K, et al. (2010) Ciguatera incidence and fish toxicity in Okinawa, Japan. Toxicon 56: 656-661.

21. Harikrishnan R, Balasundaram C, Heo M-S (2010) Molecular studies, disease status and prophylactic measures in grouper aquaculture: Economic importance, diseases and immunology. Aquaculture 309: 1-14.

22. Ottolenghi F, Silvestri C, Giordano P, Lovatelli A, New MB (2004) Capture-based aquaculture. The fattening of eels, groupers, tunas and yellowtails. Rome: FAO. 308 p.

23. Pierre S, Gaillard S, Prevot-D'Alvise N, Aubert J, Rostaing-Capailon O, et al. (2008) Grouper aquaculture: Asian success and Mediterranean trials. Aquatic Conservation-Marine and Freshwater Ecosystems 18: 297-308.

24. Sadovy de Mitcheson Y, Craig MT, Bertoncini AA, Carpenter KE, Cheung WWL, et al. (2013) Fishing groupers towards extinction: a global assessment of threats and extinction risks in a billion dollar fishery. Fish and Fisheries 14: 119-136.

25. FAO (2010) Capture Production 1950-2008. FAO Fisheries Department, Fishery Information, Data and Statistics Unit.

26. Ward RD, Hanner R, Hebert PDN (2009) The campaign to DNA barcode all fishes, FISH-BOL. Journal of Fish Biology 74: 329-356.

27. Collins RA, Armstrong KF, Meier R, Yi YG, Brown SDJ, et al. (2012) Barcoding and border biosecurity: Identifying cyprinid fishes in the aquarium trade. Plos One 7.

## Author Contributions

Conceived and designed the experiments: CS AD. Performed the experiments: CS DDH. Analyzed the data: CS DDH. Contributed reagents/materials/analysis tools: CC. Wrote the paper: CS DDH AD JIJ.

28. Hanner R, Becker S, Ivanova NV, Steinke D (2011) FISH-BOL and seafood identification: Geographically dispersed case studies reveal systemic market substitution across Canada. Mitochondrial DNA 22: 106-122.

29. Lakra WS, Verma MS, Goswami M, Lal KK, Mohindra V, et al. (2011) DNA barcoding Indian marine fishes. Molecular Ecology Resources 11: 60-71.

30. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PD (2005) DNA barcoding Australia's fish species. Philosophical transactions of the Royal Society of London Series B, Biological Sciences 360: 1847-1857.

31. Vernooy R, Haribabu E, Muller MR, Vogel JH, Hebert PDN, et al. (2010) Barcoding life to conserve biological diversity: beyond the taxonomic imperative. PLoS Biol 8.

32. Craig MT, Sadovy de Mitcheson YJ, Heemstra PC (2012) Groupers of the world: a field and market guide: NISC (Pty) Ltd, Grahamstown.

33. Ekrem T, Willassen E, Stur E (2007) A comprehensive DNA sequence library is essential for identification with DNA barcodes. Molecular Phylogenetics and Evolution 43: 530-542.

34. Smith WL, Craig MT (2007) Casting the percomorph net widely: The importance of broad taxonomic sampling in the search for the placement of serranid and percid fishes. Copeia: 35-55.

35. Nelson JS (2006) Fishes of the world: Fourth edition. 601 p.

36. Lautréod AC, Motomura H, Gallut C, Ozouf-Costaz C, Cruaud C, et al. (2013) New nuclear markers and exploration of the relationships among Serraniformes (Acanthomorpha, Teleostei): the importance of working at multiple scales. Molecular Phylogenetics and Evolution 67: 140-155.

37. Craig MT, Hastings PA (2007) A molecular phylogeny of the groupers of the subfamily Epinephelinae (Serranidae) with a revised classification of the Epinephelini. Ichthyological Research 54: 1-17.

38. Ratnasingham S, Hebert PDN (2007) The Barcode of Life Data System (www. barcodinglife.org). Molecular Ecology Notes 7: 355-364.

39. Buhay JE (2009) "Coi-like" sequences are becoming problematic in molecular systematic and DNA barcoding studies. Journal of Crustacean Biology 29: 96-110.

40. Ward RD (2009) DNA barcode divergence among species and genera of birds and fishes. Molecular Ecology Resources 9: 1077-1085.

41. Randall JE, Heemstra PC (1991) Revision of Indo-Pacific groupers (Perciformes: Serranidae: Epinephelinae), with descriptions of five new species. Indo-Pacific Fishes 20: 1-332.

42. Hall TA (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp 41: 95-98.

43. Posada D, Crandall KA (2001) Evaluation of methods for detecting recombination from DNA sequences: Computer simulations. Proceedings of the National Academy of Sciences of the United States of America 98: 13757-13762.

44. Sofford DL (2002) PAUP*. Phylogenetic Analysis Using Parsimony* and other methods. Sinauer Associates.

45. Stamatakis A (2006) RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688-2690.

46. Stamatakis A, Hoover P, Rougemont J (2008) A rapid bootstrap algorithm for the RAxML web servers. Systematic Biology 57: 758-771.

47. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574.

48. Rambaut A, Drummond A (2007) Tracer v1.4. Available: http://beast.bio.ed.ac.uk/Tracer.

49. Leache AD, Reeder TW (2002) Molecular systematics of the Eastern fence lizard (*Sceloporus undulatus*): A comparison of parsimony, likelihood, and Bayesian approaches. Systematic Biology 51: 44-68.

50. Baumann F, Bourrat M-B, Pauillac S (2010) Prevalence, symptoms and chronicity of ciguatera in New Caledonia: Results from an adult population survey conducted in Nouméa during 2005. Toxicon 56: 662-667.

51. Chateau-Degat M-L, Huin-Blondey M-O, Chinain M, Darius T, Legrand A-M, et al. (2007) Prevalence of chronic symptoms of Ciguatera disease in French Polynesian adults. American Journal of Tropical Medicine and Hygiene 77: 842-846.

52. Lehane L, Lewis RJ (2000) Ciguatera: recent advances but the risk remains. International Journal of Food Microbiology 61: 91-125.

53. Lewis ND (1986) Epidemiology and Impact of Ciguatera in the Pacific - A Review. Marine Fisheries Review 48: 6-13.

54. Olsen DA, Nellis DW, Wood RS (1984) Ciguatera in the eastern Caribbean. U S National Marine Fisheries Service Marine Fisheries Review 46: 13-18.

55. Perez-Arellano JL, Luzardo OP, Brito AP, Cabrera MH, Zumbado M, et al. (2005) Ciguatera fish poisoning, Canary Islands. Emerging Infectious Diseases 11: 1981-1982.

56. Quod JP, Turquet J (1996) Ciguatera in Réunion Island (SW Indian Ocean): Epidemiology and clinical patterns. Toxicon 34: 779-785.
57. Schlaich C, Hagelstein J-G, Burchard G-D, Schmiedel S (2012) Outbreak of ciguatera fish poisoning on a cargo ship in the port of Hamburg. Journal of Travel Medicine 19: 238-242.

58. Wong CK, Hung P, Lee KLH, Kam KM (2005) Study of an outbreak of ciguatera fish poisoning in Hong Kong. Toxicon 46: 563-571.

59. Wong C-K, Hung P, Lee KLH, Mok T, Chung T, et al. (2008) Features of ciguatera fish poisoning cases in Hong Kong 2004–2007. Biomedical and Environmental Sciences 21: 521–527.

60. Parker J, Rambaut A, Pybus OG (2008) Correlating viral phenotypes with phylogeny: Accounting for phylogenetic uncertainty. Infection Genetics and Evolution 8: 239-246.

61. R Core Team (2013) R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria. Available: http://www.R-project.org.

62. Maddison WP, Maddison DR (2011) Mesquite: a modular system for evolutionary analysis. Version 2.75. Available: http://mesquiteproject.org.

63. Skinner MP, Brewer TD, Johnstone R, Fleming LE, Lewis RJ (2011) Ciguatera fish poisoning in the Pacific islands (1998 to 2008). Plos Neglected Tropical Diseases 5.

64. Tester PA, Feldman RL, Nau AW, Faust MA, Litaker RW (2009) Ciguatera fish poisoning in the Caribbean. Smithsonian contributions to the Marine Sciences 38: 301-311.

65. Centers for disease control and prevention (2013) Ciguatera fish poisoning - New York city, 2010-2011. JAMA 309: 1102-1104.

66. Dettai A, Berkani M, Lautredou AC, Couloux A, Lecointre G, et al. (2012) Tracking the elusive monophyly of nototheniid fishes (Teleostei) with multiple mitochondrial and nuclear markers. Marine Genomics 8: 49-58.

67. Mueller RL (2006) Evolutionary rates, divergence dates, and the performance of mitochondrial genes in Bayesian phylogenetic analysis. Systematic Biology 55: 289-300.

68. Richlen ML, Morton SL, Barber PH, Lobel PS (2008) Phylogeography, morphological variation and taxonomy of the toxic dinoflagellate *Gambierdiscus toxicus* (Dinophyceae). Harmful Algae 7: 614-629.

69. Mariani S, Ellis J, O'Reilly A, Bréchon AL, Sacchi C, et al. (2014) Mass media influence and the regulation of illegal practices in the seafood market. Conservation Letters.

70. Meusnier I, Singer GAC, Landry JF, Hickey DA, Hebert PDN, et al. (2008) A universal DNA mini-barcode for biodiversity analysis. BMC Genomics 9.

71. Fricke R, Kulbicki M, Wantiez L (2011) Checklist of the fishes of New Caledonia, and their distribution in the Southwest Pacific Ocean (Pisces). Stuttgarter Beitrage zur Naturkunde A Neue Serie (Biol) 4: 341-463.

72. Vilgalys R (2003) Taxonomic misidentification in public DNA databases. New Phytologist 160: 4-5.

73. Hassanin A, Bonillo C, Bui XN, Cruaud C (2010) Comparisons between mitochondrial genomes of domestic goat (*Capra hircus*) reveal the presence of numts and multiple sequencing errors. Mitochondrial DNA 21: 68-76.

74. Steinke D, Hanner R (2011) The FISH-BOL Collaborators protocol. Mitochondrial DNA 22 (S1): 10-14.

75. Palumbi SR (1996) Nucleic acids 2: the polymerase chain reaction. Molecular systematics Second edition.pp. 205-247.

76. Chen W-J, Bonillo Cl, Lecointre G (2003) Repeatability of clades as a criterion of reliability: a case study for molecular phylogeny of Acanthomorpha (Teleostei) with larger number of taxa. Molecular Phylogenetics and Evolution 26: 262-288.

77. Lautrédou AC, Bonillo C, Denys G, Cruaud C, Ozouf-Costaz C, et al. (2010) Molecular taxonomy and identification within the Antarctic genus *Trematomus* (Notothenioidei, Teleostei): How valuable is barcoding with COI? Polar Science 4: 333-352.

78. Streelman JT, Karl SA (1997) Reconstructing labroid evolution with single-copy nuclear DNA. Proceedings of the Royal Society B-Biological Sciences 264: 1011-1020.