INTRODUCTIONTop
Neon flying squid, Ommastrephes bartramii, is widely distributed around the world’s oceans extending from the subtropics to temperate waters in the northern and southern hemispheres except in equatorial waters (Rodhouse 2005Rodhouse P.G. 2005. World squid resources. Review of the state of world marine fishery resources, FAO Tech. Rep. C2:175., Lefkaditou et al. 2011Lefkaditou E., Peristeraki P., Chartosia N., et al. 2011. Recent findings of Ommastrephes bartramii (Cephalopoda: Ommastrephidae) in the eastern Mediterranean and the implication on its range expansion. Mediterr. Mar. Sci. 12(2): 413-428.). Despite its huge potential abundance, it only supports a commercial fishery in the northwest Pacific Ocean. The exploitation of this epipelagic species was started in 1974 (Chen et al. 2009Chen X.J., Liu B.L., Wang Y.G. 2009. Cephalopod of the world. Marine Press, Beijing, China, 1064 pp.), with roughly 200000 to 300000 t landed during the 1980s by Japan, Korea and Taiwan (Burke et al. 1993Burke W.T., Freeberg M., Miles E.L. 1993. United Nations Resolutions on driftnet fishing: An unsustainable precedent for high seas and coastal fisheries management. Ocean Dev. Int. Law. 25: 127-186.). Mainland China started fishing this squid in 1993, and mainly fished in the western waters 160° east of the North Pacific. Recently, commercial squid jigging vessels of mainland China have mainly fished in the areas 150°E to 170°W and 40° to 46°N from May to October with annual catches of between 50000 and 100000 t.
The population structure of O. bartramii has been studied by many researchers for putative seasonal cohorts (Murata 1990Murata M. 1990. Oceanic resources of squids. Mar. Behav. Phys. 18: 19-71., Murata and Hayase 1993Murata M., Hayase S. 1993. Life history and biological information on flying squid (Ommastrephes bartramii) in the North Pacific Ocean. Bull. Int. Nat. North Pacific Comm. 53: l47-182.), rates of infection by helminth parasites (Bower and Margolis 1991Bower S.M., Margolis L. 1991. Potential use of helminth parasites in stock identification of flying squid, Ommastrephes bartramii, in North Pacific waters. Can. J. Zool. 69(4): 1124-1126., Nagasawa et al. 1998Nagasawa K., Mori J., Okamura H. 1998. Parasites as biological tags of stocks of neon flying squid (Ommastrephes bartramii) in the North Pacific Ocean, In: Okutani, T. (1998), Contributed papers to International Symposium on Large Pelagic Squids, July 18-19, 1996, for JAMARC’s 25th anniversary of its foundation, pp. 49-64. ) and mantle length distribution (Murata 1990Murata M. 1990. Oceanic resources of squids. Mar. Behav. Phys. 18: 19-71., Yatsu et al. 1998aYatsu A., Tanaka H., Mori J. 1998a. Population structure of the neon flying squid, Ommastrephes bartramii, in the North Pacific. In: Okutani T. (ed.), Contributed Papers to International Symposium on Large Pelagic Squids. Japan Marine Fishery Resources Research Center, Tokyo, pp. 31-48. ). Murakami et al. (1981)Murakami K., Watanabe Y., Nakata J. 1981. Growth, distribution and migration of flying squid (Ommastrephes bartrami) in the North Pacific. In: Mishima, S. (ed.), Pelagic animals and environments around the Subarctic Boundary in North Pacific. Hokkaido University, Research Institute of North Pacific Fisheries, Hakodate, pp. 161-179 (In Japanese with English abstract) identified the four stocks based on the squid body sizes extra-large (LL), large (L), small (S) and extra-small (SS). Chen et al. (2002)Chen X.J., Tian S.Q., Ye X.C. 2002. Study on population structure of flying squid in Northwestern Pacific based on gray system theory. J. Shanghai Fish. Univ. 11(4): 335-341. (In Chinese with English Abstract) found that two populations existed in the waters of 165°E westward through nine body variables based on the Grey System Theory. A geographic pattern of population genetic variability was observed in O. bartramii, with major genetic differentiation attributable to inconsistency in allele frequency distribution and in levels of genetic variation between the squid from the western and eastern parts of the species, which covers a wide area in the North Pacific Ocean (Katugin 2002Katugin O.N. 2002. Patterns of genetic variability and population structure in the North Pacific squids Ommastrephes bartramii, Todarodes pacificus and Berryteuthis magister. Bull. Mar. Sci. 71(1): 383-420.). Generally, in the North Pacific, the population of O. bartramii mainly comprises two cohorts: a) the autumn cohort hatching from September to February and b) the winter-spring cohort hatching from January to May (Katugin 2002Katugin O.N. 2002. Patterns of genetic variability and population structure in the North Pacific squids Ommastrephes bartramii, Todarodes pacificus and Berryteuthis magister. Bull. Mar. Sci. 71(1): 383-420., Chen and Chiu 2003Chen C.S., Chiu T.S. 2003. Variations of life history parameters in two geographical groups of the neon flying squid, Ommastrephes bartramii, from the North Pacific. Fish. Res. 63: 349-366., Ichii et al. 2004Ichii T., Mahapatra K., Sakai M., et al. 2004. Differing body size between the autumn and the winter–spring cohorts of neon flying squid (Ommastrephes bartramii) related to the oceanographic regime in the North Pacific: a hypothesis. Fish. Oceanogr. 13: 295-309, Bower and Ichii 2005Bower J.R., Ichii T. 2005. The red flying squid (Ommastrephes bartramii): A review of recent research and the fishery in Japan. Fish. Res. 76(1): 39-55., Chen et al. 2011Chen F., Chen X.J., Lu H.J., et al. 2011. Comparison of biological characteristics of Ommastrephes bartarmii between two different areas in the east central waters of North Pacific Ocean. J. Shanghai Fish. Univ. 20(5): 759-764. (In Chinese with English Abstract)). The autumn cohort consists of the central stock and east stock, which separate near 160°W. The winter-spring cohort also comprises the west stock and central-east stock, which separate near 170°E (Ichii et al. 2004Ichii T., Mahapatra K., Sakai M., et al. 2004. Differing body size between the autumn and the winter–spring cohorts of neon flying squid (Ommastrephes bartramii) related to the oceanographic regime in the North Pacific: a hypothesis. Fish. Oceanogr. 13: 295-309, Bower and Ichii 2005Bower J.R., Ichii T. 2005. The red flying squid (Ommastrephes bartramii): A review of recent research and the fishery in Japan. Fish. Res. 76(1): 39-55.). Although these two cohorts overlap geographically, they have been caught in different areas at separate times (Bower and Ichii 2005Bower J.R., Ichii T. 2005. The red flying squid (Ommastrephes bartramii): A review of recent research and the fishery in Japan. Fish. Res. 76(1): 39-55.). The autumn cohort (NE stock) distributed in the waters east of 170°E from May to June and the winter-spring cohort (NW stock) distributed in the waters west of 165°E from July to November have been the main fishing target for Chinese squid jigging fleets since 1998 (Chen and Chiu 2003Chen C.S., Chiu T.S. 2003. Variations of life history parameters in two geographical groups of the neon flying squid, Ommastrephes bartramii, from the North Pacific. Fish. Res. 63: 349-366., Wang and Chen 2005Wang Y.G., Chen X.J. 2005. The Resource and Biology of Economic Oceanic Squid in the World. Ocean Press, Beijing, pp. 79-295. ).
Stock identification is critical for an effective fishery management to avoid overfishing and promote the sustainable development of fisheries (Cadrin and Silva 2005Cadrin S.X., Silva V.M. 2005. Morphometric variation of yellowtail flounder. ICES J. Mar. Sci. 62: 683-694.). Stock structures are often identified and verified on the basis of their different life history strategies and genetic structures. However, morphometric traits are still often used in this field (Sajina et al. 2011Sajina A.M., Chakraborty S.K., Jaiswar A.K., et al. 2011. Stock structure analysis of Megalaspis cordyla (Linnaeus, 1758) along the Indian coast based on truss network analysis. Fish. Res. 108: 100-105). Traditional measurements are based on the conventional orthogonal method, which uses length and width to describe the variables for species with rigid body forms. Unlike fish and many crustaceans, squid have a flexible soft body without a hard surface structure. The cylindrical mantle cavity varies during locomotion and respiration, and the stretched arms and tentacles are also broken easily, often by hooks during the jigging capture process (Cabanellas-Reboredo et al. 2011Cabanellas-Reboredo M., Alós J., Palmer M., et al. 2011. Simulating the indirect handline jigging effects on the European squid (Loligo vulgaris) in captivity. Fish. Res. 110(3): 435-440., Kurosaka et al. 2012Kurosaka K., Yamashita H., Ogawa M., et al. 2012. Tentacle-breakage mechanism for the neon flying squid (Ommastrephes bartramii) during the jigging capture process. Fish. Res. 121: 9-16.). Thus, measurements based on the soft parts of squid are challenging and frequently contain errors, and an alternative structure should be used to separate stocks on the basis of body morphology. However, these measurements are fairly reliable if the approach is correct, as in the studies of Loliginidae by Pierce et al. (1994)Pierce G. J., Hastie L. C., Guerra A., et al. 1994. Morphometric variation in Loligo forbesi and Loligo vulgaris: regional, seasonal, sex, maturity and worker differences. Fish. Res. 21(1): 127-148. and of Ommastrephidae by Martínez et al. (2002)Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141: 131-143..
Hard structures, including the statolith, beak (mandibular) and gladius (pen), which contain a series of ecological information during its mysterious life history, have gradually been used for their stable and constant configuration (Bizikov and Arkhipkin 1997Bizikov V. A., Arkhipkin A. I. 1997. Morphology and microstructure of the gladius and statolith from the boreal Pacific giant squid Moroteuthis robusta (Oegopsida; Onychoteuthidae). J. Zool. 241(3): 475-492., Piatkowski et al. 2001Piatkowski U., Pütz K., Heinemann H. 2001. Cephalopod prey of king penguins (Aptenodytes patagonicus) breeding at Volunteer Beach, Falkland Islands, during austral winter 1996. Fish. Res. 52(1): 79-90., Ikeda et al.2003Ikeda Y., Arai N., Kidokoro H., et al. 2003. Strontium: calcium ratios in statoliths of Japanese common squid Todarodes pacificus (Cephalopoda: Ommastrephidae) as indicators of migratory behavior. Mar. Ecol. Prog. Ser. 251(1): 169-179., Jackson and Domeier 2003Jackson G.D., Domeier M.L. 2003. The effects of an extraordinary El Niño/La Niña event on the size and growth of the squid Loligo opalescens off Southern California. Mar. Biol. 142(5): 925-935., Guerra et al. 2010Guerra A., Rodríguez-Navarro A.B., González A.F., 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(7): 1425-1431., Ruiz-Cooley et al. 2013Ruiz-Cooley R.I., Ballance L.T., McCarthy M.D. 2013. Range Expansion of the Jumbo Squid in the NE Pacific: δ15N Decrypts Multiple Origins, Migration and Habitat Use. PloS One. 8(3): e59651.). As calcified structures embedded in cartilage, a pair of statoliths is an indispensable part of the acceleration receptor system that controls the movement and direction of the cephalopod (Hanlon and Messenger 1996Hanlon R.T., Messenger J.B. 1996. Cephalopod Behaviour. Cambridge University Press Cambridge, UK, 232 pp., Arkhipkin and Bizikov 2000Arkhipkin A.I., Bizikov V.A. 2000. Role of the statolith in functioning of the acceleration receptor system in squids and sepioids. J. Zool. 250(1): 31-55.). They are also used in studies of species identification (Clarke 1978Clarke M.R. 1978. The cephalopod statolith-introduction to its form. J. Mar. Biol. Ass. U.K. 58: 701-712., Dommergues et al. 2000Dommergues J.L., Neige P., Boletzky S.V. 2000. Exploration of morphospace using procrustes analysis in statoliths of cuttlefish and squid (Cephalopoda: Decabrachia)-evolutionary aspects of form disparity. Veliger-berkeley, 43(3): 265-276., Lombarte et al. 2006Lombarte A., Rufino M.M., Sánchez P. 2006. Statolith identification of Mediterranean Octopodidae, Sepiidae, Loliginidae, Ommastrephidae and Enoploteuthidae based on warp analyses. J. Mar. Biol. Ass. U. K. 86(04): 767-771.), age estimation (Villanueva 1992Villanueva R. 1992. Interannual growth differences in the oceanic squid Todarodes angolensis Adam in the northern Benguela upwelling system, based on statolith growth increment analysis. J. Exp. Mar. Biol. Ecol., 159(2): 157-177., Sánchez 1995Sánchez P. 1995. Age and growth of Illex coindetii. ICES J. Mar. Sci. 199: 441-444, Arkhipkin and Shcherbich 2012Arkhipkin A.I., Shcherbich Z.N. 2012. Thirty years’ progress in age determination of squid using statoliths. J. Mar. Biol. Ass. U. K. 92(06): 1389-1398.), growth pattern (Yatsu 2000Yatsu A. 2000. Age estimation of four oceanic squids, Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis, and Illex argentinus (Cephalopoda, Ommastrephidae) based on statolith Microstructure. Jpn. Agri. Res. Quart. 34: 75-80., Yatsu et al. 1997Yatsu A., Midorikawa S., Shimada T., et al. 1997. Age and growth of the neon flying squid, Ommastrephes bartrami, in the North Pacific Ocean. Fish. Res. 29: 257-270) and trace elements (such as strontium) (Durholtz et al. 1997Durholtz M.D., Lipinski M.R., Przybylowicz W.J., et al. 1997. Nuclear microprobe mapping of statoliths of chokka squid Loligo vulgaris reynaudii d’Orbigny, 1845. Biol. Bull. 193(2): 125-140., Ikeda et al. 1996Ikeda Y., Arai N., Sakamoto W., et al. 1996. PIXE analysis of trace elements in squid statoliths: comparison between Ommastrephidae and Loliginidae. International Journal of PIXE. 6: 537-542., 1997Ikeda Y., Arai N., Sakamoto W., et al. 1997. Comparison on trace elements in squid statoliths of different species’ origin as available key for taxonomic and phylogenetic study. International Journal of PIXE. 7: 141-146., Yatsu et al. 1998bYatsu A., Mochioka N., Morishita K., et al. 1998b. Strontium/calcium ratios in statoliths of the neon flying squid, Ommastrephes bartrami (Cephalopoda) in the North Pacific Ocean. Mar. Biol. 131: 275-282.) to determine life history and relevant environmental conditions. As the main feeding organ, the beak can be easily preserved, and is non-corroding, so it has been used in the studies of chemical structure (Miserez et al. 2010Miserez A., Rubin D., Waite J.H. 2010. Cross-linking Chemistry of Squid Beak. J. Biol. Chem. 285(49): 38115-38124.), aging and growth (Yatsu et al. 1997Yatsu A., Midorikawa S., Shimada T., et al. 1997. Age and growth of the neon flying squid, Ommastrephes bartrami, in the North Pacific Ocean. Fish. Res. 29: 257-270, Raya and Hernández-González 1998Raya C.P., Hernández-González C.L. 1998. Growth lines within the beak microstructure of the octopus Octopus vulgaris Cuvier, 1797. S. Afr. J. Mar. Sci. 20(1): 135-142., Yatsu and Mori 2000Yatsu A., Mori J. 2000. Early growth of the autumn cohort of neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean. Fish. Res. 45: 189-194., Raya et al. 2010Raya C.P., Bartolomé A., García-Santamaría M.T., et al. 2010. Age estimation obtained from analysis of octopus (Octopus vulgaris Cuvier, 1797) beaks: Improvement and comparisons. Fish. Res. 106: 171-176, Castanhari and Tomás 2012Castanhari G., Tomás A.R.G. 2012. Beak increment counts as a tool for growth studies of the common octopus Octopus vulgaris in southern Brazil. Bol. Inst. Pesca, São Paulo 38(4): 323-331, Perales-Raya et al. 2014Perales-Raya C., Jurado-Ruzafa A., Bartolomé A., et al. 2014. Age of spent Octopus vulgaris and stress mark analysis using beaks of wild individuals. Hydrobiologia 725(1): 105-114.), species identification (Smale et al. 1993Smale M.J., Clarke M.R., Klages N.T.W., et al. 1993. Octopod beak identification—resolution at a regional level (Cephalopoda, Octopoda: Southern Africa). S. Afr. J. Mar. Sci. 13(1): 269-293.), biomass estimation (Lu and Ickeringill 2002Lu C.C., Ickeringill R. 2002. Cephalopod beak identification and biomass estimation techniques: tools for dietary studies of southern Australian finfishes. Mus. Victoria Sci. Rep. 6: 1-65.), trophic dynamics (Cherel and Hobson 2005Cherel Y., Hobson K.A. 2005. Stable isotopes, beaks and predators: a new tool to study the trophic ecology of cephalopods, including giant and colossal squids. Proc. R. Soc. Lond. B. 272: 1601-1607., Cherel et al. 2009Cherel Y., Ridoux V., Spitz J., et al. 2009. Stable isotopes document the trophic structure of a deep-sea cephalopod assemblage including giant octopod and giant squid. Biol. Lett. 5: 364-367.) and paralarval ontogeny (Uchikawa et al. 2009Uchikawa K., Sakai M., Wakabayashi T., et al. 2009. The relationship between paralarval feeding and morphological changes in the proboscis and beaks of the neon flying squid Ommastrephes bartramii. Fish. Sci. 75(2): 317-323.).
The measurement of hard structures has been the basis of these studies. Radial measurement is a simple way to analyse the relationship between statolith shape and growth pattern (Arkhipkin 2003Arkhipkin A.I. 2003. Towards identification of the ecological lifestyle in nektonic squid using statolith morphometry. J. Mollus. Stud. 69: 171-178., Ma et al. 2009Ma J., Chen X.J., Liu B.L., et al. 2009. Morphologic Features of Statolith for Ommastrephes bartramii in the Northwest Pacific Ocean. Periodical of Ocean University of China 39(2): 215-220. (In Chinese with English abstract), Chen et al. 2010Chen X.J., Ma J., Liu B.L., et al. 2010. Effects of sexual maturity and body size on statolith shape of Ommastrephes bartramii in the Northwest Pacific Ocean. J. Fish. China. 30(6): 928-934. (In Chinese with English abstract)). The beak also has specific characteristics and has been used in species identification (Clarke 1986Clarke M.R. 1986. A handbook for the identification of cephalopod beaks. Clarendon Press, Oxford, UK, 273 pp., Ogden et al. 1998Ogden R.S., Allcock A.L., Wats P.C., et al. 1998. The role of beak shape in octopodid taxonomy. S. Afr. J. Mar. Sci. 20(1): 29-36., Lu and Ickeringill 2002Lu C.C., Ickeringill R. 2002. Cephalopod beak identification and biomass estimation techniques: tools for dietary studies of southern Australian finfishes. Mus. Victoria Sci. Rep. 6: 1-65., Xavier and Cherel 2009Xavier J.C., Cherel Y. 2009. Cephalopod beak guide for the southern ocean. British Antarctic Survey Press, Cambridge, UK, 129 pp.). It has been proved that hard structures with a stable form can perform better than soft tissues (mantle, arm, tentacle, etc.) in squid population division (Martínez et al. 2002Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141: 131-143.). Some new methods (landmarks and outline study) have also revealed that hard structures give good results in identity analysis (Dommergues 2000Dommergues J.L., Neige P., Boletzky S.V. 2000. Exploration of morphospace using procrustes analysis in statoliths of cuttlefish and squid (Cephalopoda: Decabrachia)-evolutionary aspects of form disparity. Veliger-berkeley, 43(3): 265-276., Neige and Dommergues 2002Neige P., Dommergues, J.L. 2002. Disparity of beaks and statoliths of some coleoids a morphometric approach to depict shape differentiation. Gabh. der Geol. Bun. 57: 393-399, Lombarte et al. 2006Lombarte A., Rufino M.M., Sánchez P. 2006. Statolith identification of Mediterranean Octopodidae, Sepiidae, Loliginidae, Ommastrephidae and Enoploteuthidae based on warp analyses. J. Mar. Biol. Ass. U. K. 86(04): 767-771., Neige 2006Neige P. 2006. Morphometrics of hard structures in cuttlefish. Vie Milieu 56(2): 121-127., Crespi-Abril et al. 2010Crespi-Abril A.C., Morsan E.M., Barón P.J. 2010. Analysis of the ontogenetic variation in body and beak shape of the Illex argentinus inner shelf spawning groups by geometric morphometrics. J. Mar. Biol. Ass. U. K. 90(03): 547-553.). Some researchers may be cautious about using the beak as an identification material (Xavier et al. 2007Xavier J.C., Clarke M.R., Magalhães M.C., et al. 2007. Current status of using beaks to identify cephalopods: III International Workshop and training course on Cephalopod beaks, Faial Island, Azores. Arquipélago 24: 41-48.), but it is still important for stock identification (Xavier et al. 2011Xavier J.C., Phillips R.A., Cherel Y. 2011. Cephalopods in marine predator diet assessments: why identifying upper and lower beaks is important. ICES J. Mar. Sci. 68(9): 1857-1864.). Therefore, the phenotypic characteristics, such as statoliths and beaks, can play a supportive role for stock identification and investigation.
A morphological variation of cephalopods has occurred among the species due to different genetic structures and populations, induced mainly by oceanographic environmental factors (Adkison 1995Adkison M.D. 1995. Population differentiation in Pacific salmons: local adaptation genetic drift, or the environment? Can. J. Fish. Aquat. Sci. 52(12): 2762-2777.). Multivariate analyses of morphometric data have been used in some cephalopods at intraspecies (Martínez et al. 2002Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141: 131-143.) and interspecies (Wolff 1984Wolff G.A. 1984. Identification and estimation of size from the beaks of 18 species of cephalopods from the Pacific Ocean. NOAA Technical Report NMFS 17: 1-50) levels for the evaluation of taxonomic and geographic variations. Stocks inhabiting different environments can impact their growth patterns, including the shapes of statolith and beak. It is common that the separation of stocks is only based on one material and then compared with the results obtained from different materials (Martínez et al. 2002Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141: 131-143.), and the approach of integrating two materials in the stock separation is rarely used. In this study we integrated the morphological variables of statolith and beaks to compare their differences between stocks and between the sexes and to establish a discriminant function for two seasonal cohorts in the North Pacific Ocean. Our aims were to identify different cohorts of O. bartramii by hard structures (statolith and beaks). This integrated approach can also be used for other species when different body structures are available in the stock identifications.
MATERIALS AND METHODSTop
The survey was carried out in the waters of 150°E-177°W and 38°-44°N from May to October in 2010 to 2012 by the Chinese squid jigging vessel Jinhai 827. The samples were collected randomly from the daily catches and frozen on board immediately for future analysis. The sampling station is shown in Figure 1.
|
||
|
A total of 570 individuals were collected. The mantle length (ML) was measured to the nearest 0.1 mm after thaw in the laboratory, and the sex was identified by their entirely different gonad structure. The sexual maturity stages were determined according to Lipinski and Underhill (1995)Lipinski M. R., Underhill L. G. 1995. Sexual maturation in squid: quantum or continuum. S. Afr. J. Mar. Sci. 15: 207-223. Statoliths and beaks were dissected according to Raya et al. (2010)Raya C.P., Bartolomé A., García-Santamaría M.T., et al. 2010. Age estimation obtained from analysis of octopus (Octopus vulgaris Cuvier, 1797) beaks: Improvement and comparisons. Fish. Res. 106: 171-176 and Chen et al. (2013)Chen X.J., Li J.H., Liu B.L., et al. 2013. Age, growth and population structure of Jumbo flying squid, Dosidicus gigas, off the Costa Rica Dome. J. Mar. Biol. Ass. U. K. 93:567-573. As a result, a subsample of 406 pairs of statoliths and beaks were prepared for the analysis. A pair of beaks were precisely paired with a pair of statoliths which had been picked from the same squid sample (Table 1).
Table 1. – Sample information of O. bartramii for the two stocks.
Stock structures | Sampling date | Latitude | Longitude | Number of samples | Sex (F,M) | ML (mm) |
---|---|---|---|---|---|---|
Northeast Pacific Stock | 11-21 Jun. 2010 | 39°48’N-40°09’N | 171°52’E-175°29’W | 71 | 11,60 | 212-375 |
May. to Jun. 2011 | 38°42’N-39°20’N | 172°11’E-177°30’W | 23 | 21,2 | 226-411 | |
May. to Jul. 2011 | 39°02’N-40°21’N | 174°52’E-179°58’W | 91 | 86,5 | 219-483 | |
Northwest Pacific Stock | Jul. to Oct. 2011 | 38°42’N-39°20’N | 151°23’E-159°25’E | 173 | 100,73 | 173-452 |
Aug to Nov. 2011 | 40°58’N-43°21’N | 150°21’E-156°08’E | 47 | 34,13 | 208-363 |
Photos were taken using a charge-coupled device (CCD, connecting device) for right statolith under a 50× optical microscope (Olympus) and nine parameters of statolith morphology, i.e. total statolith length (TSL), maximum width (MW), dorsal dome length (DDL), lateral dome length (LDL), dorsal lateral length (DLL), rostrum lateral length (RLL), rostrum length (RL), rostrum width (RW) and wing length (WL), were measured by using the image analysis software WT-Tiger3000 (Fig. 2A). The parameters of each statolith were measured to the nearest 0.01µm. Beaks were also measured by vernier caliper, which included upper hood length (UHL), upper crest length (UCL), upper rostrum length (URL), upper rostrum width (URW), upper lateral wall length (ULWL), upper wing length (UWL), lower hood length (LHL), lower crest length (LCL), lower rostrum length (LRL), lower rostrum width (LRW), lower lateral wall length (LLWL), and lower wing length (LWL) (Fig. 2B). All the data for beaks were measured to the nearest 0.01 mm. Measures of these variables were obtained independently by the two readers. The average was used if the range of the counts for the same sample was within 5% of the error; otherwise it was measured again and the average of the variables of each reader was used for the same sample (Francis and Mattlin 1986Francis R.I.C.C., Mattlin R.H. 1986. A possible pitfall in the morphometric application of discriminant analysis: measurement bias. Mar. Biol. 93(2): 311-313., Chen et al. 2013Chen X.J., Li J.H., Liu B.L., et al. 2013. Age, growth and population structure of Jumbo flying squid, Dosidicus gigas, off the Costa Rica Dome. J. Mar. Biol. Ass. U. K. 93:567-573).
|
||
|
A Student t-test was conducted to compare differences between the geographic stock and sexes. All of the variables were subjected to normal distribution (Kolmogorov-Smirnov test, P>0.05). Considering the impact of allometric growth (Moltschaniwskyj 1995Moltschaniwskyj N.A. 1995. Changes in shape associated with growth in the loliginid squid Photololigo sp.: a morphometric approach. Can. J. Zool. 73(7): 1335-1343., Lombarte et al. 1997Lombarte A., Sanchez P., Morales-Nin B. 1997. Intraspecific shape variability in statoliths of three cephalopod species. Vie Milieu 47: 165-169., O’Dor and Hoar 2000O’Dor R.K., Hoar J.A. 2000. Does geometry limit squid growth? ICES J. Mar. Sci. 57(1): 8-14), raw data standardization should be done before the analysis. The normalization method, as introduced by Lleonart (2000)Lleonart J., Salat J., Torres G.J. 2000. Removing allometric effects of body size in morphological analysis. J. Thero. Biol. 205: 85-93., was used to standardize morphological variables of statoliths and beaks. The accuracy of this method has been demonstrated in related investigations (Pineda et al. 2002Pineda S.E., Hernández D.R., Brunetti N.E., et al. 2002. Morphological identification of two Southwest Atlantic Loliginid squids: Loligo gahi and Loligo sanpaulensis. Rev. Invest. Desarr. Pesq. 15: 67-84., Vega et al. 2002Vega M.A., Rocha F.J., Guerra A., et al. 2002. Morphological difference between the Patagonian squid Loligo gahi populations from the Pacific and Atlantic Oceans. Bull. Mar. Sci. 71(2): 903-913, Lefkaditou and Bekas 2004Lefkaditou E., Bekas P. 2004. Analysis of beak morphometry for the horned octopus Eledone cirrhosa (Cephalopoda: Octopoda) from the Thracian Sea (NE Mediterranean). Mediterr. Mar. Sci. 5(1): 143-149., Chen et al. 2012Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481.). The TSL in statolith and UHL in beak were chosen as the independent variables to standardize the other variables of beaks (Chen et al. 2012Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481.). The standardized morphometric variables were represented by adding a lower case letter “s” after each variable, i.e. MWs, DDLs, LDLs, DLLs, RLLs, RLs, RWs, WLs; or UCLs, URLs, ULWLs, UWLs, LHLs, LCLs, LRLs, LLWLs and LWLs.
A stepwise discriminant analysis (SDA) was performed to select the significant standardized morphological variables based on the statolith, beak and the combination of the two structures (P<0.05; Rencher 2002Rencher A.C. 2002. Methods of Multivariate Analysis, 2nd edition. John Wiley & Sons, Inc, New York.), and the classification functions were developed for the three different materials (i.e., statolith, beak, and their combinations). Finally, a leave-one-out cross-validation (the Jackknife method) was used to determine rates of successful classification of squid from the two stocks for different uses of the hard structures.
RESULTSTop
Sexual dimorphism and variation of different cohorts in hard structure sizes
The morphometric variables of beaks and statoliths in both of the two cohorts are shown in Table 2 and Table 3. Significant differences were found in the beak variables between sexes for the NE cohort (P<0.05). The NW cohort also shows significant differences for all variables (P<0.05) except for UWL (P>0.05). Apparently, sexual dimorphism of statolith variables is similar to that of the beak for the NE cohort, but t-tests showed no significant differences between males and females in the NW cohort (P>0.05). The above results indicated that between-sex differences were greater for the NE stock than for the NW stock (Tables 2 and 3).
Table 2. – Beak morphological variables and P values (t-tests) of O. bartramii in the North Pacific Ocean; ***, P significant at α=0.05; ns, nonsignificant.
Variable | NE (mean±sth, mm) | NW (mean±sth, mm) | P (females between stocks) | P (males between stocks) | ||||
---|---|---|---|---|---|---|---|---|
Females | Males | P (between sexes) | Females | Males | P (between sexes) | |||
UHL | 24.76±5.15 | 17.16±0.86 | *** | 19.01±4.05 | 17.80±2.56 | *** | *** | *** |
UCL | 30.42±6.24 | 21.11±1.09 | *** | 23.13±4.96 | 21.45±3.22 | *** | *** | ns |
URL | 8.00±1.76 | 5.87±1.55 | *** | 6.27±1.35 | 5.94±0.94 | *** | *** | ns |
URW | 7.02±1.41 | 4.78±0.48 | *** | 5.14±1.22 | 4.68±0.80 | *** | *** | ns |
ULWL | 26.52±5.35 | 18.15±1.43 | *** | 20.02±4.29 | 18.53±2.85 | *** | *** | ns |
UWL | 8.23±1.70 | 6.40±0.92 | *** | 6.09±1.44 | 5.99±0.99 | ns | *** | *** |
LHL | 7.97±1.64 | 5.91±0.50 | *** | 6.33±1.30 | 5.87±0.71 | *** | *** | ns |
LCL | 15.43±3.37 | 10.90±0.85 | *** | 12.35±3.17 | 10.99±1.84 | *** | *** | ns |
LRL | 7.23±1.46 | 5.34±0.85 | *** | 5.46±1.28 | 5.08±0.89 | *** | *** | ns |
LRW | 7.31±1.50 | 5.28±0.79 | *** | 5.40±1.25 | 5.00±0.71 | *** | *** | *** |
LLWL | 22.83±4.59 | 15.42±0.95 | *** | 16.98±3.81 | 15.54±2.74 | *** | *** | ns |
LWL | 13.13±2.71 | 9.15±0.66 | *** | 9.84±2.33 | 9.22±1.46 | *** | *** | ns |
Table 3. – Statolith morphological variables and P values (t-tests) of O. bartramii in the North Pacific Ocean.***, P significant at α=0.05; ns, nonsignificant.
Variable | NE(mean±sth, μm) | NW(mean±sth, μm) | P (females between stocks) | P (males between stocks) | ||||
---|---|---|---|---|---|---|---|---|
Females | Males | P (between sexes) | Females | Males | P (between sexes) | |||
TSL | 1446.44±126.76 | 1272.02±58.51 | *** | 1260.71±107.13 | 1248.79±76.41 | ns | *** | *** |
MW | 844.75±96.23 | 745.62±39.20 | *** | 724.92±82.89 | 717.53±76.07 | ns | *** | *** |
DDL | 601.04±94.41 | 548.32±70.18 | *** | 550.14±86.05 | 540.142±92.63 | ns | *** | ns |
DLL | 651.57±115.43 | 515.99±72.88 | *** | 575.90±117.82 | 555.96±121.70 | ns | *** | *** |
LDL | 866.25±105.53 | 740.12±63.13 | *** | 719.50±75.55 | 711.53±65.48 | ns | *** | *** |
RLL | 872.40±114.00 | 774.33±73.05 | *** | 802.65±109.54 | 816.26±99.29 | ns | *** | *** |
RL | 447.53±57.71 | 386.66±41.72 | *** | 418.36±53.86 | 414.37±46.41 | ns | *** | *** |
RW | 232.89±33.47 | 212.50±27.05 | *** | 187.87±34.19 | 187.10±32.44 | ns | *** | *** |
WL | 1123.73±112.01 | 1017.86±59.64 | *** | 1008.41±103.69 | 1000.45±74.22 | ns | *** | ns |
For a given sex, different stocks also revealed different patterns with different hard structures. Female individuals had significant differences in beak morphology between stocks (t-test, P<0.01), but had no significant differences in male squids (P>0.05), with the exception of UHL, UWL and LRW between stocks (P<0.05). The morphometric characteristics of statoliths showed significant difference between NE and NW cohorts in beak shape for female squids (P<0.05). DDL and WL showed no differences in male squids (P>0.05), but the other variables of statoliths were significantly different between stocks (P<0.05) (Table 3).
Discriminant analysis using standardized beak variables
We divided the squid samples into four groups based on geography and sex. The beak was chosen as the only analytical material first. Stepwise discriminant analyses showed that the six variables, URWS, UWLS, LCLS, ULWLS, LRWS and URLS, could explain the morphological features among the four groups. The Wilks λ was decreased from 0.551 to 0.419 (Table 4). Three canonical functions effectively separated these four groups, explaining 82.1%, 15.6% and 2.3% of the total variance. The distribution of the four groups on function 1 overlapped except for western females, which could be more easily identified (Fig. 3). The cross-validation rate was 36.0% for the western males (NW-M), 43.6% for the western females (NW-F), 68.7% for the eastern males (NE-M) and 74.6% for the eastern females (NE-F), showing slightly lower values than the original ones (Table 4). The classification functions with coefficients are presented in Table 5.
Table 4. – The results estimated fromstepwise discriminant analysis of standardized beak variables for each sex of two populations, and a classification matrix with percentages of correctly classified individuals and cross-validation results.
Step | Variable | F to enter | Wilks λ | df 1 | df 2 | P value of Wilks λ |
---|---|---|---|---|---|---|
1 | URWS | 108.583 | 0.551 | 3 | 400 | <0.001 |
2 | UWLS | 51.138 | 0.522 | 6 | 798 | <0.001 |
3 | LCLS | 39.256 | 0.469 | 9 | 969 | <0.001 |
4 | ULWLS | 30.997 | 0.448 | 12 | 1051 | <0.001 |
5 | LRWS | 25.828 | 0.433 | 15 | 1094 | <0.001 |
6 | URLS | 22.366 | 0.419 | 18 | 1118 | <0.001 |
Group | Number of specimens classified | Original (%) | Cross-validation (%) | |||
NW-M | NW-F | NE-M | NE-F | |||
NW-M | 33 | 29 | 20 | 4 | 38.4 | 36.0 |
NW-F | 32 | 59 | 21 | 21 | 44.4 | 43.6 |
NE-M | 8 | 12 | 47 | 0 | 70.1 | 68.7 |
NE-F | 6 | 10 | 14 | 88 | 74.6 | 74.6 |
Total | 86 | 133 | 67 | 118 | 56.9 | 55.7 |
|
||
|
Table 5. – Coefficients of parameters in classification functions using beaks as a material.
Parameters | Group | |||
---|---|---|---|---|
NW-M | NW-F | NE-M | NE-F | |
URLS | –8.442 | –9.151 | –12.161 | –13.027 |
URWS | –62.884 | –61.306 | –59.247 | –53.328 |
ULWLS | 217.729 | 216.815 | 209.660 | 220.197 |
UWLS | –16.961 | –20.450 | –12.756 | –17.187 |
LCLS | –21.760 | –17.860 | –22.924 | –28.883 |
LRWS | –72.472 | –72.042 | –66.745 | –66.229 |
Constant | –207.284 | –211.796 | –198.656 | –232.122 |
Discriminant analysis using standardized statolith variables
When referring to statoliths, RLS, LDLS, RLLS, RWS and WLS were chosen as the most important variables from ten morphological parameters for the population discrimination (Table 6). The total Wilks λ was 2.626 for the five variables. Canonical functions 1 and 2, explaining 80.9% and 18.5%, respectively, together explained almost 100% of all the groups (99.4%). The distribution of the NE-M revealed little overlap with the other three groups on function 1 (Fig. 4). Therefore, the highest cross-validation rate was 86.6% for the NE-M, and the remainder was 36.0% for the NW-M, 37.6% for the NW-F and 50.8% for the NE-F. The classification result was similar to that of the beak (Table 6). The classification functions with coefficients are presented in Table 7.
Table 6. – The results estimated from stepwise discriminant analysis of standardized statolith variables for each sex of two populations, and a classification matrix with percentages of correctly classified individuals and cross-validation results.
Step | Variable | F to enter | Wilks λ | df 1 | df 2 | P value of Wilks λ |
---|---|---|---|---|---|---|
1 | RLS | 87.087 | 0.605 | 3 | 400 | <0.001 |
2 | LDLS | 45.585 | 0.555 | 6 | 798 | <0.001 |
3 | RLLS | 34.740 | 0.506 | 9 | 969 | <0.001 |
4 | RWS | 27.079 | 0.490 | 12 | 1051 | <0.001 |
5 | WLS | 22.910 | 0.470 | 15 | 1094 | <0.001 |
Group | Number of specimens classified | Original (%) | Cross-validation (%) | |||
NW-M | NW-F | NE-M | NE-F | |||
NW-M | 35 | 31 | 3 | 17 | 40.7 | 36.0 |
NW-F | 42 | 55 | 7 | 29 | 41.4 | 37.6 |
NE-M | 2 | 0 | 58 | 7 | 86.6 | 86.6 |
NE-F | 13 | 18 | 26 | 61 | 51.7 | 50.8 |
Total | 86 | 133 | 67 | 118 | 55.1 | 52.8 |
|
||
|
Table 7. – Coefficients of parameters in classification functions using statoliths as a material.
Parameters | Group | |||
---|---|---|---|---|
NW-M | NW-F | NE-M | NE-F | |
LDLS | 53.722 | 55.337 | 55.841 | 63.644 |
RLLS | 118.483 | 116.221 | 111.051 | 110.942 |
WLS | 372.294 | 372.568 | 367.833 | 365.621 |
RLS | 114.364 | 115.273 | 105.825 | 112.025 |
RWS | 36.761 | 36.775 | 39.616 | 40.201 |
Constant | –4572.500 | –4579.0 31 | –4371.268 | –4524.986 |
Discriminant analysis with combination of standardized beak and statolith variables
When combined with the two hard structures, the SDA results showed that eight variables (URWS, MWS, RLLS, UWLS, LLWLS, DLLS, LCLS and RLS) were identified effectively among all the four groups within geography and sex. Total Wilks λ was 1.630, decreasing sharply from 0.551 to 0.127 (Table 8). Canonical function 1 showed the highest rate of 86.5%, followed by 12.3% and 1.2% for canonical functions 2 and 3, respectively (Table 8). The different geographical stocks were effectively separated by function 1, and there was little overlap between sexes in the NW stock (Fig. 5). Overall, the successful classification rate was 71.7% with 53.5% for the NW-M, with the remainder for the NW-F at 56.4%, 95.5% for the NE-M and 81.4% for the NE-F respectively, over 15% higher than in the classification with either hard structure alone. The classification functions with coefficients are presented in Table 9.
Table 8. – The results estimated from stepwise discriminant analysis of combined standardized statolith and beak variables for each sex of two populations, and a classification matrix with percentages of correctly classified individuals and cross-validation results.
Step | Variable | F to enter | Wilks λ | df 1 | df 2 | P value of Wilks λ |
---|---|---|---|---|---|---|
1 | URWS | 108.583 | 0.551 | 3 | 400 | <0.001 |
2 | MWS | 158.336 | 0.208 | 6 | 798 | <0.001 |
3 | RLLS | 114.037 | 0.172 | 9 | 969 | <0.001 |
4 | UWLS | 88.776 | 0.157 | 12 | 1051 | <0.001 |
5 | LLWLS | 73.143 | 0.147 | 15 | 1094 | <0.001 |
6 | LDLS | 63.396 | 0.137 | 18 | 1118 | <0.001 |
7 | LCLS | 55.599 | 0.131 | 21 | 1132 | <0.001 |
8 | RLS | 49.309 | 0.127 | 24 | 1140 | <0.001 |
Group | Number of specimens classified | Original (%) | Cross-validation (%) | |||
NW-M | NW-F | NE-M | NE-F | |||
NW-M | 50 | 33 | 3 | 0 | 58.1 | 53.5 |
NW-F | 46 | 78 | 7 | 2 | 58.6 | 56.4 |
NE-M | 0 | 1 | 65 | 1 | 97.0 | 95.5 |
NE-F | 0 | 1 | 21 | 96 | 81.4 | 81.4 |
Total | 86 | 133 | 67 | 118 | 73.8 | 71.7 |
|
||
|
Table 9. – Coefficients of parameters in classification functions using beaks and statoliths as materials.
Parameters | Group | |||
---|---|---|---|---|
NW-M | NW-F | NE-M | NE-F | |
MWS | 1250.352 | 1246.126 | 1227.866 | 1218.402 |
DLLS | 150.988 | 150.077 | 143.442 | 144.485 |
RLLS | 508.075 | 505.310 | 492.537 | 493.931 |
RLS | –3.536 | –2.393 | –7.247 | –1.873 |
URWS | –169.509 | –168.181 | –164.857 | –159.222 |
UWLS | 28.458 | 24.553 | 29.583 | 26.374 |
LCLS | –216.664 | –212.422 | –212.382 | –219.119 |
LLWLS | –282.051 | –280.431 | –274.622 | –263.275 |
Constant | –10956.687 | –10885.578 | –10404.056 | –10417.180 |
DISCUSSIONTop
Since the 1960s, the structures and morphological characteristics of the statolith and beak of cephalopods have attracted much interest (Clarke 1962Clarke M.R. 1962. The identification of cephalopod “beaks” and the relationship between beak size and total body weight. Bull. Br. Mus. nat. Hist. Zool. 8: 419-480, 1978Clarke M.R. 1978. The cephalopod statolith-introduction to its form. J. Mar. Biol. Ass. U.K. 58: 701-712., 2003Clarke M.R. 2003. Potential of statoliths for interpreting coleoid evolution: a brief review. Ber. Palä. Abh. 3: 37-47., Villanueva 1992Villanueva R. 1992. Interannual growth differences in the oceanic squid Todarodes angolensis Adam in the northern Benguela upwelling system, based on statolith growth increment analysis. J. Exp. Mar. Biol. Ecol., 159(2): 157-177., Arkhipkin 2003Arkhipkin A.I. 2003. Towards identification of the ecological lifestyle in nektonic squid using statolith morphometry. J. Mollus. Stud. 69: 171-178., Chen et al. 2012Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481.). For this reason we chose the statolith and beak as materials in this study for their rigid characteristics and wide usage. Meanwhile, data standardization can effectively remove the influence of allometric growth. There are also other sorts of data standardization, such as log-transformation. Some studies have already proved that the method of data standardization used in this study is more effective than the original data, and there are other methods for identification with a low rate of misclassification (Chen et al. 2012Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481., Fang et al. 2012Fang Z., Chen X.J., Lu H.J., et al. 2012. Morphological differences in statolith and beak between two spawning stocks for Illex argentinus. Acta. Ecol. Sin. 32(19): 5986-5997. (In Chinese with English Abstract))
Sexual dimorphism usually occurs in cephalopods (Mercer et al. 1980Mercer M.C., Misra R.K., Hurley G.V. 1980. Sex determination of the Ommastrephid squid Illex illecebrosus using beak morphometric. Can. J. Fish. Aquat. Sci. 37: 283-286, Pineda et al. 2002Pineda S.E., Hernández D.R., Brunetti N.E., et al. 2002. Morphological identification of two Southwest Atlantic Loliginid squids: Loligo gahi and Loligo sanpaulensis. Rev. Invest. Desarr. Pesq. 15: 67-84., Vega et al. 2002Vega M.A., Rocha F.J., Guerra A., et al. 2002. Morphological difference between the Patagonian squid Loligo gahi populations from the Pacific and Atlantic Oceans. Bull. Mar. Sci. 71(2): 903-913, Chen et al. 2012Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481.). Our study also found this difference in the NE stock in both the statolith and beak (P<0.05), and this was mainly affected by the sex-segregated migration whereby males are separated from northward migrating females, which stay at the spawning/nursery ground to avoid cannibalism (Yatsu 1997Yatsu A., Midorikawa S., Shimada T., et al. 1997. Age and growth of the neon flying squid, Ommastrephes bartrami, in the North Pacific Ocean. Fish. Res. 29: 257-270, O’Dor and Dawe 1998O’Dor R.K., Dawe E.G. 1998. Illex illecebrosus. In: Rodhouse P.G., Dawe E.G., O’Dor R.K. (eds), Squid recruitment dynamics. FAO Fish Tech Pap 376, Rome, pp. 77-104., Ichii et al. 2009Ichii T., Mahapatra K., Sakai M., et al. 2009. Life history of the neon flying squid: effect of the oceanographic regime in the North Pacific Ocean. Mar. Ecol. Prog. Ser. 378: 1-11.). This meant that we could not find more male individuals in our study areas. With the disparity of latitude distribution, environmental conditions can eventually affect the morphology of hard structures. This was also reflected in the results from the discriminant analysis (Tables 2 and 3).
The NW stock revealed a different situation in morphometric characteristics. Most beak characteristics varied according to sex (P<0.05); however, none of the statolith variables showed differences between the sexes (P>0.05). Both sexes of the NW stock had lived in the same area, and shared a similar oceanographic environment with a single migration pattern. Therefore, the difference in the beaks might have been caused by asynchronous maturation, which is a type of reproductive strategy commonly found in other cephalopods (Rocha et al. 2001Rocha F., Guerra A., González A.F. 2001. A review of reproductive strategies in cephalopods. Biol. Rev. 76: 291-304.). However, this sexual dimorphism was small, as is shown by the overlapped dots in the discriminant analysis (Figs 3 and 4).
The population structure of O. bartramii in the Pacific Ocean has been discussed in previous studies (Murakami et al. 1981Murakami K., Watanabe Y., Nakata J. 1981. Growth, distribution and migration of flying squid (Ommastrephes bartrami) in the North Pacific. In: Mishima, S. (ed.), Pelagic animals and environments around the Subarctic Boundary in North Pacific. Hokkaido University, Research Institute of North Pacific Fisheries, Hakodate, pp. 161-179 (In Japanese with English abstract), Murata 1990Murata M. 1990. Oceanic resources of squids. Mar. Behav. Phys. 18: 19-71., Bower and Margolis 1991Bower S.M., Margolis L. 1991. Potential use of helminth parasites in stock identification of flying squid, Ommastrephes bartramii, in North Pacific waters. Can. J. Zool. 69(4): 1124-1126.), and the two main seasonal cohorts (the NE stock with large size located east of 170°E and the NW stock with small size located west of 170°E, Chen and Chiu 2003Chen C.S., Chiu T.S. 2003. Variations of life history parameters in two geographical groups of the neon flying squid, Ommastrephes bartramii, from the North Pacific. Fish. Res. 63: 349-366.) lived in different habitats. One of the reasons for the difference in statolith (both sexes) and beak (female) morphology is a separated migration trajectory. The NW stock occurs in the subtropical frontal zone (STFZ) with rich productivity, whereas NE stock occurs in the subtropical domain which is less productive because it is far from the transition zone chlorophyll front (TZCF). When the TZCF shifts northward in spring, the migration pattern of the NW stock moves with the change of the TZCF, but the NE stock remains southward to the north of the TZCF until summer or autumn (Ichii et al. 2009Ichii T., Mahapatra K., Sakai M., et al. 2009. Life history of the neon flying squid: effect of the oceanographic regime in the North Pacific Ocean. Mar. Ecol. Prog. Ser. 378: 1-11.). Therefore, the NW stock is of large size in the enhanced productivity water, feeding on myctophid (Symbolophorus californiensis, Ceratoscopelus warmingi) and squid (Onychoteuthis borealijaponica, Gonatus berryi). The small-sized NE stock feed on euphausiids, amphipods and fish (Maurolicus imperatorius) (Ichii et al. 2004Ichii T., Mahapatra K., Sakai M., et al. 2004. Differing body size between the autumn and the winter–spring cohorts of neon flying squid (Ommastrephes bartramii) related to the oceanographic regime in the North Pacific: a hypothesis. Fish. Oceanogr. 13: 295-309). This disparity of feeding may induce variation in the beak morphology (Uchikawa et al. 2009Uchikawa K., Sakai M., Wakabayashi T., et al. 2009. The relationship between paralarval feeding and morphological changes in the proboscis and beaks of the neon flying squid Ommastrephes bartramii. Fish. Sci. 75(2): 317-323.). The high rate of classification also demonstrated the difference between the two stocks (Fig. 5).
Many studies have focused on interspecies identification, with relatively high classification rates using beaks (Martínez et al. 2002Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141: 131-143., Pineda et al. 2002Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141: 131-143., Chen et al. 2012Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481.). The total correct classification rate was nearly 50% for both statoliths and beaks, but increased to about 70% for the combination of both hard parts. Thus, adding more relevant variables can increase the success of classification, using SDA to improve the accuracy of classification by retaining some highly correlated variables although the results of correct classification are still low compared with other studies of cephalopods (Martínez et al. 2002Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: Ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141: 131-143., Chen et al. 2012Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481.). The result was credible as the first three variables entered the combination group, which all entered the beak-based and statolith-based SDA (Table 8).
There are still two sub-stocks in the northwest Pacific (central stock and east stock) and northeast Pacific (west stock and central-east stock) (Bower and Ichii 2005Bower J.R., Ichii T. 2005. The red flying squid (Ommastrephes bartramii): A review of recent research and the fishery in Japan. Fish. Res. 76(1): 39-55.). The relatively high successful classification rates demonstrate that morphology of hard parts in sub-stocks showed less variation and it is hard to separate them using only morphology of hard structures.
In conclusion, hard structures, such as statolith and beak in cephalopod, are credible materials to identify at an intraspecies level. The NE stock of O. bartramii showed sexual dimorphism in hard parts, but the NW stock showed no significant difference by sex. The among-stock differences were greater in the female group than in the male group. The classification approach could separate the two stocks, especially the NE stock, with more than 80% accuracy. Meanwhile, data standardization and variable combination could help improve the accuracy of the classification. As the overlapped dots showed in the results, we must be cautious of over-interpretation of phenotype-based analysis, which is more reliable when combined with genetic methods.
ACKNOWLEDGEMENTSTop
This work was funded by the National Science Foundation of China (NSFC41306127 and NSFC41276156), the National Science Foundation of Shanghai (13ZR1419700), the Innovation Programme of Shanghai Municipal Education Commission (13YZ091), the State 863 projects (2012AA092303), the Funding Program for Outstanding Dissertations in Shanghai Ocean University, the Funding Scheme for Training Young Teachers in Shanghai Colleges and the Shanghai Leading Academic Discipline Project (Fisheries Discipline). Y Chen’s involvement was supported by SHOU International Center for Marine Studies and Shanghai 1000 Talent Program.
REFERENCESTop
Adkison M.D. 1995. Population differentiation in Pacific salmons: local adaptation genetic drift, or the environment? Can. J. Fish. Aquat. Sci. 52(12): 2762-2777.
http://dx.doi.org/10.1139/f95-865
Arkhipkin A.I. 2003. Towards identification of the ecological lifestyle in nektonic squid using statolith morphometry. J. Mollus. Stud. 69: 171-178.
http://dx.doi.org/10.1093/mollus/69.3.171
Arkhipkin A.I. 2005. Statolith as black boxes (life recorders) in squid. Mar. Freshwater Res. 56(5): 573-583.
http://dx.doi.org/10.1071/MF04158
Arkhipkin A.I., Bizikov V.A. 2000. Role of the statolith in functioning of the acceleration receptor system in squids and sepioids. J. Zool. 250(1): 31-55.
http://dx.doi.org/10.1111/j.1469-7998.2000.tb00575.x
Arkhipkin A.I., Shcherbich Z.N. 2012. Thirty years’ progress in age determination of squid using statoliths. J. Mar. Biol. Ass. U. K. 92(06): 1389-1398.
http://dx.doi.org/10.1017/S0025315411001585
Bizikov V. A., Arkhipkin A. I. 1997. Morphology and microstructure of the gladius and statolith from the boreal Pacific giant squid Moroteuthis robusta (Oegopsida; Onychoteuthidae). J. Zool. 241(3): 475-492.
http://dx.doi.org/10.1111/j.1469-7998.1997.tb04839.x
Bower J.R., Ichii T. 2005. The red flying squid (Ommastrephes bartramii): A review of recent research and the fishery in Japan. Fish. Res. 76(1): 39-55.
http://dx.doi.org/10.1016/j.fishres.2005.05.009
Bower S.M., Margolis L. 1991. Potential use of helminth parasites in stock identification of flying squid, Ommastrephes bartramii, in North Pacific waters. Can. J. Zool. 69(4): 1124-1126.
http://dx.doi.org/10.1139/z91-158
Burke W.T., Freeberg M., Miles E.L. 1993. United Nations Resolutions on driftnet fishing: An unsustainable precedent for high seas and coastal fisheries management. Ocean Dev. Int. Law. 25: 127-186.
http://dx.doi.org/10.1080/00908329409546030
Cabanellas-Reboredo M., Alós J., Palmer M., et al. 2011. Simulating the indirect handline jigging effects on the European squid (Loligo vulgaris) in captivity. Fish. Res. 110(3): 435-440.
http://dx.doi.org/10.1016/j.fishres.2011.05.013
Cadrin S.X., Silva V.M. 2005. Morphometric variation of yellowtail flounder. ICES J. Mar. Sci. 62: 683-694.
http://dx.doi.org/10.1016/j.icesjms.2005.02.006
Castanhari G., Tomás A.R.G. 2012. Beak increment counts as a tool for growth studies of the common octopus Octopus vulgaris in southern Brazil. Bol. Inst. Pesca, São Paulo 38(4): 323-331.
Chen C.S. 2010. Abundance trends of two neon flying squid (Ommastrephes bartramii) stocks in the North Pacific. ICES J. Mar. Sci. 67(7): 1336-1345.
Chen C.S., Chiu T.S. 2003. Variations of life history parameters in two geographical groups of the neon flying squid, Ommastrephes bartramii, from the North Pacific. Fish. Res. 63: 349-366.
http://dx.doi.org/10.1016/S0165-7836(03)00101-2
Chen X.J., Tian S.Q., Ye X.C. 2002. Study on population structure of flying squid in Northwestern Pacific based on gray system theory. J. Shanghai Fish. Univ. 11(4): 335-341. (In Chinese with English Abstract)
Chen X.J., Liu B.L., Wang Y.G. 2009. Cephalopod of the world. Marine Press, Beijing, China, 1064 pp.
Chen X.J., Ma J., Liu B.L., et al. 2010. Effects of sexual maturity and body size on statolith shape of Ommastrephes bartramii in the Northwest Pacific Ocean. J. Fish. China. 30(6): 928-934. (In Chinese with English abstract)
Chen F., Chen X.J., Lu H.J., et al. 2011. Comparison of biological characteristics of Ommastrephes bartarmii between two different areas in the east central waters of North Pacific Ocean. J. Shanghai Fish. Univ. 20(5): 759-764. (In Chinese with English Abstract)
Chen X.J., Lu H.J., Liu B.L., et al. 2012. Species identification of Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis and Illex argentinus (Ommastrephidae) using beak morphological variables. Sci. Mar. 76(3): 473-481.
http://dx.doi.org/10.3989/scimar.03408.05B
Chen X.J., Li J.H., Liu B.L., et al. 2013. Age, growth and population structure of Jumbo flying squid, Dosidicus gigas, off the Costa Rica Dome. J. Mar. Biol. Ass. U. K. 93:567-573
http://dx.doi.org/10.1017/S0025315412000422
Cherel Y., Hobson K.A. 2005. Stable isotopes, beaks and predators: a new tool to study the trophic ecology of cephalopods, including giant and colossal squids. Proc. R. Soc. Lond. B. 272: 1601-1607.
http://dx.doi.org/10.1098/rspb.2005.3115
Cherel Y., Ridoux V., Spitz J., et al. 2009. Stable isotopes document the trophic structure of a deep-sea cephalopod assemblage including giant octopod and giant squid. Biol. Lett. 5: 364-367.
http://dx.doi.org/10.1098/rsbl.2009.0024
Clarke M.R. 1962. The identification of cephalopod “beaks” and the relationship between beak size and total body weight. Bull. Br. Mus. nat. Hist. Zool. 8: 419-480
Clarke M.R. 1978. The cephalopod statolith-introduction to its form. J. Mar. Biol. Ass. U.K. 58: 701-712.
http://dx.doi.org/10.1017/S0025315400041345
Clarke M.R. 1986. A handbook for the identification of cephalopod beaks. Clarendon Press, Oxford, UK, 273 pp.
Clarke M.R. 2003. Potential of statoliths for interpreting coleoid evolution: a brief review. Ber. Palä. Abh. 3: 37-47.
Crespi-Abril A.C., Morsan E.M., Barón P.J. 2010. Analysis of the ontogenetic variation in body and beak shape of the Illex argentinus inner shelf spawning groups by geometric morphometrics. J. Mar. Biol. Ass. U. K. 90(03): 547-553.
http://dx.doi.org/10.1017/S0025315409990567
Dommergues J.L., Neige P., Boletzky S.V. 2000. Exploration of morphospace using procrustes analysis in statoliths of cuttlefish and squid (Cephalopoda: Decabrachia)-evolutionary aspects of form disparity. Veliger-berkeley, 43(3): 265-276.
Durholtz M.D., Lipinski M.R., Przybylowicz W.J., et al. 1997. Nuclear microprobe mapping of statoliths of chokka squid Loligo vulgaris reynaudii d’Orbigny, 1845. Biol. Bull. 193(2): 125-140.
http://dx.doi.org/10.2307/1542758
Fang Z., Chen X.J., Lu H.J., et al. 2012. Morphological differences in statolith and beak between two spawning stocks for Illex argentinus. Acta. Ecol. Sin. 32(19): 5986-5997. (In Chinese with English Abstract)
http://dx.doi.org/10.5846/stxb201109101331
Francis R.I.C.C., Mattlin R.H. 1986. A possible pitfall in the morphometric application of discriminant analysis: measurement bias. Mar. Biol. 93(2): 311-313.
http://dx.doi.org/10.1007/BF00508269
Guerra A., Rodríguez-navarro A.B., González A.F., 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(7): 1425-1431.
Hanlon R.T., Messenger J.B. 1996. Cephalopod Behaviour. Cambridge University Press Cambridge, UK, 232 pp.
Ichii T., Mahapatra K., Sakai M., et al. 2004. Differing body size between the autumn and the winter–spring cohorts of neon flying squid (Ommastrephes bartramii) related to the oceanographic regime in the North Pacific: a hypothesis. Fish. Oceanogr. 13: 295-309
http://dx.doi.org/10.1111/j.1365-2419.2004.00293.x
Ichii T., Mahapatra K., Sakai M., et al. 2009. Life history of the neon flying squid: effect of the oceanographic regime in the North Pacific Ocean. Mar. Ecol. Prog. Ser. 378: 1-11.
http://dx.doi.org/10.3354/meps07873
Ikeda Y., Arai N., Sakamoto W., et al. 1996. PIXE analysis of trace elements in squid statoliths: comparison between Ommastrephidae and Loliginidae. International Journal of PIXE. 6: 537-542.
http://dx.doi.org/10.1142/S0129083596000594
Ikeda Y., Arai N., Sakamoto W., et al. 1997. Comparison on trace elements in squid statoliths of different species’ origin as available key for taxonomic and phylogenetic study. International Journal of PIXE. 7: 141-146.
http://dx.doi.org/10.1142/S0129083597000175
Ikeda Y., Arai N., Kidokoro H., et al. 2003. Strontium: calcium ratios in statoliths of Japanese common squid Todarodes pacificus (Cephalopoda: Ommastrephidae) as indicators of migratory behavior. Mar. Ecol. Prog. Ser. 251(1): 169-179.
http://dx.doi.org/10.3354/meps251169
Jackson G.D., Domeier M.L. 2003. The effects of an extraordinary El Ni-o/La Ni-a event on the size and growth of the squid Loligo opalescens off Southern California. Mar. Biol. 142(5): 925-935.
Katugin O.N. 2002. Patterns of genetic variability and population structure in the North Pacific squids Ommastrephes bartramii, Todarodes pacificus and Berryteuthis magister. Bull. Mar. Sci. 71(1): 383-420.
Kurosaka K., Yamashita H., Ogawa M., et al. 2012. Tentacle-breakage mechanism for the neon flying squid (Ommastrephes bartramii) during the jigging capture process. Fish. Res. 121: 9-16.
http://dx.doi.org/10.1016/j.fishres.2011.12.016
Lefkaditou E., Bekas P. 2004. Analysis of beak morphometry for the horned octopus Eledone cirrhosa (Cephalopoda: Octopoda) from the Thracian Sea (NE Mediterranean). Mediterr. Mar. Sci. 5(1): 143-149.
http://dx.doi.org/10.12681/mms.219
Lefkaditou E., Peristeraki P., Chartosia N., et al. 2011. Recent findings of Ommastrephes bartramii (Cephalopoda: Ommastrephidae) in the eastern Mediterranean and the implication on its range expansion. Mediterr. Mar. Sci. 12(2): 413-428.
http://dx.doi.org/10.12681/mms.41
Lipinski M. R., Underhill L. G. 1995. Sexual maturation in squid: quantum or continuum. S. Afr. J. Mar. Sci. 15: 207-223
http://dx.doi.org/10.2989/02577619509504844
Lleonart J., Salat J., Torres G.J. 2000. Removing allometric effects of body size in morphological analysis. J. Thero. Biol. 205: 85-93.
http://dx.doi.org/10.1006/jtbi.2000.2043
Lombarte A., Sanchez P., Morales-Nin B. 1997. Intraspecific shape variability in statoliths of three cephalopod species. Vie Milieu 47: 165-169.
Lombarte A., Rufino M.M., Sánchez P. 2006. Statolith identification of Mediterranean Octopodidae, Sepiidae, Loliginidae, Ommastrephidae and Enoploteuthidae based on warp analyses. J. Mar. Biol. Ass. U. K. 86(04): 767-771.
http://dx.doi.org/10.1017/S0025315406013683
Lu C.C., Ickeringill R. 2002. Cephalopod beak identification and biomass estimation techniques: tools for dietary studies of southern Australian finfishes. Mus. Victoria Sci. Rep. 6: 1-65.
Ma J., Chen X.J., Liu B.L., et al. 2009. Morphologic Features of Statolith for Ommastrephes bartramii in the Northwest Pacific Ocean. Periodical of Ocean University of China 39(2): 215-220. (In Chinese with English abstract)
Martínez P., Sanjuan A., Guerra A. 2002. Identification of Illex coindetii, I. illecebrosus and I. argentines (Cephalopoda: ommastrephidae) throughout the Atlantic Ocean by body and beak characters. Mar. Biol. 141:131-143
http://dx.doi.org/10.1007/s00227-002-0796-7
Mercer M.C., Misra R.K., Hurley G.V. 1980. Sex determination of the Ommastrephid squid Illex illecebrosus using beak morphometric. Can. J. Fish. Aquat. Sci. 37: 283-286
http://dx.doi.org/10.1139/f80-035
Miserez A., Rubin D., Waite J.H. 2010. Cross-linking Chemistry of Squid Beak. J. Biol. Chem. 285(49): 38115-38124.
http://dx.doi.org/10.1074/jbc.M110.161174
Moltschaniwskyj N.A. 1995. Changes in shape associated with growth in the loliginid squid Photololigo sp.: a morphometric approach. Can. J. Zool. 73(7): 1335-1343.
http://dx.doi.org/10.1139/z95-157
Murakami K., Watanabe Y., Nakata J. 1981. Growth, distribution and migration of flying squid (Ommastrephes bartrami) in the North Pacific. In: Mishima, S. (ed.), Pelagic animals and environments around the Subarctic Boundary in North Pacific. Hokkaido University, Research Institute of North Pacific Fisheries, Hakodate, pp. 161-179 (In Japanese with English abstract)
Murata M. 1990. Oceanic resources of squids. Mar. Behav. Phys. 18: 19-71.
http://dx.doi.org/10.1080/10236249009378779
Murata M., Hayase S. 1993. Life history and biological information on flying squid (Ommastrephes bartramii) in the North Pacific Ocean. Bull. Int. Nat. North Pacific Comm. 53: l47-182.
Nagasawa K., Mori J., Okamura H. 1998. Parasites as biological tags of stocks of neon flying squid (Ommastrephes bartramii) in the North Pacific Ocean, In: Okutani, T. (1998), Contributed papers to International Symposium on Large Pelagic Squids, July 18-19, 1996, for JAMARC’s 25th anniversary of its foundation, pp. 49-64.
Neige P. 2006. Morphometrics of hard structures in cuttlefish. Vie Milieu 56(2): 121-127.
Neige P., Dommergues, J.L. 2002. Disparity of beaks and statoliths of some coleoids a morphometric approach to depict shape differentiation. Gabh. der Geol. Bun. 57: 393-399
O’Dor R.K., Dawe E.G. 1998. Illex illecebrosus. In: Rodhouse P.G., Dawe E.G., O’Dor R.K. (eds), Squid recruitment dynamics. FAO Fish Tech Pap 376, Rome, pp. 77-104.
O’Dor R.K., Hoar J.A. 2000. Does geometry limit squid growth? ICES J. Mar. Sci. 57(1): 8-14.
http://dx.doi.org/10.1006/jmsc.1999.0502
Ogden R.S., Allcock A.L., Wats P.C., et al. 1998. The role of beak shape in octopodid taxonomy. S. Afr. J. Mar. Sci. 20(1): 29-36.
http://dx.doi.org/10.2989/025776198784126476
Perales-Raya C., Jurado-Ruzafa A., Bartolomé A., et al. 2014. Age of spent Octopus vulgaris and stress mark analysis using beaks of wild individuals. Hydrobiologia 725(1): 105-114.
http://dx.doi.org/10.1007/s10750-013-1602-x
Piatkowski U., Pütz K., Heinemann H. 2001. Cephalopod prey of king penguins (Aptenodytes patagonicus) breeding at Volunteer Beach, Falkland Islands, during austral winter 1996. Fish. Res. 52(1): 79-90.
http://dx.doi.org/10.1016/S0165-7836(01)00232-6
Pierce G. J., Hastie L. C., Guerra A., et al. 1994. Morphometric variation in Loligo forbesi and Loligo vulgaris: regional, seasonal, sex, maturity and worker differences. Fish. Res. 21(1): 127-148.
http://dx.doi.org/10.1016/0165-7836(94)90100-7
Pineda S.E., Hernández D.R., Brunetti N.E., et al. 2002. Morphological identification of two Southwest Atlantic Loliginid squids: Loligo gahi and Loligo sanpaulensis. Rev. Invest. Desarr. Pesq. 15: 67-84.
Raya C.P., Hernández-González C.L. 1998. Growth lines within the beak microstructure of the octopus Octopus vulgaris Cuvier, 1797. S. Afr. J. Mar. Sci. 20(1): 135-142.
http://dx.doi.org/10.2989/025776198784126368
Raya C.P., Bartolomé A., García-Santamaría M.T., et al. 2010. Age estimation obtained from analysis of octopus (Octopus vulgaris Cuvier, 1797) beaks: Improvement and comparisons. Fish. Res. 106: 171-176
http://dx.doi.org/10.1016/j.fishres.2010.05.003
Rencher A.C. 2002. Methods of Multivariate Analysis, 2nd edition. John Wiley & Sons, Inc, New York.
http://dx.doi.org/10.1002/0471271357
Rocha F., Guerra A., González A.F. 2001. A review of reproductive strategies in cephalopods. Biol. Rev. 76: 291-304.
http://dx.doi.org/10.1017/S1464793101005681
Rodhouse P.G. 2005. World squid resources. Review of the state of world marine fishery resources, FAO Tech. Rep. C2:175.
Ruiz-Cooley R.I., Ballance L.T., McCarthy M.D. 2013. Range Expansion of the Jumbo Squid in the NE Pacific: δ15N Decrypts Multiple Origins, Migration and Habitat Use. PloS One. 8(3): e59651.
http://dx.doi.org/10.1371/journal.pone.0059651
Sajina A.M., Chakraborty S.K., Jaiswar A.K., et al. 2011. Stock structure analysis of Megalaspis cordyla (Linnaeus, 1758) along the Indian coast based on truss network analysis. Fish. Res. 108: 100-105
http://dx.doi.org/10.1016/j.fishres.2010.12.006
Sánchez P. 1995. Age and growth of Illex coindetii. ICES J. Mar. Sci. 199: 441-444
Smale M.J., Clarke M.R., Klages N.T.W., et al. 1993. Octopod beak identification—resolution at a regional level (Cephalopoda, Octopoda: Southern Africa). S. Afr. J. Mar. Sci. 13(1): 269-293.
http://dx.doi.org/10.2989/025776193784287338
Uchikawa K., Sakai M., Wakabayashi T., et al. 2009. The relationship between paralarval feeding and morphological changes in the proboscis and beaks of the neon flying squid Ommastrephes bartramii. Fish. Sci. 75(2): 317-323.
http://dx.doi.org/10.1007/s12562-008-0036-2
Vega M.A., Rocha F.J., Guerra A., et al. 2002. Morphological difference between the Patagonian squid Loligo gahi populations from the Pacific and Atlantic Oceans. Bull. Mar. Sci. 71(2): 903-913
Villanueva R. 1992. Interannual growth differences in the oceanic squid Todarodes angolensis Adam in the northern Benguela upwelling system, based on statolith growth increment analysis. J. Exp. Mar. Biol. Ecol., 159(2): 157-177.
http://dx.doi.org/10.1016/0022-0981(92)90034-8
Wang Y.G., Chen X.J. 1998. The current exploitation of cephalopod resources in the world and the development of Chinese distant-water squid jigging fisheries. J. Shanghai Fish. Univ. 7: 285-287 (In Chinese with English abstract)
Wang Y.G., Chen X.J. 2005. The Resource and Biology of Economic Oceanic Squid in the World. Ocean Press, Beijing, pp. 79-295.
Wolff G.A. 1984. Identification and estimation of size from the beaks of 18 species of cephalopods from the Pacific Ocean. NOAA Technical Report NMFS 17: 1-50
Xavier J.C., Cherel Y. 2009. Cephalopod beak guide for the southern ocean. British Antarctic Survey Press, Cambridge, UK, 129 pp.
Xavier J.C., Clarke M.R., Magalhães M.C., et al. 2007. Current status of using beaks to identify cephalopods: III International Workshop and training course on Cephalopod beaks, Faial Island, Azores. Arquipélago 24: 41-48.
Xavier J.C., Phillips R.A., Cherel Y. 2011. Cephalopods in marine predator diet assessments: why identifying upper and lower beaks is important. ICES J. Mar. Sci. 68(9): 1857-1864.
http://dx.doi.org/10.1093/icesjms/fsr103
Yatsu A. 2000. Age estimation of four oceanic squids, Ommastrephes bartramii, Dosidicus gigas, Sthenoteuthis oualaniensis, and Illex argentinus (Cephalopoda, Ommastrephidae) based on statolith Microstructure. Jpn. Agri. Res. Quart. 34: 75-80.
Yatsu A., Mori J. 2000. Early growth of the autumn cohort of neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean. Fish. Res. 45: 189-194.
http://dx.doi.org/10.1016/S0165-7836(99)00112-5
Yatsu A., Midorikawa S., Shimada T., et al. 1997. Age and growth of the neon flying squid, Ommastrephes bartrami, in the North Pacific Ocean. Fish. Res. 29: 257-270
http://dx.doi.org/10.1016/S0165-7836(96)00541-3
Yatsu A., Tanaka H., Mori J. 1998a. Population structure of the neon flying squid, Ommastrephes bartramii, in the North Pacific. In: Okutani T. (ed.), Contributed Papers to International Symposium on Large Pelagic Squids. Japan Marine Fishery Resources Research Center, Tokyo, pp. 31-48.
Yatsu A., Mochioka N., Morishita K., et al. 1998b. Strontium/calcium ratios in statoliths of the neon flying squid, Ommastrephes bartrami (Cephalopoda) in the North Pacific Ocean. Mar. Biol. 131: 275-282.
http://dx.doi.org/10.1007/s002270050320