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

HAL
open science

# Derivation of economic values for production traits in aquaculture species

Kasper Janssen, Paul Berentsen, Mathieu Besson, Hans Komen

► To cite this version:

Kasper Janssen, Paul Berentsen, Mathieu Besson, Hans Komen. Derivation of economic values for production traits in aquaculture species. Genetics Selection Evolution, 2017, 49 (1), pp.5. 10.1186/s12711-016-0278-x. http://hal-01479138

HAL Id: hal-01479138

https://hal.science/hal-01479138

Submitted on 28 Feb 2017

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

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

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

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

# Derivation of economic values for production traits in aquaculture species

Kasper Janssen¹*, Paul Berentsen², Mathieu Besson¹,³ and Hans Komen¹

**Abstract**

**Background:** In breeding programs for aquaculture species, breeding goal traits are often weighted based on the desired gains but economic gain would be higher if economic values were used instead. The objectives of this study were: (1) to develop a bio-economic model to derive economic values for aquaculture species, (2) to apply the model to determine the economic importance and economic values of traits in a case-study on gilthead seabream, and (3) to validate the model by comparison with a profit equation for a simplified production system.

**Methods:** A bio-economic model was developed to simulate a grow-out farm for gilthead seabream, and then used to simulate gross margin at the current levels of the traits and after one genetic standard deviation change in each trait with the other traits remaining unchanged. Economic values were derived for the traits included in the breeding goal: thermal growth coefficient (TGC), thermal feed intake coefficient (TFC), mortality rate (M), and standard deviation of harvest weight ($\sigma_{HW}$). For a simplified production system, improvement in TGC was assumed to affect harvest weight instead of growing period. Using the bio-economic model and a profit equation, economic values were derived for harvest weight, cumulative feed intake at harvest, and overall survival.

**Results:** Changes in gross margin showed that the order of economic importance of the traits was: TGC, TFC, M, and $\sigma_{HW}$. Economic values in € (kg production)$^{-1}$ (trait unit)$^{-1}$ were: 0.40 for TGC, -0.45 for TFC, -7.7 for M, and -0.0011 to -0.0010 for $\sigma_{HW}$. For the simplified production system, similar economic values were obtained with the bio-economic model and the profit equation. The advantage of the profit equation is its simplicity, while that of the bio-economic model is that it can be applied to any aquaculture species, because it can include any limiting factor and/or environmental condition that affects production.

**Conclusions:** We confirmed the validity of the bio-economic model. TGC is the most important trait to improve, followed by TFC and M, and the effect of $\sigma_{HW}$ on gross margin is small.

**Background**

In Europe, over 80% of aquaculture production originates from breeding programs, which in most cases apply family selection with the aim of improving multiple traits simultaneously [1]. Breeding goal traits are often weighted based on desired gains rather than economic values [2], which compromises economic gain [3, 4].

In aquaculture species, economic values are available
only for a few species, although their importance has

repeatedly been underlined, e.g. [5, 6]. Profit equations
have been used to derive economic values for Nile tilapia
(*Oreochromis niloticus*) [7], common carp (*Cyprinus carpio*) [8], Australian abalones (*Haliotis rubra* and *H. laevigata*) [9], and crayfish (*Cherax tenuimanus*) [10]. Besson
et al. [11, 12] used a bio-economic model to derive eco-
nomic values for African catfish (*Clarias gariepinus*) that
were produced in a land-based aquaculture system in
which water is treated and recirculated and for growth
rate in European seabass (*Dicentrarchus labrax*) under
varying temperature conditions, respectively.

For livestock species with simple and highly controlled
production systems, such as pig production, economic
values can be derived from a profit equation, e.g. [13]. For

<footnote>*Correspondence: kasper.janssen@wur.nl
¹Animal Breeding and Genomics, Wageningen University and Research, Droevendaalsesteeg 1, 6708 PB Wageningen, The Netherlands
Full list of author information is available at the end of the article</footnote>
production systems with a higher degree of complexity, partly due to seasonal variation, such as in dairy cattle and sheep farming, profit equations may fail to provide an adequate description of the farming system and bio-economic models are required [14–16]. In general, bio-economic models provide a more accurate description of farming systems than profit equations and are, therefore, increasingly used to estimate economic values [2].

Fish farms are complex production systems for two reasons. First, fish are kept outdoors in most farming systems and, thus, are exposed to fluctuating environmental conditions. Seasonal variation in temperature causes variation in growth rate of fish, because of their ectothermic nature. Fish are harvested at a constant weight rather than at a constant age, hence the length of a production cycle depends on the stocking date. Second, production output of a farm is determined by constraints such as oxygen availability [11] and stocking density. Stocking density constrains production output for many important aquaculture species, including Atlantic salmon (*Salmo salar*), European seabass, and gilthead seabream (*Sparus aurata*). Thus, bio-economic models could prove useful to derive economic values for aquaculture species.

The objectives of this study were: (1) to develop such a bio-economic model, (2) to apply the model to determine the economic importance and economic values for various traits in a case-study on gilthead seabream, and (3) to validate the model by comparison with a profit equation for a simplified production system.

# Methods

## Traits

The breeding goal considered here includes growth rate, feed intake rate, mortality rate, and uniformity in harvest weight. Growth rate affects revenues, feed costs and juvenile costs; feed intake rate affects feed costs; mortality rate affects feed costs and juvenile costs. Feed and juveniles are major costs in production [17], thus including growth and feed intake in the breeding goal is common practice in livestock [18, 19]. Uniformity, i.e. size variation around the mean harvest weight, determines the distribution of fish over price categories at harvest, and thus affects revenues via the average sales price of fish.

Economic values are specific for the unit in which a trait is expressed [20]. Here, growth rate is expressed in units of thermal growth coefficient (TGC) [21]. TGC is a standardized measure of growth in fish that takes stocking weight and temperature variation over the lifespan of a fish into account. TGC is widely used [22], is more accurate than other measures of growth rate [23], and is relatively robust to differences in temperature regimes [24, 25]. Feed intake is assumed to be determined by the same variables as bodyweight and gain in bodyweight,

because energy requirement is largely determined by bodyweight and gain in bodyweight [26]. In this study, feed intake rate is, therefore, expressed in units of thermal feed intake coefficient TFC, a TGC analogue. TFC takes stocking weight and temperature variation over the lifespan of a fish into account. The TFC model is independent of the TGC model, i.e. a modelled change in growth rate will not affect modelled feed intake rate and vice versa, which is a prerequisite to derive economic values. Mortality rate (M) is expressed as % of mortality per day. Uniformity is expressed as the standard deviation of harvest weight ($\sigma_{HW}$) in grams.

## Bio-economic model

The bio-economic model developed by Besson et al. [12] was adapted to simulate production systems with seasonal variation in temperature and in which density constrains production output, in this case a typical grow-out farm for gilthead seabream in Greece. The model is a deterministic simulation model that is programmed in R version 2.12.2. [27]. The model simulates operation of the farm during an average year. The farm consists of 20 cages of 2800 m³ each and produces about 550 tons annually. As shown in Fig. 1, the model consists of three hierarchical parts: a fish model, a cage model, and a farm model. Inputs into the fish model are: stocking date, temperature coefficients, TGC, and TFC. Outputs of the fish model for each stocking date are: bodyweight per day per fish, feed consumption per day per fish, and harvest date. Inputs into the cage model are: outputs of the fish model, M, cage volume, and feed prices. Outputs of the cage model for each stocking date are per production cycle of a cage: fish production, number of juveniles stocked, feed consumption, and feed costs. A production cycle is the period between stocking and harvesting a cage. Inputs of the farm model are: outputs of the cage model, number of cages, price of juveniles, price of packing, and sales prices. Outputs of the farm model are total per year: fish production, number of juveniles stocked, feed consumption, feed costs, juvenile costs, packing costs, revenues from fish sales, and gross margin.

The model was used to derive economic values of the traits mentioned above. Economic values give the expected change in profit from a small change in trait level, keeping the level of all other traits constant. Genetic change does not affect fixed costs, hence change in profit due to genetic change equals change in gross margin. When change in trait level equals the additive genetic standard deviation ($\sigma_A$) [15], the resulting change in gross margin indicates how important that trait is, because $\sigma_A$ is indicative of the rate at which breeding values can be improved [28]. To determine the relative importance of the traits and to derive economic values, the model was
## Farm model

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Schematic overview of the bio-economic model. TGC thermal growth coefficient, TFC thermal feed intake coefficient, BW bodyweight, CFI cumulative feed intake, DFI daily feed intake, M mortality rate</figcaption></figure>

run under two situations: first before genetic change and second, after a change of one $\sigma_A$ in one trait with the other traits kept constant. For the trait ‘uniformity’, minimum and maximum values of the possible range of $\sigma_A$ were used, because the actual value was unknown. Trait levels were changed in the desired direction of the genetic change. Economic values were expressed per kg of fish produced in the situation before genetic change [29] and were calculated as:

Economic value

$$= \left( \frac{\text{gross margin}_A - \text{gross margin}_B}{\text{trait level}_A - \text{trait level}_B} \right) / \text{fish production}_B, \quad (1)$$

where subscripts indicate before (B) and after (A) genetic change.

## Model equations

The “Estimation of model coefficients” section describes the derivation of model coefficients from farm data. These coefficients (see Table 1) are used as the equations’ coefficients of the bio-economic model. The three parts of the bio-economic model are described in the following three subsections.

### Estimation of model coefficients

Coefficients in the equations to describe temperature, fish growth, feed intake, and number of fish per cage were

derived from recent farm data of the company Andromeda S.A., which is hereafter referred to as ‘data’. The data included daily records of temperature, feed provided, and mortality, and regular records of bodyweight from 15 cages of a farm located in Vonitsa, Northern Greece during the period 2013 through 2015.

Seasonal variation in daily water temperature throughout the year followed a sinusoidal pattern. Therefore, the equation to describe daily temperature ($T_{s,a}$) in °C was [30]:

$$T_{s,a} = TM - TA \cdot \sin(2\pi \cdot (t_{s,a} - t_A)/365), \quad (2)$$

where TM is the average annual temperature (°C), TA is the range of temperatures around TM (°C), $t_A$ is the time of the year at which $T_{s,a}$ equaled TM, and $t_{s,a}$ represents the date defined as:

$$t_{s,a} = s + a, \quad (3)$$

where stocking date $s$ ($s = 1, \dots, n$) equals 1 on January 1 2013 and $a$ ($a = 0, \dots, n$) is the age of the fish (days). To estimate TM, TA, and $t_A$, Eq. 2 was fitted to the data by means of non-linear least-squares regression in R. Table 1 shows the resulting coefficients.

Bodyweight in seabream can be predicted from the stocking weight and the sum of daily effective temperatures. Daily effective temperature is the daily temperature minus 12 °C, where 12 °C represents the minimum temperature for seabream growth [31]. Therefore, the
<table><caption>Table 1 Estimated coefficients used in model equations</caption><thead><tr><th>Symbol</th><th>Meaning</th><th>Value</th><th>Standard error</th><th>Unit</th></tr></thead><tbody><tr><td>TM</td><td>Annual mean temperature</td><td>19.57</td><td>0.0119</td><td>°C</td></tr><tr><td>TA</td><td>Amplitude of temperature</td><td>-4.806</td><td>0.0167</td><td>°C</td></tr><tr><td>t<sub>A</sub></td><td>Date at which temperature equals TM</td><td>-32.49</td><td>0.2057</td><td>Day</td></tr><tr><td>TGC</td><td>Thermal growth coefficient</td><td>12.6</td><td>0.0847</td><td>g<sup>2/3</sup>/(day degrees · 1000)</td></tr><tr><td>TFC</td><td>Thermal feed intake coefficient</td><td>8.25</td><td>0.157</td><td>g<sup>0.544</sup>/(day degrees · 1000)</td></tr><tr><td>p</td><td>Weight exponent to predict cumulative feed intake</td><td>0.544</td><td>0.00282</td><td>–</td></tr><tr><td>M</td><td>Mortality rate</td><td>0.0300</td><td>0.000164</td><td>%/day</td></tr></tbody></table>

equation to describe bodyweight (in g) at $t_{s,a}$ ($BW_{t_{s,a}}$) was [32]:

$$BW_{t_{s,a}} = \left( BW_{t_{s,0}}^{2/3} + \frac{TGC}{1000} \cdot \sum_{i=t_{s,0}}^{(t_{s,a})-1} (T_i - 12) \right)^{3/2}, \quad (4)$$

where $BW_{t_{s,0}}$ is bodyweight at stocking (in g), and $\sum_{i=t_{s,0}}^{(t_{s,a})-1} (T_i - 12)$ represents the sum of effective temperatures (day degrees) over the lifespan of a fish excluding $t_{s,a}$. To estimate $TGC$, Eq. 4 was fitted to the data by means of non-linear least-squares regression in R. Instead of fixing exponents 2/3 and 3/2, fitting these to the data resulted in values of 0.612 and 1/0.612 but this barely improved accuracy of the model. Values of 2/3 and 3/2 were preferred, because standardization of growth models allows for a better comparison of growth rate across studies. Analogous to $BW_{t_{s,a}}$, the model to describe cumulative feed intake (in g) at $t_{s,a}$ ($CFI_{t_{s,a}}$) was:

$$CFI_{t_{s,a}} = \left( BW_{t_{s,0}}^{p} + \frac{TFC}{1000} \cdot \sum_{i=t_{s,0}}^{t_{s,a}} (T_i - 12) \right)^{1/p} - BW_{t_{s,0}}, \quad (5)$$

where $p$ is a weight exponent, and $\sum_{i=t_{s,0}}^{t_{s,a}} (T_i - 12)$ represents the sum of effective temperatures (day degrees) over the lifespan of a fish, including $t_{s,a}$. The term $BW_{t_{s,0}}$ was subtracted from $(BW_{t_{s,0}}^p + \frac{TFC}{1000} \cdot \sum_{i=t_{s,0}}^{t_{s,a}} (T_i - 12))^{1/p}$ to force the model through the intercept. To estimate $TFC$ and $p$, Eq. 5 was fitted to the data by means of non-linear least-squares regression in R. Parameter $M$ (%/day) was assumed to be constant over time, hence the number of fish alive decreased exponentially in time. The model to describe the number of fish at $t_{s,a}$ ($N_{t_{s,a}}$) was:

$$N_{t_{s,a}} = N_{t_{s,0}} \cdot \left( 1 - \frac{M}{100\%} \right)^{(t_{s,a})-1}, \quad (6)$$

where $N_{t_{s,0}}$ is the number of fish stocked. To estimate $M$, Eq. 6 was fitted to the data by means of non-linear least-squares regression in R.

### Fish model

Date ($t_{s,a}$) was modelled as in Eq. 3. Temperature ($T_{t_{s,a}}$) was modelled as in Eq. 2. Bodyweight ($BW_{t_{s,a}}$) was modelled as in Eq. 4, where stocking weight ($BW_{t_{s,0}}$) was 4.4 g, equal to the average in the data. Equation 4 was rewritten to calculate the harvest date ($t_{s,h}$) as:

$$\sum_{i=t_{s,0}}^{(t_{s,h})-1} (T_i - 12) = \frac{BW_{t_{s,h}}^{2/3} - BW_{t_{s,0}}^{2/3}}{TGC/1000}, \quad (7)$$

where $BW_{t_{s,h}}$ is the average harvest weight, here set to the desired market weight of 400 g. Solving the right hand side of the equation, yields $\sum_{i=t_{s,0}}^{(t_{s,h})-1} (T_i - 12) = 4084$ day degrees. Cumulative feed intake ($CFI_{t_{s,a}}$) was modelled as in Eq. 5. $CFI_{t_{s,h}}$ was set equal to $CFI_{(t_{s,h})-1}$, because fish are not fed on the day that they are harvested. Daily feed intake ($DFI_{t_{s,a}}$) was modeled as:

$$DFI_{t_{s,a}} = CFI_{(t_{s,a})+1} - CFI_{t_{s,a}}, \quad (8)$$

### Cage model

To maximize production, standing stock in a cage reaches the maximum allowable density of 15 kg/m³ at harvest. For the 2800 m³ cages, per production cycle fish production is thus 42,000 kg or 100,500 fish. To compensate for mortality, a larger number of fish is stocked than harvested. The number of fish in a cage at $t_{s,a}$ ($N_{t_{s,a}}$) was modeled as:

$$N_{t_{s,a}} = 100,500 \cdot \left( 1 - \frac{M}{100\%} \right)^{(t_{s,h} - t_{s,a})}, \quad (9)$$

The number of juveniles stocked per cage ($N_{t_{s,0}}$) was calculated by substituting $t_{s,a}$ for $t_{s,0}$. Daily feed intake per cage (kg) at $t_{s,a}$ ($DFI_{cage_{t_{s,a}}}$) was modeled as:
$$DFIcage_{t_{s,a}} = N_{t_{s,a}} \cdot DFI_{t_{s,a}}/1000. \qquad (10)$$

Total feed consumption per cage (kg) stocked at date $t_{s,0}$ ($TFIcage_{t_{s,0}}$) was calculated as:

$$TFIcage_{t_{s,0}} = \sum_{i=t_{s,0}}^{t_{s,h}} DFIcage_i. \qquad (11)$$

Depending on bodyweight, fish are fed different feed types. Daily feed costs per cage ($DFCcage_{t_{s,a}}$) were calculated as the product of $DFIcage_{t_{s,a}}$ and feed price per size category. Table 2 shows the price of feed types for the fish size categories. Total feed costs per cage (€) stocked at date $t_{s,0}$ ($TFCcage_{t_{s,0}}$) were calculated as:

$$TFCcage_{t_{s,0}} = \sum_{i=t_{s,0}}^{t_{s,h}} DFCcage_i. \qquad (12)$$

### Farm model

For the farm model, the cage model was run repeatedly over the whole range of stocking dates, from 1 (January 1th) to 365 days (December 31th) by one-day steps. Per stocking date, age at harvest ($t_{s,h} - t_{s,0}$), $TFIcage_{t_{s,0}}$, $TFCcage_{t_{s,0}}$ and $N_{t_{s,0}}$ were calculated. These results were averaged over all stocking dates to compute the average production cycle of a cage. The period between two successive production cycles is three days [Personal communications Andromeda S.A., 2015]. For the 20 cages that are present on the farm, the number of production cycles per year was calculated as:

$$Production cycles per year = 20 \cdot \frac{365}{average(t_{s,h} - t_{s,0}) + 3}. \qquad (13)$$

Average results per production cycle were multiplied by the number of production cycles per year to compute outputs at the farm level per year: fish production, number of juveniles stocked, feed consumption, and feed costs. Juvenile costs at the farm level were calculated by multiplying the number of juveniles stocked by the price of juveniles of €0.20 per piece. Packing costs at the farm

level were calculated by multiplying fish production by the price of packing of €0.33 per kg fish. Revenues from fish sales at the farm level were calculated as the product of fish production and average sales price. Average sales price was computed as the proportion of fish of each category in Table 3 multiplied by its corresponding sales price. Three percent of the fish harvested are deformed. The size distribution of the remaining fish was calculated from a normal probability density function with $\mu = 400$ and $\sigma_{HW} = 60$ [Personal communications Andromeda S.A., 2015].

### Additive genetic standard deviation of traits

The economic importance of each trait in the breeding goal depends on the change in gross margin, which itself depends on the change in trait level of one $\sigma_A$. The genetic coefficient of variation (CVA) can be used to estimate $\sigma_A$ from the mean trait level ($\mu$) [28]:

$$\sigma_A = \frac{CV_A}{100\%} \cdot \mu. \qquad (14)$$

For TFC, $\sigma_A$ can be estimated from the genetic variation in $BW_{t_{s,h}}$. For this study, $BW_{t_{s,h}}$ was set equal to 400 g, and $\sum_{i=t_{s,0}}^{(t_{s,h})-1} (T_i - 12)$ was 4084 day degrees (Eq. 7). CVA of bodyweight was estimated to be 10.6% based on data from Navarro et al. [33]. For $BW_{t_{s,h}}$, $\sigma_A$ was thus 42.4 g. The distribution of $BW_{t_{s,h}}$ was simulated in R as $BW_{t_{s,h}} = \mu + z_n \cdot \sigma_A$, where $\mu = 400$ and $\sigma_A = 42.4$, and $z_n$ is a standard normal distribution ($z_n \sim N(0, 1)$) with $n = 1, \dots, 10^6$. From this simulation, $\sigma_A$ of TGC was estimated as (Appendix 1):

$$\sigma_A \text{ of TGC} \approx \sqrt{\text{Var}(\text{TGC}_n)} \approx \sqrt{\left(\frac{1000}{4084}\right)^2 \cdot \text{Var}(BW_{t_{s,h}}^{2/3})} \\ = \frac{1000}{4084} \sqrt{\frac{1}{10^6 - 1} \cdot \sum_{i=1}^{n} (BW_i^{2/3} - 400^{2/3})^2} \\ = 0.95 \text{ g}^{2/3}/(\text{day degrees} \cdot 1000). \qquad (15)$$

<caption>Table 3 Sales price per fish size category in 2014 [Personal communications Andromeda S.A., 2015]</caption>

<table><caption>Table 2 Feed price per fish size category in 2014 [Personal communications Andromeda S.A., 2015]</caption><thead><tr><th>Fish size (g)</th><th>Price (€/kg)</th></tr></thead><tbody><tr><td>&lt;7</td><td>2.21</td></tr><tr><td>7-13</td><td>1.97</td></tr><tr><td>13-30</td><td>1.65</td></tr><tr><td>30-80</td><td>1.26</td></tr><tr><td>80-300</td><td>1.12</td></tr><tr><td>&gt;300</td><td>1.17</td></tr></tbody></table>

<table><thead><tr><th>Category (g)</th><th>Price (€/kg)</th></tr></thead><tbody><tr><td>&lt;100</td><td>0</td></tr><tr><td>100-200</td><td>1.65</td></tr><tr><td>200-300</td><td>4.15</td></tr><tr><td>300-400</td><td>4.52</td></tr><tr><td>400-600</td><td>4.63</td></tr><tr><td>&gt;600 g</td><td>5.27</td></tr><tr><td>Deformed</td><td>2.52</td></tr></tbody></table>
For TFC, $\sigma_A$ can be estimated from the genetic variation in both $BW_{t,0}$ and $CFI_{t,s,h}$. For $CFI_{t,s,h}$, $\sigma_A$ can be approximated by (Appendix 2):

$$ \sigma_A \text{ of } CFI_{t,s,h} \approx \frac{1.75}{r_A} \cdot \sigma_A \text{ of } BW_{t,s,h}, \quad (16) $$

where $r_A$ is the genetic correlation between $BW_{t,s,h}$ and $CFI_{t,s,h}$, which was assumed to be 0.90 [34]. Solving Eq. (16), $\sigma_A$ of $CFI_{t,s,h}$ was equal to 82 g. Based on an average $CFI_{t,s,h}$ of 713 g in our study, the CVA of $CFI_{t,s,h}$ was 12%, which is close to values reported for other species [34, 35]. Genetic variances for $CFI_{t,s,h}$ and $BW_{t,s,h}$ were simulated to calculate $\sigma_A$ of TFC. In our study, the average $BW_{t,s,h}$ was 4.4 g. Based on a CVA of 10.6% for bodyweight, $\sigma_A$ of $BW_{t,s,h}$ was equal to 0.45 g. The distribution of $CFI_{t,s,h}$ was simulated in R as $CFI_{t,s,h} = \mu + z_n \cdot \sigma_A$, where $\mu = 713$, $\sigma_A = 82$, and $z_n$ is a standard normal distribution ($z_n \sim N(0, 1)$) with $n = 1, \dots, 10^6$. The distribution of $BW_{t,s,h}$ was simulated in R as $BW_{t,s,h} = \mu + z_n \cdot \sigma_A$, where $\mu = 4.4$, $\sigma_A = 0.45$ and $z_n$ is a standard normal distribution ($z_n \sim N(0, 1)$) with $n = 1, \dots, 10^6$. A covariance of zero was assumed between $CFI_{t,s,h}$ and $BW_{t,s,h}$. Based on the simulations, $\sigma_A$ of TFC was estimated as:

$$ \begin{aligned} \sigma_A \text{ of TFC} &\approx \sqrt{\text{Var}(\text{TFC}_n)} \\ &= \sqrt{\text{Var}\left(1000 \cdot \frac{(CFI_{t,s,h} + BW_{t,s,h})^p - BW_{t,s,h}^p}{4084}\right)} \\ &= \frac{1000}{4084} \cdot \sqrt{\frac{1}{10^6 - 1} \cdot \sum_{i=1}^{n} \left( \frac{(CFI_i + BW_i)^{0.544} - BW_i^{0.544}}{-(713 + 4.4)^{0.544} - 4.4^{0.544}} \right)^2} \\ &= 0.55 \text{ g}^{0.544}/(\text{degree days} \cdot 1000). \end{aligned} \quad (17) $$

For $M$, $\sigma_A$ can be estimated from the genetic variation in cumulative mortality at harvest ($CM_{t,s,h}$). The average $CM_{t,s,h}$ was 14.9%. In animal breeding, an underlying liability scale is commonly used to analyze mortality and survival [36]. Heritability of $CM_{t,s,h}$ on the liability scale was assumed to be 0.17 [37] and by definition $\sigma_P$ is equal to 1, hence $\sigma_A = \sqrt{h^2} = \sqrt{0.17}$. Before genetic change, the deviation of the threshold from the mean ($x_B$) was calculated from the quantile function of a normal distribution in R as $x_B = -qnorm(0.149) = 1.04$. After genetic change by one $\sigma_A$, the deviation from the threshold from the mean ($x_A$) becomes: $x_A = x_B + \sigma_A = 1.04 + \sqrt{0.17} = 1.45$. After genetic change, $CM_{t,s,h}$ was calculated from the distribution function of a normal distribution in R as:

$$ \begin{aligned} CM_{t,s,h} &= (1 - pnorm(x_A)) \cdot 100\% \\ &= (1 - pnorm(1.45)) \cdot 100\% = 7.34\%. \end{aligned} \quad (18) $$

Average age at harvest was equal to 539 days (Table 4). $M$ after genetic change was calculated as:

$$ \begin{aligned} M &= \left(1 - \left(1 - \frac{CM_{t,s,h}}{100}\right)^{\frac{1}{539}}\right) \cdot 100 \\ \left(1 - \left(1 - \frac{7.34}{100}\right)^{\frac{1}{539}}\right) \cdot 100 &= 0.014\%/\text{day}. \end{aligned} \quad (19) $$

The difference in $M$ before and after genetic change was 0.016, which was treated as the $\sigma_A$ of $M$.

Genetic improvement of uniformity reduces the environmental variance of bodyweight. For environmental variance of bodyweight, $CV_A$ was calculated as [38]:

$$ CV_A = \frac{SD(\sigma_E^2)}{\sigma_E^2} \cdot 100\%, \quad (20) $$

where $SD(\sigma_E^2)$ is the genetic standard deviation of environmental variance and $\sigma_E^2$ is the mean environmental variance. Environmental variance equals phenotypic variance minus genetic variance [39]. The $CV_A$ of environmental variance of bodyweight is about 20% in rainbow trout [40] and 41.7% in Atlantic salmon (Salmo salar) [41]. For seabream, the actual value was unknown, hence a minimum of 20% and maximum of 40% were used to represent both extremes of the possible range of $\sigma_A$. In this study, the trait uniformity was expressed on the standard deviation scale instead of the variance scale. On the standard deviation scale, the $CV_A$ is half as large as on the variance scale [42, 43]. For $BW_{t,s,h}$, $\overline{\sigma_E} = \sqrt{\sigma_{HW}^2 - \sigma_A^2} = \sqrt{60^2 - 42.4^2} = 42.45$ g. For the minimum $CV_A$ of the environmental standard deviation of bodyweight of 10%, $\sigma_A$ equals 4.2 g, and for the maximum $CV_A$ of the environmental standard deviation of bodyweight of 20%, $\sigma_A$ equals 8.5 g.

## Validation of the bio-economic model

To validate the bio-economic model, a simplified production system was assumed for which a profit equation can be developed. In this simplified production system, fish were harvested at a constant sum of effective temperatures instead of constant bodyweight. The sum of effective temperatures at harvest was assumed to be unaffected by genetic change. This allowed a profit equation to be set up as a function of the traits: harvest weight ($BW_{t,s,h}$), cumulative feed intake at harvest ($CFI_{t,s,h}$), and survival at harvest ($S_{t,s,h}$). In the bio-economic model, $S_{t,s,h} = \frac{N_{t,s,h}}{N_{t,s,0}} \cdot 100\%$. The bio-economic model was adapted by changing the harvest criterion from a bodyweight of 400 g to a sum of effective temperatures of 4084 day degrees. Thus, an increase in TGC resulted in a greater harvest weight instead of a shorter growing period. One $\sigma_A$ change in TGC led to
change in $BW_{t_{s,h}}$; one $\sigma_A$ change in $TFC$ led to change in $CFI_{t_{s,h}}$; one $\sigma_A$ change in $M$ led to change in $S_{t_{s,h}}$. Economic values were derived from the bio-economic model using Eq. 1 and trait levels of $BW_{t_{s,h}}$, $CFI_{t_{s,h}}$ and $S_{t_{s,h}}$.

In the profit equation, profit at the farm level was described as:

$$ \text{Profit} = \frac{1000 \cdot Q}{BW_{t_{s,h}}} \cdot \left( BW_{t_{s,h}} \cdot \frac{\text{sales price} - \text{packing costs}}{1000} - CFI_{t_{s,h}} \cdot \frac{1}{0.5 + S_{t_{s,h}}/200} \cdot \left( \frac{\text{feed price}}{1000} \right) - \frac{\text{juvenile price}}{S_{t_{s,h}}/100} \right) - \text{fixed costs}, \quad (21) $$

where $Q$ represents production output of the farm (kg) and was calculated as the product of maximum stocking density, cage volume and number of production cycles per year. $Q$ is not affected by genetic change when the harvest criterion is a sum of effective temperatures of 4084 day degrees. Economic values were calculated as partial derivatives of the profit equation and divided by $Q$ to express them per kg fish production. For $BW_{t_{s,h}}$, the economic value was calculated as:

$$ \text{Economic value}_{BW_{t_{s,h}}} = \frac{\delta \text{Profit}}{\delta BW_{t_{s,h}}} \cdot \frac{1}{Q} = \frac{1000}{BW_{t_{s,h}}^2} \cdot \left( \frac{CFI_{t_{s,h}} \cdot (\text{feed price}/1000)}{0.5 + S_{t_{s,h}}/200} + \frac{\text{juvenile price}}{S_{t_{s,h}}/100} \right). \quad (22) $$

For $CFI_{t_{s,h}}$, the economic value was calculated as:

$$ \text{Economic value}_{CFI_{t_{s,h}}} = \frac{\delta \text{Profit}}{\delta CFI_{t_{s,h}}} \cdot \frac{1}{Q} = - \frac{\text{feed price}}{BW_{t_{s,h}} \cdot (0.5 + S_{t_{s,h}}/200)}. \quad (23) $$

For $S_{t_{s,h}}$, the economic value was calculated as:

$$ \text{Economic value}_S = \frac{\delta \text{Profit}}{\delta S} \cdot \frac{1}{Q} = \frac{CFI_{t_{s,h}} \cdot \text{feed price}}{BW_{t_{s,h}} \cdot (50 + S_{t_{s,h}} + S_{t_{s,h}}^2/200)} + \frac{100,000 \cdot \text{juvenile price}}{BW_{t_{s,h}} \cdot S_{t_{s,h}}^2}. \quad (24) $$

# Results

## Production results before genetic change

Tables 4 and 5 show the results of the model for key production variables and costs, respectively. The annual fish production was about 565 tons and gross margin about 759,000€. Average feed costs were €1.18/kg feed and average sales price was €4.49/kg fish.

<table><caption>Table 4 Key production variables of the gilthead seabream farm before genetic change</caption><thead><tr><td>Item FCR</td><td>Value</td></tr></thead><tbody><tr><td>Number of juveniles stocked (year<sup>-1</sup>)</td><td>1,659,945</td></tr><tr><td>Feed consumption (kg/year)</td><td>1,070,177</td></tr><tr><td>Fish production (kg/year)</td><td>564,661</td></tr><tr><td>Cages stocked (year<sup>-1</sup>)</td><td>13.4</td></tr><tr><td>Average age at harvest (day)</td><td>539</td></tr><tr><td>Survival (%)</td><td>85.1</td></tr><tr><td>Biological FCR<sup>a</sup> (kg feed/kg fish)</td><td>1.80</td></tr><tr><td>Economic FCR<sup>b</sup> (kg feed/kg fish)</td><td>1.92</td></tr></tbody></table>

## Production results after genetic change and economic values

The effect of the genetic change on production results is illustrated in Table 6. Changes in gross margin show that the order of economic importance of traits was: TGC, TFC, $M$, and $\sigma_{HW}$. The effect on gross margin of one $\sigma_A$ change for each trait relative to the effect of a 8.5 g decrease in $\sigma_{HW}$ was 43-fold for TGC, 28-fold for TFC, and 12-fold for $M$. Non-linearity was strongest for $\sigma_{HW}$ (results not presented for TGC, TFC, and $M$), for which a doubling of change in trait level from -4.2 to -8.5 g led to 7.7% overestimation of the increase in gross margin.

The mechanisms by which changes in trait levels determined changes in gross margin were as follows. An increase in TGC resulted in a lower age at harvest (Eq. 7) and consequently, the number of production cycles per year increased (Eq. 13) and at the farm level, the annual number of juveniles stocked and annual fish production increased. An increase in TGC did not affect daily feed consumption (Eq. 4) and consequently, cumulative feed intake at harvest decreased because the sum of effective temperatures at harvest decreased (Eq. 5) and at the

<table><caption>Table 5 Economic results for the gilthead seabream farm before genetic change</caption><thead><tr><td>Item</td><td>Farm level (€)</td><td>Fish level (€/kg)</td></tr></thead><tbody><tr><td>Feed costs</td><td>1,259,917</td><td>2.23</td></tr><tr><td>Juvenile costs</td><td>331,989</td><td>0.59</td></tr><tr><td>Packing costs</td><td>186,338</td><td>0.33</td></tr><tr><td>Total variable costs</td><td>1,778,244</td><td>3.15</td></tr><tr><td>Total revenues</td><td>2,537,166</td><td>4.49</td></tr><tr><td>Gross margin</td><td>758,922</td><td>1.34</td></tr></tbody></table>

<footnote><sup>a</sup> Biological FCR = feed consumption/(fish production + biomass mortality - biomass juveniles)</footnote>

<footnote><sup>b</sup> Economic FCR = feed consumption/fish production</footnote>
<table><caption>Table 6 Effect of genetic change on production results relative to the situation without genetic change</caption><thead><tr><th>Trait</th><th>Genetic change (trait unit)</th><th>Δ Juveniles stocked (year<sup>-1</sup>)</th><th>Δ Feed consumption (kg/year)</th><th>Δ Fish production (kg/year)</th><th>Δ Gross margin (€)</th></tr></thead><tbody><tr><td>TGC<sup>a</sup></td><td>+0.95</td><td>89,400</td><td>-68,409</td><td>36,309</td><td>213,131</td></tr><tr><td>TFC<sup>b</sup></td><td>-0.55</td><td>0</td><td>-120,300</td><td>0</td><td>140,891</td></tr><tr><td>M<sup>c</sup></td><td>-0.016</td><td>-137,294</td><td>-34,446</td><td>0</td><td>69,531</td></tr><tr><td>σ<sub>HW</sub><sup>d</sup></td><td>-4.2</td><td>0</td><td>0</td><td>0</td><td>2636</td></tr><tr><td>σ<sub>HW</sub></td><td>-8.5</td><td>0</td><td>0</td><td>0</td><td>4952</td></tr></tbody></table>

farm level, the annual feed consumption decreased. A decrease in TFC decreased total feed consumption per production cycle (Eqs. 5, 11) but the number of juveniles stocked per production cycle and the fish production per production cycle remained unaltered, thus at the farm level, only annual feed consumption decreased. A decrease in M reduced the number of juveniles stocked per production cycle (Eq. 9) and consequently, daily feed intake per cage decreased because the average number of fish per cage per day was smaller (Eq. 10), but fish production per production cycle was unaltered. Thus at the farm level, the annual number of juveniles stocked and annual feed consumption decreased. The effect of σ<sub>HW</sub> on the average sales price is illustrated in Fig. 2: less variation led to more sales in size category 300 to 400 g (€4.52/kg) at the expense of sales in size category 200 to 300 g (€4.15/kg). Production results were unaltered by a change in σ<sub>HW</sub>.

Economic values are in Table 7, which shows that the economic value of σ<sub>HW</sub> was similar for both levels of genetic change.

<figure><img src="image_5.png" /><figcaption>Fig. 2 Distribution of harvest weight over sales price categories at different standard deviations of harvest weight (σ<sub>HW</sub>)</figcaption></figure>

<table><caption>Table 7 Economic values of traits for gilthead seabream</caption><thead><tr><th>Trait</th><th>Baseline trait level (trait unit)</th><th>Genetic change (trait unit)</th><th>Economic value [€ (kg production)<sup>-1</sup> (trait unit)<sup>-1</sup>]</th></tr></thead><tbody><tr><td>TGC<sup>a</sup></td><td>12.6</td><td>+0.95</td><td>0.40</td></tr><tr><td>TFC<sup>b</sup></td><td>8.25</td><td>-0.55</td><td>-0.45</td></tr><tr><td>M<sup>c</sup></td><td>0.0300</td><td>-0.016</td><td>-7.7</td></tr><tr><td>σ<sub>HW</sub><sup>d</sup></td><td>60</td><td>-4.2</td><td>-0.0011</td></tr><tr><td>σ<sub>HW</sub></td><td>60</td><td>-8.5</td><td>-0.0010</td></tr></tbody></table>

## Comparison of economic values from the bio-economic model and the profit equation

Table 8 shows that, for the simplified production system, the economic values derived from the bio-economic model and the profit equation were similar.

## Discussion

### Validity of the bio-economic model

For the simplified production system, the profit equation and the bio-economic model return similar economic values, which confirm the validity of the bio-economic model. To further validate the bio-economic model, production results were compared to those of other studies. FCR (Table 4) was within the range of 1.5 to 2 reported by Sola et al. [44] but considerably lower than the 2.3 value reported by EAS-EATiP [45]. Overall survival was 85%, which is within the range reported by EAS-EATiP [45]. A comparison with a cost-breakdown for large-scale production of gilthead seabream and European seabass (Dicentrarchus labrax) is in Table 9 [17]. In the FAO data, variable costs are higher, largely because labor, energy, and medicines and veterinary services were not considered to be variable costs in our study. Trends in the increase in productivity per person [46] support the

<footnote>a Thermal growth coefficient [g<sup>2/3</sup>/(day degrees · 1000)]</footnote>

<footnote>b Thermal feed intake coefficient [g<sup>0.544</sup>/(day degrees · 1000)]</footnote>

<footnote>c Mortality rate (%/day)</footnote>

<footnote>d Standard deviation of harvest weight (g)</footnote>

<footnote>a Thermal growth coefficient [g<sup>2/3</sup>/(day degrees · 1000)]</footnote>

<footnote>b Thermal feed intake coefficient [g<sup>0.544</sup>/(day degrees · 1000)]</footnote>

<footnote>c Mortality rate (%/day)</footnote>

<footnote>d Standard deviation of harvest weight (g)</footnote>
<table><caption>**Table 9** Cost-breakdown for gilthead seabream production</caption><thead><tr><th rowspan="2">Item</th><th colspan="2">Proportion of total costs (%)</th></tr><tr><th>Our study<sup>a</sup></th><th>Barazi-Yeroulanos [17]</th></tr></thead><tbody><tr><td>Feed</td><td>50</td><td>48</td></tr><tr><td>Juveniles</td><td>13</td><td>11</td></tr><tr><td>Marketing (incl. packing)</td><td>7</td><td>18</td></tr><tr><td>Labor</td><td>-</td><td>3</td></tr><tr><td>Energy</td><td>-</td><td>4</td></tr><tr><td>Medicines and veterinary services</td><td>-</td><td>2</td></tr><tr><td>Other</td><td>-</td><td>4</td></tr><tr><td>Total variable costs</td><td>70</td><td>89</td></tr><tr><td>Total fixed costs</td><td>-</td><td>11</td></tr></tbody></table>

assumption that labor should be treated more as a fixed than a variable cost. Medicine costs may vary, but veterinary costs are likely to be fixed per farm. Energy costs are to a larger extent determined by farm layout than by realized production and thus can be considered as fixed. Altogether, total variable costs may have been slightly underestimated in our study, but FCR and overall survival matched well to current industry standards.

## Breeding goal

In the breeding goal, TGC, TFC, and M are equivalent to respectively $BW_{t_{s,h}}$, $CFI_{t_{s,h}}$, and $S_{t_{s,h}}$ when $BW_{t_{s,0}}$ is much smaller than $BW_{t_{s,h}}$ (Appendixes 1, 2). When the sum of effective temperatures is the harvest criterion, one $\sigma_A$ change in TGC, TFC, and M led to changes in $BW_{t_{s,h}}$, $CFI_{t_{s,h}}$, and $S_{t_{s,h}}$ that were very similar to the $\sigma_A$ of these traits (Table 8), which demonstrates their equivalence. If the economic values of TGC, TFC, and M were calculated for the sum of effective temperatures instead of harvest weight as the harvest criterion, they would be slightly lower for TGC [0.34€ (kg production)$^{-1}$ ($g^{0.544}/(\text{day degrees} \cdot 1000)$)$^{-1}$] and unaltered for TFC and M. In agreement with Wilton and Goddard [47], economic values were similar for both harvest criteria. Although

both sets of traits are equivalent in the breeding goal, there are pros and cons to each one. $BW_{t_{s,h}}$ is commonly used as a selection criterion and thus its use in the breeding goal is straightforward. However, $BW_{t_{s,h}}$ is a management parameter that is strongly influenced by the growing period and temperature regime. TGC corrects for heterogeneity in stocking weight, growing period, and temperature regime, and, therefore, allows for a better comparison of breeding values across conditions than $BW_{t_{s,h}}$ [11, 24, 25].

FCR could be used as an alternative to TFC in the breeding goal. An advantage of feed intake compared to FCR is that it relates directly to feed costs [48]. An advantage of FCR is that it illustrates the effect of improvement in efficiency on, for example, environmental impacts, as in Besson et al. [49]. Feed intake is often considered more appropriate as a breeding goal trait than FCR [18, 19], with a common argument that traits expressed as ratio's are disadvantageous in animal breeding [19]. Selection for a ratio, e.g. FCR, results in a lower selection response than selection for both components of the ratio, e.g. feed intake and growth [50]. However, in the same way that FCR is a ratio, growth is the ratio of feed intake to FCR, thus a breeding goal that includes both growth and FCR is equivalent to a breeding goal that includes growth and feed intake. The economic value of growth depends on which other trait, feed intake or FCR, is included in the breeding goal [48].

## Economic values

Bio-economic models and profit equations are both suitable to derive economic values. An advantage of a profit equation compared to a bio-economic model is its simplicity. However, its applicability is limited to specific situations, because environmental conditions are ignored. For example, the profit equation cannot be used to derive economic values for a range of temperature regimes, as was done in Besson et al. [11] by using the bio-economic model. Such properties may be of particular interest for breeding programs that aim at supplying many farms. In addition, alternative constraints on production output

<table><caption>**Table 8** Economic values derived from the bio-economic model and profit equation</caption><thead><tr><th rowspan="2">Trait</th><th rowspan="2">Baseline trait level (trait unit)</th><th rowspan="2">Genetic change (trait unit)</th><th colspan="2">Economic value [€ (kg production)<sup>-1</sup> (trait unit)<sup>-1</sup>]</th></tr><tr><th>Bio-economic model</th><th>Profit equation</th></tr></thead><tbody><tr><td><i>BW</i><sub>t<sub>s,h</sub></sub><sup>a</sup></td><td>400</td><td>43.6</td><td>0.0074</td><td>0.0072</td></tr><tr><td><i>CFI</i><sub>t<sub>s,h</sub></sub><sup>b</sup></td><td>713</td><td>-80.0</td><td>-0.0031</td><td>-0.0032</td></tr><tr><td><i>S</i><sub>t<sub>s,h</sub></sub><sup>c</sup></td><td>85.1</td><td>7.66</td><td>0.016</td><td>0.019</td></tr></tbody></table>

both sets of traits are equivalent in the breeding goal, there are pros and cons to each one. $BW_{t_{s,h}}$ is commonly used as a selection criterion and thus its use in the breeding goal is straightforward. However, $BW_{t_{s,h}}$ is a management parameter that is strongly influenced by the growing period and temperature regime. TGC corrects for heterogeneity in stocking weight, growing period, and temperature regime, and, therefore, allows for a better comparison of breeding values across conditions than $BW_{t_{s,h}}$ [11, 24, 25].

FCR could be used as an alternative to TFC in the breeding goal. An advantage of feed intake compared to FCR is that it relates directly to feed costs [48]. An advantage of FCR is that it illustrates the effect of improvement in efficiency on, for example, environmental impacts, as in Besson et al. [49]. Feed intake is often considered more appropriate as a breeding goal trait than FCR [18, 19], with a common argument that traits expressed as ratio's are disadvantageous in animal breeding [19]. Selection for a ratio, e.g. FCR, results in a lower selection response than selection for both components of the ratio, e.g. feed intake and growth [50]. However, in the same way that FCR is a ratio, growth is the ratio of feed intake to FCR, thus a breeding goal that includes both growth and FCR is equivalent to a breeding goal that includes growth and feed intake. The economic value of growth depends on which other trait, feed intake or FCR, is included in the breeding goal [48].

## Economic values

Bio-economic models and profit equations are both suitable to derive economic values. An advantage of a profit equation compared to a bio-economic model is its simplicity. However, its applicability is limited to specific situations, because environmental conditions are ignored. For example, the profit equation cannot be used to derive economic values for a range of temperature regimes, as was done in Besson et al. [11] by using the bio-economic model. Such properties may be of particular interest for breeding programs that aim at supplying many farms. In addition, alternative constraints on production output

<footnote><sup>a</sup> Relative to revenues</footnote>

<footnote><sup>a</sup> Harvest weight (g)</footnote>

<footnote><sup>b</sup> Cumulative feed intake at harvest (g)</footnote>

<footnote><sup>c</sup> Survival at harvest (%)</footnote>
such as oxygen availability cannot be dealt with by the profit equation but can be incorporated in the bio-economic model, as discussed later. Furthermore, the profit equation is rigid in terms of trait definition, which has led to the false assumption that harvest weight changes following genetic improvement, whereas in the bio-economic model genetic improvement of growth rate leads to a reduction in the growing period.

From a profit function, economic values can be computed from either its partial derivative with respect to trait level, or from an increase or decrease in trait level relative to the current mean. In this study, simulated changes in trait levels correspond to desired directions of genetic change of one genetic standard deviation. However, for a non-linear profit function, Goddard [51] demonstrated that economic values that maximize profit in the next generation may depend on selection responses. Dekkers et al. [52] showed that economic gain is slightly higher when economic values are derived as the partial derivative of a non-linear profit equation at the genetic level of the next generation than at the genetic level of the current generation. This implies two things:

(1) Economic values are closer to optimum when the simulated change in trait level resembles its expected rather than its desired direction.

(2) Economic values are closer to optimum when the simulated change in trait level equals the difference between trait levels in the current and next generation than when it is the partial derivative at current genetic levels.

In our study, expected and desired directions of change in trait levels were identical, except for *TFC*, which may increase in practice due to its genetic correlation with *TGC* [34]. A genetic standard deviation generally provides a better proxy for the difference between trait levels in the current and next generation than an infinitesimal change, and hence will result in an economic value that is closer to optimum than the conventional partial derivative at the current trait level.

This is the first time that economic values have been derived for uniformity in aquaculture, here expressed as $σ_{HW}$. In recent years, there has been increasing interest to improve uniformity [40, 53–55]. Improvement in uniformity affects average sales price and reduces the need for size-grading. However, for seabream production, reducing the need for size-grading would not result in major cost-savings, because seabream is size-graded only once during grow-out. Thus, a potential effect on grading frequency was excluded from the economic value of uniformity. Furthermore, uniformity has been suggested to affect feed intake, growth rate and mortality [56, 57].

Economic consequences of changes in other traits were accounted for in their respective economic values. By exploiting genetic correlations, selection for uniformity may be used to improve the other traits in the breeding goal.

## Application to other aquaculture species

In its current form, the bio-economic model can be easily applied for the derivation of economic values for other species produced in systems where stocking density limits its production output, such as cages and flow-through tanks. This would require different values for the coefficients of Table 1, maximum stocking density, stocking and harvest weight, and input and output prices. Equations 4 and 5 require some species-specific modifications, such as alternative values for exponents 2/3 and 3/2 in Eq. 4 [58] or a different minimum temperature for growth.

Adaptations to the model are required for species that are reared in production systems for which constraints on production output are different, such as recirculating aquaculture systems and ponds. When the constraint on production output is different from stocking density, the number of fish stocked per production cycle (Eq. 9) is determined by other parameters. For recirculating aquaculture systems, treatment capacity of the biofilter can be a constraint on production output [12]. In this case, daily nitrogen excretion by fish is the parameter that determines the number of fish stocked per production cycle. Daily nitrogen excretion by fish can be predicted from the difference between daily feed consumption and daily gain in bodyweight, as described in Besson et al. [12]. In both cages and ponds, oxygen availability can be a constraint on production output. In this case, daily oxygen consumption per fish is the parameter that determines the number of fish stocked per production cycle. Daily oxygen consumption per fish can be predicted from daily feed consumption and daily gain in bodyweight, as described in Besson et al. [11]. With the above modifications, the same bio-economic model was applied for the derivation of economic values for African catfish produced in recirculating aquaculture systems [12], European seabass produced in cages [11], gilthead seabream produced in cages (this study), turbot produced in tanks (unpublished results), and Nile tilapia produced in ponds (unpublished results).

## Conclusions

We developed a bio-economic model to derive economic values for a wide range of aquaculture species. Its validity was confirmed by the comparison to a profit equation for a simplified production system and by comparison of the production results to those of other studies. Application
of the bio-economic model to gilthead seabream resulted in economic values for TGC, TFC, M, and σHW. TGC was the most important trait to improve, followed by TFC and M. The effect of σHW on gross margin was small.

## Authors' contributions

KJ developed the bio-economic model, performed the analysis, and wrote the manuscript. HK and PB contributed to discussions and writing of the manuscript. MB helped in developing the bio-economic model and contributed to discussions. All authors read and approved the final manuscript.

## Author details

¹ Animal Breeding and Genomics, Wageningen University and Research, Droevendaalsesteeg 1, 6708 PB Wageningen, The Netherlands. ² Business Economics Group, Wageningen University and Research, Hollandseweg 1, 6706 KN Wageningen, The Netherlands. ³ Génétique Animale Biologie Intérative, INRA, AgroParisTech, Université Paris-Saclay, 78350 Jouy-en-Josas, France.

## Acknowledgements

We are grateful for the extensive production data provided by Andromeda S.A. and their feedback on model assumptions. We thank Piter Bijma for his help in the computation of genetic standard deviations.

## Competing interests

The authors declare that they have no competing interests.

## Availability of data and material

Production data of Andromeda S.A. will not be shared, because confidentiality was agreed.

## Funding

The research leading to these results has received funding from the European Union's Seventh Framework Programme (KBBE.2013.1.2-10) under grant agreement no. 613611.

# Appendices

## Appendix 1: Genetic variation in TGC

Genetic variation in TGC depends on genetic variation both in BW<sub>t<sub>s,0</sub></sub> and BW<sub>t<sub>s,h</sub></sub>:

$$Var(A_{TGC}) = Var \left( 1000 \cdot \frac{A_{BW_{t_{s,h}}}^{2/3} - A_{BW_{t_{s,0}}}^{2/3}}{\sum_{i=t_{s,0}}^{(t_{s,h})-1} (T_i - 12)} \right), \quad (25)$$

where A<sub>TGC</sub> is the genotype for TGC, A<sub>BW<sub>t<sub>s,h</sub></sub></sub> is the genotype for harvest weight, and A<sub>BW<sub>t<sub>s,0</sub></sub></sub> is the genotype for stocking weight. Equation 25 can be rewritten as:

$$Var(A_{TGC}) = \left( \frac{1000}{\sum_{i=t_{s,0}}^{(t_{s,h})-1} (T_i - 12)} \right)^2 \cdot \left( Var(A_{BW_{t_{s,h}}}^{2/3}) + Var(A_{BW_{t_{s,0}}}^{2/3}) - 2 \cdot cov(A_{BW_{t_{s,h}}}^{2/3}, A_{BW_{t_{s,0}}}^{2/3}) \right). \quad (26)$$

Because BW<sub>t<sub>s,0</sub></sub> is much smaller than BW<sub>t<sub>s,h</sub></sub>, Var(A<sub>BW<sub>t<sub>s,0</sub></sub></sub><sup>2/3</sup>) - 2 · cov(A<sub>BW<sub>t<sub>s,h</sub></sub></sub><sup>2/3</sup>, A<sub>BW<sub>t<sub>s,0</sub></sub></sub><sup>2/3</sup>) is much smaller than Var(A<sub>BW<sub>t<sub>s,h</sub></sub></sub><sup>2/3</sup>) [59]. Thus, Eq. 26 can be reduced to Eq. 15.

## Appendix 2: Genetic variation in cumulative feed intake at harvest

The regression coefficient of the genotype for CFI<sub>t<sub>s,h</sub></sub> (A<sub>CFI<sub>t<sub>s,h</sub></sub></sub>) on the difference between A<sub>BW<sub>t<sub>s,h</sub></sub></sub> and A<sub>BW<sub>t<sub>s,0</sub></sub></sub> can be calculated as:

$$b(A_{CFI_{t_{s,h}}}, A_{BW_{t_{s,h}}} - A_{BW_{t_{s,0}}}) = r_A \cdot \sqrt{\frac{Var(A_{CFI_{t_{s,h}}})}{Var(A_{BW_{t_{s,h}}} - A_{BW_{t_{s,0}}})}}, \quad (27)$$

where b(A<sub>CFI<sub>t<sub>s,h</sub></sub></sub>, A<sub>BW<sub>t<sub>s,h</sub></sub></sub> - A<sub>BW<sub>t<sub>s,0</sub></sub></sub>) is the regression coefficient, and r<sub>A</sub> is the genetic correlation coefficient. For the regression of A<sub>CFI<sub>t<sub>s,h</sub></sub></sub> on A<sub>BW<sub>t<sub>s,h</sub></sub></sub> - A<sub>BW<sub>t<sub>s,0</sub></sub></sub>, the intercept corresponds to feed consumption to meet maintenance energy requirements (CFI<sub>maintenance</sub>) at zero growth. Assuming a digestible energy content of the diet of 17 kJ/g, CFI<sub>maintenance</sub> is calculated as [26]:

$$CFI_{maintenance} = \sum_{i=t_{s,0}}^{t_{s,h}} 47.89 \cdot (BW_i/1000)^{0.80}/17 = 20 \text{ g.} \quad (28)$$

For A<sub>BW<sub>t<sub>s,h</sub></sub></sub> - A<sub>BW<sub>t<sub>s,0</sub></sub></sub> = 400 - 4.4 = 395.6 g, A<sub>CFI<sub>t<sub>s,h</sub></sub></sub> equals 1.80 · 395.6 = 713 g, where 1.80 is the biological FCR (Table 4). The regression coefficient can thus be approximated as (713 - 20)/395.6 = 1.75 and Eq. 27 can be rewritten as:

$$1.75 = r_A \cdot \sqrt{\frac{Var(A_{CFI_{t_{s,h}}})}{Var(A_{BW_{t_{s,h}}}) + Var(A_{BW_{t_{s,0}}}) - 2 \cdot cov(A_{BW_{t_{s,h}}}, A_{BW_{t_{s,0}}})}}. \quad (29)$$

Because BW<sub>t<sub>s,0</sub></sub> is much smaller than BW<sub>t<sub>s,h</sub></sub>, Var(A<sub>BW<sub>t<sub>s,0</sub></sub></sub>) - 2 · cov(A<sub>BW<sub>t<sub>s,h</sub></sub></sub>, A<sub>BW<sub>t<sub>s,0</sub></sub></sub>) is much smaller than Var(A<sub>BW<sub>t<sub>s,h</sub></sub></sub>) [59]. Thus Eq. 29 can be reduced to Eq. 16.

Received: 30 May 2016 Accepted: 5 December 2016

Published online: 05 January 2017

## References

1. Janssen K, Chavanne H, Berentsen P, Komen H. Impact of selective breeding on European aquaculture. Aquaculture. 2016.; doi:10.1016/j. aquaculture.2016.03.012.

2. Nielsen HM, Amer PR, Byrne TJ. Approaches to formulating practical breeding objectives for animal production systems. Acta Agric Scand A Anim. 2014;64:2–12.
3. Shook GE. Major advances in determining appropriate selection goals. *J Dairy Sci.* 2006;89:1349–61.

4. Gibson JP, Kennedy BW. The use of constrained selection indexes in breeding for economic merit. *Theor Appl Genet.* 1990;80:801–5.

5. Gjedrem T, Baranski M. Selective breeding in aquaculture: an introduction. *Dordrecht: Springer;* 2009.

6. Gjedrem T, Thodesen J. Selection. In: Gjedrem T, editor. Selection and breeding programs in aquaculture. *Dordrecht: Springer;* 2005. p. 89–111.

7. Ponzoni RW, Nguyen NH, Khaw HL. Investment appraisal of genetic improvement programs in Nile tilapia (*Oreochromis niloticus*). *Aquaculture.* 2007;269:187–99.

8. Ponzoni RW, Nguyen NH, Khaw HL, Ninh NH. Accounting for genotype by environment interaction in economic appraisal of genetic improvement programs in common carp *Cyprinus carpio*. *Aquaculture.* 2008;285:47–55.

9. Zuniga-Jara S, Marin-Riffo MC. A bioeconomic model of a genetic improvement program of abalone. *Aquacult Int.* 2014;22:1533–62.

10. Henryon M, Purvis IW, Berg P. Definition of a breeding objective for commercial production of the freshwater crayfish, marron (*Cherax tenuimanus*). *Aquaculture.* 1999;173:179–95.

11. Besson M, Komen H, Aubin J, de Boer IJM, Poelman M, Quillet E, et al. Economic values of growth and feed efficiency for fish farming in recirculating aquaculture system with density and nitrogen output limitations: a case study with African catfish (*Clarias gariepinus*). *J Anim Sci.* 2014;92:5394–405.

12. Besson M, Vandeputte M, van Arendonk JAM, Aubin J, de Boer IJM, Quillet E, et al. Influence of water temperature on the economic value of growth rate in fish farming: the case of sea bass (*Dicentrarchus labrax*) cage farming in the Mediterranean. *Aquaculture.* 2016. doi:10.1016/j.aquaculture.2016.04.030.

13. Knap PW. Breeding robust pigs. *Austr J Exp Agric.* 2005;45:763–73.

14. Hietala P, Wolfova M, Wolf J, Kantanen J, Juga J. Economic values of production and functional traits, including residual feed intake, in Finnish milk production. *J Dairy Sci.* 2014;97:1092–106.

15. van Middelaar CE, Berentsen PBM, Dijkstra J, van Arendonk JAM, de Boer IJM. Methods to determine the relative value of genetic traits in dairy cows to reduce greenhouse gas emissions along the chain. *J Dairy Sci.* 2014;97:5191–205.

16. Byrne TJ, Amer PR, Fennessy PF, Cromie AR, Keady TWJ, Hanrahan JP, et al. Breeding objectives for sheep in Ireland: a bio-economic approach. *Livest Sci.* 2010;132:135–44.

17. Barazi-Yeroulanos L. Synthesis of Mediterranean marine finfish aquaculture—a marketing and promotion strategy. In: *Studies and reviews. General Fisheries Commission for the Mediterranean*. no. 88. Rome: FAO; 2010. p. 1–198.

18. Emmerson DA. Commercial approaches to genetic selection for growth and feed conversion in domestic poultry. *Poult Sci.* 1997;76:1121–5.

19. Veerkamp RF, Pryce JE, Spurlock D, Berry D, Coffey M, Lovendahl P, et al. Selection on feed intake or feed efficiency: a position paper from gDMI breeding goal discussions. *Interbull Bull.* 2013;47:15–22.

20. Wolfová M, Wolf J. Strategies for defining traits when calculating economic values for livestock breeding: a review. *Animal*. 2013;7:1401–13.

21. Iwama GK, Tautz AF. A simple growth model for Salmonids in hatcheries. *Can J Fish Aquat Sci.* 1981;38:649–56.

22. Jobling M. The thermal growth coefficient (TGC) model of fish growth: a cautionary note. *Aquacult Res.* 2003;34:581–4.

23. Cho CY. Feeding systems for rainbow trout and other salmonids with reference to current estimates of energy and protein requirements. *Aquaculture.* 1992;100:107–23.

24. Sae-Lim P, Kause A, Mulder HA, Martin KE, Barfoot AJ, Parsons JE, et al. Genotype-by-environment interaction of growth traits in rainbow trout (*Oncorhynchus mykiss*): a continental scale study. *J Anim Sci.* 2013;91:5572–81.

25. Trong TQ, Mulder HA, van Arendonk JAM, 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. *Aquaculture.* 2013;384–387:119–27.

26. Lupatsch I, Kissil GW, Sklan D. Comparison of energy and protein efficiency among three fish species gilthead sea bream (*Sparus aurata*), European sea bass (*Dicentrarchus labrax*) and white grouper (*Epinephelus aeneus*): energy expenditure for protein and lipid deposition. *Aquaculture.* 2003;225:175–89.

27. R Core Team. *R: A Language and environment for statistical computing*. 2015. http://www.R-project.org/.

28. Houle D. Comparing evolvability and variability of quantitative traits. *Genetics*. 1992;130:195–204.

29. Groen AF. Economic values in cattle-breeding. 2. Influences of production circumstances in situations with output limitations. *Livest Prod Sci.* 1989;22:17–30.

30. Cacho OJ. Protein and fat dynamics in fish: a bioenergetic model applied to aquaculture. *Ecol Model.* 1990;50:33–56.

31. Mayer P, Estruch V, Blasco J, Jover M. Predicting the growth of gilthead sea bream (*Sparus aurata* L.) farmed in marine cages under real production conditions using temperature- and time-dependent models. *Aquacult Res.* 2008;39:1046–52.

32. Mayer P, Estruch VD, Jover M. A two-stage growth model for gilthead sea bream (*Sparus aurata*) based on the thermal growth coefficient. *Aquaculture.* 2012;358–359:6–13.

33. Navarro A, Zamorano MJ, Hildebrandt S, Ginés R, Aguilera C, Afonso JM. Estimates of heritabilities and genetic correlations for growth and carcass traits in gilthead seabream (*Sparus auratus* L.), under industrial conditions. *Aquaculture.* 2009;289:225–30.

34. Quinton CD, Kause A, Koskela J, Ritola O. Breeding salmonids for feed efficiency in current fishmeal and future plant-based diet environments. *Genet Sel Evol.* 2007;39:431.

35. Kause A, Tobin D, Dobly A, Houlihan D, Martin S, Mantysaari EA, et al. Recording strategies and selection potential of feed intake measured using the X-ray method in rainbow trout. *Genet Sel Evol.* 2006;38:389–409.

36. Falconer DS, Mackay TFC. Introduction to quantitative genetics. 3rd ed. Harlow: Longman; 1989.

37. Vehvilainen H, Kause A, Quinton C, Koskinen H, Paananan T. Survival of the currently fittest: genetics of rainbow trout survival across time and space. *Genetics*. 2008;180:507–16.

38. Mulder HA, Bijma P, Hill WG. Prediction of breeding values and selection responses with genetic heterogeneity of environmental variance. *Genetics.* 2007;175:1895–910.

39. Hazel LN. The genetic basis for constructing selection indexes. *Genetics.* 1943;28:476–90.

40. Sae-Lim P, Kause A, Janhunen M, Vehvilainen H, Koskinen H, Gjerde B, et al. Genetic (co)variance of rainbow trout (*Oncorhynchus mykiss*) body weight and its uniformity across production environments. *Genet Sel Evol.* 2015;47:46.

41. Sonesson AK, Odegard J, Ronnegard L. Genetic heterogeneity of within-family variance of body weight in Atlantic salmon (*Salmo salar*). *Genet Sel Evol.* 2013;45:41.

42. Hill WG, Mulder HA. Genetic analysis of environmental variation. *Genet Res (Camb).* 2010;92:381–95.

43. Sell-Kubiak E, Bijma P, Knol EF, Mulder HA. Comparison of methods to study uniformity of traits: application to birth weight in pigs. *J Anim Sci.* 2015;93:900–11.

44. Sola L, Moretti A, Crosetti D, Karaiskou N, Magoulas A, Rossi AR, et al. Gilthead seabream—*Sparus aurata*. In: D Crossetti, S Lapègue, I Olesen, T Svaasand, editors. Genetic effects of domestication, culture and breeding of fish and shellfish, and their impacts on wild populations. *Genimpact Final Scientific Report. Viterbo;* 2007. p. 47–54.

45. EAS-EATIP. Performance of the sea bass and sea bream sector in the Mediterranean. In Minutes of a Workshop held within Aquaculture Europe: 16 October 2014; San Sebastian; 2014.

46. University of Stirling. Study of the market for aquaculture produced seabass and seabream species; 2004. http://ec.europa.eu/fisheries/documentation/studies/aquaculture_market_230404_en.pdf. Accessed 12 April 2016.

47. Wilton JW, Goddard ME. Selection for carcass and feedlot traits considering alternative slaughter end points and optimized management. *J Anim Sci.* 1996;74:37–45.

48. Goddard ME. Consensus and debate in the definition of breeding objectives. *J Dairy Sci.* 1998;81:6–18.

49. Besson M, Aubin J, Komen H, Poelman M, Quillet E, Vandeputte M, et al. Environmental impacts of genetic improvement of growth rate and feed conversion ratio in fish farming under rearing density and nitrogen output limitations. *J Clean Prod.* 2016;116:100–9.

50. Gunsett FC. Linear index selection to improve traits defined as ratios. *J Anim Sci.* 1984;59:1185–93.
51. Goddard ME. Selection indexes for non-linear profit-functions. Theor Appl Genet. 1983;64:339-44.

52. Dekkers JCM, Birke PV, Gibson JP. Optimum linear selection indexes for multiple generation objectives with nonlinear profit-functions. Anim Sci. 1995;61:165-75.

53. Mulder HA, Bijma P, Hill WG. Selection for uniformity in livestock by exploiting genetic heterogeneity of residual variance. Genet Sel Evol. 2008;40:37.

54. Khaw HL, Ponzoni RW, Yee HY, Aziz MA, Mulder HA, Marjanovic J, et al. Genetic variance for uniformity of harvest weight in Nile tilapia (Oreochromis niloticus). Aquaculture. 2016;451:113-20.

55. Marjanovic J, Mulder HA, Khaw HL, Bijma P. Genetic parameters for uniformity of harvest weight in the GIFT strain of Nile tilapia estimated using double hierarchical generalized linear models. Genet Sel Evol. 2016;48:41.

56. Jobling M. Simple indices for the assessment of the influences of social environment on growth performance, exemplified by studies on Arctic charr. Aquacult Int. 1995;3:60-5.

57. Gilmour KM, DiBattista JD, Thomas JB. Physiological causes and consequences of social status in salmonid fish. Integr Comp Biol. 2005;45:263-73.

58. Dumas A, France J, Bureau DP. Evidence of three growth stanzas in rain-bow trout (*Oncorhynchus mykiss*) across life stages and adaptation of the thermal-unit growth coefficient. Aquaculture. 2007;267:139-46.

59. Rutten MJM, Komen H, Bovenhuis H. Longitudinal genetic analysis of Nile tilapia (*Oreochromis niloticus* L.) body weight using a random regression model. Aquaculture. 2005;246:101-13.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central