Predicting the occurrence of Atlantic bluefin tuna (Thunnus thynnus) larvae in the northern Gulf of Mexico: building a classification model from archival data

Although bluefin tuna are found throughout the Atlantic Ocean, spawning in the western Atlantic has been recorded predominantly in the Gulf of Mexico (GOM) in spring. Larval bluefin tuna abundances from the northern GOM are formulated into an index used to tune the adult stock assessment, and the variability of this index is currently high. This study investigated whether some of the variability in larval bluefin tuna abundances was related to environmental conditions, by defining associations between larval bluefin tuna catch locations, and a suite of environmental variables. We hypothesized that certain habitat types, as defined by environmental variables, would be more likely to contain bluefin tuna larvae. Favorable habitat for bluefin tuna larvae was defined using a classification tree approach. Habitat within the Loop Current was generally less favorable, as were warm‐core rings, and cooler waters on the continental shelf. The location and size of favorable habitat was highly variable among years, which was reflected in the locations of larval bluefin tuna catches. The model successfully placed bluefin tuna larvae in favorable habitat with nearly 90% accuracy, but many negative stations were also located within theoretically favorable habitat. The probability of collecting larval bluefin tuna in favorable habitat was nearly twice the probability of collecting bluefin tuna larvae across all habitats (35.5 versus 21.0%). This model is a useful addition to knowledge of larval bluefin tuna distributions; however, the incorporation of variables describing finer‐scale features, such as thermal fronts, may significantly improve the model’s predictive power.     

Whole-body heat transfer coefficient and body temperature change of juvenile Pacific bluefin tuna <Emphasis Type="Italic">Thunnus orientalis</Emphasis> according to growth

Bluefin tuna maintain a higher body temperature than ambient sea water. Body heat is derived mainly from metabolic heat to elevate and maintain regional body temperature that is higher than the ambient, while heat loss is caused by heat transfer throughout the whole body surface and gills. Retention of high body temperature is thought to differ at each growth stage, so that a larger body mass maintains a higher body temperature. We evaluated the whole-body heat transfer coefficient, thermal difference between each tissue and water temperature, and metabolic heat in tissues during swimming of juvenile bluefin tuna as a function of fork length (FL) using a small thermometer and a treadmill-type flow tank. A system for maintaining high body temperature was well developed in fish with FL greater than 20.0 cm. Whole-body heat transfer coefficient was fitted to a −0.695 power of mass. Juvenile bluefin tuna showed a transition speed of 3.0 FL/s at which they switched from aerobic to anaerobic motion.     

Inclusion of prey data improves prediction of bluefin tuna (Thunnus thynnus) distribution

We examined the distribution of Atlantic bluefin tuna ( Thunnus thynnus ) in the Gulf of Maine, Northwest Atlantic Ocean, from 17 to 23 August 1995, in relation to physical and biological parameters. Specifically, we fit a binomial GLM to the bluefin tuna presence–absence data and predictor variables that include: sea surface temperature (SST), ocean depth, distance to an SST front, time-lagged density of SST fronts, and an interpolated surface of Atlantic herring ( Clupea harengus ) density. In addition, we use simple and partial Mantel tests to examine whether bluefin tuna presence–absence data are significantly associated with these predictors, once spatial autocorrelation is accounted for. Results suggest that the distribution of bluefin tuna significantly correlated with herring density ( z  =   3.525, P  =   0.000424), and that inclusion of biological variables results in a more parsimonious model. Mantel tests results indicate that bluefin tuna abundance is significantly correlated with herring density after the effect of spatial structure is removed (Mantel r  =   0.043, P  <   0.019).     

Horizontal and vertical movements of juvenile bluefin tuna (Thunnus orientalis) in relation to seasons and oceanographic conditions in the eastern Pacific Ocean

Electronically tagged juvenile Pacific bluefin, Thunnus orientalis, were released off Baja California in the summer of 2002. Time‐series data were analyzed for 18 fish that provided a record of 380 ± 120 days (mean ± SD) of ambient water and peritoneal cavity temperatures at 120 s intervals. Geolocations of tagged fish were estimated based on light‐based longitude and sea surface temperature‐based latitude algorithms. The horizontal and vertical movement patterns of Pacific bluefin were examined in relation to oceanographic conditions and the occurrence of feeding events inferred from thermal fluctuations in the peritoneal cavity. In summer, fish were located primarily in the Southern California Bight and over the continental shelf of Baja California, where juvenile Pacific bluefin use the top of the water column, undertaking occasional, brief forays to depths below the thermocline. In autumn, bluefin migrated north to the waters off the Central California coast when thermal fronts form as the result of weakened equatorward wind stress. An examination of ambient and peritoneal temperatures revealed that bluefin tuna fed during this period along the frontal boundaries. In mid‐winter, the bluefin returned to the Southern California Bight possibly because of strong downwelling and depletion of prey species off the Central California waters. The elevation of the mean peritoneal cavity temperature above the mean ambient water temperature increased as ambient water temperature decreased. The ability of juvenile bluefin tuna to maintain a thermal excess of 10°C occurred at ambient temperatures of 11–14°C when the fish were off the Central California coast. This suggests that the bluefin maintain peritoneal temperature by increasing heat conservation and possibly by increasing internal heat production when in cooler waters. For all of the Pacific bluefin tuna, there was a significant correlation between their mean nighttime depth and the visible disk area of the moon.     

Morphological features and functions of bluefin tuna change with growth

Morphological features of Pacific bluefin tuna Thunnus orientalis have important functions for its fast swimming. Morphological features of tuna change with growth; therefore, morphological functions may develop during this process. In this study, we precisely quantified the morphology of bluefin tuna with growth from juvenile to young adult using a three-dimensional laser profiler and evaluated the fluid dynamic characteristics of tuna using computational fluid dynamics (CFD) and an accurate model based on measured data. As results of measurement of morphological features, the aspect ratio and sweepback angle, which are indices of hydrodynamic characteristics for a hydrofoil, suggested that the lift force of caudal fin was increased as the tuna grows. The results of CFD analysis showed that the coefficient of drag force gradually decreased with growth. Pectoral fins generated lift force, and the ratio of lift force to submerged weight (FL/SW) increased as the tuna grew to 0.2 m total length (TL). After the tuna exceeded 0.2 m TL, FL/SW changed to a wider range in angle of attack as the tuna grew. These results suggest that the morphological function of bluefin tuna develops to enhance its swimming ability as it grows from juvenile to young adult.     

Environmental influences on Atlantic bluefin tuna (Thunnus thynnus) catch per unit effort in the southern Gulf of St. Lawrence

The Atlantic bluefin tuna (Thunnus thynnus) population in the western Atlantic supports substantial commercial and recreational fisheries. Despite quota establishment and management under the auspices of the International Commission for the Conservation of Atlantic Tunas, only small increases in population growth have been estimated. In contrast to other western bluefin tuna fisheries indices, contemporary estimates of catch per unit effort (CPUE) in the southern Gulf of St. Lawrence have increased rapidly and are at record highs. This area is characterized by the Cold Intermediate Layer (CIL) that is defined by waters <3°C and located at depths of 30–40 m in September. We investigated the influence of several in situ environmental variables on the bluefin tuna fishery CPUE using delta‐lognormal modelling and relatively extensive and consistent oceanographic survey data, as well as dockside monitoring and mandatory logbook data associated with the fishery. Although there is considerable spatial and temporal variation of water mass characteristics, the amount of available habitat in the southern Gulf of St. Lawrence (assuming a > 3°C thermal ambit) for bluefin tuna has been increasing. The percentage of the water column occupied by the CIL was a significant environmental variable in the standardization of CPUE estimates. There was also a negative relationship between the spatial extents of the CIL and the fishery. Properties of the CIL account for variation in the bluefin tuna CPUE and may be a factor in determining the amount of available feeding habitat for bluefin tuna in the southern Gulf of St. Lawrence.     

Why do young Pacific bluefin tuna repeatedly dive to depths through the thermocline?

Young Pacific bluefin tuna Thunnus orientalis with archival tags were released in the East China Sea. Time-series data for depth, and ambient and peritoneal temperatures for nine fish, recorded every 128 s, were analyzed. Our objectives were to describe monthly changes in diving patterns in relation to the ambient thermal structure and the occurrence of feeding events during March–June, and to discuss possible reasons why bluefin tuna repeatedly dived to depths below the thermocline in terms of their thermoconservation mechanisms. It was found that the fish repeatedly dived through the thermocline at intervals of 1.2 h on average, and the dive frequency was high during March–May. However, the dive frequency and periodicity decreased in June, when the gradient became steeper. In contrast, feeding events increased in June. These results indicate that from March to May, bluefin tuna repeatedly dive because food biomass is inadequate at the surface, and they stop undertaking repeated dives in June when food becomes more readily available at the surface, in addition to low visibility caused by low solar radiation. Further, the range of heat transfer times for these fish was so long that their peritoneal temperature was probably maintained by engaging in brief dives. The periodicity of dives may lead to a lower fluctuation in the peritoneal temperature, suggesting that the dives are a kind of behavioral thermoregulation.     

Tracking bluefin tuna reproductive migration into the Mediterranean Sea with electronic pop‐up satellite archival tags using two tagging procedures

Thirteen adult bluefin tuna were tracked with electronic pop‐up satellite tags during their reproductive migration towards Mediterranean spawning grounds as they entered the Strait of Gibraltar. Fish were caught in tuna traps and tagged either underwater, with the aid of a modified spear gun, or on the deck of the boat. Fish tagged on board initially showed a shallower behavior than those tagged in the water. The pattern of horizontal movements was also different between both groups. Shortly after tagging, the eight fish tagged in the water entered the Mediterranean Sea. Six of these fish reached the spawning ground located southwest of the Balearic archipelago before heading back for the Atlantic, whereas the other two traveled farther east, reaching its easternmost longitudes between Formentera and Sardinia and the South Tyrrhenian Sea, respectively. In contrast, two out of the five fish tagged on board never entered the Mediterranean Sea, and another one did enter the Mediterranean when the reproductive season was already over. These results suggest an impact of the tagging procedure on the post‐release behavior of bluefin tuna. Excluding the tags that popped‐off east of the Strait of Gibraltar, bluefin tuna stayed in the Mediterranean Sea for 22–28 days. Analysis of the median depth indicated a shallow behavior during both day and nighttime throughout the return phase of the fish from the Mediterranean Sea to the Atlantic Ocean with the exception of the area around the Strait of Gibraltar, where they showed a deeper behavior that coincided with a marked vertical gradient in the currents.     

Bluefin tuna (Thunnus thynnus) distribution in relation to sea surface temperature fronts in the Gulf of Maine (1994–96)

Fishery‐linked aerial surveys for bluefin tuna (Thunnus thynnus) were conducted in the Gulf of Maine (GOM) from July through October, 1994–96. Each year, from 507 to 890 surface schools were detected and their locations examined in relation to oceanographic conditions. Correlations between bluefin tuna presence and environmental variables were explored for sea surface temperature (SST), distance to a SST front, frontal density (relative density of all SST fronts seen in a given 1 km area for 2 weeks prior to each tuna sighting), and bottom depth and slope. Mean SST associated with bluefin schools was 18.1°C (±2.8). Schools were located at a mean distance of 19.7 km (±19.6) from SST fronts, and in water masses with an average frontal density of 28.2 m km^?2 (±35.7). Mean bottom depth of detected schools was 139.0 m (±70.3), and mean bottom slope was 0.7% rise (±0.7). A binomial generalized linear model fit to these variables indicated that bluefin are seen closer to fronts than locations in which no tuna were seen. Using simple and partial Mantel tests, we investigated the spatial correlation between bluefin tuna presence and the environmental variables, controlling for spatial autocorrelation. For each day that schools were sighted, we performed 24 Mantel tests, on a combination of response and predictor variables. The spatial relationship between bluefin tuna and SST fronts was inconsistent. Our analysis identified significant spatial structure in the bluefin school locations that had no significant correlation with any of the measured environmental features, suggesting that other untested features, such as prey density, may be important predictors of bluefin distribution in the GOM.     

Turbulence effect on survival and feeding of pacific bluefin tuna <Emphasis Type="Italic">Thunnus orientalis</Emphasis> larvae,on the basis of a rearing experiment

ABSTRACT:   Physical conditions such as oceanic turbulence related to food availability are considered to be important factors affecting fish larval survival. Rearing experiments were conducted to elucidate the effects of turbulence on the survival and feeding rates during the initial feeding period of Pacific bluefin tuna Thunnus orientalis . Six levels of turbulence intensity were provided by changing flow rates from pipes set on the bottom of rearing tanks. The result showed a dome-shaped relationship between turbulence level and survival rate, in which the feeding rate appeared higher at a logged turbulence energy dissipation rate of −6.32, and decreased at both higher and lower turbulence levels. Compared with the turbulence intensity in the ocean, the optimal turbulence level for Pacific bluefin tuna larvae corresponded to the turbulence caused by sea surface winds with speeds of 4–12.5 m/s. The estimated optimal turbulence intensity for Pacific bluefin tuna larvae is comparable to that for yellowfin tuna Thunnus albacares .     

Retinal ganglion cell topography in juvenile Pacific bluefin tuna Thunnus orientalis (Temminck and Schlegel)

The retinal ganglion cell distribution, which is known to reflect fish feeding behavior, was investigated in juvenile Pacific bluefin tuna Thunnus orientalis. During the course of examination, regularly arrayed cells with a distinctive larger soma, which may be regarded as motion-sensitive cells, were found. The topographical distribution of ordinary-sized ganglion cells, which is usually utilized to estimate fish visual axis and/or visual field characteristics, showed that the highest-density area, termed the area centralis, was localized in the ventral-temporal retina. The retinal topography of ordinary-sized ganglion cells seems to reflect the bluefin tuna’s foraging behavior; while cruising, cells in the area centralis may signal potential prey, such as small schooling pelagic fishes or squids, that are present in the upward-forward direction. Judging from morphological characteristics, the large ganglion cells localized in the small temporal retinal area seem to be equivalent to physiologically categorized off-center Y-cells of cat, which are stimulated by a transient dark spot in a bright visual field. It was inferred that presumed large off-center cells in the temporal retina detect movements of agile prey animals escaping from bluefin tuna as a silhouette against environmental light.     

Effects of open- and closed-system temperature changes on blood O2-binding characteristics of Atlantic bluefin tuna (Thunnus thynnus)

We investigated the effects of open- and closed-system temperature changes on the O₂ affinity of Atlantic bluefin tuna (Thunnus thynnus) blood using in vitro methods essentially identical to those previously employed on tropical tuna species. Bluefin tuna blood has a general O₂ affinity (P ₅₀ = 2.6–3.1 kPa or 19–23 mm Hg at 0.5% CO₂) similar to that of skipjack tuna, yellowfin tuna, and kawakawa blood (P ₅₀ = 2.8–3.1 kPa at 0.5% CO₂) but significantly above that of bigeye tuna blood (P ₅₀ = 1.6–2.0 kPa at 0.5% CO₂). We therefore hypothesize that bluefin tuna are less tolerant of hypoxia than bigeye tuna. Further, we found the P ₅₀ of bluefin tuna blood to be slightly reduced by a 10°C open-system temperature increase (e.g., from 4.83 kPa at 15°C to 3.95 kPa at 25°C) and to be completely unaffected by a 10°C closed-system temperature change. Bluefin tuna blood, therefore, had a significantly reduced Bohr effect when subjected to the inevitable changes in P CO ₂ and plasma pH that accompany closed-system temperature shifts (0.04–0.09 Δlog P₅₀ΔpH⁻¹) compared with the effects of changes in plasma pH accomplished by changing P CO ₂ alone (0.81–0.94 Δlog P₅₀ Δ pH⁻¹). This response is similar to that of skipjack tuna blood, but different from yellowfin or bigeye tuna blood. During closed-system temperature changes at oxygen levels above P ₅₀, however, bluefin tuna blood showed a reversed temperature effect (i.e., P O ₂ decreased in response to an increase in temperature). Unlike in other tuna species, temperature effects on O₂ affinity of bluefin tuna whole blood were similar to those previously reported for hemoglobin solutions, suggesting that red cell-mediated ligand changes are not involved.     

Variability of preferred environmental conditions for Atlantic bluefin tuna (Thunnus thynnus) larvae in the Gulf of Mexico during 1993–2011

The Gulf of Mexico (GOM) is the primary spawning ground for western Atlantic bluefin tuna (Thunnus thynnus). In this work, information reported by previous studies about the preferred environmental conditions for the occurrence of bluefin tuna larvae in the GOM is integrated into a dimensionless index, the BFT_Index. This index is used to evaluate the spatial and temporal variability of areas with favorable environmental conditions for larvae within the GOM during 1993–2011. The main findings of this work are that: (i) the proposed index successfully captures the spatial and temporal variability in the in situ occurrence of bluefin tuna larvae; (ii) areas with favorable environmental conditions for larvae in the GOM exhibit year‐to‐year spatial and temporal variability linked with mesoscale ocean features and sea surface temperature; and (iii) comparison of the BFT_Index‐derived variability with recruitment of age‐0 fish estimated from recent stock assessment indicates that changes in environmental conditions may drive a relevant component (~58%) of the recruitment variability. The comparison with the recruitment dataset further revealed the existence of key regions linked with recruitment in the central/northern GOM, and that the Loop Current may function as a trap for larvae, possibly leading to low survival rates. Above (below) average conditions for occurrence of larvae in the GOM during spring were observed in 2000, 2001, 2002, 2006–2008, and 2011 (1994, 1996, 1998, 1999, 2003 and 2010). Results reported here have potential applications to assessment of bluefin tuna.     

美济礁附近海域3种金枪鱼肌肉成分检测与营养评价

为分析中国南海美济礁附近海域大目金枪鱼(Thunnus obesus)、蓝鳍金枪鱼(T.thynnus)、黄鳍金枪鱼(T.albacares)的营养成分与营养价值,该研究采用国标法检测了3种金枪鱼的水分、蛋白质、灰分、脂肪酸和氨基酸成分及其结构,并进行营养评价。结果显示,3种金枪鱼中黄鳍金枪鱼蛋白质含量最高,蓝鳍金枪鱼脂肪和灰分含量最高;以氨基酸评分(AAS)和必需氨基酸指数(EAAI)为评分标准,3种金枪鱼氨基酸评分均大于等于1;3种金枪鱼各种必需氨基酸构成均优于联合国粮食及农业组织(FAO)/世界卫生组织(WHO)模式;蓝鳍金枪鱼的油酸含量远高于其他2种金枪鱼,且差异显著(P<0.05);蓝鳍金枪鱼的单不饱和脂肪酸(MUFA)和多不饱和脂肪酸(PUFA)种类最多且含量最高。结果表明,3种金枪鱼均具有较丰富的营养成分、较完美的营养比例和较高的营养价值;3种金枪鱼相比,黄鳍金枪鱼可提供更多的蛋白质,蓝鳍金枪鱼提供更多的不饱和脂肪酸。美济礁附近海域3种金枪鱼具有很高的开发价值,成分检测对研究设计这3种金枪鱼专用驯化养殖饲料有重要的参考价值。     

Movement and behaviour of large southern bluefin tuna (Thunnus maccoyii) in the Australian region determined using pop‐up satellite archival tags

Pop‐up satellite archival tags (PSATs) were deployed on 52 large (156–200 cm length to caudal fork) southern bluefin tuna (Thunnus maccoyii) in the western Tasman Sea during the austral winters of 2001–2005. Southern bluefin tuna (SBT) were resident in the Tasman Sea for up to 6 months with movements away from the tagging area occurring at highly variable rates. The data indicated a general tendency for SBT to move south from the tagging area in the Western Tasman Sea. Four individuals migrated west along the southern continental margin of Australia and into the Indian Ocean. Three individuals moved east into the central Tasman Sea, with one individual reaching New Zealand. We also describe the first observed migration of an SBT from the Tasman Sea to the Indian Ocean spawning grounds south of Indonesia. Individuals spent most of their time relatively close to the Australian coast, with an estimated 84% of time spent in the Australian Fishing Zone. SBT favored temperatures between 19 and 21°C, adjusting their depth to the vertical temperature distribution. Distinct diurnal diving patterns were observed and adjustment of depth to maintain constant ambient light levels over a 24‐h period. The findings of this study are a significant advance toward greater understanding of the spatial dynamics of large SBT and understanding the connectivity between distant regions of their distribution.     

Possibility for decreasing of mercury content in bluefin tuna <Emphasis Type="Italic">Thunnus orientalis</Emphasis> by fish culture

ABSTRACT:   The bluefin tuna tested were reared for 22 months from eggs before the beginning of the experiment, and sampling was performed every 3 months over the following year. The experimental results showed that the mercury concentration in the muscle ranged from 0.32 to 0.85 μg/g, which is lower than that found in wild bluefin tuna of a similar size. Increase in the mercury concentration corresponding to the increase in body weight was not shown, and it was quite different with wild bluefin tuna. Furthermore, no significant relationship was found between the lipid concentration and the mercury concentration in muscle. Among the internal organs of cultured bluefin tuna, the heart (0.32–0.66 μg/g), liver (0.43–0.99 μg/g) and spleen (0.59–1.0 μg/g) contained higher concentrations of mercury. It was estimated that the full-cycle cultured bluefin tuna had been fed small fish containing lower concentrations of mercury, and that the mercury concentration of tuna would be almost equal throughout the year because the effect of mercury accumulation would be weakened by body growth. Therefore, it was concluded that selecting diet fish species might decrease mercury contents in cultured bluefin tuna.     

Recruitment abundance index of Pacific bluefin tuna using fisheries data on juveniles

ABSTRACT:   The recruitment abundance index of Pacific bluefin tuna Thunnus orientalis was estimated from 1980 to 2003 fishing year by using the troll fishery data in Nagasaki Prefecture, western Japan. It has been shown that the troll fishery in Nagasaki Prefecture operates with good time–area coverage of the species habitat, and that the fishing power slightly changed during the period analyzed, based on fisheries statistics, published information, and interviews with the fishers. Average catch per unit effort (CPUEs) were standardized by a generalized linear model (GLM) considering the effects of fishing year, season and landing area. Standardized CPUE of age-0 bluefin tuna showed larger fluctuations year by year than the nominal CPUE combined for all ages. High CPUEs in fishing years of 1981, 1994, 1996 and 1999 were observed. Data from these years agreed with the higher recruitments estimated by virtual population analysis (VPA) or higher catch of age-0 fish reported for the Pacific side. The age-specific standardized CPUE of age-0 bluefin tuna in this study was judged to be a useful indicator of recruitment.     

Diving behavior of immature, feeding Pacific bluefin tuna (Thunnus thynnus orientalis) in relation to season and area: the East China Sea and the Kuroshio–Oyashio transition region

Twenty‐four archival tags were recovered from Pacific bluefin tuna previously released off Tsushima Island in the East China Sea. By analysis of the time‐series data of the pressure and the ambient and internal temperature from the 24 tags, we examined the relationship between the tuna's pattern of diving and the thermocline depth. In the East China Sea, diving and feeding events occurred throughout almost the entire day in both winter and summer, suggesting that the purpose of diving is for feeding. In summer, the feeding frequency was greater than that in winter, which corresponds to the fact that growth is more rapid in summer than in winter. During summer in the Kuroshio–Oyashio transition region, on the other hand, feeding events were much more frequent than those in the East China Sea, in spite of a lower diving frequency. The mean horizontal distance traveled was also significantly higher and it seems that in this area they may move horizontally to feed on prey accumulated at the surface. We conclude that, in addition to the ambient temperature structure, the vertical and horizontal distribution of prey species plays an important role in the feeding behavior of Pacific bluefin tuna. One bluefin tuna migrated to the Oyashio frontal area, where both the horizontal and the vertical thermal gradients are much steeper. The fish spent most of the time on the warmer side of the front and often traveled horizontally to the colder side during the day, perhaps to feed. This implies that there is a thermal barrier effect, in this case from the Oyashio front, on their behavior. The frequency of feeding events was low, although all the monitored fish dived every dawn and dusk, irrespective of the seasons or location. It is possible that these twice‐daily diving patterns occurred in response to the change in ambient light at sunrise and sunset.     

Characterization and ontogenetic development of digestive enzymes in Pacific bluefin tuna Thunnus orientalis larvae

The major digestive enzymes in Pacific bluefin tuna Thunnus orientalis larvae were characterized, and the physiological characteristics of the enzymes during early ontogeny were clarified using biochemical and molecular approaches. The maximum activity of trypsin (Try), chymotrypsin (Ct) and amylase (Amy) was observed at pH 6–11, 8–11 and 6–9, respectively. Maximum activity of Try, Ct and Amy occurred at 50 °C, that of lipase (Lip) was at 60 °C and that of pepsin (Pep) was at 40–50 °C. These pH and thermal profiles were similar to those for other fish species but differed from those previously reported for adult bluefin tuna. Enzyme activity for all enzymes assayed was found to decrease at high temperatures (Try, Ct, Amy and Pep: 50 °C; Lip: 40 °C), which is similar to findings for other fish species with one marked exception—increased Try activity was observed at 40 °C. Lip activity appeared to be dependent on bile salts under our assay conditions, resulting in a significant increase in activity in the presence of bile salts. Ontogenetic changes in pancreatic digestive enzymes showed similar gene expression patterns to those of other fish species, whereas marked temporal increases in enzyme activities were observed at 10–12 days post hatching (dph), coinciding with previously reported timing of the development of the pyloric caeca in bluefin tuna larvae. However, complete development of digestive function was indicated by the high pep gene expression from 19 dph, which contradicts the profile of Pep activity and previously reported development timing of the gastric gland. These findings contribute to the general knowledge of bluefin tuna larval digestive system development.