Buffalos milk yield analysis using random regression models

Data comprising 1,719 milk yield records from 357 females (predominantly Murrah breed), daughters of 110 sires, with births from 1974 to 2004, obtained from the Programa de Melhoramento Genetic de Bubalinos (PROMEBUL) and from records of EMBRAPA Amazonia Oriental - EAO herd, located in Belem, Para, Brazil, were used to compare random regression models for estimating variance components and predicting breeding values of the sires. The data were analyzed by different models using the Legendre's polynomial functions from second to fourth orders. The random regression models included the effects of herd-year, month of parity date of the control; regression coefficients for age of females (in order to describe the fixed part of the lactation curve) and random regression coefficients related to the direct genetic and permanent environment effects. The comparisons among the models were based on the Akaike Infromation Criterion. The random effects regression model using third order Legendre's polynomials with four classes of the environmental effect were the one that best described the additive genetic variation in milk yield. The heritability estimates varied from 0.08 to 0.40. The genetic correlation between milk yields in younger ages was close to the unit, but in older ages it was low.

Estimates of direct and maternal genetic effects for weights from birth to 600 days of age in Nelore cattle

Estimates of direct and maternal variance and heritability for weights at each week (up to 280 days of age) and month of age (up to 600 days of age) in Zebu cattle are presented. More than one million records on 200 000 animals, weighed every 90 days from birth to 2 years of age, were available. Data were split according to week (data sets 1) or month (data sets 2) of age at recording, creating 54 and 21 data sets, respectively. The model of analysis included contemporary groups as fixed effects, and age of dam (linear and quadratic) and age of calf (linear) effects as covariables. Random effects fitted were additive direct and maternal genetic effects, and maternal permanent environmental effect. Direct heritability estimates decreased from 0.28 at birth, to 0.12-0.13 at about 150 days of age, stayed more or less constant at 0.14-0.16 until 270 days of age and increased with age after that, up to 0.25-0.26. Maternal heritability estimates increased from birth (0.01) to a peak of 0.14 for data sets 1 and 0.07-0.08 for data sets 2 at about 180-210 days of age, before decreasing slowly to 0.07 and 0.05, respectively, at 300 days, and then rapidly diminished after 300 days of age. Permanent environmental effects were 1.5 to four times higher than genetic maternal effects and showed a similar trend.

Comparison of random regression models to estimate genetic parameters for milk production in Guzerat (Bos indicus) cows

Random regression models have been widely used to estimate genetic parameters that influence milk production in Bos taurus breeds, and more recently in B. indicus breeds. With the aim of finding appropriate random regression model to analyze milk yield, different parametric functions were compared, applied to 20,524 test-day milk yield records of 2816 first-lactation Guzerat (B. indicus) cows in Brazilian herds. The records were analyzed by random regression models whose random effects were additive genetic, permanent environmental and residual, and whose fixed effects were contemporary group, the covariable cow age at calving (linear and quadratic effects), and the herd lactation curve. The additive genetic and permanent environmental effects were modeled by the Wilmink function, a modified Wilmink function (with the second term divided by 100), a function that combined third-order Legendre polynomials with the last term of the Wilmink function, and the Ali and Schaeffer function. The residual variances were modeled by means of 1, 4, 6, or 10 heterogeneous classes, with the exception of the last term of the Wilmink function, for which there were 1, from 0.20 to 0.33. Genetic correlations between adjacent records were high values (0.83-0.99), but they declined when the interval between the test-day records increased, and were negative between the first and last records. The model employing the Ali and Schaeffer function with six residual variance classes was the most suitable for fitting the data. © FUNPEC-RP.

Genetic evaluation using multi-trait and random regression models in Simmental beef cattle

The Brazilian Association of Simmental and Simbrasil Cattle Farmers provided 29,510 records from 10,659 Simmental beef cattle; these were used to estimate (co)variance components and genetic parameters for weights in the growth trajectory, based on multi-trait (MTM) and random regression models (RRM). The (co)variance components and genetic parameters were estimated by restricted maximum likelihood. In the MTM analysis, the likelihood ratio test was used to determine the significance of random effects included in the model and to define the most appropriate model. All random effects were significant and included in the final model. In the RRM analysis, different adjustments of polynomial orders were compared for 5 different criteria to choose the best fit model. An RRM of third order for the direct additive genetic, direct permanent environmental, maternal additive genetic, and maternal permanent environment effects was sufficient to model variance structures in the growth trajectory of the animals. The (co)variance components were generally similar in MTM and RRM. Direct heritabilities of MTM were slightly lower than RRM and varied from 0.04 to 0.42 and 0.16 to 0.45, respectively. Additive direct correlations were mostly positive and of high magnitude, being highest at closest ages. Considering the results and that pre-adjustment of the weights to standard ages is not required, RRM is recommended for genetic evaluation of Simmental beef cattle in Brazil. ©FUNPEC-RP.

Multiple-trait random regression models for the estimation of genetic parameters for milk, fat, and protein yield in buffaloes

In this study, genetic parameters for test-day milk, fat, and protein yield were estimated for the first lactation. The data analyzed consisted of 1,433 first lactations of Murrah buffaloes, daughters of 113 sires from 12 herds in the state of São Paulo, Brazil, with calvings from 1985 to 2007. Ten-month classes of lactation days were considered for the test-day yields. The (co)variance components for the 3 traits were estimated using the regression analyses by Bayesian inference applying an animal model by Gibbs sampling. The contemporary groups were defined as herd-year-month of the test day. In the model, the random effects were additive genetic, permanent environment, and residual. The fixed effects were contemporary group and number of milkings (1 or 2), the linear and quadratic effects of the covariable age of the buffalo at calving, as well as the mean lactation curve of the population, which was modeled by orthogonal Legendre polynomials of fourth order. The random effects for the traits studied were modeled by Legendre polynomials of third and fourth order for additive genetic and permanent environment, respectively, the residual variances were modeled considering 4 residual classes. The heritability estimates for the traits were moderate (from 0.21-0.38), with higher estimates in the intermediate lactation phase. The genetic correlation estimates within and among the traits varied from 0.05 to 0.99. The results indicate that the selection for any trait test day will result in an indirect genetic gain for milk, fat, and protein yield in all periods of the lactation curve. The accuracy associated with estimated breeding values obtained using multi-trait random regression was slightly higher (around 8%) compared with single-trait random regression. This difference may be because to the greater amount of information available per animal. © 2013 American Dairy Science Association.

Genetic parameters for body weight in buffaloes (Bubalus bubalis) in Colombia using random regression models

The objective of this research was to estimate (co) variance functions and genetic parameters for body weight in Colombian buffalo populations using random regression models with Legendre polynomials. Data consisted of 34,738 weight records from birth to 900 days of age from 7815 buffaloes. Fixed effects in the model were contemporary group and parity order of the mother. Random effects were direct and maternal additive genetic, as well as animal and maternal permanent environmental effects. A cubic orthogonal Legendre polynomial was used to model the mean curve of the population. Eleven models with first to sixth order polynomials were used to describe additive genetic direct and maternal effects, and animal and maternal permanent environmental effects. The residual was modeled considering five variance classes. The best model included fourth and sixth order polynomials for direct additive genetic and animal permanent environmental effects, respectively, and third-order polynomials for maternal genetic and maternal permanent environmental effects. The direct heritability increased from birth until 120 days of age (0.32 +/- 0.05), decreasing thereafter until one year of age (0.18 +/- 0.04) and increased again, reaching 0.39 +/- 0.09, at the end of the evaluated period. The highest maternal heritability estimates (0.11 +/- 0.05), were obtained for weights around weaning age (weaning age range is between 8 and 9.5 months). Maternal genetic and maternal permanent environmental variances increased from birth until about one year of age, decreasing at later ages. Direct genetic correlations ranged from moderate (0.60 +/- 0.060) to high (0.99 +/- 0.001), maternal genetic correlations showed a similar range (0.41 +/- 0.401 and 0.99 +/- 0.003), and all of them decreased as time between weighings increased. Direct genetic correlations suggested that selecting buffalos for heavier weights at any age would increase weights from birth through 900 days of age. However, higher heritabilities for direct genetic weights effects after 600 days of age suggested that selection for these effects would be more effective if done during this age period. A greater response to selection for maternal ability would be expected if selection used maternal genetic predictions for weights near weaning. (C) 2013 Elsevier B.V. All rights reserved.

Random regression models to estimate genetic parameters for milk production of Guzerat cows using orthogonal Legendre polynomials

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Random regression models using different functions to model milk flow in dairy cows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic associations between weight at maturity and maturation rate with ages and weights at first and second calving in Canchim beef cattle

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Análise genético-quantitativa de pesos aos 8, 12, 18 e 24 meses de idade em um rebanho da raça Guzerá

O objetivo deste trabalho foi obter estimativas de herdabilidades, correlações genéticas, tendências genéticas e correlações de rank dos touros, para os pesos aos 8, 12, 18 e 24 meses de idade, no rebanho Guzerá do Campus da UNESP, Ilha Solteira, SP. As herdabilidades e os valores genéticos dos animais foram estimados por modelo animal, usando o programa computacional MTDFREML. As correlações genéticas (Pearson) e de rank dos touros (Spearman) foram obtidas pelo procedimento CORR do SAS, utilizando os valores genéticos dos animais, enquanto as tendências genéticas foram calculadas pelo procedimento REG do SAS, utilizando a média anual dos referidos valores genéticos. O modelo estatístico para a análise genética incluiu os efeitos fixos de grupo contemporâneo e a covariável idade da vaca ao parto (efeitos linear e quadrático), e os efeitos aleatórios genético aditivo direto, genético aditivo materno, de ambiente permanente da vaca e residual. As estimativas obtidas para a herdabilidade direta foram 0,14; 0,08; 0,08 e 0,13 e para materna, 0,01; 0,02; 0,02 e 0,05, respectivamente, para os pesos aos 8, 12, 18 e 24 meses de idade. As estimativas dos coeficientes de correlação genética foram positivas e de alta magnitude entre todos os pesos estudados. As tendências genéticas anuais foram baixas e significativas apenas para os pesos aos 8 e 18 meses de idade. As correlações de rank dos touros foram moderadas e significativas, implicando em razoável manutenção de posição de classificação dos touros, quando se comparam, dois a dois, os pesos estudados.

Efeitos genéticos e ambientais sobre o perímetro escrotal de animais da raça Caracu

Neste trabalho foram avaliados os efeitos de idade (IDS) e peso aos 15 meses (P15) sobre medidas de perímetro escrotal ao sobreano de 1.892 machos da raça Caracu e estimada a herdabilidade desta característica. Utilizaram-se dois modelos de análise: em um dos modelos, foram incluídos o efeito de grupo de contemporâneos (GC) e, como covariável, a idade ao sobreano (efeitos linear e quadrático); e, no outro, o efeito de GC e, como covariável, o peso ao sobreano (efeitos linear e quadrático). A idade não teve efeito significativo no perímetro escrotal, o que está relacionado à homogeneidade dos animais nos grupos de contemporâneos, enquanto os efeitos linear e quadrático de P15 foram significativos sobre o perímetro escrotal, indicando grande influência do peso sobre a variação desta característica. A herdabilidade do perímetro escrotal foi estimada pelo método Bayesiano utilizando-se um modelo animal. O modelo incluiu os efeitos de GC e P15 (linear e quadrático) e os efeitos genéticos aditivo direto e residual. O valor médio estimado de herdabilidade do PE ao sobreano foi 0,38 e comprova que esta característica pode ser utilizada como critério de seleção para precocidade sexual em programas de melhoramento genético de animais da raça Caracu.

Análise da variabilidade genética aditiva de características de crescimento na raça Nelore

Foram utilizados dados de cinqüenta e um rebanhos participantes do Programa de Melhoramento Genético da Raça Nelore (PMGRN), distribuídos nos estados de Goiás (GO), Mato Grosso do Sul (MS), Mato Grosso (MT), Minas Gerais (MG), São Paulo (SP), Maranhão (MA) e Bahia (BA). Foram obtidas estimativas de parâmetros genéticos para os pesos padronizados aos 120 (P120), 455 (P455) e 550 (P550) dias de idade. Análises unicaráter e bicaráter foram realizadas por modelo animal usando o aplicativo MTDFREML. Para P120 foi utilizado um modelo que incluiu como efeitos fixos, grupo de contemporâneos e classe de idade da vaca ao parto, e como aleatórios, os efeitos genéticos direto, materno e de ambiente permanente da vaca. Para P455 e P550, o modelo utilizado incluiu os mesmos efeitos fixos e o efeito genético direto do animal. ANas análises unicaráter, as estimativas de herdabilidade direta foram 0,29, 0,51 e 0,47 para P120, P455 e P550, respectivamente. Nas análises bicaráter, observaram-se coeficientes de herdabilidade direta de 0,50 e 0,58 para P120, 0,50 e 0,53 para P455 e 0,44 e 0,49 para P550. As correlações genéticas estimadas entre P120 e P455, P120 e P550 e P455 e P550, foram 0,92, 0,93 e 0,96, respectivamente. As estimativas de herdabilidade obtidas para P455 e as correlações genéticas deste peso com P120 e P550 sugerem que a avaliação genética pode ser feita aos 15 meses de idade em substituição aos 18 meses.

Herdabilidade e correlações quanto a peso, perímetro escrotal e escores visuais à desmama, em bovinos Canchim

O objetivo deste trabalho foi estabelecer o modelo mais adequado para avaliação genética de bovinos Canchim e estimar os parâmetros genéticos de características produtivas à desmama. Foram utilizados dados de: 12.103 animais, quanto ao peso (PD); 5.278, quanto ao perímetro escrotal (PE); 8.343, quanto ao escore visual da conformação frigorífica (CF); 9.111, quanto ao escore de umbigo (UM); e 7.986, quanto ao escore de pelame (PEL). Os modelos estatísticos incluíram os efeitos fixos e os efeitos aleatórios genéticos aditivos direto, materno e de ambiente permanente materno, em diferentes combinações. As análises foram feitas pelo método da máxima verossimilhança restrita livre de derivadas. O modelo completo foi o mais adequado para PD, PE, CF e UM, enquanto o modelo com uso apenas dos efeitos genéticos aditivos direto e materno foi o mais adequado para PEL. As estimativas de herdabilidade direta foram 0,17, 0,13, 0,20, 0,18, e 0,52 para PD, PE, CF, UM e PEL, respectivamente, o que indica a possibilidade de se obter progresso genético por meio da seleção para essas características, principalmente para PEL. As correlações genéticas aditivas diretas entre as características variaram de -0,16 a 0,61. As correlações entre PD e PE e entre PD e CF indicam que a seleção para PD deve proporcionar ganho genético em PE e CF

Modelagem e estimação de parâmetros genéticos e fenotípicos para pesos do nascimento à seleção (378 dias) de machos Nelore

Genetic parameters for weights (17, 942 records), obtained in intervals of 60 days, from the birth to selection (378 days of age), of 2,582 males of the Nellore breed was estimated in univariate analyses by the Maximum Restricted Likelihood method. The models of analysis models included the fixed effects of contemporary groups, month of birth, mother age and age when the weights were collected as covariate. Three models differing in random effects were tested: the model 1 (M1) was adjusted for the direct and maternal addictive genetic effects and maternal permanent environment; in model 2 (M2) the maternal genetic effect was excluded; and the model 3 (M3) was only adjusted for the direct addictive genetic effect. The test of likelihood (LRT) detected significant differences, for all the ages, of M2 and M1 in relation to the simple model (M3), evidencing the importance of the maternal effects. Except for the birth weight (0.40), low values (0.05 to 0.12) of h(2) were found for M1 and M2 until 8 months of age and, after this period, reasonable increase could be observed, reaching 0.28 until 13 months of age. The estimates of the total variance fraction, due to the effect of maternal permanent environment, were high and practically became unaffected between the models 1 and 2. Maternal effects, not necessarily decomposed (in genetic addictive and permanent environment), affected the Nellore males growth. Models that contemplate maternal effects, besides the genetic addictive direct effects, are more realistic to describe the trajectory of the variances in the initial phases of growth of Nelore male calves.

Seleção para peso pós-desmama em um rebanho Gir. 2. Estimativas de variâncias e parâmetros genéticos dos efeitos direto e materno para características de crescimento

Componentes de variância e parâmetros genéticos para características de crescimento foram estimados usando diferentes modelos em um rebanho da raça Gir. Utilizou-se o método da máxima verossimilhança restrita sob modelo animal univariado. Os modelos de análise incluíram os efeitos fixos de mês de nascimento, grupo contemporâneo e idade da vaca. Cinco modelos diferindo quanto aos efeitos aleatórios foram testados. Para todas as características da fase pré-desmama, o teste de razão de verossimilhança (LRT) indicou o modelo com efeito genético aditivo direto e efeitos maternos (genético e de ambiente permanente) como o de melhor ajuste. As estimativas de herdabilidade direta para peso ao nascer (PN), peso aos quatro meses corrigido para 120 dias (P120), peso à desmama corrigido para 210 dias (P210) e ganho diário na fase pré-desmama (GPRE) foram, respectivamente, 0,31± 0,07; 0,14 ± 0,06; 0,23 ± 0,07 e 0,22 ± 0,07. Para as características da fase pós-desmama, o modelo que forneceu o melhor ajuste aos dados incluiu apenas o efeito genético aditivo direto. As estimativas de herdabilidade direta para peso de machos ao final da prova de ganho de peso (P378), peso de fêmeas corrigido para 550 dias (P550), ganho diário na prova de ganho de peso (G112), altura aos 378 dias em machos (AM) e altura aos 550 dias, em fêmeas (AF) foram, respectivamente: 0,45 ± 0,11; 0,29 ± 0,11; 0,37 ± 0,11; 0,79 ± 0,13 e 0,36 ± 0,0. Os efeitos maternos, tanto o genético quanto o de ambiente permanente, foram fontes de variação importantes para as características da fase pré-desmama, não sendo verificada influência desses efeitos sobre as características da fase pós-desmama.