Avaliação genética de bovinos da raça holandesa usando a produção de leite no dia do controle

Utilizaram-se 153.963 controles mensais de produção de leite e 13.273 primeiras lactações de vacas da raça Holandesa, com partos entre 1989 a 1998, com o objetivo de estimar parâmetros genéticos, fenotípicos e de meio ambiente para produção de leite no dia do controle (PLDC) e produção até 305 dias de lactação (P305), bem como verificar a conveniência de se utilizar a PLDC em avaliações genéticas, em substituição à P305. Foram utilizados quatro modelos (modelo animal). Em dois, modelos 1 e 2, consideraram-se os controles mensais de produção, coletados ao longo da lactação, como medidas repetidas de um animal, diferenciados segundo o critério de formação de grupos de animais contemporâneos. No modelo 1 (PLDCM01), as produções dos animais foram agrupadas segundo o rebanho-ano-estação de controle de produção, enquanto, no modelo 2 (PLDCM02), o agrupamento considerou o rebanho-ano-estação de parto. No modelo 3, analisaram-se os controles mensais de produção, como características individuais (C01 até C10); e no modelo 4, analisou-se a tradicional P305. As análises foram realizadas utilizando-se o método da Máxima Verossimilhança Restrita, por meio do sistema MTDFREML. As estimativas de herdabilidade para PLDC, com o uso do modelo 1, modelo 2 e para P305, foram de 0,27, 0,15 e 0,25, respectivamente. As herdabilidades para os controles mensais variaram de 0,11±0,02 (C01) a 0,21±0,03 (C08), e os maiores valores ocorreram a partir do quarto mês. A correlação de ordem entre os valores genéticos obtidos, para P305 e para PLDC (modelo 1), foi de 0,62, para touros, e de 0,78, para vacas. Concluiu-se que é viável a utilização da PLDC em estudos envolvendo a produção de leite e que os controles do meio da lactação, se usados para seleção, podem apresentar vantagens em relação à P305.

Palavras-chave: controle mensal, método REML, modelo animal, produção de leite, raça Holandesa

153,963 test day milk yield records and 13,273 first lactations of Holstein cows calving between 1989 and 1998, were used with the objective of estimating genetic, phenotypic and environmental parameters for test day milk yield (PLDC) and 305 day milk yield (P305) and to study the convenience of using test day yields in genetic evaluations to replace P305. Four models were used. Models 1 and 2 differed according contemporary grouping and monthly milk records were considered as repeated measures. In model 1 (PLDCM01) records were grouped by herd-year-season of test day yield and in model 2 (PLDCM02) by herd-year-season of calving. In a third (model 3), monthly yield records were analyzed as individual traits (C01 to C10); and the fourth (model 4) was the traditional 305-day model. Restricted Maximum Likelihood methodology was used with the MTDFREML system. The estimates of heritability for PLDC, using model 1, model 2 and for P305 were 0.27, 0.15 and 0.25, respectively. Heritabilities for monthly milk records ranged from 0.11+0.02 (C01) to 0.21+0.03 (C08), with the largest values occurring beginning in the fourth month. Genetic correlation estimates between monthly records and P305 ranged from 0.76 to 1.00, with the highest correlations occurring in the middle lactation. It was concluded that using test day milk yield is promising and that selecting for middle lactation records could have advantages over 305-day milk yield.

Keywords: animal model, Holstein, milk yield, REML methodology, test day yields

A medida padrão de produção de leite mais usada nas avaliações genéticas de vacas e touros é a produção até 305 dias de lactação (P305). A produção de cada mês, calculada com base na produção de leite no dia do controle (PLDC), é acumulada com o total de meses anteriores, para estimação da produção até 305 dias. Mesmo em países de pecuária leiteira mais tecnificada, até recentemente, uma vez estimada esta produção, os controles mensais eram descartados, em virtude da falta de recursos computacionais para armazenamento (Ptak & Schaeffer, 1993; Wiggans & Goddard, 1997). Com os recentes avanços na área de informática, este descarte não é mais necessário. Os controles podem, então, ser usados em estudos da curva de lactação de fatores que afetam cada controle e para avaliações genéticas dos animais.

Atualmente, nas avaliações genéticas em gado de leite, modelos de repetibilidade, em que a PLDC é analisada como medida repetida, têm sido usados em substituição à modelos para P305 (Swalve, 1995a,b; Schaeffer et al., 2000). Países como os Estados Unidos, Canadá, Austrália, Nova Zelândia, entre outros, já utilizam modelos para a PLDC, em avaliações de características produtivas e contagem de células somáticas (Meyer et al., 1989; Reents et. al., 1995; Swalve, 2000). Têm-se observado muitas vantagens de tal mudança. As herdabilidades estimadas para a PLDC têm sido similares às encontradas para P305 (Meyer et al., 1989; Danell, 1982; Ptak & Schaeffer, 1993). Além disso, os modelos para PLDC podem considerar a forma da curva de lactação, o efeito do número de dias em produção, o efeito peculiar do dia do controle para todas as vacas e os efeitos específicos para cada vaca no dia do controle, tais como número de ordenhas, duração do período seco, duração do período de serviço anterior e corrente, prenhez ou doença. As vacas podem ser avaliadas antes que suas lactações sejam encerradas, sem necessidade da utilização de fatores de ajustamento ou projeção da produção. Também, os touros podem ser avaliados com maior confiabilidade, uma vez que poderá haver maior número de filhas com produção (Ptak & Schaeffer, 1993; Swalve, 1995b). Com a PLDC é possível avaliar os animais mais cedo, permitindo a seleção dos melhores para a reprodução no início da vida útil, possibilitando redução no intervalo de gerações.

O uso da PLDC permite também flexibilidade nos esquemas de controle leiteiro, possibilitando redução do número de controles por lactação (Swalve, 2000). Para tanto, são necessários estudos sobre a utilização de menor número de controles nas avaliações genéticas em gado de leite.

O agrupamento rebanho-ano-estação de parto (RAE) é muito usado em avaliações genéticas em gado de leite para P305, para comparar as produções de animais sujeitos aos mesmos efeitos de meio ambiente. Nos estudos sobre a PLDC, o agrupamento poderá ser feito também, considerando a estação em que o controle de produção de leite é realizado, em vez da estação de parto. A substituição de RAE por rebanho-ano-estação de controle (RDC), em modelos para PLDC, possibilita a comparação de animais que têm produção controlada no mesmo período, enfatizando os efeitos de ambiente na data de controle (Swalve, 1995b). Além disso, na estimação de componentes de variância, têm-se observado redução da variância residual, aumento da variância genética e conseqüente aumento da herdabilidade para produção de leite, quando classes de RDC são usadas nos modelos para PLDC.

Assim, objetivou-se com este trabalho definir um critério de formação de grupos de animais contemporâneos em modelo de repetibilidade para PLDC, e estimar parâmetros genéticos, fenotípicos e de meio ambiente para PLDC e P305, e ainda verificar a conveniência de se utilizar a PLDC em avaliações genéticas, em substituição à P305.

Foram analisados dados de produção de leite no dia do controle, coletados ao longo da lactação, como medida repetida de um mesmo animal e como características individuais e, também, a produção até 305 dias de lactação. Os dados usados nas análises foram provenientes do Serviço de Controle Leiteiro da Associação de Criadores de Gado Holandês de Minas Gerais.

As lactações nas quais ocorreram doenças, abortos, mortes ou vendas de vacas, as lactações em andamento e também aquelas encerradas antes de 305 dias sem registro da data ou causa de secagem, não foram usadas. Também foram eliminados controles inferiores a 3 kg ou superiores a 45 kg e vacas com idade ao parto menor que 20 ou maior que 48 meses. Após as eliminações, foram calculadas as produções até 305 dias de lactação ajustadas para duas ordenhas, conforme Teixeira (1998) e as produções no primeiro e no último controle, de acordo com o método oficial regulamentado pelo Ministério da Agricultura e do Abastecimento (Brasil, 1986). Não foram calculadas as produções até 305 dias de lactação de vacas com menos de três controles mensais de produção ou com intervalo de controles maior que 75 dias, antes de completar 305 dias de lactação. Foram ainda descartadas as produções de leite até 305 dias de lactação inferiores a 1.200 kg ou superiores a 13.000 kg.

Após estas restrições, restaram 153.963 controles mensais de produção e 13.273 primeiras lactações de vacas com partos registrados no período de 1989 a 1998.

Os controles mensais de produção de leite foram divididos em intervalos de, aproximadamente, 30 dias. Uma vez que o controle leiteiro nem sempre é realizado em intervalos regulares e com o intuito de evitar eliminação dos mesmos, considerou-se um período de, no máximo, 15 dias de atraso nas realizações das pesagens mensais. Assim, como primeiro controle, foram consideradas as pesagens efetuadas no período de quatro a 45 dias após o parto; o segundo, de 30 a 75; o terceiro, de 60 a 105; e assim sucessivamente. Quando havia mais de um controle no mesmo intervalo, um deles era transferido para o intervalo seguinte, desde que este estivesse dentro dos limites do próximo intervalo; caso contrário, um deles era eliminado aleatoriamente.

Foram utilizados quatro modelos nas análises. O modelo 1 (PLDCM01), para observações repetidas, também usado por Ptak & Schaeffer (1993), foi:

em que y_ijkl = produção de leite no dia do controle; RDC_i = efeito do rebanho-ano-estação i em que o controle foi realizado (estações: águas, de outubro a março e seca, de abril a setembro); CG_j = efeito da composição genética j (animais puros de origem e animais de composição genética igual ou superior a 31/32 Holandês); b₁ = coeficiente de regressão linear da característica y_ijkl, sobre a idade da vaca no parto; b₂ a b₅ = coeficientes de regressão da característica y_ijkl, sobre o número de dias em lactação (DIM); X₁ = idade da vaca no parto, em dias; X₂ = DIM/c, em que c = 305 é uma constante positiva utilizada na redução da magnitude das covariâncias; X₃ = (DIM/c)²; X₄ = ln(c/DIM); X₅ = (ln(c/DIM))²; ak = efeito aleatório do animal k; pe_k = efeito aleatório de ambiente permanente sobre a vaca k, durante a lactação; e_ijkl = erro aleatório associado a cada observação.

Admitiu-se que o efeito de ambiente permanente e o efeito do erro aleatório são independentemente e identicamente distribuídos, com média zero e variância e , respectivamente, e não-correlacionados com o efeito de animal.

Nesta análise, foram considerados somente os dados de produção de vacas com, no mínimo, três controles coletados ao longo da lactação. Além disso, os touros deveriam ter, pelo menos, cinco filhas e cada classe de RDC deveria possuir, no mínimo, três controles de filhas de, pelo menos, dois touros.

O modelo 2 (PLDCM02) é semelhante ao modelo 1, exceto que o efeito fixo de RDC foi substituído por RAE (rebanho-ano-estação de parto). Também, as restrições para os dados foram as mesmas para o modelo 1.

O modelo 3, apresentado em seguida, foi usado nas análises dos controles mensais de produção de leite, como características individuais.

em que y_ijk = controle mensal de produção de leite (C01 até C10); b₁ = coeficiente de regressão linear da característica y_ijk, sobre a idade da vaca no parto; b₂ = coeficiente de regressão linear da característica y_ijk, sobre o intervalo do parto ao primeiro controle.

Além das restrições quanto ao número de filhas por touro e de observações por classe de RAE, descritas nos modelos 1 e 2, foram excluídos os animais sem o primeiro controle.

O modelo 4, tradicionalmente usado para analisar a produção até 305 dias de lactação, foi:

em que y_ijk = produção de leite até 305 dias de lactação; b₁ e b₂ = coeficientes de regressão linear e quadrático da característica y_ijk, sobre idade da vaca no parto.

Para esta análise também, os touros deveriam ter, pelo menos, cinco filhas e cada classe de RAE deveria possuir, no mínimo, três registros de produção de filhas de, pelo menos, dois touros.

Além das análises de característica única, foram feitas as de características múltiplas (análise com duas características), entre os controles mensais de produção de leite (C01 até C10) e a produção até 305 dias.

Em todas as análises foi utilizado um arquivo de "pedigree", contendo identificação de animal, pai e mãe, de 46.369 indivíduos, resultando em 62.885 animais diferentes na matriz de parentesco.

Os componentes de variância e covariância, necessários à estimação dos parâmetros genéticos e à predição dos valores genéticos dos animais, para as características em estudo, foram obtidos pelo método da Máxima Verossimilhança Restrita (REML), utilizando-se o aplicativo MTDFREML (Boldman et al., 1995).

A confiabilidade dos valores genéticos, correlação de ordem entre eles e a coincidência na seleção e descarte de animais foram usadas como critério para comparação das alternativas de avaliação genética. Análises de característica única foram utilizadas para obtenção dos valores genéticos e confiabilidades. A base de comparação foram as predições dos valores genéticos para P305. Os coeficientes de correlação de ordem de Spearman foram obtidos por meio do sistema SAS (1990).