1Vol.:(0123456789) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports Heritable variation in swimming performance in Nile tilapia (Oreochromis niloticus) and negative genetic correlations with growth and harvest weight Samuel Bekele Mengistu1,2,3, Arjan P. Palstra1*, Han A. Mulder1, John A. H. Benzie2,4, Trong Quoc Trinh2, Chantal Roozeboom1 & Hans Komen1 Nile tilapia is predominantly produced in smallholder ponds without aeration. We hypothesize that Nile tilapia with high oxygen uptake efficiency (O2UE) may perform better under these conditions than Nile tilapia with low O2UE. Critical swimming speed (Ucrit, in cm s−1) is a potential indicator for O2UE. Our objectives were to estimate variance components for Ucrit and fish size at swim testing early in life, and genetic correlations (rg) between Ucrit with harvest weight (HW) and daily growth coefficient (DGC) later after grow-out in a non-aerated pond. Substantial heritability was found for absolute Ucrit (0.48). The estimated rg between absolute Ucrit and fish size at testing were all strong and positive (range 0.72–0.83). The estimated rg between absolute Ucrit and HW, and absolute Ucrit and DGC were − 0.21 and − 0.63 respectively, indicating that fish with higher absolute Ucrit had lower growth in the non- aerated pond as compared to fish with lower absolute Ucrit. These results suggest a juvenile trade-off between swimming and growth performance where fish with high Ucrit early in life show slower growth later under conditions of limited oxygen availability. We conclude that Ucrit in Nile tilapia is heritable and can be used to predict growth performance. Nile tilapia (Oreochromis niloticus) is predominantly produced in smallholder tilapia ponds without aeration. In non-aerated ponds, dissolved oxygen (DO) drops below critical level (3 mg l−11) during the night. Low DO in smallholder farms negatively affects Nile tilapia growth2. It may be expected, therefore, that Nile tilapia with high oxygen uptake efficiency may grow better under these conditions than Nile tilapia with low oxygen uptake efficiency. As critical swimming speed (Ucrit) may reflect the oxygen uptake efficiency, the hypothesis is that fish with high Ucrit will grow better under conditions where oxygen is limiting. A high throughput method to assess the individual variation in oxygen uptake efficiency is by subjecting fish to exhaustive exercise in a critical swimming challenge test. In this test, swimming speeds are incrementally increased at prescribed intervals until fish stop swimming and fatigue3,4. Individual fish fatigue when swimming at a specific speed interval for a certain period, from which the Ucrit3 can be determined. Recently we have devel- oped and applied such tests for gilthead seabream (Sparus aurata) and Atlantic salmon (Salmo salar)5. Oxygen uptake is maximal at Ucrit, although the anaerobic component by fast skeletal muscle increases when nearing Ucrit6. Near Ucrit, the metabolic demand for oxygen is becoming greater than can be provided by ventilatory and circulatory systems7. Fish that are able to consume more oxygen can swim faster, or reverse for the connection that we are interested in: faster swimming fish have higher oxygen uptake efficiency. Particularly for tilapia, the link between Ucrit and maximal oxygen consumption may be strong because tilapia has a high Ucrit (4.94 ± 0.45 BL s−1 for ~ 15 cm fish) and a very high maximum metabolic rate8. Hence, Ucrit could be a good indicator of oxygen uptake efficiency of individual tilapia. The heredity of athletic performance has received considerable research attention in dog9, horse10,11 and human12. Genetic parameter estimates for swimming performance in fish are scarce, but suggest that swimming OPEN 1Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH Wageningen, The Netherlands. 2WorldFish, Jalan Batu Maung, Batu Maung, 11960 Bayan Lepas, Penang, Malaysia. 3School of Animal and Range Sciences, College of Agriculture, Hawassa University, P. O. Box 5, Hawassa, Ethiopia. 4School of Biological Earth and Environmental Sciences, University College Cork, Cork, Ireland. *email: arjan.palstra@wur.nl 2Vol:.(1234567890) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ performance has a heritable component. Broad sense heritabilities (i.e. not corrected for dominance and epistatic interaction effects)13 of swimming performance were estimated by Garenc et al.14 in stickleback (Gasterosteus aculeatus), by Hurley and Schom15 in Atlantic salmon and by Nicoletto16 in guppy (Poecilia reticulata). More recently, Vandeputte et al.17 estimated the additive genetic variance component for relative Ucrit (Ucrit divided by standard length) in European sea bass (Dicentrarchus labrax) and found a heritability of 0.55, with a negative genetic correlation with body weight. We therefore aimed first to estimate variance components for swimming performance in Nile tilapia expressed as Ucrit and to estimate the genetic correlation between Ucrit and fish size at swim testing early in life. Next, tested fish were stocked in a non-aerated pond and grown to harvest weight, to determine the genetic correlations between Ucrit early in life and harvest weight (HW) and daily growth coefficient (DGC) later in life. Materials and methods Ethics statement. This study utilised phenotypic data collected as part of the GIFT selective breeding pro- gram managed by WorldFish at the Aquaculture Extension Centre of the Malaysian Department of Fisheries at Jitra, Kedah State, Malaysia (6° 15′ 32° N; 100° 25′ 47° E). This study was approved by the internal WorldFish eth- ics committee. All fish in the GIFT breeding population are managed in accordance with the Guiding Principles of the Animal Care, Welfare and Ethics Policy of WorldFish. Experimental fish. Nile tilapia of the Genetically Improved Farmed Tilapia (GIFT) strain from genera- tion 18 was used in this experiment. The 60 full sib and half sib families were produced using 31 males and 58 females, of which two females were used twice with different males. The planned mating ratios were one male to at least two females. However, the successful matings were: 12 males each mated with one female (resulting in 12 full sib families), 12 males each mated with two females (resulting in 12 half sib groups equivalent to 24 full sib families), 4 males each mated with 3 females (four half sib groups equivalent to 12 full sib families) and 3 males each mated with 4 females (three half sib groups equivalent to 12 full sib families). Each full sib family was reared separately in a hapa (fine mesh net enclosure) set up in an earthen pond. The image analysis was done as described previously by Mengistu et al.18. In total 1500 photographs were loaded into tpsUtil software19 and digitized for six landmarks using tpsDig 2.3020. Landmarks 1 and 2 were on the 0 and 20 cm marks on the ruler which was photographed together with the fish for scaling. The landmarks 3 and 4 were used to measure standard length, the distance between the tip of the snout to the base of caudal fin. The landmarks 5 and 6 were the dorsal and ventral landmarks where the distance is maximum. These landmarks were used to calculate height, the maximum dorso-ventral distance (Fig. 1). To obtain the distance between the Cartesian coordinates, these landmarks were analysed in R software using Geomorph package version 3.0.721 and the true distance in cm was computed based on the reference scale. Swim test experiment. The swim test was done in 30 working days, one swim test per day. Thirty to 35 relatively bigger fingerlings from 60 full sib families were selected, PIT tagged and housed in a tank. Three weeks after PIT tagging, 25 fish in a range from 5 to 10 cm standard length at swim testing (SLtest, in cm) from each of the 60 full sib families were measured using a ruler with a centimetre scale, weighed (Wtest, in g) and photo- graphs were made, one day before the swim test. The SLtest and height at swim testing (Htest, in cm) of the fish used in our analysis were obtained from the photographs of each fish using image analysis. The number of fish tested per family was 25 and the number of fish per test was 50 fish. Therefore, we tested either 10 fish from 5 families or 5 fish from 10 families which resulted in all 25 fish from each family being tested in three consecutive days. To determine the Ucrit, a Brett type (rectangular oval shape raceway) swim flume of 230 cm length and 90 cm width with a water depth of 40 cm was used22. Water current was created using a Minn Kota Terrova 80 lbs propeller. The propeller has 10 speed settings, in this experiment speed levels from 2 to 10 were used. Supple- mentary Table 1 provides the flow speeds measured at each of the settings. As the assessment of Ucrit requires Figure 1. Nile tilapia with landmarks 1:6. Landmarks 1 and 2 marks a reference scale of 20 cm length, landmarks 3 and 4 represent the snout and base of the caudal fin, respectively, landmarks 5 and 6 were used to measure height (maximum dorso-ventral length) of the experimental fish. 3Vol.:(0123456789) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ all fish to fatigue, this experimental set-up could be applied for early life testing at small size and not for older and larger fish. Feeding was stopped 24 h before the beginning of the swim testing. The fish were acclimatised for one hour in the swimming flume without flow. After acclimation, the propeller was turned on to induce swimming at the second setting. The time at each setting was fixed at 30 min and flow increments continued until all fish fatigued. At each setting, the average water flow velocity was recorded using a FP111 Global Water Flow Probe (FP111, Global Water, USA). The swim test could take maximally 4.5 h, with 9 propeller speed levels. A fish fatigued when it touched the back fence and could not be stimulated to continue swimming. Each fatigued fish was scooped out immediately and PIT tag number and time at fatigue were recorded. The mean DO in the tank just before resuming the swim test was 5.6 ± 0.4 mg l−1 (71.2% saturation), ranging from 4.9 to 6.4 mg l−1, and during the swim test it was 7.6 ± 0.4 mg l−1 (97.6% saturation), ranging from 6.4 to 8.8 mg l−1. The mean water temperature in the tank just before resuming the swim test was 27.7 ± 0.6 °C, ranging from 26.5 to 28.5 °C, and during the swim test it was 28.3 ± 0.6 °C, ranging from 26.5 to 29.9 °C. Calculation of critical swimming performance and surface area. Absolute and relative critical swimming speed (Ucrit) was used as a measure of swimming performance and calculated according to Brett3: where U−1 is the highest velocity maintained for the prescribed period in cm s−1 , U is velocity increment in cm s−1 , t is time to fatigue at final velocity level in minutes, t is the time each velocity level is maintained at (= 30 min) and SLtest is standard length of fish at swim testing in cm. Figures were produced using Minitab software23. Surface area at swim testing (SAtest) of Nile tilapia is similar to the area of an ellipse and was calculated as: Grow-out in the non-aerated pond. Swim tested fish were stocked in a non-aerated pond for grow-out. The pond size was 500 m2 and the stocking density was 3 fish per m2 . During the grow-out period, DO was above 5 mg l−1 except from 9:00 p.m. to 9:00 a.m. when DO would drop below 3 mg l−1 . Fish were weighed and photographed before stocking into the non-aerated pond. The mean weight at cultivation start (Wstart) was 10.8 g and the coefficient of variation (CV) was 23.7. The fish were fed commercial feed at a rate of 3 to 5% of their body weight depending on their sizes, with the percentage of feed decreasing with size. The fish were harvested after 145 or 146 days of grow-out. Each fish was weighed at harvest. At harvest the sex of a random half of the fish (763 fish) were determined. Daily growth coefficient (DGC)24,25 was computed as: where HW is harvest weight and Wstart is stocking weight. Statistical analysis. Phenotypic and genetic parameters were estimated using ASReml version 4.126. The following animal model was used: where y is a vector of either absolute Ucrit, or relative Ucrit in the univariate model, b is the vector of fixed effects, that is test day and sex fitted as class variable for relative Ucrit while for absolute Ucrit three different models were fitted with: (1) test day and sex fitted as class variables, (2) test day and sex as class variables and Wtest as a covariate and (3) test day and sex as class variables and SLtest as a covariate, sex was not significant in all the three models; therefore, sex was removed from the models; a is a vector of additive genetic effects, c is a vector of environmental effects common to full sibs (‘hapa effect’), and e a vector of residual effects. The X, Z1 and Z2 are design matrices assigning phenotypic values to the levels of fixed effect, additive genetic and common environmental effects, respectively. The effect of sex was also not significant when subset of the data with only 763 sexed fish was analysed. Bivariate models were used to estimate the phenotypic and genetic correlations between absolute Ucrit, and traits such as Wtest, SLtest, Htest, SAtest, HW and DGC. In the bivariate models test day and sex were fitted as a class variable for absolute Ucrit, age at harvest was fitted as a covariate for HW and sex was fitted as a class vari- able for DGC. Common environmental effect was fitted as a random variable to all the traits except for DGC in the bivariate model absolute Ucrit and DGC. The bivariate model with absolute Ucrit and DGC did not converge when a common environmental effect was fitted as a random effect on both traits. The additive genetic effects were normally distributed as N = ([ 0 0 ] ,A⊗ [ σ 2a,1 ra,12σa,1σa,2 ra,21σa,2σa,1 σ 2 a,2 ]) , where A is the numerator genetic (1)Absolute Ucrit = U−1 + ( t t ) U (2)Relative Ucrit = ( U−1 + ( t �t ) �U ) /SLtest (3)SA = 1 4 π ∗ SLtest ∗Htest (4)DGC = [ 3 √ HW − 3 √ Wstart time in days ] × 100 (5)y = Xb+ Z1a + Z2c + e 4Vol:.(1234567890) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ relationships matrix and σ2a,1 ( σ2a,2 ) being the additive genetic variance of trait 1(2) and ra,12(21) being the genetic correlation between trait 1 and 2. The pedigree depth was 18 generations, i.e. from the current generation G18 all the way back to the first generation of GIFT in WorldFish, Malaysia. The common environmental effects were normally distributed as N = ([ 0 0 ] , I⊗ [ σ 2c,1 rc12σc,1σc,2 rc21σc,2σc,1 σ 2 c,2 ]) , where I being an identity matrix and σ2c,1 ( σ2c,2 ) being the common environmental variance of trait 1(2) and rc,12(21) being the common environmental correlation between trait 1 and 2. The residual effects were normally distributed as N = ([ 0 0 ] , I⊗ [ σ 2 e,1 re, 12σ 2 e,1σ 2 e,2 re,21σ 2 e,2σ 2 e, 1 σ 2 e,2 ]) , where σ2e,1 ( σ2e,2 ) being the residual variance of trait 1(2) and re,12(21) being the residual correlation between trait 1 and 2. Heritability ( h2 ) and the ratio of common environmental variance (c2 ) to phenotypic variance ( σ 2p ) of each trait was computed as h2 = σ2a/σ2p and c2 = σ2c/σ2p , respectively. The significance of the random effects were tested using loglikelihood ratio test with one degree of freedom27. To test whether the genetic correlation is larger than zero, a model without constraining the covariance was tested against a model where the covariance was constrained to zero. The full model, i.e. a model with both common environmental effects and additive genetic effects as random effects, was tested against a reduced model, i.e. a model with either only common environmental effect or additive genetic effects as a random effect. The common environmental variances were not significantly different from zero (P > 0.05) except for relative Ucrit (P = 0.006). The most likely reason that the common environmental effect was not significant in most cases was because of the almost complete confounding of sire genetic, dam genetic and common environmental effects in the experiment. This reflected the fact that 24 of the males were mated to one or two females resulting in 12 families with no half sib families and 12 families with only one half sib family (40% of the total families), making genetic and common environmental effects difficult to disentangle. Although, common environmental effects were not significant for most traits, common environmental effects explained a substantial part of the phenotypic variance and were kept in the model, to prevent overestimation of the additive genetic variance. The loglikelihood for the bivariate model with absolute Ucrit and DGC did not converge when common environmental effect was fitted as a random effect on both traits. Therefore, the common environmental effect was fitted as a random effect only on absolute Ucrit in the bivariate model with absolute Ucrit and DGC. Results Biometric data. In total 1500 fish were swim tested and stocked in the non-aerated pond. Out of the swim tested 1500 fish, the swimming performance data of seven fish were missing and resulted in 1493 Ucrit records. The descriptive statistics for age at swim testing (Agetest), Wtest, SLtest, Htest, SAtest, weight at cultivation start, and HW and DGC later in life are presented in Table 1. Out of the stocked 1500 fish, ultimately 1199 were har- vested which is equivalent to 79.9% survival. Swimming performance. The mean absolute Ucrit and relative Ucrit were 69.1 ± 5.5 cm s−1 and 9.7 ± 0.9 SL s−1, respectively (Table 1). Absolute Ucrit and relative Ucrit values showed normal distributions (Fig. 2). There was substantial variation in swimming performance between family means (Fig. 3), indicating existence of genetic variation. Genetic parameters. Variances, heritability and the ratio of common environmental variance to the phe- notypic variance (c2) effect for absolute and relative Ucrit are presented in Table 2. The heritability for absolute Ucrit was 0.48 ± 0.17 when Wtest or SLtest was not fitted in the model as a covariate. The heritability for absolute Ucrit was 0.41 ± 0.16 when SLtest was fitted as a covariate and 0.44 ± 0.16 when Wtest was fitted as a covariate. The Table 1. Number of fish (N), mean, standard deviation (SD), coefficient of variation (CV) and minimum and maximum values for critical swimming speed (Ucrit), absolute and relative, age at swim testing (Agetest), body weight at swim testing (Wtest), standard length at swim testing (SLtest), and body height at swim testing (Htest), surface area at swim testing (SAtest), weight at cultivation start(Wstart), harvest weight (HW) and daily growth coefficient (DGC). N Mean SD CV Min Max Absolute Ucrit (cm s−1) 1493 69.1 5.5 7.9 50.6 83.8 Relative Ucrit (SL s−1) 1493 9.7 0.9 9.8 6.9 13.3 Agetest (days) 1500 86.8 12.1 13.9 65 139 Wtest (g) 1500 10.8 2.6 23.7 4.8 20.1 SLtest (cm) 1500 7.2 0.6 8.0 5.3 8.9 Htest (cm) 1500 2.7 0.2 9.0 2.1 3.5 SAtest (cm2) 1500 15.3 2.5 16.4 9.2 24.0 Wstart (g) 1199 27.4 15.1 47.8 7.3 94.2 HW (g) 1199 417.7 88.1 21.1 153.4 778.9 DGC 1199 3.1 0.3 10.8 1.7 3.8 5Vol.:(0123456789) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ Figure 2. Distribution of absolute Ucrit (cm s−1) and relative Ucrit (SL s−1). Figure 3. Histogram of family average absolute Ucrit (cm s−1) for each of the 60 families. Table 2. Additive genetic variance ( σ 2a ), common environmental variance ( σ 2c ), phenotypic variance ( σ 2p ), heritability and common environmental effect ( c2 ) of absolute and relative critical swimming speed (Ucrit). *Absolute Ucrit without body weight or standard length at swim testing in the model. **Absolute Ucrit when standard length at swim testing was included in the model as covariate. ***Absolute Ucrit when body weight at swim testing was included in the model as covariate. σ 2 a σ 2 c σ 2 p Heritability c2 Absolute Ucrit* 8.90 0.43 18.45 0.48 ± 0.17 0.02 ± 0.05 Absolute Ucrit** 6.71 0.59 16.20 0.41 ± 0.16 0.04 ± 0.05 Absolute Ucrit*** 6.79 0.55 16.09 0.44 ± 0.16 0.03 ± 0.05 Relative Ucrit 0.08 0.07 0.55 0.15 ± 0.13 0.13 ± 0.06 6Vol:.(1234567890) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ heritability for relative Ucrit (0.15 ± 0.13) was low. The common environmental effect explained a small propor- tion of the phenotypic variance (0.02 to 0.04) for absolute Ucrit, while the contribution was substantial for the phenotypic variance of relative Ucrit (0.13). The analyses with absolute Ucrit with Wtest or SLtest in the model as a covariate showed that Ucrit contained considerable heritable variation even when corrected for body size. The additive genetic variance contributed a significant proportion to the phenotypic variance of absolute Ucrit (P = 0.000) and absolute Ucrit when either Wtest (P < 0.001) or SLtest (P = 0.001) was fitted in the model as covariate, while the contribution was not significant for relative Ucrit (P = 0.175). The contribution of common environmental effect to the phenotypic variance of absolute Ucrit (P = 0.584) and for absolute Ucrit when either Wtest (P = 0.384) or SLtest were fitted as covariates (P = 1.000) were not significant, while the contribution to the phenotypic variance of relative Ucrit (P = 0.007) was significant. The estimated genetic correlations (rg) and phenotypic correlations (rp) between absolute Ucrit and Wtest, SLtest, Htest, SAtest, HW and DGC are presented in Table 3. The genetic correlations were significant (P < 0.05) except for the genetic correlation between Ucrit early in life and HW later in life (P = 0.507) based on likelihood ratio test27. The estimated rg and rp correlations between absolute Ucrit and Wtest were 0.78 and 0.44, respectively. The less than one rg between absolute Ucrit and Wtest indicates the presence of genetic variance in absolute Ucrit that is not explained by Wtest. Genetic and phenotypic correlations with the other size measurements SLtest, Htest and SAtest were very similar. Fish with higher absolute Ucrit had lower HW and DGC after grow-out in a non-aerated pond. The estimated rg and rp between absolute Ucrit and HW were − 0.21 and − 0.04, respectively and the estimated rg and rp between absolute Ucrit and DGC were − 0.63 and − 0.24, respectively. The negative genetic correlations between Ucrit and HW and between absolute Ucrit and DGC indicate that fish with higher absolute Ucrit perform less in terms of HW and DGC compared to fish with lower absolute Ucrit. Discussion Our objectives were to estimate variance components for swimming performance in Nile tilapia, assessed as critical swimming speed (Ucrit) early in life, and to estimate the genetic correlation between Ucrit and body size early in life, and harvest weight (HW) and Daily Growth Coefficient (DGC) later in life after a grow-out period in a non-aerated pond. For the first time, we show with a large-scale experiment that swimming performance is heritable in Nile tilapia, and that the genetic correlation with harvest weight is strongly negative, even when corrected for body size at testing. Heritabilities, the genetic correlations, methodology and the practical applica- tion of a swimming performance test in breeding programs are discussed. This study shows the existence of heritable variation in critical swimming performance with a moderate heritability of 0.41–0.48. Our heritability estimate for Ucrit early in life is in the same range as reported previ- ously for other species and for similar traits (for summary see Table 4). Of the four studies that estimated genetic parameters for swimming performance in fish, only the study that assessed the burst swimming performance trait is not comparable with Ucrit in our study14. Our heritability estimate for relative Ucrit (0.15) was not significantly different from zero, which is different from the heritability of 0.55 for relative maximum swimming speed in European sea bass Dicentrarchus labrax17. The difference in heritability of relative Ucrit might be due to a species specific difference, particularly reflecting the high or long body shape of tilapia and seabass, respectively. Species specific differences also exist in the relation between Ucrit and body size. Absolute Ucrit was genetically strongly correlated with body weight at swim testing (0.78). This is higher than the estimated genetic correla- tion between swimming stamina and body weight in Atlantic salmon (0.23)15. The genetic correlation between absolute Ucrit and standard length (0.83) was also different from the estimated rg between swimming stamina and fork length in Atlantic salmon (− 0.14)15. To the best of our knowledge there are no studies on rg between absolute Ucrit and traits such as Htest, SAtest, HW and DGC with which to compare our results. In our study, the rg estimates between absolute Ucrit and Htest, between absolute Ucrit and SAtest were 0.72 and 0.83, respectively. These strong genetic correlations between absolute Ucrit and SLtest and Wtest early in life show that larger fish swim faster in absolute terms. The estimated rg values between absolute Ucrit and HW and absolute Ucrit and DGC were − 0.21 and − 0.63, respectively, meaning that fish with high Ucrit at testing had lower growth rate (DGC) and harvest weight (HW) later in life. These negative genetic correlations do not support our hypothesis that Nile tilapia with higher Ucrit, reflecting higher oxygen uptake efficiency, are those that perform better in terms of weight increase in non- aerated ponds where hypoxia is frequent. Instead, the negative rg shows that fish with higher Ucrit early in life show Table 3. Genetic and phenotypic correlations between absolute critical swimming speed (Ucrit) and body weight at swim testing (Wtest), standard length at swim testing (SLtest) height at swim testing (Htest), surface area at swim testing (SAtest), harvest weight (HW) and daily growth coefficient (DGC). Ucrit was estimated in a bivariate model without Wtest or SLtest as covariate. rg rp Wtest 0.78 ± 0.18 0.44 ± 0.05 SLtest 0.83 ± 0.19 0.43 ± 0.05 Htest 0.72 ± 0.22 0.37 ± 0.05 SAtest 0.83 ± 0.18 0.42 ± 0.05 HW − 0.21 ± 0.29 − 0.04 ± 0.06 DGC − 0.63 ± 0.15 − 0.24 ± 0.07 7Vol.:(0123456789) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ less body weight increase later in life. These data do not provide insight on fish body shape and composition at slaughter size. For example, it may be that fish with higher Ucrit are the leaner fish later as compared to fish with lower Ucrit. Fish with lower Ucrit may be heavier but not necessarily have more fillet mass. Results of a Ucrit test in Gilthead seabream (Sparus aurata), also a high bodied fish, showed that the (residual) Ucrit was negatively correlated with fillet mass suggesting that fast swimmers build lower fillet mass later in life5. A plausible expla- nation for our results may be the existence of a juvenile trade-off between swimming and growth performance where fish with high Ucrit early in life show slower growth later. Young juveniles may choose to either swim fast or grow fast, which may represent, for instance, two anti-predator strategies: to be able to escape predators or to become too large to be eaten rapidly. Studies have shown that a trade-off between growth rate and locomotor performance can exist28, for instance during accelerated growth29 which can negatively influence muscle cel- lularity and development30,31. Indeed, fast-growing growth hormone (GH) transgenic carp32 had lower critical swimming performance than non-transgenic controls. Fast-growing GH transgenic salmon had similar critical swimming speeds than non-transgenic controls but were also able to consume considerable more oxygen33 and may thus have compensated for lower critical swimming performance. In our study, 1,493 fish were used to estimate genetic parameters. The mating ratio used to produce the experimental fish was 1 male to 1 – 4 females, which gave full sib and half sib families. The previous studies that estimated genetic parameters used a much lower number of fish (range 96–129) as compared to our study and Table 4. Summary of the studies that estimated heritability for different swimming performance traits, genetic phenotypic correlation between swimming performance and body weight, and between swimming performance and body length. Trait Species Comments Heritability Genetic (rg) and phenotypic (rp) correlations References Critical Swimming speed Guppy (Poecilia reticulata) Measured by increasing the water velocity every 3 min until the fish fatigued 16 full sib families were used (96 fish in total) 0.24 ± 0.19 Not given 16 Swimming stamina (similar trait with critical swimming speed) Atlantic salmon (Salmo salar) Measured as the total time the fish swam until fatigue by increasing water velocity incrementally every 4 min 11 full sib families were used (129 fish in total) 0.24 ± 0.16 rg = 0.23 and r P = 0.85 (between stamina and body weight) rg = − 0.14 and r P = 0.18 (between stamina and body length) 15 Absolute burst swimming (cm/s) (not comparable with critical swim- ming speed) Threespine stickleback (Gasterosteus aculeatus) Measured as distance swam in 160 ms using video recording, 2 months old 193 fish from 25 full sib families were used 0.41* Not given 14 Relative burst swimming (body length/s) (not comparable with criti- cal swimming speed) Threespine stickleback Measured as distance swam in 160 ms using video recording 2 months old 193 fish from 25 full sib families were used 0.37 Not given 14 Absolute burst swimming (cm/s) (not comparable with critical swim- ming speed) Threespine stickleback Measured as distance swam in 160 ms using video recording 3.6 months old 181 fish from 25 full sib families were used 0.02 Not given 14 Relative burst swimming (body length/s) (not comparable with criti- cal swimming speed) Threespine stickleback Measured as distance swam in 160 ms using video recording 3.6 months old 181 fish from 25 full sib families were used 0.00 Not given 17 Relative maximum sustained speed (similar trait with relative critical swimming speed) European sea bass (Dicentrarchus labrax) Measured as the last fully accom- plished water velocity 547 fish from 366 full sib families, paternal and maternal half sib families were used 0.55 rg = − 0.64 and r P = − 0.56 between relative maximum sustained speed and body weight 17 Absolute critical swimming speed Nile tilapia (Oreochromis niloticus) Explained in Sect. 2.3 of this paper 1493 fish, full sib and half sib families Absolute Ucrit without including body weight/standard length in the model as covariate 0.48 rg = 0.87 and r P = 0.44 between absolute Ucrit and body weight This study Absolute critical swimming speed Nile tilapia 1493 fish, full sib and half sib families Absolute Ucrit when body weight was included in the model as covariate 0.42 This study Absolute critical swimming speed Nile tilapia 1493 fish, full sib and half sib families Absolute Ucrit when standard length was included in the model as covariate 0.41 This study Relative critical swimming speed Nile Tilapia 1493 fish, full sib and half sib families 0.15 This study 8Vol:.(1234567890) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ estimated broad sense heritability using full sib families (Table 4)15,16. The much larger sample size gave a much higher precision of estimates of narrow-sense heritability. Furthermore, broad-sense heritability estimates are biased estimates of narrow-sense heritabilities, because broad-sense heritabilities contain non-additive genetic variation due to dominance and epistasis that is not heritable from parent to offspring and may contain com- mon environmental effects, because in such full-sib designs estimation of common environmental effects is not feasible34. Narrow-sense heritability, however, is the ratio of additive genetic variance to phenotypic variance13 and therefore a better indication of the proportion of genetic variation that is transmitted to the next generation. In our study, we used half sib families that enabled us to estimate a narrow sense heritability. Similarly, Vandeputte et al.17 estimated narrow sense heritability using half sib families based on 547 fish. The main difference in the swimming performance trait between our study and Vandeputte et al.17 was that these authors did not include the last water velocity level that the fish did not fully complete. Besides the species difference mentioned earlier, also the number of fish and the way the swimming performance was calculated could provide additional explanation for the difference in the parameter estimates between our study and Vandeputte et al.17. Critical swimming speed can be calculated in four different ways: as absolute Ucrit, with or without Wtest or SLtest as covariate in the model, as relative Ucrit, or as residual Ucrit which is the difference in Ucrit of an individual fish with the predicted value on basis of its length5. Analysing absolute Ucrit without a covariate for either Wtest or SLtest, has the highest additive genetic variance, but part of that genetic variance is due to genes affecting body size. The use of fish with similar body weight at similar SLtest is practically difficult as the variation is consider- able; in our experiment the Wtest was from 4.8 to 20.1 g for fish from 5.4 to 10 cm SL. Therefore, it is important to account for Wtest or SLtest in the analysis to be able to estimate heritable variation in Ucrit independent of body size. Relative Ucrit is a ratio of Ucrit to SLtest for which the estimated heritability was not significantly different from zero in our study. Relative Ucrit is a ratio trait and therefore the genetic variance becomes a complex function of absolute Ucrit and SLtest. Ratio traits are generally not recommended in animal breeding35. For instance, the herit- ability of a ratio trait cannot be used to predict the genetic change for the ratio trait36. Therefore, we recommend using the absolute Ucrit and to fit either Wtest or SLtest as a covariate in a model when estimating heritability. Such an analysis shows the existence of heritable variation in Ucrit beyond body size. The less than unity genetic correlation between absolute Ucrit and Wtest indicates the presence of genetic variation in Ucrit, independently of Wtest. A genetic correlation of unity between two traits means that the two traits are controlled by the same genes while a genetic correlation of less than unity indicates that there are additional genes that are not common for the two traits and only control one of the two traits. The negative rg between absolute Ucrit and HW, and between Ucrit and DGC, clearly indicates that selection for high harvest weight will favour faster growing animals with lower Ucrit. Whether this is desirable needs to be determined. One can speculate that under conditions of hypoxia, as frequently encountered in non-aerated ponds or ponds with algal blooms, smaller, more active fish will have a higher chance of survival. In optimal management con- ditions, however, growth rate can be further increased by including Ucrit at testing in the breeding goal, next to harvest weight. Fish with higher Ucrit may also be more resilient: swimming exercise improves physiological fitness; cardiovascular and respiratory performance, and increases mitochondrial densities and muscle tissue capillarization37. Also the immune system capacity appears to be linked to swimming performance as Castro et al.38 found 21 virus-responsive genes with significantly higher transcript abundance in phenotypically poor swimmers as compared to good swimmers in Atlantic salmon. In conclusion, including absolute Ucrit in a breeding goal in addition to HW and DGC could be beneficial if the aim is to select for fitter fish, especially in environments where oxygen is limiting. Absolute Ucrit can be measured at an early stage on the selection candidates themselves, high throughput and non-invasively although size of the tested fish may be restricted due to difficulty in reaching sufficiently high flow speed. However, selection on Ucrit with 10% selection intensity for the highest values of Ucrit could lead to a 19% reduction in mean harvest weight of the offspring, compared to direct selection on harvest weight. In practice we recommend a two-stage selection scheme, where selection in the first stage is on retaining 90% of the fittest fish in terms of Ucrit, followed by a second stage selection on harvest weight. This study showed for the first time the existence of significant additive genetic variance for critical swimming speed in Nile tilapia. Favourable rg between Ucrit and traits such as Wtest, SLtest, Htest and SAtest early in life were found. The main finding demonstrated a negative rg between Ucrit and HW later in life, and between Ucrit and DGC later in life. Including Ucrit in the breeding goal may help to improve resilience of Nile tilapia. Data availability Data are available from the corresponding author upon reasonable request. Received: 15 December 2020; Accepted: 6 May 2021 References 1. Stickney, R. R. Aquaculture: An Introductory Text (CABI, Wallingford, 2017). 2. Mengistu, S. B., Mulder, H. A., Benzie, J. A. H. & Komen, H. A systematic literature review of the major factors causing yield gap by affecting growth, feed conversion ratio and survival in Nile tilapia (Oreochromis niloticus). Rev. Aquacult. 12, 524–541. https:// doi. org/ 10. 1111/ raq. 12331 (2020). 3. Brett, J. R. The respiratory metabolism and swimming performance of young sockeye salmon. J. Fish. Res. Board Can. 21, 1183– 1226. https:// doi. org/ 10. 1139/ f64- 103 (1964). 4. Plaut, I. Effects of fin size on swimming performance, swimming behaviour and routine activity of zebrafish Danio rerio. J. Exp. Biol. 203, 813–820 (2000). 9Vol.:(0123456789) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ 5. Palstra, A. P., Kals, J., Böhm, T., Bastiaansen, J. & Komen, H. Swimming performance and oxygen consumption as non-lethal indicators of production traits in Atlantic salmon and Gilthead seabream. Front. Physiol. https:// doi. org/ 10. 3389/ fphys. 2020. 00759 (2020). 6. Videler, J. J. Fish Swimming (Springer, New York, 1993). 7. Jones, D. R. & Randall, D. J. Respiration and circulation during exercise in fish. In Fish Physiology Vol. 7 (eds Hoar, W. S. & Randall, D. J.) 425–501 (Academic Press, New York, 1978). 8. McKenzie, D. J. et al. Effects of growth hormone transgenesis on metabolic rate, exercise performance and hypoxia tolerance in tilapia hybrids. J. Fish Biol. 63, 398–409. https:// doi. org/ 10. 1046/j. 1095- 8649. 2003. 00162.x (2003). 9. Kim, J. et al. Genetic selection of athletic success in sport-hunting dogs. Proc. Nat. Acad. Sci. USA 115, E7212–E7221 (2018). 10. Hellsten, T. et al. Review of genetic parameters estimated at stallion and young horse performance tests and their correlations with later results in dressage and show-jumping competition. Livest. Sci. 103, 1–12 (2006). 11. Hill, E. W., Ducro, B. J., van Weeren, R., Barneveld, A. & Back, W. Genetic contributions to exercise and athletic performance. In Equine Locomotion 2nd edn (eds Back, W. & Clayton, H. M.) 245–266 (Elsevier, New York, 2013). 12. Issurin, V. B. Evidence-based prerequisites and precursors of athletic talent: A review. Sports Med. 47, 1993–2010. https:// doi. org/ 10. 1007/ s40279- 017- 0740-0 (2017). 13. Falconer, D. S. & Mackay, T. F. C. Introduction to Quantitative Genetics 4th edn, 464 (Longman, Harlow, 1996). 14. Garenc, C., Silversides, F. G. & Guderley, H. Burst swimming and its enzymatic correlates in the threespine stickleback (Gasterosteus aculeatus): Full-sib heritabilities. Can. J. Zool. 76, 680–688. https:// doi. org/ 10. 1139/ z97- 236 (1998). 15. Hurley, S. M. & Schom, C. B. Genetic control of swimming stamina in Atlantic salmon (Salmo salar). Can. J. Genet. Cytol. 26, 57–61. https:// doi. org/ 10. 1139/ g84- 010 (1984). 16. Nicoletto, P. F. Offspring quality and female choice in the guppy, Poecilia reticulata. Anim. Behav. 49, 377–387. https:// doi. org/ 10. 1006/ anbe. 1995. 0050 (1995). 17. Vandeputte, M. et al. Quantitative genetic variation for post-stress cortisol and swimming performance in growth-selected and control populations of European sea bass (Dicentrarchus labrax). Aquaculture 455, 1–7. https:// doi. org/ 10. 1016/j. aquac ulture. 2016. 01. 003 (2016). 18. Mengistu, S. B. et al. Genotype by environment interaction between aerated and non-aerated ponds and the impact of aeration on genetic parameters in Nile tilapia (Oreochromis niloticus). Aquaculture 529, 735704. https:// doi. org/ 10. 1016/j. aquac ulture. 2020. 735704 (2020). 19. Rohlf, F. J. tpsUtil, Version 1.74. (Department of Ecology and Evolution, State University of New York at Stony Brook, 2017). 20. Rohlf, F. J. tpsDig2, Version 2.30. (Department of Ecology and Evolution, State University of New York at Stony Brook, 2017). 21. Adams, D. C., Collyer, M. L., Kaliontzopoulou, A. Geomorph: Software For Geometric Morphometric Analyses. R package Version 3.0.6. (2018). 22. Palstra, A. P. Physiological Testing of Fish in WorldFish Facilities in Malaysia. Design of Swim Flume and Standard Operating Pro- cedure for High Throughput Testing of Juvenile Tilapia. Confidential Report 516 pp. 16 (Wageningen Livestock Research, 2016). 23. Minitab Statistical Software Version 17. (2010). https:// www. minit ab. com 24. Iwama, G. K. & Tautz, A. F. A simple growth model for salmonids in hatcheries. Can. J. Fish. Aquat. Sci. 38, 649–656. https:// doi. org/ 10. 1139/ f81- 087 (1981). 25. Trọng, T. Q., Mulder, H. A., van Arendonk, J. A. M. & Komen, H. Heritability and genotype by environment interaction estimates for harvest weight, growth rate, and shape of Nile tilapia (Oreochromis niloticus) grown in river cage and VAC in Vietnam. Aqua- culture 384–387, 119–127. https:// doi. org/ 10. 1016/j. aquac ulture. 2012. 12. 022 (2013). 26. ASReml User Guide Release 4.1 Functional Specification (2015). 27. Lynch, M. & Walsh, B. Genetics and Analysis of Quantitative Traits 980 (Sunderland, Sinauer, 1997). 28. Billerbeck, J. M., Lankford, T. E. & Conover, D. O. Evolution of intrinsic growth and energy acquisition rates. I. Tradeoffs with swimming performance in Menidia menidia. Evolution 55, 1863–1872 (2001). 29. Lee, W.-S., Monaghan, P. & Metcalfe, N. B. The trade-off between growth rate and locomotor performance varies with perceived time until breeding. J. Exp. Biol. 213, 3289–3298 (2010). 30. Galloway, T. F., Kjorsvik, E. & Kryvi, H. Muscle growth and development in Atlantic cod larvae (Gadus morhua L.) related to dif- ferent somatic growth rates. J. Exp. Biol. 202, 2111–2120 (1999). 31. Johnston, I. A. Muscle metabolism and growth in Antarctic fishes (suborder Notothenioidei): Evolution in a cold environment. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 136, 701–713 (2003). 32. Li, D. et al. Rapid growth cost in “all-fish” growth hormone gene transgenic carp: Reduced critical swimming speed. Chin. Sci. Bull. 52, 1501–1506 (2007). 33. Stevens, E. D., Sutterlin, A. & Cook, T. Respiratory metabolism and swimming performance in growth hormone transgenic Atlantic salmon. Can. J. Fish. Aquat. Sci. 55, 2028–2035 (1998). 34. Lozano-Jaramillo, M., Komen, H., Wientjes, Y. C. J., Mulder, H. A. & Bastiaansen, J. W. M. Optimizing design to estimate genetic correlations between environments with common environmental effects. J. Anim. Sci. https:// doi. org/ 10. 1093/ jas/ skaa0 34 (2020). 35. Zetouni, L., Henryon, M., Kargo, M. & Lassen, J. Direct multitrait selection realizes the highest genetic response for ratio traits. J. Anim. Sci. 95, 1921–1925. https:// doi. org/ 10. 2527/ jas. 2016. 1324 (2017). 36. Gunsett, F. C. Merit of utilizing the heritability of a ratio to predict the genetic change of a ratio. J. Anim. Sci. 65, 936–942. https:// doi. org/ 10. 2527/ jas19 87. 65493 6x (1987). 37. Palstra, A. P. & Planas, J. V. Fish under exercise. Fish Physiol. Biochem. 37, 259–272 (2011). 38. Castro, V. et al. Disease resistance is related to inherent swimming performance in Atlantic salmon. BMC Physiol. 13, 1–12. https:// doi. org/ 10. 1186/ 1472- 6793- 13-1 (2013). Acknowledgements This publication was made possible through support provided principally by the Koepon Foundation, by the International Fund for Agricultural Development (IFAD), the CGIAR Research Program on Fish Agri-Food Systems (FISH) led by WorldFish. The program is supported by contributors to the CGIAR Trust Fund. We thank the WorldFish breeding team, and the Department of Fisheries Malaysia at Jitra for their technical sup- port throughout the project. Author contributions S.M.: Writing—Original draft preparation, Designing the experiment, Investigation, Data curation, Formal analysis. A.P.: Writing—Review & Editing, Conceptualization, Methodology, Designing the experiment. H.M.: Writing—Review & Editing, Conceptualization, Designing the experiment, Software, Supervision. J.B.: Writing— Review & Editing, Funding acquisition, Resources. T.Q.T.: Writing—Review & Editing, Investigation, Resources. C.R.: Writing -Review & Editing, analysis. H.K.: Writing—Review & Editing, Conceptualization, Designing the experiment, Funding acquisition, Supervision, Project administration. 10 Vol:.(1234567890) Scientific Reports | (2021) 11:11018 | https://doi.org/10.1038/s41598-021-90418-w www.nature.com/scientificreports/ 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- 021- 90418-w. Correspondence and requests for materials should be addressed to A.P.P. 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. 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:// creat iveco mmons. org/ licen ses/ by/4. 0/. © The Author(s) 2021