Stock discrimination and vertical migration of two gadoid species based on otolith stable isotope analysis

Abstract

Pacific cod (Gadus macrocephalus) and walleye pollock (Gadus chalcogrammus) are commercially important species in bottom trawl fisheries in Japanese waters. However, their first-year life history processes and whether current management units align with the scale of their biological processes remain unclear. To investigate this, we conducted otolith macrostructural analysis on age-1 Pacific cod and walleye pollock collected from the Pacific Ocean off Hokkaido, Aomori Prefecture, and Sendai Bay (Honshu). Stable isotope analysis (δ¹³C, δ¹⁸O) was performed on both the inside and outside of the otolith check using isotope ratio mass spectrometry (IRMS) and secondary ion mass spectrometry (SIMS). The δ¹⁸O shift across the check, combined with higher δ¹³C and δ¹⁸O values on the outside of the check compared to the inside, suggests that both species undergo a habitat transition during their first-year settlement. Stock discrimination was more successful for individuals from Hokkaido and Sendai Bay than for those from Aomori Prefecture. When isotope data from the two southern regions (Aomori Prefecture and Sendai Bay) were combined into a single Honshu stock, the overall classification success increased from 71 to 76%, with Cohen’s Kappa (κ) rising significantly from 0.31 to 0.60. These results indicate that managing stocks based on geographic distribution is both scientifically valid and practical though inter-stock admixture should be taken into account. Overall, otolith stable isotope analysis proves to be a valuable tool for identifying key life history events, stock structure, and migration patterns in these gadoid species.

Keywords

Gadus macrocephalus / Gadus chalcogrammus / Otolith / Stable isotopes / Stock structure / Vertical migration

References

[1] Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001, 26: 32-46

[2] Artetxe-Arrate I, Fraile I, Clear N, Darnaude AM, Dettman DL, Pécheyran C, Farley J, Murua H. Discrimination of yellowfin tuna Thunnus albacares between nursery areas in the Indian Ocean using otolith chemistry. Mar Ecol Prog Ser. 2021, 673: 165-181.

[3] Artetxe-Arrate I, Fraile I, Farley J, Darnaude AM, Clear N, Rodríguez-Ezpeleta N, Dettman DL, Pécheyran C, Krug I, Médieu A. Otolith chemical fingerprints of skipjack tuna (Katsuwonus pelamis) in the Indian Ocean: first insights into stock structure delineation. PLoS ONE. 2021, 16. e0249327

[4] Bastow T, Jackson G, Edmonds J. Elevated salinity and isotopic composition of fish otolith carbonate: stock delineation of pink snapper, Pagrus auratus, in Shark Bay, Western Australia. Mar Biol. 2002, 141: 801-806.

[5] Becker M, Andersen N, Erlenkeuser H, Humphreys MP, Tanhua T, Körtzinger A. An internally consistent dataset of δ¹³C-DIC in the North Atlantic Ocean–NAC13v1. Earth Syst Sci Data. 2016, 8: 559-570.

[6] Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.

[7] Brophy D, Rodríguez-Ezpeleta N, Fraile I, Arrizabalaga H. Combining genetic markers with stable isotopes in otoliths reveals complexity in the stock structure of Atlantic bluefin tuna (Thunnus thynnus). Sci Rep. 2020, 10: 14675.

[8] Campana SE. How reliable are growth back-calculations based on otoliths?. Can J Fish Aquat Sci. 1990, 472219-2227.

[9] Campana SE. Chemistry and composition of fish otoliths: pathways, mechanisms and applications. Mar Ecol Prog Ser. 1999, 188: 263-297.

[10] Campana SE, Neilson JD. Microstructure of fish otoliths. Can J Fish Aquat Sci. 1985, 42: 1014-1032.

[11] Campana SE, Thorrold SR. Otoliths, increments, and elements: keys to a comprehensive understanding of fish populations?. Can J Fish Aquat Sci. 2001, 58: 30-38.

[12] Casas JM, Piñeiro C. Growth and age estimation of greater fork-beard (Phycis blennoides Brünnich, 1768) in the north and northwest of the Iberian Peninsula (ICES Division VIIIc and IXa). Fish Res. 2000, 47: 19-25.

[13] Chung M-T, Trueman CN, Godiksen JA, Grønkjær P. Otolith δ¹³C values as a metabolic proxy: approaches and mechanical underpinnings. Mar Freshwater Res. 2019, 70: 1747-1756.

[14] Currey LM, Heupel MR, Simpfendorfer CA, Williams AJ. Inferring movement patterns of a coral reef fish using oxygen and carbon isotopes in otolith carbonate. J Exp Mar Biol Ecol. 2014, 456: 18-25.

[15] Darnaude AM, Sturrock A, Trueman CN, Mouillot D, Campana SE, Hunter E. Listening in on the past: what can otolith δ¹⁸O values really tell us about the environmental history of fishes?. PLoS ONE. 2014, 9. e108539

[16] Degens ET, Deuser WG, Haedrich RL. Molecular structure and composition of fish otoliths. Mar Biol. 1969, 2: 105-113.

[17] DiMaria R, Miller J, Hurst T. Temperature and growth effects on otolith elemental chemistry of larval Pacific cod, Gadus macrocephalus. Environ Biol Fishes. 2010, 89: 453-462.

[18] Dufour E, Gerdeaux D, Wurster CM. Whitefish (Coregonus lavaretus) respiration rate governs intra-otolith variation of δ¹³C values in Lake Annecy. Can J Fish Aquat Sci. 2007, 64: 1736-1746.

[19] Elsdon TS, Wells BK, Campana SE, Gillanders BM, Jones CM, Limburg KE, Secor DH, Thorrold SR, Walther BD. Otolith chemistry to describe movements and life-history parameters of fishes: hypotheses, assumptions, limitations and inferences. Oceanography and Marine Biology. 2008CRC Press303336

[20] Fukuda S. The mark and release of adult Pacific cod Gadus macrocephalus in the entrance of Mutsu Bay, Aomori prefecture, Japan. Saibai Giken. 1985, 14: 71-77

[21] Funamoto T, Yamamura O, Shida O, Itaya K, Mori K, Hiyama Y, Sakurai Y. Comparison of factors affecting recruitment variability of walleye pollock Theragra chalcogramma in the Pacific Ocean and the Sea of Japan off northern Japan. Fish Sci. 2014, 80: 117-126.

[22] Hane Y, Kimura S, Yokoyama Y, Miyairi Y, Ushikubo T, Ishimura T, Ogawa N, Aono T, Nishida K. Reconstruction of temperature experienced by Pacific bluefin tuna Thunnus orientalis larvae using SIMS and microvolume CF-IRMS otolith oxygen isotope analyses. Mar Ecol Prog Ser. 2020, 649: 175-188.

[23] Hane Y, Ushikubo T, Yokoyama Y, Miyairi Y, Kimura S. Natal origin of Pacific bluefin tuna Thunnus orientalis determined by SIMS oxygen isotope analysis of otoliths. PLoS ONE. 2022, 17. e0272850

[24] Hashimoto R. Investigation of feeding habits and variation of inhabiting depths with cod (Gadus macrocephalus) distributed on the northeastern fishing ground in Japan. Bull Tohoku Reg Fish Res Lab. 1974, 33: 51-67

[25] Helser T, Kastelle C, Crowell A, Ushikubo T, Orland IJ, Kozdon R, Valley JW. A 200-year archaeozoological record of Pacific cod (Gadus macrocephalus) life history as revealed through ion microprobe oxygen isotope ratios in otoliths. J Archaeol Sci Rep. 2018, 21: 1236-1246

[26] Helser TE, Kastelle CR, McKay JL, Orland IJ, Kozdon R, Valley JW. Evaluation of micromilling/conventional isotope ratio mass spectrometry and secondary ion mass spectrometry of δ¹⁸O values in fish otoliths for sclerochronology. Rapid Commun Mass Spectrom. 2018, 32: 1781-1790.

[27] Høie H, Andersson C, Folkvord A, Karlsen Ø. Precision and accuracy of stable isotope signals in otoliths of pen-reared cod (Gadus morhua) when sampled with a high-resolution micromill. Mar Biol. 2004, 144: 1039-1049.

[28] Høie H, Otterlei E, Folkvord A. Temperature-dependent fractionation of stable oxygen isotopes in otoliths of juvenile cod (Gadus morhua L.). ICES J Mar Sci. 2004, 61: 243-251.

[29] Honda S, Oshima T, Nishimura A, Hattori T. Movement of juvenile walleye pollock, Theragra chalcogramma, from a spawning ground to a nursery ground along the Pacific coast of Hokkaido, Japan. Fish Oceanogr. 2004, 13: 84-98.

[30] Hüssy K, Mosegaard H, Jessen F. Effect of age and temperature on amino acid composition and the content of different protein types of juvenile Atlantic cod (Gadus morhua) otoliths. Can J Fish Aquat Sci. 2004, 61: 1012-1020.

[31] Ihssen P, Booke H, Casselman J, McGlade J, Payne N, Utter F. Stock identification: materials and methods. Can J Fish Aquat Sci. 1981, 38: 1838-1855.

[32] Jones JB, Campana SE. Stable oxygen isotope reconstruction of ambient temperature during the collapse of a cod (Gadus morhua) fishery. Ecol Appl. 2009, 19: 1500-1514.

[33] Kalish JM. ¹³C and ¹⁸O isotopic disequilibria in fish otoliths: metabolic and kinetic effects. Mar Ecol Prog Ser. 1991, 75: 191-203.

[34] Kastelle CR, Helser TE, McKay JL, Johnston CG, Anderl DM, Matta ME, Nichol DG. Age validation of Pacific cod (Gadus macrocephalus) using high-resolution stable oxygen isotope (δ¹⁸O) chronologies in otoliths. Fish Res. 2017, 185: 43-53.

[35] Kastelle CR, Helser TE, Laurel BJ, Copeman LA, Stone KR, McKay JL. Oxygen isotope fractionation in otoliths: experimental results from four North Pacific and Arctic gadid species. Mar Ecol Prog Ser. 2022, 686: 159-175.

[36] Kerr LA, Hintzen NT, Cadrin SX, Clausen LW, Dickey-Collas M, Goethel DR, Hatfield EM, Kritzer JP, Nash RD. Lessons learned from practical approaches to reconcile mismatches between biological population structure and stock units of marine fish. ICES J Mar Sci. 2017, 74: 1708-1722.

[37] Kim S-T, O'Neil JR. Equilibrium and nonequilibrium oxygen isotope effects in synthetic carbonates. Geochim Cosmochim Acta. 1997, 61: 3461-3475.

[38] Kitagawa D. Monitoring on the demersal fish resources in the Tohoku area. Kaiyo Mon. 2002, 34: 793-798

[39] Kritzer JP, Liu OR. Cadrin SX, Kerr LA, Mariani S. Fishery management strategies for addressing complex spatial structure in marine fish stocks. Stock identification methods. 20142New York, Academic Press2957.

[40] Kroopnick P. The distribution of ¹³C of ΣCO₂ in the world oceans. Deep-Sea Res a: Oceanogr Res Pap. 1985, 32: 57-84.

[41] Legendre P, Legendre L (1998) Numerical ecology. Second English edition. Elsevier, Amsterdam.

[42] LeGrande AN, Schmidt GA. Global gridded data set of the oxygen isotopic composition in seawater. Geophys Res Lett. 2006, 33: 1-5.

[43] Li J, Jiang F, Wu R, Zhang C, Tian Y, Sun P, Yu H, Liu Y, Ye Z, Ma S. Tidally induced temporal variations in growth of young‐of‐the‐year Pacific cod in the Yellow Sea. J Geophys Res Oceans. 2021, 126. e2020JC016696

[44] Madea T, Takahashi T, Ueno M. Annual life period of the adult Alaska pollack in the adjacent waters of the Funka Bay, Hokkaido (in Japanese). Bull Jap Soc Sci Fish. 1981, 47: 741-746.

[45] Maeda T. Life cycle and behavior of adult pollock (Theragra chalcogramma [Pallas]) in water adjacent to Funka Bay, Hokkaido Island. Bull Int North Pac Fish Commn. 1986, 45: 39-65

[46] Maeda T, Takahashi T, Ueno M. Behaviour in each life period of adult Alaska pollack in the adjacent waters of Funka Bay, Hokkaido [Japan, Theragra chalcogramma]. Bull Jap Soc Sci Fish. 1983, 49: 577-585.

[47] Martino JC, Doubleday ZA, Chung M-T, Gillanders BM. Experimental support towards a metabolic proxy in fish using otolith carbon isotopes. J Exp Biol. 2020, 223. jeb217091

[48] McQueen K, Hrabowski J, Krumme U. Age validation of juvenile cod in the Western Baltic Sea. ICES J Mar Sci. 2019, 76: 430-441.

[49] Mercier L, Darnaude AM, Bruguier O, Vasconcelos RP, Cabral HN, Costa MJ, Lara M, Jones DL, Mouillot D. Selecting statistical models and variable combinations for optimal classification using otolith microchemistry. Ecol Appl. 2011, 21: 1352-1364.

[50] Mishima S. Stock and biological characteristics of the cod Gadus macrocephalus TILESIUS around Japan (in Japanese). Annu Rep INPFC. 1989, 42: 172-179

[51] Morioka T, Kuwada H. The upper limit of inhabiting temperature and the diet of juvenile Pacific cod Gadus macrocephalus in the northern part of Nanao Bay and its vicinity, Japan. Bull Jpn Soc Sci Fish. 2002, 68: 345-350.

[52] Nakamura T, Tanaka T, Kimura O, Fujii Y, Haraguchi K, Endo T. Comparison of radiocesium and stable isotope ratios of carbon and nitrogen among three stocks of Pacific cod (Gadus macrocephalus) around Hokkaido, Japan. Mar Pollut Bull. 2018, 127: 39-44.

[53] Nakatani T. Studies on the early life history of walleye pollock Theragra chalcogramma in Funka Bay and vicinity, Hokkaido. Mem Fac Fish Hokkaido Univ. 1988, 35: 1-46

[54] Nakatani T. Year class strength and early life hisonory of the Pacific population of walleye pollock, Gadus chalcogrammus, in Japan. Mem Fac Fish Sci Hokkaido Univ. 2016, 58: 1-11

[55] Nakatani T, Maeda T. Transport process of the Alaska pollack eggs in Funka Bay and the adjacent waters, Hokkaido (in Japanese). Bull Jap Soc Sci Fish. 1981, 47: 1115-1118.

[56] Narimatsu Y, Ueda Y, Okuda T, Hattori T, Fujiwara K, Ito M. The effect of temporal changes in life-history traits on reproductive potential in an exploited population of Pacific cod, Gadus macrocephalus. ICES J Mar Sci. 2010, 67: 1659-1666.

[57] Narimatsu Y, Kakehi S, Ito S-i, Okazaki Y, Inagawa R, Yano T. Impact of the Great East Japan Earthquake tsunami on growth and survival of Pacific cod (Gadus macrocephalus). Can J Fish Aquat Sci. 2015, 72: 1629-1638.

[58] Narimatsu Y, Sohtome T, Yamada M, Shigenobu Y, Kurita Y, Hattori T, Inagawa R. Nakata K, Sugisaki H. Why do the radionuclide concentrations of Pacific cod depend on the body size?. Impacts Fukushima nuclear accident on fish and fishing grounds. 2015, Tokyo, Springer123137.

[59] Nelson CS, Northcote TG, Hendy CH. Potential use of oxygen and carbon isotopic composition of otoliths to identify migratory and non-migratory stocks of the New Zealand common smelt: a pilot study. N Z J Mar Freshw Res. 1989, 23: 337-344.

[60] Panfili J, de Pontual H, Troadec H, Wrigh PJ. Manual of fish sclerochronology. 2002, France, Ifremer-IRD coedition

[61] Park HJ, Park TH, Lee C-I, Kang C-K. Ontogenetic shifts in diet and trophic position of walleye pollock, Theragra chalcogramma, in the western East Sea (Japan Sea) revealed by stable isotope and stomach content analyses. Fish Res. 2018, 204: 297-304.

[62] Popper AN, Ramcharitar J, Campana SE. Why otoliths? Insights from inner ear physiology and fisheries biology. Mar Freshw Res. 2005, 56: 497-504.

[63] Rehberg-Haas S, Hammer C, Hillgruber N, Hüssy K, Temming A. Otolith microstructure analysis to resolve seasonal patterns of hatching and settlement in western Baltic cod. ICES J Mar Sci. 2012, 69: 1347-1356.

[64] Reiss H, Hoarau G, Dickey-Collas M, Wolff WJ. Genetic population structure of marine fish: mismatch between biological and fisheries management units. Fish Fish. 2009, 10: 361-395.

[65] Riccioni G, Landi M, Ferrara G, Milano I, Cariani A, Zane L, Sella M, Barbujani G, Tinti F. Spatio-temporal population structuring and genetic diversity retention in depleted Atlantic bluefin tuna of the Mediterranean Sea. Proc Natl Acad Sci U S A. 2010, 107: 2102-2107.

[66] Sakuma K, Yoshikawa A, Goto T, Fujiwara K, Ueda Y. Delineating management units for Pacific cod (Gadus macrocephalus) in the Sea of Japan. Estuar Coast Shelf Sci. 2019, 229. 106401

[67] Sakurai Y, Hattori T. Reproductive behavior of Pacific cod in captivity. Fish Sci. 1996, 62: 222-228.

[68] Schade FM, Weist P, Krumme U. Evaluation of four stock discrimination methods to assign individuals from mixed-stock fisheries using genetically validated baseline samples. Mar Ecol Prog Ser. 2019, 627: 125-139.

[69] Secor DH. Cadrin SX, Kerr LA, Mariani S. The unit stock concept: bounded fish and fisheries. Stocks identify methods. 20142Academic Press728.

[70] Shephard S, Trueman C, Rickaby R, Rogan E. Juvenile life history of NE Atlantic orange roughy from otolith stable isotopes. Deep Sea Res Part I Oceanogr Res Pap. 2007, 54: 1221-1230.

[71] Sherwood GD, Rose GA. Stable isotope analysis of some representative fish and invertebrates of the Newfoundland and Labrador continental shelf food web. Estuar Coast Shelf Sci. 2005, 63: 537-549.

[72] Shimada A. Seasonal movements of Pacific cod, Gadus macrocephalus, in the eastern Bering Sea and adjacent waters based on tag-recapture data. Fish Bull. 1994, 92: 800-816

[73] Shimizu M, Isoda Y. The transport process of walleye pollock eggs into Funka Bay in winter. Bull Jpn Soc Fish Oceanogr. 1997, 61: 134-143

[74] Stanley RR, Bradbury IR, DiBacco C, Snelgrove PV, Thorrold SR, Killen SS. Environmentally mediated trends in otolith composition of juvenile Atlantic cod (Gadus morhua). ICES J Mar Sci. 2015, 72: 2350-2363.

[75] Suda A, Nagata N, Sato A, Narimatsu Y, Nadiatul H, Kawata M. Genetic variation and local differences in Pacific cod Gadus macrocephalus around Japan. J Fish Biol. 2017, 90: 61-79.

[76] Tagliabue A, Bopp L. Towards understanding global variability in ocean carbon‐13. Glob Biogeochem Cycles. 2008, 22. GB1025

[77] Takatsu T, Nakatani T, Mutoh T, Takahashi T. Feeding habits of Pacific cod larvae and juveniles in Mutsu Bay, Japan. Fish Sci. 1995, 61: 415-422.

[78] Takatsu T, Yoshida Y, Kooka K, Sugimoto K, Takahashi T. Spatial and temporal distribution of Pacific cod Gadus macrocephalus juveniles in Mutsu Bay, Japan. Bull Jap Soc Fish Oceanogr. 2001, 65: 6-14

[79] Tanaka H, Hamatsu T, Mori K. Comparison of potential fecundity models for walleye pollock Gadus chalcogrammus in the Pacific waters off Hokkaido, Japan. J Fish Biol. 2017, 90: 236-248.

[80] Thorrold SR, Campana SE, Jones CM, Swart PK. Factors determining δ¹³C and δ¹⁸O fractionation in aragonitic otoliths of marine fish. Geochim Cosmochim Acta. 1997, 61: 2909-2919.

[81] Titus K, Mosher JA, Williams BK. Chance-corrected classification for use in discriminant analysis: ecological applications. Am Midl Nat. 1984, 111: 1-7.

[82] Trueman CN, MacKenzie K, Palmer M. Identifying migrations in marine fishes through stable-isotope analysis. J Fish Biol. 2012, 81: 826-847.

[83] Trueman CN, Chung M-T, Shores D. Ecogeochemistry potential in deep time biodiversity illustrated using a modern deep-water case study. Phil Transac R Soc b: Biol Sci. 2016, 371. 20150223

[84] Tsuji S. Alaska pollack population, Theragra chalcogramma, of Japan and its adjacent waters, I: Japanese fisheries and population studies. Mar Freshw Behav Physiol. 1989, 15: 147-205.

[85] Ueno M, Nakamura H, Maeda T, Hirakawa H. Ecological studies on the alaska pollack in the adjacent waters of the Funka Bay, Hokkaido [Japan], 1: feeding season. Bull Jap Soc Sci Fish. 1975.

[86] Weidman CR, Millner R. High-resolution stable isotope records from North Atlantic cod. Fish Res. 2000, 46: 327-342.

[87] Wright JD. Global climate change in marine stable isotope records. Quaternary geochronology methods and applications. 2000, Washington D.C, American Geophysical Union671682

[88] Wu R, Li J, Zhang C, Pan X, Jiang T, Yang J, Liu S, Tian Y. Fine-scale variability in otolith chemistry: application to the life history analysis of Pacific cod (Gadus macrocephalus) in the Yellow Sea. Estuar Coast Shelf Sci. 2021, 257. 107392

[89] Wu R, Zhu Q, Katayama S, Tian Y, Li J, Fujiwara K, Narimatsu Y. Early life history affects fish size mainly by indirectly regulating the growth during each stage: a case study in a demersal fish. Mar Life Sci Technol. 2023, 5: 75-84.

[90] Yamamura O, Watanabe K, Shimazaki K. Feeding habits of Pacific cod, Gadus macrocephalus, off eastern Hokkaido, north Japan. Proc NIPR Symp Polar Biol. 1993, 6: 44-54

[91] Yamamura O, Funamoto T, Chimura M, Honda S, Oshima T. Interannual variation in diets of walleye pollock in the Doto area, in relation to climate variation. Mar Ecol Prog Ser. 2013, 491: 221-234.

[92] Yang Y-S, Kang S-K, Kim S-A, Kim S-S. Relationship between oxygen isotopic composition of walleye pollock (Theragra chalcogramma) otoliths and seawater temperature. Ocean Polar Res. 2008, 30: 249-258.

[93] Yoon T. Reproductive cycle of female walleye pollock, Theragra chalcogramma (Pallas), in the adjacent waters of Funka Bay, Hokkaido. Bull Fac Fish Hokkaido Univ. 1981, 32: 22-38

[94] Yoseda K. Reproduction of adult Pacific cod and distribution of its eggs, larvae and juveniles off the Noto-Island, Ishikawa Prefecture. Saibai Giken. 1992, 21: 21-30

RIGHTS & PERMISSIONS

The Author(s)