- Open Access
Energy acquisition strategy for reproduction in a semelparous squid
Frontiers in Zoology volume 19, Article number: 28 (2022)
Energy demand for reproduction leads to a wide diversity of foraging and life-history strategy among wild animals, linking to a common objective to maximize reproductive success. Semelparous squid species in particular can use up to 50% of the total energy intake for reproduction. However, the energy acquisition strategy for reproduction is still a controversial issue regarding whether the squid shift in diet ontogenetically. Here we used Argentinean shortfin squid (Illex argentinus) as a case study to investigate the strategy of energy acquisition for reproduction, by analyzing energy density of the squid’s reproductive tissues including ovary, nidamental glands and oviduct eggs, and stable isotopes and fatty acids of the squid’s ovary.
The reproductive energy (the sum of the energy accumulated in ovary, nidamental glands and oviduct eggs) increased significantly with maturation. The ovary nitrogen stable isotopes (δ15N) showed a significant increase with maturation, but the increase by maturity stage was not equal to the typical enrichment of about 3‰ per trophic level. Isotopic niche width showed an increasing trend with maturation, and isotopic niche space exhibited greater overlap at advanced maturity stages. The relative amounts of 16:0, 20:5n3 and 20:4n6 in the ovary, tracing for carnivores and top predators, increased after the onset of maturation. The overall fatty acid profiles of the ovary showed significant differences among maturity stages, but obvious overlaps were found for mature squids. Mixed-effects model results revealed that reproductive energy was positively correlated with δ15N values. The reproductive energy was also positively related to the relative amounts of 18:0 and 20:4n6, respectively tracing for herbivores and top predators.
Our results validate that the squid shifts to feed on higher trophic prey for reproduction as energy demand increases once maturation commences. However, the squid does not shift feeding habits at a trophic level but instead broadens prey spectrum, coupled with increasing intake of higher trophic prey items, to meet the energy demand for reproduction. Such energy acquisition strategy may be selected by the squid to maximize reproductive success by balancing energy intake and expenditure from foraging, warranting future studies that aim to clarify such strategy for reproduction among semelparous species.
Energy demands scale differently with species growth and diet, and particularly the follow of reproductive energy leads to a wide diversity of life-history among wild animals . Behavioral and life-history strategies for an animal are selected for or against over evolutionary time to acquire energy from prey and manage energy budgets . At a coarse scale, these strategies are centered on maximizing reproduction to leave as many descendants as possible [3, 4]. Indeed, the most substantial metabolic expenditure for many species occurs during reproduction [5, 6]. This is particularly a case for semelparous animals, which on average, invest much more energy into a single reproductive episode before death than iteroparous animals that reproduce in multiple cycles over the course of their lifetime [1, 7]. However, due to the conservation of matter and energy where resources allocated to one function are not available for others, animals must face a fundamental trade-off in allocating limited resources to a basic life-history function such as survival and reproductive investment . Driven by the trade-off, energy demand and allocation to reproduction appear to interact with the physiological and metabolic constraints of an animal to give rise to the remarkable diversity of foraging and life-history strategies observed in nature .
Squid are one of the most important cephalopod species , playing an essential role in the trophic webs of marine ecosystems  and global production of marine fisheries . They are well-known for having a monocyclic life-history, characterized by a short lifespan, fast growth and semelparous reproduction . During maturation they usually allocate large amounts of energy for reproduction that can amount to 50% of their body mass . This allocation of energy influences the number and success of offspring, and ultimately determines population size and stability over time [13, 14]. To meet energy demands for reproduction, these species likely evolved an ontogenetic shift in diet, coupled with opportunistic and voracious foraging behavior [10, 15]. For example, Ommastrophid squids shift in diet from crustaceans to fish as they grow, in which there is a marked increase in the proportion of fish [16,17,18]. Since prey items such as fish occupy higher trophic positions, they have greater nutrient content  and greater energy density [20, 21]. It is expected that the uptake of such prey items would improve metabolic energy and support better life-history functions . However, higher trophic animals generally have better predator avoidance, which in turn increases the energy spent by the consumers on hunting . Such energy expenditures should theoretically be saved when the consumer forages for small prey species to gain a comparative energy resource . It is also assumed that the reason why squid species shift feeding habits ontogenetically is because they become better at hunting as they grow rather than the increased need of energy for reproduction [15, 25]. It is still unclear how reproductive energy needs are met for semelparous animals, which is important for understanding their life-history.
We used Argentinean shortfin squid Illex argetinus as the model system to investigate the energy acquisition strategy for the semelparous squids during maturation. This species represents a typical semelparous squid species, and is well documented for its fast growth, short life span and large amount of energy drawn for reproduction once maturation commences . Furthermore, this species not only plays a key role as a transient “biological pump” in the southwest Atlantic ecosystem , but it is also one of the most important commercial species in the world . In this study, we estimated the reproductive energy of I. argentinus by applying the technique of tissue energy density, and evaluated the energy acquisition strategy for reproduction using biochemical markers of stable isotopes and fatty acid profiles. The technique of tissue energy density has been successfully applied to investigate reproductive energy accumulation during maturation for I. argentinus [28, 29], as well as squid species Dosidicus gigas  and Sthenoteuthis oualaniensis . Generally, carbon (δ13C) and nitrogen (δ15N) stable isotopes are transferred from dietary sources to consumers in a predictable manner , where δ15N is typically enriched by about 3‰ per trophic level, and δ13C shows little change due to trophic transfer [32, 33]. Thus, the δ15N of a given consumer can be used to reflect the variations of prey, and the isotopic space can be used to indicate feeding niche widths of the predator . Fatty acids are essential components of all living organisms and passed with little to no modification from prey to predator, which make them useful tracers of diets and marine food-web structure . Analyzing the stable isotopes and fatty acids of reproductive tissues can infer the energy acquisition strategy for reproduction for the squid.
We randomly collected the Argentinean shortfin squid I. argetinus specimens from the catches of commercial jiggers on the high seas area of the Patagonian shelf in the southwest Atlantic Ocean (61° 09′ W ~ 62° 53′ W, 46° 08′ S ~ 47° 51′ S; Fig. 1) from January to March 2019. We aimed to investigate the strategy of energy acquisition when energy demand for reproduction is the greatest for the squid [28, 35]. We particularly aimed to clarify whether the squid exclusively feeds on prey items at a higher trophic or lower trophic position, or uses a combination, to meet reproductive energy demands. To further clarify the energy acquisition strategy, we also used linear mixed-effects models  to evaluate the relationship between reproductive energy and stable isotopic values and fatty acid profiles. These analyses will allow for a more comprehensive understanding of the life-history strategy in squid species and demonstrate the potential application of these methods to other semelparous animals in nature.
The reproductive system of I. argetinus increases dramatically once maturation commences , therefore we selected 68 female specimens at maturity stages from physiologically maturing (stage III) to spawning (stage VI) after the scales proposed by ICES  and Lin et al. , for reproductive energy analysis, carbon (δ13C) and nitrogen (δ15N) stable isotope analysis, and fatty acids analyses. The sample size was 12 at stage III (physiologically maturing), 15 at stage IV (early physiologically mature), 15 at stage V (late physiologically mature), 15 at stage VI (functionally mature), and 11 at stage VII (spawning). The squids' measured dorsal mantle length ranged from 184 to 262 mm and weighed body mass from 129 to 374 g. Both mantle length and body weight exhibited unimodal distribution (Additional file 1: Fig. S1a–b), indicating that the squids were from a single population stock . Body weight was significantly correlated with mantle length, with a power function W = 0.00013 × L2.66 (r2 = 0.67, p < 0.05; Additional file 1: Fig. S1c). There were no significant differences in mantle length or body weight among maturity stages (mantle length, F = 0.92, p = 0.46; body weight, F = 1.10, p = 0.37; Additional file 1: Fig. S2).
The energy density of ovary, nidamental glands and oviduct eggs (henceforth “eggs”) was determined using an automatic isoperibol calorimeter. The average energy density over all maturity stages was 23.44 ± 0.89 kJ/g dry mass for ovaries, 25.38 ± 0.40 kJ/g dry mass for eggs, and 19.93 ± 0.48 kJ/g dry mass for nidamental glands, respectively. Energy density in the ovary increased significantly with maturation (F = 12.57, p < 0.05). Greater energy density within ovaries was observed at stages IV-VI; while both eggs and nidamental glands remained relatively constant (Fig. 2a–c). The reproductive energy, estimated as the sum of the energy accumulated in the ovary, nidamental glands and eggs, increased significantly with maturation (F = 30.31, p < 0.05), with a five-fold increase seen from stage III (74.12 ± 31.07 kJ) to stage VI (372.79 ± 82.85 kJ) (Fig. 2d).
The δ15N values of the ovary were determined from 13.03 to 16.64 ‰, and exhibited a significant increase trend with maturation (F = 7.77, p < 0.05), reaching a peak value of 15.05 ± 0.68 ‰ at stage VI (Table 1; Fig. 3a). However, the increase of δ15N values by maturity stage was not equal to a typical enrichment of about 3‰ per trophic level . The δ13C values varied from − 18.03‰ to − 15.71‰, but did not vary significantly between maturity stages (F = 0.71, p = 0.58), with an overall mean value of − 17.26 ± 0.59 ‰ (Table 1; Fig. 3b).
Both δ15N and δ13C values showed greater variability (standard deviations [SD]) with maturation (Table 1). The corrected standard ellipse area (SEAc), calculated as the standard ellipse area of isotopic covariance matrix corrected by sample size using Bayesian approach , ranged from 0.57 at stage III to 1.32 at stage VI (Table 1) and revealed the increase of trophic niche width with maturation. In addition, the Bayesian standard ellipse area (SEAb), a posterior estimate of isotopic matrix using the Bayesian approach , was significantly different between maturity stages (Kruskal–Wallis, χ2 = 95.51, p < 0.05; Fig. 3c), further revealing the differences in the isotopic niches for the squids between maturity stages. The overlap in the isotope niches was estimated as the ratio of overlap proportion and the sum of the non-overlapping area between each consecutive SEAc . This ellipse plots with 40% Bayesian credible intervals was almost separated between stages III and IV however, was obvious overlap after reaching stage IV, with the largest overlap occurring between stages IV and V (Table 1). This observation was confirmed by the overlap of the SEAc of isotopic data for each maturity stage (Fig. 3d).
Fatty acid profiles
A total of 28 fatty acids (FAs) were determined in the squid’s ovaries, 14 of which had relative mean values greater than 0.5% (Additional file 1: Table S1). Total content of fatty acids (total FAs) increased significantly from stages III–IV, and the 14 fatty acids accounted for 90.06% ~ 97.97% of total FAs, with an average 97.04 ± 1.34%. Nine out of the 14 FAs showed significant differences in the relative amounts between maturity stages, where an increasing trend was found for 14:0 (F = 20.41, p < 0.05), 20:2 (F = 4.17, p = 0.006), 20:3n3 (F = 3.08, p = 0.024), 20:4n6 (F = 4.48, p = 0.004) and 20:5n3 (F = 3.20, p = 0.028), a decreasing trend for 18:0 (F = 4.31, p = 0.005), 18:1n9 (F = 3.53, p = 0.013) and 20:1 (F = 18.37, p < 0.05), and an increase for 16:0 from stages III–VI (F = 18.18, p < 0.05), followed by a decrease at stage VII (Fig. 4a; Additional file 1: Table S1). The main fatty acid classes, namely saturated fatty acids (SFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA), were also significantly different in the relative amounts between maturity stages. SFA increased from stages III–VI (F = 8.87, p < 0.05), PUFA showed an increasing trend with maturation (F = 3.56, p = 0.035), while MUFA decreased from stages III–VI (F = 14.15, p < 0.05) (Fig. 4b; Additional file 1: Table S1).
The analysis of similarity (ANOSIM) revealed that the overall FA profiles showed significant differences among maturity stages (global R-value = 0.35, p = 0.001). However, the significant differences were only found between stages III and IV (global R-value = 0.30, p = 0.004) and between stages VI and VII (global R-value = 0.37, p = 0.001) (Table 2). Non-metric multidimensional scaling (NMDS) results showed a large amount of overlap in the overall FA profiles between stages IV, V and VI, and obvious segregations between stages III and IV and between stages VI and VII (Fig. 4c). FAs 18:0, 18:1n9 and 20:1 exhibited the highest values at stage III, 14:0 and 16:0 at stages V and VI, and 20:2, 20:3n3 and 20:4n6 at stage VII (Fig. 4c; Additional file 1: Table S1).
Relationship between reproductive energy and nitrogen isotopes and fatty acids
The relationships between reproductive energy and δ15N and FAs were separately carried out by applying linear mixed-effects models (LMMs), where mantle length and sea surface chlorophyll-a concentration were included as the explanatory variables, and maturity stage was included as a random effect to account for the potential effect caused by sexual development. Results from the LMMs revealed that the reproductive energy was positively correlated with δ15N values (t-value = 3.02, p = 0.004; Fig. 5a; Additional file 1: Table S2). Similarly, the reproductive energy was also positively related to the relative amounts of 14:0 (t-value = 2.97, p = 0.005), 18:0 (t-value = 2.13, p = 0.04), 20:4n6 (t-value = 2.60, p = 0.01) (Fig. 5b–d, Additional file 1: Table S3), where 18:0 is assumed to be the trophic marker for herbivores , and 20:4n6 for top predators . These findings suggest that the larger amount of reproductive energy is obtained by increasing intake of higher trophic prey items. However, there was no significant relationship between reproductive energy and mantle length (p > 0.05; Additional file 1: Table S2, S3). The reproductive energy was also not significantly correlated with the sea surface chlorophyll-a concentration (p > 0.05; Additional file 1: Table S2, S3), although the sea surface chlorophyll-a concentration in the sampling area varied significantly between weeks (F = 850.69, p < 0.05; Additional file 1: Figure S3).
The energy density of the squid ovary increased as maturation progressed, while both eggs and nidamental glands were relatively stable throughout sexual maturation. This pattern is consistent with previous research conducted on the same species [28, 29] and other squid species (e.g. D. gigas ; S. oualaniensis ). The increase of energy density in the ovary may be related to yolk production during ovarian development . Stable energy density of eggs is an indication of similar quality egg production, which must be of great importance for these semelparous species to maximize reproductive success . Consequently, the reproductive energy increased significantly with maturation (Fig. 1d), most likely in relation to the enlargement of reproductive systems . More importantly, the increase of reproductive energy should be considered as the increase of energy demand for reproduction. Illex argentinus has been found to significantly increase energy allocation to reproduction once maturation commences [28, 29]. Many monocyclic species have also been reported to direct large amounts of energy to reproduction at the onset of maturation [12, 43, 44]. This is most likely an evolutionary selection given the semelparous nature of the species and the need to maximize reproductive success.
Energy acquisition strategy for reproduction
The squid adopts an energy acquisition strategy where individuals feed on higher trophic prey items leading to an increase of energy in reproductive organs during maturation. Such phenomenon is supported by the results of LMMs that were performed the reproductive energy on either the nitrogen isotopes or the selected FA profiles, where the reproductive energy is positively related to δ15N values (Fig. 5a), and the relative amounts of carnivorous fatty acid tracer 18:0 and top predator fatty acid tracer 20:4n6 (Fig. 5b–c). Given the fact that energy allocation to reproduction is dramatically increasing once sexual development commences , the intake of higher trophic prey items is likely necessary to meet the energetic demand following the hypothesis that prey species at higher trophic position have a greater nutrient content .
The squid shifts feeding habits with maturation, especially during the course from stages III–IV, evidenced by the increase trend of δ15N values (Fig. 3a). Further evidence can be provided by the fatty acids analyses (Fig. 4). The relative amounts of 16:0 and 20:5n3, dietary tracers for carnivores [40, 45], increase along with maturation, while the relative amounts of 18:0 and 20:1, tracing for herbivores , decrease after the onset of maturation. Moreover, the relative amounts of 20:4n6, a tracer for top predator , also shows a statistically significant increasing trend throughout sexual maturation. Since consumers are subject to biochemical limitations in biosynthesizing dietary fatty acids , it is reasonable to expect that the squid may decrease the intake of prey species at lower trophic level, and turn to higher trophic preys after the onset of maturation. This finding is consistent with the results of stomach content analyses where the squids change in diet composition as they grow, switching from crustaceans at small sizes to fishes and cephalopods at large sizes [47,48,49,50].
However, the squid does not shift prey items at a trophic level completely. The enrichment of δ15N values between each two consecutive maturity stages is not typically enriched by about 3‰ per trophic level . By contrast, the variance of δ15N values at each maturity stage increased with maturation, being the greatest variance occurred just before spawning stage (Table 1). The Bayesian isotopic niche analyses confirmed that the squid at more advanced maturity stage has greater isotopic niche width, indicated by the increase of SEAc and SEAb (Table 1; Fig. 3c). Given the isotopic niche linking to prey , these observations indicate that the squid increases prey spectrum to meet energy demands for reproduction . The maturing squids showed larger amounts of herbivorous FA tracers 18:0 and 20:1 , while mature and spawning squids significantly increased the relative amounts of 16:0 and 20:4n6, both which have been used as indicators of carnivores and top predators respectively [40, 41]. The results might expect that maturing squids primarily feed on herbivores, and mature and spawning ones switch to feed on carnivorous prey items. It is noted that the mature squids may prey on similar prey items, evidenced by the overlaps in the isotopic niche (Fig. 3d) and FA profiles (Fig. 4c). Comparing to small prey species, the one at higher trophic position has greater energy density [20, 21], and would improve the metabolic energy of consumers . However, prey species at higher trophic position may perform better to avoid predation . Consumers must trade-off energy acquisition and the cost of capturing prey . Therefore, the squid may balance energy expenditure and gain to maximize energy intake through ontogenetic shift in diet, coupled with broadening prey spectrum.
The strategy of energy acquisition for reproduction is independent from body size that has been considered as an important factor linked to feeding ability [15, 25]. LMMs results indicated that reproductive energy was not correlated with mantle length (Additional file 1: Tables S2, S3). This seems contrast to previous reports that larger individuals invest much more reproductive energy , and many species do exhibit reproductive energy output scaling with body size . However, most of the studies linearized reproductive investment to body size without considering the potential effects of maturity state. In squid species, notably, somatic growth is gradually decreasing after reaching maturation , while the reproductive system expands dramatically, especially during physiologically maturing stage . Thus, the influence of body size on reproductive energy acquisition may become less important after the onset of maturation in squids.
Resource availability is also an important factor that influences consumer’s feeding habits [52, 53]. For squid species, many of them have been reported that feeding behavior is related to the availability of food resources [17, 53], potentially linking to their opportunistic life-history and voracious foraging behavior . In the present study, we observed significant variations in the weekly surface chlorophyll-a concentration in the sampling area (Additional file 1: Fig. S1), suggesting that the primary production are temporal variable . Unexpectedly, the reproductive energy was not related to the surface chlorophyll-a concentration (Additional file 1: Tables S2, S3). The southwestern Atlantic system has high biological production because of the high variability of physio-chemical and biological attributes occurring at the shelf and slope water columns . Therefore, the prey community can be resilient to the temporal change of productivity and maintain a status of dynamic equilibrium . Indeed, the sampling stations were apparently not situated in water columns with high chlorophyll-a concentration (Additional file 1: Fig. S3a). This case can indirectly provide evidence to explain the non-significant relationship between reproductive energy and chlorophyll-a concentration.
The semelparous Argentinean shortfin squid invests large amounts of energy for reproduction during maturation. In order to meet the energy demand, the squid ontogenetically shifts feeding habits from lower to higher trophic prey items during maturation. However, the shift is not equal to a complete trophic level. Instead, the squid broadens the prey spectrum, coupled with increasing intake of higher trophic prey. Energy acquisition for reproduction was independent of body size and primary production. This strategy of energy acquisition for reproduction may be evolutionarily selected by the squid to maximize reproductive values by balancing energy intake and expenditure from foraging. To our best knowledge, this work is the first study that uses multiple techniques including energy density, stable isotopes and fatty acid biomarkers to investigate energy acquisition strategy among semelparous species. The results put forward our understanding of the squid’s life-history in terms of energy acquisition for reproduction. Future studies should aim to confirm whether the strategy identified in this study holds for other semelparous species, which have common monocyclic life-histories but can vary in breeding patterns and therefore, reproductive strategies.
The squid specimens were collected over the Patagonian shelf of the southwest Atlantic Ocean (Fig. 1). It is well known that the southwest Atlantic is characterized by the presence of the Brazil Current (BC), a southward-flowing warm western boundary current, and the northward-flowing Malvinas Current (MC), in which both currents encounter at approximately 38°S . Meanwhile, intense frontal transitions at various near shore locations and along the shelf break promote vertical circulations that inject nutrients into the upper layer, leading to the enhanced growth of phytoplankton . Consequently, the waters along and across the shelf and slope are inhabited by diverse organisms, including mammal, bird, fish and cephalopod species in all life stages .
Specimens and process
Squid specimens were obtained from the catches of commercial jigging fishery that operated in the high seas of southwest Atlantic (61° 09′ W ~ 62° 53′ W, 46° 08′ S ~ 47° 51′ S; Fig. 1) from January 8th to March 4th in 2019. Specimens were randomly collected onboard, and frozen immediately at – 30 °C for laboratory experiments. After defrosted at room temperature in laboratory, each specimen was taken dorsal mantle length (ML) to the nearest 1.0 mm and body weight (BW) to the nearest 1.0 g. Each specimen was dissected to identify sex and assigned a macro-scale maturity stage. The macro-scale maturity stage was after the scales proposed by ICES  and Lin et al. : I immature, II developing, III physiologically maturing, IV–V, physiologically mature, VI, functionally mature, VII, spawning, VIII, spent. 152 specimens were sexed and assigned the maturity stages. Preliminary analysis indicated that about one fifth of the specimens were at stages I–II, another one fifth at stage III, and the remaining at stages IV–VII. No spent individuals were found, probably because spent individuals die shortly after spawning .
Since the squid develops reproductive system from stage III , we randomly selected 68 female specimens at maturity stages from III to VII for energy density analysis, stable isotope analysis and fatty acids analyses (12 at stage III, 15 at stage IV, 15 at stage V, 15 at stage VI, and 11 at stage VII). For each selected specimen, the ovary, nidamental glands and eggs were taken and lyophilized to constant weight in a freeze-dried chamber (Scientz-10 N lab lyophilizer, Ningbo Scientz Biotechnology Co., LTD.). The dry weight of each tissue was weighed to the nearest 0.1 mg, and then ground into fine powder using Scientz-48 grinder (Ningbo Scientz Biotechnology Co., LTD.).
An approximately 1.5–4.0 g of powdered tissue was used to determine the energy density (kJ g−1) of the ovary, nidamental glands and eggs of each squid specimen, separately, using an automatic isoperibol calorimeter (Model 6400, Parr Instrument Company, Moline, IL, USA). The powdered tissue was gently added to a capsule and then placed into the capsule holder of the calorimeter, which allows for automatically determining the energy density within several minutes.
The absolute energy (kJ) of a given tissue was calculated as the energy density multiplied by the dry tissue weight. Then, the reproductive energy was estimated as the sum of absolute energy of the ovary, nidamental glands and eggs.
Fatty acids analyses
An approximately 1 g of powdered ovary was used to extract lipids with a methanol: dichloromethane: water solvent mixture (20:10:8 by volume) according to Bligh and Dyer method . The extracted lipids were stored at − 20 °C in methanol: dichloromethane (2:1 by volume) with 0.01% butylated hydroxytoluene (BHT) as antioxidant for fatty acid methyl esters (FAME) analysis within 24–48 h. The lipid extracted tissue was recycled and lyophilized to constant weight for stable isotope analysis.
A modification of GAQSIQ method  was used to analyze the FAME of the squid’s ovary. The method has been successfully tested in soma, ovary and digestive gland samples of I. argentinus  and Dosidicus gigas . FAME were determined using an Agilent 7890B Gas Chromatography (GC) coupled to a 5977A series Mass Spectrometer Detector (MSD, Agilent Technologies, Inc. USA), equipped with a fused silica 60 m × 0.25 nm open tubular column (HB-88: 0.20 μm, Agilent Technologies, Inc. USA). Individual FAME was identified by the retention times and mass spectra with a known concentration internal standard 19:0 (GLC 37, Nu-Chek Prep, Inc.). The separation was carried out with helium as the carrier gas, and a thermal gradient programed from 125 to 250 °C, with the auxiliary heater at 280 °C.
Total content of fatty acids (total FAs) was determined as milligram per gram of dry mass (mg/g dry mass), based on the mass fraction of FAMEs relative to the internal standard 19:0 . Individual fatty acid (FA) was expressed as percentage of total FAs. Individual fatty acid that accounted for less than 0.5% of the total FAs was excluded from statistical analyses below. The individual fatty acid was also grouped into saturated fatty acids (SFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA), which were also expressed as percentage of total FAs.
Stable isotope analysis
Lipid fraction in soft tissues tends to be depleted in 13C isotope, and lipid extraction has become a standard procedure in stable isotope analysis [62, 63]. The lipid-extracted ovary samples (see “Fatty acids analysis”) were used to determined carbon (δ13C) and nitrogen (δ15N) stable isotope values. A powdered sub-sample of ovary tissue (≈ 0.3 mg) was placed into tin capsule and analyzed in a SerCon Integra 2 integrated elemental analyzer and an isotope ratio mass spectrometer (EA-IRMS) at the Stable Isotope Core Laboratory in Third Institute of Oceanography (Ministry of Natural Resources, China). Isotopic values are reported in standard δ-notion in parts per thousand (‰), where δ13C or δ15N = [(Rsample/Rstandard) − 1] × 1000, with Rsample and Rstandard representing the ratios of 13C/12C and 15 N/14 N of the sample and the standard reference material, respectively. The reference material was Vienna Pee Dee Belemnite (VPDB) for carbon and atmospheric nitrogen (N2) for nitrogen. The measurement errors were approximately 0.02‰ and 0.02‰ for δ13C and δ15N, respectively.
Significant differences (p < 0.05) among maturity stages were tested for the tissue energy density, reproductive energy, stable isotopes, total fatty acids content and relative amount of individual fatty acid. All data were checked for normality of distribution with the one-sample Kolmogorov–Smirnov test and for homogeneity of the variances with the Levene’s test . When the normality was satisfied, one-way analyses of variance (ANOVA) was applied to determine the significant difference between maturity stages, followed by Tukey’s honestly significant post-hoc test. When normal distribution and/or homoscedasticity were not achieved, data were subjected to a Kruskall–Wallis nonparametric test and a Games–Howell post hoc test was performed .
To determine the energy acquisition strategy during maturation, we evaluated the ovary isotopic niche based on Bayesian inference, and assessed the similarity of ovary fatty acid profiles. The ovary isotopic niche can reflect the dietary niche characteristics of the squid , and the ovary fatty acid profiles can be used to evaluate the potential similarity of dietary items . We calculated the isotopic niche widths for each maturity stage using Stable Isotope Bayesian Ellipses in R (SIBER, ), including the Bayesian estimate of standard ellipse area (SEAb), the corrected standard ellipse area (SEAc) and the isotopic niche overlap estimated as the ratio of overlap proportion and the sum of non-overlap areas between each consecutive SEAc (non-overlap SEAc proportion) based on 1000 replications.
We applied analysis of similarity (ANOSIM) to examine differences in the overall fatty acid profiles between maturity stages. In addition, non-metric multidimensional scaling (NMDS) was applied to visualize the differences in the overall FA profiles. As individual FA was expressed as percentage of total FAs, a square-rooted transformation was used to avoid over-emphasis of extreme values . A Bray–Curtis dissimilarity matrices was employed in the ANOSIM and NMDS . The analyses were performed in the package ‘vegan’ in R platform .
We further applied linear mixed-effects models (LMMs)  to investigate the energy acquisition strategy for reproduction in I. argentinus. LMMs were respectively performed the reproductive energy on δ15N values and selected fatty acid profiles, using maturity stage to account for sexual development effects and unexplained differences among individuals from different maturity stages. The selected FA profiles were those FAs that were found significant differences between maturity stages. Log-transformed mantle length of the squid and sea surface chlorophyll-a concentration (Chla) were also considered as explanatory variables to account for the effects of body size and primary production, respectively. The Chla data were downloaded from NOAA (https://oceanwatch.pifsc.noaa.gov/erddap/). To reduce the bias of daily variation, we used weekly Chla (Dataset ID: noaa_snpp_chla_weekly). LMMs were performed in the package ‘lme4’ in R platform .
Availability of data and materials
All data that support the findings of this study are included in the manuscript and supplementary material. Data are however used under license from the Distant Squid Fisheries Sci-Tech Group (SHOU), and the usage is only permitted from the authors and Distant Squid Fisheries Sci-Tech Group (SHOU) upon reasonable request.
McBride RS, Somarakis S, Fitzhugh GR, Albert A, Yaragina NA, Wuenschel MJ, et al. Energy acquisition and allocation to egg production in relation to fish reproductive strategies. Fish Fish. 2015;16:23–57.
Llodra ER. Fecundity and life-history strategies in marine invertebrates. Adv Mar Biol. 2002;43:87–170.
Marshall DJ, Keough MJ. The evolutionary ecology of offspring size in marine invertebrates. Adv Mar Biol. 2007;53:1–60.
Smith C, Wootton RJ. The remarkable reproductive diversity of teleost fishes. Fish Fish. 2015;17:1208–15.
English S, Bonsall MB. Physiological dynamics, reproduction-maintenance allocations, and life history evolution. Ecol Evol. 2019;9:9312–23.
Clarke A. Energy flow in growth and production. Trends Ecol Evol. 2019;34:502–9.
Johns ME, Warzybok P, Bradley RW, Jahncke J, Lindberg M, Breed GA. Increased reproductive investment associated with greater survival and longevity in Cassin’s auklets. Proc Royal Soc B. 2018;285:20181464.
Stearns SC. The evolution of life histories. London: Oxford University Press; 1992.
Stephens PA, Houston AI, Harding KC, Boyd IL, McNamara JM. Capital and income breeding: the role of food supply. Ecology. 2014;95:882–96.
Boyle P, Rodhouse P. Cephalopods: ecology and fisheries. Oxford: Wiley-Blackwell; 2005.
FAO. FAO yearbook: fishery and aquaculture statistics 2018/FAO annuaire. Rome: FAO; 2020.
Jackson GD, Semmens JM, Phillips KL, Jackson CH. Reproduction in the deepwater squid Moroteuthis ingens, what does it cost? Mar Biol. 2004;145:905–16.
Calow P. Fact and theory: an overview. In: Boyle PR, editor. Cephalopod life cycles vol II comparative reviews. London: Academic Press; 1987. p. 351–65.
Boggs CL. Dynamics of reproductive allocation from juvenile and adult feeding: radiotracer studies. Ecology. 1997;78:192–202.
Rosas-Luis R, Navarro J, Martínez-Baena F, Sánchez P. Differences in the trophic niche along the gladius of the squids Illex argentinus and Doryteuthis gahi based on their isotopic values. Reg Stud Mar Sci. 2017;11:17–22.
Ruiz-Cooley RI, Markaida U, Gendron D, Aguíñiga S. Stable isotopes in jumbo squid (Dosidicus gigas) beaks to estimate its trophic position: comparison between stomach contents and stable isotopes. J Mar Biolog Assoc. 2006;86:437–45.
Martínez-Baena F, Navarro J, Albo-Puigserver M, Palomera I, Rosas-Luis R. Feeding habits of the short-finned squid Illex coindetii in the western Mediterranean Sea using combined stomach content and isotopic analysis. J Mar Biolog Assoc. 2016;96:1235–42.
Gong Y, Ruiz-Cooley RI, Hunsicker ME, Li Y, Chen X. Sexual dimorphism in feeding apparatus and niche partitioning in juvenile jumbo squid Dosidicus gigas. Mar Ecol Prog Ser. 2018;607:99–112.
Marshall KN, Essington TE. Energetic conditions promoting top-down control of prey by predators. PLoS ONE. 2011;6: e29723.
Schaafsma FL, Cherel Y, Flores H, van Franeker JA, Lea M-A, Raymond B, et al. Review: the energetic value of zooplankton and nekton species of the southern ocean. Mar Biol. 2018;165:129.
Dessier A, Dupuy C, Kerric A, Mornet F, Authier M, Bustamante P, et al. Variability of energy density among mesozooplankton community: new insights in functional diversity to forage fish. Prog Oceanogr. 2018;166:121–8.
García S, Domingues P, Navarro JC, Hachero I, Garrido D, Rosas C. Growth, partial energy balance, mantle and digestive gland lipid composition of Octopus vulgaris (Cuvier, 1797) fed with two artificial diets. Aquac Nutr. 2011;17:e174–87.
Burger JR, Hou C, Brown JH. Toward a metabolic theory of life history. PNAS. 2019;116:26653–61.
Bhat U, Kempes CP, Yeakel JD. Scaling the risk landscape drives optimal life-history strategies and the evolution of grazing. PNAS. 2020;117:1580–6.
Gong Y, Li Y, Chen X, Yu W. trophic niche and diversity of a pelagic squid (Dosidicus gigas): a comparative study using stable isotope, fatty acid, and feeding apparatuses morphology. Front Mar Sci. 2020;7:642.
Rodhouse PGK, Arkhipkin AI, Laptikhovsky V, Nigmatullin C, Waluda CM. Illex argentinus Argentine shortfin squid. In: Rosa R, Pierce G, O’Dor R, editors. Advances in squid biology ecology and fisheries part II - oegopsid squids. New York: Nova Science Publishers; 2013. p. 109–48.
Arkhipkin AI. Squid as nutrient vectors linking Southwest Atlantic marine ecosystems. Deep-Sea Res. 2013;II(95):7–20.
Lin D, Chen X, Wei Y, Chen Y. The energy accumulation of somatic tissue and reproductive organs in post-recruit female Illex argentinus and the relationship with sea surface oceanography. Fish Res. 2017;185:102–14.
Lin D, Chen X, Wei Y, Xu M. Energy accumulation of both soma and reproductive organs and its allocation to reproduction in Argentinean short-fin squid Illex argentinus. J Fish China. 2017;41:70–80.
Han F, Chen X, Lin D. Energy accumulation and its relation to sea surface environments in Dosidicus gigas from the equatorial water of the Eastern Pacific. J Fish Sci China. 2020;27:427–37.
Zhu K, Yao J, Chen X, Liu W, Lin D. Energy accumulation of both somatic and reproductive tissues and the allocation to reproduction in the dwarf form individuals of Sthenoteuthis oualaniensis in the South China Sea. J Shanghai Ocean Univ. 2020;29:910–20.
Layman CA, Araujo MS, Boucek R, Hammerschlag-Peyer CM, Harrison E, Jud ZR, et al. Applying stable isotopes to examine food-web structure: an overview of analytical tools. Biol Rev. 2012;87:545–62.
Post DM. Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology. 2002;83:703–18.
Iverson SJ. Tracing aquatic food webs using fatty acids: from qualitative indicators to quantitative determination. In: Kainz M, Brett MT, Arts MT, editors. Lipids in aquatic ecosystems. New York: Springer; 2009. p. 281–308.
Lin D, Han F, Xuan S, Chen X. Fatty acid composition and the evidence for mixed income–capital breeding in female Argentinean short-fin squid Illex argentinus. Mar Biol. 2019;166:90.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
ICES. Report of the workshop on sexual maturity staging of cephalopods, 8–11 November 2010, ICES CM 2010/ACOM:49. Livorno, Italy: ICES, 2010.
Lin D, Chen X, Fang Z. Preliminary study on reproductive biology of summer spawning stock of Illex argentinus in the southwestern Atlantic Ocean. J Fish China. 2014;38:843–52.
Jackson AL, Inger R, Parnell AC, Bearhop S. Comparing isotopic niche widths among and within communities: SIBER – stable isotope Bayesian ellipses in R. J Anim Ecol. 2011;80:595–602.
El-Sabaawi R, Dower JF, Kainz M, Mazumder A. Characterizing dietary variability and trophic positions of coastal calanoid copepods: insight from stable isotopes and fatty acids. Mar Biol. 2009;156:225–37.
Colombo SM, Wacker A, Parrish CC, Kainz MJ, Arts MT. A fundamental dichotomy in long-chain polyunsaturated fatty acid abundance between and within marine and terrestrial ecosystems. Environ Rev. 2017;25:163–74.
Lin D, Chen X, Chen Y, Zhu G, Yan J. Ovarian development in Argentinean shortfin squid Illex argentinus: group-synchrony for corroboration of intermittent spawning strategy. Hydrobiologia. 2017;795:327–39.
McGrath B, Jackson G. Egg production in the arrow squid Nototodarus gouldi (Cephalopoda: Ommastrephidae), fast and furious or slow and steady? Mar Biol. 2002;141:699–706.
Bowerman TE, Pinson-Dumm A, Peery CA, Caudill CC. Reproductive energy expenditure and changes in body morphology for a population of Chinook salmon Oncorhynchus tshawytscha with a long distance migration. J Fish Biol. 2017;90:1960–79.
Pethybridge HR, Nichols PD, Virtue P, Jackson GD. The foraging ecology of an oceanic squid, Todarodes filippovae: the use of signature lipid profiling to monitor ecosystem change. Deep-Sea Res. 2013;II(95):119–28.
Schmidt-Nielsen K. Animal physiology: adaptation and environment. 5th ed. Cambridge, UK: Cambridge University Press; 1997.
Ivanovic M, Brunetti NE. Food and feeding of Illex argentinus. Antarct Sci. 1994;6:185–93.
Santos RA, Haimovici M. Food and feeding of the short-finned squid Illex argentinus (Cephalopoda: Ommastrephidae) off southern Brazil. Fish Res. 1997;33:139–47.
Rosas-Luis R, Sánchez P, Portela JM, del Rio JL. Feeding habits and trophic interactions of Doryteuthis gahi, Illex argentinus and Onykia ingens in the marine ecosystem off the Patagonian Shelf. Fish Res. 2014;152:37–44.
Chang S-W, Chen R-G, Liu T-H, Lee Y-C, Chen C-S, Chiu T-S, Ko C-Y. Dietary Shifts and Risks of Artifact Ingestion for Argentine Shortfin Squid Illex argentinus in the Southwest Atlantic. Front Mar Sci. 2021;8:675560. https://doi.org/10.3389/fmars.2021.675560.
Barneche DR, Robertson DR, White CR, Marshall DJ. Fish reproductive-energy output increases disproportionately with body size. Science. 2018;360:642–5.
Houston AI, Higginson AD, McNamara JM. Optimal foraging for multiple nutrients in an unpredictable environment. Ecol Lett. 2011;14:1101–7.
Portner EJ, Markaida U, Robinson CJ, Gilly WF. Trophic ecology of Humboldt squid, Dosidicus gigas, in conjunction with body size and climatic variability in the Gulf of California. Mexico Limnol Oceanogr. 2020;65:732–48.
Lin D, Chen X. Top predator reveals the stability of prey community in the western subarctic Pacific. PLoS ONE. 2020;15: e0234905.
Piola AR, Palma ED, Bianchi AA, Castro BM, Dottori M, Guerrero RA, et al. Physical oceanography of the SW Atlantic shelf: a review. In: Hoffmeyer MS, Sabatini ME, Brandini FP, Calliari DL, Santinelli NH, editors., et al., Plankton ecology of the Southwestern Atlantic: from the subtropical to the Subantarctic realm. Cham: Springer International Publishing; 2018. p. 37–56.
Cepeda GD, Temperoni B, Sabatini ME, Viñas MD, Derisio CM, Santos BA, et al. Zooplankton communities of the Argentine continental shelf (SW Atlantic, ca. 34°–55°S), An Overview. In: Hoffmeyer MS, Sabatini ME, Brandini FP, Calliari DL, Santinelli NH, editors., et al., Plankton ecology of the southwestern Atlantic: from the subtropical to the Subantarctic realm. Cham: Springer International Publishing; 2018. p. 171–99.
Campos EJD, Miller JL, Müller TJ, Peterson RG. Physical oceanography of the southwest Atlantic Ocean. Oceanography. 1995;8:87–91.
Laptikhovsky VV, Nigmatullin CM. Egg size, fecundity, and spawning in females of the genus Illex (Cephalopoda: Ommastrephidae). ICES J Mar Sci. 1993;50:393–403.
Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37:911–7.
GAQSIQ. Determination of total fat, saturated fat, and unsaturated fat in foods: hydrolytic extraction-Gas chromatography (GB/T 22223-2008). Beijing: Standards Press of China; 2008.
Chen X, Han F, Zhu K, Punt AE, Lin D. The breeding strategy of female jumbo squid Dosidicus gigas: energy acquisition and allocation. Sci Rep. 2020;10:9639.
Logan JM, Jardine TD, Miller TJ, Bunn SE, Cunjak RA, Lutcavage ME. Lipid corrections in carbon and nitrogen stable isotope analyses: comparison of chemical extraction and modelling methods. J Anim Ecol. 2008;77:838–46.
Ruiz-Cooley RI, Garcia KY, Hetherington ED. Effects of lipid removal and preservatives on carbon and nitrogen stable isotope ratios of squid tissues: Implications for ecological studies. J Exp Mar Bio Ecol. 2011;407:101–7.
Zar JH. Biostatistical analysis. 4th ed. New Jersey: Prentice Hall; 1999.
Zuur A, Ieno EN, Smith GM. Analyzing ecological data. 1st ed. New York: Springer-Verlag; 2007.
R Core Team. R: a language and environment for statistical computing. In: R Core Team, (eds).Vienna, Austria, R Foundation for Statistical Computing; 2020.
This work is a contribution of the Distant Squid Fisheries Sci-Tech Group, SHOU. We thank the staff members of the Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai Ocean University for providing assistance at the laboratory. We are grateful to technician Shaoqin Wang for the assistance with fatty acids determination. We are also grateful to Dr. Roshni Subramaniam for his help with editing the revised MS, and to the Editor Ulrich Technau and the anonymous reviewers for their constructive comments and suggestions on the manuscript.
This work was supported by National Natural Science Foundation of China (41876144) and Shanghai Talent Development Funding (2020107) to DL, and National Key Research and Development Project of China (2019YFD0901404), National Natural Science Foundation of China (41876141) and Shanghai Science and Technology Innovation Program (19DZ1207502) to XC.
Ethics approval and consent to participate
Specimens were collected as dead squids from the commercial jigging fisheries landings, during the fishing season from January to March 2019. The specimens were analyzed in the laboratory using methods that are in line with current Chinese national standards, namely Laboratory Animals—General Requirements for Animal Experiment (GB/T 35823-2018). As all material sampled in this work was obtained from commercial jigging fisheries landings and already dead, there was no requirement for ethical approval of sampling protocols as it did not include live organisms.
Consent for participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Fatty acid composition of ovary of Argentinean shortfin squid. Table S2. Linear mixed-effects models results for reproductive energy relation to nitrogen stable isotope ratios, mantle length and chlorophyll-a concentration. Table S3. Linear mixed-effects models results for reproductive energy relation to selected fatty acids, mantle length and chlorophyll-a concentration. Fig. S1. Mantle length and body weight of Argentinean shortfin squid for this work. Fig. S2. Body size distribution of Argentinean shortfin squid with sexual maturation. Fig. S3. Sea surface chlorophyll-a (Chla, mg g-1) concentration around the sampling stations in the Southwest Atlantic Ocean.
About this article
Cite this article
Lin, D., Zang, N., Zhu, K. et al. Energy acquisition strategy for reproduction in a semelparous squid. Front Zool 19, 28 (2022). https://doi.org/10.1186/s12983-022-00473-w
- Energy demand
- Acquisition strategy