Idade ao primeiro parto e intervalo de partos de cinco grupos genéticos Holandês x Gir no Brasil

Age at first calving and calving interval of five Holstein x Gir genetic groups in Brazil

MELHORAMENTO, GENÉTICA E REPRODUÇÃO

Idade ao primeiro parto e intervalo de partos de cinco grupos genéticos Holandês x Gir no Brasil

Age at first calving and calving interval of five Holstein x Gir genetic groups in Brazil

Olivardo Facó^I; Raimundo Nonato Braga Lôbo^II; Raimundo Martins Filho^III; Francisco de Assis Melo Lima^IV

^IPesquisador Embrapa Caprinos, Doutorando em Zootecnia – PDIZ/UFC; Fazenda Três Lagoas - Estrada Sobral Groaíras, Km 4, Caixa Postal D10 - CEP 62011-970 – Sobral - CE ([e-mail]; [e-mail])

^IIPesquisador Embrapa Caprinos – Sobral - CE ([e-mail]). Bolsista CNPq

^IIIProfessor colaborador PDIZ/UFC e Bolsista DCR/CNPq/FAPEPI/UFPI – Fortaleza - CE ([e-mail])

^IVProfessor do Departamento de Zootecnia - CCA – UFC ([e-mail])

RESUMO

Tendo em vista a grande importância dos bovinos mestiços Holandês x Gir (H x G) nos sistemas de produção de leite no Brasil e a íntima relação entre a reprodução e a produtividade destes sistemas, foram obtidos registros de partos e datas de nascimento junto à Associação Brasileira dos Criadores de Girolando, no intuito de avaliar o desempenho reprodutivo de vários grupos genéticos HxG. Os dados, num total de 2.600 registros de idade ao primeiro parto (IPP) e 3.763 de intervalos de partos (IDP), foram analisados por meio do procedimento GLM do pacote estatístico SAS (1999). As comparações do desempenho dos grupos genéticos foram feitas por meio de modelos classificatórios incluindo o efeito de grupo genético. Foram feitas estimativas dos efeitos genéticos, isto é, da diferença aditiva entre Holandês e Gir, do efeito de dominância e do efeito epistático aditivo x aditivo, utilizando-se modelos de regressão. O grupo meio-sangue apresentou desempenho superior, ou seja, menores IPP e IDP. Foi observada tendência de elevação do IDP com o aumento da proporção de genes da raça Holandesa. O efeito de dominância influenciou significativamente as duas características estudadas, enquanto o efeito de interação epistática do tipo aditivo x aditivo influenciou significativamente apenas a IPP. Os resultados indicaram que a perda de interações epistáticas favoráveis contribuiu para a alta IPP no grupo 5/8 Holandês.

Palavras-chave: cruzamentos, efeitos epistáticos, Girolando, heterose, desempenho reprodutivo

ABSTRACT

Reproduction efficiency of Holstein x Gir crossbred animals in Brazilian dairy farms were evaluated using data provided by the Brazilian Association of Girolando. A total of 2,600 records of age at first calving (IPP) and 3,763 records of calving intervals (IDP) were analyzed using the SAS^® GLM procedure. Estimates of genetics effects: additive difference between Holstein and Gir breeds, dominance and additive x additive epistatic effects were obtained by regression models including genetic groups to model the performance of crossbred groups. The genetic group 1/2 showed better performance: lower averages for IPP and IDP than the other groups. It was observed a trend of increasing IDP as the proportion of Holstein genes increased. The dominance effect influenced significantly both traits, while the additive x additive epistatic effects significantly influenced only IPP. The results indicated that loss of favorable epistatic interactions contributed to high IPP average in the 5/8 group.

Key Words: crossbreeding, epistatic effects, Girolando, heterosis, reproductive efficiency

Introdução

Na tentativa de melhorar a produtividade dos sistemas de produção de leite sob condições tropicais, tem-se utilizado em larga escala o cruzamento de raças zebuínas, que apresentam excelente adaptação às condições tropicais, com raças de origem européia, especializadas em produção de leite. Isso ocorre, geralmente, em razão de sérios problemas de adaptação dos animais puros de raças especializadas sob condições tropicais, como estresse térmico, baixa qualidade dos alimentos, manejo inadequado, parasitas, entre outros, que, em muitos casos, inviabilizam o sistema de produção.

Entre os fatores que mais influenciam a produtividade em rebanhos leiteiros, destaca-se a eficência reprodutiva, especialmente em rebanhos mestiços Europeu x Zebu, nos quais geralmente ocorrem lactações mais curtas nos animais dos grupos genéticos com menor proporção de genes da raça Holandesa (Facó et al., 2002). Estas lactações mais curtas precisam ser compensadas com menores intervalos de partos para que o período seco não seja excessivamente extenso.

Outro fator relacionado à eficência reprodutiva é a idade ao primeiro parto, indicativo de precocidade sexual e, portanto, de elevada importância econômica, uma vez que marca o início da vida produtiva de uma fêmea leiteira e influencia os custos de reposição das matrizes (Lemos et al., 1992).

No Brasil, a maior parte da produção de leite é obtida com a utilização de animais mestiços zebuínos, sobretudo os Holandês x Gir. Considerando a importância deste tipo racial no panorama da produção de leite nacional, em 1989, o Ministério da Agricultura, juntamente com as associações representativas, traçou as normas (Brasil, 1992) para a formação da Raça Girolando - Gado Leiteiro Tropical (5/8 Holandês + 3/8 Gir - bi-Mestiço). Para criar uma ampla base genética para a formação da raça sintética, a princípio, animais de vários grupos genéticos Holandês x Gir estão sendo controlados. Assim, a partir dos dados de controle de partos, fornecidos pela Associação Brasileira dos Criadores de Girolando, torna-se possível o estudo do desempenho reprodutivo dos vários grupos genéticos.

Este trabalho foi realizado com o objetivo de comparar o desempenho reprodutivo de cinco grupos genéticos Holandês x Gir, por meio da análise dos diversos fatores genéticos e de ambiente que influenciam a expressão das características idade ao primeiro parto (IPP) e intervalo de partos (IDP) e, a partir dos resultados, contribuir para maior conhecimento acerca das potencialidades dos diversos sistemas de cruzamento para melhorar o desempenho reprodutivo de rebanhos leiteiros sob condição tropical.

Material e Métodos

Foram obtidos registros de partos e datas de nascimento junto à Associação Brasileira dos Criadores de Girolando, a partir dos quais foi composto o arquivo de dados que possibilitou o cálculo das idades ao primeiro parto (IPP) e dos intervalos de parto (IDP).

As idades ao primeiro parto foram calculadas pela diferença em dias entre a de nascimento e a data do primeiro parto, enquanto os intervalos de partos foram obtidos pela diferença em dias entre dois partos consecutivos.

As análises descritivas dos dados foram realizadas por meio dos procedimentos estatísticos contidos no programa SAS (1999), verificando-se as restrições, as limitações e a consistência das informações.

O total de registros de idade ao primeiro parto e de intervalo de partos era inicialmente de 2.799 e 4.058, respectivamente. Foram excluídos das análises os registros considerados atípicos, ou seja, com IPP <600 ou >1460 e IDP <300. Foram ainda excluídos os registros provenientes de fazendas com menos de cinco observações, restando 2.600 e 3.763 registros de IPP e IDP, respectivamente. Os registros de IPP corresponderam a animais nascidos entre 1989 e 1995, enquanto os registros de IDP, a partos ocorridos entre 1991 e 1997.

Mais de 90% dos registros estudados foram provenientes de fazendas localizadas no estado de Minas Gerais.

Os dados foram analisados pelo método dos quadrados mínimos por meio do procedimento GLM (SAS, 1999). Para o estudo dos IDPs, foram definidos os seguintes modelos:

Modelo 1

Y_jklmno = m + g x_1j + d x_2j + gg x_3j + MP_k + AP_l + OP_m + FAZ_n + e_jklmno

Modelo 2

Y_jklmno = m + g x_1j + d x_2j + MP_k + AP_l + OP_m + FAZ_n + e_jklmno

Modelo 3

Y_jklmno = m + GG_j + MP_k + AP_l + OP_m + FAZ_n + e_jklmno

em que Y_jklmno = intervalo de partos; m = intercepto; GG_j = efeito do j^ésimo grupo genético (1/4, 1/2, 5/8, 3/4 e = 7/8 de genes holandês); MP_k = efeito do k^ésimo mês de parto; AP_l = efeito do l^ésimo ano de parto; OP_m = efeito da m^ésima ordem de parição (1ª. à 9ª.); FAZ_n = efeito da n^ésima fazenda; g = diferença aditiva entre as raças (Holandês e Gir); x_1j = proporção esperada de genes Holandês nos indivíduos do j^ésimo grupo genético; d = efeito de dominância; x_2j = proporção esperada de locos ocupados por um alelo de cada raça nos indivíduos do j^ésimo grupo genético; gg = efeito epistático aditivo x aditivo; x_3j = fração esperada de efeitos epistáticos do tipo aditivo x aditivo nas raças parentais recuperadas nos indivíduos do j^ésimo grupo genético; e_jklmno = efeito residual.

Todos os efeitos, exceto residual (e), foram considerados fixos.

Para a característica IPP, os modelos foram os mesmos utilizados para a IDP, com as seguintes modificações: (1) em vez de mês e ano de parto, foram considerados mês e ano de nascimento, respectivamente, e (2) não foi incluído o efeito da ordem de parição.

Os valores utilizados para x₁, x₂ e x₃ encontram-se na Tabela 1 e sua derivação foi realizada segundo Dickerson (1973) e Eisen (1989), citados por Lemos et al. (1992). Em função das diferentes formas de se obterem animais 5/8 e da indisponibilidade de informação acerca do genótipo dos pais, o valor de x₂ utilizado para este grupo genético foi de 0,57333 (média de 0,75; 0,50 e 0,47), conforme sugerido por Martinez et al. (1988).

A comparação do desempenho dos grupos genéticos para as duas características foi feita utilizando-se o modelo 3. Em razão do pequeno número de registros de animais dos grupos genéticos 7/8 e Holandês (>15/16), esses animais foram agrupados como grupo genético >7/8.

Segundo Robison et al. (1981), o uso de um modelo que contenha o efeito de grupo genético contabilizaria todos os efeitos genéticos, inclusive aqueles não-lineares e os epistáticos. Portanto, a comparação dos resultados obtidos pelo modelo 3 com aqueles obtidos pelos modelos de regressão 1 e 2 deve gerar informação acerca da importância dos efeitos genéticos não-lineares e epistáticos.

A qualidade do ajuste promovido pelos modelos 1 e 2 foi avaliada por meio do teste F, em que o numerador foi a diferença entre a soma de quadrados do resíduo dos modelos de regressão 1 ou 2 e a soma de quadrados do resíduo do modelo 3, dividida pela diferença entre os graus de liberdade do resíduo do modelo de regressão considerado e os graus de liberdade do resíduo do modelo 3, enquanto o denominador foi o quadrado médio residual do modelo 3 (Robison et al., 1981).

As estimativas para os efeitos genéticos (g, d e gg) foram produzidas a partir do modelo (1 ou 2) que apresentou melhor ajuste quando comparado ao modelo 3.

Resultados e Discussão

Todos os efeitos de ambiente incluídos nos modelos 1, 2 e 3 influenciaram significativamente (P<0,05) as características estudadas e não são discutidos aqui, pois foram incluídos apenas para aumentar a precisão das análises genéticas.

Na Tabela 2 são apresentadas as médias estimadas pelo método dos quadrados mínimos para cada grupo genético. O grupo genético 1/2 teve desempenho superior, isto é, menor IPP que os demais, exceto o grupo genético >7/8 de genes holandês. Ainda para IPP, o grupo 5/8 apresentou desempenho inferior aos demais grupos, exceto ao 1/4. As vacas do grupo 3/4 apresentaram IPPs intermediárias, tendo o desempenho inferior às meio-sangue, mas semelhante às do grupo >7/8.

Esses resultados confirmam aqueles encontrados por Lemos et al. (1992), em que as F1 tiveram primeiro parto mais jovens e as 5/8 em idade mais avançada. Rege et al. (1994), estudando o desempenho de animais mestiços Jersey x Ghana Shorthorn (Bos taurus) e Jersey x Sokoto Gudali (Bos indicus), observaram que animais F₁ apresentaram menores idades ao primeiro parto que os puros (Ghana Shorthorn ou Sokoto Gudali), porém, ao contrário do observado neste estudo, os animais F1 apresentaram IPPs semelhantes às dos animais produzidos por retrocruzamento (3/4 de genes da Raça Jersey).

Para a característica IDP, o grupo genético 1/2 apresentou desempenho superior, ou seja, menor IDP, aos grupos 3/4 e >7/8, porém, semelhante aos grupos 1/4 e 5/8. O grupo >7/8 apresentou desempenho inferior aos demais, exceto ao 5/8, enquanto o 3/4 não diferiu dos grupos 1/4 e 5/8. Verificou-se tendência de elevação do IDP à medida que aumentou a participação de genes da raça Holandesa. Esta tendência está provavelmente ligada à tendência de maior produção de leite dos grupos genéticos com maior proporção de genes holandeses (Facó et al., 2002).

Madalena et al. (1990) observaram diferenças de IDP entre seis grupos genéticos Holandês x Guzerá apenas no primeiro intervalo de partos, em que os animais 5/8 apresentaram menores IDP's. No entanto, esses autores destacaram que esse bom desempenho reprodutivo dos 5/8 não foi suficiente para compensar a baixa produção de leite exibida por esses animais.

Por outro lado, Guimarães et al. (2002) não observaram diferença significativa quanto IDP de vários grupos genéticos Holandês x Gir, variando desde o Gir puro até o Holandês puro por cruza. Também Lemos et al. (1997), comparando estratégias de cruzamentos, como F₁, 5/8 bi-mestiça, rotacionado com repetição do Holandês e absorção pelo Holandês, e Madalena et al. (1983), comparando o desempenho de animais 3/4, 7/8 e Holandês puro, observaram IDPs semelhantes.

Rege et al. (1994), estudando o desempenho de animais mestiços Jersey x Ghana Shorthorn, observaram que animais F₁ apresentaram menores intervalos de partos que animais puros Ghana Shorthorn e 3/4 Jersey, mas não observaram diferenças significativas entre animais com proporções genéticas semelhantes no cruzamento Jersey x Sokoto Gudali.

Na Tabela 3 estão apresentadas as comparações dos ajustes dos modelos epistático e aditivo-dominante em relação àquele obtido pelo modelo classificatório para as duas características.

Verificou-se que para a IPP o modelo aditivo-dominante não ajustou adequadamente os dados, uma vez que a variação adicional explicada pelo modelo classificatório foi estatisticamente significativa (P<0,01). O modelo epistático, no entanto, garantiu um ajuste tão eficiente quanto o modelo classificatório que inclui o efeito de grupo genético.

Lemos et al. (1992) encontraram resultado semelhante sob baixo nível de manejo, porém, sob alto nível manejo, verificaram que tanto o modelo aditivo-dominante quanto o epistático ajustaram os dados tão bem quanto o modelo classificatório.

Para a característica IDP, tanto o modelo aditivo-dominante quanto o modelo epistático promoveram ajustes semelhantes àquele obtido pelo modelo classificatório. No entanto, o modelo aditivo-dominante apresentou um ajuste mais próximo do modelo classificatório, evidenciado pelo menor valor de F e pelos menores erros-padrão das estimativas dos efeitos genéticos (Tabelas 3 e 4).

Na Tabela 4 são apresentadas estimativas da diferença genética aditiva entre as raças Holandês e Gir (g), dos efeitos de dominância (d) e de interações epistáticas do tipo aditivo x aditivo (gg). Para a IPP, estas estimativas foram geradas pelo modelo epistático, uma vez que o aditivo-dominante não promoveu ajuste adequado. Para a IDP, no entanto, foi utilizado o modelo aditivo-dominante, em razão do menor erro-padrão das estimativas e do menor valor de F obtido na comparação com o modelo classificatório.

As interações epistáticas aditiva x aditiva influenciaram significativamente (P<0,001) a expressão da idade ao primeiro parto, evidenciando a importância dos efeitos epistáticos para essa característica nas condições deste estudo. Este fato, somado à inadequação do modelo aditivo-dominante, é indicativo de que o fraco desempenho do grupo genético 5/8 em relação à IPP decorreu não apenas da perda de heterose em relação ao F₁, mas também à perda de interações epistáticas favoráveis.

Os sinais negativos das estimativas dos efeitos genéticos (d e gg) para a IPP indicam que tanto o efeito de dominância (d) quanto das interações epistáticas aditiva x aditiva (gg) foram benéficos no sentido de contribuírem para a redução da idade ao primeiro parto, corroborando os resultados obtidos por Lemos et al. (1992) sob baixo nível de manejo.

Quanto ao IDP, apenas o efeito de dominância (d) foi significativo (P<0,05), contribuindo para a redução deste intervalo. Embora a diferença genética aditiva entre as raças não tenha apresentado significância, verificou-se uma tendência de elevação do IDP, à medida que se elevou a participação dos genes da raça Holandesa, o que está reforçado pelo sinal positivo da estimativa de g.

Mackinnon et al. (1996) também não encontraram significância dos efeitos de interações epistáticas para o IDP.

Nas condições deste estudo, verifica-se que a heterose para a IPP resultou tanto dos efeitos de dominância quanto dos de interações epistáticas, enquanto, para o IDP, os efeitos de dominância foram os principais responsáveis pela heterose.

Conclusões

A adequação do modelo genético de regressão aditivo x dominante deve ser sempre testada, pois esse modelo nem sempre se ajusta de maneira conveniente para os efeitos genéticos, ou seja, nem sempre os efeitos de interação epistáticas, dentre outros efeitos genéticos, são negligenciáveis.

Os resultados observados indicam que a utilização de cruzamentos para produção de animais F1 constitui-se em uma importante alternativa para elevar a precocidade sexual e a fertilidade dos rebanhos leiteiros brasileiros.

O baixo desempenho dos animais 5/8 Holandês demonstra a necessidade de um rigoroso e bem conduzido processo seletivo durante a formação da raça Girolando. Embora sejam características cujas herdabilidades variam de moderada a baixa magnitude, sendo esperado pouco ganho genético, sua elevada importância econômica justifica adequada atenção.

Agradecimento

À Associação Brasileira dos Criadores de Girolando, que cedeu atenciosamente os dados para este estudo e esteve sempre disponível para esclarecimentos.

Literatura Citada

Recebido em: 04/02/04

Aceito em: 15/07/05

• BRASIL. Departamento Nacional de Produção Agropecuária – Coordenação de Produção Animal. Normas para formação da raça Girolando Brasília: DF, 1992. 31p.
• FACÓ, O.; LÔBO, R.N.B.; MARTINS FILHO, R. et al. Análise do desempenho produtivo de diversos grupos genéticos Holandês x Gir no Brasil. Revista Brasileira de Zootecnia, v.31, n.5, p.1944-1952, 2002.
• GUIMARÃES, J.D.; ALVES, N.G.; COSTA, E.P. et al. Eficiência reprodutiva e produtiva em vacas das raças Gir, Holandês e cruzadas Holandês x Zebu. Revista Brasileira de Zootecnia, v.31, n.2, p.641-647, 2002.
• LEMOS, A.M.; MADALENA, F.E.; TEODORO, R.L. et al. Comparative Performance of six Holstein-Friesian x Guzera grades in Brazil. 5. Age at first calving. Revista Brasileira de Genética, v.15, n.1, p.73-83, 1992.
• LEMOS, A.M.; VERNEQUE, R.S.; TEODORO, R.L. et al. Efeito da estratégia de cruzamentos sobre características produtivas e reprodutivas em vacas do sistema mestiço do CNPGL-EMBRAPA. Revista Brasileira de Zootecnia, v.26, n.4, p.704-708, 1997.
• MACKINNON, M.J.; THORPE, W.; BAKER, R.L. Sources of genetic variation for milk production in a crossbred herd in the tropics. Animal Science, v.62, p.5-16, 1996.
• MADALENA, F.E.; LEMOS, A.M.; TEODORO, R.L. et al. Dairy production and reproduction in Holstein-Friesian and Guzera crosses. Journal of Dairy Science, v.73, p.1872-1886, 1990.
• MADALENA, F.E.; VALENTE, J.; TEODORO, R.L. et al. Produção de leite e intervalo entre partos de vacas HPB e mestiças HPB:Gir num alto nível de manejo. Pesquisa Agropecuária Brasileira, v.18, n.2, p.195-200, 1983.
• MARTINEZ, M.L.; LEE, A.J.; LIN, C.Y. Age and Zebu-Holstein additive and heterotic effects on lactation performance and reproduction in Brazil. Journal of Dairy Science, v.71, n.3, p.800-808, 1988.
• REGE, J.E.O.; ABOAGYE, G.S.; AKAH, S. et al. Crossbreeding Jersey with Ghana Shorthorn and Sokoto Gudali cattle in a tropical environment: additive and heterotic effects for milk production, reproduction and calf growth traits. Animal Production, v.59, p.21-29, 1994.
• ROBISON, O.W.; McDANIEL, B.T.; RINCON, E.J. Estimation of direct and maternal additive and heterotic effects from crossbreeding experiments in animals. Journal of Animal Science, v.52, n.1, p.44-50, 1981.
• STATISTICAL ANALYSIS SYSTEMS - SAS. User's guide Cary: 1999.

Datas de Publicação

• Publicação nesta coleção
  13 Jan 2006
• Data do Fascículo
  Dez 2005

Histórico

• Aceito
  15 Jul 2005
• Recebido
  04 Fev 2004