# The relationships between growth rate and mitochondrial metabolism varies over time

Jean-Baptiste Quéméneur, Morgane Danion, Joëlle Cabon, Sophie Collet, José-Luis Zambonino-Infante, Karine Salin

► To cite this version:

Jean-Baptiste Quéméneur, Morgane Danion, Joëlle Cabon, Sophie Collet, José-Luis Zambonino-Infante, et al.. The relationships between growth rate and mitochondrial metabolism varies over time. Scientific Reports, 2022, 12 (1), pp.16066. 10.1038/s41598-022-20428-9 . anses-03980511

HAL Id: anses-03980511

https://anses.hal.science/anses-03980511

Submitted on 9 Feb 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# The relationships between growth rate and mitochondrial metabolism varies over time

Jean-Baptiste Quéméneur<fnref n="1" />, Morgane Danion², Joëlle Cabon², Sophie Collet<fnref n="1" />,
José-Luis Zambonino-Infante<fnref n="1" /> & Karine Salin<fnref n="1" />✉

Mitochondrial metabolism varies significantly between individuals of the same species and can influence animal performance, such as growth. However, growth rate is usually determined before the mitochondrial assay. The hypothesis that natural variation in mitochondrial metabolic traits is linked to differences in both previous and upcoming growth remains untested. Using biopsies to collect tissue in a non-lethal manner, we tested this hypothesis in a fish model (Dicentrarchus labrax) by monitoring individual growth rate, measuring mitochondrial metabolic traits in the red muscle, and monitoring the growth of the same individuals after the mitochondrial assay. Individual variation in growth rate was consistent before and after the mitochondrial assay; however, the mitochondrial traits that explained growth variation differed between the growth rates determined before and after the mitochondrial assay. While past growth was correlated with the activity of the cytochrome c oxidase, a measure of mitochondrial density, future growth was linked to mitochondrial proton leak respiration. This is the first report of temporal shift in the relationship between growth rate and mitochondrial metabolic traits, suggesting an among-individual variation in temporal changes in mitochondrial traits. Our results emphasize the need to evaluate whether mitochondrial metabolic traits of individuals can change over time.

Growth performance is closely related to an individual's fitness by influencing body size and, in turn, age at maturity, fecundity rate, foraging ability, and predator escape<sup>1,2</sup>. The physiological processes involved in somatic growth such as cell division and protein synthesis require ATP, which is mainly produced by the mitochondria<sup>3</sup>. Interestingly, individuals that are most efficient at converting nutrients into ATP will typically have faster growth than those that do so less efficiently<sup>4,5</sup>. This can occur either directly via higher mitochondrial efficiency of ATP production (experimentally tracked as the amount of ATP generated per molecule of oxygen consumed, ATP/O ratio) or indirectly through higher rates of oxygen consumption (OXPHOS) for ATP production or more mitochondria<sup>6,7</sup>. Conversely, variation in the availability in ATP between individuals might also occur because individuals differ in the proportion of the oxygen that is consumed by mitochondria to offset the proton leakage (LEAK respiration) across the inner mitochondrial membrane<sup>8</sup>. In each case, the somatic growth can be upregulated from more cell division and protein synthesis and hence biomass production with availability in ATP.

Our current understanding of the relationship between mitochondrial properties and growth performance is largely based on studies of growth trials that occur over a period of time anterior to mitochondrial measurement<sup>4,5,9</sup>, hereafter called past growth. The main reason for this caveat is that studying mitochondrial metabolic traits usually involves terminal sampling. If mitochondrial properties are only determined toward the end of the growth trial, it is the variation in growth performance that might actually promote differences in mitochondrial metabolic traits among individuals. Hence, ascertaining the causality in a relationship between mitochondrial metabolic traits and growth performance in the absence of a clear temporal relationship can be a source of debate. Previous work in tadpoles has demonstrated that experimental variation in mitochondrial metabolic traits lead to variation in future growth<sup>10</sup>. However, it is also likely that mitochondrial variation induced by an artificial substance, the 2,4-Dinitrophenol, is not representative of natural variation in mitochondrial metabolic traits. Natural variation in mitochondrial metabolic traits might be related to individual variation in future growth, and, in turn, predict the rate of growth subsequent to the mitochondrial assay, but that has not been documented previously.

<footnote n="1">¹Ifremer, Laboratory of Environmental Marine Sciences, University Brest, CNRS, IRD, 29280 Plouzané, France. ²Anses, Ploufragan-Plouzané Niort Laboratory, VIMEP Unit, Technopôle Brest-Iroise, 29280 Plouzané, France. ✉email: Karine.salin@ifremer.fr</footnote>
Investigations of the relationship between mitochondrial metabolic traits and growth performance have generally focused on a growth trial that represents a relatively short period of time. For example, individual brown trout and broilers that had lower rates of LEAK respiration exhibited better growth performance, where growth was measured over 11 to 20 days in trout¹¹ and 7 days in broilers⁵ before the measurement of mitochondrial function. Similarly, study in juvenile brown trout demonstrates a positive association between mitochondrial ATP/O ratio and growth, but growth was determined over a trial of 14 days for a trout that can live for years⁴. Growth rate is generally repeatable over time when measured under constant conditions¹². However, metabolic traits of mitochondria are flexible and can change dramatically in as little as few hours¹³. Experimental studies now show that mitochondrial metabolic traits can change in less than a month in different taxa, including insects¹⁴, bivalves¹⁵, fishes¹⁶, birds¹⁷ and mammals¹⁸. However, individuals can differ in the degree to which their mitochondrial traits change as a function of time and therefore the link between mitochondrial traits and growth performance may decline over time.

In this study, we measured mitochondrial metabolic traits using a non-terminal technique for tissue sampling in fish, which can then be used to measure performance following the sampling time. Using European sea bass (Dicentrarchus labrax), we studied whether individual variations in past and future growth rates were related to individual variations in mitochondrial metabolic traits. We tested the predictions that (i) individual variation in growth rate, both before and after the mitochondrial assay, is related to individual variation in similar mitochondrial traits, and (ii) the relationship between growth rate and mitochondrial traits declines over time. To test these hypotheses, we examined the growth rates up to 20 weeks before and 12 weeks after mitochondrial measurement, respectively. We analyzed mitochondrial metabolic traits (rates of OXPHOS and LEAK respiration, ATP production, and cytochrome c oxidase (COX) activity, a measure of mitochondrial density) in the red muscle. We chose this tissue because it can be sampled using biopsy punches to collect a muscle tissue plug in a non-lethal manner¹⁹,²⁰ and has, in contrast to white muscle, a particularly high mass-specific mitochondrial metabolism²¹,²².

## Methods

The experiments were approved by the French Ethics Committee in charge of Animal Experimentation (no.2019072411491441) and were in accordance with institutional and ARRIVE guidelines.

### Animal collection and husbandry.
In May 2019, juvenile European sea bass, Dicentrarchus labrax (Linnaeus 1758) (6 months old, mass 5 g), were transferred from a fish farm (Turbot Ichtus, Trédarzec, France) to the Ifremer rearing facility (Plouzané, France). Fish were kept in a common tank for 5 months, maintained under a 12 L: 12 D photoperiod, and fed at satiety three times a week using commercial pellets (Neo Start, Le Gouëssant, Lamballe, France).

In October 2019, fish (n = 40) were anaesthetized (Tricaïne; 125 mg L⁻¹), weighed (41.5 ± 1.8 g, MCE11201S-2S00-0, Sartorius, Göttingen, Germany), and implanted subcutaneously with an identification tag (RFID; Biologid, Bernay, France). The fish were then randomly allocated to ten replicate 400 L tanks supplied with open-flow, fully aerated seawater (oxygen saturation > 95%, salinity 32 ppt), thermo-regulated during winter to avoid falling below 13 °C, and fed at satiety three times a week. Temperature was recorded weekly. To account for the potential effect of temperature variation over the duration of the trial (15.5 ± 0.5 °C, range: 13.1–17.9 °C) on growth, a correlations analysis was performed between temperature and specific growth rate (SGR). No statistical relationship was found between SGR and temperature (Spearman R² = 0.060, P = 0.596). Additional fish (n = 40) were present in the tanks (final density: n = 8 per tank) for the need of another project.

### Growth measurements.
Body mass (BM) was measured about every four weeks from October 2019 to June 2020. The fish were fasted for 48 h and anesthetized before each BM measurement (± 0.1 g). The specific growth rate (% day⁻¹) was estimated as follows:

$$ \text{Specific Growth Rate} = \frac{\ln (\text{final BM}) - \ln (\text{initial BM})}{\text{days elapsed}} \times 100 $$

In March 2020, a red muscle biopsy sample was collected from fish to measure the mitochondrial metabolic traits. Past growth was defined as specific growth rates before the analysis of mitochondrial metabolic traits (past specific growth rate, SGR<sub>past</sub>). SGR<sub>past</sub> were calculated using the BM at the muscle biopsy as the final BM and the BM at 7, 11, 16, and 20 weeks before the muscle biopsy as the initial BM (Fig. 1). Future growth was defined as specific growth rates after analysis of mitochondrial metabolic traits (future specific growth rates, SGR<sub>future</sub>). SGR<sub>future</sub> were calculated using the BM at 4, 8, and 12 weeks after the muscle biopsy as the final BM and the BM at the muscle biopsy as the initial BM. In European sea bass, most of the somatic growth occur within the first 3 to 5 years of life, so several months of growth measurement at the juvenile stage might be representative of the overall growth of the animal.

### Muscle biopsy procedure.
Muscle biopsy was performed as a non-lethal means of sampling tissue for the mitochondrial assay while allowing us to determine future growth rate. Fish were anaesthetized with tricaïne (as above), weighed (76.7 ± 3.6 g), and biopsied. A skin incision (< 10 mm in length) was made with a scalpel below the lateral line and between the dorsal and caudal fins. Then, a core of the red muscle was collected using a biopsy punch (2 mm LCH-PUK-20, Kai Medical, Solingen, Germany). The core of red muscle was immediately cleaned of white muscle with a scalpel. The red muscle sample was then weighed (5.6 ± 0.2 mg; AC210P-0F1, Sartorius, Göttingen, Germany), and transferred to ice-cold respiration buffer (20 mmol L⁻¹ Taurine, 10 mmol L⁻¹
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Experimental design. Juvenile European sea bass (n=40) were weighted about every four weeks over a 32-week period. At week 20, a biopsy of red muscle was used for mitochondrial assay. Specific growth rates (SGR) were calculated relative to the time of the biopsy. Past growth rate corresponds to SGR calculated before the biopsy, and future growth rate corresponds to SGR calculated after the biopsy.</figcaption></figure>

KH₂PO₄, 20 mmol L⁻¹ HEPES, 110 mmol L⁻¹ D-sucrose, 60 mmol L⁻¹ K-lactobionate, 1 g L⁻¹ BSA fatty acid free, and pH 7.0 at 13 °C). The incision was disinfected (Vétédine® Solution, Vetoquinol, Magny-Vernois, France) and filled with powdered bandage (ORAHESIVE, ConvaTec®, Deeside, UK). The fish were placed in a recovery tank before returning to their original tank.

The newly developed sampling procedure for mitochondrial assay required several preliminary works. A pilot experiment on a different set of fish evaluated the immunological consequences of biopsy on inflammation resulting from bacterial infection or damaged tissue. There was no significant effect of the biopsy on the immune system (see the supplementary materials and methods for details of all assay protocols and results Fig. S1b). This experiment also demonstrated that biopsy slightly reduced the growth rate but fish continued to gain mass (Fig. S1a). Finally, fish needed anesthesia for the biopsy procedure, and tricaine is the most judicious to study mitochondria in fish muscle. Tricaine acts by blocking voltage-gated sodium channels and blocks neural action potentials²³. Muscle sodium channels are relatively insensitive to tricaine²³, which suggest that tricaine are unlikely to affect the physiology of muscle cells. However, it would be relevant to test for potential effect of tricaine on mitochondrial function in fish muscle.

**Tissue homogenate for mitochondrial assay.** Immediately after biopsy, muscle samples were shredded using micro-dissecting scissors to obtain a homogenous solution with a particle size less than 0.5 mm (tested by pipetting through a 1 mL tip), an homogenization procedure adapted from²⁴. The tissue preparation by shredding used here allowed a rapid and efficient muscle permeabilization, with no loss of tissue, while maintaining the quality of the mitochondria. Validations of the permeabilization method are described in¹¹,²⁵. The homogenates were diluted further in respirometry buffer to obtain a final concentration of 1 mg mL⁻¹ (mean ± s.e = 1.0 ± 0.2 mg mL⁻¹). Red muscle homogenizations were carried out on ice. Only four fish could be run simultaneously for the mitochondrial assay (one measurement per fish), and two runs per day were performed; therefore, five days were required to analyze the forty fish.

**Measurement of mitochondrial metabolic traits.** Oxygen consumption and magnesium green fluorescence were detected simultaneously in four respirometry chambers (2 mL, a chamber per fish) equipped with fluorescent sensors and recorded using DatLab software (Oroboros Instruments, Innsbruck, Austria) at 13 °C with continuous stirring. Immediately after homogenisation and dilution, 2.2 mL of the homogenate was added to one of the four chambers. Mitochondrial metabolic traits were determined using the method described by Chinopoulos, et al.²⁶. First, an adenylate kinase inhibitor (P1,P5-Di(adenosine-5') pentaphosphate pentasodium salt, 0.1 mmol L⁻¹), complex I-linked substrates (pyruvate 5 mmol L⁻¹ and malate 0.5 mmol L⁻¹), and MgGreen (2.2 µmol L⁻¹) were added to each well. After addition of EGTA (0.1 mmol L⁻¹) and EDTA (15 µmol L⁻¹), the fluorescent signal of MgGreen was calibrated using ten successive injections of MgCl₂ (1 mmol L⁻¹). Succinate (10 mmol L⁻¹) was added to supply complex II with the energy substrates. OXPHOS respiration and ATP production were measured by adding ADP (2 mmol L⁻¹). Leak respiration was measured by inhibiting ATP synthesis by adding carboxyattractyloside (4 µmol L⁻¹). The rate of ATP disappearance owing to ATPase activity was also measured under these conditions. This measurement was subtracted from ATP production measured during OXPHOS respiration. An inhibitor of coenzyme Q-cytochrome c reductase (CIII), antimycin A (2.5 µmol L⁻¹), was then added. This measure was subtracted from other mitochondrial respiratory rates to correct for oxygen consumption unrelated to the respiratory chain activity. Finally, mitochondrial density was estimated with COX activity by adding ascorbate (8 mmol L⁻¹) followed by N,N,N,N-tetramethyl-p-phenylenediamine (TMPD; 0.5 mmol L⁻¹) to the chambers. The rate of ATP production and ATP/O ratio were calculated as described in Salin, et al.²⁷. Briefly, we converted the fluorescent signal of the free magnesium concentration to the ATP concentration. We used the same binding affinity (K_d) for Mg²⁺ bound to ATP and ADP as in Thoral, et al.²⁸ because the mitochondrial analyses were performed under the same conditions of temperature and homogenate concentration (K_d-ATP = 0.266 mmol L⁻¹, K_d-ADP = 1.803 mmol L⁻¹). The ATP/O ratio was calculated as the ratio of ATP production to OXPHOS respiration. Rates of oxygen consumption and ATP production were expressed as pmol s⁻¹ mg⁻¹ wet tissue.

In an additional experiment, we determined the technical repeatability of mitochondrial metabolic traits. Duplicated measurements of mitochondrial metabolic traits of two biopsies from the same individual were significantly reproducible (Intraclass correlation coefficients: OXPHOS respiration r = 0.41, *p* < 0.001; LEAK respiration r = 0.48, *p* < 0.001; COX activity: r = 0.40, *p* < 0.001; n = 19; ATP production r = 0.44, *p* = 0.001; n = 14).
<table><caption>**Table 1.** Descriptive analyses of specific growth rates (SGR) and mitochondrial metabolic traits (OXPHOS respiration, ATP production, ATP/O ratio, LEAK respiration, and cytochrome c oxidase [COX] activity) in red muscle of juvenile European sea bass (n = 40).</caption><thead><tr><th></th><th>N</th><th>Min</th><th>Max</th><th>Mean ± SE</th></tr></thead><tbody><tr><td colspan="5"><strong>Past growth rates</strong></td></tr><tr><td>SGR<sub>past</sub> 20 weeks (% day<sup>-1</sup>)</td><td>40</td><td>0.21</td><td>0.30</td><td>0.37 ± 0.08</td></tr><tr><td>SGR<sub>past</sub> 16 weeks (% day<sup>-1</sup>)</td><td>40</td><td>0.22</td><td>0.49</td><td>0.36 ± 0.07</td></tr><tr><td>SGR<sub>past</sub> 11 weeks (% day<sup>-1</sup>)</td><td>40</td><td>0.21</td><td>0.54</td><td>0.34 ± 0.07</td></tr><tr><td>SGR<sub>past</sub> 7 weeks (% day<sup>-1</sup>)</td><td>40</td><td>0.12</td><td>0.42</td><td>0.30 ± 0.07</td></tr><tr><td colspan="5"><strong>Future growth rates</strong></td></tr><tr><td>SGR<sub>future</sub> 4 weeks (% day<sup>-1</sup>)</td><td>40</td><td>0.03</td><td>0.40</td><td>0.26 ± 0.08</td></tr><tr><td>SGR<sub>future</sub> 8 weeks (% day<sup>-1</sup>)</td><td>40</td><td>0.12</td><td>0.37</td><td>0.29 ± 0.05</td></tr><tr><td>SGR<sub>future</sub> 12 weeks (% day<sup>-1</sup>)</td><td>40</td><td>0.20</td><td>0.48</td><td>0.33 ± 0.06</td></tr><tr><td colspan="5"><strong>Mitochondrial metabolic traits</strong></td></tr><tr><td>OXPHOS respiration (pmol O<sub>2</sub> s<sup>-1</sup> mg<sup>-1</sup> red muscle)</td><td>40</td><td>37.16</td><td>98.37</td><td>67.21 ± 15.03</td></tr><tr><td>ATP production (pmol ATP s<sup>-1</sup> mg<sup>-1</sup> red muscle)</td><td>40</td><td>65.02</td><td>387.67</td><td>176.60 ± 63.20</td></tr><tr><td>ATP/O ratio</td><td>40</td><td>0.53</td><td>2.44</td><td>1.28 ± 0.33</td></tr><tr><td>LEAK respiration (pmol O<sub>2</sub> s<sup>-1</sup> mg<sup>-1</sup> red muscle)</td><td>40</td><td>4.27</td><td>11.30</td><td>7.40 ± 1.48</td></tr><tr><td>COX activity (pmol O<sub>2</sub> s<sup>-1</sup> mg<sup>-1</sup> red muscle)</td><td>40</td><td>72.53</td><td>118.01</td><td>92.13 ± 10.46</td></tr></tbody></table>

**Statistical analysis.** Intra-class correlation (ICC) was used to test the consistent differences in growth rates among individuals before and after the mitochondrial assay. Individual consistency of growth rate was tested between the most similar duration of the growth trial, that is, SGR<sub>past</sub> at 7 and 11 weeks was correlated to SGR<sub>future</sub> at 8 and 12 weeks, respectively. Linear mixed models (LMMs) were used to determine the relationships between fish growth rate and mitochondrial metabolic traits. Each SGR (SGR<sub>past</sub> over 7, 11, 16, and 20 weeks, and SGR<sub>future</sub> over 4, 8, and 12 weeks) was run in separate models. The models included SGR as the dependent variable, OXPHOS respiration, ATP production, LEAK respiration, and COX activity as continuous predictors, with the fish tank as a random factor. To control for the effect of initial body mass on growth rate, initial body mass was included as a covariate in each model. The day and the processing run of the mitochondrial assay, as well as the Oroboros device, and the respirometry chamber used for the analysis of mitochondrial assay were initially included as potential random factors in each model but were not significant and were subsequently removed. We also tested whether the degree of mitochondrial efficiency in generating ATP, as determined by the ATP/O ratio, explained individual variation in growth, which had to be examined in a separate linear model to prevent problems associated with multicollinearity between ATP/O, OXPHOS respiration, and ATP production.

Since mitochondrial density can influence mass-specific rates of mitochondrial metabolic traits (Pearson correlations: all r > 0.37 and p < 0.05) and SGR. Estimate of COX-independent mitochondrial traits were used in preliminary LMM analyses of SGR. Residuals from separate linear regressions relating mitochondrial metabolic traits (e.g. OXPHOS and LEAK respiration and ATP production) to COX activity were calculated. As before, we used LMM with SGR as the dependent variable and COX-independent OXPHOS and LEAK respiration, COX-independent ATP production, and COX activity as continuous predictors. Since the patterns of the results of growth rate analyses were the same whether mitochondrial metabolic traits were measured in terms of mass-specific or COX-independent tissues, only those for mass-specific mitochondrial metabolic traits are reported here.

All analyses were based on a sample size of 40 fish and all statistical analyses were performed using R (4.0.3; package lme4; package partR2). Data are presented as means ± SEM, the significance level was set to *p* < 0.05, and inclusive R² (IR²) was added to illustrate the variance explained by a predictor in a model²⁹.

## Results

Individuals varied considerably in their past and future growth rates (Table 1), but as expected, individual variation in growth was significantly repeatable before and after the mitochondrial assay (ICC: all r > 0.31 and p < 0.05). The consistency in the rank order of an individual growth rate between before and after the biopsy procedure shows that the biopsy procedure did not impact the inter-individual differences in growth. Future growth rate was negatively correlated with LEAK respiration (Fig. 2, Table S1). Indeed, variations in LEAK respiration predicted the SGR<sub>future</sub> at 4 weeks (t<sub>34.00</sub> = -2.23, *p* = 0.033, IR² = 0.10; Table S1). The relationship between SGR<sub>future</sub> and LEAK respiration was no longer significant after 8 and 12 weeks of growth (Table S1). Future growth rates were unrelated to OXPHOS respiration, ATP production, or COX activity in the muscle mitochondria (Table S1).

Individuals with the fastest past growth had lower COX activity in their red muscle mitochondria (Fig. 3, Table S2). Indeed, variations in SGR<sub>past</sub> at 7, 11, 16, and 20 weeks were negatively correlated with COX activity (7 weeks: t<sub>31.07</sub> = -3.68, *p* < 0.001, IR² = 0.09; 11 weeks: t<sub>31.96</sub> = -3.98, *p* < 0.001, IR² = 0.10; 16 weeks: t<sub>33.82</sub> = -3.46, *p* = 0.001, IR² = 0.05; 20 weeks: t<sub>33.58</sub> = -3.79, *p* < 0.001, IR² = 0.12; Table S2). Past growth rates were unrelated to OXPHOS respiration, LEAK respiration, or ATP production in the muscle mitochondria (Table S2).
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Future growth rates over 4 weeks was negatively related to the mitochondrial LEAK respiration of red muscle. Future Specific Growth Rates (SGR<sub>future</sub>) as a function of variation in Leak respiration in red muscle of juvenile European seabass (n = 40). SGR<sub>future</sub> were calculated between body mass measured in the mitochondrial assay and posterior body masses. Continuous lines show significant effect (*p* < 0.05) and inclusive R<sup>2</sup> (IR<sup>2</sup>) were added to illustrate the variance explained by LEAK respiration in the variation of SGR<sub>future</sub>. See Table S1 for statistical analyses.</figcaption></figure>

Discussion

We accepted our hypothesis that individual variations in both past and future growth rates are related to mito-
chondrial metabolic traits. However, the mitochondrial traits that explained the growth variation differed
between the growth rates determined before and after the mitochondrial assay. While past SGR was correlated
with COX activity, future SGR was linked to LEAK respiration in the red muscle mitochondria. Nonetheless, as
predicted, individuals with the fastest future growth had lower LEAK respiration rates than fish that grew slowly
after the mitochondrial assay. However, contrary to expectations, individuals with the fastest past growth had
lower COX activity than fish that grew slowly before the mitochondrial assay. We partially accept our hypothesis
that the relationship between growth rate and mitochondrial properties declines over time, with the decline
being especially marked for future growth rates (after four weeks), while that relationship persists over 20 weeks
of past growth rates.

Low LEAK respiration may confer an advantage for future growth. Our fish showed a negative relationship between the rates of LEAK respiration and future growth, calculated as the gain in mass between the time of the mitochondrial assay and four weeks later. This suggests that individuals with higher proton leakage dissipate more energy to offset leakage in place of growth, as previously shown in fish<sup>11</sup>, birds<sup>9</sup> and mammals<sup>30</sup>. But, we also found that this relationship did not persist after eight weeks of growth. The mitochondrial metabolic traits of an individual can change over time<sup>17,31,32</sup>. However, plastic changes are often not similar across individuals<sup>33,34</sup>. Individuals might differ in mitochondrial plasticity simply as a result of changes in mitochondrial traits over time as animals grow or interact in their social environment.

Individual variations in past growth periods were not a function of mitochondrial LEAK respiration. Rather, past SGRs were negatively related to COX activity. This negative relationship is inconsistent with previous findings that individuals with the fastest growth have higher COX activities than fish that grow slowly<sup>6,7,35</sup>. Although the relationships between growth and metabolic traits have generally been studied in the mitochondria of white muscle<sup>6,35</sup>, liver and intestine<sup>35</sup> , and the whole body<sup>36</sup>, one would expect the rank order of an individual’s red and white muscle mitochondria traits, (i.e. its mitochondrial metabolic traits in red and white muscle relative to other individuals in the population) to be consistent between tissues. However, there is also increasing evidence of the absence of covariation in mitochondrial metabolic traits between tissues of the same individual<sup>11,25</sup> as well as evidence of tissue-specific relationships between mitochondrial metabolic traits and whole-animal performance<sup>4,25</sup>. Since our data are correlative, we cannot rule out that intrinsic differences in growth may drive variations in tissue composition if some individuals have more fat reserves and, therefore, lower mass-specific mitochondrial density. An alternative explanation might lie in the fact that activity of the COX is representative of rates of substrate oxidation, where correlation of markers of mitochondrial density can be poor<sup>37</sup> (but see<sup>38</sup>). While measurements of several markers of mitochondrial density were beyond the scope of the present study, citrate synthase activity can be performed from the homogenate after the respirometry protocol. In this way, mitochondrial respiration and CS activity can be determined on the same tissue preparation. Clearly, more research is needed to ascertain whether individual variations in growth are expected to have important consequences for mitochondrial density in some tissues.

Our study questions the idea that mitochondrial metabolism is the underlying mechanism of individual
variation in growth measured before the mitochondrial assay. The differential effect of mitochondrial metabolic
traits on past and future growth is a complex but highly relevant issue, since almost all studies that have exam-
ined the effect of mitochondrial metabolism on animal performance have determined mitochondrial metabolic
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Past growth rates were negatively related to cytochrome c oxidase activity of red muscle mitochondria. Past Specific Growth Rates (SGR<sub>past</sub>) as a function of variation in Cytochrome C Oxidase (COX) activity in red muscle of juvenile European seabass (n=40). SGR<sub>past</sub> were calculated between body mass measured before the mitochondrial assay (at 20, 16, 11, and 7 weeks before) and the body mass at the mitochondrial assay. Duration of the growth trials is indicated on the top right of each panel. For 20 weeks period, plotted represent SGR<sub>past</sub> since initial body mass has not significant effect on SGR. For 16, 11 and 7 week periods, SGR is plotted as partial residuals of SGR<sub>past</sub> evaluated at mean initial body mass (46.6 g, 53.8 g, 60.8 g, respectively) since initial body masses had significant effect on SGR. Continuous lines show significant effect (*p* < 0.05) and inclusive R<sup>2</sup> (IR<sup>2</sup>) were added to illustrate the part of past SGR variation explained by COX activity. See Table S2 for statistical analyses.</figcaption></figure>

traits after, and not before, animal performance. These differential effects occurred despite our experimental animals being of similar age and nutritional state, and having been maintained under identical environmental conditions before and after mitochondrial analyses. Temporal changes in mitochondrial phenotype might occur within weeks³⁹, days¹⁷ and hours³¹. Since the fish growth were determined over 5 months, we cannot rule out that internal clocks of the sea bass, despite relatively constant rearing conditions, played a role in influencing mitochondrial metabolic traits. Little is known about the temporal variation of the mitochondrial phenotype, but Stier, et al.¹⁷ found repeatability in mitochondrial metabolic traits over a period of 10 days in the red blood cells of pied flycatchers (*Ficedula hypoleuca*). Knowledge regarding how mitochondrial traits change over time is currently a crucial step for future research on the relevance of mitochondrial metabolic traits for organismal performance.
Our data emphasize the importance of information regarding the temporal repeatability and variation of mitochondrial metabolic traits. Such data are still lacking because of the need to cull animals when sampling tissues for mitochondrial assays. Although it is possible to make cross-sectional comparisons of animals, these do not reveal individual variations in mitochondrial changes. We and others have partially solved this issue by measuring mitochondrial metabolic traits from biopsy¹⁹ or blood sampling⁴⁰,⁴¹, allowing repeated measurements. However, whether measurements in skeletal muscle or blood cells are representative of mitochondrial metabolic traits in other tissues is unclear, as mitochondrial metabolic traits can differ among tissues²², ⁴² but see⁴³. Combining longitudinal and comprehensive analyses across tissues of mitochondrial metabolic traits will be more effective in increasing the understanding of among- and within-individual variations in whole-organism performance⁴⁴.

## Data availability

The dataset supporting this article is made available in the Supplementary Material.

Received: 17 June 2022; Accepted: 13 September 2022

Published online: 27 September 2022

## References

1. Armstrong, D. P., Keevil, M. G., Rollinson, N. & Brooks, R. J. Subtle individual variation in indeterminate growth leads to major variation in survival and lifetime reproductive output in a long-lived reptile. *Funct. Ecol.* **32**, 752–761. https://doi.org/10.1111/1365-2435.13014 (2018).

2. Sibly, R. M., Kodric-Brown, A., Luna, S. M. & Brown, J. H. The shark-tuna dichotomy: why tuna lay tiny eggs but sharks produce large offspring. *R. Soc. Open Sci.* **5**, 180453. https://doi.org/10.1098/rsos.180453 (2018).

3. Nicholls, D. G. & Ferguson, S. J. in *Bioenergetics (Third Edition)* (eds David G. Nicholls & Stuart J. Ferguson) 249–270 (Academic Press, 2003).

4. Salin, K. et al. Differences in mitochondrial efficiency explain individual variation in growth performance. *Proc. R. Soc. B: Biol. Sci.* **286**, 20191466. https://doi.org/10.1098/rspb.2019.1466 (2019).

5. Bottje, W. G. & Carstens, G. E. Association of mitochondrial function and feed efficiency in poultry and livestock species. *J. Anim. Sci.* **87**, E48–E63. https://doi.org/10.2527/jas.2008-1379 (2009).

6. Nathanailides, C. & Stickland, N. C. Activity of cytochrome c oxidase and lactate dehydrogenase in muscle tissue of slow growing (lower modal group) and fast growing (upper modal group) Atlantic salmon. *J. Fish Biol.* **48**, 549–551. https://doi.org/10.1111/j.1095-8649.1996.tb01449.x (1996).

7. Perello-Amoros, M. et al. Mitochondrial adaptation to diet and swimming activity in gilthead seabream: Improved nutritional efficiency. *Front. Physiol.* https://doi.org/10.3389/fphys.2021.678985 (2021).

8. Brand, M. D. The efficiency and plasticity of mitochondrial energy transduction. *Biochem. Soc. Trans.* **33**, 897–904 (2005).

9. Toyomizu, M. et al. Meat-type chickens have a higher efficiency of mitochondrial oxidative phosphorylation than laying-type chickens. *Comp. Biochem. Physiol. A-Mol. Integr. Physiol.* **159**, 75–81. https://doi.org/10.1016/j.cbpa.2011.01.020 (2011).

10. Salin, K., Luquet, E., Rey, B., Roussel, D. & Voituron, Y. Alteration of mitochondrial efficiency affects oxidative balance, development and growth in frog (*Rana temporaria*) tadpoles. *J. Exp. Biol.* **215**, 863–869. https://doi.org/10.1242/geb.062745 (2012).

11. Salin, K., Auer, S. K., Anderson, G. J., Selman, C. & Metcalfe, N. B. Inadequate food intake at high temperatures is related to depressed mitochondrial respiratory capacity. *J. Exp. Biol.* **219**, 1356–1362. https://doi.org/10.1242/jeb.133025 (2016).

12. Auer, S. K., Salin, K., Rudolf, A. M., Anderson, G. J. & Metcalfe, N. B. Flexibility in metabolic rate confers a growth advantage under changing food availability. *J. Anim. Ecol.* **84**, 1405–1411. https://doi.org/10.1111/1365-2656.12384 (2015).

13. Brown, J. C. L. & Staples, J. F. Mitochondrial metabolism during fasting-induced daily torpor in mice. *Biochimica et Biophysica Acta (BBA) Bioenergetics* **1797**, 476–486. https://doi.org/10.1016/j.bbabio.2010.01.009 (2010).

14. Jorgensen, L. B., Overgaard, J., Hunter-Manseau, F. & Pichaud, N. Dramatic changes in mitochondrial substrate use at critically high temperatures: a comparative study using *Drosophila*. *J. Exp. Biol.* https://doi.org/10.1242/jeb.240960 (2021).

15. Ouillon, N., Sokolov, E. P., Otto, S., Rehder, G. & Sokolova, I. M. Effects of variable oxygen regimes on mitochondrial bioenergetics and reactive oxygen species production in a marine bivalve, *Mya arenaria*. *J. Exp. Biol.* https://doi.org/10.1242/jeb.237156 (2021).

16. Le Roy, A., Mazue, G. P. F., Metcalfe, N. B. & Seebacher, F. Diet and temperature modify the relationship between energy use and ATP production to influence behavior in zebrafish (*Danio rerio*). *Ecol. Evol.* **11**, 9791–9803. https://doi.org/10.1002/ece3.7806 (2021).

17. Stier, A., Bize, P., Hsu, B. Y. & Ruuskanen, S. Plastic but repeatable: rapid adjustments of mitochondrial function and density during reproduction in a wild bird species. *Biol. Lett.* **15**, 5. https://doi.org/10.1098/rslb.2019.0536 (2019).

18. Granata, C., Oliveira, R. S., Little, J. P., Renner, K. & Bishop, D. J. Training intensity modulates changes in PGC-1α and p53 protein content and mitochondrial respiration, but not markers of mitochondrial content in human skeletal muscle. *Faseb J.* **30**, 959–970. https://doi.org/10.1096/fj.15-276907 (2016).

19. Porter, C., Reidy, P. T., Bhattarai, N., Sidossis, L. S. & Rasmussen, B. B. Resistance exercise training alters mitochondrial function in human skeletal muscle. *Med. Sci. Sports Exerc.* **47**, 1922–1931. https://doi.org/10.1249/mss.0000000000000605 (2015).

20. Jeffries, K. M. et al. The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes. *Comp. Biochem. Physiol. B: Biochem. Mol. Biol.* **256**, 110629. https://doi.org/10.1016/j.cbpb.2021.110629 (2021).

21. Picard, M., Hepple, R. T. & Burelle, Y. Mitochondrial functional specialization in glycolytic and oxidative muscle fibers: tailoring the organelle for optimal function. *Am. J. Physiol.-Cell Physiol.* **302**, C629–C641. https://doi.org/10.1152/ajpcell.00368.2011 (2012).

22. Leary, S. C. et al. Fiber-type differences in muscle mitochondrial profiles. *Am. J. Physiol.-Regulat. Integr. Compar. Physiol.* **285**, R817–R826. https://doi.org/10.1152/ajpregu.00058.2003 (2003).

23. Attili, S. & Hughes, S. M. Anaesthetic tricaine acts preferentially on neural voltage-gated sodium channels and fails to block directly evoked muscle contraction. *Plos One* https://doi.org/10.1371/journal.pone.0103751 (2014).

24. Kuznetsov, A. V. et al. Evaluation of mitochondrial respiratory function in small biopsies of liver. *Anal. Biochem.* **305**, 186–194. https://doi.org/10.1006/abio.2002.5658 (2002).

25. Salin, K. et al. Variation in metabolic rate among individuals is related to tissue-specific differences in mitochondrial leak respiration. *Physiol. Biochem. Zool.* **89**, 511–523. https://doi.org/10.1086/688769 (2016).

26. Chinopoulos, C., Kiss, G., Kawamata, H. & Starkov, A. A. in *Methods in Enzymology* Vol. 542 (eds Galluzzi Lorenzo & Kroemer Guido) 333–348 (2014).

27. Salin, K. et al. Simultaneous measurement of mitochondrial respiration and ATP production in tissue homogenates and calculation of effective P/O ratios. *Physiol. Rep.* **4**, e13007 (2016).

28. Thoral, E., Roussel, D., Chinopoulos, C., Teulier, L. & Salin, K. Low oxygen levels can help to prevent the detrimental effect of acute warming on mitochondrial efficiency in fish. *Biol. Lett.* **17**, 20200759. https://doi.org/10.1098/rsbl.2020.0759 (2021).
29. Stoffel, M. A., Nakagawa, S. & Schielzeth, H. partR2: partitioning R-2 in generalized linear mixed models. PeerJ https://doi.org/10.7717/peerj.11414 (2021).

30. Caldeira da Silva, C. C., Cerreira, F. M., Barbosa, L. F., Medeiros, M. H. & Kowaltowski, A. J. Mild mitochondrial uncoupling in mice affects energy metabolism, redox balance and longevity. Aging Cell 7, 552–560 (2008).

31. Hoppel, F. et al. Effects of ultramarathon running on mitochondrial function of platelets and oxidative stress parameters: A pilot study. Front. Physiol. https://doi.org/10.3389/fphys.2021.632664 (2021).

32. Jacques, M. et al. Mitochondrial respiration variability and simulations in human skeletal muscle: The Gene SMART study. *Faseb J* **34**, 2978–2986. https://doi.org/10.1096/fj.201901997RR (2020).

33. Killen, S. S., Adriaenssens, B., Marras, S., Claireaux, G. & Cooke, S. J. Context dependency of trait repeatability and its relevance for management and conservation of fish populations. *Conserv. Physiol.* **4**, 19. https://doi.org/10.1093/conphys/cow007 (2016).

34. Mathot, K. J. & Dingemanse, N. J. Energetics and behavior: unrequited needs and new directions. *Trends Ecol. Evol.* **30**, 199–206. https://doi.org/10.1016/j.tree.2015.01.010 (2015).

35. Eya, J. C., Ashame, M. F., Pomeroy, C. F., Manning, B. B. & Peterson, B. C. Genetic variation in feed consumption, growth, nutrient utilization efficiency and mitochondrial function within a farmed population of channel catfish (*Ictalurus punctatus*). *Comp. Biochem. Physiol. B: Biochem. Mol. Biol.* **163**, 211–220. https://doi.org/10.1016/j.cbpb.2012.05.019 (2012).

36. Mathers, E. M., Houlihan, D. F., McCarthy, I. D. & Burren, L. J. Rates of growth and protein synthesis correlated with nucleic acid content in fry of rainbow trout, *Oncorhynchus mykiss* - effects of age and temperature. *J. Fish Biol.* **43**, 245–263. https://doi.org/10.1111/j.1095-8649.1993.tb00426.x (1993).

37. Scott, G. R., Guo, K. H. & Dawson, N. J. The mitochondrial basis for adaptive variation in aerobic performance in high-altitude deer mice. *Integr. Comparat. Biol.* https://doi.org/10.1093/icb/icy056 (2018).

38. Larsen, S. et al. Biomarkers of mitochondrial content in skeletal muscle of healthy young human subjects. *J. Physiol.-London* **590**, 3349–3360. https://doi.org/10.1113/jphysiol.2012.230185 (2012).

39. Bouchard, P. & Guderley, H. Time course of the response of mitochondria from oxidative muscle during thermal acclimation of rainbow trout Oncorhynchus mykiss. *J. Exp. Biol.* **206**, 3455–3465. https://doi.org/10.1242/jeb.00578 (2003).

40. Braganza, A., Annarapu, G. K. & Shiva, S. Blood-based bioenergetics: An emerging translational and clinical tool. *Mol. Aspects Med.* https://doi.org/10.1016/j.mam.2019.100835 (2020).

41. Stier, A. et al. How to measure mitochondrial function in birds using red blood cells: A case study in the king penguin and perspectives in ecology and evolution. *Methods Ecol. Evol.* **8**, 1172–1182. https://doi.org/10.1111/2041-210X.12724 (2017).

42. Fernández-Vizarra, E., Enriquez, J. A., Pérez-Martos, A., Montoya, J. & Fernández-Silva, P. Tissue-specific differences in mitochondrial activity and biogenesis. *Mitochondrion* **11**(1), 207–213. https://doi.org/10.1016/j.mito.2010.09.011 (2011).

43. Tyrrell D. J., Bharadwaj M. S., Jorgensen M. J., Register, T. C. & Molina, A. J. Blood cell respirometry is associated with skeletal and cardiac muscle bioenergetics: Implications for a minimally invasive biomarker of mitochondrial health. *Redox Biology* **10**, 65–77. https://doi.org/10.1016/j.redox.2016.09.009 (2016).

44. Rebecca, E. et al. Integrating Mitochondrial Aerobic Metabolism into Ecology and Evolution. *Trends in Ecology & Evolution* **36**(4), 321–332. https://doi.org/10.1016/j.tree.2020.12.006 (2021).

## Acknowledgements

We thank Christophe Lebigre and Hélène Ollivier for their suggestions on an earlier version of the manuscript, and two referees for helpful comments that improved the manuscript.

## Author contributions

J.B.Q., S.C., J.C., M.D. and K.S. conceived of the ideas and designed the methodology. J.B.Q., J.C. and M.D. collected data. J.B.Q., J.C., M.D. and K.S. analyzed the data. J.B.Q., M.D. and K.S. wrote the manuscript. J.L.Z.I. revised the manuscript and added comments. All authors approved the final manuscript for publication.

## Funding

This work was supported by the Ifremer and Brittany regions to KS.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-20428-9.

**Correspondence** and requests for materials should be addressed to K.S.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022