Growth of Mediterranean young-of-the-year bluefin tuna Thunnus thynnus (Scombridae): regional differences and hatching periods

INTRODUCTION

The Atlantic bluefin tuna (Thunnus thynnus) is a highly migratory species with a wide distribution throughout the North Atlantic Ocean and the Mediterranean Sea. The species comprises two different populations, one spawning in the Mediterranean and one in the Gulf of Mexico (Rooker et al. 2008Rooker J.R., Secor D.H., De Metrio G., et al. 2008. Natal homing and connectivity in Atlantic bluefin tuna populations. Science 322: 742-744. https://doi.org/10.1126/science.1161473 , Rodríguez-Ezpeleta et al. 2019Rodríguez-Ezpeleta N., Díaz-Arce N., Walter J.F., et al. 2019. Determining natal origin for improved management of Atlantic bluefin tuna. Front. Ecol. Environ. 17: 439-444 https://doi.org/10.1002/fee.2090 ). These populations share common feeding grounds and exhibit a wide spatial overlap throughout the Atlantic, as revealed by satellite tags and genetic and microchemistry analyses (Boustany et al. 2008Boustany A.M., Reeb C.A., Block B.A. 2008. Mitochondrial DNA and electronic tracking reveal population structure of Atlantic bluefin tuna (Thunnus thynnus). Mar. Biol. 156: 13-24. https://doi.org/10.1007/s00227-008-1058-0 , Rooker et al. 2014Rooker J.R., Arrizabalaga H., Fraile I., et al. 2014. Crossing the line: migratory and homing behaviors of Atlantic bluefin tuna. Mar. Ecol. Prog. Ser. 504: 265-276. https://doi.org/10.3354/meps10781 , Rodríguez-Ezpeleta et al. 2019Rodríguez-Ezpeleta N., Díaz-Arce N., Walter J.F., et al. 2019. Determining natal origin for improved management of Atlantic bluefin tuna. Front. Ecol. Environ. 17: 439-444 https://doi.org/10.1002/fee.2090 ). Despite the complex population structure of Atlantic bluefin tuna (ABFT), the International Commission for the Conservation of Atlantic Tunas (ICCAT) manages it as two different stocks separated by the 45°W meridian. The complexity of the ABFT structure is becoming increasingly conspicuous as knowledge of it progresses. The assumption that there are only two spawning areas has been questioned (Lutcavage et al. 1999Lutcavage M.E., Brill R.W., Skomal G.B., et al. 1999. Results of pop-up satellite tagging of spawning size class fish in the Gulf of Maine: do North Atlantic bluefin tuna spawn in the mid-Atlantic? Can. J. Fish. Aquat. Sci. 56: 173-177. https://doi.org/10.1139/f99-016 ). Indirect evidence based on reproductive studies (Goldstein et al. 2007Goldstein J., Heppell S., Cooper A., et al. 2007. Reproductive status and body condition of Atlantic bluefin tuna in the Gulf of Maine, 2000-2002. Mar. Biol. 151: 2063-2075. https://doi.org/10.1007/s00227-007-0638-8 ) and tag data (Block et al. 2005Block B.A., Teo S.L., Walli A., et al. 2005. Electronic tagging and population structure of Atlantic bluefin tuna. Nature 434: 1121-1127. https://doi.org/10.1038/nature03463 , Walli et al. 2009Walli A., Teo S.L., Boustany A., et al. 2009. Seasonal movements, aggregations and diving behavior of Atlantic bluefin tuna (Thunnus thynnus) revealed with archival tags. PLoS ONE 4: e6151. https://doi.org/10.1371/journal.pone.0006151 , Galuardi et al. 2010Galuardi B., Royer F., Golet W., et al. 2010. Complex migration routes of Atlantic bluefin tuna (Thunnus thynnus) question current population structure paradigm. Can. J. Fish. Aquat. Sci. 67: 966-976. https://doi.org/10.1139/F10-033 ) point to the existence of undocumented spawning grounds which, at least for the Slope Sea, was definitively confirmed by Richardson et al. (2016)Richardson D.E., Marancik K.E., Guyon J.R., et al. 2016. Discovery of a spawning ground reveals diverse migration strategies in Atlantic bluefin tuna (Thunnus thynnus). Proc. Nat. Acad. Sci. 113: 3299-3304. https://doi.org/10.1073/pnas.1525636113 . In addition, several authors have pointed out the possibility of additional spawning areas in the Atlantic, such as in the Canary and Azores Islands (Lutcavage et al. 1999Lutcavage M.E., Brill R.W., Skomal G.B., et al. 1999. Results of pop-up satellite tagging of spawning size class fish in the Gulf of Maine: do North Atlantic bluefin tuna spawn in the mid-Atlantic? Can. J. Fish. Aquat. Sci. 56: 173-177. https://doi.org/10.1139/f99-016 , Di Natale and Idrissi 2012Di Natale A., Idrissi M. 2012. 2011 ICCAT Atlantic Wide Research Programme for bluefin tuna (GBYP). Coordination Detailed Activity Report For Phase 2. Collect. Vol. Sci. Pap. ICCAT 68: 176-207), but the lack of direct studies in these areas prevents us from drawing conclusions on the subject for the time being. However, a recent identification of ABFT larvae in the Cantabrian Sea (Rodríguez et al., 2019Rodriguez J.M., Johnstone C. Lozano-Peral D. 2019. Is the Southern Bay of Biscay a spawning ground for Thunnus thynnus and Sarda sarda? ICCAT SCRS/2019/201) is an additional piece of evidence of the high degree of complexity of the ABFT spawning structure.

The complexity of ABFT population structure is further exacerbated when the Mediterranean component is considered; this has been a subject of discussion for years, as it could have important management implications (Viñas et al. 2011Viñas J., Gordoa A., Fernández-Cebrián R., et al. 2011. Facts and uncertainties about the genetic population structure of Atlantic bluefin tuna (Thunnus thynnus) in the Mediterranean. Implications for fishery management. Rev. Fish Biol. Fish. 21: 527-541. https://doi.org/10.1007/s11160-010-9174-6 , Arrizabalaga et al. 2019Arrizabalaga H., Arregui I., Medina A., et al. 2019. Life History and Migrations of Mediterranean Bluefin Tuna. In: Block B.A. (ed), The Future of Bluefin Tunas: Ecology, Fisheries Management, and Conservation, pp 67-93). The potential existence of discrete bluefin tuna subpopulations or contingents within the Mediterranean has been widely discussed (Renzoni et al. 1978Renzoni A., Bernhard M., Sara R., et al. 1978. Comparison between the Hg body burden of Thunnus thynnus from the Mediterranean and the Atlantic. IV Journées Etud. Pollutions, Morales-Nin and Fortuño 1990Morales-Nin B., Fortuño J.M. 1990. Mercury body burden and otolith characteristics of bluefin tuna from the northwest Mediterranean (Balearic Sea). Sci. Mar. 54: 277-285, Fromentin and Powers 2005Fromentin J.M., Powers J.E. 2005. Atlantic bluefin tuna: population dynamics, ecology, fisheries and management. Fish Fish 6: 281-306. https://doi.org/10.1111/j.1467-2979.2005.00197.x ), because the presence of bluefin tuna in the Mediterranean is not seasonal as in the Gulf of Mexico but permanent throughout the year. Genetic studies on bluefin tuna within the Mediterranean have shown different views on the existence of genetic subpopulations. Some studies have found regional differences within the Mediterranean (Viñas et al. 2003Viñas J., Pla C., Tawil M., et al. 2003. Mitochondrial genetic characterization of bluefin tuna (Thunnus thynnus) from three Mediterranean (Libya, Malta, Tunisia); and one Atlantic locations (Gulf of Cadiz). ICCAT Col Vol Sci Pap 55: 1282-1288, Carlsson et al. 2004Carlsson J., McDowell J., Díaz-Jaimes P., et al. 2004. Microsatellite and mitochondrial DNA analyses of Atlantic bluefin tuna (Thunnus thynnus thynnus) population structure in the Mediterranean Sea. Mol. Ecol. 13: 3345-3356. https://doi.org/10.1111/j.1365-294X.2004.02336.x , Riccioni et al. 2010Riccioni G., Landi M., Ferrara G., et al. 2010. Spatio-temporal population structuring and genetic diversity retention in depleted Atlantic bluefin tuna of the Mediterranean Sea. Proc. Nat. Acad. Sci. 107: 2102-2107. https://doi.org/10.1073/pnas.0908281107 ), but the most recent ones based on Single nucleotide polymorphisms and using larvae and age 0 fish as reference samples did not detect any genetic structure (Rodríguez-Ezpeleta et al. 2019Rodríguez-Ezpeleta N., Díaz-Arce N., Walter J.F., et al. 2019. Determining natal origin for improved management of Atlantic bluefin tuna. Front. Ecol. Environ. 17: 439-444 https://doi.org/10.1002/fee.2090 ). In addition, e-tagging studies have identified both resident individuals (Fromentin and Lopuszanski 2013Fromentin J.-M., Lopuszanski D. 2013. Migration, residency, and homing of bluefin tuna in the western Mediterranean Sea. ICES J. Mar. Sci. 71: 510-518. https://doi.org/10.1093/icesjms/fst157 , Cermeño et al. 2015Cermeño P., Quílez-Badia G., Ospina-Alvarez A., et al. 2015. Electronic tagging of Atlantic bluefin tuna (Thunnus thynnus, L.) reveals habitat use and behaviors in the Mediterranean Sea. PLoS ONE 10: e0116638. https://doi.org/10.1371/journal.pone.0116638 ) and fish that just use the Mediterranean to reproduce (e.g. Aranda et al. 2013Aranda G., Medina A., Santos A., et al. 2013. Evaluation of Atlantic bluefin tuna reproductive potential in the western Mediterranean Sea. J. Sea Res. 76: 154-160. https://doi.org/10.1016/j.seares.2012.08.005 , Abascal et al. 2016Abascal F.J., Medina A., De la Serna J.M., et al. 2016. Tracking bluefin tuna reproductive migration into the Mediterranean Sea with electronic pop-up satellite archival tags using two tagging procedures. Fish. Oceanogr. 25: 54-66. https://doi.org/10.1111/fog.12134 ). A recent review suggests that more populations or contingents might exist than was previously thought, and it seems more likely that ABFT that originated in the western Mediterranean migrate to the Atlantic for feeding more intensively than those that originated in the eastern Mediterranean (Arrizabalaga et al. 2019Arrizabalaga H., Arregui I., Medina A., et al. 2019. Life History and Migrations of Mediterranean Bluefin Tuna. In: Block B.A. (ed), The Future of Bluefin Tunas: Ecology, Fisheries Management, and Conservation, pp 67-93).

Additional information also reveals differences between Mediterranean regions. The spawning grounds in the Mediterranean, determined through histological analysis and larval findings, are well known: the Balearic waters in the western Mediterranean (Medina et al. 2002Medina A., Abascal F., Megina C., et al. 2002. Stereological assessment of the reproductive status of female Atlantic northern bluefin tuna during migration to Mediterranean spawning grounds through the Strait of Gibraltar. J. Fish. Biol. 60: 203-217. https://doi.org/10.1111/j.1095-8649.2002.tb02398.x ), the Levantine Sea in the eastern Mediterranean (Karakulak et al. 2004Karakulak S., Oray I., Corriero A., et al. 2004. Evidence of a spawning area for the bluefin tuna (Thunnus thynnus L.) in the eastern Mediterranean. J. Appl. Ichthyol. 20: 318-320. https://doi.org/10.1111/j.1439-0426.2004.00561.x ), and the South Tyrrhenian Sea and the waters around Malta and off the Tunisian coast and off the eastern coast of Sicily in the central Mediterranean (e.g. Nishida et al. 1998Nishida T., Tsuji S., Segawa K. 1998. Spatial data analyses of Atlantic bluefin tuna larval surveys in the 1994 ICCAT BYP. ICCAT Col. Vol. Sci. Pap. 48: 107-110, Corriero et al. 2003Corriero A., Desantis S., Deflorio M., et al. 2003. Histological investigation on the ovarian cycle of the bluefin tuna in the western and central Mediterranean. J. Fish. Biol. 63: 108-119. https://doi.org/10.1046/j.1095-8649.2003.00132.x , Giovanardi and Romanelli 2010Giovanardi O., Romanelli M. 2010. Preliminary note on tuna larvae in samples from the coasts of the southern-central Mediterranean Sea collected by the MV Arctic Sunrise in June/July 2008. ICCAT Col Vol Sci Pap 65: 740-743). The spawning migration path from the Atlantic to the Mediterranean has shown a strong connection with all Mediterranean spawning grounds except the Levantine Sea (Block et al. 2005Block B.A., Teo S.L., Walli A., et al. 2005. Electronic tagging and population structure of Atlantic bluefin tuna. Nature 434: 1121-1127. https://doi.org/10.1038/nature03463 , Walli et al. 2009Walli A., Teo S.L., Boustany A., et al. 2009. Seasonal movements, aggregations and diving behavior of Atlantic bluefin tuna (Thunnus thynnus) revealed with archival tags. PLoS ONE 4: e6151. https://doi.org/10.1371/journal.pone.0006151 , Arrizabalaga et al. 2019Arrizabalaga H., Arregui I., Medina A., et al. 2019. Life History and Migrations of Mediterranean Bluefin Tuna. In: Block B.A. (ed), The Future of Bluefin Tunas: Ecology, Fisheries Management, and Conservation, pp 67-93), which places the origin of the eastern spawners into question. In addition, several studies have reported differences in the fish size structure between Mediterranean spawning regions (Karakulak et al. 2004Karakulak S., Oray I., Corriero A., et al. 2004. Evidence of a spawning area for the bluefin tuna (Thunnus thynnus L.) in the eastern Mediterranean. J. Appl. Ichthyol. 20: 318-320. https://doi.org/10.1111/j.1439-0426.2004.00561.x , Heinisch et al. 2008Heinisch G., Corriero A., Medina A., et al., 2008. Spatial-temporal pattern of bluefin tuna (Thunnus thynnus L. 1758) gonad maturation across the Mediterranean Sea. Mar. Biol. 154: 623-630. https://doi.org/10.1007/s00227-008-0955-6 , Addis et al. 2016Addis P., Secci M., Biancacci C., et al. 2016. Reproductive status of Atlantic bluefin tuna, Thunnus thynnus, during migration off the coast of Sardinia (western Mediterranean). Fish. Res. 181: 137-147. https://doi.org/10.1016/j.fishres.2016.04.009 ). The spawners in the eastern Mediterranean seem to be smaller than those entering the Mediterranean and those from the western Mediterranean spawning grounds, while sizes in the central region had a wider distribution, which might be representative of a mixing spawning ground of fish coming from the Atlantic and eastern Mediterranean (Heinisch et al. 2008Heinisch G., Corriero A., Medina A., et al., 2008. Spatial-temporal pattern of bluefin tuna (Thunnus thynnus L. 1758) gonad maturation across the Mediterranean Sea. Mar. Biol. 154: 623-630. https://doi.org/10.1007/s00227-008-0955-6 ).

Differences in spawning periods among Mediterranean regions are known: spawning takes place in June-July in the western basin (e.g. Susca et al. 2001Susca V., Corriero A., Deflorio M., et al. 2001. New results on the reproductive biology of the bluefin tuna (Thunnus thynnus) in the Mediterranean. ICCAT Col Vol Sci Pap 52: 745-751, Medina et al. 2002Medina A., Abascal F., Megina C., et al. 2002. Stereological assessment of the reproductive status of female Atlantic northern bluefin tuna during migration to Mediterranean spawning grounds through the Strait of Gibraltar. J. Fish. Biol. 60: 203-217. https://doi.org/10.1111/j.1095-8649.2002.tb02398.x , Corriero et al. 2003Corriero A., Desantis S., Deflorio M., et al. 2003. Histological investigation on the ovarian cycle of the bluefin tuna in the western and central Mediterranean. J. Fish. Biol. 63: 108-119. https://doi.org/10.1046/j.1095-8649.2003.00132.x ) and one month earlier, May-June, in the eastern basin (Duclerc et al. 1974Duclerc J., Piccinetti C., Piccinetti Manfrin G., et al. 1974. Nouvelles donne sur la reproduction du thon rouge et d’autres espèces de thonidés en Méditerranée. Rapp. Comm. Int. Mer Médit 22: 37, Karakulak et al. 2004Karakulak S., Oray I., Corriero A., et al. 2004. Evidence of a spawning area for the bluefin tuna (Thunnus thynnus L.) in the eastern Mediterranean. J. Appl. Ichthyol. 20: 318-320. https://doi.org/10.1111/j.1439-0426.2004.00561.x , Oray and Karakulak 2005Oray I., Karakulak F. 2005. Further evidence of spawning of bluefin tuna (Thunnus thynnus L., 1758) and the tuna species (Auxis rochei Ris., 1810, Euthynnus alletteratus Raf., 1810) in the eastern Mediterranean Sea: preliminary results of TUNALEV larval survey in 2004. J. Appl. Ichthyol. 21: 236-240. https://doi.org/10.1111/j.1439-0426.2005.00658.x ). These differences are attributed to temporal differences in water warming between regions. In the Mediterranean Sea the differences in temperature and chlorophyll between the eastern and western basins are permanent throughout the year (d’Ortenzio and Ribera d’Alcalà 2009d’Ortenzio F., Ribera d’Alcalà M. 2009. On the trophic regimes of the Mediterranean Sea: a satellite analysis. Biogeosciences 6: 139-148. https://doi.org/10.5194/bg-6-139-2009 , Skliris et al. 2012Skliris N., Sofianos S., Gkanasos A., et al. 2012. Decadal scale variability of sea surface temperature in the Mediterranean Sea in relation to atmospheric variability. Ocean Dyn. 62: 13-30. https://doi.org/10.1007/s10236-011-0493-5 , Shaltout and Omstedt 2014Shaltout M., Omstedt A. 2014. Recent sea surface temperature trends and future scenarios for the Mediterranean Sea. Oceanologia 56: 411-443. https://doi.org/10.5697/oc.56-3.411 ). Temperature and food availability can affect growth and survival during the early stages of fish development (Anderson 1988Anderson J.T. 1988. A review of size dependent survival during pre-recruit stages of fishes in relation to recruitment. J. Northwest Atl. Fish. Sci. 8: 55-66. https://doi.org/10.2960/J.v8.a6 ), as evidenced in Pacific and Atlantic bluefin tuna larvae (García et al. 2013García A., Cortés D., Quintanilla J., et al. 2013. Climate-induced environmental conditions influencing interannual variability of Mediterranean bluefin (Thunnus thynnus) larval growth. Fish. Oceanogr. 22: 273-287. https://doi.org/10.1111/fog.12021 , Satoh et al. 2013Satoh K., Tanaka Y., Masujima M., et al. 2013. Relationship between the growth and survival of larval Pacific bluefin tuna, Thunnus orientalis. Mar. Biol. 160: 691-702. https://doi.org/10.1007/s00227-012-2124-1 , Ishihara et al. 2019Ishihara T., Watai M., Ohshimo S., et al. 2019. Differences in larval growth of Pacific bluefin tuna (Thunnus orientalis) between two spawning areas, and an evaluation of the growth-dependent mortality hypothesis. Environ Biol Fishes 102: 581-594. https://doi.org/10.1007/s10641-019-00855-w ). Consequently, the oceanographic differences within the Mediterranean may cause differences in the growth rate during the earlier stages of development of the species. In spite of this, no comparative studies on the larval or young-of-the-year (YOY) growth within the Mediterranean have yet been carried out. So far, studies on ABFT larval growth are limited to the larvae collected in the western spawning ground in the Balearic waters (García et al. 2006García A., Cortés D., Ramírez T., et al. 2006. First data on growth and nucleic acid and protein content of field-captured Mediterranean bluefin (Thunnus thynnus) and albacore (Thunnus alalunga) tuna larvae: a comparative study. Sci. Mar. 70: 67-78. https://doi.org/10.3989/scimar.2006.70s267 , 2013García A., Cortés D., Quintanilla J., et al. 2013. Climate-induced environmental conditions influencing interannual variability of Mediterranean bluefin (Thunnus thynnus) larval growth. Fish. Oceanogr. 22: 273-287. https://doi.org/10.1111/fog.12021 ), whereas the studies on YOY are performed with juveniles caught in the central Mediterranean (Santamaria et al. 2009Santamaria N., Bello G., Corriero A., et al. 2009. Age and growth of Atlantic bluefin tuna, Thunnus thynnus (Osteichthyes: Thunnidae), in the Mediterranean Sea. J. Appl. Ichthyol. 25: 38-45. https://doi.org/10.1111/j.1439-0426.2009.01191.x La Mesa et al. 2005La Mesa M., La Mesa G., Micalizzi M. 2005. Age and growth of madeira scorpionfish, Scorpaena maderensis Valenciennes, 1833, in the central Mediterranean. Fish. Res. 74: 265-272. https://doi.org/10.1016/j.fishres.2005.01.018 ).

The objectives of the present study were to analyse the growth rates of ABFT YOY in the Mediterranean basin in 2013 and 2016 and investigate potential differences in growth rates between eastern, central and western regions. We analysed the microstructure of ABFT otoliths using the daily increment counts to estimate the age and growth of juveniles collected in the framework of the Grand Bluefin Tuna Year Programme (GBYP). The GBYP was officially adopted by the ICCAT Commission in 2008 and aims to improve basic data collection through data mining, understanding of key biological and ecological processes and assessment models, and provision of scientific advice on stock status.

MATERIALS AND METHODS

Sampling

 ⌅

The otoliths of ABFT YOY were taken from the GBYP biological data bank; they were sampled within the GBYP biological sampling programme in several Mediterranean regions. During the sampling, catch date was recorded and fish fork length (FL) was measured to the nearest 0.1 cm. Each sagittal otolith was carefully extracted and cleaned. After drying, the otoliths were stored in plastic vials and kept in the GBYP data bank. The samples were requested to the GBYP for two years, 2013 and 2016, and were selected according to the specific objectives for each year. In 2013, the objective was to analyse the growth rate of YOY in the whole Mediterranean basin; a total of 60 otoliths were selected, 20 per region (eastern, western and central) and from samples caught close together in time (Fig. 1). In 2016, the objective was to analyse the growth rates for each Mediterranean region, which required a wide range of sizes in each region; consequently, for each region the otoliths were selected from two subsets of samples separated in time. The length frequency distribution of each sample, recorded in the GBYP data bank, was examined to ensure that the otoliths selected for each location represented the whole size range of the fish caught (Table 1). The juveniles collected in the Tyrrhenian Sea in August 2016 displayed a wider size range, and consequently the sample size was larger than that of those from the other regions. Another particularity of 2016 was that only 27 juvenile individuals were collected in the eastern region, so the analysed subsample (20) represented practically all the fish available in the GBYP data bank for that year and region.

Fig. 1.  Sampling sites. Western locations (diamonds), central locations (squares) and eastern locations (circles). Symbols filled with grey and black outlines represent the years 2013 and 2016, respectively.

Table 1.  Summary of collected data: dates, locations, number of individuals and size ranges.

Year	Area	Location	Catch date	Sample size	Fork length (cm)
2013	Western	Balearic Islands	19 October	20	32.5 - 41.0
	Central	Malta	10-23 September	20	19.4 - 27.6
	Eastern	North Cyprus	15-20 August	20	21.0 - 26.4
2016	Western	Balearic Islands	10-30 September	10	20.4 - 33.1
		Balearic Islands	2-18 November	11	35.0 - 48.0
	Central	South Tyrrhenian	19-29 August	20	18.4 - 52.6
		South Tyrrhenian	12 October	2	29.5 - 32.6
		South Tyrrhenian	4-17 December	13	43.0 - 57.0
	Eastern	North Cyprus	22-29 July	5	14.7 - 17.1
		North Cyprus	13-31 August	5	16.1 - 30.4
		North Cyprus	2-10 September	10	20.5 - 35.8

RESULTS

The otolith microstructure showed a concentric pattern of increments from the core region until around the sixth increment, and then became larger and elongated along the postero-anterior axis, which was gradually reduced (Fig. 2). In 2013, 83.3% of the readings differed from each other by less than 5%, whereas 16.7% of the readings differed by 5% to 10% and no age estimate differed more than 10%. In 2016, 87% of the readings differed by less than 5% and the differences between the remaining percentage were below 9%.

A total of 134 specimens between 14.7 and 57 cm FL were aged. This is the largest sample size on YOY bluefin tuna otolith microstructure analysis so far. Age estimates ranged from 45 to 192 days. An overview of the results (Fig. 3 A) showed some abnormal observations: young individuals (<80 days) with extraordinarily large sizes, even above 50 cm. These observations corresponded to YOY caught in the Tyrrhenian Sea in August 2016. The catch was comprised of individuals with a wide length range, spanning 18.4 to 52.6 cm FL, but with a narrow age range (59 to 78 days). The age and length of the exceptional observations were double-checked. First, the age was re-read and the new readings confirmed the previous age estimates. Next, the reported fork lengths of the YOY captured in 2016 were examined with the size of their otoliths based on the known linear relationship between fish length (FL) and otolith length (e.g. Jenkins and Davis 1990Jenkins G.P., Davis T.L. 1990. Age, growth rate, and growth trajectory determined from otolith microstructure of southern bluefin tuna Thunnus maccoyii larvae. Mar. Ecol. Prog. Ser. 63: 93-104. https://doi.org/10.3354/meps063093 , La Mesa et al. 2005La Mesa M., La Mesa G., Micalizzi M. 2005. Age and growth of madeira scorpionfish, Scorpaena maderensis Valenciennes, 1833, in the central Mediterranean. Fish. Res. 74: 265-272. https://doi.org/10.1016/j.fishres.2005.01.018 , Gunn et al. 2008Gunn J.S., Clear N.P., Carter T.I., et al. 2008. Age and growth in southern bluefin tuna, Thunnus maccoyii (Castelnau): direct estimation from otoliths, scales and vertebrae. Fish. Res. 92: 207-220. https://doi.org/10.1016/j.fishres.2008.01.018 ). The radius length (LR) was measured for each otolith and the linear regression between FL and LR were estimated. The model fitted with the data holding the exceptional observations only explained 46% of FL variability but increased to 92% when the sample was excluded (Fig. 3 B). These results indicated that the extraordinarily large sizes in that specific sample were wrong. Consequently, every FL observation outside the 99% prediction interval was discarded from the growth analysis; a total of 13 fish were excluded from the growth analysis because their sizes were outside the prediction interval but their ages were included in the birth date estimations.

The birth dates found for each region varied from 24 May to 11 July in the eastern region, from 31 May to 4 August in the central region and from 2 June to 29 July in the western region. The extension of the birth intervals in the eastern and western regions was close to two months (Fig. 4). The main differences between regions (Fig. 4) were: 1) the eastern region consistently showed earlier ending and earlier hatching, with the peak around the first fortnight of June for both sampled years, 2) the hatching interval was wider in the central region than in the other regions and showed two peaks, in mid-June and in mid-July, and 3) the peak observed in the western region in late June took place between the two peaks observed in the central region.

The growth models fitted for 2013 and 2016 observations explained more than 90% of the daily growth variability of YOY in the whole Mediterranean basin (Fig. 5). The t-test analysis showed no significant differences between the growth rates of 2013 (3.25 mm d^-1) and 2016 (2.70 mm d^-1). In contrast, the results of the growth models fitted for each region separately displayed high variability, from 2.52 mm d^-1 to 4.01 mm d^-1 (Table 2). While the western and central regions showed no significant differences between each other, the daily growth rate of the eastern region (4.01 mm d^-1) was significantly faster. These differences in the growth rates between basins remain after the additional analysis in the western and central basin, in which the largest fish (>380 mm) were excluded (Table 2).

DISCUSSION

Spawning onset of T. thynnus seems to be temperature-related for all tuna species (Schaefer 2001Schaefer K.M. 2001. Reproductive biology of tunas. Fish Physiol.19: 225-270. https://doi.org/10.1016/S1546-5098(01)19007-2 ). In the Mediterranean, the reported temperatures during the ABFT spawning period range from 19.5°C to 26.5°C (García et al. 2005García A., Alemany F., De la Serna J., et al. 2005. Preliminary results of the 2004 bluefin tuna larval surveys off different Mediterranean sites (Balearic Archipelago, Levantine Sea and the Sicilian Channel). ICCAT Col. Vol. Sci. Pap. 58: 1261-1270, Teo et al. 2007Teo S.L., Boustany A.M., Block B.A. 2007. Oceanographic preferences of Atlantic bluefin tuna, Thunnus thynnus, on their Gulf of Mexico breeding grounds. Mar. Biol. 152: 1105-1119. https://doi.org/10.1007/s00227-007-0758-1 , Gordoa and Carreras 2014Gordoa A., Carreras G. 2014. Determination of Temporal Spawning Patterns and Hatching Time in Response to Temperature of Atlantic Bluefin Tuna (Thunnus thynnus) in the Western Mediterranean. PLoS ONE 9: e90691. https://doi.org/10.1371/journal.pone.0090691 ), but neither spawning nor water warming are simultaneous along its almost 4000 km of length. The time differences in spawning between regions assumed to be related to the warming time lag range from June to July in the western basin (e.g. Susca et al. 2001Susca V., Corriero A., Deflorio M., et al. 2001. New results on the reproductive biology of the bluefin tuna (Thunnus thynnus) in the Mediterranean. ICCAT Col Vol Sci Pap 52: 745-751, Medina et al. 2002Medina A., Abascal F., Megina C., et al. 2002. Stereological assessment of the reproductive status of female Atlantic northern bluefin tuna during migration to Mediterranean spawning grounds through the Strait of Gibraltar. J. Fish. Biol. 60: 203-217. https://doi.org/10.1111/j.1095-8649.2002.tb02398.x , Corriero et al. 2003Corriero A., Desantis S., Deflorio M., et al. 2003. Histological investigation on the ovarian cycle of the bluefin tuna in the western and central Mediterranean. J. Fish. Biol. 63: 108-119. https://doi.org/10.1046/j.1095-8649.2003.00132.x ) and from May to June in the eastern basin (Duclerc et al. 1974Duclerc J., Piccinetti C., Piccinetti Manfrin G., et al. 1974. Nouvelles donne sur la reproduction du thon rouge et d’autres espèces de thonidés en Méditerranée. Rapp. Comm. Int. Mer Médit 22: 37, Karakulak et al. 2004Karakulak S., Oray I., Corriero A., et al. 2004. Evidence of a spawning area for the bluefin tuna (Thunnus thynnus L.) in the eastern Mediterranean. J. Appl. Ichthyol. 20: 318-320. https://doi.org/10.1111/j.1439-0426.2004.00561.x , Oray and Karakulak 2005Oray I., Karakulak F. 2005. Further evidence of spawning of bluefin tuna (Thunnus thynnus L., 1758) and the tuna species (Auxis rochei Ris., 1810, Euthynnus alletteratus Raf., 1810) in the eastern Mediterranean Sea: preliminary results of TUNALEV larval survey in 2004. J. Appl. Ichthyol. 21: 236-240. https://doi.org/10.1111/j.1439-0426.2005.00658.x ). These spawning windows in the Mediterranean spawning regions are to a certain extent in accordance with the hatching time intervals shown in this study. However, they do not support the hypothesis of a 2016 winter spawning event based on the anomalously large YOY caught in the Tyrrhenian sea (Di Natale et al. 2017Di Natale A., Tensek S., Celona A., et al. 2017. Another peculiar situation for YOY of bluefin tuna (Thunnus thynnus) in the Mediterranean Sea in 2016. ICCAT Collect. Vol. Sci. Pap. 74: 2839-2849), which have been shown here to be errors.

It should be highlighted that the estimated hatching time intervals shown here may be only indicative of the minimum spawning window expected in each region. Notwithstanding this limitation, the observed hatching interval in the central Mediterranean, which spanned from the end of May to the beginning of August, was wider than the spawning window reported by previous studies. The sampling periods reported by these studies, either based on gonadal development or larval survey, were limited to May-July (Corriero et al. 2003Corriero A., Desantis S., Deflorio M., et al. 2003. Histological investigation on the ovarian cycle of the bluefin tuna in the western and central Mediterranean. J. Fish. Biol. 63: 108-119. https://doi.org/10.1046/j.1095-8649.2003.00132.x , García et al. 2005García A., Alemany F., De la Serna J., et al. 2005. Preliminary results of the 2004 bluefin tuna larval surveys off different Mediterranean sites (Balearic Archipelago, Levantine Sea and the Sicilian Channel). ICCAT Col. Vol. Sci. Pap. 58: 1261-1270, Heinisch et al. 2008Heinisch G., Corriero A., Medina A., et al., 2008. Spatial-temporal pattern of bluefin tuna (Thunnus thynnus L. 1758) gonad maturation across the Mediterranean Sea. Mar. Biol. 154: 623-630. https://doi.org/10.1007/s00227-008-0955-6 ). We must go back to 1932 to find the first and last indication of bluefin tuna spawning in August, when some larvae were found around Sicilian waters (Sanzo 1932Sanzo L. 1932. Uova e primi stadi larvali di tonno (Orcynus thynnus Ltkn.). Mem. R. Com. Talass. Ital. 189: 1-16. https://doi.org/10.1080/11250003209429197 ). Currently, there is no information to infer whether spawning events in August are permanent or occasional. Similarly, in the eastern region (Levantine Sea), the hatching interval lasted until mid-July, which is well beyond the previously reported spawning period (May-June).

The results of the central region covered a wide hatching period (late May to early August) with a bimodality in the hatching frequency in mid-June and mid-July. The observed pattern could be an artefact of a weak sampling coverage. However, this is unlikely given that the expected peak under a progressive spawning process would be found in late June and/or early July, which is precisely when we observed the minimum frequency. Furthermore, these results are consistent with the length frequency bimodality of YOY found by Relini et al. (1995)Relini M., G. Palandri, Torchia G. 1995. Tagging of Thunnus thynnus juveniles in the Ligurian Sea, Autumn 1994. ICCAT Col Vol Sci Pap ICCAT 44: 378, who suggested two different spawning pulses in the central region. The existence of different spawning pulses might be due to the existence of different bluefin tuna contingents (e.g. residents and Atlantic migrants) using the central spawning region at slightly different times, but further research is needed to support this hypothesis.

The sampled location of the YOY does not necessarily correspond to their hatching place. However, based on current knowledge, it is likely that sampling of YOYs occurred close to their hatching places in the western and eastern regions. The Mediterranean is comprised of two basins that are connected by the strait of Sicily, but each of them has a counter-clockwise surface current (Millot 1999Millot C. 1999. Circulation in the western Mediterranean Sea. J. Mar. Syst. 20: 423-442. https://doi.org/10.1016/S0924-7963(98)00078-5 , Hamad et al. 2005Hamad N., Millot C., Taupier-Letage I. 2005. A new hypothesis about the surface circulation in the eastern basin of the Mediterranean Sea. Prog. Oceanogr. 66: 287-298. https://doi.org/10.1016/j.pocean.2005.04.002 ) which, in addition to the long distance between the western and eastern spawning regions, makes the exchange of YOY between western and eastern regions unlikely. The central region represents a different scenario due to its large extension, from Sicily to the Aegean Sea, which could receive some YOY from the Levantine Sea. Unfortunately, little is currently known about YOY migration capability because of the lack of scientific tagging studies for this age class.

The few studies on age and growth of YOY in the Mediterranean (Table 2) showed much faster growth rates than those reported for YOY in the western Atlantic (Brothers et al. 1983Brothers E.B., Prince E.D., Lee D.W. 1983. Age and growth of young-of-the-year bluefin tuna, Thunnus thynnus, from otolith microstructure. NOAA Tech. Report NMFS 8: 49-59, Arai et al. 2020Arai K., Graves J.E., Secor D.H.. 2020. Sub-annual cohort representation among young-of-the-year recruits of the western stock of Atlantic bluefin tuna. Fish. Res. 225: 105476. https://doi.org/10.1016/j.fishres.2019.105476 ), which were around 1.1 to 1.5 mm d^-1. The results of the growth model for the central Mediterranean showed no major differences from previous estimates in this region (La Mesa et al. 2005La Mesa M., La Mesa G., Micalizzi M. 2005. Age and growth of madeira scorpionfish, Scorpaena maderensis Valenciennes, 1833, in the central Mediterranean. Fish. Res. 74: 265-272. https://doi.org/10.1016/j.fishres.2005.01.018 , Santamaria et al. 2009Santamaria N., Bello G., Corriero A., et al. 2009. Age and growth of Atlantic bluefin tuna, Thunnus thynnus (Osteichthyes: Thunnidae), in the Mediterranean Sea. J. Appl. Ichthyol. 25: 38-45. https://doi.org/10.1111/j.1439-0426.2009.01191.x ), with the exception of the data reported by Megalofonou (2006)Megalofonou P. 2006. Comparison of otolith growth and morphology with somatic growth and age in young-of-the-year bluefin tuna. J. Fish. Biol. 68: 1867-1878. https://doi.org/10.1111/j.1095-8649.2006.01078.x , which differed from the rest. In the latter study, the model showed faster growth: even faster than our estimation for the eastern region. The YOY collected in that study came from a very large area, which included various seas of the central Mediterranean: the south Thyrrhenian, the Ionian, the Adriatic and also the Aegean Sea. However, it is unlikely all these fish have their origin in the central Mediterranean spawning ground. In particular, the YOY collected in the Aegean Sea may have come from the Levantine Sea because this is the closest spawning ground and the anti-clockwise direction of the surface current from the Levantine to the Aegean Sea would favour the transport of these specimens.

The present study has dealt with YOY information, and, although only individuals surviving the recruitment phase may be analysed, one might assume that they are the best adapted for each scenario. There is evidence for the Pacific bluefin tuna that offspring survival depends largely on the growth rates during the late larval phase (Tanaka et al. 2006Tanaka Y., Satoh K., Iwahashi M., et al. 2006. Growth-dependent recruitment of Pacific bluefin tuna Thunnus orientalis in the northwestern Pacific Ocean. Mar. Ecol. Prog. Ser. 319: 225-235. https://doi.org/10.3354/meps319225 , Watai et al. 2017Watai M., Ishihara T., Abe O., et al. 2017. Evaluation of growth-dependent survival during early stages of Pacific bluefin tuna using otolith microstructure analysis. Mar. Freshw. Res. 68: 2008-2017. https://doi.org/10.1071/MF16337 ). Food availability and water temperature are the critical factors affecting the growth (Heath 1992Heath M. 1992. Field investigations of the early life stages of marine fish. Pages 1-174 Advances in Mar. Biol. Elsevier. https://doi.org/10.1016/S0065-2881(08)60039-5 ) and survival rates (Campana 1996Campana S.E. 1996. Year-class strength and growth rate in young Atlantic cod Gadus morhua. Mar. Ecol. Prog. Ser. 135: 21-26. https://doi.org/10.3354/meps135021 , Meekan and Fortier 1996Meekan M.G., Fortier L. 1996. Selection for fast growth during the larval life of Atlantic cod Gadus morhua on the Scotian Shelf. Mar. Ecol. Prog. Ser. 137: 25-37. https://doi.org/10.3354/meps137025 ) of marine fish larvae. For bluefin tuna, high water temperature seems to promote larval growth rates (García et al. 2013García A., Cortés D., Quintanilla J., et al. 2013. Climate-induced environmental conditions influencing interannual variability of Mediterranean bluefin (Thunnus thynnus) larval growth. Fish. Oceanogr. 22: 273-287. https://doi.org/10.1111/fog.12021 , Satoh et al. 2013Satoh K., Tanaka Y., Masujima M., et al. 2013. Relationship between the growth and survival of larval Pacific bluefin tuna, Thunnus orientalis. Mar. Biol. 160: 691-702. https://doi.org/10.1007/s00227-012-2124-1 ). However, the results of a recent study on Pacific bluefin tuna pointed out that both temperature and food availability have an effect on larval growth but at different stages of development (Ishihara et al. 2019Ishihara T., Watai M., Ohshimo S., et al. 2019. Differences in larval growth of Pacific bluefin tuna (Thunnus orientalis) between two spawning areas, and an evaluation of the growth-dependent mortality hypothesis. Environ Biol Fishes 102: 581-594. https://doi.org/10.1007/s10641-019-00855-w ).

The findings of the current study showed a significantly faster YOY growth in the eastern Mediterranean, and the results ruled out the possibility that this was due to the size range of the YOY in this region. The warmer temperature in this region compared with the western region (Shaltout and Omstedt 2014Shaltout M., Omstedt A. 2014. Recent sea surface temperature trends and future scenarios for the Mediterranean Sea. Oceanologia 56: 411-443. https://doi.org/10.5697/oc.56-3.411 ) could be an important factor determining the observed differences in growth. It is worth remembering that the results of growth rates are skewed toward higher values because only surviving individuals are analysed (Le Pape and Bonhommeau 2015Le Pape O., Bonhommeau S. 2015. The food limitation hypothesis for juvenile marine fish. Fish Fish 16: 373-398. https://doi.org/10.1111/faf.12063 ), because key survival factors might change between regions.

The population structure of ABFT has been the object of debate for decades, particularly regarding the eastern Mediterranean population. The presence of fish in the Mediterranean all year around (Cermeño et al. 2015Cermeño P., Quílez-Badia G., Ospina-Alvarez A., et al. 2015. Electronic tagging of Atlantic bluefin tuna (Thunnus thynnus, L.) reveals habitat use and behaviors in the Mediterranean Sea. PLoS ONE 10: e0116638. https://doi.org/10.1371/journal.pone.0116638 ), along with the potentially smaller connectivity between the Atlantic and the eastern Mediterranean (Walli et al. 2009Walli A., Teo S.L., Boustany A., et al. 2009. Seasonal movements, aggregations and diving behavior of Atlantic bluefin tuna (Thunnus thynnus) revealed with archival tags. PLoS ONE 4: e6151. https://doi.org/10.1371/journal.pone.0006151 , Arrizabalaga et al. 2019Arrizabalaga H., Arregui I., Medina A., et al. 2019. Life History and Migrations of Mediterranean Bluefin Tuna. In: Block B.A. (ed), The Future of Bluefin Tunas: Ecology, Fisheries Management, and Conservation, pp 67-93), raises the question of the existence of a single eastern Mediterranean population or contingent. Thus, the extent to which the faster growth in the eastern Mediterranean can be attributed just to the warmer environmental conditions, excluding biological potential differences, should be further investigated.

In synthesis, we can conclude that the growth rate of ABFT YOY is faster in the eastern Mediterranean, possibly because of the warmer conditions of the region, but ontogenetic differences cannot be excluded. We suggest the possibility of occurrence of two spawning pulses in the central Mediterranean, which might correspond to the existence of two different bluefin tuna contingents (e.g. residents and Atlantic migrants) using the central spawning region at slightly different times. However, further research is needed to support this hypothesis.