# Sequential Isotopic Signature Along Gladius Highlights Contrasted Individual Foraging Strategies of Jumbo Squid (Dosidicus gigas)

Anne Lorrain, Juan Argüelles, Ana Alegre, Bertrand Arnaud, Jean-Marie Munaron, Pierre Richard, Yves Cherel

► To cite this version:

Anne Lorrain, Juan Argüelles, Ana Alegre, Bertrand Arnaud, Jean-Marie Munaron, et al.. Sequential Isotopic Signature Along Gladius Highlights Contrasted Individual Foraging Strategies of Jumbo Squid (Dosidicus gigas). PLoS ONE, 2011, 6 (7), pp.e22194. 10.1371/journal.pone.0022194. hal-00613959

HAL Id: hal-00613959

https://hal.science/hal-00613959v1

Submitted on 26 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Sequential Isotopic Signature Along Gladius Highlights Contrasted Individual Foraging Strategies of Jumbo Squid (*Dosidicus gigas*)

Anne Lorrain¹*, Juan Argüelles², Ana Alegre², Arnaud Bertrand²,³, Jean-Marie Munaron¹, Pierre Richard⁴, Yves Cherel⁵

1 IRD, LEMAR UMR CNRS/UBO/IRD 6539, Plouzané, France, 2 IMARPE, Callao, Lima, Perú, 3 IRD, UMR212 EME IFREMER/IRD/UM2, Sète, France, 4 Littoral, Environnement et Sociétés, UMR 6250 CNRS-Université de La Rochelle, La Rochelle, France, 5 Centre d'Etudes Biologiques de Chizé, UPR 1934 du CNRS, Villers-en-Bois, France

## Abstract

**Background:** Cephalopods play a major role in marine ecosystems, but knowledge of their feeding ecology is limited. In particular, intra- and inter-individual variations in their use of resources has not been adequately explored, although there is growing evidence that individual organisms can vary considerably in the way they use their habitats and resources.

**Methodology/Principal Findings:** Using δ¹³C and δ¹⁵N values of serially sampled gladius (an archival tissue), we examined high resolution variations in the trophic niche of five large (>60 cm mantle length) jumbo squids (*Dosidicus gigas*) that were collected off the coast of Peru. We report the first evidence of large inter-individual differences in jumbo squid foraging strategies with no systematic increase of trophic level with size. Overall, gladius δ¹³C values indicated one or several migrations through the squid's lifetime (~8–9 months), during which δ¹⁵N values also fluctuated (range: 1 to 5‰). One individual showed an unexpected terminal 4.6‰ δ¹⁵N decrease (more than one trophic level), thus indicating a shift from higher- to lower-trophic level prey at that time. The data illustrate the high diversity of prey types and foraging histories of this species at the individual level.

**Conclusions/Significance:** The isotopic signature of gladii proved to be a powerful tool to depict high resolution and ontogenic variations in individual foraging strategies of squids, thus complementing traditional information offered by stomach content analysis and stable isotopes on metabolically active tissues. The observed differences in life history strategies highlight the high degree of plasticity of the jumbo squid and its high potential to adapt to environmental changes.

**Citation:** Lorrain A, Argüelles J, Alegre A, Bertrand A, Munaron J-M, et al. (2011) Sequential Isotopic Signature Along Gladius Highlights Contrasted Individual Foraging Strategies of Jumbo Squid (*Dosidicus gigas*). PLoS ONE 6(7): e22194. doi:10.1371/journal.pone.0022194

**Editor:** Tamara Natasha Romanuk, Dalhousie University, Canada

**Received** February 22, 2011; **Accepted** June 17, 2011; **Published** July 14, 2011

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

**Funding:** This work is a contribution to the cooperation agreement between IMARPE and IRD, and of the LMI DISCOH. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

* E-mail: anne.lorrain@ird.fr

## Introduction

Pelagic squids play a major role in the trophic structure of marine ecosystems as they constitute large pelagic biomasses, and are voracious predators and important prey for other organisms [1]. However, knowledge of their life history traits, including feeding strategies and diets is limited by biases resulting from net feeding and subsequent stomach content analyses [2]. Individual variation (intra- and inter-) in their use of resources is not well known, although there is growing evidence that individuals in a population vary considerably in the way they use their habitats and resources [3,4].

A new tool to investigate long-term trophic ecology at both the population and individual levels has been recently developed by combining stable isotope measurements and sequential sampling of hard tissues that remain inert after synthesis [5–8]. Carbon isotopes (δ¹³C) can be used to investigate the animals’ habitats and migration patterns, while δ¹⁵N values depict the associated trophic

variations [5]. Indeed, in contrast to δ¹⁵N values which show a stepwise enrichment from prey to predators [9], carbon isotope values vary little with each trophic level. They therefore mostly reflect spatial variations of the environment and can indicate offshore vs. inshore, benthic vs. pelagic feeding, or even latitudinal variations in foraging habitats [10,11]. In squids, stable isotope profiles along the gladius, their internal chitinous shell, have been recently analyzed and seem to be a promising tool to produce a chronological record of dietary information over their lifetime [12,13]. The gladius grows continuously by accretion of new molecules of chitin and proteins at the proximal part of the gladius (near the head), with no metabolic turnover after synthesis. Consequently, it retains molecules sequentially laid down throughout the lives of cephalopods, and sequential sampling along this structure gives a dietary record of the organism throughout its life.

The jumbo squid (*Dosidicus gigas*, Ommastrephidae) has a strong economic and ecological importance in the eastern Pacific Ocean.

<figure><img src="image_2.png" /></figure>
It is a fast growing species, from 1 mm mantle length (ML) at birth to more than 1 m ML when 1–2 years old [14]. In the last decade, jumbo squids have shown a significant increase in abundance and distribution in the Eastern Pacific [15,16]. Both climatic and anthropogenic factors impact jumbo squids and marine organisms, but a current hypothesis is that individual feeding histories of jumbo squid strongly control the population structure [17]. Here, we measured sequential isotopic values along the gladius of five large individual jumbo squids that were sampled the same day at two nearby locations (within 50 km) off Peru, thus examining temporal variations in habitat and resource use through their life span. Our objective was to look at intra- and inter-individual isotopic variations to investigate the overall isotopic niche of jumbo squids and their degree of specialization. As this species is thought to be highly migratory [14] and to have ontogenic shifts in its diet [13,18], we expected strong variations in its isotopic carbon and nitrogen gladius profiles. Finally, variability in food habits of large jumbo squids was also investigated by the traditional method of stomach content analyses of several individuals sampled along the northern coast of Peru. These stomach analyses represent a snapshot in animal's diet and contrast with life history reconstructions allowed by longitudinal investigations along the gladius, but they nicely complement our isotopic analyses.

## Materials and Methods

### Ethics statement

Animals in this study were cared for in accordance with the guidelines of the ethics committee of the Institute of Research for Development and the Sea Institute of Peru (agreement 006481/01).

### Squid sampling

Gladii and statoliths of large (>60 cm ML) jumbo squids were sampled from 5 adults collected in offshore waters of Northern Peru on November 19, 2008 (Table 1). They were stored in ethanol and then frozen at -20°C. Age was estimated by counting the number of increments per statolith following [19]. Gladii were cleaned in distilled water for 5 min in an ultrasonic bath. They were then measured, dried and cut into 1 cm consecutive sections respecting the V shape of the growth lines following [13], Fig. 1. The sections sampled at the proximal end of the gladius (near the head) correspond to the most recently deposited material produced when the animal is an adult (oldest age), while distal sections of the gladius correspond to older material deposited when the animal was young. Only the proostracum was sampled,

resulting in 50 to 62 sections per individual for gladii ranging between 65 and 80 cm length. Each cm of gladius averages different amounts of time (time-averaging, see 20) because of the 3-dimensional nature of gladius growth and varying growth rates along squids' life [20,21]. Consequently, proximal (youngest) cm-wide sections represent shorter time periods than distal (oldest) sections. Sections (n = 274) were ground into powder, weighed with a micro-balance, and packed in tin containers for isotopic analysis. The main constituents of gladius are chitin and proteins [22,23]. Hence, with chitin being depleted in ¹⁵N compared to proteins and diet, muscle isotope values of the corresponding individuals were also analyzed to control for these differences. Muscles were dried, ground to a fine powder and lipids were removed using cyclohexane [24]. Finally, the gladius and lipid-free muscle of six small (4–5 cm ML) jumbo squid individuals sampled along the coast of Peru in February 2009 at 3 and 9°S were analyzed to track latitudinal variations in baseline stable isotope values. Their gladius was cut into 1 cm consecutive sections following the method used for larger squids.

### Isotopic analysis

Samples were analyzed using an elemental analyzer (Flash EA 1112, Thermo Scientific, Milan, Italy) coupled to an isotope ratio mass spectrometer (Delta V Advantage with a Conflo IV interface, Thermo Scientific, Bremen, Germany). Results are expressed in standard δ notation based on international standards (Vienna Pee Dee Belemnite for δ¹³C and N₂ in air for δ¹⁵N) following the formula: δ¹³C or δ¹⁵N = [(Rsample/Rstandard) - 1] × 10³ (in ‰), where R is ¹³C/¹²C or ¹⁵N/¹⁴N. Reference gas calibration was done using reference materials (USGS-24, IAEA-CH6, IAEA-600 for carbon; IAEA-N1, -N2, -N3, -600 for nitrogen). Analytical precision based on isotope values of the acetanilide (Thermo Scientific) used to estimate C and N content for each sample series was <0.1‰ both for carbon and nitrogen. Values are means ± s.d.

### Stomach content analysis

Stomach were collected from 37 large jumbo squids (range 65–85 cm ML) caught during 2008 along a latitudinal gradient (from 4°S to 11°S). They were kept frozen (–20°C) until analysis. In the laboratory, fresh remains were divided into the main prey items, which were weighed to calculate their proportion by mass in the diet. Prey items were identified to the minimum possible taxon using published keys and descriptions and by comparison with material held in our own reference collection, using fish otoliths and bones, cephalopod beaks and crustacean exoskeletons.

<table><caption>Table 1. Age, size, coordinates, isotope values and C/N mass ratios of the five large jumbo squids.</caption><thead><tr><th rowspan="2">Size (cm)</th><th rowspan="2">Coordinates</th><th rowspan="2">Age (days)</th><th colspan="3">Muscle</th><th colspan="6">Gladius</th><th colspan="6">C/N</th></tr><tr><th>δ¹⁵N</th><th>δ¹³C</th><th>C/N</th><th>N</th><th>Mean</th><th>Max</th><th>Min</th><th>Range</th><th>Mean</th><th>Max</th><th>Min</th><th>Range</th><th>Mean</th><th>Max</th><th>Min</th><th>Range</th></tr></thead><tbody><tr><td>65.2</td><td>8°37'S, 80°42'W</td><td>258</td><td>12.5</td><td>-15.9</td><td>3.2</td><td>50</td><td>8.2</td><td>9.2</td><td>7.5</td><td>0.7</td><td>-15.8</td><td>-14.9</td><td>-16.6</td><td>1.7</td><td>3.9</td><td>4.1</td><td>3.8</td><td>0.3</td></tr><tr><td>72.2</td><td>ND</td><td>14.8</td><td>-16.0</td><td>3.2</td><td>53</td><td>10.0</td><td>12.5</td><td>7.4</td><td>2.7</td><td></td><td>-16.1</td><td>-15.0</td><td>-17.0</td><td>2.0</td><td>3.8</td><td>4.0</td><td>3.8</td><td>0.2</td></tr><tr><td>71.2</td><td>242</td><td>11.5</td><td>-16.5</td><td>3.2</td><td>58</td><td>6.9</td><td>7.5</td><td>6.4</td><td>0.5</td><td></td><td>-16.9</td><td>-16.0</td><td>-17.6</td><td>1.7</td><td>3.8</td><td>4.0</td><td>3.7</td><td>0.3</td></tr><tr><td>81.2</td><td>8°14'S, 80°32'W</td><td>276</td><td>17.7</td><td>-16.2</td><td>3.2</td><td>62</td><td>12.5</td><td>14.6</td><td>9.6</td><td>2.9</td><td>-16.4</td><td>-15.5</td><td>-16.8</td><td>1.3</td><td>3.8</td><td>4.0</td><td>3.8</td><td>0.2</td></tr><tr><td>66.5</td><td></td><td>252</td><td>12.4</td><td>-15.6</td><td>3.2</td><td>51</td><td>7.4</td><td>8.4</td><td>6.5</td><td>0.9</td><td>-15.7</td><td>-14.9</td><td>-16.7</td><td>1.8</td><td>3.9</td><td>4.0</td><td>3.8</td><td>0.2</td></tr></tbody></table>

Bold values correspond to individuals represented on figs 2 and 3. N corresponds to the number of serial samples along the gladius.

doi:10.1371/journal.pone.0022194.t001

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Jumbo squid gladius and isotopic sampling. Schematic dorsal view of *Dosidicus gigas* gladius (adapted from Perez et al. 1996). Sampling was done through 1 cm sections along the gladius proostracum (a) without sampling the conus (b) following the V-shape of the growth lines. Sections close to the proximal end correspond to the most recent formed material and to the oldest age of squid.
doi:10.1371/journal.pone.0022194.g001</figcaption></figure>

## Results

Assuming a deposition rate of one increment per day in statoliths, the lifespan of the five large individuals was ~8–9 months (257±14 days), and ~2 months (57±7 days) for the six small individuals (4–5 cm ML) (Table 1, 2). C/N ratios were consistent among and within individuals (range: 3.8–4.1), thus indicating similar biochemical composition of the samples allowing inter- and intra-individual comparisons (Table 1). Gladii isotopic signatures spread over a relatively small range of δ¹³C values (from −17.6 to −14.9‰, a 2.7‰ difference), but over a much larger range of δ¹⁵N values (from 6.4 to 14.6‰, a 8.2‰ difference). The same patterns were observed with muscle isotopic values, with δ¹³C ranging from −16.5 to −15.6‰ and δ¹⁵N from 11.5 to 17.7‰ (Table 1). Mean muscle and gladius isotope values were highly correlated for both δ¹³C (δ¹³C<sub>muscle</sub> = 0.67×δ¹³C<sub>gladius</sub>−5.22, R² = 0.94, n = 5) and δ¹⁵N (δ¹⁵N<sub>muscle</sub> = 1.09×δ¹⁵N<sub>gladius</sub>+3.96, R² = 0.99, n = 5). All data are provided in Table S1. Overall, the isotopic profiles along the gladius showed no consistency (i.e. they had different trajectories) both between and within the five large individuals, with all isotopic records shifting through time (Fig. S1).

In all squids, intra-individual gladius δ¹³C variations were substantial (from 1.3 to 2.0‰), indicating foraging along isotopic gradients.

Nitrogen isotopes also showed temporal variations through time (from 1.1 to 5.0‰) with no systematic increase through size/age. Two of the five individuals (individuals A and B) adopted contrasting foraging strategies (Figs. 2, 3), i.e. they showed similar δ¹³C values, but their mean δ¹⁵N values (8.2 vs. 12.5‰) and ranges were different (Table 1). This strong isotope difference in gladius δ¹⁵N values among individuals A and B was also observed in their muscle isotopic values (12.5 vs. 17.7‰). Noticeably, individual B showed a strong δ¹⁵N decrease (4.6‰) at the end of its life (just before capture).

Isotopic values of small squid gladii (4–5 cm) sampled at 3°S and 9°S (Table 2) revealed strong latitudinal differences in isotopic baseline levels for both δ¹³C (3.4‰ range) and δ¹⁵N values (5.2‰ range). Indeed, trophic level differences alone cannot account for these large variations in such small animals. As for large squids, mean muscle and gladius isotope values were positively correlated (n = 6, δ¹³C<sub>muscle</sub> = 0.88×δ¹³C<sub>gladius</sub>−1.14, R² = 1.00; δ¹⁵N<sub>muscle</sub> = 0.99×δ¹⁵N<sub>gladius</sub>+4.52, R² = 0.99).

Overall, stomach content analysis across northern Peru revealed that large jumbo squids feed on a large diversity of prey (Table 3). Food was dominated by fish (in particular *Vinciguerria* sp., Nomeidae) and cephalopods (*Dosidicus gigas* and other species) but crustaceans such as euphausiids can also make a large contribution to the diet (up to 45% by mass).

<table><caption>**Table 2.** Age, size, isotopic values (± s.d) and C/N mass ratios of three small individuals per latitude.</caption><thead><tr><th rowspan="3">Size (cm)</th><th rowspan="3">Coordinates</th><th rowspan="3">Age (days)</th><th colspan="4">Muscle</th><th colspan="6">Gladius</th><th colspan="6">C/N</th></tr><tr><th colspan="3">δ¹⁵N</th><th rowspan="2">N</th><th colspan="3">δ¹⁵N</th><th colspan="3">δ¹³C</th><th colspan="3">C/N</th></tr><tr><th>δ¹⁵N</th><th>δ¹³C</th><th>C/N</th><th>Mean</th><th>Max</th><th>Min</th><th>Range</th><th>Mean</th><th>Max</th><th>Min</th><th>Range</th><th>Mean</th><th>Max</th><th>Min</th><th>Range</th></tr></thead><tbody><tr><td>6.4</td><td>3°29'S, 81°59'W</td><td>68</td><td>9.5</td><td>-18.9</td><td>3.3</td><td>4</td><td>5.1</td><td>5.9</td><td>4.3</td><td>1.6</td><td>-19.2</td><td>-19.0</td><td>-19.4</td><td>0.4</td><td>4.1</td><td>4.2</td><td>3.9</td><td>0.3</td></tr><tr><td>5.2</td><td></td><td>66</td><td>9.6</td><td>-19.0</td><td>3.3</td><td>3</td><td>4.9</td><td>5.1</td><td>4.7</td><td>0.4</td><td>-19.3</td><td>-19.2</td><td>-19.4</td><td>0.2</td><td>4.3</td><td>4.3</td><td>3.9</td><td>0.4</td></tr><tr><td>4.1</td><td></td><td>56</td><td>9.1</td><td>-18.8</td><td>3.3</td><td>3</td><td>4.8</td><td>5.2</td><td>4.5</td><td>0.7</td><td>-19.3</td><td>-19.2</td><td>-19.6</td><td>0.6</td><td>4.0</td><td>4.2</td><td>3.9</td><td>0.3</td></tr><tr><td>4.2</td><td>9°12'S, 79°10'W</td><td>56</td><td>14.5</td><td>-16.3</td><td>3.3</td><td>3</td><td>10.5</td><td>11.0</td><td>10.0</td><td>1.0</td><td>-16.3</td><td>-16.2</td><td>-16.4</td><td>0.2</td><td>4.0</td><td>4.0</td><td>4.0</td><td>0.0</td></tr><tr><td>4.0</td><td></td><td>53</td><td>14.6</td><td>-16.3</td><td>3.3</td><td>3</td><td>10.1</td><td>11.0</td><td>9.6</td><td>1.4</td><td>-16.2</td><td>-16.0</td><td>-16.4</td><td>0.4</td><td>4.0</td><td>4.0</td><td>3.9</td><td>0.1</td></tr><tr><td>4.4</td><td></td><td>52</td><td>14.7</td><td>-16.2</td><td>3.2</td><td>3</td><td>9.8</td><td>10.0</td><td>9.3</td><td>0.7</td><td>-16.3</td><td>-16.0</td><td>-16.6</td><td>0.6</td><td>4.1</td><td>4.2</td><td>4.1</td><td>0.1</td></tr></tbody></table>

N corresponds to the number of serial samples along the gladius.

doi:10.1371/journal.pone.0022194.t002

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** High resolution isotopic profile along individual A gladius. δ¹³C (filled circles) and δ¹⁵N (filled diamonds) values along the length of the gladius of a 65.2 cm ML jumbo squid. Grey symbols represent periods of migration, while black symbols illustrate a more fixed isotopic habitat.</figcaption></figure>

doi:10.1371/journal.pone.0022194.g002

Discussion

We report here the first results on jumbo squid trophic isotope ecology off Peru using sequential stable isotope values along the gladius of five large individuals. Previous work [13] reported that the jumbo squid gladius is a powerful tool to determine the geographic origin of squids and their ontogenic variations, but they did not investigate variations within and between individuals at a given location. Our results reveal strong intra- and inter-individual variations in stable isotope values along gladii of jumbo squids, underlining the potential of the isotopic tool to depict individual variations in foraging strategies.

In marine ecosystems, δ¹³C values vary greatly with latitude and/or inshore/offshore gradients [10,25]. Carbon isotope values of small individuals caught at 3 and 9°S showed a 3‰ δ¹³C variation, confirming the presence of latitudinal gradients of the isotopic baseline along the coast of Peru. We cannot however discriminate between this latitudinal and a potential neritic/oceanic gradient in our δ¹³C values. The δ¹³C values of the most recent gladius sections of the five large individuals captured on the same day (B, C and A, D, E) were similar (−15.5, −15.7 and −16.1, −16.4, −16.3‰, respectively), suggesting foraging within the same habitat before capture. Intra-individual δ¹³C variations on gladius isotopic profiles (Figs. 2, 3) were interpreted as reflecting

<figure><img src="image_7.png" /><figcaption>**Figure 3.** High resolution isotopic profile along individual B gladius. δ¹³C (filled circles) and δ¹⁵N (filled diamonds) values along the length of the gladius of a 81.2 cm ML jumbo squid. Grey symbols represent period of migration, while black symbols illustrate a more fixed isotopic habitat.</figcaption></figure>

doi:10.1371/journal.pone.0022194.g003

either migrations or more resident periods. Flat lying (invariant) sections of the carbon isotope profiles are assumed to represent resident locations or spatially limited movements (i.e., habitat 1, 2, 3 on Figs. 2, 3), while sections where the δ¹³C profiles vary (grey symbols) represent a migration between habitats. With resident and migration period defined, intra-individual δ¹⁵N variations during resident periods (i.e., within habitats 1, 2 or 3; Figs. 2, 3) was interpreted as a change in trophic position. Variability in δ¹⁵N profiles during migrations (changing δ¹³C), however, could not be interpreted since baseline δ¹⁵N variations may occur during squid migrations. Indeed, the δ¹⁵N values of the six small squids sampled from 3°S to 9°S off Peru suggested such baseline isotopic changes with a 5.2‰ δ¹⁵N variation along a 6° latitudinal gradient. The eastern South Pacific off Peru encompasses one of the most intense and shallow oxygen minimum zones (OMZ) of the world's oceans [26,27,28] and these OMZs are generally the site of intense denitrification [29] that is known to increase the baseline δ¹⁵N values [30,31]. In such productive regions, denitrification preferentially removes ¹⁵N-depleted NO₃⁻ and leaves residual nitrate ¹⁵N-enriched [32]. Variable anoxic conditions along the Peruvian coast [e.g., strong anoxic conditions occurring south of 5–10°S, 28] can therefore contribute to latitudinal gradients in the δ¹⁵N baseline level that are reflected in organisms at higher trophic levels.

<figure><img src="image_8.png" /></figure>
<table><caption>**Table 3.** Diet (% by mass) of large jumbo squids (65–85 cm ML) collected off Peru (4–11°S) during 2008.</caption><thead><tr><th rowspan="2"></th><th colspan="3">Latitude</th></tr><tr><th>4-5°S</th><th>8-9°S</th><th>10-11°S</th></tr></thead><tbody><tr><td><strong>Prey</strong></td><td></td><td></td><td></td></tr><tr><td><em>Dosidicus gigas</em></td><td><strong>53.8</strong></td><td>6.2</td><td></td></tr><tr><td>Other Cephalopods</td><td></td><td><strong>62.2</strong></td><td>0.1</td></tr><tr><td>Euphausiids</td><td></td><td><strong>45.4</strong></td><td></td></tr><tr><td><em>Pleuroncodes monodon</em></td><td>0.1</td><td></td><td></td></tr><tr><td>Myctophidae</td><td>0.3</td><td></td><td>0.1</td></tr><tr><td><em>Vinciguerria sp.</em></td><td></td><td><strong>30.6</strong></td><td>2.1</td></tr><tr><td>Nomeidae</td><td></td><td></td><td><strong>96.6</strong></td></tr><tr><td>Other Fishes</td><td>0.4</td><td>1.0</td><td>1.1</td></tr><tr><td><strong>Other information</strong></td><td></td><td></td><td></td></tr><tr><td>Number of stomachs</td><td>14</td><td>16</td><td>7</td></tr><tr><td>% full stomachs</td><td>85.7</td><td>62.5</td><td>85.7</td></tr></tbody></table>

doi:10.1371/journal.pone.0022194.t003

The observed intra-individual δ¹³C variations in jumbo squid gladii revealed a complex life history (i.e. several migrations). Noticeably, overlapping δ¹³C values among individuals were observed, but every squid is unique and exhibits a different isotopic pattern through its life. Using biologging, Bazzino et al. [33] showed that 4 large jumbo squids (>70 cm ML) displayed different horizontal movements within 2–4 weeks of life [33]. They proposed that these changes reflect alterations in foraging behaviour, and searches for optimal foraging locations [33]. Their findings of short time scale migrations [33] are consistent with the long-term differences in gladius isotopic variations we observed among individuals. However, the relatively small intra-individual variations in δ¹³C values together with overlapping values among the five individuals suggest that all individual sampled probably foraged in a restricted area at the scale of northern Peru. Furthermore, the observed intra and inter-individual differences never reached the 3‰ δ¹³C difference observed between small jumbo squids sampled at 3°S and 9°S.

Concerning δ¹⁵N variations, the profiles of the two individuals A and B illustrate the variability in foraging strategies (Figs. 2, 3). While the individual A showed two periods of dietary change that can be interpreted by increases in trophic levels, individual B showed a very different and distinct pattern. A 2.2‰ δ¹⁵N increase early in life was followed by a strong 4.6‰ drop, which illustrates a decrease in its trophic level (probably from fish to euphausiids) at the end of its life. Euphausiids have been found to make a large contribution to the diet of large jumbo squids off Chile [>38% in mass, 2]. Stomach contents collected off Peru also showed that they forage for a large variety of prey, including euphausiids (Table 3). Our results therefore reveal the high variability of trophic patterns of large jumbo squids, and refute the systematic increase of trophic level through ontogeny reported by most studies in squids [e.g. 8,12,13,18], except by [34] who did not find ontogenic differences in jumbo squids (14–87 cm) diet in the Gulf of California.

Interestingly, muscle δ¹⁵N values of individuals A and B (12.5 vs. 17.7‰) would have lead us to incorrectly conclude that they had different trophic levels, with squid B probably being cannibalistic. However, gladius analyses showed that δ¹⁵N values were different throughout almost all their life, suggesting that the two individuals probably came from different locations with

different isotopic baselines. Individual B, which presented a much higher δ¹⁵N, most probably came from further south where the OMZ is more intense (and therefore so is denitrification) and shallower [28,35], both which would increase the δ¹⁵N values of the baseline level (as is suggested by the highest δ¹⁵N values of the 3 small jumbo squids sampled at 9°S, Table 2). Furthermore, the strong decrease of trophic level at the end of the life of individual B seems to preclude a cannibalistic comportment and even to suggest a lower trophic level than individual A.

Finally, simultaneous analysis of gladius and muscle in both small and large squids showed a strong correlation between average muscle and gladius isotopic values, with similar values for δ¹³C (0.2±0.2‰ difference) but consistently lower δ¹⁵N values for gladius (−4.5±0.3‰) compared to muscles. It must be kept in mind however that gladius and muscles have different metabolic turnover rates. Muscle values correspond to an integration of a few weeks of life whereas average values provided for gladius correspond to the squid's entire life. The correlations and average differences between isotopic values of muscles and gladius are therefore only valid for the small individuals (2 months old) where muscle and gladius match the same time period. These consistently lower δ¹⁵N values in gladius are due to the presence of chitin, being depleted in ¹⁵N (but not in ¹³C) compared to proteins which are the main constituents of muscles [22,23]. The knowledge of this 4.5‰ difference between gladii and muscles can be useful for comparison purposes in future studies where only one tissue is available, for example to calculate trophic positions through gladius samples.

In summary, these results show that jumbo squids living in the same environment at a given time can have completely different historical backgrounds. The data suggest that (i) jumbo squid actively migrate, (ii) they do not have a fixed diet and foraging strategy, and (iii) they do not systematically increase their trophic level with age. These differences in life history strategies confirm the high degree of plasticity of the species and its high potential to rapidly colonize new areas and adapt to environmental variability. The next step will be to use stable isotopes in combination with tagging studies to better assess the time-changes in habitats and their associated resources that characterize the life history of jumbo squids. Stable isotopes in squid gladii therefore offer an excellent opportunity to depict high resolution and ontogenic variations in individual foraging strategies of squids, resolving some bias inherent to stomach content analysis and stable isotopes on metabolically active tissues.

## Supporting Information

**Figure S1** δ¹³C (upper panel) and δ¹⁵N (lower panel) values along the length of the gladius of the five large jumbo squids. Grey and black symbols represent sampling at different places.

(TIF)

**Table S1** Gladius isotopic values and C/N mass ratios of individuals A to E.

(DOC)

## Acknowledgments

The authors thank A. Arkhipkin for statolith analysis, G. Guillou for stable isotope analysis, N. Coic for figure design and D.P. Gillikin and two anonymous reviewers for their useful comments on the manuscript.

## Author Contributions

Conceived and designed the experiments: AL YC JA. Performed the experiments: AL JA AA JMM. Analyzed the data: AL AB YC. Contributed reagents/materials/analysis tools: PR JMM. Wrote the paper: AL AB YC.

<figure><img src="image_9.png" /></figure>
References

1. Rodhouse PG, Nigmatullin CM (1996) Role as consumers. Phil Trans R Soc B 351: 1003-1022. doi:10.1098/rstb.1996.0090.

2. Ibañez CM, Arancibia H, Cubillos LA (2008) Biases in determining the diet of jumbo squid *Dosidicus gigas* (D'Orbigny 1835) (Cephalopoda: Ommastrephidae) off southern-central Chile (34°S-40°S). Helgol Mar Res 62: 331-338. doi: 10.1007/s10152-008-0120-0.

3. Bolnick DI, Svanback R, Fordyce JA, Yang LH, Davis JM, et al. (2003) The ecology of individuals: incidence and implications of individual specialization. Am Nat 161: 1-28. doi:10.1086/343878.

4. Bearhop S, Adams CE, Waldron S, Fuller RA, Macleod H (2004) Determining trophic niche width: a novel approach using stable isotope analysis. J Anim Ecol 73: 1007-1012. doi:10.1111/j.0021-8790.2004.00861.x.

5. Cherel Y, Kernaléguen L, Richard P, Guinet G (2009) Whisker isotopic signature depicts migration patterns and multi-year intra- and inter-individual foraging strategies in fur seals. Biol Lett 5: 830-832. doi:10.1098/rsbl.2009.0552.

6. Newsome SD, Tinker MT, Monson DH, Oftedal OT, Ralls K, et al. (2009). Using stable isotopes to investigate individual diet specialization in California sea otters (*Enhydra lutris nereis*). Ecology 90: 961-974. doi:10.1890/07-1812.1.

7. Vander Zanden HB, Bjorndal KA, Reich KJ, Bolten AB (2010) Individual specialists in a generalist population: results from a long-term stable isotope series. Biol Lett 6: 711-714. doi:10.1098/rsbl.2010.0124.

8. Guerra A, Rodriguez-Navarro AB, Gonzalez AF, Romanek CS, Alvarez-Lloret P, et al. (2010) Life-history traits of the giant squid *Architeuthis dux* revealed from stable isotope signatures recorded in beaks. ICES J Mar Sci 67: 1425-1431. doi: 10.1093/icesjms/fsq091.

9. Vanderklift MA, Ponsard S (2003) Source of variation in consumer- diet δ¹⁵N enrichment: a meta-analysis. Oecologia 136: 169-182. doi: 10.1007/s00442-003-1270-z.

10. Cherel Y, Hobson KA (2007) Geographical variation in carbon stable isotope signatures of marine predators: a tool to investigate their foraging areas in the Southern Ocean. Mar Ecol Prog Ser 329: 281-287. doi:10.3354/meps329281.

11. Jaeger A, Lecomte VJ, Weimerskirch H, Richard P, Cherel Y (2010) Seabird satellite tracking validates the use of latitudinal isoscapes to depict predators' foraging areas in the Southern Ocean. Rapid Commun Mass Spectrom 24: 3456-3460. doi: 10.1002/rcm.4792.

12. Cherel Y, Fontaine C, Jackson GD, Jackson CH, Richard P (2009) Tissue, ontogenic and sex-related differences in δ¹³C and δ¹⁵N values of the oceanic squid *Todarodes filippovae* (Cephalopoda: Ommastrephidae). Mar Biol 156: 699-708. doi:10.1007/s00227-008-1121-x.

13. Ruiz-Cooley RI, Villa EC, Gould WR (2010) Ontogenetic variation of δ¹³C and δ¹⁵N recorded in the gladius of the jumbo squid *Dosidicus gigas*: geographic differences. Mar Ecol Prog Ser 399: 187-198. doi: 10.3354/meps08383.

14. Nigmatullin CHM, Nesis KN, Arkhipkin AI (2001) A review of the biology of the jumbo squid *Dosidicus gigas* (Cephalopoda: Ommastrephidae). Fish Res 54: 9-19.

15. Bograd SJ, Castro CG, Di Lorenzo E, Palacios DM, Bailey G, et al. (2008) Oxygen declines and the shoaling of the hypoxic boundary in the California Current. Geophys Res Lett 35: L12607. doi:10.1029/2008GL034185.

16. Zeidberg LD, Robison BH (2008) Invasive range expansion by the Humboldt squid, *Dosidicus gigas*, in the eastern. North Pacific Proc Natl Acad Sci 104: 12948-12950.

17. Argüelles J, Tafur R, Taipe A, Villegas P, Keyl F, et al. (2008) Size increment of jumbo flying squid *Dosidicus gigas* mature females in Peruvian waters, 1989-2004. Prog Oceanogr 79: 308-312.

18. Ruiz-Cooley RI, Markaida U, Gendron D, Aguiñiga S (2006) Stable isotopes in jumbo squid (*Dosidicus gigas*) beaks to estimate its trophic position: comparison

between stomach content and stable isotopes. J Mar Biol Ass of the UK 86: 437-445. doi:10.1017/S0025315406013324.

19. Arkhipkin AI (1991) Methods for cephalopod age and growth studies with emphasis on statolith ageing techniques. In: Jereb P, Ragonese S, Boletzky SV, eds. Proceedings of the International Workshop on Squid Age Determination Using Statoliths. NTR-ITPP Publ. Especial No. 1. Istituto di Tecnologia della Pesca e del Pescato, Mazara del Vallo, Sicily, Italy. pp 11-17.

20. Goodwin DH, Schöne BR, Dettman DL (2003) Resolution and delity of oxygen isotopes as paleotemperature proxies in bivalve mollusk shells: Models and observations. Palaios 18: 110-125.

21. Perez JAA, O'Dor RK, Beck P, Dawe EG (1996) Evaluation of gladius dorsal surface structure for age and growth studies of the short-finned squid, *Illex illecebrosus* (Teuthoidea: Ommastrephidae). Can J Fish Aquat Sci 53: 2837-2846.

22. DeNiro MJ, Epstein S (1978) Influence of diet on the distribution of carbon isotopes in animals. Geochim Cosmochim Acta 42: 495-506.

23. Webb S, Hedges REM, Simpson SJ (1998) Diet quality influences the δ¹³C and δ¹⁵N of locusts and their biochemical components. J Exp Biol 201: 2903-2911.

24. Kojidanovik J, Richard P, Le Corre M, Cosson RP, Bustamante P (2008) Effects of lipid extraction on δ¹³C and δ¹⁵N values in seabird muscle, liver and feathers. Waterbirds 31: 169-178.

25. Rau GH, Sweeney RE, Kaplan IR (1982) Plankton ¹³C:¹²C ratio changes with latitude: differences between northern and southern oceans. Deep-Sea Res I 29: 1035-1039.

26. Chavez F, Bertrand A, Guevara-Carrasco R, Soler P, Csirke J (2008) The northern Humboldt Current System: Brief history, present status and a view towards the future. Prog Oceanogr 79: 95-105.

27. Paulmier A, Ruiz-Pino D (2009) Oxygen minimum zones (OMZs) in the modern ocean. Prog Oceanogr 80: 113-128.

28. Bertrand A, Ballón M, Chaigneau A (2010) Acoustic observation of living organisms reveals the oxygen minimum zone. PLoS ONE 5: e10330. doi:10.1371/journal.pone.0010330.

29. Naqvi SWA, Jayakumar DA, Narvekar PV, Naik H, Sarma VVSS, et al. (2000) Increased marine production of N2O due to intensifying anoxia on the Indian continental shelf. Nature 408: 346-349.

30. Sigman DM, Altabet MA, McCorkle DC, Francois R, Fischer G (1999) The δ¹³N of nitrate in the Southern Ocean: nitrate consumption in surface waters. Global Biogeochem Cycles 13: 1149-1166.

31. Graham BS, Koch PL, Newsome SD, McMahon KW, Auroioles D (2010) Using isoscapes to trace movements and foraging behaviour of top predators in oceanic ecosystems. In: West JB, Bowen GJ, Dawson TE, Tu KP, eds. Isoscapes: Understanding Movement, Pattern, and Process on Earth through Isotope Mapping. Berlin: Springer-Verlag. pp 299-318.

32. Voss M, Dippner JW, Montoya JP (2001) Nitrogen isotope patterns in the oxygen-deficient waters of the Eastern Tropical North Pacific Ocean. Deep-Sea Res 48: 1905-1921.

33. Bazzino G, Gilly WF, Markaida U, Salinas-Zavala CA, Ramos-Castillejos R (2010) Horizontal movements, vertical-habitat utilization and diet of the jumbo squid (*Dosidicus gigas*) in the Pacific Ocean off Baja California Sur, Mexico. Prog Oceanogr 86: 59-71. doi:10.1016/j.pocean.2010.04.017.

34. Markaida U, Sosa-Nishizaki O (2003) Food and feeding habits of jumbo squid *Dosidicus gigas* (Cephalopoda: Ommastrephidae) from the Gulf of California, Mexico. J Mar Biol Assoc UK 83: 1-16.

35. Fuenzalida R, Schneider W, Garcés-Vargas J, Bravo L, Lange C (2009) Vertical and horizontal extension of the oxygen minimum zone in the eastern South Pacific Ocean. Deep Sea Res 56: 992-1003.

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