Genetic evaluation of European quails by random regression models

BREEDING, GENETIC AND REPRODUCTION

Genetic evaluation of European quails by random regression models¹

Flaviana Miranda Gonçalves^I; Aldrin Vieira Pires^II,V; Idalmo Garcia Pereira^III; Eduardo Silva Cordeiro Drumond^I; Vivian Paula Silva Felipe^IV; Sandra Regina Freitas Pinheiro^II

^IPrograma de Pós-Graduação em Zootecnia - Departamento de Zootecnia/UFVJM, 39100-000, MG, Brazil

^IIDepartamento de Zootecnia/UFVJM, 39100-000, Diamantina, MG, Brazil

^IIIDepartamento de Zootecnia/Escola de Veterinária, 30123-970,UFMG, MG, Brazil

^IVPrograma de Pós-Graduação em Zootecnia - Departamento de Zootecnia/Escola de Veterinária, 30123-970- UFMG, MG, Brazil

^VSupported by CNPq and FAPEMIG

ABSTRACT

The objective of this study was to compare different random regression models, defined from different classes of heterogeneity of variance combined with different Legendre polynomial orders for the estimate of (co)variance of quails. The data came from 28,076 observations of 4,507 female meat quails of the LF1 lineage. Quail body weights were determined at birth and 1, 14, 21, 28, 35 and 42 days of age. Six different classes of residual variance were fitted to Legendre polynomial functions (orders ranging from 2 to 6) to determine which model had the best fit to describe the (co)variance structures as a function of time. According to the evaluated criteria (AIC, BIC and LRT), the model with six classes of residual variances and of sixth-order Legendre polynomial was the best fit. The estimated additive genetic variance increased from birth to 28 days of age, and dropped slightly from 35 to 42 days. The heritability estimates decreased along the growth curve and changed from 0.51 (1 day) to 0.16 (42 days). Animal genetic and permanent environmental correlation estimates between weights and age classes were always high and positive, except for birth weight. The sixth order Legendre polynomial, along with the residual variance divided into six classes was the best fit for the growth rate curve of meat quails; therefore, they should be considered for breeding evaluation processes by random regression models.

Key Words:Coturnix coturnix coturnix, genetic parameters, heterogeneity of variance

Introduction

The growing market for products of animal origin with high protein content and good acceptance by the consumer has caused the meat quail production to increase as well. However, meat quail production faces a lack of national genetic material. Currently, genetic material with less productive potential available in the market is used, which leaves the national market dependent on genetic material from other countries.

Therefore, breeding programs are important, since they estimate breeding parameters and establish selection strategies to renew the herds with animals of higher breeding potential at every new generation.

The main piece of information used in processes of genetic evaluation and selection is body weight during the growth period. Meat quail studies determine these weights at regular intervals, so the same animal has several weights determined during its life, thus characterizing repeated measurements or longitudinal data. This repeated measurement of the same animal over time has received increasing interest of animal breeding researchers. Repeating traits have been studied using covariance functions. Among them, random regression models have received more attention from animal breeding researchers (Tholon & Queiroz, 2008).

Random regression models enable the prediction of breeding values for the interval at which measurements are taken and to estimate individual growth rate curves considering the relationship between animals, so they promote better use of the data and potential increase in the selection accuracy (Dionello et al., 2008).

According to El Faro & Albuquerque (2003), multiple regression analyses that use covariance functions also allow one to consider the heterogeneity of residual variances on the control day. Modeling can improve the partition of total variance, which increases accuracy of estimated genetic parameters, maximizing the response to selection.

Therefore, the objective of this study was to compare different random regression models defined from different classes of variance heterogeneity combined with different orders of Legendre polynomial for the genetic evaluation of meat quails.

Material and Methods

Female meat quails of the LF1 lineage (Coturnix coturnix coturnix) from Programa de Melhoramento Genético de Codornas of the Department of Animal Science, Universidade Federal dos Vales do Jequitinhonha e Mucuri, in Diamantina, MG, Brazil, were utilized in this experiment. For the analyses, a total of 28,076 records of body weights from 4,510 meat quails and 4,681 animals in the relationship matrix, taken every seven days (at birth, 7, 14, 21, 28, 35 and 42 days of age), collected from June 2010 to May 2011 were used. Three generations from the base population were studied; the first with 793 birds, the second with 2,189 and the third with 1,528 quails.

Quails with at least three weight records were considered for the analysis. They belonged to 174 contemporary groups, consisting of generation, hatching and sex. Quails were reared from birth to 42 days of age in pens with concrete floor and wood shaving beds, equipped with heating. Birds received feed and water ad libitum. From birth to 21 days of age, the diet supplied contained 25% CP and 2,900 kcal ME/kg, and from 22 to 42 days, the diet contained 24% CP and 2,925 kcal ME/kg, according to nutritional requirements specified in the literature (Oliveira et al., 2002; Fridrich et al., 2005; Corrêa et al., 2007; Veloso et al., 2010).

At birth, quails were identified, weighed and transferred to the pens. At 35 days of age, 210 females and 105 males were selected according to body weight and then transferred to individual cages to form the next generation on. The selected quails remained in the galvanized wire cage (95 cm long, 28.5 cm wide and 21.5 cm high) equipped with egg collector (15.5 cm in front of the cages), repeating the entire cycle every generation.

The three hatchings of the first generation totaled 793 birds; the five hatchings of the second generation totaled 2,189 birds, while the four hatchings of the third generation totaled 1,528 birds. The hatches of the same generation were collected every seven days and the quails were reared under the same management practices described before.

Body weights were analyzed by random regression models. The fixed and random regression models were represented by continuous functions, whose ages were described by the Legendre orthogonal polynomials. The different order (2nd, 3rd, 4th, 5th and 6th order) models adjusted for the random effects can be represented as follows:

in which y_ij = body weight on day j of quail i; F = set of fixed effects, constituted by the contemporary group (174 subclasses); β_m = fixed regression weight coefficient of the Legendre polynomial m; k_β = polynomial order for the mean curve, varying from two to six (k_β = 2, 3, ... , 6) to obtain the average growth rate curve of the quail population; α_im and ρ_im are animal additive genetic and permanent environmental regression coefficients, respectively, for quail i; k_α and k_ρ = Legendre adjusted polynomial orders; φ_m = Legendre polynomial function of the standardized age m (-1 < age < 1); and ε_ij = residual random effect.

The model above can be described by the following matrix:

Var(a) = k_aA; Var(p) = k_pI_Nd; Var(e) = R in which y = vector of N observations on 4,510 animals; β = fixed effect vector that includes the fixed regression coefficients used to model the mean trajectory of the population; α = vector k_a × N_d of additive genetic random regression coefficients, in which N_d = total number of animals in the relationship matrix (4,681); p = vector k_p × N of permanent environmental random regression coefficient; e = random error vector; X, Z₁ and Z₂ = incident matrices of fixed regression coefficients, animal additive genetic and permanent environmental random regression coefficients, respectively; k_a and k_p = covariance matrices between the random regression coefficients for animal additive genetic and permanent environmental effects, respectively. A is the matrix of numerators of coefficient of relationship between individuals; I_Nd = size identity matrix N_d; = direct product operator; and R = diagonal matrix containing the residual variances.

The residual variance model considered classes of residual variance from homogeneous to heterogeneous, according to different ages: class 1, homogeneous residual variance (1 to 42 days of age); class 2, heterogeneous variances in two different periods (1 to 7 days and 8 to 42 days of age); class 3, heterogeneous variances in three periods (1 to 7, 8 to 14 and 15 to 42 days of age); class 4, heterogeneous variances in four periods (1 to 7, 8 to 14, 15 to 21 and 22 to 42 days of age); class 5, heterogeneous variances in five periods (1 to 7, 8 to 14, 15 to 21, 22 to 28 and 29 to 42 days of age); and, class 6: heterogeneous variances in six periods (1 to 7, 8 to 14, 15 to 21, 22 to 28, 29 to 35 and 36 to 42 days of age), totaling six different classes of variance heterogeneity to be used in the data analysis. Each of the six classes combined with the different orders of Legendre polynomial (2nd, 3rd, 4th, 5th and 6th order) resulted in 30 models to be tested.

The models were compared by the logarithm of the likelihood function (Log_e L), the likelihood ratio, Akaike Information Criterion (AIC = -2Log_eL +2p) and the Schwarz's Bayesian Information Criterion (BIC = -2Log_eL + pLog_e(N - r)) tests, in which: p = number of model parameters; N = total number of observations; r = rank of the incidence matrix for the fixed effects; and Log_e L = logarithm of the likelihood function. For the likelihood ratio test (LRT), 1% probability was adopted and the statistics was given by the equation: LRT_ij = 2Log_eL_i - 2Log_eL_j, in which Log_eL_i = maximum likelihood function for the complete model i and Log_eL_j is the maximum likelihood function for the reduced model. The LRT estimate was compared to the given Chi-square value with d degrees of freedom at 1% significance level, in which d is the difference between the number of parameters estimated by the complete and reduced models.

The AIC and BIC tests allow to compare non-nested models and both penalize models with more parameters, but BIC is the most stringent and tends to favor more parsimonious models (Wolfinger, 1993; Nunez-Antón & Zimmerman, 2000). Low values for AIC and BIC and high values for Log_eL indicate good fit.

The covariances between random regression coefficients for animal additive genetic and permanent environmental effects, according to the best fit model by the chosen criteria, were estimated by restricted maximum likelihood (REML) using the software Wombat (Meyer, 2006).

After estimating the covariance components necessary to obtain the solutions of the system of equations of the mixed model, the estimates of the matrices G and P, both 6 × 6, became known and therefore possible to represent the covariance structure effects throughout the studied interval using the covariance functions. After structuring the covariances, one can estimate the variance or covariance values for any point or combination points within the age range considered, to determine heritability and correlations as well.

The estimates of covariances for animal additive genetic and permanent environmental effects for animal weight on day t were obtained by:

Weight heritability of the birds can be estimated at any age within the range of the interval studied. Direct heritability represents the fraction of phenotypic variance explained by variance of the additive genetic value of studied birds and was obtained as follows:

The genetic and permanent environmental correlations with weights of the different weighing days, t_i and t_j , were estimated by:

The genetic value for each animal was obtained as follows:

VG = b_0j + b_1jt + b_2jt² + b_3jt³ + b_4jt⁴ + b_5jt⁵ in which b₀ is the regression coefficient of the intercept and b₁, b₂, b₃, b₄ and b₅are the regression coefficients of the n-order polynomial regression model, and t is the bird age represented by the Legendre polynomial.

The genetic trend of additive effect for generations was estimated by the linear regression coefficient of mean genetic values of birds in the year of birth, using PROC REG of the software SAS (Statistical Analysis System, version 9.0).

Results and Discussion

Random regression modeling using Legendre polynomial requires the definition of the best fit order for each random effect considered in the analysis model. Thus, the polynomial order was gradually increased (ranging from 2 to 6) in order to determine the minimum order able to describe the covariance structures as a function of time (Table 1).

As the order of the Legendre polynomial increased, the values of Log L, BIC and AIC decreased. For all the studied residual variance classes, the sixth order Legendre polynomial yielded the lowest values of BIC and AIC (Table 1). According to Log L, the highest polynomial order of the model yielded significant values (P<0.01) for the likelihood ratio test (LRT), which indicates that the model with the most parameters (model 66 with 48 parameters), in which K_a= 6 and K_p= 6 for additive genetic and permanent environmental effects, respectively, was the best fit for the experimental data. Therefore, it is recommended to use six classes of variance heterogeneity (one for each week of life) and to adopt the sixth order Legendre polynomial (model 66) when studying the growth rate curves of meat quail when considering a total of seven weight data.

Similar results were reported by Bonafé et al. (2011), who also reported the sixth order Legendre polynomial as the best fit for growth rate curves of meat quail. Akbas et al. (2004), when assessing the weight of meat quails, also reported best fit for the sixth order Legendre polynomial for random effects and constant residual variance at all ages.

The comparison of the fitting residue model between the six classes of variance residual heterogeneity showed increased Log L significantly (P<0,01) for the likelihood ratio test (LRT) as the number of heterogeneous classes increased. All Heterogeneous Variance Models had better fit compared with the Homogeneous models (Class 1), given the lower values of BIC and AIC. The values of Log L, BIC and AIC changed considerably and significantly (P<0.01) for LRT up to Class 6 model, thus suggesting this to be the best fit.

Sarmento (2007), while studying the growth rate curve of Santa Inês sheep, found that the Homogeneous Residual Variance model was inadequate, while the Residual Variance model in five classes was a good fit for the data. Other authors also emphasized the need to consider Heterogeneous Residual Variance models to assess growth rate curves (Meyer, 2000; Sousa et al., 2008), for milk production on the control day (El faro & Albuquerque, 2003; Sarmento et al., 2008) and sensitivity studies of genotypes in meat quails (Felipe, 2010).

The random regression coefficient of the intercept (b₀) displays the largest variance of all random effects. For all other coefficients, variance decreased as the function order increased (Table 2). The correlations between intercept and linear coefficient were positive and high for all random effects. The correlations between intercept and the quadratic and cubic coefficients were always negative. For other coefficients, the correlations changed between positive and negative, with small values close to one.

The intercept justifies the high variance for the random effects in the model (Table 2). Despite being the model with the most parameters (model 66) and the best fit according the criteria adopted in the study (LRT, BIC and AIC), this result demonstrates that the regression coefficients of higher order tend to add less variance to the random effects.

Additive genetic variance increased up to 28 days of age and decreased slightly up to 42 days (Table 3). It should be noted that animal permanent environmental variance increased with time and dropped slightly in the last week.

Similar behavior was reported by Bonafé (2008), whose estimates considering heterogeneous variance ranged from 0.33 to 290.30 for the UFV1 meat quail and 0.33 to 288.30 for the lineage UFV2. Similar results were also reported by Akbas et al. (2004), while studying meat quails considering homogeneous residual variance, in which permanent environment variance also increased with age.

The phenotypic and animal permanent environmental variances displayed similar behavior, increased in the first week, decreased from 7 to 14 days and increased again until 35 days of age due to scale effect, since the weights increased with age. However, animal permanent environmental variance displayed lower values (Figure 1).

Similar results were reported by Akbas et al. (2004), Dionello et al. (2008) and Bonafé et al. (2011), while studying meat quails, in which both variances also increased with age.

Additive genetic variance increased from birth to 35 days of age and dropped slightly in the last week. The residual variance also decreased at 21 and 35 days of age; however, it increased up to 42 days of age. The drop of additive genetic variance at 42 days of age may be related to animal behavior patterns, since they become more agitated in the last rearing stages, females reach sexual maturity and males develop competitive behavior to determine patterns of social dominance, such as competition for territory, food and water. These factors may lead to an increase in residual variance that possibly reflects on the declining additive genetic variance.

Heritability decreased along the growth rate curve, ranging from 0.51 (1 day) to 0.16 (42 days of age). These results contradict Bonafé et al. (2011), who reported that heritability decreased up to 14 days, increased from 14 to 35 days and decreased again after 35 days (Figure 2).

According to Akbas et al. (2004), the heritability estimated for body weight of meat quails increased until the 5th week (0.61), with a slight drop during the last week (0.44), but remained high. Dionello et al. (2008) reported increasing heritability from birth (0.01) to 42 days of age (0.50) for the EV1 lineage, but low values (0.01 to 0.10) for the EV2 lineage.

Silva et al. (2010), while studying meat quails using multiple trait model, also reported that estimates decreased as age increased, varying from 0.43 at 28 days to 0.27 at 147 days of age.

The downward trend in heritability from 35 to 42 days may be related to the increase of the residual variance in the last two weeks, in which caused a steeper decrease of the additive genetic variance.

Estimates of genetic correlations (Table 4) between weights at different ages were high and positive, except for the correlation between birth weight and life weight, thus suggesting that selection should not be based on birth weight. From 7 days of age, all genetic correlations between body weights were high and positive, indicating that selection can start at this age. The highest genetic correlation (0.972) was observed between 21 and 28 days of age, suggesting that selection at 28 days of age can be effective, considering that correlation between 28 and 35 days of age (selection age) was also high (0.965). These results are consistent with those reported by Dionello et al. (2008) and Bonafé et al. (2011), who also reported high genetic correlation from 7 days of age.

The permanent environmental correlations were, in general, similar to the genetic correlations, which highlights the importance of strong permanent effect on body weight, indicating that environmental factors may influence body weight at all ages.

The genetic trend (Figure 3) displayed an increasing pattern over generations (1.2586 g/generation) for the genetic value for body weight at 35 days of age.

A fact that may explain the low breeding value of the second generation is the high number of birds that were born from the five hatchings. As a consequence, the higher number of birds per lot accentuated the competitive behavior among them. With increasing age, the birds became more agitated, tending to establish hierarchy for food and territory, which increased the residual variance, causing the breeding value of the animals to drop. This shows that animals raised out of an ideal condition can express their genetic potential differently; it would be a classic case of genotype × environment interaction.

Conclusions

The use of sixth order Legendre polynomial along with residual variance divided into six periods was the best fit for the growth rate curve of female LF1 meat quails, and they should be considered for breeding evaluation processes by random regression models. Heritability estimates decreased along the growth rate curves; therefore, genetic gains are higher when younger birds are selected. The genetic trend of breeding values increased over generations, thus symbolizing genetic progress of the studied population.

Received September 22, 2011 and accepted April 18, 2012.

Corresponding author: [e-mail]

• AKBAS, Y.; TAKMAN, Ç.; YAYLAC, E. Genetic parameters for quail body weights using a random regression model. South African Journal of Animal Science, v.34, n.2, p.104-109, 2004.
• BONAFÉ, C.M. Avaliação do crescimento de codornas de corte utilizando modelos de regressão aleatória 2008. 58f. Dissertação (Mestrado em Genética e Melhoramento) - Universidade Federal de Viçosa, Viçosa, MG.
• BONAFÉ, C.M.; TORRES, R.A.T.; SARMENTO, J.L.R.S. et al. Modelos de regressão aleatória para descrição da curva de crescimento de codornas de corte. Revista Brasileira de Zootecnia, v.40, n.4, p.765-771, 2011.
• CORRÊA, G.S.S.; SILVA, M.A.; CORRÊA, A.B. et al. Exigências em proteína bruta para codornas de corte EV1 em crescimento. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.59, n.5, p.1278-1286, 2007.
• DIONELLO, N.J.L.; CORREA, G.S.S.; SILVA, M.A. et al. Estimativas da trajetória genética do crescimento de codornas de corte utilizando modelos de regressão aleatória. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.60, n.2, p.454-460, 2008.
• EL FARO, L.; ALBUQUERQUE, L.G. Utilização de modelos de regressão aleatória para produção de leite no dia do controle, com diferentes estruturas de variâncias residuais. Revista Brasileira de Zootecnia, v.32, n.5, p.1104-1113, 2003.
• FELIPE, V.P.S. Estudo da interação genótipo-ambiente em codornas de corte 2010. 59f. Dissertação (Mestrado em Genética e Melhoramento Animal) - Universidade Federal de Minas Gerais, Belo Horizonte.
• FRIDRICH, A.B.; VALENTE, B.D.; FELIPE SILVA, A.S. et al. Exigência de proteína bruta para codornas européias no período de crescimento. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.57, n.2, p.261-265, 2005.
• MEYER, K. "WOMBAT" - Digging deep for quantitative genetic analyses by restricted maximum likelihood. In: WORLD CONGRESS ON GENETIC APPLIED TO LIVESTOCK PRODUCTION, 8., 2006, Belo Horizonte. Proceedings... Belo Horizonte, [2006] (CD-ROM).
• MEYER, K. Random regressions to model phenotypic variation in monthly weights of Australian beef cows. Livestock Production Science, v.65, n.1, p.19-38, 2000.
• NUNEZ-ANTÓN, V.N.; ZIMMERMAN, D.L. Modelling nonstationary longitudinal data. Biometrics, v.56, p.699-705, 2000.
• OLIVEIRA, E.G.; ALMEIDA, M.I.M.; MENDES, A.A. et al. Desempenho produtivo de codornas de ambos os sexos para corte alimentadas com dietas com quatro níveis protéicos. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.7, n.2, p.75-80, 2002.
• SARMENTO, J.L.R. Modelos de regressão aleatória para avaliação genética da curva de crescimento de ovinos da raça Santa Inês 2007. 101f. Tese (Doutorado em Genética e Melhoramento) - Universidade Federal de Viçosa, Viçosa, MG.
• SARMENTO, J.L.R.; ALBUQUERQUE, L.G.; TORRES, R.A. et al. Comparação de modelos de regressão aleatória para estimação de parâmetros genéticos em caprinos leiteiros. Revista Brasileira de Zootecnia, v.37, n.10, p.1788-1796, 2008.
• SILVA, L.P.; LEITE, C.D.S.; SOUSA, F.M. et al. Parâmetros genéticos para pesos corporais de matrizes de codornas de corte. In: SIMPÓSIO BRASILEIRO DE MELHORAMENTO ANIMAL, 8., 2010, Maringá. Anais... Maringá, [2010]. (CD-ROM).
• SOUSA, J.E.R.; SILVA, M.A.; SARMENTO. et al. Homogeneidade e heterogeneidade de variância residual em modelos de regressão aleatória sobre o crescimento de caprinos Anglo-Nubianos. Pesquisa Agropecuária Brasileira, v.43, n.12, p.1725-1732, 2008.
• THOLON, P.; QUEIROZ, S.A.; Utilização de diferentes estruturas de variância residual em modelos de regressão aleatória para descrição da curva de crescimento de perdizes (Rhynchotus rufescens) criadas em cativeiro. Revista Caatinga, v.21, n.2, p.37-47, 2008.
• VELOSO, R.C.; DRUMOND, E.S.C.; VASCONCELOS, R.C. et al. Proteína bruta e energia metabolizável para codornas de corte: 1 - Fase inicial. In: SIMPÓSIO INTERNACIONAL, 4.; CONGRESSO BRASILEIRO DE CODORNAS DE CORTE, 3., 2010, Lavras. Anais... Lavras, [2010]. (CD-ROM).
• WOLFINGER, R. Covariance structure selection in general mixed models. Community of Statistics - Simulation, v.22, n.4, p.1079-1106, 1993.

1

Supported by CNPq, CAPES and FAPEMIG.

Publication Dates

• Publication in this collection
  02 Oct 2012
• Date of issue
  Sept 2012

History

• Received
  22 Sept 2011
• Accepted
  18 Apr 2012