Ontogenetic Trophic Shifts by Ommastrephes bartramii in the North Paciﬁc Ocean Based on Eye Lens Stable Isotopes

: The neon ﬂying squid ( Ommastrephes bartramii ) plays an important ecological role in the North Paciﬁc. An analysis of stable isotopes in eye lenses was conducted to investigate the inter- and intravariation of the trophic ecology of stocks at the eastern and western North Paciﬁc throughout the life cycle. δ 13 C and δ 15 N values gradually increased with ontogenetic growth of the squid, which was associated with geographic migrations and increased the trophic level. For both stocks, from the paralarval to the juvenile stage, the trophic niche breadth increased, which might be the reason that the swimming and feeding ability improved as they entered the juvenile stage. Meanwhile, interactions between different ecosystems led to a greater diversity of food sources; thus, their feeding targets were no longer limited to plankton but shifted toward small ﬁsh and other cephalopods. Then, from the juvenile to the subadult–adult stage, the trophic niche breadth decreased, which can be explained by that O. bartramii had a selective preference for certain prey as ontogenetic growth proceeded, and they seemed to focus more on larger prey in the subadult–adult stage. Furthermore, the small amount of overlap between early and later life cycles suggested a signiﬁcant trophic niche separation among different trophic ecologies and spatial ecologies. This study provides an understanding of diet shifts in neon ﬂying squid in the North Paciﬁc Ocean, primarily including diet shifts during their individual development and differences in trophic variation between the two stocks.


Introduction
The neon flying squid (Ommastrephes bartramii) is an oceanic squid of large size that is widespread in the temperate and subtropical waters of the Pacific, Indian and Atlantic Oceans [1,2]. The species is usually found at depths ranging from the sea surface to about 1500 m and takes part in diel vertical migration, and its vertical and geographical migration varies with its growth [3][4][5]. O. bartramii is an economically and ecologically important species; it is also an opportunist and is short-lived [6,7]. It plays an important role in the structure and functioning of oceanic ecosystems [8] and is a voracious predator that feeds on fish, cephalopods and crustaceans such as the California lanternfish (Symbolophorus californiensis) and the Berry armhook squid (Gonatus berryi) [9]. In turn, O. bartramii is an important prey species for both large fish such as swordfish (Xiphias gladius) [10] and marine mammals such as sperm whales [11]. Currently, O. bartramii is fished commercially on a large scale in the North Pacific Ocean [12]. The fishing grounds are generally considered to be divided into the eastern and western North Pacific [13], with the 170 • E line of latitude acting as the dividing line [14]. Arkhipkin (2005) [45] noted that statoliths are calcareous structures, and they are confirmed to be more suitable for determining individual age based on the growth bands in the statoliths compared with beak, gladius and radula [46]. The growth of these bands is mainly controlled by the organism's internal circadian rhythm [47] and effectively reflects its early growth and development [48]. Growth bands can also be used to determine the time of occurrence of particular time nodes such as the hatching and initial feeding periods [49]. Among others, Arkhipkin and Perez [43] (1998) and Arkhipkin (2005) [45] also indicated that the postnuclear zone, dark zone and peripheral zone in squid statoliths correspond to the paralarval, juvenile and subadult-adult stages, respectively. Furthermore, research on eye lenses can be combined with age determination using statoliths.
Currently, there are extensive studies related to the trophic ecology and habitat of O. bartramii in the North Pacific Ocean. For example, Parry (2008) [50] used muscle nitrogen stable isotopes to study the trophic variation in the range of mantle lengths (MLs) during the ontogenetic stage of O. bartramii individuals; Wang et al. (2022) [51] selected the beaks for carbon and nitrogen stable isotope analysis to study the interannual and ontogenetic variation in the trophic ecology of O. bartramii. However, fewer studies have used eye lens stable isotopes combined with age by using statolith for retrospective analysis of the entire life cycle. The aim of this study is to investigate the ontogenetic changes in terms of trophic ecology and habitat of the squid O. bartramii and compare the difference in the feeding ecology of two squid stocks throughout the life cycle. The findings will help researchers better understand the life history of the neon flying squid in the North Pacific Ocean.  Table 1). A total of 40 O. bartramii samples (23 females and 17 males) were selected, taking into account ML and sex. All of the samples were frozen whole at −20 • C onboard the vessels and immediately transported to the Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai, China, after landing. Measurements were then conducted to acquire biological data, including dorsal mantle length to the nearest 1 mm, body weight (BW) to the nearest 1 g, sex and maturity [52]. The samples were dissected in the lab after defrosting. A total of 40 pairs of statoliths were extracted from the statocysts in the squid heads and stored in 75% ethanol for further analysis. Scalpels and forceps were used to remove the whole of the eye lens from each individual and separate the two parts of the lens. Only the larger, posterior part was kept, which was then stored in 75% ethanol for isotope analysis; the anterior lens was discarded.

Statolith Processing and Increment Counting
In this study, the statoliths were treated by using standard techniques and aging methodology [53]. Uozumi and Ohara (1993) [54] basically supported the "one incrementone day" hypothesis by the continuous sampling method. Following this, their growth increments were then counted to determine age. Longitudinal sections of the statoliths were used so that the growth increments were more distinct and easier to distinguish [55]. The polished statolith slices were placed under an optical microscope (×400 magnification) and photographed using a CCD (Connecting Device, Shanghai, China). Photoshop 8.0 software was then used to superimpose the images. The number of growth increments from the nucleus to postnuclear zones, dark zone and peripheral zone were used to determine the age of paralarval, juvenile and subadult-adult squid, respectively [45]. Two independent counts were made for each statolith and the mean count was accepted as the final value if there was a difference of less than 5% between these counts. Finally, the date of hatching of each squid was back-calculated from its capture date [56].

Lens Delamination and Processing
Deionized water was used to clean posterior lenses to remove any residual tissue, and then a vernier caliper (that measured to the nearest 0.01 mm) was used to measure the lens diameters (LDs). The lenses were then placed into glass Petri dishes using two forceps to delaminate them, that is, separate them into successive segments starting from the outer edge. Before each peeling, vernier calipers were used to remeasure the diameter and width of the last peeling. To make sure that each segment that was to be used in the isotope analysis contained sufficient tissue, the width of each segment was kept to about 0.5 mm. As the squid eyes varied in size, the number of extracted segments varied from three to nine. The extracted segments were placed into 1.5 mL pointed centrifuge tubes and numbered, ready for further analysis.
All of the excised segments were rinsed with Milli-Q ultrapure water for 5 min to remove possible contaminants; they were then freeze-dried using German Christ Alpha 1-4 LSG freeze-dryers at −50 • C for at least 24 h and ground to a fine homogeneous powder using a RETSCH MM200 automatic ball mill. Powder weighing between 1 mg and 2 mg was subsampled and placed into 0.3-mg tin capsules for testing.

Stable Isotope Analysis
The stable isotope values were determined at the Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources of the Ministry of Education, Shanghai 201306, China, by using an ISOPRIME 100 isotope ratio mass spectrometer (Isoprime Corporation, Cheadle, UK). Vienna Pee Dee Belemnite (V-PDB) was used as the standard substance of δ 13 C; atmospheric nitrogen (N 2 ) was used for δ 15 N. In this study, the lipid content was not corrected because of the low C:N ratios (~3.5) [57]. During the sample testing, in order to ensure the stability of the instrument and the accuracy of the results, three standard isotope samples were inserted every tenth sample. The analytical precision of both the δ 13 C and δ 15 N measurements was ±0.2‰. The samples were combusted in the element analyzer to generate CO 2 and N 2 , and the mass spectrometer detected the ratio of 13 C to 12 C in the Fishes 2022, 7, 295 5 of 18 CO 2 and of 15 N to 14 N in the N 2 . By comparing these ratios with those of the standard substances, the δ 13 C and δ 15 N values were then calculated using the formula where X is either 13 C or 15 N, and δ is the sample's heavy isotope to the light isotope ratio. The atomic ratios of 13 C/ 12 C or 15 N/ 14 N in the sample and standard substance are represented by R sample and R standard , respectively.

Data Analysis
Linear regression analysis was used to estimate the ratio of δ 15 N to δ 13 C for the series of eye segments to determine whether this differed from 3, which was the value usually considered to indicate a typical trophic fractionation. The relationship between age and LD was used to back-calculate the lens sizes for individuals at different ontogenetic stages. The SIBER 2.1.6 package [58] in R [59] was used to draw the trophic niches corresponding to different ontogenetic stages, to calculate the niche areas and to estimate niche breadths and calculate the amount of trophic niche overlap. The standard ellipse area corrected (SEAc) contained about 40% of the carbon and nitrogen stable isotope data; thus, it could be used to represent the core niche of each individual [60,61] and estimate the breadth of trophic niches and the amount of trophic niche overlap between different ontogenetic stages. A scale proposed by Langton (1982) [62] was used, and its range of overlap values was from 0 to 1: values of 0-0.29 indicated low overlap, values of 0.30-0.60 indicated medium overlap and values higher than 0.60 corresponded to a high degree of overlap between groups. Significance analysis of differences in δ 13 C and δ 15 N values between different groups and within the groups at different growth stages was conducted by using an analysis of variance (ANOVA). Normal test and homogeneity of variance test were performed before data analysis, and if the above tests were not satisfying, nonparametric tests were performed to analyze the corresponding levels of significant differences, after which multiple comparisons were performed for different growth stages of the same stock and different stocks of the same growth stage, and Bonferroni corrections were performed. In addition, the relationship between LD and ML was then constructed by using simple linear regression, and the Locally Weighted Regression Scatterplot Smoother (LOESS) in R [59] was used to fit the δ 13 C and δ 15 N values of the two squid groups to LD by using a smoothing parameter of 0.5, and correlation analysis was performed between δ 13 C and δ 15 N values and LD, respectively.

Ontogenetic Growth Distribution of the Squid
The mantle lengths of the sampled squid ranged from 176 mm to 534 mm. The dominant ranges were 501-550 mm for the eastern stock and 251-350 mm for the western stock ( Figure 2). Overall, the eye lens diameters ranged from 4.26 mm to 16.43 mm, and the dominant ranges were 12-16 mm and 10-12 mm for the eastern and western stocks, respectively ( Figure 3).

Ontogenetic Growth Distribution of the Squid
The mantle lengths of the sampled squid ranged from 176 mm to 534 mm. The dominant ranges were 501-550 mm for the eastern stock and 251-350 mm for the western stock ( Figure 2). Overall, the eye lens diameters ranged from 4.26 mm to 16.43 mm, and the dominant ranges were 12-16 mm and 10-12 mm for the eastern and western stocks, respectively ( Figure 3).  A total of 29 out of 40 statoliths were processed successfully and their growth increments were counted. The ages of the sampled squids ranged from 130 d to 276 d. For females, the age ranges were 167-276 d for the eastern stock and 142-259 d for the western stock. For males, the age ranges were 130-174 d for the eastern stock and 140-202 d for the western stock. The age corresponding to each ontogenetic stage was determined from the number of growth increments in the statolith. Based on these results, the embryos were found to be aged 0 days, the paralarval 15-41 days, the juvenile 64-110 days and the subadult-adult 116-276 days ( Table 2). Based on the age results of the individuals and collection dates, we back-calculated the hatching date of individuals. The relationship between age and the lens diameter was as follows (

Ontogenetic Growth Distribution of the Squid
The mantle lengths of the sampled squid ranged from 176 mm to 534 mm. The dominant ranges were 501-550 mm for the eastern stock and 251-350 mm for the western stock ( Figure 2). Overall, the eye lens diameters ranged from 4.26 mm to 16.43 mm, and the dominant ranges were 12-16 mm and 10-12 mm for the eastern and western stocks, respectively ( Figure 3).  A total of 29 out of 40 statoliths were processed successfully and their growth increments were counted. The ages of the sampled squids ranged from 130 d to 276 d. For females, the age ranges were 167-276 d for the eastern stock and 142-259 d for the western stock. For males, the age ranges were 130-174 d for the eastern stock and 140-202 d for the western stock. The age corresponding to each ontogenetic stage was determined from the number of growth increments in the statolith. Based on these results, the embryos were found to be aged 0 days, the paralarval 15-41 days, the juvenile 64-110 days and the subadult-adult 116-276 days ( Table 2). Based on the age results of the individuals and collection dates, we back-calculated the hatching date of individuals. The relationship between age and the lens diameter was as follows (  A total of 29 out of 40 statoliths were processed successfully and their growth increments were counted. The ages of the sampled squids ranged from 130 d to 276 d. For females, the age ranges were 167-276 d for the eastern stock and 142-259 d for the western stock. For males, the age ranges were 130-174 d for the eastern stock and 140-202 d for the western stock. The age corresponding to each ontogenetic stage was determined from the number of growth increments in the statolith. Based on these results, the embryos were found to be aged 0 days, the paralarval 15-41 days, the juvenile 64-110 days and the subadult-adult 116-276 days ( Table 2). Based on the age results of the individuals and collection dates, we back-calculated the hatching date of individuals. The relationship between age and the lens diameter was as follows ( Figure 4): LD = 0.0504 × (Age) + 1.3966 (R 2 = 0.5214, p < 0.001, n = 29). The linear relationship between the eye lens diameter and mantle length was significant (p < 0.05, R 2 = 0.782), and the mantle length increased continuously as the eye lens diameter increased ( Figure 5), which meant that the LD, instead of the ML, could be used in the study of δ 13 C and δ 15 N values variation throughout the growth of the squid.

Distribution of Stable Isotope Values for the Eye Lenses
A total of 297 segments were extracted from 40 eye lenses. δ C values ranged from −21.04‰ to −18.89‰, and the δ N values from 3.67‰ to 12.27‰. Among them, there was a significant correlation between δ N and δ C values of successive segments of 19 individuals with positive correlation (p < 0.05). Only six specimens had a value of 3 for the ratio of δ N to δ C for the series of segments (Table 3).

Distribution of Stable Isotope Values for the Eye Lenses
A total of 297 segments were extracted from 40 eye lenses. δ 13 C values ranged from −21.04‰ to −18.89‰, and the δ 15 N values from 3.67‰ to 12.27‰. Among them, there was a significant correlation between δ 15 N and δ 13 C values of successive segments of 19 individuals with positive correlation (p < 0.05). Only six specimens had a value of 3 for the ratio of δ 15 N to δ 13 C for the series of segments (Table 3).
With the comparison of the differences between the δ C and δ N values for the eastern and western squid stocks, it can be seen that there were no significant differences in the δ C and δ N values at the stage of paralarval and juvenile (p > 0.05), and at the stage of subadult-adult, the δ N values had significant difference between two stocks (F = 8.181, p < 0.05), while the opposite was true for δ C values ( Figure 6). In addition, ANOVA showed no significant difference in ML between the two stocks (p > 0.05).   For the western squid stock, the δ 13 C values in the lens segments were found to range from −20.63‰ to −19.36‰ (−19.90‰ ± 0.31‰), −20.55‰ to −19.14‰ (−19.88‰ ± 0.35‰) and −20.79‰ to −18.89‰ (−19.55‰ ± 0.26‰) for the paralarval, juvenile and subadultadult stages, respectively. The corresponding ranges of δ 15 N values were 3.67‰ to 6.34‰ (4.56‰ ± 0.64‰), 3.76‰ to 9.85‰ (6.33‰ ± 1.39‰) and 5.33‰ to 10.06‰ (8.80‰ ± 0.82‰), respectively. The results of ANOVA showed significant differences in δ 13 C and δ 15 N values of different growth stages (F C = 43.769, p C < 0.05; F N = 102.968, p N < 0.05) ( Figure 6; Table 4). With the comparison of the differences between the δ 13 C and δ 15 N values for the eastern and western squid stocks, it can be seen that there were no significant differences in the δ 13 C and δ 15 N values at the stage of paralarval and juvenile (p > 0.05), and at the stage of subadult-adult, the δ 15 N values had significant difference between two stocks (F = 8.181, p < 0.05), while the opposite was true for δ 13 C values ( Figure 6). In addition, ANOVA showed no significant difference in ML between the two stocks (p > 0.05).
For the eastern stock, the SEAc values for the paralarval, juvenile and subadult-adult squid were 0.65‰ 2 , 1.06‰ 2 and 0.43‰ 2 , respectively. There was a trophic separation between the paralarval and subadult-adult squid. There was a low amount of trophic overlap between the paralarval and juvenile, juveniles and subadult-adults, the corresponding values were 0.24 and 0.14, respectively, whereas there was a trophic separation between the paralarval and subadult-adult squid (Figure 7; Table 5). trophic separation between the paralarval and the subadult-adult squid. The trophic overlap ratio between paralarval and juveniles and that between juveniles and subadultadults were 0.29 and 0.21, respectively ( Figure 7; Table 5).

Assessment of the Suitability of Using Squid Eye Lenses to Determine the Ontogenetic Variation in Stable Isotope Values
Loess curves were constructed to show the variations in the stable carbon and nitrogen isotope values during squid ontogenesis based on the LDs and the stable carbon and nitrogen isotope values that were found for each eye lens segment (Figure 8). It can be  For the western stock, the SEAc values for the paralarval, juvenile and subadult-adult squid were found to be 0.64‰ 2 , 1.52‰ 2 and 0.59‰ 2 , respectively. There was a trophic separation between the paralarval and the subadult-adult squid. The trophic overlap ratio between paralarval and juveniles and that between juveniles and subadult-adults were 0.29 and 0.21, respectively (Figure 7; Table 5).

Assessment of the Suitability of Using Squid Eye Lenses to Determine the Ontogenetic Variation in Stable Isotope Values
Loess curves were constructed to show the variations in the stable carbon and nitrogen isotope values during squid ontogenesis based on the LDs and the stable carbon and nitrogen isotope values that were found for each eye lens segment (Figure 8). It can be seen that, for both squid stocks, both the δ 13 C and δ 15 N values increased significantly as LD increased (δ 13 C: R 2 West = 0.142, p < 0.05; R 2 East = 0.082, p < 0.05; δ 15 N: R 2 West = 0.608, p < 0.05; R 2 East = 0.542, p < 0.05; Figure 8).

Discussion
In the study of the changes in isotopic values over the life cycle of the individual squid, trophic effects should be considered; however, it is also necessary to consider whether variations in baseline values due to geographic migration influence the values. Moreno et al. [63] considered that ignoring the baseline values or using a single baseline for multiple geographical regions would distort the result. Actually, a value of 3:1 for the ratio of δ N/δ C with increases in trophic level is widely used in trophic stable isotope studies [64]; also, assuming that the squid does not move to a different geographic area during their lifetime, the ratio of trophic increase should remain close to 3:1 theoretically [41]. However, in this study, only a few specimens were found that had a δ N/δ C ratio that was similar to the theoretical value (Table 3), which meant that, for these squids, the ontogenetic isotopic variation was related not only to trophic growth but also to geographic movement. In addition, based on the feeding analysis of O. bartramii by Watanabe (2004) [9] and Watanabe (2008) [65], the winter-spring and fall cohorts fed mainly on fishes of Myctophidae and the feeding proportion of fish increased gradually as the individuals grew. Moreover, the winter-spring cohort also consumed fishes such as Engraulis japonicas and Maurolicus imperatorius, and cephalopods such as Watasenia scillallans, while the fall cohort consumed cephalopods such as Berryteuthis anonychus, Gonatus spp., Abraliopsis spp. and Onychoteuthis borealijaponica. Among them, Myctophidae, E. japonicas and cephalopods have geographic migratory behavior. For example, the larvae and fry of Myctophidae in the western North Pacific are mainly distributed in the western transient zone, and when the length reaches 40mm, fishes migrate northward to the Oyashio Current and involve in the Western Subarctic Gyre [66,67]. Therefore, migratory prey maybe have an effect on the δ N/δ C ratio variation during the growth and migration of individual predators.
In this study, both the δ C and δ N values of the eastern and western squid stocks presented an increasing trend throughout the life cycle (Figure 8). Similar results have Loess curves showed that the δ 13 C values of the western stock of squid were higher than those of the eastern stock. For both stocks, the growth rates of δ 13 C values were similar in the early stage of the life cycle with the increase in LD. When LD was larger than about 8.95 mm, the δ 13 C values of the eastern stock tended to be flattened, while the growth rate of δ 13 C values decreased for the western stock after LD was larger than about 8.34 mm, but was still larger than that of the fall cohort ( Figure 8).
In contrast, when the LD was less than about 6.03 mm, the δ 15 N values of the eastern stock were greater than that of the western stock, while the δ 15 N values of the western stock increased at a greater rate compared with the eastern stock in the early stage of individual development; when the LD was greater than about 7.82 mm and 8.87 mm, respectively, both the δ 15 N values growth rate of the western and eastern stock decreased. However, the growth rate of the western squid stock remains higher than that of the eastern stock, which meant that the difference in the δ 15 N values between the two stocks gradually widened (Figure 8).

Discussion
In the study of the changes in isotopic values over the life cycle of the individual squid, trophic effects should be considered; however, it is also necessary to consider whether variations in baseline values due to geographic migration influence the values. Moreno et al. [63] considered that ignoring the baseline values or using a single baseline for multiple geographical regions would distort the result. Actually, a value of 3:1 for the ratio of δ 15 N/δ 13 C with increases in trophic level is widely used in trophic stable isotope studies [64]; also, assuming that the squid does not move to a different geographic area during their lifetime, the ratio of trophic increase should remain close to 3:1 theoretically [41]. However, in this study, only a few specimens were found that had a δ 15 N/δ 13 C ratio that was similar to the theoretical value (Table 3), which meant that, for these squids, the ontogenetic isotopic variation was related not only to trophic growth but also to geographic movement. In addition, based on the feeding analysis of O. bartramii by Watanabe (2004) [9] and Watanabe (2008) [65], the winter-spring and fall cohorts fed mainly on fishes of Myctophidae and the feeding proportion of fish increased gradually as the individuals grew. Moreover, the winter-spring cohort also consumed fishes such as Engraulis japonicas and Maurolicus imperatorius, and cephalopods such as Watasenia scillallans, while the fall cohort consumed cephalopods such as Berryteuthis anonychus, Gonatus spp., Abraliopsis spp. and Onychoteuthis borealijaponica. Among them, Myctophidae, E. japonicas and cephalopods have geographic migratory behavior. For example, the larvae and fry of Myctophidae in the western North Pacific are mainly distributed in the western transient zone, and when the length reaches 40mm, fishes migrate northward to the Oyashio Current and involve in the Western Subarctic Gyre [66,67]. Therefore, migratory prey maybe have an effect on the δ 15 N/δ 13 C ratio variation during the growth and migration of individual predators.
In this study, both the δ 13 C and δ 15 N values of the eastern and western squid stocks presented an increasing trend throughout the life cycle ( Figure 8). Similar results have been obtained for other hard structures found in O. bartramii, including beaks [68] and also obtained for gladii in jumbo squid (Dosidicus gigas) [69]. The δ 13 C values increased continuously for both the eastern and western squid stocks during the ontogenetic growth of O. bartramii, which implied that they tended to migrate from lower baseline δ 13 C values to higher baseline δ 13 C values. Bower and Ichii [70] found that the North Pacific population of O. bartramii could be further divided into fall cohort (hatching period mainly from September to February) and winter-spring cohort (hatching period mainly from January to May, while it could be extended to August). In addition, the fall cohort mainly spread east of 170 • E and near 160 • W, while the majority of the winter-spring cohort was located west of 170 • E [71]. In this study, based on the catching sites and hatching dates of the two stocks, it was finally determined that the western stock was the winter-spring cohort, while the eastern stock was the fall cohort. In fact, in the early stage of individual development, the STFZ where the fall cohort was located was close to the transition zone chlorophyll front (TZCF), thus bringing higher winter productivity for them, while the STD where the spawning ground of the winter-spring population was located was less productive. As the TZCF moved northward, the chlorophyll concentration in the TZ and SAFZ was higher, thus bringing greater productivity, so the migration of the winter-spring and fall cohorts from south to north during individual development was also a gradual movement to the high productivity zone [72]. Previous studies have suggested that male individuals of the fall cohort tended to stay in subtropical waters and did not migrate northward, but in this study, the sites where males were collected were located within the SAFZ, which meant that the males had some migration strategies [73]. However, due to the small sample size of male individuals, they cannot represent the whole stock. In addition, as noted above, migratory prey tends to consume different carbon sources, which could also contribute to the increase.
Furthermore, during ontogenesis, the δ 13 C values of the western stock were higher than the eastern stock ( Figure 8). The possible reason was that large-and mesoscale climatic and oceanic characteristics such as the Kuroshio Current and Oyashio Current dominate the influence over population dynamics of O. bartramii in the western North Pacific [74] and the Kuroshio and Oyashio Current are also significant sources of primary production for O. bartramii fishing grounds [75]; by contrast, in the eastern North Pacific, due to the offshore trade winds, upwelling occurs in the California Current and in the region, most of the year it occurs strong offshore surface transport [76]. In addition, Fang (2016) [77] showed that compared with the western stock, the eastern stock was farther from shore, which was assumed in this study to result in lower δ 13 C values in the eastern stock of O. bartramii during migration than the western stock.
Similar to the trends in the δ 13 C values, in this study, the δ 15 N values of both western and eastern stocks of O. bartramii were found to continuously increase from the paralarval to the subadult-adult stage (Figure 8). In the squid-occupied ocean, if the effects of trophic variation during ontogenetic growth were excluded, during the progressive formation of eye lenses, there would be a systematic decrease in δ 15 N values, which was often associated with habitat shifts during migration and was evidenced by a decrease in baseline δ 15 N values, according to Somes et al. (2010) [78]. However, it is contrary to the results of this study. In addition, there was a significant correlation between δ 15 N values and LD for both stocks in this study. Here, we suggest that as the squids undergo ontogenetic growth, their beak sizes get gradually larger, which means that they can feed on prey with greater body length and from higher trophic levels, including fish and other cephalopods [79,80]; thus, there is a shift from lower to higher trophic levels as they grow. Fang et al. (2016) [68] noted that the differences between two stocks are often related to differences in population growth rates, composition of feeding food and migration routes. When the LD was less than about 6.03 mm, the δ 15 N values of the eastern stock were greater than that of the western stock, and individuals in the LD range were basically in the paralarval-juvenile stage. Ichii et al. [19] showed that in the early stage of the life cycle, the growth rate of the fall cohort was faster compared to winter-spring populations, which therefore meant that they were able to feed on higher trophic level prey earlier and thus had higher δ 15 N values. In the later stages of individual development, the difference in δ 15 N values between the two stocks increased. While in the later stages of growth, the fall cohort grew less rapidly than the winter-spring cohort, which had an effect on δ 15 N values. In addition based on the diet composition consumed by the two stocks described above, it was found that although their trophic levels were essentially similar, the differences in feeding composition brought about by geographical differences between the two stocks could also have some influence on δ 15 N values, as O. bartramii are opportunistic, voracious predators. In addition, the variation in baseline values of different geographic regions during the migration of the two stocks would also be reflected in the variation of δ 15 N values.
Multiple comparisons showed that δ 13 C and δ 15 N values were significantly different between the paralarval-juvenile and subadult-adult stages of both stocks, but not between the paralarval and juvenile stage (Table 4). We suggest that the significant difference in δ 13 C values could be explained by geographic differences, as the habitat during the paralarval-juvenile stage was less diverse, squid were initially concentrated in the southern spawning ground and then started to move northward, while the significant differences in δ 15 N values were often related to differences in food composition in different areas and the size of ML. For eastern and western stocks, within each growth stage, we observed that the variability in the δ 15 N values was higher in the juvenile stage (western stock: 6.09‰; eastern stock: 3.84‰) ( Figure 6). It is usually considered that the range of δ 15 N values (NR) reflects the trophic length that is occupied by the squid: the larger the NR, the more trophic levels that are occupied [58]. O. bartramii is a voracious, opportunistic predator [81]; thus, at the juvenile stages, for the western stock, the Kuroshio Current and Oyashio Current bring rich and different trophic levels of food resources for the squid to obtain; for the eastern stock, the spawning ground has abundant feeding resources due to their proximity to TZCF [82]. However, later in the life cycle of the squid, as a result of the increase in mantle length and migration to its feeding grounds, O. bartramii may have a preference for larger prey from higher trophic levels. The breadth of the food resources then decreases, which causes the NR of O. bartramii in the later stages of the life cycle to decrease. This phenomenon can also be observed for some fish and other cephalopod species [83,84].
Niches not only reflect the spatial position occupied by a species within a certain environmental range at a certain time, but also reflect the influence of environmental and other factors on the species [85]. The niche breadth reflects the diversity or uniformity of species' feeding resources: the greater the breadth, the greater the environmental adaptability of the species. The niche overlap is another important indicator to measure the degree of competition and the similarity of resource utilization between or within species at different periods [86,87].
In this study, we focused on the niche variation of squid exogenous feeding after hatching. For both the eastern and western stocks, the SEAc values peak at the juvenile stage, which suggests that there is some variation in the resources used for ontogenetic growth [88]. Many species use similar resources when small [89]. At the paralarval stage, the swimming ability of O. bartramii is poor and their main source of food is plankton; thus, given the predation risk, the individual squid prefers to rely on similar resources and has a low adaptive capacity to the environment. Later, as a result of their ontogenetic growth, the squid tends to migrate to more open habitats, and the amount of diversity in the isotopic niches that they occupy increases [90]. According to previous studies, we found that, after the entry of the juvenile stage, the swimming ability of the squid improves as they migrate into the open water [91]. Their food sources then become more diverse due to the interaction between different ecosystems: their feeding targets are no longer limited to plankton and begin to shift toward small fish and other cephalopods [9,65]. This explains the larger SEAc values found for the juvenile stage of O. bartramii. For both stocks, the SEAc values first increase and then decline. Similar to the trends in NR discussed above, the SEAc values decrease in the later life stages. In their investigation of shark niches, Bethea et al. [92] pointed out that a reduction in the size of the niche may be due to larger individuals focusing more on specialized prey. Given the migratory activity of O. bartramii, this strengthens the conclusion that, in the later stages of their life cycle when they migrate to their feeding grounds, O. bartramii with larger ML focus more on larger prey (Figure 7).
The δ 13 C and δ 15 N values found for eye lenses suggest that for both the eastern and western stocks in the early stage of their life cycle, the amount of overlap in the SEAc between squid at the paralarval and juvenile stages was low (0.29 for the western stock and 0.24 for the eastern stock) ( Table 5). This was followed by a continuous decrease in the overlap between the juvenile and subadult-adult stages (western stock: juvenile/subadultadult: 0.21; eastern stock: juvenile/subadult-adult: 0.14), indicating that there is a gradual trophic niche separation between paralarval-juvenile and subadult-adult squid. This separation is related to changes in the O. bartramii habitat and influencing factors throughout its life cycle, exactly as Camarillo-Coop et al. [93] suggested that, as the size of the squid increases, so does the amount and variety of their prey, and that a gradual shift from small to large prey also leads to changes in the niche occupied by the squid [94].

Conclusions and Prospects
In this study, in the light of combined information from statolith determined age and stable eye lens isotopes, together with details of the squid life cycle and interactions with environment, we conclude the ontogenetic changes in their trophic preference and the O. bartramii collected were sexually immature, which meant that these individuals could not fully demonstrate the trophic variation that O. bartramii exhibited during the life cycle. In future studies, we will aim to carry out more comprehensive sampling to reduce the influence of this factor on the study results, and we will also aim to explore the trophic influence of North Pacific O. bartramii on other fishery organisms in the North Pacific Ocean, as well as the influence on the North Pacific ecosystem based on climatic and fishery conditions.

Data Availability Statement:
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.