Os controles individuais de produção de leite de vacas Gir, utilizados neste trabalho, foram coletados mensalmente pela Associação Brasileira de Criadores, com sede em São Paulo, SP, e Associação Brasileira de Criadores de Zebu, com sede em Uberaba, MG.

O arquivo base, do qual foram extraídas as informações utilizadas neste trabalho, continha 55.717 registros de controle leiteiro mensal, obtidos em duas ordenhas diárias, de partos ocorridos de 1973 a 1999, de diversas ordens de parto em 17 rebanhos e estavam armazenados no banco de dados do Arquivo Zootécnico Nacional, gerenciado pela Embrapa Gado de Leite. Este arquivo possuía igualmente a produção de leite em 305 dias de lactação (PL305) calculada de acordo com o método oficial reconhecido pelo Ministério da Agricultura (Brasil, 1986).

As características estudadas foram a produção de leite no dia do controle leiteiro (PLDC1 a PLDC10) e a produção de leite em 305 dias de lactação (PL305), das três primeiras lactações. Do arquivo base, foram gerados sub-arquivos de trabalho para verificações e consistências dos dados, utilizando para tal fim procedimentos disponíveis no software SAS (1990). Outras consistências e cálculos, bem como tabelas e gráficos, foram efetuadas utilizando os programas Word e Excel, da Microsoft (1997).

Inicialmente, foram eliminadas informações de lactações com menos de três controles mensais, rebanhos e anos com menos de 100 lactações. Posteriormente foram incluídas as variáveis idade da vaca ao parto, intervalo parto-primeiro controle leiteiro (ipc) e estações do parto e do controle leiteiro (estação 1 - meses de janeiro, fevereiro e março; Estação 2 - meses de abril, maio e junho; Estação 3 - meses de julho, agosto e setembro; Estação 4 - meses de outubro, novembro e dezembro).

As estações (do parto ou do controle) foram agrupadas dentro de rebanho e do ano (do parto ou do controle), criando-se as sub-classes de grupos contemporâneos de rebanho-ano-estação do parto (raep305) e rebanho-ano-estação do controle (raec1 a raec10).

Foram então impostas novas restrições aos dados, eliminando registros de vacas com idade ao parto inferior a 24 meses e superior a 120 meses, produções do primeiro controle com intervalo parto-primeiro controle leiteiro inferior a 5 dias e superior a 35 dias, grupos contemporâneos com menos de 10 observações e touros com menos de 2 filhas.

Foi gerado, finalmente, o arquivo para as análises, contendo 32.779 controles mensais de 3.605 lactações de 2.082 vacas, filhas de 281 touros, com partos ocorridos de 1987 a 1999 em 11 rebanhos, contendo as seguintes variáveis: número da vaca, número do pai, número da mãe, ordem do parto, raec1 a raec10, raep305, idade da vaca ao parto, ipc, duração da lactação, PLDC1 a PLDC10 e PL305.

Foi criado, para formação da matriz de numeradores dos coeficientes de parentesco (NRM), um arquivo de pedigrees, utilizado em todas as análises, contendo a identificação da vaca, do pai e da mãe, resultando em 2.082 pedigrees decifrados, com 3.407 animais diferentes e 9 animais que apresentaram endogamia média de 0,26.

Os componentes de variância, os parâmetros genéticos e os valores genéticos preditos dos animais foram estimados pelo método de máxima verossimilhança restrita (REML), sob modelo animal, utilizando o aplicativo MTDFREML (Multiple Trait Derivative-Free Restricted Maximum Likelihood), descrito por Boldman et al. (1995).

No método de máxima verossimilhança restrita, desenvolvido formalmente por Patterson & Thompson (1971), para modelos mistos balanceados ou não, o vetor de observações é modelado considerando duas partes diferentes, uma referente aos efeitos fixos e outra aos aleatórios, de maneira que a função de verossimilhança restrita das observações seja dada pela soma das funções de cada parte. A maximização da função de verossimilhança restrita da parte referente aos efeitos aleatórios, em relação aos componentes de variância, elimina o viés (tendência) resultante da perda de graus de liberdade na estimação dos efeitos fixos do modelo (Lopes et al., 1993).

O modelo animal, utilizado nas avaliações genéticas de bovinos leiteiros, é um substituto da metodologia de comparação com contemporâneos (MCC). A principal diferença entre as metodologias é a ênfase no parentesco entre os animais, em lugar do enfoque no touro, em que o desempenho das filhas de um touro é utilizado na determinação de seu mérito genético. No modelo animal, as avaliações genéticas são baseadas no animal e no seu parentesco com outros animais que estão sendo avaliados e, machos e fêmeas são considerados simultaneamente, isto é, cada animal irá influenciar as avaliações de seus parentes, pois há uma equação para cada indivíduo.

Nestas análises, o critério de convergência considerado foi a variância do simplex menor que 10^-9. A cada convergência, o programa era reiniciado, usando como valores de (co)variância iniciais aqueles obtidos na análise anterior, o que permitia checar possíveis 'máximos locais'. Este procedimento era repetido várias vezes até que os valores da função de verossimilhança entre duas tentativas repetidas não alterassem mais do que 10^-5 unidades.

O modelo proposto para o estudo, em termos matriciais, foi:

As pressuposições usuais para este modelo são:

As estimativas de herdabilidade foram calculadas como:

em que:

As estimativas dos valores genéticos preditos dos touros nas PLDC e PL305 foram organizados em arquivos para obtenção da porcentagem de touros em comum em diferentes níveis de seleção e para calcular as correlações de Spearman no conjunto global de dados.

As médias observadas, os desvios-padrão e os coeficientes de variação para produções de leite no dia do controle leiteiro (PLDC) e para produção de leite em 305 dias de lactação (PL305), podem ser visualizadas na Tabela 1.

Os coeficientes de variação (CV) foram elevados e aumentaram com o avançar da lactação, com valor mais elevado na PL305. Seriam, conforme Sampaio (1998), características muito instáveis. Os CV foram maiores que os citados nos trabalhos de Albuquerque (1984), Gadini (1985, 1997), Machado (1997) e Ferreira (1999), os quais apresentaram valores de 20 a 30%.

Os altos CV encontrados refletem, principalmente, desuniformidade de produções entre os rebanhos estudados. A diferença entre as produções dos rebanhos chegou a 66, 64, 61, 59, 58, 56, 55, 53, 51 e 52% nas PLDC1 a PLDC10, respectivamente, e de 63% na PL305.

As estimativas dos componentes de variância e de herdabilidade das PLDC e PL305 estão apresentadas na Tabela 2.

Os valores de herdabilidade estimadas para PLDC2 até PLDC4 e para PL305 tiveram valores próximos. Do quinto controle até o décimo, estes valores foram muito reduzidos em relação àqueles já citados. No primeiro controle de produção de leite, a estimativa de herdabilidade foi mais elevada.

Os menores valores para a herdabilidade, a partir do primeiro controle, se deveram à maior redução da variância genética aditiva (diminuição de até 70%, de 2,10 kg² para 0,62 kg²) em relação à variância fenotípica (queda de até 35%, de 9,15 kg² para 6,15 kg²). O maior valor da estimativa de herdabilidade, ocorrido no primeiro controle, pode ser explicado pela provável padronização de tratamento nutricional pré e pós- parto e pelo ajustamento da produção de leite para o efeito do intervalo parto-primeiro controle.

Foram observadas na literatura consideráveis desigualdades para as estimativas de herdabilidade, devido às disparidades entre as populações estudadas e aos diferentes métodos de análise utilizados.

Strabel & Szwaczkowski (1997) encontraram herdabilidade para PLDC muito superior à obtida para PL305, enquanto Swalve (1995) e Firat et al. (1997a) observaram maiores valores de herdabilidade para PL305.

A tendência dos valores de herdabilidade citados pela literatura, envolvendo animais das raças européias, mostram maiores valores na segunda fase da lactação (Wilmink, 1987; Pander et al., 1992; Firat et al., 1997b; Machado, 1997; Tijani et al., 1999), onde é alegado haver menor variação nas produções destes controles, devido à influência de meio ambiente ser mais expressiva no início e no final da lactação, contrariando os resultados deste trabalho. Kettunen et al. (1998) verificaram maior valor da herdabilidade no primeiro controle, confirmando os achados do presente estudo.

Sikka & Taneja (1981), em vacas Sahiwal, e Albuquerque (1984), com vacas da raça Gir, verificaram aumento nos coeficientes de herdabilidade da produção parcial e da produção parcial acumulada, respectivamente, até o terceiro controle e posterior decréscimo. Já Gadini (1985) notou tendência de diminuição da herdabilidade do primeiro ao quinto controle e posterior aumento até o décimo, em rebanho da raça Gir.

Na Tabela 3 estão expostos os resultados médios dos valores genéticos preditos de 281 touros e suas confiabilidades de estimação.

Para intensidade de seleção constante, a resposta à seleção seria, então, dependente da herdabilidade e do desvio-padrão fenotípico. Assim, quanto maior a diferença dos animais avaliados, maior poderá ser o diferencial de seleção, com maior progresso genético nas características de herdabilidade mais alta.

Se os reprodutores são avaliados pelas produções de suas filhas em diferentes controles, aqueles obtidos em classes de controles de maior desvio-padrão fenotípico, provavelmente, terão maiores valores genéticos. As médias e amplitudes de valores genéticos dos reprodutores tenderam a ser maiores nas classes do primeiro ao quarto controles. As PLDC1 a PLDC4 foram, também, aquelas com desvios-padrão (Tabela 1) e herdabilidades (Tabela 2) mais elevadas. Então, o progresso genético seria maior nas PLDC1 a PLDC4.

Na Tabela 4, são apresentados os percentuais de touros em comum, para níveis crescentes de animais selecionados pela classificação da PL305 em comparação com as das PLDC, de acordo com os valores genéticos preditos estimados.

As altas correlações de ordem entre as avaliações dos touros revelam grande associação linear entre os valores genéticos previstos usando as PLDC e PL305, indicando que a ordem dos animais nas diferentes avaliações tende a ser a mesma, no conjunto global de dados, conforme afirmaram Ptak & Schaeffer (1993) e Verneque et al. (1998).

Entretanto, a seleção de 1% dos reprodutores com base nos seus valores genéticos para PLDC e PL305, leva a divergência entre os critérios de seleção, à exceção do quarto controle. A partir deste percentual, há maior proporção de touros coincidentes quanto às classes de PLDC e PL305.

Ao se comparar, a partir dos 5% melhores reprodutores, com base nos valores genéticos da PLDC2 a PLDC5 com a PL305, a porcentagem de indivíduos comuns foi acima de 80%. Nas PLDC8 a PLDC10, ocorreu menor coincidência de touros com PL305, sendo que as confiabilidades das estimativas dos valores genéticos dos touros também foram mais baixas, bem como as correlações de ordem no conjunto global de dados.

Assim, maiores ganhos pela seleção podem ser obtidos, quando se utilizam os primeiros controles como critério de seleção, pois o uso desta resposta correlacionada seria útil por poder ser medida preco cemente durante a lactação, reduzindo o intervalo de gerações. Esta alternativa poderá reduzir o custo e o número de controles leiteiros, uma vez que serão utilizados menos controles de produção de cada vaca por lactação.

Poderão ser incluídos, nas avaliações genéticas, animais com apenas uma medida de produção, minimizando vícios por descartes de lactações incompletas (Fimland, 1983; Ptak & Schaeffer, 1993). Martinez et al. (1998) informam que muitos animais são eliminados da avaliação do teste de progênie da raça Gir por apresentarem algum problema que não permite continuidade da lactação ou porque são vendidos antes de a encerrarem. A não inclusão destas lactações pode adulterar o resultado das avaliações.

Os resultados obtidos por Ptak & Schaffer (1993), Swalve (1995) e Ferreira (1999) mostraram ocorrer, também, altas correlações de ordem entre os valores genéticos dos touros em algumas PLDC e, alteração de hierarquia de touros nos diversos controles, em relação à classificação para PL305. Swalve (1995) verificou que a classificação dos touros não foi afetada drasticamente quando usou PLDC para avaliação.

O valor genético dos animais avaliados são expressos em termos relativos e são predições, portanto, podem variar em função da característica selecionada, do modelo utilizado e do número de informações disponíveis (Verneque, 1994 e Torres, 1998).

Vários autores argumentaram que a utilização de rebanho-ano-estação do controle em vez de rebanho-ano-estação do parto, quando se trabalha com PLDC, possibilita agrupar os controles sujeitos aos mesmos efeitos ambientes e comparar os animais com produções obtidas no mesmo período (Meyer et al., 1989; Ptak & Schaeffer, 1993; Swalve, 1995, Ferreira, 1999), maximizando a variância genética e minimizando a variância residual, obtendo melhores estimativas de herdabilidade e valores genéticos preditos mais precisos.

Os modelos que utilizam as PLDC poderão, por sua vez, inserir efeitos peculiares associados ao dia do controle, melhorando as estimativas dos parâmetros genéticos, reduzindo ações que seriam considerados como circunstanciais e que não poderiam constar do modelo com a PL305 (Trus & Buttazzoni, 1990; Pander et al., 1992; Stanton et al., 1992; Van Tassel et al., 1992; Reents et al., 1995; Wiggans & Goddard, 1997; Ferreira, 1999).

As altas correlações de ordem entre as avaliações dos touros revelam grande associação linear dos valores genéticos preditos usando as PLDC e PL305, indicando que a ordem dos animais nas diferentes avaliações tende a ser a mesma, no conjunto global de dados.

Os touros em comum que seriam selecionados pelos valores genéticos preditos das PLDC2 a PLDC5, foram acima de 80%, a partir de 5% dos melhores classificados na PL305, podendo, quando se desejar, usar dos resultados destas avaliações genéticas para seleção de animais.

Novos estudos necessitam ser realizados, incluindo vacas com lactações incompletas, agregando e conjugando múltiplas PLDC nas análises com base na PL305, testando modelos mais refinados que incluam fatores específicos aos controles, visando melhorar a precisão das avaliações das PLDC, pois na PL305 estes não constariam do modelo.

ALBUQUERQUE, L.G. Parâmetros genéticos das produções de leite parciais acumuladas de um rebanho Gir no Estado de São Paulo. Ribeirão Preto: Universidade de São Paulo, 1984. 79p. Dissertação (Mestrado em Genética) - Universidade de São Paulo, 1984.

ALI, T.E.; SCHAEFFER, L.R. Accounting for covariances among test day milk yields in dairy cows. Canadian Journal of Animal Science, v.67, n.3, p.637-644, 1987.

BOLDMAN, K.G.; KRIESE, L.A.; Van VLECK, L.D. et al. A Manual for use of MTDFREML. A set of programs to obtain estimative of variances and covariances [DRAFT]. Beltsville: Department of Agriculture, Agricultural Research, 1995. 125p.

BOWMAN, J.C. Introdução ao melhoramento genético animal. São Paulo: Universidade de São Paulo, 1981. 87p.

BRASIL - Ministério da Agricultura, Abastecimento e Reforma Agrária. Normas técnicas para execução do serviço de controle leiteiro em bovídeos. DOU - Diário Oficial da República Federativa do Brasil, Brasília, DF, 15 out. 1986, n.195, Seção I, p.15.532 a 15.535.

DANELL, B. Studies on lactation yield and individual test day yields of Swedish dairy cows. II. Estimates of genetic and phenotypic parameters. Acta Agriculture Scandinavia, v.32, n.1, p.83-91, 1982.

FERREIRA, W.J. Parâmetros genéticos para produção de leite no dia do controle de vacas da raça Holandesa. Viçosa, MG: Universidade Federal de Viçosa, 1999. 103p. Dissertação (Mestrado em Zootecnia) - Universidade Federal de Viçosa, 1999.

FIMLAND, E. Some properties of using the additive genetic relationship in the mixed model technique. Zeitschrift fuer Tierzuchtungund Zuechtungsbiologie, v.100, p.361-374, 1983.

FIRAT, M.Z.; THEOBALD, C.M.; THOMPSON, R. Univariate analysis of test day milk yields of British Holstein-Friesian heifers using Gibbs Sampling. Acta Agriculture Scandinavia, v.47, n.4, p.213-220, 1997a.

FIRAT, M.Z.; THEOBALD, C.M.; THOMPSON, R. Multivariate analysis of test day milk yields of British Holstein-Friesian heifers using Gibbs sampling. Acta Agriculture Scandinavia, v.47, n.4, p.221-220, 1997b.

GADINI, C.H. Influências genéticas e ambientes sobre a curva de lactação de um rebanho Gir. Ribeirão Preto: Universidade de São Paulo, 1985. 96p. Dissertação (Mestrado em Genética) - Universidade de São Paulo, 1985.

GADINI, C.H. Genetic evaluations of test day production traits and somatic cell scores. Lincoln: University of Nebraska, 1997. 91p. Dissertation (Doctor of Philosophy) - University of Nebraska, 1997.

GODDARD, M.E. Milk recording data needed for genetic evaluation. In: Biennal Session of the International ComitTee for Animal Recording, 27., paris. Proceedings... Paris: ICAR/EEAP, 1991. p.161-166.

KEOWN, J.F.; Van VLECK, L.D. Selection on test day fat percentage and milk production. Journal of Dairy Science, v.54, n.2, p.199-203, 1971.

KETTUNEN, J.F.; MANTYSAARI, I.; STRANDÉN, I. et al. Estimation of genetic parameters for first lactation test day milk production using random regression models. In: WORLD CONGRESS ON GENETICS APLIED TO LIVESTOCK PRODUCTION, 6., 1998, Armidale. Proceedings... Armidale: U.N.E., 1998. v.23, p.307-310.

LEDIC, I.L. Aplicação prática das provas de progênie do Gir leiteiro a nível de rebanho Gir puro e em vacas mestiças. In: SEMINÁRIO INTERNACIONAL DE MONTERIA DE GANADO DE DOBLE PROPÓSITO, GYR-LECHERO Y BUFALOS, 1., 1996, Monteria. Memorias... Monteria: PRODESA, 1996. p.94-105

LOPES, P.S.; MARTINS, E.N.; SILVA, M.A. et al. Estimação de componentes de variância. Viçosa, MG: Universidade Federal de Viçosa, 1993. 61p.

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 em Zootecnia) - Universidade Estadual Paulista, 1997.

MARTINEZ, M.L.; VERNEQUE, R.S. Programa nacional de melhoramento genético. Balde branco, n. 439, 2001 (Encarte técnico: Produção & Rusticidade - Gir Leiteiro, a solução para os trópicos).

MARTINEZ, M.L.; VERNEQUE, R.S.; TEODORO, R.L. Informativo Gir leiteiro. Juiz de Fora: EMBRAPA, 1998. 11p.

MELLO, C.M.R.; GONÇALVES, T.M.; VERNEQUE, R.S. et al. Avaliação genética de touros da raça Gir usando produções em lactações completas ou parciais. II. Correlações e coincidência de ordem. In: SIMPÓSIO NACIONAL DE MELHORAMENTO ANIMAL, 3., 2000, Belo Horizonte. Anais... Belo Horizonte: FEPMVZ, 2000. p.331-332

MEYER, K.; GRASER, H.; HAMMOND, K. Estimates of genetic parameters for first lactation test day production of Australian black and white cows. Livestock Production Science, v.21, n.3, p.177-199, 1989.

MICROSOFT. Microsoft® Corporation Advanced Software, Inc. Santa Rosa, 1997.

PANDER, B.L.; HILL, W.G.; THOMPSON, R. Genetic parameters of test day records of British Holstein-Friesian heifers. Animal Production, v.55, n.1, p.11-21, 1992.

PATTERSON, H.D.; THOMPSON, R. Recovery of interblock information when block sizes are unequal. Biometrika, v.58, n.2, p.545-554, 1971.

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.

RIBAS, M.; PEREZ, B. Monthly test day milk records and yield at 244 days. II. Genetic parameters in first lactation. Cuban Journal of Agriculture Science, v.24, n.2, p.139-144, 1990.

SAMPAIO, I.B.M. Estatística aplicada à experimentação animal. Belo Horizonte: FEPMVZ, 1998. 221p.

SAS INSTITUTE. SAS/STAT^® user’s guide. version 6. 4.ed. Cary: 1990. 1290p.

SCHAEFFER, L.R.; MINDER, C.E.; McMILLAN, U. et al. Non-linear techniques for predicting 305-day lactation production of Holsteins and Jerseys. Journal of Dairy Science, n.1, v.60, p.1636-1644, 1977.

SIKKA, A.K.; TANEJA, V.K. Inheritage of part lactation records in Sahiwal cattle. Indian Journal of Animal Science, v.51, n.9, p.809-811, 1981.

STRABEL, T.; SZWACZKOWSKI, T. Additive genetic and permanent environmental variance components for test day milk traits in Black-White cattle. Livestock Production Science, v.48, n.2, p.91-98, 1997.

TIJANI, A.; WIGGANS, G.R.; Van TASSEL, C.P. et al. Use of (co)variance function to describe (co)variances for test day yield. Journal of Dairy Science, v.82, n.1, p.226, 1999.

TORRES, R.A. Efeito da heterogeneidade de variância na avaliação genética de bovinos da raça Holandesa no Brasil. Belo Horizonte: Universidade Federal de Minas Gerais, 1998. 124p. Tese (Doutorado em Melhoramento Genético Animal) - Escola de Veterinária, Universidade Federal de Minas Gerais, 1998.

Van TASSEL, C.P.; QUASS, R.L.; EVERETT, R.W. Parameter estimates for 305-day ME records and 305-day test day residual records. Journal of Dairy Science, v.75, p.251, 1992. (suppl. 1)

Van VLECK, L.D.; HENDERSON, C.R. Use of part lactation records in sire evaluation. Journal of Dairy Science, v.44, n.8, p.1511-1518, 1961.

VERNEQUE, R.S. Procedimentos numéricos e estimação de componentes de covariância em análise multivariada pelo método da máxima verossimilhança restrita: modelos mistos aplicados ao melhoramento animal. Piracicaba: Escola Superior de Agricultura "Luiz de Queiroz", 1994. 155p. Tese (Doutorado em Estatística e Experimentação Animal) - Escola Superior de Agricultura "Luiz de Queiroz", 1994.

WILMINK, J.B.M. Efficiency of selection for different cumulative milk, fat and protein yields in first lactation. Livestock Production Science, v.17, n.3, p.211-224, 1987.

¹ Parte da tese de Doutorado do primeiro autor apresentada à FCAV/UNESP de Jaboticabal. ² Pesquisador da EPAMIG, Caixa Postal 351, 38001-970, Uberaba, MG. ³ Pesquisadores da EMBRAPA Gado de Leite, R. Eugênio do Nascimento 610, 36038-330, Juiz de Fora, MG. ⁴ Pesquisadora do IZ, Av. Bandeirantes 2419, 14001-970, Ribeirão Preto, SP. ⁵ Professores da FCAV/UNESP, Via Prof. Paulo Donato Castellane s/n, 14884-900, Jaboticabal, SP. ⁶ Pesquisador da EPAMIG, Caixa Postal 351, 38001-970, Uberaba, MG.