Os trabalhos têm mostrado que, de maneira geral, as estimativas de herdabilidade das produções de lactações parciais são semelhantes àquelas obtidas para a produção até 305 dias (TANDON e HARVEY, 1984; GENGLER, KEOWN e VAN VLECK, 1995; e WILMINK, 1987).

Produções de leite de lactações em progresso, projetadas para 305 dias, são altamente correlacionadas com a produção aos 305 dias de lactação. Estimativas de correlações genética de 0,85; 0,96; 0,99; e 1,00 para a produção aos 60, 120, 180 e 240 dias de lactação, projetadas para 305 dias, foram reportadas por WILMINK (1987), analisando as primeiras lactações de vacas holandesas. Estimativas similares foram obtidas por VAN VLECK (1964), TANDON e HARVEY (1984), SETHI e JAIN (1993), GENGLER, KEOWN e VAN VLECK (1995) e ZAHED, KHALIL e SOLIMAN (1997).

Correlações fenotípicas foram estimadas de registros mensais de várias lactações de vacas da raça Holandês por KEOWN e VAN VLECK (1971), que observaram estimativas variando de 0,65 a 0,90, sendo os maiores valores observados entre 120 e 150 dias de lactação. Vários outros estudos constataram altas estimativas de correlações fenotípicas das produções de leite no dia do controle com a produção padronizada aos 305 dias (WILMINK, 1987; PANDER, HILL e THOMPSON, 1992; e PANDER e HILL, 1993) e aos 244 dias de lactação (RIBAS, PEREZ e GUZMAN, 1994).

Os resultados acima demostram que o momento mais adequado para projeção da lactação se encontra entre o quinto e o sexto controle, parecendo, assim, não serem necessários mais que 150 dias de produção para que a lactação possa ser usada na avaliação do mérito genético dos animais.

O objetivo deste trabalho foi estimar parâmetros genéticos e fenotípicos para produção de leite até 305 dias e produções parciais projetadas.

Os dados utilizados neste estudo são provenientes de rebanhos da raça Gir, que participam do programa de melhoramento dessa raça, coordenado pela Embrapa Gado de Leite, em convênio com a Associação Brasileira dos Criadores de Gir Leiteiro (ABCGIL).

Após eliminar informações indesejáveis ao estudo (lactações consideradas anormais; lactações com produção zero em controles mensais; lactações com duração da lactação menor que 90 dias; lactações cujas vacas apresentaram idade ao primeiro parto menor que 18 meses; vacas que não apresentavam a primeira lactação; lactações de ordem de parto superior a cinco; e lactações de vacas sabidamente tratadas com Somatotropina Bovina Recombinante - bSTR) e considerando-se, ainda, um mínimo de cinco filhas por touro e de duas lactações por subclasse de rebanho-ano (RA), distribuídas nas duas estações de parição consideradas (1 = águas, outubro a março e 2 = seca, abril a setembro), analisaram-se 4595 lactações de 2254 vacas, filhas de 145 touros e de 1618 matrizes, distribuídas em 18 rebanhos, com partos ocorridos nos anos de 1980 a 1997.

As produções acumuladas foram calculadas até 90, 150, 210, 240 e 305 dias, usando o método do intervalo de teste, segundo EVERETT e CARTER (1968), e de uso oficial pelo Ministério da Agricultura (BRASIL, 1986).

A partir de 91, 151, 211 ou 241 dias de lactação, projetaram-se, aleatoriamente, 10, 30, 50 ou 70% das lactações, tanto para a duração da lactação observada, como para 305 dias, utilizando-se funções quadráticas logarítmicas obtidas por GONÇALVES, MARTINEZ e MILAGRES (1997).

Para classificação das vacas quanto ao nível de produção, utilizou-se a média da produção de leite na primeira lactação. Assim, o nível baixo de produção compreende todas as vacas que apresentaram produção média de leite, na primeira lactação, inferior a 2.230 kg de leite (média da produção de leite na primeira lactação), e nível alto de produção as demais.

Para projeção da lactação a partir de 91, 151 e 211 dias de lactação, utilizaram-se as seguintes equações:

, para baixo nível de produção da vaca, e , para alto nível de produção da vaca.

Para projeção da lactação a partir de 241 dias de lactação, utilizaram-se as seguintes equações:

, para baixo nível de produção da vaca, e , para alto nível de produção da vaca, sendo a estimativa da produção de leite diária, em quilogramas, e t o período de tempo desde o parto, em dias. A produção no período projetado é dada por: , em que é a produção estimada no dia i e t, como definido acima.

Além da característica produção de leite até 305 dias (P305d), foram geradas outras 32 estimativas de produção, conforme mostrado na Tabela 1.

Os dados foram preparados utilizando-se o pacote computacional "SAS^®" (STATISTICAL ANALYSIS SYSTEM) versão 6.12 (1996).

Para a análise utilizou-se o seguinte modelo estatístico:

em que

é o vetor das variáveis dependentes;

X é a matriz de incidência dos efeitos fixos;

é o vetor de efeitos fixos (rebanho-ano, época de parto, idade da vaca ao parto com termos linear e quadrático);

Z é a matriz de incidência dos efeitos genéticos diretos;

é o vetor dos efeitos aleatórios dos valores genéticos diretos do animal;

W é a matriz de incidência dos efeitos permanente de ambiente;

é o vetor dos efeitos permanentes de ambiente; e

é o vetor dos erros aleatórios associados a cada observação.

Assumiu-se que , , e possuem distribuição conjunta normal multivariada, como apresentado a seguir.

em que , e ou , e , quando se trata de análise de uma ou mais características, respectivamente, em que A é a matriz, de ordem m, de coeficientes de parentesco entre os animais avaliados; Z, a matriz, de ordem n x m, de incidência dos valores genéticos; W, a matriz, de ordem n x p, de incidência dos efeitos permanentes de ambiente; I, uma matriz identidade, de ordem n; I_p, uma matriz identidade, de ordem p; G₀, a matriz, de ordem q, de variâncias e covariâncias genética aditiva; P₀, a matriz, de ordem q, de variâncias e covariâncias dos efeitos permanentes de ambiente; R₀, a matriz, de ordem q, de variâncias e covariâncias dos resíduos; , a variância genética aditiva; , a variância de efeito permanente de ambiente; , a variância de efeito temporário de ambiente, sendo que n é o número total de observações; p, o número de indivíduos com observações; m, o número de indivíduos da população (com ou sem informações); e q, o número de características em análise.

A estimativa da média da produção de leite até 305 dias (P305d) foi de 2445,75 kg (13,35), com coeficiente de variação de 37,35%; as médias para as produções estimadas (PE) para duração da lactação observada variaram de 2381,98 (12,48) a 2458,44 (13,64) e as estimativas para as lactações projetadas (PE) para 305 dias, de 2407,43 (12,43) a 2544,68 (13,56).

Verificaram-se valores semelhantes para as médias da produção de leite das diferentes características consideradas, sendo as médias das produções estimadas pela projeção de 10 ou 30% das lactações mais semelhantes à média para P305d.

Constatou-se que, embora as PE sejam semelhantes à P305d, algumas lactações poderiam estar sendo super ou subestimadas, em valores expressivos, pela projeção, contudo estas diferenças se anularam na média, não tendo sido constatadas pelas diferenças entre as PE e P305d. A duração da lactação média nos dados em estudo foi de 292,41 (0,94) dias.

Pelas estimativas de herdabilidade (), repetibilidade () e seus respectivos erros-padrão (EP) obtidos neste estudo (Tabela 2), observa-se que os valores de e das PE, a partir de 211 ou 241 dias de lactação, foram mais semelhantes aos obtidos para P305d do que as estimativas obtidas com base nas PE a partir de 91 ou 151 dias de lactação.

As PE pela projeção de 10 ou 30% das lactações foram as que apresentaram estimativas de herdabilidade e repetibilidade mais semelhantes às obtidas para a P305d, resultado já esperado, tendo em vista que foram projetadas poucas lactações, sendo, assim insuficiente para alterar significativamente estas produções, quando comparadas à P305d.

Estimativas de 0,17 (0,047) e 0,50 (0,047) para a herdabilidade e a repetibilidade, respectivamente, foram obtidas para P305d, considerada como característica padrão, por ser utilizada em grande parte dos trabalhos de avaliação genética de bovinos leiteiros.

Estimativas mais elevadas obtidas para as PE pela projeção de 50 ou 70% das lactações podem ser explicadas, em parte, pelo fato de que maior número de registros novos foi gerado nestes casos, o que parece ter levado à redução na variância ambiental e aumento na variância genética aditiva.

As estimativas de obtidas neste trabalho assemelham-se às obtidas por TANDON e HARVEY (1984) e SETHI e JAIN (1993), para produções de lactações "truncadas" em diferentes estádios da lactação, e WILMINK (1987), em que todas as lactações foram projetadas a partir de diferentes estádios da lactação, para 305 dias.

Em animais da raça Gir, estimativas de 0,25 (0,047), 0,37 (0,056), 0,12 (0,04) e 0,24 (0,06) foram obtidas por RAMOS (1984), SOUZA (1991), VERNEQUE (1994) e VERNEQUE et al. (1998), respectivamente, sendo as duas primeiras pelo método dos quadrados mínimos e as demais pelo REML, sob modelo animal. Estimativa superior à deste trabalho, utilizando REML sob modelo de touro, foi descrita por PANDER et al. (1992), em registros de vacas da raça Holandês ( = 0,49), sendo que KEOWN e VAN VLECK (1971), AURAN (1976), DANELL (1982), WILMINK (1987), RIBAS e PEREZ (1990) e RIBAS et al. (1994) obtiveram valores mais próximos (0,25; 0,25; 0,30; 0,31; 0,29; e 0,29, respectivamente), também para vacas da raça Holandês.

A menor estimativa obtida no presente trabalho pode ter sido causada, especialmente, pela grande variabilidade nos registros de produção de leite da população estudada, sendo que a maior parte desta variabilidade, possivelmente, seja devida a fatores ambientais.

As estimativas de herdabilidade () e repetibilidade (t) e os respectivos erros-padrão (EP) da produção de leite até 305d e das produções de leite estimadas, obtidas por análise bivariada, mostraram que os valores obtidos para as PE a partir de 241 dias, foram mais semelhantes às obtidas para P305d.

As estimativas de herdabilidade e repetibilidade obtidas para P305d, oriundas das análises uni e bivariada, foram semelhantes, tendo em vista que as médias das estimativas de herdabilidade e repetibilidade, da produção de leite até 305 dias, obtidas por análise bivariada, foram de 0,19 (0,048) e 0,52 (0,047), enquanto as estimativas obtidas por análise univariada foram 0,17 (0,047) e 0,50 (0,047), respectivamente.

Estimativas de repetibilidade de 0,45 (0,05) para produção de leite até 90 dias; 0,46 (0,05) para produção até 150 dias; 0,47 (0,05) para produção até 210 dias; e 0,46 (0,05) para produção até 305 dias de lactação foram reportadas por VERNEQUE et al. (1998), obtidas pelo método REML em registros de vacas da raça Gir. Trabalhos anteriores de SOUZA (1991), pelo método dos quadrados mínimos, e VERNEQUE (1994), pelo método REML, relataram estimativas de 0,40 (0,019) e 0,46 para produção de leite até 305 dias, em registros de produção de leite de vacas da raça Gir.

É possível constatar, por intermédio das Tabelas 3 e 4, que as estimativas dos coeficientes de correlação genética entre as produções de leite estimadas pela projeção de uma porcentagem das lactações e P305d foram altas e positivas.

Quando as lactações foram projetadas para a duração da lactação observada, as estimativas das correlações genéticas entre PE e P305d variaram de 0,94 (0,483) a 1,00 (0,589) e, no caso em que as lactações foram projetadas para 305 dias de lactação, de 0,94 (0,466) a 1,00 (0,541). Pela semelhança entre as estimativas, é possível constatar que não houve diferença entre projetar para duração da lactação observada ou para 305 dias, o que pode ser explicado pelo fato de a média de duração das lactações, 292,41 dias, ser próxima a 305 dias.

Observa-se que as estimativas obtidas para a correlação entre produção de leite até 305 dias (P305d) e PE pela projeção a partir de 211 ou 241 dias foram ligeiramente mais altas (0,99 a 1) que as obtidas entre P305d e PE pela projeção a partir de 91 ou 151 dias (0,94 a 1). Estas diferenças foram maiores entre P305d e PE pela projeção de 70% das lactações a partir de 91 dias de lactação. Pode-se constatar tendência de valores pouco menores (0,94 a 0,97) nas estimativas das correlações, com o aumento da porcentagem de lactações projetadas e do período de projeção (projeção de 70% das lactações a partir de 91 ou 151 dias de lactação).

As estimativas obtidas no presente trabalho são similares às encontradas na literatura: 0,85; 0,96; 0,99 e 1 entre P305d e as produções de lactações projetadas a partir de 60, 120, 180 e 240 dias de lactação, respectivamente (WILMINK, 1987).

Estimativas também elevadas foram obtidas por TANDON e HARVEY (1984), SETHI e JAIN (1993) e VERNEQUE et al. (1998) entre a produção de leite de lactações parciais e a produção de leite até 305 dias.

Tal como as estimativas dos coeficientes de correlação genética, as estimativas dos coeficientes de correlação fenotípica foram altas e positivas (Tabelas 3 e 4).

Semelhantemente, também foram obtidas as menores estimativas entre P305d e as PE, por meio da projeção de 70% das lactações a partir de 91 dias de lactação (0,82 - entre P305d e P090d70%DL e 0,73 - entre P305d e P090d70%).

Nota-se, ainda, que quanto mais tardia for a projeção das lactações e menor a porcentagem de lactações projetadas, maiores serão as estimativas dos coeficientes de correlação, sugerindo, neste caso, resultados mais confiáveis nos testes de avaliação genética, quando lactações parciais são também utilizadas.

Pode-se observar que, à semelhança do que ocorreu com as estimativas de correlação genética, no caso da correlação fenotípica, as estimativas tendem a ser maiores entre as PE pela projeção de uma porcentagem das lactações, para a duração da lactação observada, em comparação às estimativas obtidas das PE pela projeção das lactações para 305 dias.

A literatura é escassa em estudos nos quais é estimada a correlação fenotípica entre produção de leite de lactações parciais e P305d. Contudo, os poucos relatos encontrados na literatura mostram valores de correlação relativamente altos. Estimativas variando de 0,64 a 0,99 foram obtidas por AURAN (1976), entre produções de leite acumuladas aos 30, 60, 90, ..., 300 dias e a produção de leite aos 305 dias de lactação. Semelhantemente, estimativas variando de 0,69 a 0,97; 0,65 a 0,99; e 0,60 a 0,97 foram obtidas, respectivamente, por DANELL (1982), PANDER e HILL (1993) e RIBAS et al. (1994), entre as produções de leite acumuladas e a produção aos 305 dias ou aos 244 dias de lactação, no caso do último estudo. Estas estimativas são inferiores às do presente estudo, o que se deve, possivelmente, ao fato destes autores não terem projetado as lactações, e tão somente utilizado as produções acumuladas nos diferentes estágios da lactação estudados.

Os coeficientes de correlação, tanto do efeito permanente () como temporário () de ambiente, foram próximos a um, seguindo tendência similar à obtida para os coeficientes de correlação genética e fenotípica. Entretanto, é possível constatar que as correlações do efeito temporário de ambiente, embora mais precisas (com menores erros-padrão), foram inferiores às demais.

As menores estimativas de correlação dos efeitos permanente e temporário de ambiente foram obtidas para P090d70%, sendo as mesmas de 0,78 (0,205) e 0,60 (0,038), respectivamente.

Observa-se que as estimativas são mais altas entre PE pela projeção de menor porcentagem de lactações e P305d. Semelhantemente, as estimativas entre as PE pela projeção das lactações a partir de 211 ou 241 dias de lactação e a produção de leite até 305 dias foram mais altas.

As altas estimativas de correlação entre as produções obtidas pela projeção de uma porcentagem das lactações e a produção até 305 dias de lactação eram esperadas, já que a produção observada em parte da lactação está incluída na produção da lactação total.

As herdabilidades das produções estimadas pela projeção de 10% das lactações a partir de 91, 151 e 211 dias ou pela projeção de 10, 30, 50 ou 70% das lactações a partir de 241 dias foram semelhantes à obtida para a produção de leite até 305 dias.

Produções de leite estimadas pela projeção de 10, 30, 50 ou 70% das lactações a partir de 91, 151, 211 ou 241 dias são altamente correlacionadas com a produção de leite até 305 dias.

AURAN, T. 1976. Studies on monthly and cumulative monthly milk yield records. Acta Agric. Scand., (26):3-9.

BRASIL. Ministério da Agricultura. Secretaria Nacional de Produção Agropecuária. Portaria número 45 de 10 de outubro 1986. Normas Técnicas para Execução do Serviço de Controle Leiteiro em Bovídeos. Diário Oficial [da República Federativa do Brasil], Brasília, 124(198):15533-15535, 15 de out. 1986. Seção 1.

BOLDMAN, K.G., KRIESE, L.A., VAN VLECK, L.D. et al. 1995. A manual for of MTDFREML: a set of programs to obtain estimates of variances and covariances (DRAFT). Washington USDA, Agricultural Research Service. 115p.

DANELL, B. 1982. Studies on lactation yield and individual test-day yields of Swedish dairy cow. II. Estimates of genetic and phenotypic parameters. Acta Agric. Scand., (32):103-114.

EVERETT, R.W., CARTER, H.W. 1968. Accuracy of test interval method of calculating dairy herd improvement association records. J. Dairy Sci., 51(12):1936-1941.

GENGLER, N., KEOWH. J.K., VAN VLECK, L.D. 1995. Various persistency measures and relationships with total, partial and peak milk yields. Rev. Bras. Gen., 18(2):237-243.

GONÇALVES, T.M., MARTINEZ, M.L., MILAGRES, J.C. 1997. Curva de lactação na raça Gir. 2. Influência dos fatores de meio ambiente, estimativas de repetibilidade e herdabilidade para os parâmetros da curva de lactação quadrática logarítmica. R. Bras. Zootec., 26(1):88-97.

KEOWN, J.F., VAN VLECK, L.D. 1971. Selection on test-day fat percentage and milk production. J. Dairy Sci., 54(2):199-203.

PANDER, B.L., HILL, W.G. 1993. Genetic evaluation of lactation yield from test day records on incomplete lactation. Lvstck. Prod. Sci., 37(1-2):23-36.

PANDER, B.L., HILL, W.G., THOMPSON, R. 1992. Genetic parameters of test day records of British Holstein - Friesian heifers. Anim. Prod., 55(1):11-21.

RAMOS, A. A. Estudo das características reprodutivas e produtivas de zebuínos leiteiros da raça Gir nos trópicos. Botucatu: UNESP, 1984. 224p. Tese (Livre Docência na Disciplina de Bovinocultura de Leite) - Universidade Estadual Paulista, 1984.

RIBAS, M., PEREZ, B. 1990. Monthly test day milk records and yield at 244 days. II. Genetic parameters in first lactation . Cuban J. Agric. Sci., 24(2):129-144.

RIBAS, M., PEREZ, B., GUZMAN, G. 1994. Monthly test-day yields and the extension of lactations in Cuba. Cuban J. Agric. Sci., 28(2):125-139.

SETHI, I.C., JAIN, J.P. 1993. Sire evaluation on partial records in dairy cattle. Ind. J. Anim. Sci., 63(8):869-872.

SOUZA, E.M. Estudos de fatores genéticos e no meio que influenciam o desempenho de rebanhos de Gir leiteiro. Viçosa: UFV, 1991. 103p. Dissertação (Mestrado em Melhoramento Genético) - Universidade Federal de Viçosa, 1991.

TANDON, P.K., HARVEY, W.R. 1984. Best linear unbiased prediction of sire breeding values from part lactations of daughters. J. Dairy Sci., 67(10):2399-2406.

VAN VLECK, L.D. p.1-31. Animal model for bull and cow evaluation. In: Large dairy herd management symposium. Gainesville, 1992.

VAN VLECK, L.D. 1964. Genetic parameters of five-month lactation records. J. Dairy Sci., 47(4):421-425.

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: ESALQ, 1994. 157p. Tese (Doutorado em Agronomia) - Escola Superior de Agricultura "Luis de Queiroz", 1994.

VERNEQUE, R.S., MARTINEZ, M.L., TEODORO, R.L. Avaliação genética de vacas e touros com base na produção de leite em diferentes estágios da lactação. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 35, 1998, Botucatu. Anais .... Botucatu: SBZ, 1998. p.255-257.

WILMINK, J.B.M. 1987. Efficiency of selection for different cumulative milk, fat, and protein yields in first lactation. Lvstck. Prod. Sci., (17):211-224.

² Doutorando em Zootecnia - USP.