1. Introduction
The collection of polychaetes in subtidal and intertidal mudflat habitats is a widespread activity worldwide, playing a key role in coastal economies [
1,
2,
3]. Live polychaetes are used for multiple proposes, from maturation diets for crustacean and finfish broodstock in aquaculture facilities [
4], to perform bioremediation processes of marine fish farms effluents [
5,
6], and as live baits for recreational and commercial fishing [
2]. Several polychaetes species are exploited for these purposes, such as
Arenicola marina,
Atilla virens,
Diopatra neapolitana,
Glycera dibranchiata,
Halla parthenopeia,
Hediste diversicolor,
Marphysa sanguinea,
Namalycastis rhodochorde,
Ophelia neglecta,
Perinereis linea,
Scoletoma impatiens and
Sigalion squamosus [
1,
2]. Over the past decades,
D. neapolitana has been one of the main polychaete species harvested in Ria de Aveiro, a coastal lagoon in mainland Portugal, as well as in other estuaries along the Portuguese coast line (e.g. Tagus estuary, Sado estuary and the coastal lagoon Ria Formosa) [
2,
7,
8,
9,
10]. The overexploitation of this resource has prompted a decline in local populations [
1] and, in an attempt to reduce the impact of the harvesting activity, several measures have subsequently been implemented. According to the Portuguese legislation (Portaria nº 1228/2010) bait gatherers can only operate with a personal license and are only allowed to work using hand gathering or a restricted gear. A more recent Ordinance was published in January 2014 (Portaria nº 14/2014) in an attempt to define the maximum daily catch limit. The daily limit is assumed to reflect the Maximum Sustainable Yield (MSY) of digging, ensuring a sustainable income to diggers.
According to this Ordinance, the daily catch limit for annelids should be 0.5 L day−1 per digger, excluding the tube [
11]. In 1999, a commercial bait harvesting value of around 200 million euros was estimated for Europe but, the existing gaps in the supervision and regulation of this activity, could result in an underestimation of the real value [
12]. In 2001/2002, the harvesting of
D. neapolitana was quantified for the first time in Ria de Aveiro, with an estimate annual catch volume of around 45 tons and a commercial value of approximately 350 thousand euros [
13]. Nevertheless, most commercial diggers, at a national and international level, do not have valid licenses, neither properly report their harvesting, as such, they contribute to foster to a parallel economy by not declaring their sales for tax purposes and impair any management plan ruling this activity to be successfully implemented [
12,
14]. The commercial potential of live fishing bait is so high that several attempts were made to intensively culture
D. neapolitana [
15,
16]. Nonetheless, several constraints still hamper this approach and the harvesting of specimens from wild populations remains the sole supply source of this highly priced polychaete species.
The effects of bait harvesting on species density and population structure has a direct impact on the community (neglecting non-target species) and, consequently, on ecosystem functioning and processes [
17,
18]. The polychaetes population response to harvesting is influenced by the nature and extent of bait digging pressure and by the demographic characteristics of the population being exploited [
14]. The current status of this species and the sustainability of its harvest are issues of growing concern [
11]. As such, reliably tracing the geographic origin of polychaetes harvested from the wild could be paramount for a sustainable management plan for bait harvesting, that could foresee no take periods and areas, as well as expose illegal poaching activities by less scrupulous bait diggers.
At present, several tools have been applied to confirm the geographic origin of marine organisms [
19,
20]. Fatty acid (FA) profiles recorded in soft tissues (e.g. adductor muscle of bivalves), have been successfully used for this purpose in common cockles (
Cerastoderma edule) [
21,
22] and manila clams (
Ruditapes decussatus) [
23]. The particular physicochemical conditions of each ecosystem shape FA composition of marine organisms in the sense that salinity and temperature are known to modulate the structure, fluidity and thus the composition of cell membranes [
24]. Higher salinity fluctuations and/or lower water temperatures promote a decrease in the levels of saturated FA (SFA) and an increase in the concentration of polyunsaturated FA (PUFA), responsible for the stabilization of the bilayer structure [
24]. Elemental fingerprints (EF) use the elemental profile recorded in hard biogenic structures, such as shells of common cockles (
C. edule) [
25,
26], otoliths of California halibut (
Paralichthys californicus) and Garibaldi (
Hypsypops rubicundus) [
27] and bony plates of Long-Snouted Seahorse (
Hippocampus guttulatus) [
28]. Considering that elements are influenced by the environmental and chemical features of each ecosystem [
29] and that these mineral structures grow throughout the year, EF have already been successfully used to discriminate specimens originating from geographically close locations [
30,
31,
32].
The combination of two different traceability tools can be more effective to confirm the geographic origin of marine organisms, as already confirmed by Zhang, et al. [
33] and Perez, et al. [
34] that showed that the use of stable isotope ratios in combination with FA profiles could successfully discriminate scallop species (
Patinopecten yessoensis,
Chlamys farreri, and
Argopecten irradians) and warty venus (
Venus vecurrosa), respectively, from different geographic locations. This combination of tools was also able to trace the geographic origin and seasonality of the whitemouth croaker (
Micropogonias furnieri) [
35]. Moreover, Matos, et al. [
36], using stable isotopes and elemental fingerprints effectively traced the geographic origin of eastern oysters (
Crassostrea virginica).
In an attempt to contribute to a better management of D. neapolitana stocks, the present study tested, for the first time, if the combination of FA profiles of the whole polychaete body and the EF of its jaws differs between specimens originanting from different locations in Ria de Aveiro a coastal lagoon in mainland Portugal where the capture of this polychaete being used as live bait for sports fishing is an important economic and social activity.
4. Discussion
The use of FA profiles of soft tissues and EF of mineral structures in marine species has been optimized in order to put forward faster, more accurate and environmentally safer methods to discriminate the geographic origin of these organisms, namely those that feature an important commercial value [
23,
28,
43,
44,
45,
46,
47,
48]. Most of available studies to date using these methods, either when these are applied individually or combining more than one approach (e.g., FA profiles and stable isotopes [
34], elemental fingerprint and stable isotopes [
36]), are mostly focused on food safety issues, rather than in the implementation of effective management plans for endogenous marine resources that may be vulnerable to poaching. The present study showed, for the first time, that the combination of biogeochemical tools (FA profiles and EF) using the whole body and the jaws of
D. neapolitana (respectively) can successfully be used to confirm the geographic origin of bait digging with a high accuracy level (
Table 1). It is therefore legitimate to say that these natural barcodes can be successfully used towards the implementation of more effective fishery management plans, allowing the enforcement of no take zones.
The FA profiles displayed by the whole body of
D. neapolitana revealed as the most dominant FA the 16:0, followed by the PUFA 20:5n-3. In general, this trend was similar to that found for other polychaetes by Fernandes, et al. [
49] for
Hediste diversicolor and by Jerónimo, et al. [
50] for
H. diversicolor,
D. neapolitana, Sabella cf.
pavonina and
Terebella lapidaria, all collected in the same location in Ria de Aveiro. The same was found for
Alvinella pompejana [
51],
Arenicola marina [
52],
Nepthys hombergii and
Lanice conchilega [
53] sampled in others locations.
The FA belonging to the SFA and PUFA classes were responsible for most of the differences recorded among locations (p < 0.05;
Table S1 on supplementary data). These dissimilarities in FA profiles of polychaetes among locations were likely associated with a differential physiological response to changes in environmental conditions (e.g. salinity and temperature) that shape the environment on their sampling locations [
22,
49,
54,
55]. At higher temperatures, the reorganization of the membrane structure is needed to maintain membrane fluidity and homeostasis, leading to an increase in FA saturation or the prevalence of shorter-chain FA [
56]. Salinity is responsible for changes in FA profiles involved in the osmoregulation process, inducing changes in membrane permeability [
57]. Higher salinities, like those registered at M1, I and E are associated with a decrease in PUFA to reduce membrane permeability [
58].
The biogenic carbonate hard parts of marine species, such as shells, otoliths, plates, fish scales and fish bones, incorporate and retain elements from the surrounding environment throughout their lifetimes [
28,
47,
59,
60,
61,
62]. It is important to note that this study represents the first dataset reporting the EF of polychaete jaws. Similar to other mineral structures, (e.g.
Cerastoderma edule shells [
26,
63]), the levels of elements recorded in aragonitic jaws of
D. neapolitana [
64] differed among locations within Ria de Aveiro. These discrepancies are associated with both local physical conditions and elements availability in the environment. Considering that location M2 is more upstream than M1 (within the same channel and closer to the inlet), hydrological conditions differ, which result in a substantial enrichment of Ba and Mn upstream due to freshwater inputs and nutrient runoff, as suggested by Ricardo et al., [
26]. In fact, the presence of high levels of Ba and Mn in location M2 had already been previously reported for these authors in
C. edule shells [
26]. The highest concentrations of K and Na observed in M1 could be related with its geographical proximity to the inlet of the coastal lagoon, suggesting a potential association with the physicochemical properties of oceanic seawater (e.g. temperature, salinity [
65,
66]). Locations I and M2 registered the highest levels of Ni in
D. neapolitana jaws. This trend was previously reported for the body of
D. neapolitana by Pires, et al. [
67] exactly on those same locations. High levels of Ni could be associated with anthropogenic impacts [
68,
69], such as the presence of an important shipyard, commercial harbor, and industrial activities in location I [
70], along with agriculture runoff in location M2 [
71].
The use of Random Forest classifications based on FA profiles and EF enhanced the discriminations of geographic origin between specimens of
D. neapolitana (
Table 1). Indeed, the use of each tool individually was associated with some constrains. When FA profiles were used alone to determine the geographic origin of collected specimens, location E was well discriminated, contrarily to the other locations. When using EF, locations I and M1 were discriminated with a high level of accuracy, contrarily to locations E and M2 (
Table 1). Thus, the combination of FA profiles and EF showed to be a more efficient approach to successfully allocate sampled specimens of
D. neapolitana to their true geographic origin, thus increasing the success rate (
Table 1). The combination of different fingerprinting methods (e.g. FA profiles with stable isotopes or multi-elements with stable isotope) had already been successfully employed to discriminate the geographic origin of different species. When using FA profiles combined with stable isotope analysis, the geographic origin of sea cucumbers (
Apostichopus japonicus) [
72], jumbo squids (
Dosidicus gigas) and scallops (
Patinopecten yessoensis,
Chlamys farreri, and
Argopecten irradians) [
33] was determined with high success rate of correct allocations. Employing a multi-elements and stable isotope approach for different shrimp species (
Penaeus indicus, P
. merguiensis,
P. monodon,
P. notialis,
P vannamei,
Pleoticus muelleri, and
Pandalus borealis) was also highly successful when aiming to allocate their geographic origin [
73]. These approaches were also used to distinguish production methods (farmed vs. wild), namely by using FA profiles and stable isotopes on European eel (
Anguilla anguilla) [
74], and Atlantic salmon (
Salmo salar) [
75], as well as for shrimps when using multi-elements and stable isotope analysis (see above; [
73]).
Author Contributions
Conceptualization, F.R. and R.C; methodology, F.R., M.L.L., and R.M.; validation, M.R.D., C.P. and R.C.; formal analysis, F.R. and R.C.; investigation, F.R., M.R.D., C.P. and R.C.; resources, M.R.D., C.P. and R.C.; writing—original draft preparation, F.R. and R.C.; writing—review and editing, F.R., M.L.L., R.M., M.R.D., E.F.S., C.P. and R.C.; visualization, F.R.; supervision, M.R.D., E.F.S., C.P. and R.C.; project administration, R.C.; funding acquisition, C.P. and R.C. All authors have read and agreed to the published version of the manuscript.