Standard metabolic rate does not associate with age‐at‐maturity genotype in juvenile Atlantic salmon

Abstract Atlantic salmon (Salmo salar) is a species with diverse life‐history strategies, to which the timing of maturation contributes considerably. Recently, the genome region including the gene vgll3 has gained attention as a locus with a large effect on Atlantic salmon maturation timing, and recent studies on the vgll3 locus in salmon have indicated that its effect might be mediated through body condition and accumulation of adipose tissue. However, the cellular and physiological pathways leading from vgll3 genotype to phenotype are still unknown. Standard metabolic rate is a potentially important trait for resource acquisition and assimilation and we hypothesized that this trait, being a proxy for the maintenance energy expenditure of an individual, could be an important link in the pathway from vgll3 genotype to maturation timing phenotype. As a first step to studying links between vgll3 and the metabolic phenotype of Atlantic salmon, we measured the standard metabolic rate of 150 first‐year Atlantic salmon juveniles of both sexes, originating from 14 different families with either late‐maturing or early‐maturing vgll3 genotypes. No significant difference in mass‐adjusted standard metabolic rate was detected between individuals with different vgll3 genotypes, indicating that juvenile salmon of different vgll3 genotypes have similar maintenance energy requirements in the experimental conditions used and that the effects of vgll3 on body condition and maturation are not strongly related to maintenance energy expenditure in either sex at this life stage.


| INTRODUC TI ON
Due to their close links with fitness, life-history traits such as age at maturity, offspring number, and size have been a long-term focus of research in biology (Stearns, 1989). Recent advances in genomics have enabled the identification of genes associated with such traits in an increasingly broad range of species (Barson et al., 2015;Lamichhaney et al., 2016;Narum et al., 2018;Troth et al., 2018).
However, follow-up studies for determining the mechanistic basis underlying such associations remain rare, and so, a thorough understanding of the underlying biology, as well as the wider consequences of these genes, is lacking.
Age at maturity is a key life-history trait as it is the source of a trade-off between current and future reproduction (Roff, 2001;Stearns, 1989), and the different life-history strategies of the Atlantic salmon (Salmo salar, Linnaeus, 1758) are a good example of this trade-off. A considerable proportion of Atlantic salmon growth occurs during their marine migration, after which they return to freshwater as sexually mature fish. The time spent at sea (sea age at maturity) varies from 1 to 5 years, with an individual's size roughly doubling with every extra year spent feeding at sea (Hutchings & Jones, 1998). Spending more time at sea gives a potential fitness advantage, as larger, later-maturing individuals have a higher reproductive success (Fleming, 2011;Fleming & Einum, 2011;Mobley et al., 2020). However, the increased time at sea also comes with a higher risk of dying before first reproduction, resulting in a trade-off between size at maturity and risk of mortality (Stearns, 2000).
Analyzing the genetic background for this variation, a genomewide association study (GWAS) covering 57 European Atlantic salmon populations identified a large-effect locus on chromosome 25 that is significantly associated with age at maturity, explaining 39% of the variation in this trait (Barson et al., 2015). The strongest candidate gene in this region was the vestigial-like family member 3 (vgll3) gene. The two alleles of this locus, E and L, associated with either early or late age at maturity, respectively, as well as shorter or longer body length (in equally aged returning individuals).
Taking studies of this gene to the laboratory, examining body condition and maturation probability in a common garden setting, Debes et al. (2021) found that the vgll3*E allele associated with a higher body condition in males and females (but not body length), as well as a higher probability of male maturation in the first autumn after hatching. Interestingly, Debes et al. (2021) also demonstrated a connection between vgll3 and body condition that is independent of maturation status. At the molecular level, expression of vgll3 in Atlantic salmon has been linked with activation of the intracellular HIPPO signaling pathway, suggesting a role in inhibiting gonad development and promoting adipocyte differentiation (Kjaerner-Semb et al., 2018;Kurko et al., 2020;Verta et al., 2020). Furthermore, Ahi et al. (2021) reported strong regulatory effects of vgll3 genotype on the expression of a number of reproductive axis genes in juvenile male Atlantic salmon, indicating that the broad effects of vgll3 can be mediated via the endocrine system. The findings on the mechanistic basis of the vgll3-maturation pathway so far indicate that different vgll3 genotypes are linked with differing patterns of resource allocation or assimilation, where the vgll3*E allele is positively influencing juvenile body condition, enabling earlier gonadal development in males. However, many aspects of the complete molecular and physiological pathways for this process remain unclear.
Early-life performance in acquiring and efficiently using resources sets the stage for the timing of maturation: sexual maturation and reproduction are costly processes, and a key part of juvenile development is the accumulation and storage of excess energy to support future maturation (Hutchings & Jones, 1998;Simpson, 1992;Thorpe, 2007). Standard metabolic rate (SMR) is a potentially important trait with regards to this as it relates to the organism's maintenance energy requirements, capacity for nutrient assimilation, and responses to changing resource availability Auer, Bassar, et al., 2020;Bochdansky et al., 2005;Millidine et al., 2009;Rosenfeld et al., 2015). The effect of vgll3 genotype on body condition as well as its strong effect on age at maturity in Atlantic salmon raises the question of whether vgll3 might assert some of its effects via SMR, warranting an investigation on the potential influence of vgll3 on metabolic phenotype. If high SMR confers an improved ability to grow and exploit resources, or if it is indicative of a fast pace of life, we might expect high-SMR individuals to mature at an earlier time than low-SMR individuals . If vgll3 was associated with SMR, it could be one of the traits through which this gene mediates its effect on maturation timing: with early-maturing vgll3*EE individuals being predicted to have a higher metabolic rate than late-maturing vgll3*LL individuals. Additionally, beyond the effects on growth, body condition, and gonad development, little is known about the effects of vgll3 on a broader range of physiological systems. Given the range of traits found to associate with SMR, testing for an association between SMR and vgll3 could help in refining the potential traits and physiological processes in the scope of vgll3's effects, such as differences in digestive capacity (Millidine et al., 2009), enzyme activity (Norin & Malte, 2012), mitochondrial leak respiration (Salin et al., 2016(Salin et al., , 2019, and behavior (Binder et al., 2016;Biro & Stamps, 2010;Cutts et al., 1998;Metcalfe et al., 1995;Yamamoto et al., 1998).
To further our understanding of the mechanistic functions of vgll3, in this study we test if there is an association between SMR and vgll3 genotype in a common-garden environment using 150 juvenile Atlantic salmon (69 males and 81 females) from 14 families with either vgll3*EE or vgll3*LL genotypes. Based on (1) earlier studies implicating vgll3 as a resource allocation locus and (2) the association among a high SMR, pace of life, and growth performance in high-resource environments, we hypothesized that the difference in resource allocation and maturation timing between the vgll3 genotypes could stem from a difference in metabolic rate. From this, we predicted that individuals with an early-maturation genotype (vgll3*EE) would have a higher SMR than those with the latematuration genotype (vgll3*LL).

| Experimental animals and husbandry
The experiment was conducted under an animal experiment permit granted by the Finnish Project Authorisation Board (permit nr. ESAVI/4511/2020). The cohort of Atlantic salmon used in this study was established using parental individuals deriving from a first-generation hatchery broodstock of salmon originating from the river Kymijoki in Finland, managed by the Natural Resources Institute Finland (LUKE) at their hatchery in Laukaa. Eggs and milt were collected and transferred to the experimental facilities at the University of Helsinki in October 2019 and parents were crossed to create 14 vgll3 homozygote full-sib families (seven vgll3*EE families and seven vgll3*LL families). We used homozygous families to maximize sample sizes of individuals with genotypes expected to exhibit the largest differences associated with vgll3, for example, in SMR.
Fertilized eggs were incubated in darkness in vertical incubators in replicated, family-specific compartments, with a water temperature of 7.2°C until March 2020.
After hatching, on March 6 th , 2020, the alevins were transferred to the experimental facilities at Lammi Biological Station (61°04′45″N, 025°00′40″E, Lammi, Finland), several weeks before they commenced independent feeding. Each family was reared in a randomly selected separate circular 165-L tank (90 cm diameter) and supplied with a continuous flow of UV filtered water from the local lake Pääjärvi, warmed by 1°C through a heat-exchange system.
The incoming water was directed to create a slow circular flow in the tank. The photoperiod was adjusted according to the local latitude for the entire experiment. Initially, fish were fed eight times/day with commercial 0.2 mm pellet food (Vita, Veronesi, Italy), that is, ad libitum. As fish grew, the 0.2 mm feed was gradually replaced with 0.5 mm pellet food at an increasing frequency (up to 12 times/day).
The temperature was 4.5°C when the transfer to Lammi Biological Station was made, and rose gradually from an average of 4.7°C in March to 11.5°C in July (Figure 1). The average temperature during this period was 7.38°C. Tanks were flushed of uneaten food daily, and tank surfaces were carefully scrubbed clean of dirt and algae once or twice each week (depending on temperature). Tanks were checked for dead individuals daily; the mortality rate throughout the study was 3.9%. The specific vgll3 genotype of each family was not known to people participating in fish husbandry, or to those conducting the respirometry experiments, and was only revealed after the respirometry data had been finalized.

| Overall respirometry procedure
Starting on July 4 th , 2020 (at about 2150 degree days), one batch of 16 fish was tested using 16 respirometers (Appendix 1- Figure A1 submerged in the acclimation tank. Air exposure of fish was avoided both during capture from holding tanks and when moving them to respirometers by transferring fish in small plastics cups filled with water. A detailed description of the respirometry holding system is in Appendix 1-materials and methods.

| Respirometry details
Standard metabolic rate was measured in fish at rest using intermittent flow respirometry (Chabot et al., 2016;Forstner, 1983;Steffensen et al., 1984). Respirometry details and results are included as recommended by Killen et al. (2021) and a detailed description of the SMR definition, protocol, and respirometry system is in Appendix 1-materials and methods. Fish were put in their respirometers around 12:00 at noon, and the procedure was stopped at 7:00 in the morning the next day (total time 19 h).
Oxygen concentrations were logged for the entire time the fish were in the respirometers, and data from the period of acclimation were removed during analysis by the MLND calculation (Chabot et al., 2016; See extraction of respirometry data in Appendix 1-materials and methods). Temperature was kept the same as in the acclimation tanks at 10.53 ± 0.02°C (mean ± SD) during respirometry. The respirometry tank was covered with a light-blocking tarp to leave the respirometers in complete darkness, making the fish unable to see each other during the procedure. The respirometry system was in a separate room from other activities in the research building to avoid disturbances.
Respirometers were set to cycle between the open flush phase and the closed measurement phase in 5-(flush) and 15-min (measurement) intervals, resulting in a total of 50-60 measurement cycles for each batch. Sampling frequency of oxygen concentration F I G U R E 1 Average daily water temperature in holding tanks at Lammi Biological Station. Water was supplied from the local lake Pääjärvi, and temperature thus fluctuated according to the temperature of the lake. Average temperature of the holding tanks during the respirometry period and the 3 days leading up to it was 10.34 ± 0.77°C (SD). All tested fish were held at a constant temperature of 10.5°C in acclimation tanks for 40 h before respirometry was set to one sample every 2 s. The average minimum oxygen concentration reached during the closed phase was 9.52 ± 0.76 mg/L (±SD), equating an average drop in dissolved oxygen of 12 ± 7% (±SD). For the entire experiment, the lowest overall oxygen concentration reached was 4.22 mg/L, and there were 10 cycles in total where the oxygen concentration reached below 6 mg/L. The insides of respirometers were brushed daily to minimize bacterial or algal growth, and oxygen probes were gently wiped with a paper tissue moistened with a 70% ethanol solution before the measurements. After half of the batches were completed, the oxygen probes were recalibrated and the entire system (except acclimation tanks) was cleaned and disinfected with a solution of bleach, followed by thorough rinsing with fresh water. Background respiration was recorded for each batch of measurements before and after putting the fish in the respirometer. This was done individually for each chamber in each batch. The length of each background respiration measurement was 15 min (one closed cycle).
The average background respiration was 7.48 ± 4.55% (±SD) of the total respiration in the chamber with fish.
After each respirometry batch, fish were removed from their respirometers and euthanized (sodium bicarbonate-buffered methanesulfonate overdose, 250 mg/L), then carefully dried on each side using a tissue paper before their body mass was measured to the nearest 0.01 g using a precision scale (Scout STX222; Ohaus). A fin clip was taken from the caudal fin and stored in ethanol for verifying the vgll3 genotype and for sex determination. We used Kompetitive allele-specific polymerase chain reaction (KASP™, LGC, UK) assays (He et al., 2014) for the vgll3 TOP SNP and the sex-specific SDY locus as described in Sinclair-Waters et al. (2021). At the time of respirometry testing, the body mass distributions of the vgll3*EE and vgll3*LL families were slightly different. This was caused by an earlier experiment utilizing the same study population, which sampled individuals non-randomly in an attempt to size-match individuals from different families, causing significant body mass difference between vgll3 genotypes in this experiment that is not necessarily caused by the vgll3 genotype. As a consequence of this, no statistical analysis on body mass differences was done, and all analyses controlled for body mass.

| Statistical analysis
To test if vgll3 genotype associated with a change in standard metabolic rate, we used the following linear mixed model (model 1): where log 10 SMR (mg O 2 /L/h) was fitted against log 10 body mass (g; accounting for allometric scaling), including vgll3 genotype (LL/EE), sex (m/f), and the interaction between vgll3 genotype and sex as fixed effects, thus including both the direct effects of vgll3 genotype and sex on standard metabolic rate, as well as a potential sex difference in the effect of vgll3. Family (1-14) was included in the model as a random effect, accounting for variation caused by both tank and family effects as they were confounded in the experiment (one family per tank). Respirometry batch (1-10) was also included as a random effect.
To control for potential temperature effects on metabolic rate via thermal acclimation, the average holding tank temperature (°C) for the last 3 days before the fish entered the acclimation tank was included as an additional fixed effect.
Log 10 SMR ∼ Log 10 BodyMass + vgll3 + sex + vgll3:Sex + AvgTemp + Family + Batch F I G U R E 2 Overall respirometry procedure and structure. Fish were reared in 14 tanks, separated by family, with each family having either the vgll3*EE or vgll3*LL genotype. Fish were tested in 10 batches of 16, with an acclimation period of 40 h (constant temperature, no feed) before being held in respirometers for about 19 h while oxygen consumption was being measured To test if there was an association between vgll3 and body condition, we used a similar approach as above, where we modeled body condition using the following mixed-effect model (model 2): where condition is modeled as a function of vgll3 genotype, sex, and the interaction between vgll3 genotype and sex (fixed effects), as well as family and batch (random effects). Body condition was calculated as Fulton's K, which is the body weight (g) divided by the cube of the body length (mm) multiplied by 100, and represents the relationship between the body weight and length (Ricker, 1975).
p-and F-test values were calculated for fixed effects using Satterthwaite's method with type III tests in models with interactions, and type II tests in models without interactions. The full models were compared to three alternative models: A no interaction model removed the interaction between sex and vgll3 genotype; a sex-only model further removed vgll3 genotype as a fixed factor, and a simple model removed both sex and vgll3 genotype. These models were compared to the full model using Akaike's information criteria (AIC) and by testing the predictive power of each model by cross-validation. Cross-validation was performed by separating the dataset into two sets of either even-or odd-numbered respirometry batches, and then using the model parameters derived from each set try predicting the SMR of the other; the root-mean-square error was then used to compare the predicted values against the observed ones for each model. All models were visually examined and confirmed to have normally distributed residuals. Statistical summaries of both sets of models are shown in Tables A1 and A2 in Appendix 2-analysis and data.
Nine individuals were removed from the analysis because of technical issues (pump, oxygen logger, or identification failure), and one individual that entered acclimation was not tested because of unexpected mortality, leaving the total number of individuals used in analysis at N = 150.

| RE SULTS
A total of 150 fish were successfully measured for SMR, of which 76 had the vgll3*LL genotype (34 males and 42 females) and 74 had the vgll3*EE genotype (35 males and 39 females; Appendix 2- Table   A3). The number of individuals tested was similar for each family, with an average of 10  Table A2).
The vgll3 genotype had no significant effect on SMR ( Figure 3).
There was also no significant effect of sex (p = .59) or the sex:vgll3 interaction (p = .69). The models excluding vgll3, the vgll3:sex interaction, or both sex and vgll3 as fixed factors, all had a lower AIC and root-mean-square deviation (from cross-validation) than the full model, indicating that the models excluding vgll3 and sex had a higher relative quality and higher predictive power. The metabolic scaling coefficient was estimated to 0.88 ± 0.05 (95% CI) in the full model and 0.89 ± 0.24 in the simple model, and the coefficient of variation was 6.23% and 6.19% for the simple and full model, respectively. The coefficient of variation was calculated as the SD of residuals after scaling them from log 10 to log e scale, as suggested by Garland (1984). See Appendix 2- Table A1 for the full statistical summary.
Condition ∼ vgll3 + sex + vgll3: Sex + Family + Batch F I G U R E 3 Main results. (a) Log 10 -log 10 plot of standard metabolic rate (mg O 2 /L/h) and body mass (g) for every included individual (N = 150). The square and circle with a black border indicate genotype means, and horizontal error bars indicate the standard error for mean body mass. (b) Estimates of model coefficients for the full mixed linear model of log 10 standard metabolic rate (model 1). The full model summary is shown in Appendix 2- Table A1. (c) Predicted means of mass-specific standard metabolic rate (mg O 2 /L/kg/h) for male and female individuals with either vgll3*EE or vgll3*LL genotypes, fixed for a body mass of 1.08 g and an average temperature of 10.28°C (preceding the respirometry acclimation and trial) at the population level, based on the model estimates in b Using the residual standard deviation of 0.027 obtained from the full model, a post hoc power analysis was performed to estimate the minimum effect size this study should have been able to detect with a power of 80%. Repeatedly simulating datasets (Appendix 2- Figure A2) where SMR is affected by vgll3 genotype and sex as specified in the full model-and testing different values for the effect of vgll3 and the vgll3:sex interaction -, we found a minimum effect size of 4.25% and 6% on mass-adjusted (non-log) SMR for the vgll3 and the vgll3:sex effects, respectively (translating to 0.018 and 0.025 for the model estimates). The power analysis indicates that although we did not detect a vgll3 effect in this experiment, we cannot rule out effects of vgll3 and vgll3:sex on massadjusted SMR smaller than 4.25% and 6%.

| DISCUSS ION
Vgll3 genotype has a strong effect on Atlantic salmon life history (Ayllon et al., 2015;Barson et al., 2015;Debes et al., 2021;Verta et al., 2020), and recent findings have suggested that this effect might be mediated by vgll3-induced differences in body condition (Debes et al., 2021). Taking this together with the connections between standard metabolic rate (SMR) and resource acquisition and assimilation Auer, Bassar, et al., 2020;Auer, Solowey, et al., 2020;Bochdansky et al., 2005;Millidine et al., 2009;Rosenfeld et al., 2015), we hypothesized that vgll3 genotype might be asserting some of its effects through changes in energetic physiology that would be reflected in SMR. However, we did not detect any difference in either SMR or body condition between juvenile Atlantic salmon individuals of different homozygous vgll3 genotypes. The lack of difference in body condition might be related to body size or life stage; the fish in the study of Debes et al. (2021) were considerably more developmentally advanced, averaging 19 g, and had been reared for approximately 3100 degree days at stable daily temperatures with many of the males maturing the same autumn. In comparison, the fish in this study averaged just 1.08 g and had been reared at approximately 2150 degree days under a variable temperature regime, which could have further reduced growth (Imholt et al., 2011). It might then be possible that the body condition effects of vgll3 genotype only start to occur later in development when the fish has reached larger body mass, or in the months prior to maturation. Considering the SMR results, measurement error is one factor that could potentially mask vgll3 effects, but the SMR measurements in this study were generally precise, with a metabolic scaling coefficient well within the range for what has been observed in other fishes (Clarke & Johnston, 1999) and a low coefficient of variation for the model's residuals, indicating good precision of the results. Our post hoc power simulation based on our observed residual variation indicated that this experiment, with its comparably high number of individuals, should have been able to detect effects of vgll3 genotype down to an effect size as small as 4.25% on mass-adjusted SMR, as well as a sex-by-vgll3 genotype interaction effect down to 6% (Appendix 2- Figure A2). This would be a comparably weak effect considering that mass-and sex-adjusted metabolic rates among individuals have been reported to vary two to threefold within a population Burton et al., 2011;Metcalfe et al., 2016). Given this, the lack of difference in SMR prior to maturation or prior to a change in body condition does not support SMR of juvenile salmon as a mediator of vgll3 effects on age at maturity at this particular age and in these conditions. Finally, these results were supported by a parallel study using a different crossing design and juvenile salmon of larger body size (Prokkola et al., 2021) where no association between SMR and vgll3 genotype was found.
The results indicate that any potential organismal differences caused by vgll3 genotype do not affect physiological function or body composition in a way that changes SMR at this life stage in these conditions. Such differences could, for example, have been organ size (Rosenfeld et al., 2015), enzyme activity (Norin & Malte, 2012), mitochondrial leak respiration (Salin et al., 2016(Salin et al., , 2019, or capacities related to digestion or resource acquisition/assimilation (Auer et al., 2015;Millidine et al., 2009). Alternatively, but less parsimoniously, there may be tissue-dependent variation in vgll3-effects on metabolic rates which cancel each other out at the organismal level, resulting in similar whole-animal metabolic rates. For example, vgll3*EE individuals could have more metabolically active digestive systems (allowing for faster nutrient assimilation) or increased muscle mass (high MR), but then a higher investment in adipose tissue (low MR) masks this difference's effect on the mass-specific metabolic rate.
Nevertheless, our results indicate that energy-allocation effects of vgll3 are not driven by transient whole-animal maintenance energy requirements at this particular life stage and in these conditions. Metabolic phenotype is a multifaceted set of traits, covering not just standard metabolic rate, but also the maximum metabolic rate, aerobic scope, and daily energy expenditure, to mention some.
Additionally, the consequences of variation in metabolic traits are context dependent and may change under different environmental conditions and life stages (Auer, Bassar, et al., 2020;Auer, Solowey, et al., 2020;Bochdansky et al., 2005;Millidine et al., 2009;Norin & Metcalfe, 2019). For example, compared to a natural setting, the ad libitum feeding is a key feature of the environment experienced by the individuals used in this study until the 2 days fasting before the SMR measurements. The benefit of a high metabolic rate may depend on resource availability (Auer, Bassar, et al., 2020), which in the wild will be more variable both spatially and temporally compared to this experiment. We initially hypothesized that vgll3 might influence metabolic rate such that vgll3*EE individuals have a higher SMR, benefiting them under high-resource conditions. This is not what we found, and our findings rather suggest that in ad libitum (highresource) conditions, vgll3 does not influence growth via changes in SMR. Future research on this, and related topics, would benefit from investigations in low-or fluctuating resource availability environments in order to assess potential effects via, for example, phenotypic plasticity in metabolic rate (Auer, Solowey, et al., 2020;Norin & Gamperl, 2018).
We initially planned to also measure metabolic rate under exhaustive swimming to obtain data on maximum metabolic rate (MMR), but had to abandon these trials as we were unable to motivate | 7 of 14 ÅSHEIM Et al.
juvenile salmon of this age and size to do any exhaustive swimming, either by swim tunnel or by hand chasing. This behavior might be specific to this life stage, as we were eventually successfully able to exercise fish up to MMR by hand chasing in a parallel study using a cohort from the same year that was studied 1-2 summer months later when the fish had grown from 1 to 4 g (Prokkola et al., 2021).
Besides supporting our finding on the lack of an association between SMR and vgll3 genotype, Prokkola et al. (2021) found a significantly higher MMR and aerobic scope in vgll3*EE compared to vgll3*LL individuals, indicating that vgll3 genotype does affect the metabolic phenotype, but does so by affecting traits related to MMR without affecting SMR. In the wild, juvenile Atlantic salmon prefer staying in slow-velocity microhabitats where they act as sit-and-wait predators, darting out to catch suitable prey items as they pass by (Fraser et al., 1993;Metcalfe et al., 1997). It is thus possible that individuals at the life stage studied here are not well physiologically adapted to aerobically exhaustive swimming, but rather shorter bursts, due to their small relative muscle mass and glycogen stores available for exercise. The available literature reporting successful swimming respirometry or MMR measurements of salmon of ~1 g size is very sparse, and the experiments reported by Dabrowski (1986) and Cutts et al. (2002) are to our knowledge the only ones that have been successful in motivating juvenile salmon close to or below this size to swim. We therefore recommend focusing efforts on investigating MMR at the latter end of the early life stage, or carefully devising alternate ways of ethically inducing MMR.
The lack of association observed here between vgll3 genotype and SMR indicates that at this life stage and in these conditions, SMR will not necessarily be co-selected together with vgll3 genotype. Given these results, we do not expect selection on age at maturity to constrain variation in SMR, at least not via selection on vgll3. However, in addition to vgll3, other genes have also been identified as potentially important for maturation in Atlantic salmon, such as six6 (Sinclair-Waters et al., 2020), akap11 (Barson et al., 2015;Kurko et al., 2020), and tead3 (Christensen et al., 2017)(see (Mobley et al., 2021), for an overview). Furthermore, Prokkola et al. (2021) detected an interaction between six6 and vgll3 genotype which associated with a reduction in maximum metabolic rate in individuals with the vgll3*LL six6*LL genotype combination (although no effect of six6 on SMR was detected).
Thus, while no link between age at maturity genotype and metabolic rate was detected in the current study, it remains to be determined whether there are links with other age at maturity genes or metabolic traits. Studying these could be worthwhile to improve our overall understanding of the connection between energetics and genetic control of maturation timing.

| CON CLUS IONS
We found that the vgll3 age-at-maturity genotype did not significantly affect standard metabolic rate in small (0.5-2 g) Atlantic salmon juveniles (21 weeks post-hatch). Our results indicate that vgll3's effect on age at maturity and resource allocation is unlikely to be mediated through maintenance energy requirements or related traits and that juvenile individuals of different vgll3 genotypes face similar maintenance energy requirements at the developmental stage and conditions used in this study.

CO N FLI C T O F I NTE R E S T
The authors declare no competing interests.

O PE N R E S E A RCH BA D G E S
This article has earned an Open Data Badge. Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://zenodo.org/ recor d/5255061. https://doi.org/10.5281/zenodo.5255061.

DATA AVA I L A B I L I T Y S TAT E M E N T
The full datasets and the R scripts used to analyze them are available
Each chamber cap had two outlets, one for circulatory flow and one for flush flow. Inside the chamber, a flat circular 3D-printed plastic baffle partly covered each chamber cap to protect the outlets and to ensure a mixed water flow. The baffles were secured to the chamber caps using a stainless steel screw. The chamber caps and baffles were 3D printed using polyethylene terepthalate glycol (PETG, Devil Design, Mikolow, Poland) Respirometer volume was measured using the weight difference between dry respirometers and respirometers filled with water (while plugging the flush-flow and oxygen probe outlets and inlets).
Respirometers had a total volume of 47.3 ml, and the fish had an average weight of 1.08 ± 0.3 g (±SE) giving a 2.2% ratio of animal body mass to volume of respirometer. Oxygen probes were calibrated in O 2 -free water (deoxygenated using sodium sulphate) and air-saturated water.

Respirometry holding system
The full respirometry system consisted of three square 164 L tanks, two acclimation tanks and one respirometry tank, all located in a different room from the holding tanks. As opposed to the holding tanks, which had variable temperatures depending on the temperature of the incoming lake water, the respirometry system was set up to use the same water but to hold a constant temperature of 10.5°C.
This was achieved by setting up two 8°C cold-water reservoir tanks, which both received an incoming 1 L min -1 flow of filtered (UV and a 60 µm physical filter) lake water, but were cooled by two separate aquarium coolers. Each of the respirometry system's three tanks had a temperature-controlled pump which created a slow flow of 8°C water from the reservoirs into the tank whenever the tank temperature got higher than 10.5°C, achieving a highly stable temperature at 10.53 ± 0.02°C (mean ± SD) for the entirety of the experiment.
An overflow pipe drained water from the acclimation tanks back to the cold-water reservoir. The two acclimation tanks shared one coldwater reservoir while the respirometry tank used the other. Both the respirometry tank and its cold-water reservoir were continuously UV filtered to reduce the growth of algae and background respiration. Water in all three tanks was continuously aerated with airstones, and each acclimation tank had two Eheim 300 pumps (Eheim, Deizisau, Germany) to keep water mixed. Lighting was provided by one 80-cm Juwel Novolux LED-lights (JUWEL Aquarium, Rotenburg, Germany) for each acclimation tank, with photoperiod adjusted according to the local latitude.

F I G U R E A 1 Design of respirometers.
(a) Schematic (not to scale) overview of a single respirometer. (b) One group of four respirometers, in the respirometry tank, with fish, shown from above

Extraction of respirometry data
All oxygen traces were inspected visually for linearity (and blindly from vgll3 genotype), using regression diagnostic plots (i.e., residuals vs fitted, normal Q-Q, and scale-location), and 1387 slopes were excluded, leaving an average of 47.85 ± 5.1 (±SD) MO 2 measurements per individual used in the statistical analysis. Using the r 2 of the fitted traces to inspect for linearity was dropped as only 0.6% of all slopes (54 out of 8739) failed to pass the coefficient of determination threshold (r 2 < .95), despite spontaneous fish activity and occasional malfunction of circulation pumps. MO 2 was calculated from each of these slopes as described by Morozov et al. (2020) using the R package FishResp. Mean r 2 for all oxygen traces used in the statistical analysis was 0.997 ± 0.004 (±SD, N = 7311). For each individual, the standard metabolic rate (SMR) was defined as the mean of the lowest normal distribution (MLND) after having fitted all MO 2 measurements to up to four normal distributions as described by Chabot et al. (2016). This method automatically separates the MO 2 measurements into distributions representing the acclimation time, the measurements where the fish are quiescent, and measurements following potential spontaneous activity. In FishResp (Morozov et al., 2020), background respiration was modelled to increase linearly (based on pilot tests) from the pre-test measurement to the post-test measurement, and oxygen consumption was controlled for background respiration by subtracting the background respiration (and its predicted linear increase) from the oxygen concentration measurements before calculating the rate of oxygen consumed by the fish.

A PPE N D I X 2 A N A LYS I S A N D DATA
F I G U R E A 2 Power analysis. Shows the predicted power for this study (N = 150 individuals, alpha = 0.05, residual sd = 0.027) to detect an effect on SMR from vgll3 genotype and an interaction between vgll3 genotype and sex, given different effect sizes. The analysis gives an indication of the smallest effect sizes this study can detect, given that vgll3 genotype and sex affects SMR as in the model written below. Effect sizes are shown as the degree that the factors affect the mass-adjusted metabolic rate (mg O2 L-1h-1kg-1). For example, a 4% effect size would increase or decrease mass-specific metabolic rate with 4%. Each point represents the proportion of 10000 simulated datasets where a significant effect was found. Simulated datasets are based on a model of the form: log(SMR) = a + b × log(body mass) where a is the intercept; b is the metabolic scaling coefficient, set to 0.8; c and d are the effect sizes of genotype and the interaction between genotype and sex, respectively; and e is the error component, which here is normally distributed with a mean of 0 and a standard deviation of 0.027 (taken from results of this study). In the simulation for effect of genotype, the genotype-sex interaction was fixed at 0%. In the simulation for interaction of genotype and sex, the effect of genotype was fixed at 4%. The simulated datasets used body mass sampled randomly from a normal distribution with a mean of 1.08 g and a standard deviation of 0.3 (taken from results of this study), balancing vgll3 genotype and sex evenly across the 150 simulated individuals. The analysis is available in the uploaded R script (see data accessibility)

Simple model
Log 10 smr∼Log 10 mass+AvgTemp+(1|family)+(1+batch) 0.0656 −606.652 Note: F and p values are calculated using Satterthwaite's method using type III (full model) or type II (simple model) tests. Mixed models are fitted using restricted maximum likelihood. (A) Summary for the full model, including the fixed effect of sex (m/f) and vgll3 genotype (EE/LL), as well as their interaction. (B) Summary for the simple model, excluding fixed effect of sex or vgll3 genotype. (C) Model comparison between the full model, a simple model, no-interaction model and a sex-only model (excluding vgll3 but keeping sex as fixed factor). The Root-mean-square error (RMSE) is the result of a cross-validation of the models, which is obtained by splitting the dataset into two sets of odd-and even-numbered respirometry batches, and using model parameters estimated from one set to predict the results of the other set. The RMSE represents the deviation between the predictions and the real data, and lower values indicate a higher predictive ability. Batch was excluded as a random effect from the models in the cross-validations due to singularity. The Akaike information criterion (AIC) is an estimation of the relative quality of the models, taking residual variation and model complexity into account; A lower value indicates higher relative quality. Note: Body condition is calculated as Fulton's K, where the body condition is the body mass (g) divided by the cube of the body length (mm) and multiplied by 100. F and p values are calculated using Satterthwaite's method using type III (full model) or type II (sex-only model) tests. Mixed models are fitted using restricted maximum likelihood. (A) Summary for the full model, including fixed effect of sex (m/f) and vgll3 genotype (EE/LL), as well as their interaction. (B) Summary for the sex-only model, excluding effect of vgll3 genotype. (C) Model comparison between the full model, a no-interaction model, and sex-only model. The Root-mean-square error (RMSE) is the result of a cross-validation of the models, which is obtained by splitting the dataset into two sets of odd-and even-numbered respirometry batches, and using model parameters estimated from one set to predict the results of the other set. The RMSE represents the deviation between the predictions and the real data, and lower values indicate a higher predictive ability. Batch was excluded as a random effect from the models in the cross-validations due to singularity. The Akaike information criterion (AIC) is an estimation of the relative quality of the models, taking residual variation and model complexity into account; A lower value indicates higher relative quality.