Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A glimpse into the genetic diversity of the Peruvian seafood sector: Unveiling species substitution, mislabeling and trade of threatened species

  • Alan Marín ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    marin@fish.hokudai.ac.jp

    Current address: Laboratorio de Genética, Fisiología y Reproducción, Facultad de Ciencias, Universidad Nacional del Santa, Chimbote, Perú

    Affiliation Biodes Laboratorios Soluciones Integrales S.C.R.L., Tumbes, Perú

  • José Serna,

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Biodes Laboratorios Soluciones Integrales S.C.R.L., Tumbes, Perú

  • Christian Robles,

    Roles Investigation, Methodology, Writing – review & editing

    Affiliation Biodes Laboratorios Soluciones Integrales S.C.R.L., Tumbes, Perú

  • Beder Ramírez,

    Roles Investigation, Methodology, Writing – review & editing

    Affiliation Laboratorio Costero de Tumbes, Instituto del Mar del Perú-IMARPE, Tumbes, Perú

  • Lorenzo E. Reyes-Flores,

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Laboratorio de Genética, Fisiología y Reproducción, Facultad de Ciencias, Universidad Nacional del Santa, Chimbote, Perú

  • Eliana Zelada-Mázmela,

    Roles Investigation, Methodology, Writing – review & editing

    Affiliation Laboratorio de Genética, Fisiología y Reproducción, Facultad de Ciencias, Universidad Nacional del Santa, Chimbote, Perú

  • Giovanna Sotil,

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Laboratorio de Genética Molecular, Instituto del Mar del Perú-IMARPE, Lima, Perú

  • Ruben Alfaro

    Roles Conceptualization, Data curation, Investigation, Methodology, Writing – review & editing

    Affiliation Laboratorio de Biología Molecular, Facultad de Ciencias de la Salud, Universidad Nacional de Tumbes, Tumbes, Perú

Abstract

Peru is one of the world’s leading fishing nations and its seafood industry relies on the trade of a vast variety of aquatic resources, playing a key role in the country’s socio-economic development. DNA barcoding has become of paramount importance for systematics, conservation, and seafood traceability, complementing or even surpassing conventional identification methods when target organisms show similar morphology during the early life stages, have recently diverged, or have undergone processing. Aiming to increase our knowledge of the species diversity available across the Peruvian supply chain (from fish landing sites to markets and restaurants), we applied full and mini-barcoding approaches targeting three mitochondrial genes (COI, 16S, and 12S) and the control region to identify samples purchased at retailers from six departments along the north-central Peruvian coast. DNA barcodes from 131 samples were assigned to 55 species (plus five genus-level taxa) comprising 47 families, 24 orders, and six classes including Actinopterygii (45.03%), Chondrichthyes (36.64%), Bivalvia (6.87%), Cephalopoda (6.11%), Malacostraca (3.82%), and Gastropoda (1.53%). The identified samples included commercially important pelagic (anchovy, bonito, dolphinfish) and demersal (hake, smooth-hound, Peruvian rock seabass, croaker) fish species. Our results unveiled the marketing of protected and threatened species such as whale shark, Atlantic white marlin, smooth hammerhead (some specimens collected during closed season), shortfin mako, and pelagic thresher sharks. A total of 35 samples (26.72%) were mislabeled, including tilapia labeled as wild marine fish, dolphinfish and hake labeled as grouper, and different shark species sold as “smooth-hounds”. The present study highlights the necessity of implementing traceability and monitoring programs along the entire seafood supply chain using molecular tools to enhance sustainability efforts and ensure consumer choice.

Introduction

Peru is a major fishing country with rich marine biodiversity. About 1070 fish [1], 1018 molluscan [2], and 320 crustacean species [3] have been described from the Peruvian marine ecosystem, which is dominated by the cold nutrient-rich waters of the Humboldt Current [4]. The highly productive Peruvian upwelling system supports not only the world’s largest fishery for Peruvian anchovy (Engraulis ringens) [5, 6], but also other important planktivorous fish species (e.g., jack mackerel Trachurus murphyi) and their predators (e.g., bonito Sarda chiliensis, dolphinfish Coryphaena hippurus), which are valuable artisanal fishery resources [4]. The fishery sector plays a key role in the nation’s socio-economic growth with most artisanal production consumed directly through local markets [4]. In 2016, Peru was the fifth biggest producer of marine capture fisheries in the world, with a total production of 3.7 million tonnes [7]. Most (75%) of that production was due to anchovy catches, but high catches of other species such as Pacific chub mackerel (Scomber japonicus), jumbo flying squid (Dosidicus gigas), and South Pacific hake (Merluccius gayi), were also reported [8].

Furthermore, Peru has favorable conditions for fishery and aquaculture activities considering its 3080 km coastline, and 12000 lakes and lagoons [6, 9]. Peruvian aquaculture production makes up only 2% of the total seafood sector [9], and the main species include whiteleg shrimp (Penaeus vannamei), Peruvian scallop (Argopecten purpuratus), rainbow trout (Oncorhynchus mykiss), tilapia (Oreochromis niloticus), and Amazon fish paiche (Arapaima gigas) [9]. Fish consumption in Peru has increased in recent years. According to FAO statistics [10], the average annual per capita consumption of fish during 2013–15 was 21.8 kg, which was the highest in Latin America and the Caribbean. This increase was due in part to fish consumption campaigns launched by the Peruvian government [11]. The significant increase of fish consumption added to the depletion of marine resources and a high demand for quality seafood products, partially due to the growing number of Peruvian seafood restaurants (locally known as “cebicherias”), make this country’s sector a fertile area for mislabeling of species that are expensive, scarce or out of season.

Seafood fraud and species substitutions occur regularly in this sector and represent a global issue. Proper identification of food species is now a main concern not only for governments and companies, but also for consumers due to economic, regulatory, health and religious reasons [12]. However, conventional fish identification methods, based on morphological traits using field guides and taxonomic keys, may lead to misidentification when applied to morphologically similar or recently diverged species and can be utilized only when the full body is available. Mislabeling can happen at any point in the supply chain, from fisher to retailer; thus, determining how substitutions occur is complicated [13]. In a comprehensive analysis by Pardo et al. [14], which included 51 peer-reviewed seafood articles published from 2010 to 2015 and comprising 4500 samples, the average rate of reported misdescriptions was 30%. This is where the discriminatory power of DNA based tools can be successfully applied for seafood authentication, even if all morphological characters are gone after processing and cooking.

The most frequently used genetic marker for metazoan identification through DNA barcoding is a partial fragment (∼650 bp) located at the 5ʹ end of the mitochondrial cytochrome c oxidase subunit I (COI) gene. It has been successfully used to correctly identify different fresh and processed seafood samples [15]. Furthermore, in October 2011, the United States Food and Drug Administration (FDA) formally adopted DNA barcoding as the primary method for seafood identification [16]. The universal primers designed by Folmer et al. [17] are one of the most commonly used for the amplification of the “barcode COI region”; however these have failed to amplify PCR products from different marine organisms such as fishes, crabs, echinoderms, decapods, and scallops [18] (and references therein). Therefore new barcoding primer sets targeting more specific groups [19], as well as alternative mitochondrial genes such as cytochrome b, 12S rRNA, and 16S rRNA have been tested in different marine organisms [20, 21]. Another problem that challenges the amplification of full-length DNA barcoding fragment size (∼650 bp) is when dealing with samples that have been through extreme conditions such as high temperature and pressure in canning or cooking processes. Heat exposure and high pressure degrade large molecular weight genomic DNA to shorter size fragments mainly through enzymatic degradation, depurination, and hydrolysis, and sometimes highly degraded samples might display breaks or artifactual mutations [22]. To overcome this issue, shorter (100–200 bp) PCR fragments within the full-length barcode region (known as mini-barcoding markers) have proved to be an effective species identification tool when using degraded target DNA [23].

In spite of the increasing need to enforce regulations aiming for sustainable seafood industry and effective control of the trade of endangered species, few initiatives have been undertaken to evaluate the utility of molecular markers for authenticating products available from seafood retailers in South America. Most of those studies have surveyed single groups of Amazonian and Atlantic fish species from Brazil, including catfish and sawfish from supermarkets and fish markets [24, 25], Amazonian fish from local harbors and markets [26, 27], croaker filets from supplier companies [28], characiforms from street markets [29], and sharks from supermarkets [30]. From the South Pacific coast, reports include studies on Chilean species of commercial mollusks [31], commercial crabs from local markets [32], and salmon from supermarkets [33], while sharks from Peruvian fish landing sites [34] have also been reported.

Considering the scarcity of data regarding the authentication of seafood products in the Peruvian market and aiming to increase our knowledge of the species diversity along the supply chain from fish landing sites to markets and restaurants, this study assessed the utility of full and mini-barcoding markers for identifying a variety of local and imported fish and shellfish products in the seafood sector. Additionally, we evaluated the accuracy of fish labels and described the conservation status and current regulatory framework related to the threatened species detected in this study.

Materials and methods

Sample collection

A total of 143 national and imported seafood samples were collected from July 2016 to March 2018, covering a wide range of presentations including fresh, refrigerated, frozen, canned, dried, cooked, packed, dehydrated, marinated, fish burger, and fish roe. They included 48 samples from restaurants (RT), 29 from supermarket chains (SMC), 24 from markets (MK), 22 from fish landing sites (FLS), 12 from multimarket (MM), seven from wholesale fish market (WFM), and one from a grocery store (GS). Samples were collected along the north-central Peruvian coast in different cities from six departments namely Tumbes (TU, n = 26), Piura (PI, n = 3), Lambayeque (LA, n = 5), La Libertad (LL, n = 48), Ancash (AN, n = 25), and Lima (LI, n = 36). Sampling localities were chosen due to the more diverse marine ichthyofauna present in the north than in the south [6]. Packages and labels from all packed items were kept for further examination. For samples collected from restaurants, we checked menus and asked the wait staff twice about the name of the marine species served to confirm each seafood type. In some instances, when the wait staff was not well informed, we requested information directly from the restaurant manager or chef. When possible, we targeted high priced menu-listings, which are more prone to be substituted by cheaper species. In most cases, we took pictures of the menu list and served dishes.

DNA extraction, PCR amplification, and sequencing

DNA extraction.

A small section of muscle was excised from the inner part of all collected samples, except for samples SF2, SF27, SF73, and SF74 from which a piece of dorsal fin was collected. Tissues were rinsed with distilled water and preserved in 99% ethanol at -20°C for further DNA analysis. Genomic DNA was isolated using the cetyl-trimethylammonium bromide (CTAB) precipitation method [35] for muscle tissues and the standard phenol-chloroform protocol [36] for fin tissue samples. To rehydrate the tissue and to remove contaminants, the processed samples were soaked in distilled water prior to DNA extraction.

Full barcoding PCR amplification.

Aiming to identify the largest number of seafood species that came from a variety of retailers and processors including a wide range of taxonomic groups, 10 different primer sets were utilized including FB from the mitochondrial COI and 16S rRNA genes, and the control region, and MB from the COI and 12S rRNA genes (Table 1). PCR amplifications for the COI gene were carried out using Folmer primers LCO1490/HCO2198 [17] for fishes and invertebrates, FishF1/FishR1 [37] and “cocktail” [38] primer sets for fishes, and a degenerated version of Folmer primers COIF-ALT/COIR-ALT designed for Veneridae [39] was used for surf clams. For scallop samples, we used the Pectinidae family-specific primer set Pect16BC [18], targeting the 5’ end of the mitochondrial 16S rRNA gene. For marlin samples, we used the primer set A/G [40], flanking the complete mitochondrial control region. All primer sets used in this study are listed in Table 1.

thumbnail
Table 1. PCR primer sets used in the amplification of samples analyzed in this study.

https://doi.org/10.1371/journal.pone.0206596.t001

PCR reactions were performed in a TurboCycler Blue-Ray (BioTech, Taipei, Taiwan), a ProFlex PCR System (Applied Biosystems, Foster City, CA, USA) or a Veriti 96 Well thermal cycler (Applied Biosystems, Foster City, CA, USA) using Maximo Taq DNA Polymerase 2X-preMix (GeneOn GmbH, Nurnberg, Germany), HotStarTaq Plus Master Mix (QIAGEN, Hilden, Germany) or Maximo Taq DNA Polymerase (GeneOn GmbH, Nurnberg, Germany). Most PCR amplifications using either Folmer or FishF1/FishR1 primer sets were performed using the following protocol: 1–2 μl of genomic DNA, 10 μl of Maximo Taq DNA Polymerase 2X-preMix (GeneOn GmbH, Nurnberg, Germany), 0.5 μM each primer, in 20 μl of total volume. Thermocycling conditions were as follows: initial denaturation for 5 min at 95°C, followed by 35 cycles of denaturation for 30 s at 95°C, annealing for 40 s at 54°C (FishF1/FishR1) or at 45°C (Folmer), and extension for 1 min at 72°C, followed by a final extension for 10 min at 72°C. Master mix and PCR protocols using the other two abovementioned PCR kits are detailed in S1 Table. The PCR amplification using the “cocktail” primer set was performed with the HotStarTaq Plus Master Mix (QIAGEN, Hilden, Germany) using PCR conditions described previously [38] with slight modifications: 38 amplification cycles and annealing temperature from 50 to 52°C (see S1 Table). Primer sets COIF-ALT/COIR-ALT and A/G (control region) were amplified with HotStarTaq Plus Master Mix (QIAGEN, Hilden, Germany) and the amplification protocols are also detailed in S1 Table. The Pectinidae primer set Pect16BC was amplified following the amplification protocol described in Marín et al. [18]. All PCR products were electrophoresed in a 1.5% agarose gel and visualized under UV light.

Mini-barcoding PCR amplification.

Samples that failed PCR amplification with COI “full barcoding” (hereafter FB) primer sets were amplified using three mini-barcoding (hereafter MB) primer sets, namely Mini_SH-A, Mini_SH-D, and Mini_SH-E [23]. PCR amplification reactions were the same as described above. PCR amplification conditions were as follows: initial denaturation for 5 min at 94°C, followed by 35 cycles of denaturation for 30 s at 94°C, annealing for 30 s at 46°C (Mini_SH-A and Mini_SH-E) or at 50°C (Mini_SH-D), and extension for 30 s at 72°C, followed by a final extension for 10 min at 72°C. In addition, a primer set namely MiFish-U [21] designed for fish metabarcoding environmental DNA (eDNA) that targets a small fragment (163–185 bp) of the hypervariable region of the 12S rRNA gene was used as an MB marker. PCR protocols for MB sets Mini_SH-A, Mini_SH-D, Mini_SH-E, and MiFish-U using HotStarTaq Plus Master Mix (QIAGEN, Hilden, Germany) are detailed in S1 Table.

Sequencing.

Positive amplification products were sequenced in both directions at Macrogen Inc. sequencing facilities (Korea) on an ABI 3730xl Genetic Analyzer (Applied Biosystems, Foster City, CA) and at the Laboratory of Molecular Genetics of IMARPE (Peru) on an ABI 3500 (Applied Biosystems, Foster City, CA). Sequencing primers for all markers used in this study are indicated in S1 Table.

Data analyses.

All DNA sequencing electropherograms were manually checked and edited by removing ambiguous base calling and adapter “tail” sequences (when necessary) using MEGA 7 software [41]. Complementary strand sequences were aligned so that a contiguous consensus was obtained for each sample. To avoid the inclusion of putative nuclear copies of COI gene sequences (NUMTs), sequences were manually checked for indels and premature stop codons. Identification of DNA sequences at species level was accomplished using both the Barcode of Life Data System (BOLD, http://www.boldsystems.org) selecting “species level barcode records” database and the Basic Local Alignment Search Tool (BLAST) on the National Center for Biotechnology Information (NCBI, http://www.blast.ncbi.nlm.nih.gov/Blast.cgi) identification engine (BLASTn, highly similar sequences “megablast”). Since records deposited only in the BOLD database have been validated for both the DNA sequence and specimen data, we used this repository as our final criteria for identifying seafood species. Only sequences with a similarity index ≥98% were considered a valid match [42]. In cases of ambiguous results obtained from NCBI and BOLD databases, further phylogenetic analysis including DNA sequences from both databases was performed using the Neighbor-Joining (NJ) method with Kimura 2-parameter model (K2P) [43] and 1000 bootstrap replicates using MEGA 7 software, and the Bayesian analysis inference (BI) using MrBayes 3.2.6 [44]. The level of substitution saturation for COI datasets was evaluated using DAMBE 6 [45]. We used jModelTest 2 [46] under the Akaike and Bayesian information criterion (AIC and BIC) to find the best-fit model of evolution. Two runs were performed simultaneously, each with four Markov chains. The analyses were run for one or five million generations with sampling every 100 generations. The first 25% of the sampled trees were discarded as burn-in. Obtained phylogenetic trees were drawn using FigTree 1.4.2 program (http://tree.bio.ed.ac.uk/software/figtree/). The species names obtained by barcoding analyses were then compared to the corresponding common/market names included in the “List of main species from artisanal fish landings during 2017” (hereafter “FISHLANDINGS-2017 list”), kindly provided by the Artisanal Fishery Office at IMARPE, and also compared with the fish common names presented in Chirichigno and Cornejo (2001) [4]. Acceptable English market names were searched within the FDA Seafood List accessible from https://www.accessdata.fda.gov/scripts/fdcc/index.cfm?set=seafoodlist. Accepted marine scientific names were checked in The World Register of Marine Species (WORMS, available at http://www.marinespecies.org) and FishBase (available at http://www.fishbase.de) databases. For batoid classification, we followed the nomenclature proposed in Last et al. [47]. Furthermore, we checked the conservation status for each genetically identified species at the International Union for Conservation of Nature (IUCN Red List of Threatened Species) publically available from http://www.iucnredlist.org.

Results and discussion

Molecular identification performance

This report represents the first intensive effort to accurately identify a wide range of commercial seafood products across the Peruvian supply chain (from harvest to consumption) using molecular markers. Overall, 137 PCR products were obtained, which represents a PCR success rate of 96%. Unsuccessful PCR amplifications resulted from canned and cooked samples, most likely due to DNA degradation or inhibitors presented in the processed food. We would like to emphasize that our protocols identified all “cebiche” samples (n = 16) collected from restaurants. Cebiche is the Peruvian national dish and by far the most popular and the pride of the citizens of Peru (see Fig 1E), where seafood is marinated with lime, which could affect the effectiveness of DNA isolation challenging downstream applications.

thumbnail
Fig 1. Representative pictures of sampling sites analyzed in this study.

a. Fish Landing Site (FLS): guitarfish; b. Market (MK): Pacific menhaden; c. Supermarket chain (SMC); d. Multimarket (MM); e to i. Restaurant (RT): e. marinated seafood “cebiche”, f. spicy shellfish cream “picante de mariscos”, g. grilled octopus, h. fish and shellfish in “parihuela” soup, and i. fried Peruvian grunt.

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

A total of 131 of the 137 amplified PCR products resulted in high quality sequencing electropherograms (sequencing success rate 96%), enabling proper seafood identification of 121 (92.37%) and 10 (7.63%) samples to species and genus level respectively. Of the 131 DNA sequences, 128 samples (97.71%) showed sequence identity greater than the threshold value (≥98%). However, because of either unresolved taxonomy or short variability of the COI gene among congeners, 13 samples (including octopus, marlin, smooth-hound, tilapia, and tuna specimens) matched with more than one species within 98% identity cutoff. Phylogenetic analyses for octopus, marlin, and smooth-hound samples are described below in subsections “Restaurants”, S1 Appendix, and S2 Appendix, respectively. Only three samples (SF42, SF44, and SF117) fell below the species-level identification cutoff (due to lack of reference sequences) and consequently, only a genus level identification was possible. Three (SF117, SF119, and SF128) of the 115 specimens identified with FB (COI) did not match to any records in the BOLD database.

Mini-barcoding efficiency

The efficiency of all the MB primers tested in this study is presented in S2 Table. Species identities within the range of 98–100% for NCBI and 98.2–100% for BOLD databases were obtained in 15 tested fish samples, of which 14 could not be amplified with FB markers. Most of those included canned and cooked samples from SMCs and RTs such as fish fritters, tortilla, fried, and steamed. A canned tuna sample (SF68) was successfully amplified and sequenced by MB, but BLAST analysis showed multiple Thunnus species (T. atlanticus, T. orientalis, and T. albacares) matching at 98% identity. A low genetic distance among those tuna species hampers a species level identification through DNA barcoding, consequently multilocus approaches based on control region with ITS1 [48] and 12S with ND5 markers [21] have been recommended.

Among the three tested MB primer sets from Shokralla et al. [23], the primer sets SH-A and SH-E showed the highest PCR amplification (63% and 77% respectively) and sequencing success rates (100% and 90% respectively). All three primer sets (SH-A, SH-D, and SH-E) showed high sequence identity (above 98%, S2 Table), which allowed a correct identification of eight bony fishes, one shark, and one devil ray species. A shark filet (SF76) was identified only to the genus level due to lack of reference sequence. Interestingly, this study added elasmobranch species to the performance test of these MB primers with promising results, however, we tested only a limited number of shark and devil ray species. Further analyses using a wider range of elasmobranch species are needed to find out the efficacy of those MB primers in the identification of members of this fish group. Mini-barcoding primer sets SH-A, SH-D, and SH-E failed to amplify one canned herring sample SF77 (Marine Stewardship Council -MSC- certified). Therefore, we tested an additional primer set MiFish-U targeting a small fragment of the 12S rRNA gene [21]. This primer set showed a high performance by identifying sample SF77 as Clupea harengus with 100% identity (NCBI database), confirming its correct label information. Using the same primer set, a second canned sample (SF96) was also tested and identified as Sardina pilchardus with 100% identity (NCBI database). Sample SF96 was also identified with primer sets SH-A and SH-E with sequence similarities of 99.22% and 99.56%, respectively, in the BOLD database. Our results indicate that the primer set MiFish-U designed by Miya et al. [21] can be considered an alternative potential MB marker for the identification of canned fish samples. Overall, the MB approach applied in this study revealed four cases of misbranding, which are shown in S3 Table.

Species diversity

Even though our sampling was performed in a relatively small number of retailers, which included 21 RTs, eight MKs, seven stores from four SMCs, six FLSs, one WFM, one GS, and one MM, a highly diversified fauna including both marine (93.44%) and freshwater (6.56%) species was found over the 20-month sampling period. We identified 55 species (plus five genus-level taxa) covering 47 families, 24 orders, and six classes including Actinopterygii (45.03%), Chondrichthyes (36.64%), Bivalvia (6.87%), Cephalopoda (6.11%), Malacostraca (3.82%), and Gastropoda (1.53%) (Table 2). The most diverse group was Perciformes represented by 12 families, 16 genera, and 16 species, followed by Myliobatiformes with five families, six genera, and six species. All seafood DNA sequences obtained in this study have been deposited in GenBank/EMBL/DDBJ databases with accession numbers from MH194422 to MH194552 and from MK070511 to MK070524 (Table 3).

thumbnail
Table 2. Taxonomic classification of seafood diversity identified in this study.

Conservation status was retrieved from IUCN Red List. IUCN abbreviations: NE Not Evaluated, DD Data Deficient, LC Least Concern, NT Near Threatened, VU Vulnerable, EN Endangered. n: sample size. Identification to the genus level (ID% or Similarity% <98) is highlighted in bold.

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

thumbnail
Table 3. Species identification results of the 131 samples collected through the supply chain of the Peruvian fishery sector using full (FB) and mini-barcoding (MB), results are based in NCBI and BOLD databases.

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

Fish landing sites, wholesale fish markets, and markets.

A total of 51 samples were identified from FLSs (n = 21), WFM (n = 7) and MKs (n = 23), represented mainly by batoids (72%), sharks (100%), and bony fishes (35%), respectively (see S1 Fig). Among batoid specimens (rays and their relatives) collected from FLSs (n = 15) and MKs (n = 7), including one Chilean eagle ray specimen labeled as smooth-hound, we identified nine species belonging to seven families: Aetobatidae, Arhynchobatidae, Dasyatidae, Mobulidae, Myliobatidae, Narcinidae, and Rhinobatidae (Table 2); and one genus (Gymnura sp.) from Gymnuridae (see S3 Appendix). Three (Hypanus dipterurus, Myliobatis longirostris, and Sympterygia brevicaudata) of the nine identified species collected from different FLSs in TU (collection date January 2017) were not found in recent landing reports (from January 2014 to July 2018) from Tumbes (C. Luque, personal communication). Besides, the sample identified as shorttail fanskate S. brevicaudata (SF47, BOLD similarity 100%) was labeled as witch skate Rostroraja velezi; the latter species was included in landing records from Tumbes (C. Luque, personal communication). These findings highlight the necessity for implementation of periodic genetic monitoring programs across landing sites to support fisheries management and conservation efforts of batoid species.

In 2016, total “ray” landings from Peru reached 2440 metric tonnes (MT), with 72.34% of this production (1765 MT) destined to the domestic fresh fish market, while the remaining was used for cured fish production [8]. The giant devil ray M. mobular (Myliobatiformes: Mobulidae), which is the most landed mobulid species in northern Peru [49], was also the most abundant batoid and only mobulid species detected in this study, being found in three FLSs (n = 5), three MKs (n = 4), and one RT (n = 1). In Peru, the genus Mobula comprises five species including M. birostris (formerly Manta birostris), M. tarapacana, M. mobular (formerly M. japanica [47]), M. thurstoni, and M. munkiana [50]. Predominantly gillnets are used by small-scale and industrial fisheries to target mobulids, but they are also caught as bycatch in the tuna purse seine fishery [49]. In spite of their commercial value, conservation concerns, and promising management and conservation efforts targeting chondrichthyan species (i.e., PAN Tiburón-Perú and a law prohibiting M. birostris fishery, see S4 Appendix and S5 Appendix), there have not been any molecular studies on Peruvian marine rays. Our results could be used to better understand the diversity of commercially important rays, providing baseline data for further genetic studies necessary to design and implement conservation actions.

The smooth hammerhead Sphyrna zygaena (Carcharhiniformes: Sphyrnidae), which is the third most commonly landed shark species in Peru [51], was found at the three retailer categories described in this section (FLS, WFM, and MK), where it was usually labeled as thresher shark. S4 Table shows all identified samples grouped by retailer category and seafood type. Nine samples were identified as S. zygaena with MB and FB (BOLD similarity 100%). Of particular concern was the detection of six specimens collected during closed season (January 1 to March 10, RM N° 008-2016-PRODUCE), including five headless samples (SF50 to SF54 from WFM-LL) labeled as thresher sharks (Alopias sp.) and one sample (SF66 from FLS-TU) landed as whole body. Illegal, unreported and unregulated (IUU) fishing can be profitable due to the high demand for overexploited and protected species, and low risk of getting caught or being punished, especially when it takes place in countries where enforcement is weak [52]. Sphyrna zygaena is categorized as “Vulnerable” by the IUCN Red List of Threatened Species. Regrettably, the shark fishery in Peru is poorly regulated and monitored, mainly because fisheries managers put more effort in controlling small pelagic resources which dominate the fishing industry [34].

Other interesting findings at FLSs were the landings of one whale shark specimen (SF46) Rhincodon typus (Orectolobiformes: Rhincodontidae) and one sample (SF65) identified as striped marlin Kajikia audax (Perciformes: Istiophoridae), both as whole bodies. The commercial fishing of both species has been banned by the Peruvian government (see text in S4 Appendix and S6 Appendix). The whale shark DNA sequence obtained herein represents the first publically available nucleotide sequence of R. typus (GenBank accession number MH194467) caught in Peruvian waters, which could be useful for further comparative studies, therefore the specimen tissue and DNA from sample SF46 are available upon request.

Among bony fishes (class Actinopterygii) collected from MKs, the economically valuable genus Anisotremus (Perciformes: Haemulidae) was represented by two species: the Peruvian grunt A. scapularis and the burrito grunt A. interruptus. In one MK from Ancash, we bought a bag containing 10 fish specimens labeled as “Peruvian grunt”, however one specimen was larger and darker than the others. Molecular analysis showed that the “dark grunt” (SF106) was actually arnillo drum Cheilotrema fasciatum (family Sciaenidae), which strongly resemble grunts in appearance but is of lower economic value. Another sample (SF105) from the same bag was identified as the Peruvian grunt A. scapularis. Apparently, the arnillo drum, larger in size, was put there to increase the total product weight.

Grocery store, supermarket chains, and multimarket.

A total of 35 samples were identified from SMCs (n = 22), MM (n = 12), and a GS (n = 1), (S4 Table). The only sample collected from GS was a canned tuna (SF68) identified as Thunnus sp. The most abundant group detected in both SMCs and MM was Actinopterygii with 55% and 66%, respectively (S1 Fig). Among the samples collected from SMCs and MM, the Peruvian anchovy E. ringens (canned and fish burger presentations), dolphinfish C. hippurus (fish burger and fish roe presentations), and Humboldt squid D. gigas (seafood mix presentation) represented some of the most important species from landings for direct human consumption during the year 2016 [8]. Two canned anchovy products (SF69 and SF70, SMC at LL) labeled as “Peruvian sardine” were identified as E. ringens. In 2009, the Ministry of Production of Peru, aiming to promote internal consumption of anchovy as well as to conquer new international markets, adopted the name “Peruvian sardines” for processed (i.e., canned) Peruvian anchovies [53]. This marketing strategy is due to the fact that in the international market, sardine is usually in higher demand than anchovy [54]. However, the presence of the local sardine species Sardinops sagax also known as “Peruvian sardine” [55] may cause confusion among local consumers.

Imported items, representing both farmed and wild species, were found only in SMC and MM retailers. During 2017, Peru imported 145344 MT of seafood products valued at US$306 million; frozen and canned products accounted for 70% of the total [56]. We authenticated nine imported seafood products belonging to six species from five orders including Perciformes (Kajikia albida, frozen filet n = 1), Clupeiformes (Sardina pilchardus and Clupea harengus, canned n = 3), Salmoniformes (Salmo salar, vacuum packed filet n = 1), Siluriformes (Pangasianodon hypophthalmus, frozen and vacuum packed filets n = 3), and Decapoda (Metapenaeus dobsoni, instant noodle soup n = 1). In 2015, Chilean port authorities detected a shipment (valued at US$19 million) from Callao (Peru) of 37200 cans of Pacific menhaden (Ethmidium maculatum) labeled as horse mackerel (Trachurus murphyi) [57]. In this regard, the use of DNA-based technologies for seafood authentication is imperative to ensure proper label information, not only for domestic and imported products but also for Peruvian exports.

Interestingly, in one supermarket from Lima, we found different imported canned products carrying the Marine Stewardship Council (MSC) blue ecolabel. The MSC is an international non-profit organization that sets a standard (MSC Fishery Standard) used to assess sustainable fisheries all over the world. Currently, there are no Peruvian fisheries holding MSC certification or undergoing full assessment. In 2015, a molecular barcoding test of a total of 256 MSC labeled products (from 16 countries, covering 13 fish species) performed by an independent laboratory revealed that 99.6% were correctly labeled [58]. Herein, we were able to verify the correct species information from one MSC certified canned herring sample imported from Germany and labeled as Clupea harengus (SF77, bought in SMC-LI) using the 12S rRNA gene eDNA metabarcoding primer set MiFish-U designed by Miya et al. [21].

We detected only one case of mislabeling in an SMC (from LL) in which a filet sample labeled as guitarfish (SF58) was determined to actually be Chilean eagle ray M. chilensis (BOLD similarity 100%). However, it is difficult to determine whether it was an intentional case of mislabeling due to the fact that SMCs usually rely on wholesale seafood distributors. It is important to mention that SMCs employ trained personnel and safety protocols including the application of cold chain management systems, thus preserving food quality and ensuring food safety. However, the aforementioned practices make SMC seafood products more expensive than those of popular MKs; sometimes the price difference is as high as 300% [59].

Three smooth-hound filet samples (SF56, SF57, and SF85) were first identified as Mustelus henlei (BOLD similarity 98.15–98.55%). However, NJ and BI phylogenetic analyses (Fig A in S2 Appendix) clustered those samples in a unique clade with high nodal support. Therefore, samples SF56, SF57, and SF85 were assigned to Mustelus sp. Another smooth-hound filet sample (SF75, SMC-LI) was identified as Mustelus lunulatus (BOLD similarity 100%). In Peru, smooth-hounds (Mustelus spp.), houndsharks (Triakis spp.), and catsharks (Schroederichthys spp.) are usually reported under the same common name “tollo” ([34], FISHLANDING-2017 list). Peruvian shark landing statistics at species level include only three Mustelus species: M. whitneyi, M. mento, and M. dorsalis [51]. A reduced frequency in landing occurrences, combined with low taxonomic resolution at landing sites and a poorly regulated and monitored fishery [34], may have been masking or “diluting” the presence of M. lunulatus from landing reports. However, we cannot rule out the possibility that sample SF75 was imported from Ecuador (where M. lunulatus also occurs), which is an important source of shark imports to Peru [51].

Inaccurate identification of morphologically similar species in combination with poor taxonomic resolution of fisheries landing reports and the application of inaccurate names will not only cause considerable economic impacts but also lead to undesired consequences for fishery management [60] including local population depletion. M. lunulatus is not included in the “Identification guide to commercially important sharks from Peru” [61], which is a field guide for identifying most frequent shark species in landings of artisanal fisheries of Peru. In this regard, government fishery officers must undergo training in detecting not only main commercial species but also the ones that are infrequently landed. Field identification guides should consider including the “less commercial” species.

One surprising result was the mislabeling of tilapia Oreochromis sp. as “sand-perch” (SF88, collected from MM, see Fig 1D). Passing off cheap farmed tilapia as more expensive wild fish has been reported in previous studies [62, 63]. During 2013, frozen imports of tilapia from China accounted for more than 50% of total domestic market sales [64]. China is Peru’s biggest trade partner, with investments exceeding US$14.00 billion [65]. A Free Trade Agreement (FTA) between Peru and China was signed on April 28, 2009, and entered into force on March 1, 2010 [66] bringing valuable opportunities for Peruvian entrepreneurs to go through Chinese markets duty-free. Unfortunately, imported Chinese tilapia enters the Peruvian market at significantly lower prices than locally produced ones [64].

Shellfish accounted for 18% of SMC samples, comprising three crustacean (14%) (M. dobsoni, P. vannamei, and Romaleon setosum), and one mollusk (4%) species (Humboldt squid D. gigas). Only mollusks (17%) were collected from MM, represented by two bivalve species (Tagelus dombeii and Donax obesulus). Two shucked shellfish samples (SF118 from MK-AN and SF127 from MM-LI) labeled as surf clam “macha” (Mesodesma donacium) were identified as jackknife clam T. dombeii (BOLD similarity 99.47–100%). Peruvian populations of the surf clam M. donacium have been depleted and its fishery abruptly collapsed, mainly due to the combined effects of unregulated overexploitation and adverse climatic events (El Niño/Southern Oscillation-ENSO) [67]. As evidenced by our results, the great demand that still exists for this bivalve species makes this resource vulnerable to substitution by other species in different Peruvian cities.

Restaurants.

The results of the present study provide a snapshot of species availability in some seafood restaurants. Forty-five samples were identified covering 15 bony fishes, nine shellfish, three sharks, and one batoid species (S4 Table). Within the bony fishes, which represented 60% of restaurant samples (S1 Fig), we identified high market-value species such as flounder (Paralichthyidae), grunts (Haemulidae), and grouper and seabass (Serranidae). In Peru, groupers (Epinephelus spp. and Mycteroperca spp.) and grape-eye seabass (Hemilutjanus macrophthalmos) are considered “luxury” seafood species and in high demand, which makes them more prone not only to overexploitation but also to substitution by cheaper ones. Two restaurant samples (SF7 and SF17) labeled as grouper were identified as South Pacific hake M. gayi and dolphinfish C. hippurus, respectively. Similarly, two samples (SF8 and SF9) labeled as grape-eye seabass were found to be swordfish Xiphias gladius.

When mislabeling occurs on board or at landing, the error continues along the food chain to the consumer [68]. Consequently, final seafood retailers such as restaurants are more vulnerable to receive mislabeled products. Herein, 17 (38%) of the 45 identified samples bought from 21 restaurants across 10 different districts were molecularly identified as different species to those declared by the restaurant staff or menu list (S3 Table). Similar substitution levels (from 26 to 50%) of samples collected in restaurants have been reported in other studies [13, 69, 70]. We want to emphasize, however, that mislabeling should not always be considered as fraud. Instead, it could be the result of species misidentification, due to the confusion generated by the use of different vernacular names in different regions or countries, or when mistakes occur during product information management by mid-chain players. Fortunately, many restaurateurs view seafood sustainability as a requisite for future viability, and some remarkable initiatives have been undertaken by Peruvian chefs engaging directly with artisanal fishers [71, 72].

Our results revealed that in cases of mislabeling, the species most commonly used as a replacement was dolphinfish (C. hippurus), which was served as grouper (SF17), corvine (SF23), sailfish (SF99, SF100, SF121), and cachema weakfish (SF130) in five different restaurants from La Libertad and Lima. The fact that dolphinfish meat is being used in some restaurants to replace other “white flesh” species including corvine drum was reported in a previous study [73] based on visual inspections (R. Gozzer, personal communication). The dolphinfish’s white flesh makes this species a potential substitute for others high-priced species. Peru is the main producer of dolphinfish with estimated landings accounting for more than 50% of global catches [73]. Peruvian dolphinfish fishery is targeted exclusively by the artisanal fleet, representing one of the nation’s most important artisanal fishery; however it is poorly regulated with high levels of informality along its supply chain [73].

Eight shellfish species including mollusks (octopuses, squids, scallops, mangrove cockle, and sea snail) and crustaceans (crab and shrimp) were identified from restaurant samples (S4 Table). Two species (whiteleg shrimp P. vannamei and Peruvian scallop A. purpuratus) account for most of Peruvian mariculture production [8]. Both are well known and widely accepted by consumers, becoming a target product for most seafood restaurants. In Peru, shrimp production has been growing steadily at about 10 percent annually since 2008 [56]. In 2017, shrimp production reached 26768 MT, of which 80% (21400 MT, valued at US$164.1 million) was exported [56]. On the other hand, Peruvian scallop production has decreased significantly from 67694 MT in 2013 to 13137 MT in 2017. An estimated 3300 MT (valued at US$54.3 million) was exported in 2017 [56]. The decline in scallop production was driven mainly by the “coastal El Niño”, which affected up to 98% of production in northern Peru during 2016–17 [74]. This drop in Peruvian scallop production has affected not only the domestic market but also global scallop trade [75].

Another valuable shellfish species is the Gould octopus Octopus mimus (Octopoda: Octopodidae), which supports an important artisanal fishery in Peru. Landing estimates were 5405 MT in 2016 [8]. A genetic study has suggested the possible conspecificity between O. mimus and the Hubb’s octopus O. hubbsorum [76]. The distribution of O. mimus is believed to be restricted from northern Peru to Chile, whereas O. hubbsorum is found from the Gulf of California to Oaxaca in Mexico [76]. However, some molecular studies have reported the presence of O. mimus in Central America and Ecuador [76, 77] (and references therein).

We analyzed four octopus samples collected from RTs in LL, AN, and LI. Barcoding results matched samples SF72, SF81, and SF124 to O. mimus and sample SF14 to O. hubbsorum (BOLD similarity 100%). Our phylogenetic results (Fig 2) showed congruent topologies between BI and NJ trees, with samples SF72, SF81, and SF124 clustered within the O. mimus subclade (BA posterior probability 95%, NJ bootstrap support 77%) with a maximum of 0.2% (K2P) within-cluster divergence, while SF14 was within the O. hubbsorum subclade (BA posterior probability 97%, NJ bootstrap support 45%) showing a maximum within-cluster divergence of 0.3% (K2P). The minimum genetic distance (K2P) between both subclades was 0.7%. Interestingly, sample SF14 shares the same haplotype with the O. mimus specimen reported in Ecuador (GenBank accession KT335830) [77]. Data mined from customs information imports from the National Customs Superintendency of Peru (SUNAT) (http://www.aduanet.gob.pe/cl-ad-itconsultadwh/ieITS01Alias?accion=consultar&CG_consulta=2), showed that recent octopus imports were represented only by O. mimus from Chile (from 2014 to 2017) and O. vulgaris from the Philippines in 2014. Thus, without further evidence of O. hubbsorum imports or a recent range expansion towards the South Pacific, we assigned sample SF14 to Octopus mimus. Further studies must be carried out to solve the taxonomic status of the economically important O. mimus, which is still under debate.

thumbnail
Fig 2. Phylogenetic tree based on Bayesian inference (BI) and Neighbor-Joining (NJ) for the identification of samples SF14, SF72, SF81, and SF124 Octopus mimus.

Phylogenetic tree based on COI barcode sequences (576 bp) from samples SF14, SF72, SF81, and SF124 (Octopus mimus, this study) and other Octopus reference sequences available in BOLD and NCBI. Sample SF14 is highlighted in blue and shaded in yellow. Samples SF72, SF81, and SF124 are highlighted in green and shaded in orange. Bayesian consensus tree was inferred with five million generations under the GTR+I+G substitution model. NJ tree was constructed with 1000 bootstrap replicates under the Kimura-2-parameter (K2P) model. Nodal supports for Bayesian inference posterior probabilities and bootstrap values for NJ analysis (highlighted in bold) above 45% are shown. Samples from this study include identification code and GenBank accession numbers. Reference sequence labels include BOLD process ID and GenBank accession numbers. Vampire squid Vampyroteuthis infernalis and North Atlantic octopus Bathypolypus arcticus were used as outgroup.

https://doi.org/10.1371/journal.pone.0206596.g002

Mislabeling

Overall, 35 (26.72%) out of the 131 identified samples were found to be mislabeled, the majority came from markets and restaurants. As expected, most (94.28%) of the misrepresented samples were processed or cooked, where morphological features had been altered or removed. Mislabeled samples included 15 bony fishes (42.86%), 13 sharks (37.14%), four batoids (11.43%), and three mollusk (8.57%) specimens. S3 Table summarizes all mislabeled samples detected in each retailer category. Except for grocery store (GS), where only one sample was collected, we detected at least one mislabeling in each sampling site: FLS (n = 1, 4.76%), WFM (n = 5, 71.42%), MK (n = 8, 34.78%), SMC (n = 1, 4.54%), MM (n = 3, 25%), and RT (n = 17, 37.77%). However, it must be mentioned that sample sizes from WFM (n = 7) and GS (n = 1) were not representative, which prevented us from making any meaningful inferences on mislabeling rates from those retailer categories.

Nevertheless, our results could be used as a starting point to identify the major mislabeled species and the most common substitute species, as well as high priority stages for species substitution control along the supply chain. For example, among chondrichthyans, the most used species was the Chilean eagle ray M. chilensis, which was incorrectly labeled as manta (SF29-MK), guitarfish (SF58-SMC), and smooth-hound (SF62-MK). A high mislabeling rate (54.17%) was found among all shark samples (n = 24) collected across six different stages of the supply chain, with a total of 13 mislabeling cases involving two “Vulnerable” (S. zygaena and A. pelagicus) and one “Near Threatened” (Prionace glauca) species. This result could be a consequence of inaccurate species identification practices at early stages of the supply chain coupled with weak enforcement of shark regulations within the seafood sector. In Peru, the shark fishery is poorly monitored, worsened by the superficial taxonomic identification at landing sites across the country [34]. The bony fish most commonly used to replace other species was dolphinfish C. hippurus, which was served as grouper, corvine, sailfish, and cachema weakfish in six cases involving restaurant samples. Our results indicate that mislabeling is a common issue within the Peruvian seafood sector. Markets and restaurants accounted for the most cases of mislabeling, making those retailer categories potential candidates to be considered as priority control stages. However, further studies covering wider geographical areas with larger sample sizes from each supply chain stage are needed to support our mislabeling results.

The use of the same common name for different species or a single species having different vernacular names even within the same region or country is a common issue in seafood labeling [78]. In Peru, seafood commercial names used for more than one species include “mero” (e.g., groupers Epinephelus spp., Mycteroperca spp., Alphestes spp.), “lenguado” (e.g., flounders Paralichthys spp., Etropus spp.), “tollo” (e.g., smooth-hounds Mustelus spp., Triakis spp., catsharks Schroederichthys spp.), “ojo de uva” (i.e., grape-eye seabass Hemilutjanus macrophthalmos and mocosa ruff Schedophilus haedrichi), “almeja” (e.g., clams Semele spp., Gari spp.), and “barquillo” (e.g., Chiton spp., Acanthopleura spp.), just to mention a few. Despite this large number, few efforts have been made to regulate the application of standardized Peruvian commercial fish names [79]. To standardize the nomenclature used for seafood products, official guides with acceptable market names were published by the US Food and Drug Administration (FDA) in 1988 [80] and the European Union (EU) in 2001 [81]. To avoid ambiguities with accepted market names within the Peruvian seafood sector, the creation of an official standardized list of commercial fish names is strongly encouraged.

Conservation status and regulatory framework

A revision of the conservation status using the IUCN Red List of Threatened Species [82] showed that among all identified species, four samples belonged to species classified as “Endangered”: R. typus (n = 1) and P. hypophthalmus (n = 3); 13 shark and one marlin samples came from four species classified as “Vulnerable”: S. zygaena (n = 9), Isurus oxyrinchus (n = 1), A. pelagicus (n = 3), and K. albida (n = 1); and seven samples came from three species listed as “Near Threatened”: K. audax (n = 1), P. glauca (n = 5), and M. longirostris (n = 1). We should mention that the Endangered P. hypophthalmus production come from large-scale farms. The remaining identified samples to the species level (n = 86, 65.65%) correspond to species listed as “Least Concern” (n = 44), “Not Evaluated” (n = 23), and “Data Deficient” (n = 19).

As aforementioned, M. japanica and M. mobular belong to the same species [47], with nomenclatural priority given to M. mobular [83]. Mobula japanica was considered a wide-ranging circumtropical species assessed as “Near Threatened”, whereas M. mobular was considered a Mediterranean endemic with a Red List Assessment of “Endangered” [84]. The lumping of both species represents a change in “taxonomic concept” requiring a reassessment for the Red List, which is scheduled as part of the Global Shark Trend’s pelagic species project in 2018 [84].

In Peru, the Ministry of Production (PRODUCE) through its Vice-ministry of Fisheries is the entity responsible for the establishment and application of fisheries management regulations. The legal framework that regulates fishing activities aiming to ensure preservation and the sustainable exploitation of the aquatic resources comprises the General Fisheries Act (DL N° 25977), its Regulations on the General Fisheries Act (DS N° 012-2001-PRODUCE, modified by DS N° 015-2007-PRODUCE), and the Control Regulation and Sanction of the Fishing and Aquaculture Activities (DS N° 017-2017-PRODUCE).

A summary of the most important and recent regulations related to the three threatened species groups detected in this study (i.e., sharks, mobulids, and istiophorids) is presented in S4 Appendix, S5 Appendix, and S6 Appendix. The regulatory framework related to labeling of manufactured products detected herein is described in S7 Appendix.

Conclusions

This study represents the first attempt to assess the biodiversity present across different stages of the Peruvian supply chain using full and mini DNA barcoding, providing baseline data on the incidence of major mislabeled and the most common substitute species within the Peruvian seafood sector. Our results showed that full and mini-barcoding approaches are reliable and useful tools for species diversity determination, authentication and mislabeling detection of seafood products traded in the Peruvian market, which includes a wide range of taxonomic groups. A current drawback is the lack of barcoding reference sequences of some economically important Peruvian seafood species including smooth-hounds M. dorsalis and M. whitneyi, and butterfly ray Gymnura afuerae. In this light, the generation of a comprehensive Peruvian seafood barcoding library based on a mass genetic profiling of seafood biodiversity will be helpful to overcome these disadvantages. A major effort on seafood traceability must be undertaken by governmental agencies, fishery policy makers, and scientists to protect treasured marine species such as those on the IUCN Red List (e.g., endangered whale shark and vulnerable hammerhead shark) and to detect illegal fishing during closed seasons. The molecular evidence presented in this study suggests that illegal, unreported and unregulated (IUU) fishing activities are occurring in some areas of the Peruvian seafood sector as well as fraudulent actions within the supply chain. Peruvian artisanal fisheries lack of basic information for their proper management, with no good records on commercial fisheries landings, and no proper monitoring of seafood along the supply chain [73]. Albeit illegal incidental or opportunistic catches of threatened marine species have been already reported by Peruvian governmental inspectors and researchers, that is, however, only the tip of the iceberg, compared with what is really slipping through the net. Action plans for implementing standard and emerging DNA technologies, including rapid molecular detection techniques and environmental DNA (eDNA) to monitor endangered and heavily exploited species must be a priority concern. To strengthen traceability, strict enforcement of fish inspection programs based on DNA barcoding throughout the seafood industry and retailers must be conducted by government agencies. DNA barcoding will help to prevent and combat illegal or “pirate” fishing, especially in a mega-diverse country with high fish consumption such as Peru.

Supporting information

S1 Table. PCR mix composition and amplification conditions for full and mini-barcoding primer sets using two different commercial master mix brands.

PCR and sequencing primers are indicated.

https://doi.org/10.1371/journal.pone.0206596.s001

(XLSX)

S2 Table. PCR and sequencing efficiency of mini-barcoding primers.

PCR Polymerase chain reaction, SEQ Sequencing, P Positive result, N Negative result, MB mini-barcoding, FB fullbarcoding.

https://doi.org/10.1371/journal.pone.0206596.s002

(XLSX)

S3 Table. Summary list of all mislabeled samples.

Mislabeled samples are grouped by retailer category and location.

https://doi.org/10.1371/journal.pone.0206596.s003

(XLSX)

S4 Table. Molecular identification of 131 seafood samples.

Samples are grouped by origin (retailer category) seafood type, and retailer location. Samples code highlighted in bold and denoted with an asterisk (*) correspond to mislabeled samples. English/Spanish ("declared as"), identification results from BLAST and BOLD analysis, GenBank accession numbers and nucleotide consensus sequences generated in this study are given for each sample. Genera highlighted in bold correspond to samples identified to the genus level.

https://doi.org/10.1371/journal.pone.0206596.s004

(XLSX)

S1 Fig. Pie charts depicting percentage of seafood type contributions in each retailer category.

Wholesale fish markets (WFM, n = 7, only shark samples) and grocery store (GS, n = 1, tuna sample) categories are not included.

https://doi.org/10.1371/journal.pone.0206596.s005

(PDF)

S1 Appendix. Phylogenetic identification results of samples SF65 striped marlin Kajikia audax and SF131 Atlantic white marlin K. albida.

https://doi.org/10.1371/journal.pone.0206596.s006

(PDF)

S2 Appendix. Phylogenetic identification results of samples SF56, SF57, and SF85 smooth-hound Mustelus sp.

https://doi.org/10.1371/journal.pone.0206596.s007

(PDF)

S3 Appendix. Phylogenetic identification results of samples SF42 and SF44 butterfly ray Gymnura sp.

https://doi.org/10.1371/journal.pone.0206596.s008

(PDF)

S4 Appendix. Conservation status and regulatory framework of sharks.

https://doi.org/10.1371/journal.pone.0206596.s009

(PDF)

S5 Appendix. Conservation status and regulatory framework of mobulids.

https://doi.org/10.1371/journal.pone.0206596.s010

(PDF)

S6 Appendix. Conservation status and regulatory framework of istiophorids.

https://doi.org/10.1371/journal.pone.0206596.s011

(PDF)

S7 Appendix. Regulatory framework related to labeling of manufactured products.

https://doi.org/10.1371/journal.pone.0206596.s012

(PDF)

Acknowledgments

We thank John Bower (Hokkaido University, Japan) for improving the use of English in the final version of the manuscript. The authors are also grateful to Renato Gozzer (WWF, Peru) for his useful comments on the Peruvian dolphinfish fishery and to Carlos Luque (IMARPE-Tumbes) for his help during the comparison of batoid species detected in this study with the ones listed in landing records from Tumbes. We also thank Ximena Velez-Zuazo (CCS-Smithsonian Institution) for her kind support in comparing smooth-hound (Mustelus sp.) DNA sequences from her previous work with sequences determined in this study.

References

  1. 1. Chirichigno N, Cornejo RM. Catálogo comentado de los peces marinos del Perú. Lima: Instituto del Mar del Perú. 2001; 314 pp.
  2. 2. Ramírez R, Paredes C, Arenas J. Moluscos del Perú. Rev Biol Trop. 2003; 51(3): 225–284.
  3. 3. Food and Agriculture Organization of the United Nations (FAO). Perfíles de Pesca y Acuicultura por Países. Perú. Hojas de datos de perfiles de los países. In: Departamento de Pesca y Acuicultura de la FAO. Roma. Actualizado 1 May 2010. Available from: http://www.fao.org/fishery/ Accessed 3 May 2018.
  4. 4. Arellano CE, Swartzman G. The Peruvian artisanal fishery: changes in patterns and distribution over time. Fish Res. 2010; 101(3): 133–145.
  5. 5. Oyarzún D, Brierley CM. The future of coastal upwelling in the Humboldt current from model projections. Clim Dyn. 2018; 1–17.
  6. 6. Swartzman G, Bertrand A, Gutiérrez M, Bertrand S, Vasquez L. The relationship of anchovy and sardine to water masses in the Peruvian Humboldt Current System from 1983 to 2005. Prog Oceanogr. 2008; 79(2–4): 228–237.
  7. 7. Food and Agriculture Organization of the United Nations (FAO). The State of World Fisheries and Aquaculture 2018—Meeting the sustainable development goals. 2018. Rome. 210 pp.
  8. 8. PRODUCE. Anuario estadístico pesquero y acuícola 2016. 2018. Lima. 202 pp.
  9. 9. Centre for the Promotion of Imports (CBI). Peru’s potential as a seafood exporter. 10 Oct 2017. Available from: https://www.cbi.eu/market-information/fish-seafood/peru-potential-seafood-exporter/ Accessed 3 May 2018.
  10. 10. Food and Agriculture Organization of the United Nations (FAO). The State of World Fisheries and Aquaculture 2016 –Contributing to food security and nutrition for all. 2016. Rome. 200 pp.
  11. 11. Morales LE, Higuchi A. Is fish worth more than meat?–How consumers’ beliefs about health and nutrition affect their willingness to pay more for fish than meat. Food Qual Prefer. 2018; 65: 101–109.
  12. 12. Marín A, Fujimoto T, Arai K. Rapid species identification of fresh and processed scallops by multiplex PCR. Food Control. 2013; 32: 472–476.
  13. 13. Willette DA, Simmonds SE, Cheng SH, Esteves S, Kane TL, Nuetzel H, et al. Using DNA barcoding to track seafood mislabeling in Los Angeles restaurants. Conserv Biol. 2017; 31(5): 1076–1085. pmid:28075039
  14. 14. Pardo MA, Jiménez E, Pérez-Villarreal B. Misdescription incidents in seafood sector. Food Control. 2016; 62: 277–283.
  15. 15. Hajibabaei M, Singer G, Hebert PD, Hickey DA. DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics. Trends Genet. 2007; 23: 167–172. pmid:17316886
  16. 16. Naaum AM, Hanner R. Community engagement in seafood identification using DNA barcoding reveals market substitution in Canadian seafood. DNA Barcodes. 2015; 3(1): 74–79.
  17. 17. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994; 3: 294–299. pmid:7881515
  18. 18. Marín A, Fujimoto T, Arai K. The variable 5´ end of the 16S rRNA gene as novel barcoding tool for scallops (Bivalvia, Pectinidae). Fish Sci. 2015; 81: 73–81.
  19. 19. Lobo J, Costa P, Teixeira M, Ferreira M, Costa M, Costa F. Enhanced primers for amplification of DNA barcodes from a broad range of marine metazoans. BMC Ecol. 2013; 13:34. pmid:24020880
  20. 20. Vences M, Thomas M, Van der Meijden A, Chiari Y, Vieites DR. Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians. Front Zool. 2005; 2(1): 5. pmid:15771783
  21. 21. Miya M, Sato Y, Fukunaga T, Sado T, Poulsen JY, Sato K, Minamoto T, Yamamoto S, Yamanaka H, Araki M, Kondoh M, Iwasaki W. MiFish, a set of universal PCR primers for metabarcoding environmental DNA from fishes: detection of more than 230 subtropical marine species. R Soc Open Sci. 2015; 2(7): 150088. pmid:26587265
  22. 22. Teletchea F. Molecular identification methods of fish species: reassessment and possible applications. Rev Fish Biol Fish. 2009; 19(3): 265.
  23. 23. Shokralla S, Hellberg RS, Handy SM, King I, Hajibabaei M. A DNA mini-barcoding system for authentication of processed fish products. Sci Rep. 2015; 5: 15894. https://doi.org/10.1038/srep15894 pmid:26516098
  24. 24. Carvalho DC, Neto DAP, Brasil BSF, Oliveira DAA. DNA barcoding unveils a high rate of mislabeling in a commercial freshwater catfish from Brazil. Mitochondrial DNA. 2011; 22(S1): 97–105.
  25. 25. Palmeira CAM, da Silva Rodrigues-Filho LF, de Luna Sales JB, Vallinoto M, Schneider H, Sampaio I. Commercialization of a critically endangered species (largetooth sawfish, Pristis perotteti) in fish markets of northern Brazil: authenticity by DNA analysis. Food Control. 2013; 34(1): 249–252.
  26. 26. Ardura A, Linde AR, Moreira JC, Garcia-Vazquez E. DNA barcoding for conservation and management of Amazonian commercial fish. Biol Conserv. 2010; 143(6): 1438–1443.
  27. 27. Ardura A, Pola IG, Ginuino I, Gomes V, Garcia-Vazquez E. Application of Barcoding to Amazonian commercial fish labeling. Food Res Int. 2010; 43: 1549–1552.
  28. 28. de Brito MA, Schneider H, Sampaio I, Santos S. DNA barcoding reveals high substitution rate and mislabeling in croaker fillets (Sciaenidae) marketed in Brazil: the case of ‘pescada branca’ (Cynoscion leiarchus and Plagioscion squamosissimus). Food Res Int. 2015; 70: 40–46.
  29. 29. Rodrigues ADS, Brandão JHSG, Bitencourt JDA, Jucá-Chagas R, Sampaio I, Schneider H, Affonso PRADM. Molecular identification and traceability of illegal trading in Lignobrycon myersi (Teleostei: Characiformes), a threatened Brazilian fish species, using DNA barcode. ScientificWorldJournal. 2016; vol 2016 Article ID 9382613.
  30. 30. Bornatowski H, Braga RR, Vitule JRS. Shark mislabeling threatens biodiversity. Science. 2013; 340, 923.
  31. 31. Aguilera-Muñoz F, Valenzuela-Muñoz V, Gallardo-Escárate C. Authentication of commercial Chilean mollusks using ribosomal internal transcribed spacer (ITS) as specie-specific DNA marker. Gayana (Concepción). 2008; 72: 178–187.
  32. 32. Haye PA, Segovia NI, Vera R, Gallardo M, Gallardo-Escárate C. Authentication of commercialized crab-meat in Chile using DNA barcoding. Food Control. 2012; 25: 239–244.
  33. 33. Herrero B, Vieites JM, Espiñeira M. Authentication of Atlantic salmon (Salmo salar) using real-time PCR. Food Chem. 2011; 127(3): 1268–1272. pmid:25214125
  34. 34. Velez-Zuazo X, Alfaro-Shigueto J, Mangel J, Papa R, Agnarsson I. What barcode sequencing reveals about the shark fishery in Peru. Fish Res. 2015; 161: 34–41.
  35. 35. Nishiguchi MK, Doukakis P, Egan M, Kizirian D, Phillips A, Prendini L, Rosenbaum HC, Torres E, Wyner Y, De-Salle R, Giribet G. DNA isolation procedures. In: DeSalle R, Giribet G, Wheeler W, editors. Methods and Tools in Biosciences and Medicine: Techiques in molecular systematics and evolution. Birkhäuser Verlag; Basel, Switzerland. 2002. pp. 279–280.
  36. 36. Sambrook J, Russell DW. Molecular cloning: A laboratory manual (3rd ed.). Cold Spring Harbor, NY, USA: Cold Spring Harbor Press, 2001.
  37. 37. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PD. DNA barcoding Australia's fish species. Philos Trans R Soc Lond B Biol Sci. 2005; 360(1462): 1847–1857. pmid:16214743
  38. 38. Ivanova NV, Zemlak TS, Hanner RH, Hebert PD. Universal primer cocktails for fish DNA barcoding. Mol Ecol Resour. 2007; 7(4): 544–548.
  39. 39. Mikkelsen PM, Bieler R, Kappner I, Rawlings TA. Phylogeny of Veneroidea (Mollusca: Bivalvia) based on morphology and molecules. Zool J Linn Soc. 2006; 148(3): 439–521.
  40. 40. Lee WJ, Conroy J, Howell WH, Kocher TD. Structure and evolution of teleost mitochondrial control regions. J Mol Evol. 1995; 41(1): 54–66. pmid:7608989
  41. 41. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016; 33: 1870–1874. pmid:27004904
  42. 42. Cawthorn DM, Duncan J, Kastern C, Francis J, Hoffman LC. Fish species substitution and misnaming in South Africa: an economic, safety and sustainability conundrum revisited. Food chem. 2015; 185: 165–181. pmid:25952855
  43. 43. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980; 16(2): 111–120. pmid:7463489
  44. 44. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003; 19(12): 1572–1574. pmid:12912839
  45. 45. Xia X. DAMBE6: new tools for microbial genomics, phylogenetics, and molecular evolution. J Hered. 2017; 108(4): 431–437. pmid:28379490
  46. 46. Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012; 9(8): 772.
  47. 47. Last P, White W, de Carvalho M, Séret B, Stehmann M, Naylor G (Eds.). Rays of the World. Csiro Publishing. 2016.
  48. 48. Viñas J, Tudela S. A validated methodology for genetic identification of tuna species (genus Thunnus). PLoS ONE. 2009; 4(10): e7606. pmid:19898615
  49. 49. Alfaro-Cordova E, Del Solar A, Alfaro-Shigueto J, Mangel J, Diaz B, Carrillo O, Sarmiento D. Captures of manta and devil rays by small-scale gillnet fisheries in northern Peru. Fish Res. 2017, 195: 28–36.
  50. 50. Cornejo R, Velez-Zuazo X, Gonzalez-Pestana A, Kouri C, Mucientes GR. An updated checklist of Chondrichthyes from the southeast Pacific off Peru. Check List. 2015; 11(6): 1809.
  51. 51. Gonzalez-Pestana A, Kouri JC, Velez-Zuazo X. Shark fisheries in the Southeast Pacific: A 61-year analysis from Peru. F1000Research. 2016; 3:164.
  52. 52. Martini M. Illegal, unreported and unregulated fishing and corruption. Transparency International 2013; 392: 9. Available from: http://www.u4.no/publications/illegal-unreported-and-unregulated-fishing-and-corruption/ Accessed 4 May 2018.
  53. 53. Agencia Peruana de Noticias (ANDINA). Peruvian sardines promoted everywhere. 11 Set 2009. Available from: http://www.andina.com.pe/agencia/noticia-peruvian-sardines-promoted-everywhere-253257.aspx Accessed 3 May 2018.
  54. 54. Fréon P, Sueiro JC, Iriarte F, Miro OF, Landa Y, Mittaine JF, et al. Harvesting for food versus feed: a review of Peruvian fisheries in a global context. Rev Fish Biol Fish. 2013; 24: 381–398.
  55. 55. Cárdenas-Quintana G, Franco-Meléndez M, Salcedo-Rodríguez J, Ulloa-Espejo D, Pellón-Farfán J, Harris C. The Peruvian sardine, Sardinops sagax: Historical analysis of the fishery (1978–2005). Ciencias marinas. 2015; 41: 203–216.
  56. 56. PROMPERU (Comisión de Promoción del Perú para la Exportación y el Turismo). Informe anual: Desenvolvimiento del comercio exterior pesquero y acuícola en el Perú 2017. Departamento de Productos Pesqueros de la Sub Dirección de Promoción Internacional de la Oferta Exportable. 2018, 85 pp.
  57. 57. Declaration and labeling forgery detected in 37200 cans of alleged horse mackerel. Fish Info & Services Co. (FIS). 19 Oct 2015. Available from: http://fis.com/fis/worldnews/worldnews.asp?monthyear=&day=17&id=79907&l=e&special=&ndb=1%20target=Accessed 3 May 2018.
  58. 58. Marine Stewardship Council (MSC). From ocean to plate: How DNA testing helps to ensure traceable, sustainable seafood. March 2016. 2016; 1: 15. Available from: https://www.msc.org/docs/default-source/default-document-library/what-we-are-doing/msc-from-ocean-to-plate-traceability-and-dna-report-2016.pdf?sfvrsn=fbce98a4_4 Accessed 4 May 2018.
  59. 59. Gestión. Pescado en supermercados es 300% más caro que en mercados populares, afirma Elsa Galarza. 7 Feb 2014. Available from: https://gestion.pe/economia/pescado-supermercados-300-caro-mercados-populares-afirma-elsa-galarza-3444 Accessed 4 May 2018.
  60. 60. Fischer J. Fish identification tools: review and guidance for decision-makers. FAO Fisheries and Aquaculture Technical Paper No. 585. 2013 Rome, FAO. 107 pp.
  61. 61. Instituto del Mar del Perú (IMARPE). Guía para la determinación de tiburones de importancia comercial en el Perú. Serie de divulgación científica 2015. Año 1, Vol 1, N°2. Available from: http://biblioimarpe.imarpe.gob.pe:8080/handle/123456789/3007 Accessed 4 May 2018.
  62. 62. Jacquet JL, Pauly D. Trade secrets: renaming and mislabeling of seafood. Mar Policy. 2008; 32(3): 309–318.
  63. 63. Filonzi L, Chiesa S, Vaghi M, Marzano F. Molecular barcoding reveals mislabelling of commercial fish products in Italy. Food Res Int. 2010; 43(5): 1383–1388.
  64. 64. El Comercio. Más del 50% de tilapia que se vende en el Perú es importada. 6 Mar 2014. Available from: https://elcomercio.pe/economia/peru/50-tilapia-vende-peru-importada-167142 Accessed 3 May 2018.
  65. 65. ANDINA (Agencia Peruana de Noticias). China is Peru’s biggest trade partner. 26 Aug 2016. Available from: http://www.andina.com.pe/Ingles/noticia-china-is-perus-biggest-trade-partner-628207.aspx Accessed 3 May 2018.
  66. 66. Foreign Trade Information System (SICE OAS). Text Agreement China-Peru. 2018. Available from http://www.sice.oas.org/TPD/PER_CHN/PER_CHN_e.ASP
  67. 67. McLachlan A, Defeo O. Chapter 14 –Fisheries. In: McLachlan A, Defeo O, editors. The Ecology of Sandy Shores (Third Edition). Academic Press; 2018. pp. 331–374.
  68. 68. Gordoa A, Carreras G, Sanz N, Viñas J. Tuna Species Substitution in the Spanish Commercial Chain: A Knock-On Effect. PLoS ONE. 2017; 12(1): e0170809. pmid:28125686
  69. 69. Staffen CF, Staffen MD, Becker ML, Löfgren SE, Muniz YCN, de Freitas RHA, Marrero AR. DNA barcoding reveals the mislabeling of fish in a popular tourist destination in Brazil. PeerJ. 2017; 5: e4006. pmid:29201560
  70. 70. Christiansen H, Fournier N, Hellemans B, Volckaert F. Seafood substitution and mislabeling in Brussels’ restaurants and canteens. Food Control. 2018; 85: 66–75.
  71. 71. Advanced Conservation Strategies. A marine conservation assessment in Peru. Report prepared for the David & Lucile Packard Foundation & Fondation Ensemble. 2014. Available from: https://www.fondationensemble.org/wp-content/uploads/2015/02/ACS-Marine-Conservation-Assessment-of-Peru-final.pdf
  72. 72. PRODUCE: Restaurantes top de Lima compraron directamente a pescadores artesanales por mas de S/. 145000. Gestión. 11 Jul 2017. Available from: https://gestion.pe/economia/empresas/produce-restaurantes-top-lima-compraron-directamente-pescadores-artesanales-s-145-000-139167 Accessed 4 May 2018.
  73. 73. Amorós S, Gozzer R, Melgar V, Rovegno N. Peruvian mahi mahi fishery (Coryphaena hippurus): characterization and analysis of the supply chain. WWF- Marine Program of WWF-Peru. 2017. Available from: http://d2ouvy59p0dg6k.cloudfront.net/downloads/mahi_mahi_value_chain_en.pdf
  74. 74. Mereghetti M. European buyers face another year of scarce Peruvian scallops. UnderCurrentNews. 2 may 2017. Available from: https://www.undercurrentnews.com/2017/05/02/scarcity-of-peruvian-scallops-an-issue-for-european-buyers/
  75. 75. Globefish. Analysis and information on world fish trade. Less bivalves in world trade in 2016. 6 Jul 2017. Available from: http://www.fao.org/in-action/globefish/market-reports/resource-detail/en/c/903650/
  76. 76. Pliego-Cárdenas R, Hochberg FG, León FJGD, Barriga-Sosa IDLA. Close genetic relationships between two American octopuses: Octopus hubbsorum Berry, 1953, and Octopus mimus Gould, 1852. J Shellfish Res. 2014; 33(1): 293–303.
  77. 77. Pliego-Cárdenas R, Flores L, Markaida U, Barriga-Sosa I, Mora E, Arias E. Genetic evidence of the presence of Octopus mimus in the artisanal fisheries of octopus in Santa Elena Peninsula, Ecuador. Am Malacol Bull. 2016; 34(1): 51–55.
  78. 78. Barbuto M, Galimberti A, Ferri E, Labra M, Malandra R, Galli P, Casiraghi M. DNA barcoding reveals fraudulent substitutions in shark seafood products: the Italian case of “palombo” (Mustelus spp.). Food Res Int. 2010; 43(1): 376–381.
  79. 79. Rodríguez E, Paredes F, Zeballos J, Juarez L, Barreto J. Nomenclatura actualizada de peces comerciales del Perú. IMARPE. 1996; 1: 21–30.
  80. 80. US Department of Health and Human Services (USDHHS). ‘‘The fish list” FDA guide to acceptable market names for food fish sold in interstate commerce. Department of Health and Human Service, FDA Center for Food Safety and Applied Nutrition, available from Sup. Documents, US Govt. Printing Office, Washington, DC. 1988.
  81. 81. European Comission (EC). Comission Regulation (EC) N° 2065/2001. Laying down detailed rules for the application of Council Regulation (EC) No. 104/2000 as regards informing consumers about fishery and aquaculture products. Official Journal of the European Communities. 2001; L278: 6–8.
  82. 82. IUCN. The IUCN Red List of Threatened Species. Version 2018–1. Available from: http://www.iucnredlist.org Accessed 16 September 2018.
  83. 83. Poortvliet M. Genetics of manta and devil rays: Evolution, population genetics and conservation of a group of vulnerable pelagic filter-feeders. Rijksuniversiteit Groningen. 2015.
  84. 84. Kyne PM. It’s All in the Name: Shark Systematics and the IUCN Red List. IUCN SSC Shark Specialist Group Global Shark Trends Project 2018–2020. 2018. Available from: http://www.iucnssg.org/shark-systematics-and-the-iucn-red-list.html Accessed 21 August 2018.