Ontogenetic phase shifts in metabolism in a flounder Paralichthys olivaceus.

Size-scaling metabolism is widely considered to be of significant importance in biology and ecology. Thus, allometric relationships between metabolic rate (VO2) and body mass (M), VO2 = a(i)Mb, have long been a topic of interest and speculation. It has been proposed that intraspecifically metabolic rate scales isometrically or near isometrically with body mass during the early life history in fishes, invertebrates, birds and mammals. We developed a new perspective on intraspecific size-scaling metabolism through determination of metabolic rate in the Japanese flounder, Paralichthys olivaceus, during their early life stages spanning approximately four orders of magnitude in body mass. With the increase of body mass, the Japanese flounder had four distinct negative allometric phases in which three stepwise increases in scaling constants (ai, i = 1-4), i.e. ontogenetic phase shifts in metabolism, occurred with growth during its early life stages at around 0.002, 0.01 and 0.2 g, maintaining each scaling exponent constant in each phase (b = 0.831). These shifts in metabolism during the early life stages are similar to the tiger puffer, Takifugu rubripes. Our results indicate that ontogenetic phase shifts in metabolism are key to understanding intraspecific size-scaling metabolism in fishes.

Metabolic rate of living things is related to many biological traits such as body size, growth, phylogeny of species, rates of reproduction and genome evolution1234567. Body size of animals primarily constrains the metabolic rate through the size-scaling effect891011121314 and these two traits have co-evolved15. It has been established that the relationship between metabolic rate () and body mass (M) in animals is expressed by the allometric formula, where “a” is a scaling constant and “b” is the scaling exponent. However there are many cases where a simple allometric equation can not explain non-linear or curvilinear log-log metabolic scaling relationships16171819. Furthermore intraspecific size-scaling metabolism of animals may be a complex process in which several phases can be distinguished16202122.

Fishes are good candidates for intraspecific studies because, unlike other vertebrates such as mammals, birds or reptiles, most fishes develop from a small body size (mg level or less)2324 and attain a mass range approximately comparable to that between a mouse and an elephant25. In addition, it is not necessary to consider the phylogeny of species262728. Although reported scaling exponents in various species of fish are quite variable (approximately 0.4 to 1.3)2629, a scaling exponent has often been suggested as close to unity in the early stages of development of marine fish larvae30. On the other hand, in juveniles and adults, the relationship between metabolic rate and body mass is negatively allometric, and the mass exponents for simple regression of metabolic rate on body mass are usually less than 1. Consequently, a biphasic relationship where the scaling exponent changes from isometry or near isometry during the larval phase to negative allometry during the juvenile phase has been proposed in many fish species2531323334. Although these studies concluded that both the metabolic scaling exponent and scaling constant varies during ontogeny, the generality of the isometric scaling remains an open question requiring additional empirical evidence.

Recently we have demonstrated that in the larval and juvenile tiger puffer, Takifugu rubripes, oxygen consumption (), as a proxy for metabolism, scales with M as , and that three stepwise increases in scaling constants ai (i = 1–4), i.e. ontogenetic phase shifts in metabolism, occur, maintaining each scaling exponent constant in each phase22. The overall scaling exponent (0.948 ± 0.002) (estimate ± s.e.m.) was steeper than the scaling exponents for each developmental stage (0.795 ± 0.019)22. This result suggests that isometric or near-isometric metabolic scaling during early life stages may be produced by a combined effect of ontogenetic phase shifts in metabolism. Thus, ontogenetic shifts in metabolism are likely to occur in other teleost fish species. In the present paper, we examine if the ontogenetic phase shifts in metabolism occur with growth in the Japanese flounder Paralichthys olivaceus (Temminck & Schlegel, 1846).

Results

Rates of oxygen consumption ( in µl O2 fish−1 min−1) in relation to body mass (M in g) are plotted in Figure 1. Mass-specific rates of oxygen consumption (/M in µl O2 g−1 min−1) are also presented. increased daily from just after hatching to 7 days after hatching (DAH), with virtually no increase of body mass. After 7 DAH, body mass increased. There was no substantial difference between the results of the two different methods of respirometry (Figure 1). Therefore, the values obtained were used without any distinction for the relationship between oxygen consumption and body mass.

Figure 1

Ontogenetic changes in the rate of respiration (, diamonds) and the mass-specific rate of respiration (, circles) with increase of body mass (M) in Japanese flounder.

Symbols for indicate the two methods of respirometry (open: closed method and solid: semi-closed method). Symbols for signify the year in which data were collected (open: 2005 and solid: 2006). The vertical broken lines at around 0.0003 g represent and which increased daily from just after hatching to 7 days after hatching, with virtually no increase in body mass. Ranges covered by the four solid lines each for and indicate intragroup phases of negative allometry. The broken lines both on and represent the intergroup lines. Small symbols represent values during the transitional phases. Regression analysis of each line for is presented in Table 1, and ANCOVA in Table 2.

Two models were applied to compare the four negative allometric relationships: for each incidence of negative allometry, and for the overall line constituting these allometries3536. ‘ai’ represents an intragroup scaling constant of the i'th group, and ‘α’ the intergroup one of the groups. Equations (2) and (3) were rewritten as and where y is , x is log10M and ε and E represent the random intra- and intergroup variation in metabolism. To estimate and log10 α, we used the ordinary least-squares regression to minimise the sum of squares of µ, which is the vertical distance of the group mean (,) from the overall line (intergroup line)35. Because log10 ai is equal to (), Equation (4) was rewritten as follows2236: The scaling exponents of individual lines were significantly smaller than unity for all regression lines (P<0.05; two tailed t-test). Regression analysis of each group, except the transitional phases is given in rows 1 to 4 of Table 1. The slopes of the individual lines in each of the four groups were not significantly different (F3, 94 = 1.95, P = 0.127; one-way ANCOVA). The intragroup scaling exponent was estimated to be 0.831 ± 0.0026 (estimate ± S.E.M.), and logarithm of the scaling constants in each group at = 0.831 was estimated as follows: log10a1 = 0.373 ± 0.083, log10a2 = 0.425 ± 0.062, log10a3 = 0.539 ± 0.034, and log10a4 = 0.584 ± 0.018. The scaling constants were calculated as follows: a1 = 2.36, a2 = 2.66, a3 = 3.46, and a4 = 3.84. Regression analysis of the intergroup is given in rows 5–7 of Table 1. The intergroup scaling exponents was estimated to be 0.915 ± 0.0011.

Table 1

Intragroup (rows 1 to 4), intergroup (rows 5 to 7) regression analysis of the relationship between (: oxygen consumption in µl O2 fish−1 min−1) and log10 M (M: body mass in g) in the Japanese flounder

		Range of body mass	Scaling	Scaling exponent
Group	N	(g)	constant	(Mean ± S.E.M.)	P	R2
1*	24	0.00041–0.0015	0.80	0.680 ± 0.065	6.26 × 10−5	0.832
2*	26	0.0025–0.0099	2.57	0.825 ± 0.066	1.42 × 10−2	0.866
3*	36	0.016–0.14	3.92	0.875 ± 0.035	1.12 × 10−3	0.948
4*	16	0.24–0.90	3.67	0.782 ± 0.066	5.04 × 10−3	0.910
1–4†	102	0.00041–0.90	α = 4.28		6.58 × 10−91	0.9999
1–4§	102	0.00041–0.90	4.20	0.911 ± 0.0061	1.55 × 10−26	0.996
Total§,¶	124	0.00026–0.90	4.32	0.922 ± 0.0060	1.73 × 10−24	0.995

*Parameters were estimated to minimise the sum of squares of ε in each group; †parameters were estimated to minimise the sum of squares of µ; §parameters were estimated to minimise the sum of squares of E; ¶including the transitional phases. N is the number of determinations; P is the difference of scaling exponent from unity, examined using Students t-test, two tailed; R2 is squared correlation coefficient between and log10 M.

A one-way ANCOVA was carried out to clarify the validity of Equation (6), and the results are provided in Table 2. The µ's were not different from zero (F2, 97 = 1.71, P = 0.186; one-way ANCOVA), meaning that the intragroup means might lie on the intergroup (overall) line. The intragroup scaling exponent was significantly different from the intergroup exponent (F1, 97 = 9.67, P = 0.00246; one-way ANCOVA). This implies that the intragroup scaling constant ai increased significantly from a1 = 2.36 to a4 = 3.84 with increase in body mass (Figure 1). Thus, the metabolic rate is expressed by , in which ai increased three times during the transitional phases.

Table 2

ANCOVA table for respirometry in the Japanese flounder based on the model (equation (6))

		Degrees of		Mean-square
Term	Sum of squares	freedom	Mean square	ratio	P
log α	110.12189	1	110.121898	34124	2.26 × 10−125
µi	0.011030379	2	0.005515	1.71	0.186
	0.031216	1	0.031216	9.67	0.00246
	79.661712	1	79.661712	24685	1.42 × 10−118
εij	0.313031	97	0.003227
Total (about mean)	80.016990	101
Total (about zero)	190.138887	102

Discussion

The size-scaling metabolism for the Japanese flounder during its early life stages (that is, spanning approximately four orders of magnitude in body mass) showed ontogenetic phase shifts, because the scaling constant ai increased three times in the course of development, maintaining each scaling exponent constant in each phase () (Figure 1). These shifts in metabolic scaling are consistent with the findings of Yagi et al.22, who reported that a puffer fish had four distinct phases in which three stepwise increases in scaling constants, i.e. ontogenetic phase shifts in metabolism, occur with growth during its early life stages, maintaining each scaling exponent constant in each phase (). However, when comparing the scaling constants for the Japanese flounder to the tiger puffer ones, the range for the Japanese flounder (2.36 to 3.84) was narrower than what was observed in the previous study22 for the tiger puffer (2.75 to 7.25). This difference in the scaling constants between two species could be attributable to differences in the developmental trajectory. The Japanese flounder larvae initially resemble typical fish swimming in a vertical attitude as well as the tiger puffer. At ~0.01 g, they change their mode of life to swim in a horizontal attitude after metamorphosis, becoming bottom-dwelling fish. Juvenile and later stage Japanese flounder showed almost no activity in respirometers, lying on their floor. Therefore, measures of resting routine metabolic rate closely approximate resting metabolic rate after metamorphosis. The tiger puffer juveniles move to estuary and neighbouring mudflats at ~0.1 g, and grow there for a while in Hakata Bay and Ariake Sound, Japan22. After that, they show active benthic swimming to migrate to the East China Sea and the Yellow Sea. Although further research is required, our results suggest that differing patterns of development and activity status during ontogeny can cause variation in metabolic scaling between species.

Previously, the intraspecific metabolic rate of larval and juvenile fishes was thought to scale isometrically (b = 1.0) or near isometrically to body mass30323337383940. Isometric or near isometric scaling during early life stages have also been reported in invertebrates, birds and mammals (reviewed in Glazier16). These reports are not surprising, because isometric or near isometric scaling of metabolic rate may be made by the increases of ai during ontogeny. For example, when the scaling exponent was estimated from 0.00026 g (0 DAH) to 0.0055 g (30 DAH), it would be 1 (P = 0.052; two tailed t-test) for the Japanese flounder. This suggests that the ontogenetic phase shifts in metabolism during early life stages may not be specific only to the Japanese flounder and the tiger puffer. In fact, further analyses demonstrate that similar sequential shifts in scaling coefficient (a) over development have been observed in the tobacco hornworm41. However, Blossman-Myer & Burggren42 have found a different pattern to that found in our study: an overall scaling exponent of 0.82 versus isometric scaling within development stages (instars) of the silkworm. Although a segmented regression was applied in the present study, other functions may be more biologically relevant (e.g. techniques that allow gradual, rather than an abrupt change in the scaling exponent). In any case, we speculate that such hidden ontogenetic phase shifts in metabolism may play a key role in the linear intraspecific size-scaling metabolism that has been observed in animals.

Fishes develop from larvae, whose morphological features are appreciably different from adults, to juveniles which have adult-like features43. Previously, the developmental stages during early life for fishes have been distinguished by morphological features, such as yolk-sac, and pre-flexion, flexion and post-flexion larvae44. In this study, the ontogenetic phase shifts in metabolism at ~0.002, ~0.01 and ~0.2 g (~6, ~10 and ~25 mm in standard length; M. Yagi and S. Oikawa, unpublished data) during the transitional phases were accompanied by morphological and behavioural changes. At ~0.002 g larvae, rapid increment of swimming velocity was observed45. At ~0.01 g, the transformation from larval to juvenile stage occurs46. Fukuhara46 also reported that a relatively sharp increment of maximum swimming speed was observed for ~0.01 g juveniles. At ~0.2 g, juveniles are observed to migrate from their nursery grounds to offshore habitats. This migration is also associated with the shift of feeding habits from mysids to fish47. Thus, our results suggest that developmental stages of ontogenetic phase shifts in metabolism may be more appropriate to explain behavioural and eco-physiological traits during early life stages of fishes. In fact, there is evidence that, in the tiger puffer, increases in ai (metabolic rate) facilitate anti-predator adaptation, in which the majority of intracohort predation occurs on smaller fish that had a lower metabolic rate22. Individuals that enter the next metabolic phase with higher metabolic rate are considered to become stronger, with higher motility accompanying morphological and behavioural changes than extrapolated from the previous metabolic phase with lower metabolic rate22. Future studies should investigate the cellular, tissular and physiological mechanisms associated with ontogenetic phase shifts in metabolism, and their possible relevance to prey-predator interactions in nature.

Methods

Fish used

Fish (P. olivaceus) were hatched from artificially fertilized eggs obtained from wild parents in the northern part of Kyushu, Japan that were captured by fishermen. Experiments were conducted at the Fisheries Research Laboratory, Kyushu University, Fukuoka, Japan in 2005 and 2006. All larvae and juveniles were maintained in aerated 500 l tanks supplied with a constant flow of seawater and fed live rotifers, Brachionus rotundiformis twice daily between 2 and 22 DAH (~5000 l−1), live brine shrimp, Artemia sp. larvae twice daily two times between 20 and 32 DAH (~2000 l−1), and artificial diets three times daily thereafter. Live diets were fortified with essential fatty acids, EPA and DHA using Super Capsule Powder (Chlorella Industry, Tokyo, Japan), before feeding. The water temperature in the rearing tanks was held at 18°C (this temperature was chosen to approximate conditions experienced by the early life stages of the flounder in the wild and was also used when measuring the oxygen consumption). Fish used in the respirometry study were not fed for 3 to 24 h before experiments, depending on the fish body size, and so the individuals used for respirometry did not have food in their guts. All experimental procedures were approved by a KU committee and conducted in accordance with the Guideline for the Care and Use of Laboratory Animals of Japan, with registrations of M.Y. and S.O. for Animal Experiments in the Faculty of Agriculture and in the Graduate Course of KU.

Respirometry

Oxygen consumption () was measured in larval and juvenile fish ranging in size between 0.00026 g (wet body mass, 0 DAH) and 0.90 g (80 DAH). Resting routine rates of oxygen consumption, i.e. intermediate between the resting and routine activity states48, were determined in fasted larvae and juveniles at 18°C using one of two methods depending on the developmental stage of the fish, a closed method and a semi-closed method, based on previous studies2232. The larvae and juveniles were placed in the respiration chamber and allowed to settle for an appropriate acclimatisation period before measurement to remove possible stress caused by handling (1 h for the closed method and 1 to 4 h for the semi-closed method). The closed method was only used for larvae just after hatching to 8 DAH, because of their poor swimming ability. The closed method was based on depletion of oxygen in water in a sealed glass oxygen bottle (20 ml). A blank chamber without fish was used to eliminate background respiration. The semi-closed method was essentially a closed method in which the glass chamber was slowly flushed with air-saturated water before determination and closed during determination. Several respiration chambers (50 to 650 ml) were used depending on the fish body size (50 ml for 8–25 DAH, 100 ml for 18–44 DAH, 250 ml for 46–63 DAH, 350 ml for 46–80 DAH, 650 ml for 63–80 DAH). In the semi-closed method, the bottle which received water flowing out of the respiration chamber was used as the blank chamber. This bottle was sealed at the beginning of determination of oxygen consumption by fish, placed in the water bath for the respiration chamber during determination, and the oxygen concentration in the bottle was determined at the end of the measurement. By using this value as the initial oxygen concentration in the respiration chamber, background respiration was cancelled.

Wet body mass of experimental fish was directly determined immediately after respirometry, except for specimens used in the closed method. Wet body mass during the larval stage was determined as detailed by Yagi et al.22. In the closed method, body mass was indirectly estimated from other co-cultured individuals of similar body size, because they were fixed when the dissolved oxygen concentration of water in the respiration chamber was determined by the Winkler's titration method.

Data analysis

Rates of oxygen consumption () for the Japanese flounder were plotted on log10 plots against wet body mass (M), and four negative allometric relationships were fitted between and M, based on the developmental stages, which is similar to the tiger puffer Takifugu rubripes22, interposing three transitional phases of at approximately 0.002, 0.01 and 0.2 g of wet body mass. For the tiger puffer, at approximately ~0.002 g, the primordia of the fin rays, at ~0.01 g post-flexion larvae metamorphosed into juveniles, and at ~0.1 g squamation of the body proceeded22. For the Japanese flounder, at approximately ~0.002 g larvae, the ventral fin buds are present and the base of the dorsal and anal fins (primordia) appear, at ~0.01 g, development of fin rays is completed, and at ~0.2 g, squamation on the caudal peduncle is completed, consisting of eight scale rows46. Thus, in this study, three transitional phases of at approximately 0.002 (0.0015–0.0020 g), 0.01 (0.010–0.015 g) and 0.2 (0.15–0.22 g) g of wet body mass were interposed. Size-scaling analyses were performed by applying the statistical model described in Yagi et al.22.

Author Contributions

M.Y. and S.O. conceived of and designed the research. M.Y. cared for research animals. M.Y. and S.O. collected, processed and interpreted data. M.Y. drafted the manuscript. All authors read and commented on the manuscript.