J. Anim Sci.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kause, A.
Right arrow Articles by Ruohonen, K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kause, A.
Right arrow Articles by Ruohonen, K.
J. Anim. Sci. 2006. 84:807-817
© 2006 American Society of Animal Science


ANIMAL GENETICS

Feed efficiency of rainbow trout can be improved through selection: Different genetic potential on alternative diets1

A. Kause*,2, D. Tobin{dagger}, D. F. Houlihan{dagger}, S. A. M. Martin{dagger}, E. A. Mäntysaari*, O. Ritola{ddagger} and K. Ruohonen§

* MTT Agrifood Research Finland, Animal Production Research, Animal Breeding, FIN-31600 Jokioinen, Finland; and {dagger} School of Biological Sciences, Department of Zoology, University of Aberdeen, Aberdeen, AB24 2TZ, UK; and {ddagger} Finnish Game and Fisheries Research Institute, Tervo Fisheries Research and Aquaculture, FIN-72210 Tervo, Finland; and and § Finnish Game and Fisheries Research Institute, Turku Game and Fisheries Research, FIN-20520, Finland


    Abstract
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
To assess the genetic potential for selection of increased feed efficiency in rainbow trout (Oncorhynchus mykiss), we estimated the heritabilities and correlations for BW, daily weight gain (DG), and daily feed intake (DFI). Body weight was recorded 5 times, and DG and DFI 3 times during a feeding trial lasting 22 mo. To test the hypothesis that phenotypic and genetic parameters were influenced by a nutritional environment, fish were fed either a modern normal protein diet (NP, 40 to 45% protein and 30 to 33% lipid) or an alternative high protein diet (HP, 50 to 56% protein, 20 to 24% lipid) in a split-family design. Results showed that there were no large differences in heritabilities between the diets. Average heritability for DFI over both diets and different fish ages was low (average h2 = 0.10), indicating that modest genetic changes in response to selection can be obtained. Average heritabilities for BW and DG over both diets and different fish ages were 0.28 and 0.33, respectively. The NP diet enabled fish to express a wide range of BW, as shown by the increased coefficients of phenotypic variation for BW. Fish fed the HP diet showed increased phenotypic variation for DFI in > 750-g fish. On the NP diet, genetic correlations of DFI with DG and BW were very strong for 750- to 2,000-g fish. In contrast, on the HP diet, the respective correlations were moderate to low, revealing more genetic potential to change growth and feed intake simultaneously in opposite directions. An analysis of the predicted selection responses showed that selection solely for high DG improved feed efficiency as a correlated genetic response. Simultaneous selection for high DG and reduced DFI, in turn, may increase genetic gain in feed efficiency by a factor of 1.2 compared with selection solely for DG. However, variation for growth and feed intake and the relationships between these traits were different in different nutritional environments, leading to divergent genetic responses on the alternative diets.

Key Words: feed efficiency • feed intake • genotype-by-environment interaction • heritability • quantitative genetics • rainbow trout


    INTRODUCTION
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Studies on the quantitative genetics of feed efficiency in fish are lacking. Moreover, it is possible that current commercial fish feeds impose unfavorable genetic constraints that hamper selection efforts to improve feed efficiency. Here, we estimated heritabilities and genetic correlations for daily feed intake (DFI), daily weight gain (DG), and BW in rainbow trout fed 2 diets; a normal protein (NP) diet representing a modern feed, and an experimental high protein (HP) diet.

It was hypothesized that the selection potential for feed efficiency should be greater on an HP diet. First, we hypothesized that phenotypic and genetic variation in growth should be greater on an HP diet, allowing greater selection responses than on an NP diet. This is based on the idea that on an NP diet, protein growth is restricted due to the low dietary protein level (Kim and Kaushik, 1992Go), and high lipid level of NP diet facilitates high lipid deposition in the body of all individuals. This is assumed to lead to similar wet weight growth among individuals, regardless of their growth composition, thereby reducing variation.

In contrast, on an HP diet, the variation in BW is increased because of the variation in the capacity of fish to utilize dietary protein for growth (protein efficiency; Houlihan et al., 1995Go). Due to the high protein content of an HP diet, the capacity of inefficient fish to deposit digested protein as protein growth is greatly exceeded (Kim and Kaushik, 1992Go), and the low lipid content restricts their lipid growth. This results in poor overall growth of the inefficient fish, and hence, leads to increased variation on HP diet.

Second, we tested whether growth and feed intake were strongly correlated, as has been shown in farm animals (e.g., Clutter and Brascamp, 1998Go). Moreover, we tested whether this tight correlation could be uncoupled on the HP diet, enhancing the selection potential. Third, we assessed the usefulness of different selection strategies to improve feed efficiency.


    MATERIALS AND METHODS
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Population Structure
The rainbow trout (Oncorhynchus mykiss Walbaum) used in this study belonged to the Finnish national breeding program. All fish were housed at the Tervo Fisheries Research and Aquaculture station (a freshwater station) of the Finnish Game and Fisheries Research Institute for the duration of the experiment (June 2001 to November 2004). The breeding scheme was reported previously by Kause et al. (2005)Go.

To estimate genetic parameters, the 2001 generation was exposed to the 2 different diet treatments in a split-family design. The pedigree for every fish was known as far back as 4 generations. A total of 210 full/half-sib families were produced from 89 sires and 109 dams of 3 yr of age in a factorial mating design. Matings were completed during 3 d. Full-sib batches of eye-staged eggs were transferred to indoor 150-L family tanks in June 2001.

In February 2002, after 8 mo of growing in the family tanks, 2,931 fingerlings were removed from the family tanks, and individually tagged with PIT-transponders (Trovan Ltd., Köln, Germany) to enable individual identification. For 45 of the 210 full-sib families, an average of 39.6 individuals per family (range: 37 to 40) were randomly tagged. The large initial family size was used because these families were later sampled not only for BW and feed intake but also for destructive composition traits necessary for a parallel study (our unpublished observations). For the remaining 165 full-sib families, nondestructive sampling was planned, so an average of 7 fish per family (range: 4 to 7) were randomly tagged.

Fish Management
During tagging in February, each family was randomly split into 2 groups to be reared with different experimental diets. In May 2002 (at wk 20), the 2 diet treatments were initiated. Before the initiation of the diet treatments, all fish were fed with commercial rainbow trout dry food (Nutra Starter and Nutra Parr, Rehuraisio Inc., Raisio, Finland). The experimental diets were a modern diet with low protein (40 to 45% of DM) and high lipid content (30 to 33% of DM; NP diet), and an alternative diet with high protein (50 to 56% of DM) and low lipid content (15 to 24% of DM; HP diet; Table 1Go). The diets consisted of fish meal, fish oil, wheat meal, and wheat starch, and were supplemented for minerals and vitamins according to National Research Council (1993)Go recommendations.


View this table:
[in this window]
[in a new window]
 
Table 1. Proximate composition1 of normal protein (NP) and high protein diets (HP), including glass bead-labeled x-ray diets fed at measurement times 1 (bead size, 3 mm), 2 (bead size, 6 mm), and 4 (bead size, 7 mm)
 
The fish assigned to the 2 diets did not differ in weight at tagging (mean ± SD; NP = 62.4 ± 19.9 g, n = 1,355; and HP = 62.3 ± 19.4 g, n = 1,335). Each diet group was equally divided into 4 replicate fiberglass tanks and housed indoors. The fish were kept in 3-m3 tanks until wk 24, and transferred to 20-m3 tanks thereafter. The families were equally distributed among the tanks. Fish density in the tanks was under 20 kg/m3.

Feed was provided 4 h a day by computer-controlled pneumatic feeders (Arvo-Tec Inc., Huntokoski, Finland). To ensure overfeeding and to prevent restriction of the growth potential of the fish, feed consumption was controlled daily by visual observation. Feeding was regarded as overfeeding when uneaten feed remained on the bottom of the tanks after 4 h of feeding. The amount of excess feed was 10 to 20% over the recommendations of the feed company. Water temperature during the experiment was ambient, with seasonal fluctuations.

Recording Traits
The approach we adopted here follows the logic that genetic components of feed efficiency (DG:DFI) can be sufficiently described by the analysis of its constituent traits, DG and DFI (Gunsett, 1984Go; Kennedy et al., 1993Go; Cammack et al., 2005Go). The disadvantage of analyzing feed efficiency per se is that a trait defined as a ratio is often ill behaved statistically, and it is difficult to control which of its component traits is changing when selection is applied on the ratio. Moreover, genetic parameters of feed efficiency are determined by its component traits, and feed efficiency can be improved by selecting for the component traits with appropriate economic weights (Gunsett, 1984Go; Kennedy et al., 1993Go).

Body weight was recorded 5 times: in May 2002 (time 1), October 2002 (time 2), May 2003 (time 3), September 2003 (time 4), and November 2003 (time 5; traits: BW1 to BW5). Daily feed intake was recorded 3 times: at times 1, 2, and 4 (traits: DFI1,2,4). The average family sizes at different times are given in Table 2Go.


View this table:
[in this window]
[in a new window]
 
Table 2. Population structure and average family size of the 210 full/half-sib families used1
 
Recording times 1, 2, and 4 consisted of a 3-wk x-ray session with 3 repeated measurements of BW and DFI. The x-radiography was used to record individual feed consumption, using a portable x-ray unit (Todd Research 80/20, Essex, UK), as described by Talbot and Higgins (1983)Go. During each 3-wk session, each of the 8 tanks was measured once weekly (1 NP and 1 HP tank per day) in the same order each week. To avoid the potential effects of systematic feeding rhythms, the recording order of NP and HP tanks was reversed on successive days.

To initiate a recording session, all fish (x-ray and non-x-ray) were weighed during the first week of each session, and DFI was measured from randomly selected individuals from each family. In the second and the third weeks, the procedure was repeated but only the fish x-rayed in the first week were reweighed and x-rayed again. A 1-wk difference between DFI recordings was considered appropriate because fish appetite is reduced very little, even during the day following x-raying (McCarthy et al., 1993Go).

The x-rays were performed in the same way for all tanks during all sessions. Before x-ray, all fish from a given tank were fed as usual but the diet was labeled with radio-opaque ballotini glass beads (Jencons Scientific Ltd., Leighton Buzzard, UK; Table 1Go). The labeled pellets used at times 1, 2, and 4 consisted of 1, 0.5, and 0.3% beads, respectively, with a diameter of 400 to 600 µm. To avoid increased time and labor needed for counting the beads of large fish, bead density was reduced with increasing fish age. This is because larger fish consume more, and thus the reduced bead density prevents the absolute bead number within a fish from being elevated and unnecessarily large.

A minimum of 2 h after the feeders had stopped providing feed, the fish to be recorded were serially placed into anesthetic solution (buffered MS-222, Argent Laboratories, Redmond, WA), weighed, and x-rayed. Thereafter, the beads were counted from the films, and the weight of feed within a stomach was estimated using a calibration regression equation. The predictive calibration equation was constructed for each session and each diet separately by x-raying different but known amounts of feed (n = minimum of 8 feed samples), and then by regressing the number of counted beads from the films against the known feed weight (R2 = 0.90 to 0.99). Evacuation of the feed before x-ray was avoided by adjusting the duration of the feeding and the length of the recording period, using information from a separate experiment (data not shown).

At times 3 and 5, each fish was anesthetized and recorded for BW once. The fish were recorded for sex and maturity at time 3 by a visual inspection of secondary sexual characters, and at time 5 by the examination of gonads from slaughtered fish. Six sex x maturity classes were identified: males that were mature at 2, 3, and later years; females that were mature at 3 and later years; and fish with unknown sex and maturity age.

Traits Analyzed
Feed intake was measured daily. For this reason, we also calculated growth from BW records as DG. Daily gains at times 1, 2, and 4 for each individual were calculated following the methods of Iwama and Tautz (1981)Go and Cho (1992)Go, and using information on the daily water temperature and on the 2 to 3 repeated BW records obtained for each x-rayed individual during each x-ray session (traits: DG1,2,4). First, a growth rate (the regression coefficient bi) for the ith fish over a 3-wk period was calculated from a regression:


Formula 1[1]

where BWij1/3 is the cubic root of the BW record measured for individual i during day j, a is an intercept of the regression, CTEMPj is a cumulative temperature sum at day j, and eij is the residual. Then, DG for each individual was obtained using the following equation:


Formula 2[2]

where the terms in parentheses are the predicted cubic root BW during 2 successive days (j and j+1) obtained using Eq. 1. In other words, day j is the day when a fish was measured for BW (and x-rayed) and day j+1 is the consecutive day. Using this method, DG (in units of g/d) can be calculated for each individual with 2 or more BW records.

Next, for each individual at each 3-wk x-ray session, the repeated observations for DFI, DG, and BW were each compressed into a single (more accurate) mean. This approach was justified because salmonids have extremely high day-to-day variability in feed intake compared with other domesticated animals. For instance, an individual fish does not feed every day or may feed very little on some days (McCarthy et al., 1992Go; Jobling and Koskela, 1996Go). This is indicated by the low repeatability of DFI (r = 0.09 to 0.32; our unpublished observations). Thus, single-point estimates of DFI are potentially inaccurate, and averaging several repeated observations increases the measurement accuracy considerably (Falconer and Mackay, 1996Go).

Because some of the x-rayed fish lacked 1 of the 3 observations within a session and because of the differences among the average daily trait means, the raw means for each individual would have been unsuitable for the subsequent calculations. Therefore, weighted least squares means for each individual were obtained by accounting for the test tank-wise weekly performance. This was performed for each tank separately by fitting an ANOVA to the longitudinal feed intake data with a model including individual fish (all fish within a given test tank) and test week (wk 1, 2, and 3) as factors, and calculating least squares means for each individual (LSMEANS option, SAS Inst., Inc., Cary, NC). Because test tanks were treated separately, the test tank differences remained in the data and needed to be fitted in the subsequent statistical analyses. These least squares means at times 1, 2, and 4 were used as observations in the subsequent statistical and genetic analyses. The measures of BW3 and BW5 included only 1 record, so least squares means were not required.

Statistical Tests for Diet Differences in Trait Means
To examine the differences between the diets in the trait means, parametric ANOVA were performed for each trait separately (Proc Mixed; SAS Inst., Inc.). The fixed effects included in the model were diet, sex x maturity class, and the interaction of diet with sex x maturity class. The random factors included were test tank nested within diet, family (this consists of both common environmental effects and genetic effects of full-sib families), interaction of test tank with sex x maturity class, interaction of family with diet, and interaction of family with sex x maturity class.

The method of Kenward and Roger (1997)Go was used to calculate the correct degrees of freedom and F-tests for the fixed effects. All DFI values were transformed using square roots to obtain normally distributed residuals. When calculating statistical tests and least squares means for fish with equal size, BW was included as a covariate in the models of DFI and DG. We refer to these traits corrected for BW as DFI%1,2,4 and DG%1,2,4. Performing a repeated analysis of BW, DG, and DFI to compare diet differences across time was not possible because the extensive data size and complex mixed model led to memory capacity problems when using PROC MIXED.

Genetic Analyses
(Co)variance components were estimated using the DMU AI software that was developed for analyzing multivariate mixed models using the restricted maximum likelihood method (Jensen and Madsen, 2000Go).

A trait measured on the 2 diets was regarded as 2 different traits. The model for the diet-specific traits was:


Formula 3[3]

where animi is a random additive genetic effect of an animal (i = 1...number of observations), famtankj is a random family tank effect (j = 1 to 210 tanks), SEX-MATk is a fixed sex and maturity effect (k = 1 to 6), TESTTANKl is a fixed test tank effect (l = 1 to 4 tanks), {varepsilon}ijkl is the residual, and yijkl is an observation from individual i. When analyzing multitrait models with 2 traits measured on different diets, the residual covariance was set to zero.

To obtain genetic parameters for feed intake percentage (DFI adjusted for BW) and daily gain percentage (DG adjusted for BW), BW was included as a covariate into the models of raw DFI and DG. Body weight-adjusted DFI and DG are of interest in assessing the degree to which family differences (heritabilities) and ranking of families across diets (between-diet correlations) are a result of BW differences. They are not, however, needed for the analysis of feed efficiency genetics because the analysis of the components, raw DG and DFI, will suffice.

Heritabilities and genetic correlations between diets were obtained from bivariate models including the same character recorded on both diets. To obtain trait correlations within diets, a multitrait model with BW, DG, and DFI was run separately for each diet and each sampling time.

Using these models, additive genetic ({sigma}2A), common environment due to shared environment of full sibs during incubation and family tank rearing ({sigma}2fam tank), residual ({sigma}2R), and phenotypic variances ({sigma}2P = {sigma}2A + {sigma}2R + {sigma}2fam tank), as well as phenotypic (rP) and genetic correlations between traits (rA) were obtained. Heritability was calculated as h2 = {sigma}2A/{sigma}2P and common environment ratio as c2 = {sigma}2fam tank/{sigma}2P. If the common environment ratio was lower than 0.01, the family tank effect was removed from the model of that trait. To scale the phenotypic variance of traits with different means (Formula 3), coefficient of phenotypic variation was calculated as

Formula 3

Because heritability is a ratio, low heritability can result either from low genetic variation or from high residual variation, or both. Thus, we calculated coefficients of additive genetic variation

Formula 3

and of residual variation

Formula 3

(Houle, 1992Go).


    RESULTS
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Diet Effects
Neither BW nor daily gain differed between diets (P > 0.08; Table 3Go). However, from time 2 onwards, BW was 5.2 to 7.6% greater on the NP than on the HP diet but the differences were not statistically significant. Similarly, at times 2 and 4, fish on the NP diet had 14.8 to 17.2% greater DG, but the differences were not significant (Table 3Go). Growth displays such large variation within diets that the power of the statistical tests with 4 test tanks as replicates was not great enough to reach significance. Body lipid composition, in turn, was significantly greater on the NP diet, confirming that the diets significantly influenced components of growth (our unpublished observations).


View this table:
[in this window]
[in a new window]
 
Table 3. Sample sizes, means, and their SE for traits on 2 diets, and statistical tests for diet effect
 
In contrast to growth, both DFI and DFI percentage were greater, not lower, on the HP diet compared with the NP diet. At times 2 and 4, DFI was 39.3 and 23.3% greater on the HP diet, respectively (P ≤ 0.014; Table 3Go). At time 3, the diet difference was lower (15.6%; P = 0.16). The difference between diets in DFI percentage was even greater than for the raw DFI (all P ≤ 0.054), because of HP fish being smaller but feeding more than NP fish.

Variation on 2 Diets
The hypothesis of greater variation for growth on the alternative HP diet was not supported by the data. In contrast to the hypothesis, coefficients of phenotypic variation for BW were greater on the NP diet compared with the HP diet, and the difference was progressively increased from 3.7% at time 1 to 16.4% at time 5 (Table 4Go). Similarly, all CVP of DG and 2 of 3 CVP of DG percentage were greater on NP compared with HP diet.


View this table:
[in this window]
[in a new window]
 
Table 4. Phenotypic variances ({sigma}2P), coefficients of phenotypic variation (CVP), heritabilities (h2), common environment ratios (c2), and their SE for traits on 2 diets
 
Heritabilities for BW displayed only small differences between diets, the average heritabilities being 0.26 and 0.30 for NP and HP diets, respectively (Table 4Go). Heritabilities of BW were slightly smaller on the NP diet, even though the NP diet displayed greater phenotypic variation. This resulted from the fact that the residual variation was greater on the NP diet (average CVR = 18.8) than on the HP diet (average CVR = 16.7), whereas coefficients of additive genetic variation remained almost unaltered across the diets (average CVA on NP = 11.1 and on HP = 11.0).

For DG, the average heritabilities were 0.37 on the NP and 0.29 on the HP diets, but with overlapping confidence limits (Table 4Go). The average CVA for daily gain was 28.1 on NP and 21.9 on HP, showing that genetic variation was increased on the NP diet. The average CVR was 34.9 on NP and 32.4 on HP diet. For DG percentage, average heritabilities were 0.31 and 0.25 on the NP and HP diets, respectively, and CV showed trends similar to raw DG (data not shown). Consequently, the hypothesis of increased variation on HP diet was not supported by daily gain data either; rather, if any trend was visible it was in contrast to our hypothesis.

For DFI, in 2 out of 3 cases, both CVP and heritabilities were greater on the HP diet (Table 4Go). Average CVP were 41.0 on the NP and 49.2 on the HP diet, showing increased phenotypic variation on the HP diet. Moreover, average heritability for DFI was 0.07 on NP and 0.13 on HP diet, and average CVA was 10.4 on the NP and 17.0 on the HP diet, the values showing greater genetic variation on the HP diet. However, many of the heritabilities for DFI were small and not significantly different from zero. At the phenotypic level, DFI percentage showed increased variation on the HP diet (Table 4Go), whereas CVA did not show great differences between the diets (data not shown).

Correlations Between Growth and Feed Intake
As hypothesized, growth traits and DFI were strongly correlated on the NP diet. On the NP diet, DG and DFI were strongly correlated at all sampling times (Table 5Go). The phenotypic and genetic correlations ranged from 0.51 to 0.74 and 0.86 to 0.96, respectively. Similarly, the phenotypic (0.48 to 0.54) and genetic correlations (0.72 to 0.90) between BW and DFI were all high on the NP diet.


View this table:
[in this window]
[in a new window]
 
Table 5. Phenotypic (above diagonal) and genetic (below diagonal ± SE) correlations between traits1 on 2 diets
 
On the HP diet, growth traits and DFI were uncoupled with successive time measurements (Table 5Go). At time 1 on the HP diet, growth traits and DFI were strongly correlated, as was the case for the NP diet. In contrast, at the later sampling times on the HP diet, the correlations between growth traits and DFI were reduced. At time 4 on the HP diet, the phenotypic correlations of DFI with BW and DG were 0.06 and 0.27, and the genetic correlations were –0.29 and 0.28, respectively (Table 5Go). These correlations indicate better possibilities for simultaneous selection for high DG and against DFI on HP diet.

Genotype-by-Environment Interactions
Genetic correlations between the same measure recorded on both diets were estimated to investigate whether genetic improvement on one diet can be enhanced by selecting on an alternate diet (Falconer, 1952Go). Genetic correlations between the diets for BW were strong, above 0.71 (Table 6Go). Genetic correlation at time 2 was significantly different from unity, indicating a weak but existent genotype-by-environment interaction. For both DG and DG percentage, genetic correlations at times 1 and 2 were unity or close to it, whereas the correlations at time 4 differed from unity. Taken together, reranking of families for growth traits was only modest.


View this table:
[in this window]
[in a new window]
 
Table 6. Genetic correlations (rA) and their SE between diets for traits
 
For DFI and DFI percentage, the between-diet correlations were lower than for BW (Table 6Go). At time 1, the correlations for DFI and DFI percentage were greater than 0.71. However, at time 2, the correlations were close to zero or negative. At time 4, estimation of the correlations between feed intake traits on both diets was problematic due to the low heritabilities, the correlations being either 0.84 with high standard error (DFI) or nonestimable (DFI percentage). These results demonstrate that feed intake was more prone to reranking than BW, especially at time 2.

Predicted Genetic Responses
Given the estimated genetic parameters, different selection strategies were assessed to analyze the way feed efficiency can be improved by selection for DG and against DFI, or by selection for DG only.

First, 2 sets of phenotypic and genetic (co)variance matrices for BW, DG, and DFI for market-sized fish (> 700 g) were constructed, representing scenarios NP and HP. Because the standard errors of heritabilities for BW, DG, and DFI at times 2 and 4 were overlapping (Table 4Go), the means over both diets were used for both scenarios (h2BW = 0.20; h2DG = 0.41; h2DFI = 0.11). Because phenotypic variances, and phenotypic and genetic correlations between the traits were different in both diets, diet-specific estimates were used. The diet-specific mean estimates were obtained by calculating a mean estimate recorded at times 2 and 4.

Second, genetic gains in response to different selection strategies were predicted for DG and DFI by the standard selection index methodology for individual selection (Hazel, 1943Go; Cameron, 1997Go). Different selection strategies were assessed by switching step-by-step the relative economic weights from DG to DFI. A selection intensity of 1 was assumed. To calculate genetic response in feed efficiency, we first calculated mean feed efficiency (DG/DFI) in the base situation; that is, before selection (based on Table 3Go). Then, responses to selection were calculated for DG and DFI, and mean feed efficiency was recalculated.

The predicted genetic responses revealed that selection solely for rapid DG leads to a 17.6 to 18.6% increase in DG and simultaneously to 8.4 and 9.3% increases in feed efficiency in both NP and HP scenarios, respectively (Figures 1A and BGo). This confirms that increasing DG was related to increasing feed efficiency, and that selection solely for DG will increase feed efficiency. A similar result was observed when BW was selected for (data not shown). When heritability of DFI is much lower than that of DG (Table 4Go), the latter responds much more rapidly to selection, leading to an automatic improvement in feed efficiency, even if the genetic correlation between the traits is high. It should be noted that the absolute levels of genetic gains are arbitrary and greater than expected to occur in reality (Kause et al., 2005Go).


Figure 1
View larger version (19K):
[in this window]
[in a new window]
 
Figure 1. Predicted genetic responses (% of the mean before selection) in daily gain (DG), feed intake (DFI), and feed efficiency (FE) in normal protein (panel A) and high protein (panel B) scenarios after one generation of selection. The selection index includes DG and DFI, which are the 2 components of feed efficiency. The x-axis describes the relative economic weight given to DG. When x = 1, all economic weight is on DG; when x = 0.5, one-half of the weight is for DG and one-half is against DFI; and when x = 0, all weight is against DFI.

 
In the NP scenario, the genetic response of feed efficiency was not increased by switching economic weight from selecting for DG to selecting against DFI (Figure 1AGo). The reason was the very strong unfavorable genetic correlation between DG and DFI, which hampers simultaneous breeding of the traits in opposite directions. However, in the HP scenario, the genetic response in feed efficiency was increased from 9.3 to 11.4% by selecting against DFI (–0.25 economic weight) and for DG (0.75 economic weight; Figure 1BGo). This was a result of the correlations between DG and DFI being only moderate, allowing simultaneous breeding of the traits in opposite directions.


    DISCUSSION
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Both feed costs and nutrient losses into the water from aquaculture can be reduced by improving the feed efficiency of fish. We showed here that selection solely for rapid daily gain could improve feed efficiency as a correlated genetic response. On the experimental diet, simultaneous selection for rapid DG and against DFI, in turn, may increase genetic gain in feed efficiency by a factor of 1.2 compared with selection solely for rapid DG. The results showed, however, that the variation for growth and feed intake and the relationships between the traits might differ in different nutritional environments, influencing genetic responses. Heritabilities and genetic correlations are indeed known to display environment-dependent variation (Falconer, 1952Go; Charmantier and Garant, 2005Go).

Diet Differences in Variation and Correlations
The current industrial on-growing diets for salmonid fish have high lipid (> 35%) and low protein (< 40%) contents. Such feeds are preferred because they support a high growth rate (but also increased lipid deposition), and the lipid component of feed is cheaper than the protein. Our results did not support the hypothesis of decreased variation in BW and DG on the modern NP diet. In contrast, phenotypic variation in growth was greater on the NP diet, and coefficients of genetic variation for DG were greater on the NP diet. No clear trend was seen for the heritabilities of the growth traits. The original hypothesis was based on a belief that the NP diet allows all fish to grow well, leading to low phenotypic variation in growth. The experimental HP diet, in turn, was hypothesized to be a challenge for fish with inefficient protein use, leading to their reduced growth, and thus, to increased variation in growth. Obviously these statements are incorrect, and they need to be replaced with an alternative hypothesis. For example, it is likely that variation in feed intake and growth is more a function of interactions between various nutrients (Hardy, 2002Go), rather than of the individual components, as hypothesized here.

In addition, we showed that both lipid and protein BW displayed greater phenotypic variation (CVP) on the NP diet compared with the HP diet (our unpublished observations). Thus, the increased phenotypic variation in wet BW on the NP diet was a result of both lipid and protein components of growth. Overall, the NP diet increased the expression of diverse BW, and depending on whether this resulted in an increase in residual or genetic variation, trait heritabilities were either reduced or elevated accordingly.

In contrast to the growth traits, DFI displayed greater phenotypic and genetic variation on the HP diet. After an initial reduction in variation on the HP diet just after the diets were applied, the variation in DFI in market-sized fish (times 2 and 4) was greater on the HP than on the NP diet. Moreover, as hypothesized, on the HP diet, growth and DFI were uncoupled as fish grew, whereas on the NP diet, growth and DFI were always strongly positively (unfavorably) correlated. The high positive phenotypic and genetic correlations of growth traits with DFI on NP diet indicate that the individuals with highest BW and DG fed the most. This was in agreement with studies on other farm animals [rA between DG and DFI in beef cattle = 0.68, Koots et al. (1994)Go; in lambs > 0.7, Cammack et al. (2005)Go; in pigs = 0.65, Clutter and Brascamp, (1998)Go], and results mostly from the fact that fish with greater growth feed more. It was not possible without further studies to reveal the mechanisms behind the diet differences in variation and correlations. Although speculative, it is possible that the greater consumption of a normal protein diet by faster growing fish may be partly due to their greater protein needs compared with slower growing fish (e.g., Houlihan et al., 1995Go). When protein content was high (HP diet), fast-growing fish required no extra feeding.

The greater DFI by the fast-growing fish on the NP diet results in greater intake of excess lipid, and consequently, exposes these fish to elevated deposition of lipid stores, as shown by the analysis of body and muscle composition (our unpublished observations). In contrast, for fish fed the HP diet, the genetic correlation between BW and lipid deposition was close to zero or negative, consistent with the diminished correlation between BW and DFI intake on that diet.

It may be argued that differential maturing of NP and HP fish could explain both the greater variation of feed intake and the uncoupling of growth and DFI on the HP diet. However, no support was found for this explanation. First, changing the fixed sex x maturity factor to sex or removing it totally from the statistical models increased phenotypic variation for DFI (and feed intake percentage) but the diet differences in CVP were maintained (data not shown). Second, in the ANOVA results, the interaction of diet with sex x maturity was nonsignificant for all traits (data not shown), indicating that maturity process influenced the traits in a similar manner on both diets. Finally, phenotypic correlations of DFI with BW and DG for different sex x maturity classes remained low on the HP diet and high on the NP diet (data not shown); thus, maturity was not the cause of different correlations on the 2 diets.

Challenges in Genetic Analysis of Feed Efficiency in Fish
Restricted feeding may seem an appealing experimental treatment when studying feed efficiency. In farm animals, restricted feeding is used to reduce among-animal variation in feed intake. In this way, variation in growth is mostly a consequence of variation in feed efficiency, enhancing possibilities to improve feed efficiency by selection solely for rapid growth (Clutter and Brascamp, 1998Go). However, in socially structured fish populations such as rainbow trout, restricted feeding leads to increased dominance hierarchies, and thus, to strongly unequal distribution of feed within a population (McCarthy et al., 1992Go; Jobling, 1995Go; Jobling and Koskela, 1996Go). For instance, coefficient of phenotypic variation for feed intake under restricted feeding is greater than under satiation feeding in trout (McCarthy et al., 1992Go; Jobling and Koskela, 1996Go). Because of feed being unequally distributed between individuals during restricted feeding, it does not provide an efficient alternative to study feed efficiency in fish. Accordingly, when planning the current study, diets with different protein content were considered a more viable alternative.

Silverstein et al. (2001)Go, working on catfish, were the first to report on individual feed intake in a family-structured fish population. They found a broad-sense heritability of 0.37 to 0.41 for feed intake. However, the majority of genetic studies on feed intake and feed efficiency in fish have been based on the average performance of full-sib families held in family tanks (Kinghorn, 1983Go; Thodesen et al., 2001Go; Doupé and Lymbery, 2004Go; Kolstad et al., 2004Go), which does not allow reliable estimation of genetic parameters. This results in elevating heritabilities because calculating tank means removes the large within-family variance from the data. Accordingly, we found an average heritability of 0.10 for DFI, whereas in the previous studies, heritabilities or proportion of variation due to family structure have been 0.31 to 0.84 (Kinghorn, 1983Go; Thodesen et al., 2001Go; Kolstad et al., 2004Go).

The weakness of our approach was the inaccuracy of the x-ray method to describe true long-term feed intake. This results from feed intake displaying large day-today variation (McCarthy et al., 1992Go; Jobling and Koskela, 1996Go), and thus, the correspondence between the short-term x-ray method and long-term recording was low (our unpublished observations). Moreover, the repeatability of DFI, even for the mean of 3 repeated records, was low to moderate, increasing residual variation. However, the x-ray method is the only method to record individual feed intake from large numbers of fish held in a common tank, and is routinely used for other fish species as well (McCarthy et al., 1993Go; Jobling et al., 2001Go; Silverstein et al., 2001Go).

Selection Strategies for Increasing Feed Efficiency
The analysis of predicted genetic responses in feed efficiency showed that selection solely for rapid growth improves feed efficiency as a correlated genetic response. This observation was in agreement with some (e.g., Thodesen et al., 1999Go; Ogata et al., 2002Go) but not all (Mambrini et al., 2004Go) fish studies examining correlated genetic changes in feed efficiency in response to selection for high BW, and with studies on fish and farm animals showing that growth and feed efficiency are favorably genetically correlated (Kinghorn, 1983Go; Koots et al., 1994Go; Clutter and Brascamp, 1998Go, Thodesen et al., 2001Go; Henryon et al., 2002Go; Doupé and Lymbery, 2004Go; Kolstad et al., 2004Go; Cammack et al., 2005Go).

Interestingly, in the NP scenario, the genetic gain in feed efficiency was not increased when DFI was selected against, along with selection for high DG. In contrast, on the HP diet, selection against DFI considerably increased genetic gain in feed efficiency. Moreover, a potentially realistic scenario was simulated in which heritability of DFI was artificially elevated to 0.25 for the HP scenario, all other estimates remaining the same as previously described. In this scenario, genetic response in feed efficiency was tripled by selecting against DFI. This scenario may be realistic because novel methods for recording feed intake in the future may increase the accuracy of records, thus elevating heritability values.

Great profits could be gained by improving feed efficiency through selection. Regardless of the intuitive benefits, none of the existing fish-breeding programs uses feed intake or feed efficiency in their selection index. Instead, they select for rapid growth or for high BW at a fixed age. The reason for the lack of selection efforts to improve feed efficiency is that individual feed intake is challenging and costly to record simultaneously from a large population of fish. Moreover, the favorable genetic correlation between BW and feed efficiency reduces the motivation to select directly for feed efficiency, although it could be advantageous.

Could an HP diet be used as a novel selection environment to improve feed efficiency on current commercial diets? In a seminal paper, Falconer (1952)Go concluded that to maximize a selection response, it is sometimes beneficial to select in an alternative environment, even when this environment is not a commercial production environment. Falconer (1952)Go showed that the greater the heritability in the alternative environment and the greater the genetic correlation between the environments, the more beneficial is the indirect selection in the alternative environment. In our study, DFI (and feed efficiency and residual feed intake, data not shown) displayed moderate reranking of families across the diets, suggesting that selection on the alternative HP diet does not lead to parallel changes on the commercial diet. Environment-specific expression of heritabilities and genetic correlations of feed intake and growth are utilized, for example, in pig breeding (Clutter and Brascamp, 1998Go). The present evidence, however, is not strong enough that such practice could be implemented in rainbow trout.


    Footnotes
 
1 The authors acknowledge Tuija Paananen, Heli Wahlroos, Ángel Ruiz Bañón, César Luis Borreguero, and Alexandre Dobly for assistance in fish management, data collection, and preparation. This work was carried out with financial support from the Commission of the European Communities, Quality of Life and Management of Living Resources Programme, project Q5RS-2001-0994 "Protein and Growth Efficiency in Salmonid Selection (PROGRESS)." It does not reflect its views and in no way anticipates the Commission’s future policy on this area. Back

2 Corresponding author: Antti.Kause{at}mtt.fi

Received for publication September 6, 2005. Accepted for publication November 11, 2005.


    LITERATURE CITED
 Top
 Abstract
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 


Cameron, N. D. 1997. Selection Indices and Prediction of Genetic Merit in Animal Breeding. CAB International, Wallingford, UK.

Cammack, K. M., K. A. Leymaster, T. G. Jenkins, and M. K. Nielsen. 2005. Estimates of genetic parameters for feed intake, feeding behavior, and daily gain in composite ram lambs. J. Anim. Sci. 83:777–785.[Abstract/Free Full Text]

Charmantier, A., and D. Garant. 2005. Environmental quality and evolutionary potential: Lessons from wild populations. Proc. R. Soc. Lond. B, Biol. Sci. 272:1415–1425.

Cho, C. Y. 1992. Feeding systems for rainbow trout and other salmonids with reference to current estimates of energy and protein requirements. Aquaculture 100:107–123.

Clutter, A. C., and E. W. Brascamp. 1998. Genetics of performance traits. Pages 427–462 in The Genetics of the Pig. M. F. Rothschild and A. Ruvinsky, ed. CAB International, Wallingford, UK.

Doupé, R. G., and A. J. Lymbery. 2004. Indicators of genetic variation for feed conversion efficiency in black bream. Aquac. Res. 35:1305–1309.

Falconer, D. S. 1952. The problem of environment and selection. Am. Nat. 86:293–298.

Falconer, D. S., and T. F. C. Mackay. 1996. Introduction to Quantitative Genetics. 4th ed. Longman Group Ltd., London, UK.

Gunsett, F. C. 1984. Linear index selection to improve traits defined as ratios. J. Anim. Sci. 59:1185–1193.[Abstract/Free Full Text]

Hardy, R. W. 2002. Rainbow trout, Oncorhynchus mykiss. Pages 184–202 in Nutrient Requirements and Feeding of Finfish for Aquaculture. C. D. Webster and C. E. Lim, ed. CABI Publishing, Wallingford, UK.

Hazel, L. N. 1943. The genetic basis of constructing selection indices. Genetics 28:476–490.[Free Full Text]

Henryon, M., A. Jokumsen, P. Berg, I. Lund, P. B. Pedersen, N. J. Olesen, and W. J. Slierendrecht. 2002. Genetic variation for growth rate, feed conversion efficiency, and disease resistance exists within a farmed population of rainbow trout. Aquaculture 209:59–76.

Houle, D. 1992. Comparing evolvability and variability of quantitative traits. Genetics 130:195–204.[Abstract]

Houlihan, D. F., C. G. Carter, and I. D. McCarthy. 1995. Protein synthesis in fish. Pages 191–219 in Biochemistry and Molecular Biology of Fishes. Vol. 4. P. Hochachka and P. Mommsen, ed. Elsevier Science, Amsterdam, The Netherlands.

Iwama, G. K., and A. F. Tautz. 1981. A simple growth model for salmonids in hatcheries. Can. J. Fish. Aquac. Sci. 38:649–656.

Jensen, J., and P. Madsen. 2000. A User’s Guide to DMU. A Package for Analyzing Multivariate Mixed Models. National Institute of Animal Science, Research Center Foulum, Denmark.

Jobling, M. 1995. Simple indices for the assessment of the influences of social environment on growth performance, exemplified by studies on Arctic charr. Aquac. Int. 3:60–65.

Jobling, M., D. Covès, B. Damsga °rd, H. R. Kristiansen, J. Koskela, T. E. Petursdottir, S. Kadri, and O. Gudmundsson. 2001. Techniques for measuring feed intake. Pages 49–87 in Food Intake in Fish. D. Houlihan, T. Boujard, and M. Jobling, ed. Blackwell Science, Oxford, UK.

Jobling, M., and J. Koskela. 1996. Interindividual variations in feeding and growth in rainbow trout during restricted feeding and a subsequent period of compensatory growth. J. Fish Biol. 49:658–667.

Kause, A., O. Ritola, T. Paananen, H. Wahlroos, and E. A. Mäntysaari. 2005. Genetic trends in growth, sexual maturity and skeletal deformations, and rate of inbreeding in a breeding programme for rainbow trout. Aquaculture 247:177–187.

Kennedy, B. W., J. H. J. van der Werf, and T. H. E. Meuwissen. 1993. Genetic and statistical properties of residual feed intake. J. Anim. Sci. 71:3239–3250.[Abstract]

Kenward, M. G., and J. H. Roger. 1997. Small sample inference for fixed effects from restricted maximum likelihood. Biometrics 53:983–997.[Medline]

Kim, J. D., and S. J. Kaushik. 1992. Contribution of digestible energy from carbohydrates and estimations of protein/energy requirements for growth of rainbow trout (Oncorhynchus mykiss). Aquaculture 106:161–169.

Kinghorn, B. 1983. Genetic variation in food conversion efficiency and growth in rainbow trout. Aquaculture 32:141–155.

Kolstad, K., B. Grisdale-Helland, and B. Gjerde. 2004. Family differences in feed efficiency in Atlantic salmon (Salmo salar). Aquaculture 241:169–177.

Koots, K. R., J. P. Gibson, and J. W. Wilton. 1994. Analyses of published genetic parameter estimates for beef production traits. 2. Phenotypic and genetic correlations. Anim. Breed. Abstr. 62:309–338.

Mambrini, M., F. Médale, M. P. Sanchez, B. Recalde, B. Chevassus, L. Labbé, E. Quillet, and T. Boujard. 2004. Selection for growth in brown trout increases feed intake capacity without affecting maintenance and growth requirements. J. Anim. Sci. 82:2865–2875.[Abstract/Free Full Text]

McCarthy, I. D., C. G. Carter, and D. F. Houlihan. 1992. The effect of feeding hierarchy on individual variability in daily feeding of rainbow trout, Oncorhynchus mykiss (Walbaum). J. Fish Biol. 41:257–263.

McCarthy, I. D., D. F. Houlihan, C. G. Carter, and K. Moutou. 1993. Variation in individual consumption rates of fish and its implications for the study of fish nutrition and physiology. Proc. Nutr. Soc. 52:427–436.[Medline]

NRC. 1993. Nutrient Requirements of Fish. National Academic Press, Washington, DC.

Ogata, H. Y., H. Oku, and T. Murai. 2002. Growth, feed efficiency and feed intake of offspring from selected and wild Japanese flounder (Paralichthys olivaceus). Aquaculture 211:183–193.

Silverstein, J. T., B. G. Bosworth, G. C. Waldbieser, and W. R. Wolters. 2001. Feed intake in channel cat fish: Is there a genetic component? Aquac. Res. 32:199–205.

Talbot, C., and P. J. Higgins. 1983. A radiographic method for feeding studies on fish using metallic iron powder as marker. J. Fish Biol. 23:211–220.

Thodesen, J., B. Gjerde, B. Grisdale-Helland, and T. Storebakken. 2001. Genetic variation in feed intake, growth and feed utilization in Atlantic salmon (Salmo salar). Aquaculture 194:273–281.

Thodesen, J., B. Grisdale-Helland, S. J. Helland, and B. Gjerde. 1999. Feed intake, growth and feed utilization of offspring from wild and selected Atlantic salmon (Salmo salar). Aquaculture 180:237–246.


This article has been cited by other articles:


Home page
J ANIM SCIHome page
C. D. Quinton, A. Kause, K. Ruohonen, and J. Koskela
Genetic relationships of body composition and feed utilization traits in European whitefish (Coregonus lavaretus L.) and implications for selective breeding in fishmeal- and soybean meal-based diet environments
J Anim Sci, December 1, 2007; 85(12): 3198 - 3208.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kause, A.
Right arrow Articles by Ruohonen, K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kause, A.
Right arrow Articles by Ruohonen, K.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS