Acessibilidade / Reportar erro

Parâmetros genéticos para as produções de leite no dia do controle e da primeira lactação de vacas da raça Holandesa

Genetic parameters for test day and lactation milk yields of first lactation Holstein cows

Resumos

Foram utilizados 263.390 registros de produção de leite no dia do controle (PDC) de 32.448 primeiras lactações de vacas da raça Holandesa com parto no período de 1991 a 2001, para estimar componentes de variância e parâmetros genéticos utilizando um modelo animal e a metodologia REML. Os dados de produção foram coletados pelo Serviço de Controle Leiteiro da Associação Brasileira de Criadores de Bovinos da Raça Holandesa e suas afiliadas estaduais e disponibilizados pela Embrapa Gado de Leite. O modelo para ajuste da produção de leite até 305 dias (P305) incluiu os efeitos fixos de rebanho-ano de parto, época do parto e idade da vaca ao parto, com termos linear e quadrático, e os efeitos aleatórios de animal e erro. Os mesmos efeitos foram incluídos no modelo para as produções de leite no dia do controle (PDC) tanto sob modelo uni- e bi-caráter como sob modelo de repetibilidade (MRS), exceto para o efeito fixo de grupo contemporâneo, definido por rebanho-ano-mês do controle. Alternativamente, ajustou-se um segundo modelo de repetibilidade (MRF), que além dos efeitos presentes no MRS, incluiu as covariáveis que descrevem a curva da lactação: dias em lactação (DEL)/305 e ln(305/DEL) com termos linear e quadrático. As estimativas de h² para as PDC, com MRS e MRF foram 0,43 e 0,30, respectivamente. As estimativas de herdabilidade (h²) para as PDC variaram de 0,22 (PDC1) a 0,36 (PDC4) com o modelo uni-caráter. Para o modelo bi-caráter, as estimativas variaram de 0,23 (PDC1) a 0,33 (PDC3 e PDC4). Em ambos os modelos, observaram-se valores inferiores de h² no início e no fim do período de lactação. A estimativa h² para a P305 com o modelo uni-caráter foi 0,27, enquanto, pelo modelo bi-caráter, estas estimativas variaram de 0,27 a 0,30. As correlações genéticas (r g) entre as PDC e a P305 foram altas, variando de 0,86 (PDC1 e P305) a 0,99 (PDC3 e P305). Estimativas de herdabilidade maiores para as PDC que para a P305 e a alta correlação com a mesma indicam potencial de uso das PDC nas avaliações genéticas de animais da raça Holandesa no Brasil. Embora predominantemente altas, as estimativas de r g entre as PDC não foram homogêneas (0,64-1,0); entretanto as maiores freqüências foram para valores próximos ou iguais a 1. Assim, modelos de regressão aleatória devem ser também avaliados para se concluir sobre a melhor utilização das PDC da raça Holandesa no Brasil.

análise bi-caráter; componentes de variância; modelo de repetibilidade; parâmetros genéticos; produções de leite do dia do controle


Covariance components for test day records and lactation milk yield using 263.390 records of 32.448 first lactation Holstein cows, were estimated using animal models by REML. Besides the lactation model, two alternative repeatability models (RM) were analyzed. Lactation model included fixed effects of herd-year-season and age of cow with linear and quadratic terms, and random effects of animal and error. The first model for test-day yield (RMF) included the same effects, but fixed effect of contemporary group, defined as herd-year-month of test. Alternatively another model for test-day yield (RMF) used a logarithmic polynomial sub-model for the shape of the lactation curve. Heritability for lactation yield (0.27) was smaller than those estimated by RMF and RMS, 0.30 and 0.43, respectively. Heritability estimates for univariate (0.22-0.36) and bivariate models (0.23-0.33) for test day milk yields were found to be smallest during early and late lactation. Heritability estimate for lactation milk yield when estimated by univariate model (0.27) was smaller than estimates obtained by bivariate models (0.27-0.30). Genetic correlations were higher between consecutive test days than between test days in the beginning and end of lactation. Larger heritability estimates for test day models and large genetic correlations between test day and lactation yield (0.86-0.99) indicate a potential use of test day records in genetic evaluations.

genetic evaluation; genetic parameters; milk yield; multivariate analysis; selection; test day model; variance components


MELHORAMENTO, GENÉTICA E REPRODUÇÃO

Parâmetros genéticos para as produções de leite no dia do controle e da primeira lactação de vacas da raça Holandesa1 1 Pesquisa financiada pela Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), processo 99/03838-6 e pelo Prodetab/Embrapa, integrante do Projeto 054-02/99.

Genetic parameters for test day and lactation milk yields of first lactation Holstein cows

Cláudio Manoel Rodrigues de MeloI; Irineu Umberto PackerII; Cláudio Nápolis CostaIII; Paulo Fernando MachadoII

IProfessor da UFSC. Departamento de Aqüicultura - AQI, CEP: 88040-900, Florianópolis, SC, Brasil (cmrmelo@cca.ufsc.br)

IIProfessores do Departamento de Zootecnia - USP/ESALQ, CEP: 13418-900, Piracicaba, SP, Brasil (iupacker@esalq.usp.br)

IIIPesquisador da Embrapa Gado de Leite, CEP: 36038-330, Juiz de Fora, MG, Brasil (cnc8@cnpgl.embrapa.br)

RESUMO

Foram utilizados 263.390 registros de produção de leite no dia do controle (PDC) de 32.448 primeiras lactações de vacas da raça Holandesa com parto no período de 1991 a 2001, para estimar componentes de variância e parâmetros genéticos utilizando um modelo animal e a metodologia REML. Os dados de produção foram coletados pelo Serviço de Controle Leiteiro da Associação Brasileira de Criadores de Bovinos da Raça Holandesa e suas afiliadas estaduais e disponibilizados pela Embrapa Gado de Leite. O modelo para ajuste da produção de leite até 305 dias (P305) incluiu os efeitos fixos de rebanho-ano de parto, época do parto e idade da vaca ao parto, com termos linear e quadrático, e os efeitos aleatórios de animal e erro. Os mesmos efeitos foram incluídos no modelo para as produções de leite no dia do controle (PDC) tanto sob modelo uni- e bi-caráter como sob modelo de repetibilidade (MRS), exceto para o efeito fixo de grupo contemporâneo, definido por rebanho-ano-mês do controle. Alternativamente, ajustou-se um segundo modelo de repetibilidade (MRF), que além dos efeitos presentes no MRS, incluiu as covariáveis que descrevem a curva da lactação: dias em lactação (DEL)/305 e ln(305/DEL) com termos linear e quadrático. As estimativas de h2 para as PDC, com MRS e MRF foram 0,43 e 0,30, respectivamente. As estimativas de herdabilidade (h2) para as PDC variaram de 0,22 (PDC1) a 0,36 (PDC4) com o modelo uni-caráter. Para o modelo bi-caráter, as estimativas variaram de 0,23 (PDC1) a 0,33 (PDC3 e PDC4). Em ambos os modelos, observaram-se valores inferiores de h2 no início e no fim do período de lactação. A estimativa h2 para a P305 com o modelo uni-caráter foi 0,27, enquanto, pelo modelo bi-caráter, estas estimativas variaram de 0,27 a 0,30. As correlações genéticas (rg) entre as PDC e a P305 foram altas, variando de 0,86 (PDC1 e P305) a 0,99 (PDC3 e P305). Estimativas de herdabilidade maiores para as PDC que para a P305 e a alta correlação com a mesma indicam potencial de uso das PDC nas avaliações genéticas de animais da raça Holandesa no Brasil. Embora predominantemente altas, as estimativas de rg entre as PDC não foram homogêneas (0,64-1,0); entretanto as maiores freqüências foram para valores próximos ou iguais a 1. Assim, modelos de regressão aleatória devem ser também avaliados para se concluir sobre a melhor utilização das PDC da raça Holandesa no Brasil.

Palavras-chave: análise bi-caráter, componentes de variância, modelo de repetibilidade, parâmetros genéticos, produções de leite do dia do controle

ABSTRACT

Covariance components for test day records and lactation milk yield using 263.390 records of 32.448 first lactation Holstein cows, were estimated using animal models by REML. Besides the lactation model, two alternative repeatability models (RM) were analyzed. Lactation model included fixed effects of herd-year-season and age of cow with linear and quadratic terms, and random effects of animal and error. The first model for test-day yield (RMF) included the same effects, but fixed effect of contemporary group, defined as herd-year-month of test. Alternatively another model for test-day yield (RMF) used a logarithmic polynomial sub-model for the shape of the lactation curve. Heritability for lactation yield (0.27) was smaller than those estimated by RMF and RMS, 0.30 and 0.43, respectively. Heritability estimates for univariate (0.22-0.36) and bivariate models (0.23-0.33) for test day milk yields were found to be smallest during early and late lactation. Heritability estimate for lactation milk yield when estimated by univariate model (0.27) was smaller than estimates obtained by bivariate models (0.27-0.30). Genetic correlations were higher between consecutive test days than between test days in the beginning and end of lactation. Larger heritability estimates for test day models and large genetic correlations between test day and lactation yield (0.86-0.99) indicate a potential use of test day records in genetic evaluations.

Key Words: genetic evaluation, genetic parameters, milk yield, multivariate analysis, selection, test day model, variance components

Introdução

A produção de leite até 305 dias (P305) é o indicador zootécnico geralmente utilizado para as comparações entre vacas nas decisões de manejo e descarte nos rebanhos leiteiros. A P305 é calculada usando os registros de produção provenientes dos controles leiteiros, que são realizados em intervalos aproximadamente mensais.

A utilização da P305 nas comparações entre os animais requer o seu ajuste para os efeitos de rebanho, ano de parto, época de parto, duração da lactação, entre outros. Entretanto, estes efeitos são considerados em termos médios, pois são definidos com base na data do parto.

Uma alternativa à utilização da P305 seria o uso da produção de leite no dia do controle (PDC). Este procedimento possibilita um ajuste mais preciso para os efeitos de ambiente temporário, visto que ele permite ajustar para os efeitos presentes em um dado dia ou mês de controle (Ptak & Schaeffer, 1993; Swalve, 1995), viabiliza o uso de mais dados de um mesmo animal e o uso de dados de animais com lactações em andamento ou parciais, sem necessidade de cálculo da P305, possibilitando-se avaliações genéticas mais freqüentes e redução do intervalo de gerações (Olori, 1997).

Em virtude do uso de mais dados de uma mesma vaca, as avaliações com base nas PDC, quando comparados àquelas utilizando a P305, podem ser mais confiáveis, o que é evidenciado pelas maiores estimativas de herdabilidade obtidas para as PDC (Visscher & Goddard, 1995; Strabel & Szwaczkowski, 1997). Neste contexto, Ptak & Schaeffer (1993) indicam a necessidade de considerar, no modelo de ajuste, a forma da curva de lactação, que pode ser incorporada usando coeficientes de regressão do caráter, em função do número de dias em lactação (DEL).

Outra consideração sobre os modelos para o ajuste das PDC refere-se à formação dos grupos contemporâneos. Danell (1982) recomenda o agrupamento de todas as vacas de determinado rebanho-data do controle. Em vários estudos (Reents et al., 1995; Swalve, 1995; Pösö et al., 1998; Rekaya et al., 1999), as estimativas de herdabilidade foram maiores quando os grupos contemporâneos foram definidos como rebanho-data do controle que quando definidos como rebanho-estação de parto. Isto também é observado quando se compara as estimativas de herdabilidade obtidas para as PDC àquelas obtidas para a P305 (Swalve, 1995; Strabel & Szwaczkowski, 1997; Ferreira et al., 2003), sugerindo maior confiabilidade nas estimativas dos valores genéticos dos animais. Neste contexto, o uso de modelos para o ajuste das PDC em substituição aos modelos para o ajuste da P305 tem sido investigado por vários autores (Van Tassell et al., 1992; Reents, et al., 1995; Rekaya et al., 1995; Wiggans & Goddard, 1997; Vargas et al., 1998).

Entre as alternativas para estimar componentes de variância e parâmetros genéticos adotando-se as PDC, o ajuste do modelo de repetibilidade (MR) caracteriza-se como a mais simples. O MR assume variâncias genéticas aditivas e de ambiente permanente constantes ao longo da lactação e correlações genéticas e de ambiente permanente entre as produções de uma mesma vaca iguais a um. Alguns estudos têm mostrado que correlações genéticas entre as PDC, dentro da primeira lactação, são próximas a um (Gadini, 1997; Gadini et al., 1998; Machado et al., 1998, Rekaya et al., 1999).

Pode-se, também, considerá-las como caracteres correlacionados e, então, analisá-las como medidas repetidas. Embora não assuma nenhuma estrutura de covariância entre os registros sucessivos, este procedimento é computacionalmente mais exigente, e, mesmo assim, é necessário conhecer a magnitude das covariâncias entre as PDC na lactação e entre as covariâncias e a P305, para avaliar a potencialidade e a estratégia de sua utilização nos procedimentos de seleção.

Alguns estudos desta natureza (Swalve, 1995; Machado, 1997; Rekaya et al., 1999) têm permitido concluir que estimativas de herdabilidades são menores para as produções do início e fim do período de lactação, em razão de menores estimativas de variância genética aditiva e de maiores estimativas de variâncias residuais para as PDC nestes períodos da lactação. De modo geral, estes resultados sugerem que as PDC podem ser usadas na seleção dos animais em substituição à P305. Entretanto, a escolha dos controles a serem utilizados deve considerar as estimativas de herdabilidade para as PDC e as estimativas de correlação genética entre as PDC e a P305. Assim, as PDC que apresentam maiores herdabilidades e maiores correlações genéticas com a P305 seriam as preferidas nos processos de seleção (El Faro, 2002).

Os mesmos estudos têm mostrado que as correlações genéticas entre as PDC e P305 são superiores no período intermediário da lactação quando comparadas àquelas do período inicial e do período final da lactação (Swalve, 1995; Rekaya et al., 1999), embora Rekaya et al. (1995) e Machado (1997) tenham relatado estimativas de correlações genéticas superiores entre PDC e P305 no início e fim da lactação.

Objetivou-se, neste estudo, estimar componentes de covariância e parâmetros genéticos, empregando modelos de repetibilidade, uni e bi-caráter, para o ajuste das produções de leite no dia do controle e para a produção de leite até 305 dias de primeiras lactações de vacas da raça Holandesa.

Material e Métodos

Dados

Os registros de produção de leite foram coletados sob a supervisão do Serviço de Controle leiteiro (SCL) da Associação Brasileira de Criadores de Bovinos da Raça Holandesa (ABCBRH) e de suas filiadas estaduais e disponibilizados pelo Arquivo Zootécnico Nacional (AZN), gerenciado pela Embrapa Gado de Leite. Foram utilizados os dados de 39.048 lactações de 531 rebanhos dos estados de São Paulo (SP), Minas Gerais (MG), Rio Grande do Sul (RS), Santa Catarina (SC) e Espírito Santo (ES) controladas no período de 1991 a 2001.

Foram definidas quatro épocas de parto: 1 - janeiro a março; 2 - abril a junho; 3 - junho a setembro e 4 - outubro a dezembro. O grupo de contemporâneos (GC) foi definido como rebanho-ano-mês do controle nos modelos para ajuste das PDC e rebanho-ano de parto no modelo para ajuste da P305 e, em ambos, aplicou-se a restrição de que cada GC deveria conter pelo menos quatro vacas de 18 a 48 meses de idade.

Assim, no ajuste do modelo para P305, foram usados dados de produção de leite de 28.437 vacas pertencentes a 467 rebanhos. Para ajuste dos modelos para PDC,foram utilizadas somente as produções, coletadas entre cinco e 305 dias após o parto, de vacas com pelo menos três controles de produção, resultando em 32.448 primeiras lactações, com 263.390 registros de produção de leite provenientes de controles mensais, para análise.

O arquivo contendo dados de genealogia dos animais, para a estruturação da matriz de coeficientes de parentesco (A), foi caracterizado pela inclusão de 3.737 registros de touros aos 41.069 registros de vacas. Entretanto, animais não informativos (i.e., ambos pais desconhecidos e somente um descendente) foram eliminados. Assim, a matiz A incluiu 36.732 animais para o ajuste do modelo da P305 e 40.770 animais para o ajuste dos demais modelos.

Modelos

O modelo utilizado para o ajuste das PDC sob modelo uni- e bi-caráter ou para o ajuste da P305 foi:

em que: yijk é a PDC ou a P305; GCi, o efeito do grupo contemporâneo , formado pela combinação de rebanho-ano-mês do controle (PDC) ou de rebanho-ano de parto (P305); EPOCAj, a época do parto j; b1 e b2, os coeficientes de regressão linear e quadrático do caráter em função da idade da vaca ao parto; xijk, a idade da vaca ao parto, em dias; aijk, o efeito genético aditivo aleatório do animal k pertencente ao rebanho-ano de parto i e a época de parto j; e eijk, o erro aleatório associado a cada observação.

Para o ajuste das PDC sob modelo de repetibilidade (MRS), foi utilizado o seguinte modelo:

em que: yijk é a PDC da vaca k, que pariu na época j e foi controlada no grupo i; RAMi, o efeito do i-ésimo grupo de controle formado pela combinação de rebanho, ano, mês do controle; EPOCAj, a época do parto j; b1 e b2, os coeficientes de regressão linear e quadrático da PDC, em função da idade da vaca ao parto; xk, a idade da vaca k ao parto, em dias; ak, o efeito genético aditivo aleatório do animal k; pek, o efeito permanente da vaca k; e eijk, o erro aleatório associado a cada observação.

Alternativamente, ajustou-se outro MR incluindo os coeficientes de regressão fixa da produção de leite, em função de covariáveis que descrevem a curva de lactação média da população (MRF):

em que: q1 a q4 são os coeficientes de regressão da PDC em função das covariáveis znk que descrevem a curva da lactação, em que: em que: z1k = c; z2k = c2; z3k = 1n(1/c); z4k = (1n(1/c))2 e c = DEL/305. Os demais efeitos são como descritos na eq. (2).

Na forma matricial, o modelo (1) pode ser escrito como:

em que y é o vetor das observações; b é o vetor dos efeitos fixos; a é o vetor de efeito aleatório de animal; e é o vetor de efeito aleatório residual; e X e Z são as matrizes de incidência relativa às observações, para efeitos fixos e efeito aleatório de animal, respectivamente. Assumiu-se , em que: G é a matriz de covariância genética aditiva; R, a matriz de covariâncias residuais; A, a matriz de coeficientes de parentesco entre os animais; IN, uma matriz identidade de ordem igual ao número de animais e Ä, o operador de produto direto entre matrizes.

Da mesma forma, os MRF e MRS, eq. (2), podem ser descritos como:

em que y, b, a, e, X, e Z são como definido anteriormente; p é o vetor de efeitos permanentes de ambiente; e ZD é a matriz de incidência relativa ao efeito aleatório permanente de ambiente. Assumiu-se , , e, portanto, . Os componentes de covariância foram estimados utilizando-se o programa REMLF90 (Misztal, 1999), baseado na metodologia de máxima verossimilhança restrita (REML), com o algoritmo de Maximização da Esperança (EM), e um processo de aceleração da convergência. Definiu-se como critério de convergência o valor do quadrado das diferenças entre estimativas consecutivas menor que 10–9.

Procedeu-se à decomposição canônica das matrizes de covariância genética (G) e residual (R) estimadas para as PDC com intuito de definir a provável ordem necessária de ajuste das PDC sob modelos de regressão aleatória.

A eficiência relativa de seleção (ER) baseada na produção dos 10 controles, na produção dos cinco primeiros controles, na produção de controles alternados de dois em dois meses (PDC1, PDC4, PDC7 e PDC10) ou, ainda, na produção de cada controle individualmente, em relação à seleção pela P305, foi calculada para determinar o ganho/perda resultante da seleção pelas produções do dia do controle em relação a P305. Definiu-se o agregado como a = 30xPDC1+30xPDC2+...+30xPDC10, de forma a ter uma equivalência com a P305 (correlação esperada entre a e P305 igual a 1). A ER foi calculada como:

a correlação do agregado com a índice, G é a matriz de covariância genética aditiva ou parte dela (dependendo de quais controles estão sendo usados no cálculo da ER), P é a matriz de covariância fenotípica ou parte dela e 0,27, a herdabilidade estimada para P305.

Resultados e Discussão

As informações sobre número de lactações e de controles, médias de produção em até 305 dias e nos controles, e número de animais considerados no estudo são apresentadas na Tabela 1.

As PDC determinam a forma típica da curva de lactação, com ligeiro aumento da produção inicial de cerca de 22 kg até o pico da lactação (período de 55 a 87 dias de lactação) com aproximadamente 24 kg, decrescendo posteriormente até 19 kg de leite ao final da lactação. A variação entre as produções de leite foi homogênea (desvio-padrão médio de 6,65 kg de leite) nos controles da lactação.

Na Tabela 2, são apresentadas as estimativas dos componentes de variância e parâmetros genéticos para a P305 e para as PDC sob modelo de repetibilidade.

As estimativas de herdabilidade para as PDC (0,43 e 0,30) foram maiores que aquela obtida para a P305 (0,27), semelhantes às reportadas por Costa et al. (2000), Freitas et al. (2001) e Ferreira et al. (2003) e dentro dos limites (0,22-0,49) observados na literatura (Costa et al., 1999; Machado et al., 1999; Torres et al., 1999).

Valores de herdabilidades superiores para as PDC em relação aos obtidos para a P305 foram também relatados por Strabel & Szwaczkowski (1997) e Ferreira et al. (2003). Portanto, os resultados do presente estudo confirmam que os MR para o ajuste das PDC resultam em maiores estimativas de herdabilidade que os modelos para o ajuste da P305.

As estimativas dos componentes de variância diferiram entre os MR, com a maior diferença entre as variâncias genéticas (6,763), seguidas daquela entre as variâncias residuais (1,872), e uma diferença muito pequena entre os componentes do efeito permanente de ambiente (0,205). A estimativa de herdabilidade (0,43) obtida com o MRS foi maior que a (0,30) encontrada por meio do MRF. A redução na estimativa da variância residual, pela inclusão das covariáveis que descrevem a forma da curva da lactação era esperada, no entanto, não se esperava a redução na estimativa da variância genética aditiva, pois os coeficientes que descrevem a forma da curva de lactação foram tomados como fixos.

As estimativas das variâncias obtidas com o ajuste dos MR foram semelhantes às reportadas por Ferreira et al. (2003) e encontram-se dentro dos limites apresentados na literatura (Reents et al. 1995; Swalve 1995; Olori, 1997; Pösö et al. 1998; Rekaya et al. 1999) para a variância genética (12,70-2,77), para o efeito permanente de ambiente (8,60-3,26) e residual (7,40-3,80).

Os valores de herdabilidade (0,30) e de repetibilidade (0,64) obtidos com o ajuste do MRF, foram semelhantes aos reportados por Rekaya et al. (1999). Entretanto, o valor de herdabilidade de 0,43, verificado no ajuste do MRS foi superior a 0,31, obtido por Pool et al. (2000). Estes valores (0,43 e 0,30) também se assemelham às estimativas reportadas por Olori (1997) para modelos de repetibilidade que consideravam a forma da curva da lactação por meio de regressões polinomial quadrática, cúbica e de quarta ordem (MRF), ou na ausência destas regressões (MRS).

As estimativas dos componentes de variância, herdabilidades para as PDC, sob modelo de análise uni-caráter e para a P305 e ER, são apresentadas na Tabela 3.

As estimativas de variância genética aditiva para as PDC aumentaram do primeiro controle (5,01 kg2) até o quarto controle (7,88 kg2), reduzindo-se daí por diante até o último controle (4,80 kg2). As maiores estimativas foram obtidas para as PDC3 e PDC4. Ao contrário do observado para as variâncias genéticas, as estimativas das variâncias residuais foram superiores para as produções do início (PDC1) e do fim (PDC10) do período de lactação.

Os valores de herdabilidade para PDC aumentaram do início da lactação até a fase intermediária (quarto controle); a partir deste ponto, até o final da lactação, ocorreu diminuição destes valores (Figura 1). O maior valor das estimativas de herdabilidade para as PDC na fase intermediária da lactação ocorreu mais em função do aumento nas estimativas das variâncias genéticas que em função da redução das variâncias residuais. Estes resultados são similares aos obtidos por Olori (1997), Rekaya et al. (1999) e Ferreira et al. (2003), mas não aos apresentados por El Faro (2002), que obteve maior estimativa de herdabilidade para os controles do início e do fim do período da lactação. As estimativas de herdabilidade para as PDC foram geralmente superiores à obtida para a P305, exceto para as PDC1, PDC9 e PDC10 (Figura 1). Estes resultados também se assemelham aos reportados em vários outros estudos (Meyer et al., 1989; Swalve, 1995; Kettunen et al., 1998; Machado et al., 1998).


As estimativas de herdabilidade resultantes das análises bi-caráter (Tabela 4) foram semelhantes àquelas obtidas pelas análises uni-caráter e situaram-se dentro do intervalo de valores relatados por Rekaya et al. (1999) e Ferreira et al. (2003).

As estimativas das correlações genéticas e fenotípicas entre as PDC e entre estas e a P305 encontram-se na Tabela 5.

As correlações genéticas para as PDC variaram de 0,64 (entre PDC1 e PDC10) a 1 (entre PDC3 e PDC7 e entre PDC4 e PDC9), entretanto, as maiores freqüências foram para valores próximos ou iguais a 1 (71%>0,95 e 47%>0,98). As correlações entre as PDC a partir do sexto controle (>PDC6) foram todas iguais a 0,98. Entretanto, as correlações genéticas entre as PDC do início (<PDC3) e do final (>PDC8) da lactação foram inferiores a 0,91 (Tabela 5). As correlações fenotípicas variaram de 0,27 (entre PDC1 e PDC10) a 0,78 (entre PDC3 e PDC4; entre PDC4 e PDC5; entre PDC5 e PDC6; entre PDC6 e PDC7; e entre PDC7 e PDC8) e apresentaram o mesmo comportamento das correlações genéticas, com valores superiores entre as PDC adjacentes (Tabela 5). As correlações residuais (resultados não apresentados) foram inferiores às demais (genéticas e fenotípicas), com valores variando de 0,69 (entre PDC5 e PDC7; entre PDC7 e PDC8; e entre PDC8 e PDC10) a 0,15 (entre PDC1 e PDC10). Estes resultados, com valores de correlações genéticas superiores entre as produções dos controles consecutivos, assemelham-se aos obtidos por Meyer et al. (1989), Swalve (1995), Gadini (1997) e Olori et al. (1999).

As correlações genéticas entre as PDC e a P305 foram altas, variando de 0,86 (entre PDC1 e P305) a 0,99 (entre PDC3 e P305), com maiores valores entre as PDC do período intermediário da lactação e a P305 (Tabela 5). Dessa forma, se as decisões de descarte dos animais forem baseadas nas PDC da fase intermediária da lactação, que apresentam maiores valores de herdabilidade (Tabela 3) e maiores correlações genéticas com a P305 (Tabela 5), poderia se obter maiores ganhos genéticos na P305 que se decisões fossem tomadas com base nas PDC da fase inicial ou final da lactação.

Assim, maior ganho genético na P305 pode ser obtido por seleção indireta baseada nas PDC3 a PDC8, particularmente nas PDC3 e PDC4, que apresentam maior eficiência relativa que a seleção direta para a P305 (Tabela 3). A seleção realizada com base nas PDC poderia resultar em ganhos 29% superiores aos obtidos da seleção sob P305 (utilizando os 10 controles), além de possibilitar avaliações com produções parciais (cinco primeiros controles) e controles esparsos (controles bimensais) ainda com ganhos superiores aos esperados para P305 (23 e 27%, respectivamente, para os cinco primeiros controles e para controles bimensais) (Tabela 3). Ribas & Perez (1990) relataram maior ER para o PDC5 e PDC6 em estudo como PDC de primeiras lactações de vacas Holandesas.

Estimativas de covariâncias genéticas e fenotípicas entre as PDC são apresentadas na Tabela 6.

As covariâncias genéticas e fenotípicas entre as PDC foram maiores particularmente na fase intermediária da lactação, diminuindo em magnitude à medida em que os mesmos se distanciavam. Considerando-se que as estimativas de correlações genéticas entre as PDC foram altas e próximas à unidade entre a maioria das PDC (71% > 0,95), pode-se assumir que as correlações não diferem de 1 entre as PDC. Este fato sugere que um modelo de repetibilidade (Rekaya et al., 1999; Melo et al., 2002) poderia ser utilizado para o ajuste das PDC, visto que o mesmo é mais parcimonioso, pois necessita estimar apenas três componentes de covariância, e, portanto, é menos exigente em termos de demanda computacional e tempo de processamento.

Apesar de maior simplicidade e eventual facilidade de análise ao assumir variâncias constantes e correlações iguais a um entre as PDC, o modelo de repetibilidade não considera a forma da curva da lactação. Contudo, ajustar todas as produções disponíveis ao longo da lactação em um modelo multi-caráter é mais difícil em razão do maior número de parâmetros a serem estimados. Neste contexto, os modelos de regressão aleatória, com submodelos que permitem o ajuste da forma da curva da lactação e consideram a estrutura de covariância entre as PDC, caracterizam-se como uma interessante e promissora alternativa para o ajuste das PDC (Brotherstone et al., 2000; Pool et al., 2000; Kettunen et al., 2000; Pool & Meuwissen, 1999; Jakobsen et al., 2002) devendo ser investigados em estudos futuros.

Os seguintes autovalores (li) 60,06; 3,01; 0,52; 0,28; 0,21; 0,13; 0,10; -0,02, -0,22; -0,94 e 76,44; 24,50; 13,92; 8,31; 5,75; 4,81; 4,56; 4,14; 3,74; 3,67 são produzidos, respectivamente, pela decomposição canônica da matriz de covariância genética e fenotípica. Pela decomposição canônica da matriz G, tem-se que os dois primeiros autovalores contabilizam 96% da variância genética (desconsiderando o sinal dos demais). Todavia, para a matriz R, são necessários os quatro primeiros autovalores para contabilizar o mesmo percentual da variância residual. Os autovetores associados a estes autovalores (dois primeiros de G e quatro primeiros de R) são apresentados na Figuras 2 e 3, respectivamente.



Nota-se que os dois primeiros autovetores da matriz de covariância genética (Figura 2) correspondem, aproximadamente, aos efeitos constante e linear de dias em lactação e os quatros primeiros autovetores da matriz de covariância residual (Figura 3) aos efeitos constante, linear, quadrático e cúbico. Pode-se, portanto, sugerir que um modelo de regressão aleatória utilizando um polinômio de primeira ordem (efeitos constante e linear) poderia ser suficiente para ajustar o efeito genético no modelo. Todavia, ordens superiores (terceira ordem) seriam necessárias para ajustar os efeitos de ambiente.

Conclusões

A não-inclusão das covariáveis que descrevem a curva de lactação no modelo de ajuste das PDC, sob modelo de repetibilidade, resultou em estimativas inflacionadas da variância genética aditiva.

As PDC da fase intermediária da lactação podem ser utilizadas como critérios de seleção dos animais como forma de obter predição mais confiável do valor genético dos animais, particularmente aquelas que têm maiores herdabilidade (produções dos controles dois a oito) que a da produção de leite até 305 dias e são altamente correlacionadas com esta. Nesta situação, a seleção indireta baseada nas PDC pode ser mais eficiente que a seleção direta para P305. Assim, a utilização das PDC poderia possibilitar avaliações por produções parciais (e.g. PDC1 a PDC5) ou controles esparsos (e.g. bimensais) e ainda resultar em ganhos superires a 20% em relação a P305.

Estimativa de correlações genéticas altas e próximas à unidade entre a maioria das PDC sugerem que, eventualmente, pode-se assumir correlação igual a um entre as PDC e ajustar um modelo de repetibilidade, que é mais parcimonioso e menos exigente em termos de demanda computacional.

A decomposição canônica das matrizes de covariâncias genética e residual mostrou que um modelo de regressão aleatória utilizando um polinômio de primeira ordem pode ser suficiente para ajustar o efeito genético. Contudo, um polinômio de ordem três é necessário para ajustar os efeitos de ambiente. Neste contexto, é importante verificar o ajuste destes modelos em trabalhos futuros para se concluir sobre a melhor utilização das PDC na avaliação genética da raça Holandesa no Brasil.

Agradecimento

Ao professor Ignacy Misztal, da Universiy of Georgia, EUA, por disponibilizar os programas computacionais RENUMMAT3 e REMLF90; à ABCBRH, ao Prodtab; e aos revisores, pelas valiosas sugestões apresentadas.

Literatura Citada

Recebido em: 06/02/03

Aceito em: 29/10/04

  • BROTHERSTONE, S.; WHITE, I.; MEYER, K. Genetic modeling of dairy milk yield using orthogonal polynomial and parametric curves. Animal Science, n.70, p.407-415, 2000.
  • COSTA, C.N. An investigation into heterogeneity of variance for milk and fat yields of Holstein cows in Brazilian herd environments. Genetic and Molecular Biology, v.22, n.3, p.375-381, 1999.
  • COSTA, C.N.; BLAKE, R.W.; POLLAK, E.J. et al. Genetic analysis of Holstein cattle populations in Brazil and United States. Journal of Dairy Science, v.83, n.12, p.2963-2974, 2000.
  • DANELL, B. Studies on lactation yield and individual test-day yields of Swedish dairy cow. II. Estimates of genetic and phenotypic parameters. Acta Agriculturae Scandinavica, v.32, p.103-114, 1982.
  • EL FARO, L. Estimação de componentes de (co)variância para produção de leite no dia do controle de primeiras lactações de vacas Caracu, aplicando-se ";test day models"; de dimensão finita e modelos de regressão aleatória Jaboticabal: Universidade Estadual Paulista, 2002. 102p. Tese (Doutorado) - Universidade Estadual Paulista.
  • FERREIRA, W.J. Avaliação genética de bovinos da raça Holandesa usando a produção de leite no dia do controle. Revista da Sociedade Brasileira de Zootecnia, v.32, n.2, p.295-303, 2003.
  • FREITAS, A. F.; DURÃES, M.C.; VALENTE, J. et al. Parâmetros genéticos para produção de leite e gordura nas três primeiras lactações de vacas Holandesas. Revista Brasileira de Zootecnia, v.30, n.3, p.709-713, 2001.
  • GADINI, C.H. Genetic evaluation of test day production traits and somatic cell scores Licoln: University of Nebrask, 1997. Thesis (Ph.D) - University of Nebrask, 1997.
  • GADINI, H.C.; KEOWN, J.F.; Van VLECK, L.D. Genetic parameters of test day milk, fat and protein yields. Genetic parameters of test-day milk yield of Holstein cows (compact disc). In: WORLD CONGRESS OF GENETICS APPLIED LIVESTOCK PRODUCTION, 6., 2002, Armidale. Proceedings... Armidale: University of New England, 1998.
  • JAKOBSEN, J.H.; MADSEN, P.; JENSEN, J. et al. Genetic parameters for milk production and persistency for Danish Holsteins estimated in random regression models using REML. Journal of Dairy Science, v.85, n.6, p.1607-1616, 2002.
  • KETTUNEN, A.; MÄNTYSAARI, E.A.; PÖSÖ, J. Estimation of genetic parameters for daily milk yield of primiparous Ayrshire cows by random regression test-day models. Livestock Production Science, v.66, p.251-261. 2000.
  • KETTUNEN, A.; MÄNTYSAARI, E.A.; STRANDÉN, I. et al. Estimation of genetic parameters for first lactation test day milk production using random regression models (compact disc). In: WORLD CONGRESS OF GENETICS APPLIED LIVESTOCK PRODUCTION, 6., 1998, Armidale. Proceedings... Armidale: University of New England, 1998.
  • MACHADO, S.G. Parâmetros genéticos e de ambiente da produção de leite no dia do controle da primeira lactação de vacas da raça Holandesa Jaboticabal: Universidade Estadual Paulista, 1997. 76p. Dissertação (Mestrado) - Universidade Estadual Paulista, 1997.
  • MACHADO, S.G.; FREITAS, M.A.R.; GADINI, C.H. Genetic parameters of test-day milk yield of Holstein cows (compact disc). In: WORLD CONGRESS OF GENETICS APPLIED LIVESTOCK PRODUCTION, 6., Armidale, 1998. Proceedings... Armidale: University of New England, 1998.
  • MACHADO, S.G.; FREITAS, M.A.R.; GADINI, C.H. Genetic parameters of test day milk yields of Holstein cows. Genetic and Molecular Biology, v.22, n.3, p.383-386. 1999.
  • MELO, C.M.R.; COSTA, C.N.; PACKER, I.U. et al. Parâmetros genéticos para as produções de leite de controles individuais e total de primeiras lactações de vacas da raça Holandesa. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 39., 2002, Recife. Anais... Recife: Sociedade Brasileira de Zootecnia, 2002. CD-ROM
  • MEYER, K.; GRASER, H.U.; HAMMOND, K. Estimates of genetic parameters for first lactation test day production of Australian Black and White cows. Livestock Production Science, v.21, p.177-199, 1989.
  • MISZTAL, I. 1999. REMLF90: manual. ftp://nce.ads.uga.edu/pub/ignacy/blupf90/ docs/remlf90.pdf (05 de jan. 2001).
  • OLORI, V.E. Utilization of daily milk records in genetic evaluation of dairy cattle Edinburgh, 1997. 1v. Thesis (Ph.D) - University of Edinburgh, 1997.
  • OLORI, V.E.; HILL, W.G.; McGUIRK, B.J. et al. Estimating variance e components for test day milk records by restricted maximum likelihood with a random regression animal model. Livestock Production Science, v.61, n., p.53-63, 1999.
  • POOL, M.H.; JANSS, L.L.G.; MEUWISSEN, T.H.E. Genetic parameters of Legendre polynomials for first parity lactation curves. Journal of Dairy Science, v.83, n.11, p.2640-2649, 2000.
  • POOL, M.H.; MEUWISSEN, T.H.E. Prediction of daily milk yields from a limited number of test days using test day models. Journal of Dairy Science, v.82, n.7, p.1555-1564, 1999.
  • PÖSÖ, J.; MÄNTYSAARI, E.A.; LIDAUER, M. et al. mpirical bias in the pedigree indices of heifers evaluated using test day models (compact disc). In: WORLD CONGRESS OF GENETICS APPLIED LIVESTOCK PRODUCTION, 6., 1998, Armidale. Proceedings... Armidale: University of New England, 1998.
  • PTAK, E.; SCHAEFFER, L.R. Use of test day yields for genetic evaluation of dairy sires and cows. Livestock Production Science, v.34, n.1-2, p.23-34, 1993.
  • REENTS, R.; DEKKERS, J.C.M.; SCHAEFFER, L.R. Genetic evaluation for somatic cell score with a test day model for multiple lactations. Journal of Dairy Science, v.78, n.12, p.2858-2870, 1995.
  • REENTS, R.; JAMROZICK, J.; SCHAEFFER, L.R. et al. Estimation of genetic parameters for test day records of somatic cell score. Journal of Dairy Science, v.78, n.12, p.2847-2857, 1995.
  • REKAYA, R.; BEJAR, R.; CARABAÑO, M.J. et al. Genetic parameters for test day measurements in Spanish Holstein-Friesian. In: INTERBULL, 11., 1995, Prague. Proceedings... Prague: Czech Republic, 1995. p.1-8 (Bulletin 11).
  • REKAYA, R.; CARABAÑO, M.J.; TORO, M.A. Use de test day yields for the genetic evaluation of production traits in Holstein-Friesian cattle. Livestock Production Science, v.57, p.203-217, 1999.
  • RIBAS, M.; PEREZ, B. Monthly test day milk records and yield at 244 days. II. Genetic parameters in first lactation. Cuban Journal of Agricultural Science, v.24, n.2, p.129-144, 1990.
  • STRABEL, T.; SZWACZKOUSKI, T. Additive genetic and permanent environmental variance components for test day milk yields in Black-white cattle. Livestock Production Science, v.48, p.91-98, 1997.
  • SWALVE, H.H. The effect of test day models on the estimation of genetic parameters and breeding values for dairy yield traits. Journal of Dairy Science, v.78, n.4, p.929-938, 1995.
  • TORRES, R.A.; BERGMAN, J.A.G.; COSTA, C.N. et al. Ajustamento para heterogeneidade de variância para a produção de leite entre rebanhos da raça Holandesa no Brasil. Revista Brasileira de Zootecnia, v.28, n.2, p.295-305. 1999.
  • VAN TASSELL, C.P.; QUAAS, R.L.; EVERETT, R.W. Parameter estimates for 305-day residual records. Journal of Dairy Science, v.75, suppl. 1, p.251, 1992.
  • VARGAS, B.; PEREZ, E.; Van ARENDONK, J.A.M. Analysis of test day yield data of Costa Rican dairy cattle. Journal of Dairy Science, v.81, n.1, p.225-271, 1998.
  • VISSCHER, P.M.; GODDARD, M.E. Genetic parameters for milk yield, survival, workability and type traits for Australian dairy cattle. Journal of Dairy Science, v.78, n.1, p.205-220, 1995.
  • WIGGANS, G.R.; GODDARD, M.E. A computationally feasible test day model for genetic evaluation of traits in the United States. Journal of Dairy Science, v.80, n.8, p.1795-1800, 1997.
  • 1
    Pesquisa financiada pela Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), processo 99/03838-6 e pelo Prodetab/Embrapa, integrante do Projeto 054-02/99.
  • Datas de Publicação

    • Publicação nesta coleção
      13 Set 2005
    • Data do Fascículo
      Jun 2005

    Histórico

    • Aceito
      29 Out 2004
    • Recebido
      06 Fev 2003
    Sociedade Brasileira de Zootecnia Universidade Federal de Viçosa / Departamento de Zootecnia, 36570-900 Viçosa MG Brazil, Tel.: +55 31 3612-4602, +55 31 3612-4612 - Viçosa - MG - Brazil
    E-mail: rbz@sbz.org.br