Skip to main content
  • Research Article
  • Open access
  • Published:

Multivariate random regression analysis for body weight and main morphological traits in genetically improved farmed tilapia (Oreochromis niloticus)

Abstract

Background

Because of their high economic importance, growth traits in fish are under continuous improvement. For growth traits that are recorded at multiple time-points in life, the use of univariate and multivariate animal models is limited because of the variable and irregular timing of these measures. Thus, the univariate random regression model (RRM) was introduced for the genetic analysis of dynamic growth traits in fish breeding.

Methods

We used a multivariate random regression model (MRRM) to analyze genetic changes in growth traits recorded at multiple time-point of genetically-improved farmed tilapia. Legendre polynomials of different orders were applied to characterize the influences of fixed and random effects on growth trajectories. The final MRRM was determined by optimizing the univariate RRM for the analyzed traits separately via penalizing adaptively the likelihood statistical criterion, which is superior to both the Akaike information criterion and the Bayesian information criterion.

Conclusions

In the selected MRRM, the additive genetic effects were modeled by Legendre polynomials of three orders for body weight (BWE) and body length (BL) and of two orders for body depth (BD). By using the covariance functions of the MRRM, estimated heritabilities were between 0.086 and 0.628 for BWE, 0.155 and 0.556 for BL, and 0.056 and 0.607 for BD. Only heritabilities for BD measured from 60 to 140 days of age were consistently higher than those estimated by the univariate RRM. All genetic correlations between growth time-points exceeded 0.5 for either single or pairwise time-points. Moreover, correlations between early and late growth time-points were lower. Thus, for phenotypes that are measured repeatedly in aquaculture, an MRRM can enhance the efficiency of the comprehensive selection for BWE and the main morphological traits.

Background

From an economical point of view, growth and developmental characters are the most important traits in farmed fish species, and persistent efforts are made to genetically improve these traits in fish breeding. Growth and developmental traits, such as body weight (BWE) and morphological traits, are measured at different times. Faster BWE growth shortens time to market and selection on morphological traits allows fish shape and size to be standardized. Growth and developmental traits are dynamic quantitative traits because they vary spatially and temporally [1]. They are also called infinite-dimensional traits, since they are expressed on a continuous time and space scale [2]. In tilapia breeding, growth and developmental traits at a specific age (in day, week, or year), such as the time to market are the main targets of genetic improvement. As such, genetic analysis of the growth and developmental process can increase the efficiency of selection compared to that of traits measured at specific ages [3].

Considerable attention has been paid to the genetic analysis of BWE and morphological traits at a specific age in breeding tilapia. In earlier studies, growth traits at specific ages were considered as separate traits and analyzed by using univariate animal models [4,5,6,7,8,9,10,11,12,13]. Subsequently, multivariate animal models were applied to estimate both the heritabilities of growth traits at specific ages and the genetic correlations between traits measured at different specific ages [8, 11, 12, 14,15,16,17,18,19,20,21,22]. Although multivariate genetic analysis of growth traits has improved the accuracy of parameter estimation by using more records, at multiple ages, such an analysis is strongly limited by the variable and irregular timing of these measurements [23].

For growth and developmental traits that are recorded at multiple time-points during growth, the patterns of growth curves have been shown to be heritable [2, 3, 24]. The genetic analysis of dynamic quantitative traits was initially conducted by first fitting individual growth curves and then analyzing their estimated parameters within the framework of a multivariate animal model. However, with such genetic analyses, it was not possible to determine whether some individuals had too few records to fit the growth curves. Since then, random regression models (RRM) [3, 25] were developed to dissect dynamic phenotypes into different functions that describe the effects of various genetic and environmental factors. Only three papers on the application of RRM for the genetic analysis of growth curves have been published, in rainbow trout [26] and tilapia [9, 27]. The use of RRM allows heritabilities of growth traits at any age to be estimated, as well as genetic correlations between pairwise traits measured at different ages. Moreover, multiple trait RRM (MRRM) can estimate genetic correlations between pairs of traits measured at different ages and improve heritability estimates for each trait. This not only facilitates the genetic analysis of dynamic traits, but also improves the prediction of breeding values [3, 28].

The objective of this study was to construct and implement a MRRM to simultaneously model the genetic changes in growth traits and estimate the heritability of each growth trait, as well as the genetic correlations between pairs of growth traits measured at specific ages, such as the time to market for a population of genetically improved farmed tilapia in order to formulate a selection criterion for simultaneously selecting BWE and morphological traits. For this purpose, BWE, body length (BL), and body depth (BD) were measured six times on 1451 fish from 45 mixed families of full and half-sibs.

Methods

Experimental population

The base population of 1800 one-month old tilapia fingerlings (sex ratio of 1:1) was imported from WorldFish, Malaysia, in 2006, and was derived from 60 families with complete pedigree. The fish were systematically selected for seven generations at the experimental station of the Freshwater Fisheries Research Center in Wuxi, China. Based on breeding values for BWE at 120 days of age estimated using an animal model, 40 broodstocks per family were chosen at each generation, with a sex ratio of 1:1. The selection intensity in males and females was ~ 5%. Broodstocks from different families were randomly mated with each other and 100 to 125 families were retained in each generation. There was one generation of selection per year. In the seventh generation, in 2014, 120 males and 120 females from different families (one male and one female from each family) were selected as experimental parents. In May of that year, each male and two females ready to spawn were maintained in a 1-m3 fiberglass tank for one week. Next, 120 females with fertilized eggs in their mouths were separately placed into multiple hapas (1 m × 1 m × 1 m) in a large concrete pond (50 m × 7 m × 1 m) for a one-week incubation period and were then isolated from their progenies. At 50 days of age, 50 families, each consisting of no less than 1000 surviving progenies above 30 g, were used to construct the experimental population. At the same time, 32 progenies that were randomly chosen from each family were tagged with passive integrated transponder tags, and 1600 tagged fish were mixed in a larger concrete pond (50 m × 20 m × 1.8 m). Because of this, we did not include any effect of the rearing facilities on the measurements in the statistical model. In the subsequent experiment, all fish were fed on a standard commercial diet (crude protein: 28%, crude lipids: 4%, crude fiber: 15%, ash: 18%, total phosphorus: 1%, lysine: 1.2%) manufactured by the Feedstuff Incorporated company (Ningbo, China). Dissolved oxygen was maintained at 8.23 to 8.67 mg/L by air pumps and the water temperature naturally fluctuated between 24.3 and 26.7 °C during the entire experiment.

The tagged fish were weighed on electronic scales to measure BWE, while BL and BD were measured with calipers. Before each measurement, the fish were anesthetized with clove oil at a density of 100 mg/L. The fish were measured once every 15 days during the experiment, with the first recording at the time of tagging. A maximum of six records for each fish and each trait were available, since 4.6% of the individuals died before the end of the experiment. Although experimental fish were measured only six times, the traits were recorded on 20 different days of age because of different spawning times. A total of 7560 records were extracted for BWE, BL and BD.

To edit the dataset, fish with missing phenotypic records for any trait and no information on the parents were excluded and abnormal phenotypic values that were more than three standard from the phenotypic mean of the trait were removed. After editing, 1451 individuals with 7235 records remained, with each individual having at least two records per trait. Table 1 presents the descriptive statistics for the dataset for BWE, BL, and BD for the 20 age groups.

Table 1 Averages and standard deviations (in parentheses) for body weight (g), body length (mm) and body depth (mm) at 20 days of age

Pedigree data were collected by tracing back three generations, thus including 1604 individuals from 77 sires and 88 dams. Although maternal effects on growth traits may be large in the early growth period for tilapia, they were not considered in the analyses, because observations from only a single generation and from such a small number of dams without records render the maternal effects unidentifiable.

Random regression model

In the same setting of cultivation programs, the two sexes and six test days were considered as fixed effects, while family, additive genetic, and permanent environmental effects were considered as random effects. For the growth-related traits that were recorded at multiple time-points, the phenotype of the trait at \(t\) days of age was described by the following single-trait animal model:

$$y_{ijkl} \left( t \right) = \mu_{j} + b_{k} \left( t \right) + f_{l} \left( t \right) + a_{i} \left( t \right) + pe_{i} \left( t \right) + e_{ijkl} (t),$$
(1)

where \(y_{ijkl} \left( t \right)\) is an observation at \(t\) days of age, \(\mu_{j}\) is the \(j\)th test day effect, \(b_{k} \left( t \right)\) is the \(k\)th sex effect, \(f_{l} \left( t \right)\) is the \(l\)th family effect, \(a_{i} \left( t \right)\) is the additive genetic effect for the \(i\)th individual, \(pe_{i} \left( t \right)\) is the permanent environmental effect for the \(i\)th individual, and \(e_{ijkl} (t)\) is the residual error following a normal distribution with expectation 0 and variance \(\sigma_{e}^{2}\).

To analyze the growth curves of phenotypes that are measured repeatedly, Legendre polynomials are generally used to fit the changes in the fixed and random effects in Model (1) and are represented by:

$$L_{m} \left( t \right) = \beta_{0} + \beta_{1} \psi_{1} \left( t \right) + \cdots + \beta_{k} \psi_{k} \left( t \right) + \cdots + \beta_{m} \psi_{m} (t),$$

where \(m\) is the order, \(\beta_{k}\) is the \(k\)th regression coefficient, and \(\psi_{k} \left( t \right)\) is the \(k\)th covariate of the Legendre polynomial [2], defined by:

$$\psi_{k} \left( t \right) = \frac{1}{{2^{k} }}\mathop \sum \limits_{i = 0}^{int(k/2)} \frac{{\left( { - 1} \right)^{i} \left( {2k - 2i} \right)!}}{{i!\left( {k - i} \right)!\left( {k - 2i} \right)!}}\tau^{k - 2i} ,$$

where \(\tau = -\,2 \times \frac{{t - { \hbox{min} }(t)}}{{\hbox{max} \left( t \right)\,\,-\,\,{ \hbox{min} }(t)}}\).

By modeling changes with age in the effects of sex, family, additive genetics, and permanent environment of the Legendre polynomials with orders \(m_{1}\), \(m_{2}\), \(m_{3}\), and \(m_{4}\), respectively, Model (1) becomes the following RRM:

$$y_{ijkl} \left( t \right) = \mu_{j} + L_{m1}^{b} \left( t \right) + L_{m2}^{f} \left( t \right) + L_{m3}^{a} \left( t \right) + L_{m4}^{pe} \left( t \right) + e_{ijkl} \left( t \right),$$
(2)

where \(L_{\eta }^{\theta } \left( t \right) = \beta_{0}^{\theta } + \beta_{1}^{\theta } \psi_{1} \left( t \right) + \beta_{2}^{\theta } \psi_{2} \left( t \right) + \cdots + \beta_{\eta }^{\theta } \psi_{\eta } (t)\), with \(\theta \in \left[ {b f a pe} \right]\) and \(\eta \in [m_{1} m_{2} m_{3} m_{4} ]\).

In matrix notation, Model (2) can be written as:

$${\mathbf{y}}_{ijkl} = \mu_{j} + {\mathbf{x}}_{1i} {\mathbf{b}}_{l} + {\mathbf{x}}_{2i} {\mathbf{f}}_{k} + {\mathbf{x}}_{3i} {\mathbf{a}}_{i} + {\mathbf{x}}_{4i} {\mathbf{pe}}_{i} + {\mathbf{e}}_{ijkl} ,$$
(3)

where \({\mathbf{x}}_{\eta i} = [1 \psi_{1} \left( {t_{i} } \right) \ldots \psi_{\eta } (t_{i} )]\), and \({\varvec{\uptheta}} \in \left[ {{\mathbf{b}}\, {\mathbf{f}} \,{\mathbf{a}} \,{\mathbf{pe}}} \right] = [\beta_{0}^{\theta } \beta_{1}^{\theta } \ldots \beta_{\eta }^{\theta } ]^{T}\).

The model satisfied the following:

$$\begin{aligned} & E\left( {{\mathbf{y}}_{ijk} |\mu_{j} , {\mathbf{b}}_{l} , {\mathbf{f}}_{k} , {\mathbf{a}}_{i} , {\mathbf{pe}}_{i} } \right) = \mu_{j} + {\mathbf{x}}_{ij} {\mathbf{b}}_{l} + {\mathbf{z}}_{ik} {\mathbf{f}}_{k} + {\mathbf{z}}_{i} {\mathbf{a}}_{i} + {\mathbf{z}}_{i} {\mathbf{pe}}_{i} \\ & Cov\left( {{\mathbf{y}}_{ijk} } \right) = {\mathbf{I}}\,{ \otimes }\,{\mathbf{F}} + {\mathbf{A}}\,{ \otimes }\,{\mathbf{G}} + {\mathbf{I}}\,{ \otimes }\,{\mathbf{PE}} + {\mathbf{I}}\,{ \otimes }\,\sigma_{e}^{2} , \\ \end{aligned}$$

where \({\mathbf{A}}\) is the numerator relationship matrix and \({\mathbf{F}}\), \({\mathbf{G}}\), and \({\mathbf{P}}\) are the family, additive genetic, and permanent environmental covariance matrices for the random regression coefficients of the Legendre polynomials, respectively.

To estimate the genetic correlations between growth traits measured at multiple time-points, the RRM for those traits must be solved simultaneously by a multivariate genetic analysis. The covariance matrices in the MRRM were estimated by using restricted maximum likelihood (REML), as implemented in DMU Version 6 [29]. The starting values were set to 0 for each fixed effect, to the identity matrix for each random effect covariance matrix and to 1 for the residual variance. The convergence criterion for REML was set to 10−6.

After estimating the covariance matrices, heritabilities for the traits at any age and the genetic correlations between traits measured at different ages were estimated by the covariance functions [30]. The covariances for a trait between the \(i\)th and \(j\)th days of age can be calculated as \({\mathbf{z}}_{i} {\mathbf{Fz}}_{j}^{\text{T}}\) for the family effect, \({\mathbf{w}}_{i} {\mathbf{Aw}}_{j}^{\text{T}}\) for the additive genetic effect, and \({\mathbf{s}}_{i} {\mathbf{PEs}}_{j}^{\text{T}}\) for the permanent environmental effect. Therefore, genetic correlations between the \(i\)th and \(j\)th days of age for a trait were estimated by \(\frac{{{\mathbf{w}}_{i} {\mathbf{Aw}}_{j}^{\text{T}} }}{{\sqrt {({\mathbf{w}}_{i} {\mathbf{Aw}}_{i}^{\text{T}} )({\mathbf{w}}_{j} {\mathbf{Aw}}_{j}^{\text{T}} )} }}\), phenotypic correlations by \(\frac{{{\mathbf{P}}_{ij} }}{{\sqrt {{\mathbf{P}}_{ii} {\mathbf{P}}_{jj} } }}\) with phenotypic covariance \({\mathbf{P}}_{ij} = {\mathbf{z}}_{i} {\mathbf{Fz}}_{j}^{\text{T}} + {\mathbf{w}}_{i} {\mathbf{Aw}}_{j}^{\text{T}} + {\mathbf{s}}_{i} {\mathbf{PEs}}_{j}^{\text{T}} + \sigma_{e}^{2} {\mathbf{I}}(i = j)\), and heritability at the \(i\)th day of age by \(\frac{{{\mathbf{w}}_{i} {\mathbf{Aw}}_{i}^{\text{T}} }}{{{\mathbf{P}}_{ii} }}\).

Choice of model

When choosing a model for a longitudinal analysis, Stoica and Babu [31] observed that neither the Akaike information criterion nor the Bayesian information criterion had optimum properties in terms of consistency and efficiency, and therefore the authors introduced a new criterion, i.e. penalizing adaptively the likelihood (PAL) [31], that overcomes these issues. More recently, Corrales et al. [32] applied PAL to the selection of the order of the Legendre polynomial in RRM. Consider a set of \(m\) competing models \(M_{1} , L, M_{k - 1} , M_{k} , L, M_{m}\) that are ranked in ascending order by the number of parameters. The PAL criterion for selecting Model \(k\) is defined as:

$$PAL_{k} = -\,2\log \left( {ML_{k} } \right) + p_{k} { \log }(p_{m} )\frac{{{ \log }(r_{k} + 1}}{{{ \log }(\rho_{k} + 1)}},$$
(4)

where \(\log \left( {ML_{k} } \right)\) is the logarithm of the maximum likelihood value of Model \(k\); \(p_{k}\) is the number of parameters of Model \(k\); \(p_{m}\) is the largest number of parameters of competing models; and \(r_{k} = 2\log \left( {ML_{k - 1} } \right) - 2\log \left( {ML_{1} } \right)\) and \(\rho_{k} = 2\log \left( {ML_{m} } \right) - 2{ \log }(ML_{k - 1} )\), with \(\log \left( {ML_{1} } \right)\), \(\log \left( {ML_{k - 1} } \right)\) and \(\log \left( {ML_{m} } \right)\) being the logarithms of the maximum likelihood values of Models \(1\), \(k {-} 1\), and \(m\), respectively.

Results

Choice of MRRM

The best MRRM for BWE, BL, and BD was selected by separately optimizing the univariate RRM for growth traits. The Legendre polynomial for the effect of sex was generally modeled by the population’s growth curve, chosen as the Legendre polynomial of three orders according to the highest goodness-of-fit for the three analyzed traits. The three random effects, i.e. family, additive genetic, and permanent environmental effects, was characterized by using Legendre polynomials of different orders that ranged from 0 to 3 [3]. The models were designated as \({\text{LPm}}_{1} {\text{m}}_{2} {\text{m}}_{3}\), e.g., \({\text{LP}}121\) is a model with a Legendre polynomial of order 1 for the family effect, order 2 for the additive genetic effect, and order 1 for the permanent environmental effect. A total of 64 RRM for each analyzed trait were compared based on the PAL criterion to determine the best combination of orders of the Legendre polynomials for the three random effects.

When fitting RRM with constant permanent environmental effects, the variance estimates of these effects were statistically inferred as not significant, using student t statistics (estimate/standard error). In addition, RRM with dynamic permanent environmental effects converged poorly in DMU. Therefore, permanent environmental effects were not considered in competing models. In the remainder, for the nine competing models, model LP00, which had constant family and additive effects, was taken as the reference, because the \(- 2\log \left( {ML} \right)\) for the RRM with dynamic additive genetic and family effects was consistently lower than for model LP00. Based on Eq. (4), Table 2 presents the calculated PAL values for the nine competing RRM. Among these, \({\text{LP}}23\) was chosen for BWE with Legendre polynomials of order 2 for the family effect and order 3 for the additive genetic effect, while \({\text{LP}}13\) and \({\text{LP}}22\) were selected for BL and BD, respectively.

Table 2 Choice of RRM for BWE, BL and BD based on the PAL criterion

Genetic parameters for growth traits

Based on the model choices as described above, the final MRRM was constructed to jointly analyze BWE, BL, and BD. Based on Model (2), the MRRM, in which Legendre polynomials of different orders were nested, was expressed as:

$$\left\{ {\begin{array}{*{20}c} {y_{1} \left( t \right) = \mu_{1} + L_{3}^{{b_{1} }} \left( t \right) + L_{2}^{{f_{1} }} \left( t \right) + L_{3}^{{a_{1} }} \left( t \right) + e_{1} (t)} \\ {y_{2} \left( t \right) = \mu_{2} + L_{3}^{{b_{2} }} \left( t \right) + L_{1}^{{f_{2} }} \left( t \right) + L_{3}^{{a_{2} }} \left( t \right) + e_{2} (t)} \\ {y_{3} \left( t \right) = \mu_{3} + L_{3}^{{b_{3} }} \left( t \right) + L_{2}^{{f3_{{}} }} \left( t \right) + L_{2}^{{a_{3} }} \left( t \right) + e_{3} (t)} \\ \end{array} } \right.,$$
(5)

where \(y_{1} \left( t \right)\), \(y_{2} \left( t \right)\), and \(y_{3} \left( t \right)\) are the phenotypes of BWE, BL, and BD, respectively.

Covariance estimates

Estimates of the (co)variances (standard errors) of the random regression coefficients were for the family effects and the genetic additive effects as shown in Figs. 1 and 2, respectively.

Fig. 1
figure 1

Estimates of the (co)variances (standard errors) of the random regression coefficients for family effects

Fig. 2
figure 2

Estimates of the (co)variances (standard errors) of the random regression coefficients for additive genetic effects

Estimates of the (co)variances (standard errors) of the random regression coefficients were for the random errors as follows:

$${\mathbf{RE}} = \left[ {\begin{array}{*{20}c} {42.06 (1.0)} & {6.2} & {5.8} \\ {6.2 (0.3)} & {7.5 (0.2)} & {1.2} \\ {5.8 (0.2)} & {1.2 (0.1)} & {3.7 (0.1)} \\ \end{array} } \right].$$

For each estimated covariance matrix of the random regression coefficients (including β 0 ), the diagonal block matrices are the covariance matrices for BWE, BL, and BD, respectively, while the off-diagonal blocks are the covariance matrices between the three growth traits. All covariances of the random regression coefficients were inferred as significantly different from 0. Estimates of the genetic correlations between most of the random regression coefficients for the additive genetic effects were positive, and the few negative genetic correlations were mainly between the intercept and the cubic regression coefficients or between the linear and cubic regression coefficients. In addition, estimates of genetic covariances for the same order random regression coefficients between the pairwise growth traits were all positive and decreased with increasing order.

Heritability estimates

Table 3 lists the estimated ratios of family variances to phenotypic variances and heritabilities for selected days of age, using the univariate and multivariate RRM. The estimated ratios of family variances to phenotypic variances for BWE ranged from 0.259 and 0.365 and were consistently lower than these ratios for BL and BD for the period between 60 and 140 days of age. The estimated heritability obtained with the MRRM ranged from 0.080 to 0.635 for BWE, from 0.173 to 0.562 for BL, and from 0.098 to 0.610 for BD for the period between 60 and 140 days of age. Heritability estimates obtained from the MRRM were slightly lower than those obtained from univariate RRM for BWE and BL, but consistently higher for BD. Figure 3 shows that the heritabilities estimated with the MRRM increased with age for all three growth traits. Heritabilities changed most rapidly for BD and most slowly for BL. Heritabilities for BL were consistently higher than those for the other two traits from 60 to 75 days of age and those for BWE were higher than those for the other two traits from 75 days of age onward. Considering 120 days of age as the time to market, heritabilities were above moderate and similar among the three growth traits.

Table 3 Estimates of ratios of family to phenotypic variances (f 2) and of heritabilities (h 2) for BWE, BL and BD at selected days of age based on the univariate and multivariate RRM
Fig. 3
figure 3

Estimated heritabilities by age for the three growth traits. Square for BWE, diamond for BL, and triangle point-up for BD

Correlation estimates

Figure 4 presents a three-dimensional plot of estimates of genetic correlations between pairwise traits measured at different ages for each of the three growth traits from 60 to 140 days of age. Patterns of the genetic correlations over time were similar for BWE, BL, and BD. Genetic correlation estimates between measurements at adjacent days of age were close to 1 but decreased monotonously as the lag between the days of age increased. In addition, genetic correlations between measurements at early days of age were lower than those at later days of age. The lowest genetic correlations were 0.566 for BWE, 0.633 for BL, and 0.497 for BD between 60 and 140 days of age.

Fig. 4
figure 4

Estimates of genetic correlations between traits measured at different days of age (ra) for the three growth traits

Table 4 includes family and phenotypic correlations between the pairwise selected growth days for BWE. Similar to the genetic correlations, phenotypic correlations approached 0.8 for the period between 75 and 120 days of age. In the earlier growth period, genetic correlations were higher than the phenotypic correlations between the pairwise days of age. The minimum family correlations were found between growth traits at 60 and 140 days of age, which were 0.573, 0.702, and 0.664 for BWE, BL, and BD, respectively. Family and phenotypic correlations for BL and BD were similar to those for BWE (see Additional file 1: Tables S1 and S2).

Table 4 Estimates of phenotypic (lower triangle) and family (upper triangle) correlations for body weight between pairwise measurements at selected days of age

Patterns of the genetic correlations between the pairwise growth traits were similar to those between the pairwise days of age for the same trait, but they showed lower ridges, because the genetic correlation of the same trait on the same day is 1 (Fig. 5). Genetic correlations between traits at the same age ranged from 0.89 at 140 days of age to 0.95 at 60 days of age between BWE and BL, from 0.89 to 0.96 between BWE and BD, and from 0.85 to 0.95 between BL and BD. In contrast, genetic correlations between traits at different days of age were lowest between BWE at 140 days of age and BL at 60 days of age (0.50), between BWE at 140 days of age and BD at 60 days of age (0.47), and between BL at 60 days of age and BD at 140 days of age (0.44). At harvest (120 days of age), genetic correlations were equal to 0.92 between BWE and BL, 0.94 between BWE and BD, and 0.89 between BL and BD.

Fig. 5
figure 5

Estimates of genetic correlations (rg) between growth traits

Family and phenotypic correlations between the pairwise growth traits exhibited similar patterns as the genetic correlations but different magnitudes at the same compared days of age (Fig. 6). Family correlations were consistently higher than genetic correlations between BWE and BD and between BL and BD. For the same day of age, most phenotypic correlations of pairwise growth traits were lower than the corresponding genetic and family correlations (see Additional file 1: Tables S3, S4 and S5).

Fig. 6
figure 6

Estimates of family correlations (rf) between growth traits

Discussion

In this study, we used an MRRM to estimate heritabilities of BWE, BL, and BD based on repeated measurements during growth period, as well as genetic correlations between pairwise growth traits at specific days of age. Legendre polynomials were chosen to characterize the influence of fixed and random effects, i.e. family, additive genetic, and permanent environmental effects, on growth curves. The best MRRM was established by separately optimizing univariate RRM for each growth trait, according to the PAL. Based on the final MRRM, we found that heritabilities increased with age for all traits and that they were slightly lower than those obtained by using univariate RRM. In addition, for each trait, genetic correlations between measurements decreased monotonously with increased lag in days of age. For measurements at the same days of age, estimates of heritabilities and genetic correlations were close to those previously reported for tilapia [4,5,6, 8, 10,11,12,13,14,15,16, 18, 20,21,22, 33, 34].

For the selected RRM, changes in the genetic effects with day of age were optimally modeled for all analyzed traits by using Legendre polynomials of three orders. In contrast, Rutten et al. [9] and Turra et al. [27] used polynomials of two orders to fit the fixed and random effects, without selecting the orders of the polynomials for the different fixed and random effects. The non-significant variance for the permanent environmental effects showed that they have no impact on growth traits. Ratios of family variance [9] and of permanent environmental variance [27] to phenotypic variance are largely underestimated in a population with multiple families, which may be caused by the collinearity between the family and permanent environmental effects in the RRM used.

Repeated measurements during growth are required to estimate changes in both the fixed and the random effects with age when applying the RRM to genetic analysis of growth traits. More longitudinal measurements per individual would help to model growth curves with a higher goodness-of-fit and to more robustly estimate genetic parameters for random regression. However, more longitudinal measurements not only increase the experimental costs, but also affect fish growth, especially when all experimental individuals are measured together. In our experiment, the simultaneous measurement of 1600 tilapia fish incurred high labor costs and a measuring frequency of once every 15 days influenced fish growth to some extent. In future trials, the experimental population could be divided into several separately reared subpopulations, and each subpopulation could be observed in turn over more growth time-points. The extension of such measurements to more families from multiple generations would also help to identify maternal effects on early growth in tilapia within the MRRM framework.

Conclusions

Using repeated records of growth duration from multiple families in one generation, first we introduced MRRM to genetically analyze growth curves of body weight and main morphological traits in genetically improved farmed tilapia. The optimal MRRM was chosen by comparing the univariate RRM for the analyzed traits separately via PAL. By using the covariance functions of the MRRM, changes in heritabilities with days of age and genetic correlations between lags in days of age were estimated, which could be used to carry out early selection for each trait analyzed. More importantly, for phenotypes that are measured repeatedly in aquaculture, an MRRM can enhance the efficiency of the comprehensive selection for BWE and the main morphological traits at a specific age.

References

  1. Yang Y, Li R, Li S. An estimation method for genetic parameters of dynamic traits. Acta Vet Zootech Sin. 1996;27(2):136–42.

    Google Scholar 

  2. Kirkpatrick M, Heckman N. A quantitative genetic model for growth, shape, reaction norms, and other infinite-dimensional characters. J Math Biol. 1989;27:429–50.

    Article  CAS  PubMed  Google Scholar 

  3. Schaeffer LR. Application of random regression models in animal breeding. Livest Prod Sci. 2004;86:35–45.

    Article  Google Scholar 

  4. Bolivar RB, Newkirk GF. Response to within family selection for body weight in Nile tilapia (Oreochromis niloticus) using a single-trait animal mode. Aquaculture. 2002;204:371–81.

    Article  Google Scholar 

  5. Gall GAE, Bakar Y. Application of mixed-model techniques to fish breed improvement: analysis of breeding-value selection to increase 98-day body weight in tilapia. Aquaculture. 2002;212:93–113.

    Article  Google Scholar 

  6. Ponzoni RW, Hamzah A, Tan S, Kamaruzzaman N. Genetic parameters and response to selection for live weight in the GIFT strain of Nile Tilapia (Oreochromis niloticus). Aquaculture. 2005;247:203–10.

    Article  Google Scholar 

  7. Rutten MJM, Bovenhuis H, Komen H. Modeling fillet traits based on body measurements in three Nile tilapia strains (Oreochromis niloticus L). Aquaculture. 2004;231:113–22.

    Article  Google Scholar 

  8. Rutten MJM, Bovenhuis H, Komen J. Genetic parameters for fillet traits and body measurements in Nile tilapia (Oreochromis niloticus L.). Aquaculture. 2005;246:125–32.

    Article  Google Scholar 

  9. 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.

    Article  Google Scholar 

  10. Charo-Karisa H, Komen H, Rezk MA, Ponzoni RW, Arendonk JAMV, Bovenhuis H. Heritability estimates and response to selection for growth of Nile tilapia (Oreochromis niloticus) in low-input earthen ponds. Aquaculture. 2006;261:479–86.

    Article  Google Scholar 

  11. Maluwa AO, Gjerde B, Ponzoni RW. Genetic parameters and genotype by environment interaction for body weight of Oreochromis shiranus. Aquaculture. 2006;259:47–55.

    Article  Google Scholar 

  12. Eknath AE, Bentsen HB, Ponzoni RW, Rye M, Nguyen NH, Thodesen J, et al. Genetic improvement of farmed tilapias: composition and genetic parameters of a synthetic base population of Oreochromis niloticus for selective breeding. Aquaculture. 2007;273:1–14.

    Article  CAS  Google Scholar 

  13. Nguyen NH, Khaw HL, Ponzoni RW, Hamzah A, Kamaruzzaman N. Can sexual dimorphism and body shape be altered in Nile tilapia (Oreochromis niloticus) by genetic means? Aquaculture. 2007;272:S38–46.

    Article  Google Scholar 

  14. Khaw HL, Ponzoni RW, Danting MJC. Estimation of genetic change in the GIFT strain of Nile tilapia (Oreochromis niloticus) by comparing contemporary progeny produced by males born in 1991 or in 2003. Aquaculture. 2008;275:64–9.

    Article  Google Scholar 

  15. Khaw HL, Bovenhuis H, Ponzoni RW, Rezk MA, Charo-Karisa H, Komen H. Genetic analysis of Nile tilapia (Oreochromis niloticus) selection line reared in two input environments. Aquaculture. 2009;294:37–42.

    Article  CAS  Google Scholar 

  16. Rezk MA, Ponzoni RW, Khaw HL, Kamel E, Dawood T, John G. Selective breeding for increased body weight in a synthetic breed of Egyptian Nile tilapia, Oreochromis niloticus: response to selection and genetic parameters. Aquaculture. 2009;293:187–94.

    Article  Google Scholar 

  17. Luan TD, Olesen I, Ødegård J, Kolstad K, Dan NC, Elghobashy H, et al. Genotype by environment interaction for harvest body weight and survival of Nile tilapia (Oreochromis niloticus) in brackish and fresh water ponds. In: Proceedings of the 8th international symposium on Tilapia in aquaculture: 12–14 October 2008, Cairo; 2008. p. 133–141.

  18. Nguyen NH, Ponzoni RW, Abu-Bakar KR, Hamzah A, Khaw HL, Hong YY. Correlated response in fillet weight and yield to selection for increased harvest weight in genetically improved farmed tilapia (GIFT strain), Oreochromis niloticus. Aquaculture. 2010;305:1–5.

    Article  Google Scholar 

  19. Thodesen J, Rye M, Wang YX, Yang KS, Bentsen HB, Gjedrem T. Genetic improvement of tilapias in China: genetic parameters and selection responses in growth of Nile tilapia (Oreochromis niloticus) after six generations of multi-trait selection for growth and fillet yield. Aquaculture. 2011;322–323:51–64.

    Article  Google Scholar 

  20. Bentsen HB, Gjerde B, Nguyen NH, Rye M, Ponzoni RW, Vera MSPD, et al. Genetic improvement of farmed tilapias: genetic parameters for body weight at harvest in Nile tilapia (Oreochromis niloticus) during five generations of testing in multiple environments. Aquaculture. 2012;338–341:56–65.

    Article  Google Scholar 

  21. Akar M, Baylan M, Sangün L. A quantitative genetic analysis of six phenotypes in Oreochromis niloticus in different olds. Biyoloji Bilimleri Araştırma Dergisi. 2013;6:62–6.

    Google Scholar 

  22. Luan TD, Olesen I, Kolstad K. Genetic parameters and genotype by environment interaction for growth of Nile tilapia in low and optimal temperature. In: Proceedings of the 9th World congress on genetics applied to livestock production, 1–6 August 2006, Leipzig; 2010.

  23. Valente BD, Silva MA, Silva LOC, Bergmann JAG, Pereira JCC, Fridrich AB, et al. Estruturas de covariância de peso em função da idade de animais Nelore das regiões Sudeste e Centro-Oeste do Brasil. Arq Bras Med Vet Zootec. 2008;60:389–400.

    Article  Google Scholar 

  24. Kirkpatrick M, Lofsvold D, Bulmer M. Analysis of the inheritance, selection and evolution of growth trajectories. Genetics. 1990;124:979–93.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Henderson CR. Application of linear models in animal breeding. Guelph: University of Guelph; 1984.

    Google Scholar 

  26. Mckay LR, Schaeffer LR, Mcmillan I. Analysis of growth curves in Rainbow trout using random regression. In: Proceedings of 7th World congress on genetics applied to livestock production, 19–23 August 2002, Montpellier; 2002.

  27. Turra EM, de Oliveira DAA, Valente BD, Teixeira EDA, Prado SDA, Melo DCD, et al. Estimation of genetic parameters for body weights of Nile tilapia Oreochromis niloticus using random regression models. Aquaculture. 2012;354–355:31–7.

    Article  Google Scholar 

  28. Meyer K. Scope for a random regression model in genetic evaluation of beef cattle for growth. Livest Prod Sci. 2004;86:69–83.

    Article  Google Scholar 

  29. Madsen P, Jensen J. A user’s guide to DMU. A package for analysing multivariate mixed models. Foulum: National Institute of Animal Science; 1997.

    Google Scholar 

  30. Meyer K, Hill WG. Estimation of genetic and phenotypic covariance functions for longitudinal or ‘repeated’ records by restricted maximum likelihood. Livest Prod Sci. 1997;47:185–200.

    Article  Google Scholar 

  31. Stoica P, Babu P. Model order estimation via penalizing adaptively the likelihood (PAL). Signal Process. 2013;93:2865–71.

    Article  Google Scholar 

  32. Corrales JD, Munilla S, Cantet RJ. Polynomial order selection in random regression models via penalizing adaptively the likelihood. J Anim Breed Genet. 2015;132:281–8.

    Article  CAS  PubMed  Google Scholar 

  33. Guo HU, Wang F, Zhou J, Pan Y, Baohai LI, Wei GU, et al. Study on genetic diversity and germplasm sources of wild Yadong salmon population. J Northeast Agric Univ. 2016;1:58–65.

    Google Scholar 

  34. Tibshirani R. Regression shrinkage and selection via the lasso: a retrospective. J R Stat Soc Series B. 2011;73:273–82.

    Article  Google Scholar 

Download references

Authors’ contributions

RY proposed the method and supervised JH, JZ and JG for writing the manuscript. JH, YZ, JZ and JG wrote computer codes and conducted real data analysis. DH and PX took part in the interpretation of results. JH and PX constructed the breeding population and observed traits. All authors read and approved the final manuscript.

Acknowledgements

This study was funded by Special Scientific Research Funds for Central Non-profit Institutes, Chinese Academy of Fishery Sciences (2014B004 and 2014A04XK01) and the National Natural Science Foundations of China (31201989).

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Runqing Yang.

Additional file

12711_2017_357_MOESM1_ESM.docx

Additional file 1: Table S1. Phenotypic (Lower triangle) and family (Upper triangle) correlations of body lengths between pairwise the selected days of age. Table S2 Phenotypic (Lower triangle) and family (Upper triangle) correlations of body depths between pairwise the selected days of age. Table S3 Phenotypic correlations between body weights and body lengths at the selected days of age. Table S4 Phenotypic correlations between body weights and body depths at the selected days of age. Table S5 Phenotypic correlations between body lengths and body depths at the selected days of age.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, J., Zhao, Y., Zhao, J. et al. Multivariate random regression analysis for body weight and main morphological traits in genetically improved farmed tilapia (Oreochromis niloticus). Genet Sel Evol 49, 80 (2017). https://doi.org/10.1186/s12711-017-0357-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-017-0357-7

Keywords