Commercial traceability of Arapaima spp. ﬁsheries in the Amazon basin: can biogeochemical tags be useful?

. The development of analytical tools to determine the origin of ﬁshes is useful to better understand patterns of habitat use and to monitor, manage, and control ﬁsheries, including certiﬁcation of food origin. The application of isotopic analyses to study calciﬁed structures of ﬁshes (scales, vertebrae, and otoliths) may provide robust information about the ﬁsh geographic origin and environmental living conditions. In this study, we used Sr and C isotopic markers recorded in otoliths of wild and farmed commer-cialized pirarucu ( Arapaima spp.) to evaluate their prediction potential to trace the ﬁshes origin. Wild and farmed ﬁsh specimens, as well as food used for feeding pirarucu in cap-tivity, were collected from different sites. Isotope composition of otoliths performed by isotope-ratio mass spectrometry (IRMS; δ 13 C) and femtosecond laser ablation multi-collector inductively coupled plasma mass spectrometry (LAfs-MC-ICPMS; 87 Sr / 86 Sr) were compared to the isotopic composition of water and of the food given to the ﬁshes in the farms. Wild ﬁsh specimens that lived in environments with the largest ﬂuctuation of river water Sr isotope ratios over time presented the largest Sr isotope variations in otoliths. A quadratic discriminant analysis on otolith isotopic composition provided 58 % of correct classiﬁcation for ﬁsh production (wild and farmed) and 76 % of correct classiﬁcation for the ﬁsh region. Classiﬁcation accuracy for region varied between 100 % and 29 % for the Madeira and the Lower Amazon ﬁshes, respectively. Overall, this preliminary trial is not yet fully developed to be applied as a commercial traceability tool. However, given the importance of Arapaima spp. for food security and the generation of economic resources for millions of people in the Amazon basin, further analyses are needed to increase the discrimination performance of these biogeographical tags.

Abstract. The development of analytical tools to determine the origin of fishes is useful to better understand patterns of habitat use and to monitor, manage, and control fisheries, including certification of food origin. The application of isotopic analyses to study calcified structures of fishes (scales, vertebrae, and otoliths) may provide robust information about the fish geographic origin and environmental living conditions. In this study, we used Sr and C isotopic markers recorded in otoliths of wild and farmed commercialized pirarucu (Arapaima spp.) to evaluate their prediction potential to trace the fishes origin. Wild and farmed fish specimens, as well as food used for feeding pirarucu in captivity, were collected from different sites. Isotope composition of otoliths performed by isotope-ratio mass spectrometry (IRMS; δ 13 C) and femtosecond laser ablation multi-collector inductively coupled plasma mass spectrometry (LAfs-MC-ICPMS; 87 Sr/ 86 Sr) were compared to the isotopic composition of water and of the food given to the fishes in the farms. Wild fish specimens that lived in environments with the largest fluctuation of river water Sr isotope ratios over time presented the largest Sr isotope variations in otoliths. A quadratic discriminant analysis on otolith isotopic composition provided 58 % of correct classification for fish production (wild and farmed) and 76 % of correct classification for the fish region. Classification accuracy for region varied between 100 % and 29 % for the Madeira and the Lower Amazon fishes, respectively. Overall, this preliminary trial is not yet fully developed to be applied as a commercial traceability tool. However, given the importance of Arapaima spp. for food security and the generation of economic resources for millions of people in the Amazon basin, further analyses are needed to increase the discrimination performance of these biogeographical tags.

Introduction
Food production is becoming increasingly associated with sustainable practices ensuring environmental preservation goals. Food origin and production conditions have become important issues in national and international trade to attest to adequate practices. For instance, consumers want to know whether the fishes belong to an endangered or vulnerable species and whether they were grown in natural or farmed conditions (Baffi and Trincherini, 2016;Kim et al., 2015;Pracheil et al., 2014). In order to address some of these questions, and further improve wild and farming fish management, recent studies have used biogeochemical tracers to better understand fish population dynamics, their ecological strategies, and stock origin (Brennan and Schindler, 2017;Kennedy et al., 2005;Kerr and Campana, 2013;Pracheil et al., 2014;Rojas et al., 2007;Thresher, 1999). These tools have also been used to control national and international fish trade (Pracheil et al., 2014) and to identify the geographical origin and conditions under which fishes have been raised (Barnett-Johnson et al., 2008;Bell et al., 2007;Rojas et al., 2007;Turchini et al., 2009). A sustainable community-based management is important to ensure secure and fair food production, thus valorizing the local economy, respecting the ecosystem functioning, and maintaining ecosystem services.
Determination of fish geographic origin by stable isotope analyses has been investigated by different authors aiming to characterize food origin, manage fisheries and fish stocks, build knowledge of species life history, identify critical habitats for conservation, and avoid overexploitation of fish stocks (Baffi and Trincherini, 2016;Duponchelle et al., 2016;Garcez et al., 2015;Hauser, 2018;Hegg et al., 2015;Jordaan et al., 2016;Pouilly et al., 2014;Pracheil et al., 2014). The methodological assumptions to discriminate fish stocks are based on (i) the isotopic heterogeneity among stocks and (ii) the low mobilization rates in the analyzed tissue (Kerr and Campana, 2013). Isotopic information can be extracted from calcified structures of fishes (scale, otolith, or vertebrae) in order to preserve fish integrity for commercial use and eventually reconstruct fish life history (Campana, 1999;Pouilly et al., 2014;Pracheil et al., 2014). Strontium isotopes have been used as an origin tracer of food products because of their robust response in terms of origin authenticity and fraud detection (Baffi and Trincherini, 2016). On the other hand, carbon isotopes have been used to distinguish farmed and wild fishes according to feeding patterns (Rojas et al., 2007;Turchini et al., 2009). Fish otoliths, or ear bones, are calcified structures that grow continually and record ambient conditions throughout a fish's life (Campana, 1999). Since Sr isotopes in otoliths are not reabsorbed and do not fractionate during biological uptake, the isotopic ratio of this element is a robust geographic marker (Kennedy et al., 2000;Kerr and Campana, 2013;Pouilly et al., 2014). Most studies using Sr isotopes in fish otoliths were performed on marine and freshwater ecosystems of temperate regions (Comyns, 2008;Gillanders, 2002;Kennedy et al., 1997Kennedy et al., , 2000Kennedy et al., , 2002Woodhead et al., 2005). Only a few studies have focused on fish migration and living conditions in tropical river systems (Duponchelle et al., 2016;Garcez et al., 2015;Hauser, 2018;Hegg et al., 2015;Pouilly et al., 2014;Sousa et al., 2016;Walther et al., 2011).
The Amazon basin has the largest rainforest on the planet and constitutes a complex system of rivers, lakes, and wetlands (Oliveira, 1996). It has a variety of different river water types (white, black, or clear) that drains complex and heterogeneous geologic formations (Santos et al., 2015). Because these rivers drain rocks with different origins and ages, their waters present contrasting Sr isotopic compositions, thus providing an adequate scenario for the application of Sr isotopes as a geographic tracer (Duponchelle et al., 2016;Hauser, 2018;Pouilly et al., 2014;Santos et al., 2015). The region is also known to support a large diversity of fish species, many of which play an important economic role in the region, such as the Arapaima spp., known as one of the largest freshwater fish genus (Hrbek et al., 2007;Queiroz, 2000;Stone, 2007). The four described species (A. agassizii, A. mapae, A. leptosoma, and A. gigas) of this genus are endemic to the Amazon basin, where they are popularly called pirarucu or paiche (Arantes et al., 2010;Stewart, 2013a, b). This genus is socially, economically, and ecologically important in the region because it constitutes one of the main food sources for the local community, providing important economic resources on local and regional scales (farming, fishing, trading). Due to overexploitation, Arapaima spp. have been classified as vulnerable by CITES (Convention on International Trade in Endangered Species of Wild Fauna and Flora) and fishing is subject to legal restriction, such as seasonal fishing prohibition, minimum size of capture, and, most importantly, Arapaima's commercialization is restricted to fishes originated from management areas or aquaculture farms (Feio and Mendes, 2017). Paradoxically, Arapaima spp. are considered exotic invasive species in the upper Madeira watershed in Bolivia and Peru, after being introduced in the region in the 1970s (Van Damme et al., 2011;Figueiredo, 2013;Miranda-Chumacero et al., 2012).
Synergic actions initiated in 1989 and involving communitarian lake management, governmental conservation policies, non-governmental organization projects, and aquaculture production of Arapaima have allowed the recovery of overexploited stocks (Figueiredo, 2013;McGrath et al., 2015). For instance, more than 500 groups have permission for fishing Arapaima under an annual quota stipulated by the Brazilian Institute of Environment and Renewable Natural Resources (IBAMA), which have limited resources to monitor, control, and regulate fishery stocks (McGrath et al., 2015). Even though the "sustainable fishery" is certified with tags in order to allow traceability, these tags are easily counterfeited and illegally reused. In this context, the development of an isotopic tag to track back the origin of Arapaima fishery could reinforce the actual system of traceability and combat illegal exploitation of Arapaima stocks. Some aspects of Arapaima spp. biology qualify this group as a model for isotopic certification. This genus is described as sedentary, although individuals migrate locally from lateral lakes to rivers during flooding pulses in order to complete their life cycle (Arantes et al., 2010;Araripe et al., 2013;Castello, 2008). Also, because this genus is air breath dependent, individuals come out of the surface to breath regularly and local fishermen can estimate population sizes and realize the certification of the communitarian management. These characteristics lead to a spatial distinction between stocks of the major Amazonian regions and allow us to use biogeochemical tags as a tracer of origin. Moreover, trophic patterns of Arapaima vary among populations and over the ontogeny from omnivorous (Watson et al., 2013) to carnivorous or piscivorous, predominantly based on C 3 sources (Carvalho et al., 2018;Domingues et al., 2006;Queiroz, 2000). Because of this variation in feeding sources (Carvalho et al., 2018;Castello, 2008), a common strategy used in aqua-culture and strongly recommended by governmental manuals of Arapaima farming is the nutritional training with three rations during different life stages (Ono and Kehdi, 2013). As a reference, δ 13 C values of the food used on farming activity vary between −19.3 ‰ and −14.9 ‰, thus corresponding to a diet based mostly on C 4 macrophytes, corn, and soya beans.
In this study, we hypothesize that farmed and wild fishes of different subbasins would present different δ 13 C and 87 Sr/ 86 Sr values, depending on their food sources (δ 13 C) and geographical origins ( 87 Sr/ 86 Sr). Hence, the main objective is to test whether δ 13 C and 87 Sr/ 86 Sr measured on fish otoliths can be used as biogeochemical tags of the geographic origin and provenance (wild vs. farmed) of Arapaima specimens and, consequently, evaluate if these isotopes can be used as a traceability tool. We analyzed Sr and C isotopic composition of Arapaima's otoliths from farmed and wild fishes proceeding from four different Amazonian regions (Madeira, Solimões, Central Amazon, Lower Amazon, Fig. 1). Sr and C isotope data were analyzed across a transect in Arapaima's otoliths in order to identify habitat change or geographic mobility during the fish lifetime. These data were also compared to the C and Sr isotopic composition of fish food supplied by the farmers.

Study area
The Amazon basin represents a dynamic and heterogeneous ecosystem extending over more than 45 % of the surface area of South America. The Amazon River and its huge network of tributaries drain different geological formations (Gaillardet et al., 1997;Gibbs, 1967;Santos et al., 2015;Stallard, 1980;Stallard and Edmond, 1983) covered by primary forests, chaparral savannas, floodplains, and swamps. These habitats are therefore some of the most biodiverse in the world, particularly in regard to the Amazonian freshwater fish fauna, which are under ecological pressure from degradation by dams building, mining, changing land cover, and global climate change (Anderson et al., 2018;Carnicer et al., 2015;Castello et al., 2013;Castello and Macedo, 2016;Finer and Jenkins, 2012;Forsberg et al., 2017;Latrubesse et al., 2017;Lees et al., 2016;Winemiller et al., 2016).
The Amazon basin is a geomorphological depression located between two old and stable geological regions: the Guiana shield in the north and the Brazilian shield in the south. While the Andes mountain chain limits the western border of the basin, cratonic terrains and the Atlantic Ocean limit its eastern border. Owing to its complex geological history, rivers of the Amazon basin drain rocks with a wide range of Sr isotope compositions (Santos et al., 2015). For example, the Madeira and Negro rivers drain old rock formations, such as Precambrian and Ordovician rocks, that imprint a strong radiogenic Sr isotope signature in their waters (respectively 0.7168 ± 0.0007 and 0.7318 ± 0.0074; Santos et al., 2015). In contrast, the Solimões river drains younger formations as well as carbonate rocks, so that their water is characterized by less radiogenic Sr isotope ratios (0.7091±0.0002; Gaillardet et al., 1997;Santos et al., 2015). Because of this heterogeneity, Sr isotopes in the Amazon river waters may be used as a robust biogeographic marker for aquatic fauna (Duponchelle et al., 2016;Garcez et al., 2015;Hauser, 2018;Hegg et al., 2015;Pouilly et al., 2014;Sousa et al., 2016).
The carbon isotopic composition of an organism depends primarily on the isotopic composition of the primary producer that constitutes the basis of the trophic chain it feeds on. In particular, plants present with two main photosynthetic pathways, C 3 and C 4 , of which the δ 13 C range of values are contrasting (−32 ‰ to −24 ‰ for C 3 plants; −14 ‰ to −9 ‰ for C 4 plants; De Niro and Epstein, 1978). In general, wild Amazon fishes feed dominantly on a trophic chain derived from a C 3 carbon source (Araujo-Lima et al., 1986;Forsberg et al., 1993;Jepsen and Winemiller, 2007;Marshall et al., 2008;Mortillaro et al., 2015;Watson et al., 2013).

Fish otolith sampling
Thirty-eight otolith samples of Arapaima spp. were collected in different sites of the four main regions ( Fig. 1): 22 were obtained from professional fishermen for a commercial purpose, 6 were obtained from farmers, and 10 others were collected at Manaus and Santarém markets ( Table 1). The sagitta otoliths of each specimen were extracted by head dissection after capture. Afterward, they were washed, dried, and kept until laboratory analyses. None of the fishes had precise location of capture, and the informed origin in four main regions were ( Fig. 1): Solimões (Mamirauá reserve at the Solimões river), Madeira (Yata, Beni, and Madre de Dios rivers and Ariquemes farm), Central Amazon (Manacapuru farm, at Solimões river; Itacoatiara, at Amazon River; Novo Airão at Negro river) and Lower Amazon (Santarém, Amazon river). These regions have a complex drainage system that may include different sources of water with different Sr isotopic compositions in each site (Table 1). For example, Mamirauá reserve in Solimões region is located at the confluence of five major rivers, whose water types range from white (lower 87 Sr/ 86 Sr values) to black water (higher 87 Sr/ 86 Sr values). Similarly, Itacoatiara is located at the confluence zone of the Amazon and Madeira rivers, which also have quite distinct Sr isotopic compositions. Therefore, the lack of precise information of collection and the variability in regional water chemistry may not be exactly matched by literature data and may be important to explain the unexpected variance in fish otolith data. Table 1. δ 13 C (mean) and 87 Sr/ 86 Sr (mean) measured in sagittae otolith for 38 wild and farmed Arapaima gigas proceeding from four Amazonian regions with correct classification and correspondent water 87

Samples preparation and analytical methods
The otoliths were sonicated in distilled water, dried, and mounted in Araldite epoxy resin at MARBEC laboratory in Montpellier University (France). Afterward, they were transversally cut with a low-speed saw (Isomed Buehler, Düsseldorf, Germany, 2009) to obtain a dorsoventral slice including the otolith core. The slices were then finely polished until the core could be seen, sonicated in distilled water, and mounted on glass using Crystalbond glue. Sr isotope analyses were performed at the Laboratoire de Chimie Analytique Bio-inorganique et Environnement (LCABIE) from the Institut Pluridisciplinaire de Recherche sur l'Environnement et les Matériaux (IPREM), Université de Pau et des Pays de l'Adour and in the laboratory PSO-IFREMER (Pole Spectrometrie Océan, Brest, France). Interlaboratory crosscalibration was performed to confirm the repeatability and comparability of the analysis (see Hauser, 2018, for details). The isotope ratios were measured using an La-fs-MC-ICPMS following the procedure detailed by (Claverie et al., 2009;Tabouret et al., 2010). Laser ablation conditions were 500 Hz and 20 µJ pulse energy until the depth limit of ablation (< 30 µm); the beam spot size was 10 µm and its velocity was 5 µm s −1 . The Sr isotope ratios were obtained by transects (200 µm width; see Tabouret et al., 2010) along the major otolith axes, perpendicular to otolith growth lines. The laser-ablated material was carried with He gas to a double torch chamber in which the laser aerosol was mixed with a 2 % HNO 3 solution before introduction into the plasma (Barats et al., 2007). These conditions were adjusted to obtain the maximal plasma sensibility and stability. Interferent 87 Rb signal was monitored by 85 Rb, and 87 Sr/ 86 Sr was corrected following the procedure of Barnett-Johnson et al. (2010). Similarly, 83 Kr was measured to control 84 Kr and 86 Kr impact in 84 Sr and 88 Sr values, respectively. Finally, the ratio 88 Sr/ 86 Sr was used to correct 87 Sr/ 88 Sr and mass bias using the exponential law (Walther and Thorrold, 2008). Internal pattern 87 Sr/ 86 Sr ratio (NIESS 22, certificated by the National Institute of Japan Environmental Studies) was analyzed at the beginning and end of each session of analysis to check the repeatability of the 87 Sr/ 86 Sr measurements.
Complete transects from core to edge were performed on 10 wild otolith samples, all of which presented a flat 87 Sr/ 86 Sr ratio pattern along the transect. For the remaining samples, the transect was performed only on the final onethird of the otoliths, which records the environmental condition during the last life period of the adult fish (Fig. 2). The results presented hereafter corresponds to the final part of the otolith for all individuals.
To obtain dietary information based on carbon isotopes, the same last one-third of each otolith was micro-drilled on the same slice preparations of otolith used to gather the Sr isotope data. The slice preparations were drilled with intervals of 6.8 mm using a "New Wave Microdrill" at the Universidade de Brasília. The drilled carbonate powder samples were placed directly into vials for isotope analysis. Carbon isotopes were measured using a Delta Plus V Thermo-Fisher mass spectrometer connected to a Finnigan GasBench II at the Laboratório de Isótopos Estáveis (LAIS), Instituto de Geociências Rede de Estudos Geocronológicos, Geodinâmicos e Ambientais (GEOCHRONOS), Universidade de Brasília, Brazil. The results were validated against reference standards NBS 18 and 19 (respectively δ 13 C = −5.0 ‰ and 1.9 ‰).

Statistical analysis
ANOVA was applied to test 87 Sr/ 86 Sr and the δ 13 C mean difference in otoliths among (1) wild fish proceeding from the four regions, and (2) farmed fishes from the four regions. A t test was also applied to evaluate the mean difference between all wild vs. farmed fishes. To evaluate the use of 87 Sr/ 86 Sr and δ 13 C as a predictive tracer of fish origin (farmed or wild) and subbasin/region of capture (Upper Amazon: Madeira and Solimões, Central Amazon, Lower Amazon), a quadratic discriminant analysis (QDA; Anderson et al., 2010;Li et al., 2016) was carried out using a crossvalidation by Jackknifed (leave one out) predictions procedure. All the statistical analyses were performed in R freeware (http://www.r-project.org/, last access: 10 April 2019).
Except for two specimens, the wild fishes presented a narrow range of 87 Sr/ 86 Sr across the otoliths (Fig. 4) and their average 87 Sr/ 86 Sr values are comparable to 87 Sr/ 86 Sr of the river waters in which they were living (Table 1). The first exception is a fish from Central Amazon-Itacoatiara (I2) that exhibited 87 Sr/ 86 Sr similar to other individuals from the same site but also displayed a peak of 87 Sr/ 86 Sr value up to 0.7259 (Fig. 4a). The second exception is a specimen from the Solimões-Mamirauá area (M2) that presented higher 87 Sr/ 86 Sr values (0.7223 ± 0.0001) in comparison to the river values (0.7090-0.7100) and to other individuals from the same site (0.709-0.7110, Fig. 4c).
The data also revealed differences in the variability in 87 Sr/ 86 Sr values among specimens from the same region (Fig. 3). Individuals from Solimões, Central Amazon (0.7090 to 0.7096), and Lower Amazon (0.7086 to 0.7087) presented   Table A1).
to the Negro river waters. In contrast, 87 Sr/ 86 Sr values of the other specimens of the same region presented isotope ratios in the range of Solimões river waters during all their life (Fig. 5a). One specimen from the Lower Amazon -Santarém area (CS1) presented important fluctuations of 87 Sr/ 86 Sr values across the otolith, corresponding to values in the range of Tapajós river; although the two other specimens collected in the same area (CS2, CS3) presented flat profiles with values intermediate between the Amazon and Tapajós river waters (Fig. 5b). Finally, one of the farmed specimens of the Madeira River (R1) also showed a fluctuating 87 Sr/ 86 Sr profile (Fig. 5c), although the other four fishes showed a flat and completely overlapping profile.

Food and otolith carbon isotopic composition
The δ 13 C values of otoliths were significantly different among wild and farming specimens (ANOVA, F = 124.44, p< 0.01). Most samples of wild fish present δ 13 C consistent with C 3 sources (mean −28.9 ± 1.2 ‰). The exceptions were all samples from Itacoatiara (Central Amazon), which display a mean δ 13 C value of −18.4 ‰ ± 1.8, that fall between the C 3 and C 4 signatures (Fig. 6). In contrast, otoliths of farmed fish presented a wide range of δ 13 C (mean −17.1 ‰±7.7, min = 26.0 ‰, max = 4.8 ‰). The otoliths of the fishes from the Madeira region farms presented a mean δ 13 C value of −8.5 ‰ ± 0.1, indicating a strong contribution of C 4 plants in their feeding source. Farmed fish from the market of Santarém (Lower Amazon) presented a mean δ 13 C value of −14.4 ‰ ± 0.8 in their otoliths, thus revealing a contribution of both C 3 and C 4 plants in their feeding source. Otoliths of farmed fish from the market of Manaus (Central Amazon) presented a lower mean δ 13 C value (−24.7 ‰ ± 0.8), indicating a main C 3 feeding source.

QDA discriminant analysis
The 87 Sr/ 86 Sr and δ 13 C isotopes biplot shows that all otolith samples fall within four main groups (Fig. 6). The carbon isotopic composition combined with the average 87 Sr/ 86 Sr ratio of the fish otolith allows us to partially distinguish the different fish origins (farmed or wild) as well as their geographical region.
Wild fishes presented more negative δ 13 C values (mean −24.3 ‰ ± 3.2) corresponding to a diet more influenced by a C 3 macrophyte carbon source, whereas farmed fishes presented less negative δ 13 C values (mean −17.1 ‰ ± 7.7), corresponding to a higher influence of a C 4 carbon source. How-  Table A1). Figure 6. Biplot of mean δ 13 C and 87 Sr/ 86 Sr for 38 wild (circle symbols) and farmed (triangle symbols) Arapaima otoliths from four geographic regions. Large blue shading represent farmed and wild fish from the Madeira, and yellow squares wild and farmed fishes from Lower Amazon, Central Amazon, and Solimões. The only two samples that are out of their group are tagged with * . ever, farmed fishes presented a higher variability indicating different sources of food depending on the farm. This variability led to low predictability of fish origin (58 % of correct classification, Table 2).
On the contrary, QDA analysis gave a higher score of correct classification of a fish's region (76 %, Table 3). This per-centage varied between subbasins, from 100 % (Madeira) to 29 % (Lower Amazon).

Discussion
The data presented show major differences in otolith's C and Sr isotopic composition among the studied populations of Arapaima. In general, the 87 Sr/ 86 Sr measured in otoliths of wild specimens were similar to the 87 Sr/ 86 Sr reported in the dissolved fraction of the river water in which they were caught. Furthermore, the pattern of variation of the 87 Sr/ 86 Sr ratio throughout the life of these individuals was mostly flat, indicating that studied Arapaima predominantly stayed in water with the same chemical composition (Araripe et al., 2013;Castello, 2004Castello, , 2008Hermann et al., 2016;Núñez-Rodríguez et al., 2015;Queiroz, 2000;Viana et al., 2007). However, some variations were recorded and are related to the lack of precise information of the fishery site and the intricate mosaic of the Amazon water chemistry. They may either correspond to the seasonal variation of water composition, or to movements between the river and adjacent lagoons, which have been shown to have higher 87 Sr/ 86 Sr values than the river (Pouilly et al., 2014). On the other hand, most farmed fishes also presented a flat profile, but some specimens presented abrupt variations in 87 Sr/ 86 Sr values along the otolith. These variations may be produced by changes in the water isotopic composition, due, for example, to a transfer of the fish to another pond, or to a modification of the water source filling the pond. However, the variation could also be a consequence of food change. Most strontium otolith studies indicated the role of ambient water in the control of strontium uptake. Controversially, the role of food in the strontium otolith uptake is debated, as revealed by Sturrock et al. (2012) that reviewed the significance of food up taking processes on the Sr isotopic composition of fish. For example, Kennedy et al. (2000) suggested that food consumption in adult hatchery-reared salmon is preponderant in the uptake of strontium, although Walther and Thorrold (2006) indicated that water chemistry is the dominant factor for marine fishes. Finally, recent advances highlighted that physiological factors may also contribute to Sr control in the otolith, because these elements are transferred into the blood plasma via branchial or intestinal uptake, before reaching the endolymph fluid, and finally the otolith (Payan et al., 2004;Sturrock et al., 2014). Although the importance of food in strontium uptake is not clear, we may speculate that in natural conditions strontium composition of fish and food sources are in equilibrium with river water, so that fish and water are directly correlated, with no significant relative contribution of food in the uptake. In artificial condition, however, strontium composition of exogenous food source could be different from water and may result in a gap between fish and water strontium composition corresponding to the relative importance of food in the uptake.

Isotope record of wild Arapaima
The relationship between Sr isotope ratios in water and fish otoliths or scales has revealed to be a robust tool to study fish migration and geographical origin of population in the Amazon basin (Duponchelle et al., 2016;Garcez et al., 2015;Hauser, 2018;Hegg et al., 2015;Pouilly et al., 2014;Sousa et al., 2016). The 87 Sr/ 86 Sr values of wild fishes from the Lower Amazon and Central Amazon and from the upper Madeira and Solimões rivers presented differences according to the regions (Figs. 3 and 4). Fishes from Lower Amazon and Central Amazon presented the narrowest variation in 87 Sr/ 86 Sr values during their life. These values were around 0.7100, which agrees with the reported values for this river's waters (Santos et al., 2015). One fish presented a peak of the 87 Sr/ 86 Sr value of > 0.7250 (Fig. 4a), which could correspond to a period during which this fish has lived in a habitat with water 87 Sr/ 86 Sr values close to the Madeira River waters or some granitic shield tributaries. Because of the sharpness of the peak, it could also be interpreted as an irregularity in the otolith. In contrast, fishes from the upper Madeira region, including the Beni, Mamoré, Yata, and Madre de Dios rivers, presented a higher degree of variability in 87 Sr/ 86 Sr values (Fig. 4d, from 0.7150 to 0.7350). Fishes from this region also presented a higher 87 Sr/ 86 Sr variation across each otolith profile when compared to fishes from other sites, which is consistent with the natural seasonal variation in these older geological regions. In the Madeira waters and its upper tributaries, Santos et al. (2015) observed dissolved 87 Sr/ 86 Sr data with a high seasonal variation, owing to the nature of rocks being eroded during the rainy and dry sea-sons. The data presented here suggest that fish otoliths also record these seasonal variations.
The wild fishes from the Solimões river were caught in the Mamirauá reserve. This reserve does not correspond to the main channel of the Solimões channel, but to lateral lakes that developed in a mixing zone with other tributaries, some of which may be of black waters that generally present 87 Sr/ 86 Sr similar to that reported for the Negro river (Santos et al., 2015). One of the five fishes analyzed presented higher 87 Sr/ 86 Sr values (> 0.720), suggesting it may have lived part of its life in such a black water tributary. Besides this regional variation, the data presented revealed that in general individual fish lived in waters with a limited range of 87 Sr/ 86 Sr values, suggesting a resident behavior. The pattern of 87 Sr/ 86 Sr along the otolith from the Madeira fishes presented higher variations. This could either result from movements between habitats with contrasted 87 Sr/ 86 Sr water signatures (e.g., adjacent lakes and lagoon, as shown by Pouilly et al., 2014, for the Beni River) or from the integration of the important natural seasonal variations in 87 Sr/ 86 Sr signature in the Madeira waters described by Santos et al. (2015).
Previous studies also concluded a resident behavior of Arapaima species (Araripe et al., 2013;Castello, 2004Castello, , 2008Hermann et al., 2016;Queiroz, 2000;Viana et al., 2007), including a study of individual behavior of restocked and wild Arapaima using radio telemetry (Núñez-Rodríguez et al., 2015). A flat 87 Sr/ 86 Sr profile along the otolith does not directly implicate an absence of movement. Indeed, if a fish moves across two habitats presenting the same isotopic signature, the movement would not be revealed by otolith microchemistry analyses. On the other hand, Castello (2008) demonstrated lateral migration of pirarucu between the Solimões river and the floodplain during water pulse in the Mamirauá reserve with other observation methods. Based on our results, we can conclude that studied Arapaima did not show movements across contrasted habitats (for example white vs. black water systems; Santos et al., 2015). We cannot, however, exclude lateral movements across habitats with similar water signatures. We argue that the 87 Sr/ 86 Sr variations observed in each otolith can be a combination of (1) small changes in the isotopic composition of water due to diverse tributary sources in the hydrological seasonal cycle and/or (2) a lateral migration. As an example, the M3 fish from the Mamirauá showed ripples that might be interpreted as lateral movements between the Solimões white waters and adjacent lagoons or lakes with a slightly higher signature (see Pouilly et al., 2014). Due to the weakness of the pattern and the absence of 87 Sr/ 86 Sr seasonal data from lakes and rivers in the Amazon basin, more detailed studies would be necessary to confirm one or the other hypothesis of movement behavior, which are probably complementary.
Strontium and carbon isotopes in fish otolith record different parameters during a specimen's life. As 87 Sr/ 86 Sr values could be used as a robust fish geographical indicator, even in small scales (Pouilly et al., 2014), carbon isotope compositions (δ 13 C) are related to the feeding source of the fishes.
Most wild fishes analyzed presented δ 13 C values between −24 ‰ and −30 ‰ (Fig. 6). Hence, wild Arapaima in this study had δ 13 C mostly derived from C 3 plants (−28.9 ± 1.2 ‰), as also observed in previous studies (Domingues et al., 2006;Forsberg et al., 1993;Watson et al., 2013). However, a higher contribution of C 4 plants could be observed in some specimens; in particular, fishes from the Central Amazon region (Itacoatiara site) that presented more positive δ 13 C values (−15.4 ‰ and −20.1 ‰). Forsberg et al. (1993) showed that C 3 may account for 82.4 % to 97.5 % of the fish diet in the Amazon basin. Nonetheless, it is worth mentioning that these studies analyzed muscle tissue and there may exist isotopic differences between carbon in muscle tissue and in bone structures because of physiological pathway incorporation. Although this fractionation is not known in Amazonian fishes, available data from Atlantic cod indicate that δ 13 C values can be 15.9 ‰ higher in otoliths than in body tissues (Radtke et al., 1996). If this same fractionation were applied to the otoliths of Arapaima, the carbon source would have an unlikely value of −44.8 ‰. This fractionation difference is likely lower for Amazonian fishes as also suggested by the carbon isotopic composition of other calcified tissues, such as scales reported by Domingues et al. (2006). They showed that these calcified tissues have δ 13 C values between −18.0 ‰ and −29.2 ‰, which are in the same range as samples from the present study. The more positive δ 13 C values observed in fishes from Itacoatiara in the Central Amazon may reflect environmental heterogeneity related to water types (white, black, and clear), channel formations in the dry season, and other hydrologic seasonality related to the flood pulse concept (Domingues et al., 2006;Junk et al., 1989;Oliveira et al., 2006). Moreover, these isotopic values may be related to seasonal resource availability, such as Schizodon fasciatus that presents major digestibility of C 4 macrophytes in the várzea areas (Forsberg et al., 1993;Mortillaro et al., 2015;Oliveira et al., 2006).

Isotopic variations in farmed fish otoliths
In general, farmed fishes also presented a flat 87 Sr/ 86 Sr profile, except for three fishes (CS1, FMn2, and R1) that showed a larger 87 Sr/ 86 Sr profile variation when compared to wild fishes from the same region (Fig. 5). These variations could be related to the physical conditions of the water pond in which they have been raised and/or to abrupt changes in water type, such as seasonal pond transfer. We argue that these fishes were probably used as breeders and that the changes in Sr isotope ratio indicate that they were transferred between different ponds with different Sr isotope compositions. Indeed, fish farming manuals indicate that changing the breeders from one pond to another is an important strategy to increase reproduction (Ono and Kehdi, 2013;SEBRAE, 2010). On the other hand, four fishes (R2, R3, R4, and R5) of the Madeira farm had the exact same Sr isotopic profile, suggesting they have lived the last part of their lives in the same common pond.
Compared to wild fish, farmed fish also showed a higher variation in δ 13 C, thus indicating more diversified food sources. Fishes from the Madeira (Ariquemes farm) presented less negative δ 13 C values (from −8.7 ‰ to −4.8 ‰), which could be related to C 4 -based food (DeNiro and Epstein, 1978;Sant'Ana et al., 2010), probably containing a large proportion of corn. Farmed fishes from the Lower Amazon (Santarém farms) presented intermediate values (from −15.6 ‰ to −13.7 ‰) and those from the central Amazon (Manaus farms) had more negative values (from −26.0 ‰ to −23.2 ‰) more related to C 3 -based food. Therefore, we conclude that our hypothesis of an artificial alimentation based on C 4 plants is not always verified, and that food farming seems to depend from local or regional production or from feeding strategies used by the farm.

Combining 87 Sr/ 86 Sr and δ 13 C signatures
We aimed at verifying if the combination of Sr and C isotopes may be a powerful tool to distinguish between farmed and wild specimens from different Amazonian regions. The quality of information concerning fish origin is an important parameter for sustainable fish commercialization. The lack of a precise fishing site and the heterogeneity of the water system may bring uncertainties to the data but does not compromise the scope of this study and our conclusions. This is also a sensitive question when considering that commercialization of Arapaima is only allowed from farming or management areas. In this sense, the isotope tool applied in this study can be improved to better control commerce in the actual system of traceability tracking back the origin of fish and combat the illegal reuse of tags on illicit fisheries.
The QDA presented in this study gives the proportion of correct prediction of the origin of fishes (production method: farm or wild, geographic regions). The results showed a good but not sufficient enough (> 75 %) percentage of correct classification of the geographic origin (mainly based on 87 Sr/ 86 Sr values). This percentage is downgraded by the overlaps of 87 Sr/ 86 Sr values of some regions (Solimões, Central Amazon, and Lower Amazon). The lack of contrast in 87 Sr/ 86 Sr values between Lower Amazon, Central Amazon, and Solimões regions leads to a higher confusion: four fishes from the Lower Amazon (of a total of seven) and two fishes from Solimões (of a total of five) were misclassified in the Central Amazon region, most of them prevenient from farmed sources. On the contrary, it is upgraded by some clear contrasts existing in different Amazonian subbasins, such as the Madeira, but we can also indicate the Tapajós or Negro rivers that also presented specific values (Santos et al., 2015;review in Hauser, 2018).
On the one hand, results showed low predictability (58 %) of fish origin (farmed or wild). This is mainly due to the va-riety of food sources used to feed the farmed fishes. We hypothesized that farms used food based on a mixture of C 3 and C 4 plants (soya bean, corn) but some farms apparently used food based on C 3 plants, generating confusion with the food of wild fishes. On the other hand, all fishes sold in Manaus marked as farmed fishes presented C 3 -based δ 13 C signatures. This could mean that these supposedly farmed fish actually came from wild provenance, which is illegal and contributes to the over exploration of this natural resource. False information on the fish provenance would also hamper the precision of our approach.
As a preliminary intent, the method gave some interesting results that emphasize the potential of such analyses to obtain a well performing tool. In only a few cases, the 87 Sr/ 86 Sr values recorded in wild fish otoliths were not in agreement with the 87 Sr/ 86 Sr water of the reported origin. For instance, the 87 Sr/ 86 Sr values of wild Arapaima obtained from Santarém market, Lower Amazon, were similar to those observed in the Solimões river (CS3). Hence, it is possible that these wild specimens were caught in the Solimões river (e.g., Mamirauá reserve) and not in the Santarém area as reported by the fish seller. Nonetheless, because of the scarcity of 87 Sr/ 86 Sr water baseline in this area, a Santarém origin cannot be completely ruled out. Some farmed fishes may also have 87 Sr/ 86 Sr values that are not in agreement with the expected values of the reported origin. For example, farmed fishes from the Lower Amazon (Santarém) were probably raised in a pond filled with water from both the Amazon and Tapajós river. Thus, the farming conditions are likely to interfere with the two tracers used in this study.

Conclusions
The expected differences in δ 13 C between farmed and wild fishes (related to artificial vs. natural food sources) could not be confirmed, owing mainly to the C 4 macrophyte contribution to the natural alimentation of wild fishes and to the use of C 3 -based food sources for farmed fishes. False information on the fish provenance in markets may also have contributed to decreasing the precision of the approach and market sampling should be avoided in future studies. Another weakness of our approach is in the 87 Sr/ 86 Sr data overlapping among Amazon subbasins and the lack of a more extensive 87 Sr/ 86 Sr water baseline. Hence, this preliminary result is not yet fully sufficient to be applied as a commercial traceability tool and further analyses are needed to increase the discrimination performance because millions of people rely on Arapaima spp. for subsistence and income. Nonetheless, these initial results encourage a more detailed seasonal 87 Sr/ 86 Sr water sampling in lakes and rivers in all four regions analyzed, and especially in the Madeira and in the Mamirauá reserve, in order to refine the spatial water base and consequently to understand the causes of the otolith profile variation in wild Arapaima spp. Results also suggest fur-ther aspects of the investigation, such as controlled physiological experiment to clarify the sources (water and food) for 87 Sr/ 86 Sr otolith assimilation pathways in farmed conditions and investigating the actual importance of C 4 macrophyte influence to both farmed and wild Arapaima according to season.
Data availability. Isotopic data are available at https://doi.org/10.5281/zenodo.2646436 (Pereira et al., 2019). Data includes (1) mean carbon isotopic signature and mean Sr isotopic ratio 87 Sr / 86 Sr of Arapaima spp. (see Table 1 of the paper); (2) water 87 Sr / 86 Sr (see Table S1 in the Supplement); (3) wild fish Sr profile: data used to draw Fig. 5 (Variation of 87 Sr / 86 Sr values measured by LAfs-MC-ICPMS on wild fish otolith cores -edge transects); and (4) Farmed fish Sr profile: data used to draw Fig. 6 (Variation of 87 Sr / 86 Sr measured by LAfs-MC-ICPMS on farmed fish otolith cores -edge transects). 87 Sr / 86 Sr profile data correspond to the last 30 % of each otolith transects. For some fishes tagged with * , less than the entire column was used in the graphic of Figs. 5 and 6.
Author contributions. In this work, MP and RVS designed the overall project ideas; RVS and LAP developed the sampling plan. For the project execution, LAP collected samples from the Central Amazon and Lower Amazon and FC and MH sampled in Bolivia and Madeira. Then, the sample preparation, δ 13 C otolith analysis, and 87 Sr/ 86 Sr water analysis were made by LAP with the supervision of RVS, while the 87 Sr/ 86 Sr otolith analysis was performed by MH, CP, and SB with the supervision of MP and FDL. MP and LAP performed statistical and data analysis. Finally, LAP prepared the manuscript with the contributions from all co-authors.
Competing interests. The authors declare that they have no conflict of interest.