Análise do desempenho produtivo de diversos grupos genéticos Holandês x Gir no Brasil

Analysis of productive performance of different Holstein x Gir genetic groups in Brazil

Análise do desempenho produtivo de diversos grupos genéticos Holandês x Gir no Brasil¹

Analysis of productive performance of different Holstein x Gir genetic groups in Brazil

Olivardo FacóI, ²; Raimundo Nonato Braga LôboII, ³; Raimundo Martins FilhoIII, ⁴; Arlindo de Alencar Araripe MouraIV, ⁵

^IMestre em Zootecnia, CCA, UFC; Rua: Estado do Rio, 55/323. Bairro: Bela Vista, CEP: 60441-150 Fortaleza, CE

^IIDS., Pesquisador da Embrapa Caprinos, Sobral, CE

IIIDS., Professor do Departamento de Zootecnia, CCA, UFC

IVDS., Professor do Departamento de Zootecnia, CCA, UFC

^{Endereço para correspondência}

RESUMO

Foram obtidos dados de controle leiteiro mensal junto à Associação Brasileira dos Criadores de Girolando, a partir dos quais foram calculadas as produções de leite por lactação (PL), as produções de leite até os 305 dias de lactação (PL305) e as durações das lactações (DL). Foram analisadas 3.574 lactações, das quais 10,8% foram registradas sob regime alimentar extensivo (RAEX), 67,9% sob regime alimentar semi-intensivo (RASI) e 21,3% sob regime alimentar intensivo (RAI). Os dados foram analisados pelo método dos quadrados mínimos por meio do procedimento GLM (SAS, 1996). Análises preliminares indicaram forte interação entre grupo genético e regime alimentar. Em função disso, a comparação do desempenho dos diversos grupos genéticos foi feita separadamente para cada regime alimentar. Foram também estudados os efeitos genéticos da diferença genética aditiva (g) entre as raças Holandesa e Gir, da dominância (d) e da interação epistática do tipo aditiva x aditiva (gg). Verificou-se que nem sempre os efeitos das interações epistáticas são negligenciáveis. Não foi verificado qualquer benefício em elevar a proporção de genes da Raça Holandesa sob condições de ambiente hostis. Por outro lado, para manejos mais aprimorados, tal elevação deu indícios de ser benéfica para o aumento da produção.

Palavras-chave: cruzamentos, duração da lactação, heterose, interações epistáticas, produção de leite

ABSTRACT

Total milk yield per lactation (PL), milk yield in 305 days (PL305) and lactation length (DL) were estimated based on records of the Brazilian Association of Girolando Breeders (3,574 lactations). Milk production was recorded monthly and 10.8% of the herds were raised under an extensive feeding system, 67.9% under a semi-intensive feeding system and 21.3% under an intensive feeding system. Data were analyzed through the least square method by GLM procedure (SAS, 1996). Preliminary analysis indicated a high interaction between genetic group and feeding system. Thus, the performance of genetic groups was evaluated within each feeding system. Genetic effects of additive difference between (g) Holstein and Gir breeds, dominance (d) and additive x additive epistatic interactions (gg) were also studied. It was observed that, in some cases, the effects of epistatic interactions should be taken into account. Moreover, increasing the proportion of Holstein-Friesian genes did not have any significant effect on milk yield when herds were raised under poor management conditions. However, as such conditions improved, crosses with higher proportion of Holstein tended to show better performance.

Key Words: crossbreeding, epistatic interactions, heterosis, lactation length, milk yield

Introdução

A produtividade dos sistemas de produção de leite em áreas de clima tropical é caracteristicamente baixa em todo o mundo, quando comparada aos sistemas de clima temperado. Tal ineficiência se dá pelos aspectos, que se somam e interagem, seguintes: inadequados manejos nutricional, reprodutivo e sanitário, limitado potencial genético dos rebanhos e condições climáticas hostis.

Na tentativa de melhorar a produtividade destes sistemas, tem-se utilizado em larga escala o cruzamento de raças zebuínas (ou nativas adaptadas), que apresentam excelente adaptação às condições tropicais, com raças de origem européia especializadas para produção de leite. Isto ocorre, geralmente, devido aos sérios problemas de adaptação dos animais puros de raças especializadas sob condições tropicais (estresse térmico, baixa qualidade dos alimentos, manejo inadequado, parasitas etc), que, em muitos casos, inviabilizam o sistema de produção. Fundamentalmente, o principal objetivo deste tipo de cruzamento é utilizar-se da expressão da heterose e da complementaridade entre estes tipos zootécnicos para a obtenção de animais mais adaptados e produtivos sob tais condições.

No Brasil, a maior parte da produção de leite é oriunda da utilização de mestiços zebuínos. Dentro deste universo, ocupam posição de destaque os mestiços Holandês x Gir. Dada a importância deste tipo racial no panorama da produção de leite nacional, em 1989, o Ministério da Agricultura (Brasil, 1992), juntamente com as Associações representativas, traçaram as normas para a formação da raça Girolando - Gado Leiteiro Tropical (5/8 Holandês + 3/8 Gir-bi-Mestiço).

No intuito de 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. Desta forma, a partir dos dados de controle leiteiro, fornecidos pela Associação Brasileira dos Criadores de Girolando, foi possível estudar o desempenho produtivo dos vários grupos genéticos.

Assim, os objetivos deste trabalho foram comparar o desempenho produtivo de vários 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 produção de leite por lactação (PL), produção de leite até os 305 dias de lactação (PL305) e duração da lactação (DL), e, a partir dos resultados obtidos, contribuir para um maior conhecimento acerca das potencialidades dos diversos sistemas de cruzamentos para produção de leite sob condição tropical.

Material e Métodos

Foram obtidos dados de controle leiteiro e registros de partos junto à Associação Brasileira dos Criadores de Girolando, a partir dos quais foi composto o arquivo de dados que possibilitou o cálculo das produções de leite por lactação (PL), as produções de leite até os 305 dias de lactação (PL305) e as durações das lactações (DL).

As análises do estudo foram realizadas no Laboratório de Informática do Departamento de Zootecnia do Centro de Ciências Agrárias da Universidade Federal do Ceará, utilizando o programa Statistical Analysis System - SAS (1996).

A PL foi calculada multiplicando-se a média dos registros de controle leiteiro (kg / dia) pela duração da lactação (dias). O cálculo da PL305, quando a duração da lactação era igual ou inferior a 305 dias, foi igual ao da PL. Porém, quando as lactações excederam os 305 dias de duração, os registros de controle leiteiro foram truncados, sendo computados apenas os registros feitos dentro dos primeiros 305 dias de lactação, para efeito do cálculo da média destes registros, quando então esta média foi multiplicada por 305. A DL foi calculada pela diferença, em dias, entre as datas de parto e de secagem.

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

O total de lactações era inicialmente de 5.834, no entanto, foram impostas restrições de maneira que fossem analisadas apenas aquelas lactações com no máximo um registro de controle leiteiro mensal perdido, de modo que uma lactação de 300 dias de duração deveria ter pelo menos nove registros de controle leiteiro mensal. Além disso, foram ainda excluídos os registros de lactação que não continham a informação acerca do regime alimentar, bem como lactações provenientes de fazendas com menos de dez registros e de fazendas sob o regime alimentar intensivo com produção média inferior a 2.000 kg de leite por lactação (não se pode considerar intensivo um sistema de produção com média de produção de leite tão baixa). Por fim, as lactações sob regime de três ordenhas foram excluídas em função do pequeno número de registros sob tais condições, restando 3.574 lactações a serem analisadas, as quais foram registradas no período de 1992 a 1997, sendo todas sob regime de duas ordenhas.

Das 3.574 lactações analisadas, mais de 93,4% foram provenientes de fazendas localizadas no estado de Minas Gerais, sendo o restante de fazendas localizadas em Goiás (3,4%), São Paulo (2,1%), Mato Grosso do Sul (0,5%) e Alagoas (0,5%). 10,8% destas lactações foram registradas sob regime alimentar extensivo, 67,9% sob regime alimentar semi-intensivo e 21,3% sob regime alimentar intensivo.

Os dados foram analisados pelo método dos quadrados mínimos por meio do procedimento GLM (SAS, 1996). Com o objetivo de definir os efeitos fixos de ambiente que comporiam os modelos genéticos, foi estabelecido o modelo geral abaixo:

Y_ijklmno = m + RA_i + GG_j + RAGG_ij + MP_ik + AP_il + OP_im + FAZ_in + e_ijklmno

em que: Y_ijklmno = característica em questão (PL, PL305 ou DL); m = intercepto; RA_i = efeito do i^ésimo regime alimentar; GG_j = efeito do j^ésimo grupo genético; RAGG_ij = interação do regime alimentar i com o grupo genético j; MP_ik = efeito do k^ésimo mês de parto dentro do i^ésimo regime alimentar; AP_il = efeito do l^ésimo ano de parto dentro do i^ésimo regime alimentar; OP_im = efeito da m^ésima ordem de parição dentro do i^ésimo regime alimentar; FAZ_in = efeito da n^ésima fazenda dentro do i^ésimo regime alimentar; e_ijklmn = efeito residual.

Todos os efeitos, exceto o efeito residual (e), foram considerados fixos e influenciaram significativamente (P<0,05) as características em estudo. Verificou-se, neste ponto, o efeito altamente significativo (P<0,01) da interação entre Regime Alimentar e Grupo Genético. Em função disso, os grupos genéticos foram comparados dentro de cada regime alimentar. Da mesma forma, o estudo dos efeitos genéticos ficaria mais apropriado dentro de cada regime alimentar.

Assim, para o estudo da comparação do desempenho dos diversos grupos genéticos e dos efeitos genéticos em cada regime alimentar foram definidos os seguintes modelos:

Modelo 1

Yjklmno = 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 = característica em questão (PL, PL305 ou DL); m = intercepto; GG_j = efeito do j^ésimo grupo genético; 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; 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 aditivo x aditivo nas raças parentais que são recuperadas nos indivíduos do j^ésimo grupo genético; e_jklmno = efeito residual.

Todos os efeitos foram considerados fixos, exceto o efeito residual (e). Os valores utilizados para x₁, x₂ e x₃ estão listados na Tabela 1, sendo a derivação dos mesmos dada por Dickerson (1973) e Eisen (1989), citados por Lemos et al. (1992). Em função das diferentes formas de se obterem animais 5/8, 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 em cada regime alimentar foi procedida utilizando o modelo 3. Em função do pequeno número de registros de animais dos grupos genéticos 7/8 e Holandês (³ 15/16), estes 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. Desta forma, a comparação dos resultados obtidos a partir do modelo 3 com aqueles obtidos a partir dos 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 de testes F, em que: numerador do valor F = (soma de quadrados do resíduo do modelo 1 ou 2 menos soma de quadrados do resíduo do modelo 3) / (graus de liberdade do resíduo do modelo de regressão considerado menos graus de liberdade do resíduo do modelo 3); denominador do valor F = quadrado médio residual do modelo 3 (Robison et al., 1981).

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

Resultados e Discussão

Na Tabela 2 são apresentadas as estatísticas básicas das 3.574 lactações estudadas, enquanto, na Tabela 3, estas estatísticas são apresentadas separadamente para cada regime alimentar.

As médias observadas na Tabela 2 para PL e PL305 foram superiores àquelas observadas por Chaudhry et al. (1993), Freitas et al. (1995) e Freitas et al. (1998). Por outro lado, foram semelhantes àquelas obtidas por Mackinnon et al. (1996) e Lemos et al. (1997).

Com relação à duração da lactação, a média observada no presente estudo foi inferior às observadas por Mackinnon et al. (1996), Lemos et al. (1997) e Freitas et al. (1998).

A média de produção total por lactação encontrada no presente estudo foi igual ou superior àquelas observadas pelos demais autores citados. Por outro lado, a duração de lactação média registrada foi inferior às observadas por estes últimos. Logo, verifica-se que a produção média por dia de lactação registrada no presente estudo foi bastante superior (13,05 kg de leite por dia de lactação) às observadas por Mackinnon et al. (1996) (10,14 kg / dia), Lemos et al. (1997) (10,97 kg / dia) e Freitas et al. (1998) (6,94 kg/dia). Assim, percebe-se que, embora a média de produção de leite por dia de lactação observada tenha sido bastante superior àquelas registradas pelos demais autores citados, em função da baixa média de duração da lactação aqui observada, esta superioridade nem sempre se traduziu em maior produção de leite por lactação. Estas observações evidenciam a importância da característica duração da lactação para produção de leite sob condições tropicais com animais mestiços Europeu-Zebu e também o grande potencial produtivo dos mestiços Holandês x Gir sob tais condições.

Destacam-se, na Tabela 2, os elevados coeficientes de variação, mostrando a grande instabilidade das características estudadas. Além disso, percebe-se que não houve grande diferença entre as médias de PL e de PL305, em função da baixa média de duração das lactações.

Na Tabela 3, verifica-se que as médias das características estudadas se elevam consistentemente à medida que as condições de manejo (regime alimentar) melhoram. Da mesma forma, há uma tendência de diminuição da instabilidade das respostas.

Pode-se perceber claramente que, sob os Regimes Alimentares Extensivo (RAEX) e Semi-Intensivo (RASI), os animais não conseguiram, de uma maneira geral, manter lactações com duração adequada (305 dias). Conseqüentemente, não foram observadas grandes diferenças entre as médias de PL e as de PL305 nestes regimes alimentares. Já, sob o Regime Alimentar Intensivo (RAI), pode-se perceber uma diferença mais expressiva entre as médias de PL e PL305, pois, sob esta condição de manejo, os animais encontraram condições de estender suas lactações acima dos 305 dias de duração. Além disso, observa-se, a partir da tabela 4, que, neste regime alimentar, os animais com proporção de genes da Raça Holandesa maior ou igual a 3/4 são encontrados em maior proporção que nos demais regimes alimentares (57,9% no RAI contra 37,1 e 33,8% nos RASI e RAEX).

A comparação do desempenho dos vários grupos genéticos, para as três características estudadas, pode ser feita a partir dos resultados apresentados na Tabela 4.

Observa-se que, no Regime Alimentar Extensivo (RAEX), não houve diferença significativa entre os grupos genéticos estudados para qualquer das características em estudo (Tabela 4).

No Regime Alimentar Semi-Intensivo (RASI), os animais do grupo ³7/8 tiveram desempenho superior aos demais tanto na PL quanto na PL305, enquanto que os animais do grupo 1/4 tiveram desempenho inferior. Os animais do grupo 3/4 apresentaram maior PL do que aqueles do grupo 1/2, enquanto o grupo 5/8 apresentou PL semelhante aos anteriores. Não houve diferença significativa entre o desempenho dos animais dos grupos genéticos 3/4, 5/8 e 1/2 para a PL305 (Tabela 4).

Já, com relação à duração da lactação no RASI, o grupo genético ³7/8 foi superior aos 3/4, 1/2 e 1/4. O grupo 1/4 foi inferior aos demais, exceto ao 1/2. O grupo 3/4 apresentou maior DL que o 1/2. Já o grupo 5/8 teve DL semelhante aos ³7/8, 3/4 e 1/2 (Tabela 4).

No Regime Alimentar Intensivo (RAI), os animais dos grupos genéticos ³7/8, 3/4 e 1/2 tiveram desempenho semelhante nas três características estudadas. Os animais do grupo 5/8 tiveram o pior desempenho nas três características, apesar de não haver diferença significativa com relação aos 1/2 e 1/4 para a DL e ao 1/4 para as demais. Em função do pequeno número de observações e do elevado erro padrão na estimativa das médias por quadrados mínimos (LSMs) para o grupo genético 1/4, não houve diferença significativa entre este grupo e os 3/4 e 1/2, para a PL e a PL305, apesar das LSMs estimadas para aquele grupo genético serem bastante inferiores às estimadas para estes (Tabela 4).

Observa-se que a condição de meio mais hostil no RAEX impediu que os genótipos mais produtivos expressassem seu potencial, nivelando o desempenho dos diversos grupos genéticos. Outro fator que pode ter contribuído para este nivelamento é o valor elevado do erro-padrão das médias estimadas no RAEX, em função, principalmente, do pequeno número de registros para alguns grupos, notadamente para os ³7/8, 5/8 e 1/4 (Tabela 4).

É interessante perceber que não houve diferença estatística entre os ³7/8 e os 3/4 e 1/2 no RAI, apesar desta diferença ter ocorrido no RASI. Esperar-se-ia que esta diferença ocorresse, uma vez que no RAI foram atingidas médias superiores às do RASI, portanto, os animais com maior proporção de genes da raça Holandesa deveriam expressar com maior facilidade a sua maior habilidade leiteira. No entanto, observa-se que, no RAI, os erros-padrão das médias estimadas são bem superiores àqueles observados no RASI, levando a intervalos de confiança mais elevados e, conseqüentemente, à não significância das diferenças encontradas entre estes grupos genéticos. Um dos fatores que contribuiu para este maior intervalo de confiança das médias estimadas no RAI foi o menor número de observações. Porém, provavelmente, o fator mais importante deve estar na inadequação do modelo matemático (modelo 3) para o RAI onde o coeficiente de determinação (R²) foi bastante inferior (em torno de 0,20) em relação àqueles observados nos demais regimes alimentares (em torno de 0,40). Dentre os fatores que podem ter contribuído para esta inadequação, deve-se destacar o fato de que os dados neste regime alimentar foram provenientes de apenas cinco fazendas, sendo aproximadamente 70% destes oriundos de uma única delas.

Ao se compararem as LSMs para PL de cada grupo genético entre os três regimes alimentares, percebe-se que os animais ³7/8, 3/4 e 1/2 apresentaram aumentos expressivos nos seus desempenhos à medida que melhoraram as condições de manejo, enquanto que para os 5/8 e 1/4 estes aumentos ocorreram de maneira menos intensa. Isso significa que os animais 5/8 e 1/4 parecem demonstrar não ter potencial genético para responder com a mesma intensidade dos demais às melhorias nas condições ambientais, o que precisa de investigações futuras, a partir de maior número de dados.

O baixo desempenho produtivo verificado para os animais 5/8 concorda com os resultados de Teodoro & Verneque (1999). Este desempenho, provavelmente, deu-se em função das perdas de heterose e interações epistáticas favoráveis. Apesar disso, o desenvolvimento de uma raça sintética 5/8 Holandês 3/8 Gir pode ser justificada pela simplicidade de uso por parte dos produtores. No entanto, um rigoroso programa de seleção deve ser implementado com o objetivo de minimizar os efeitos das perdas de heterose e interações epistáticas favoráveis.

Já o baixo desempenho observado para os animais do grupo genético 1/4, justifica-se em função do baixo nível de heterose e da menor proporção de genes da raça Holandesa (especializada para produção de leite) verificada neste grupo de animais.

Ainda com relação à Tabela 4, nota-se uma grande variação no número de registros (n) para cada grupo genético nos três regimes alimentares, sendo os animais dos grupos 1/2 e 3/4 de holandês mais numerosos. Este aspecto justifica-se pelo fato de existirem três formas de se chegar ao animal 5/8 (Girolando): (1) cruzamento entre holandês puro e mestiço 1/4; (2) mestiçagem entre 1/2 e 3/4 e; (3) mestiçagem entre animais 5/8. Logo, sendo o Girolando uma raça em formação, entende-se como normal a existência de maiores proporções de animais de grupos genéticos oriundos dos primeiros cruzamentos, ou seja, 1/2, 3/4 e 1/4. No entanto, há um reduzido número de animais do grupo 1/4, provavelmente em função da preferência dos produtores (que têm como objetivo aumentar a produção de leite) em fazer o retrocruzamento utilizando a raça Holandesa.

A comparação entre os ajustes promovidos pelos modelos de regressão 1 e 2 em relação àquele promovido pelo modelo 3 pode ser feita a partir dos resultados apresentados na Tabela 5.

Observa-se que, no RAEX, a variação adicional obtida pelo modelo classificatório (modelo 3) em relação àquelas obtidas pelos dois modelos genéticos de regressão (n^os 1 e 2) não foi significativa (P>0,05). Portanto, os dois modelos de regressão ajustaram os dados tão bem quanto o modelo classificatório (n^º 3), para as três características estudadas no RAEX (Tabela 5).

Resultados semelhantes foram obtidos no RASI, exceto pelo fato de a variação adicional obtida pelo modelo 3 ter sido significativa em relação ao modelo 2 para PL305 (Tabela 5).

Já no RAI apenas o modelo 1 ajustou convenientemente os dados para a PL, enquanto que nenhum dos dois modelos alcançou tal ajuste para a PL305. Por outro lado, tanto no modelo 1 quanto no 2 foi obtido bom ajuste para a DL. É evidente a maior dificuldade dos modelos de regressão em ajustar os dados adequadamente no RAI, particularmente para as características de produção de leite. Como mencionado anteriormente, o maior nível de imprecisão, com relação aos resultados obtidos para o RAI, pode ter contribuído para esta maior dificuldade de ajuste.

Em duas situações, o modelo 2 não foi adequado para ajustar os dados, enquanto que o modelo 1 mostrou tal adequacidade. Isto mostra que em certas situações o modelo aditivo-dominante pode não ser adequado, e mais, nem sempre o efeito epistático pode ser desprezível.

Ainda na Tabela 5, observa-se que os valores de F obtidos para o modelo 1 são consistentemente inferiores àqueles obtidos para o modelo 2, exceto para a DL nos RAEX e RASI, indicando que o ajuste produzido pelo modelo 1 foi mais aproximado daquele produzido pelo modelo 3. Em função disto e também para obter estimativas do efeito das interações epistáticas aditiva x aditiva, as estimativas de g, d e gg foram obtidas a partir do modelo de regressão 1 (Tabela 6).

O efeito de dominância (d) foi significativo nos três regimes alimentares para a PL e a PL305, enquanto para a DL foi significativo apenas no RAI (Tabela 6).

As diferenças aditivas entre as raças (Holandês x Gir) (g) foram significativas para as três características, exceto no RAEX.

Já os efeitos das interações epistáticas aditiva x aditiva (gg) foram significativos apenas no RAI para a PL e a PL305 e, no RASI para a PL305, portanto, exatamente naquelas condições em que o modelo 2 (efeitos aditivos e de dominância) falhou em ajustar os dados adequadamente (Tabela 6).

As estimativas d e gg, no RAI, foram bastante discrepantes em relação às obtidas nos demais regimes alimentares, remetendo-se ao problema anteriormente citado no que diz respeito aos resultados obtidos neste regime alimentar (pequeno número de dados, sendo estes provenientes de poucas fazendas). Assim, parece mais conveniente concentrar a discussão apenas em relação às estimativas feitas nos demais regimes alimentares.

Comparando-se as estimativas de g e d (Tabela 6) nos regimes alimentares extensivo e semi-intensivo, percebe-se que as estimativas de g foram maiores no RASI do que no RAEX, enquanto as estimativas de d foram maiores no RAEX do que no RASI, para as três características estudadas, conseqüentemente a relação d/g foi maior no RAEX do que no RASI. Estes resultados concordam com aqueles obtidos por Madalena et al. (1990), em que as estimativas do efeito genético direto da raça Holandesa foram maiores nas melhores condições de manejo, enquanto que o contrário foi observado com relação às estimativas de heterose, indicando ser a heterose mais importante sob condições de meio menos favoráveis.

Em relação à duração da lactação, o efeito de dominância não foi significativo nos dois regimes alimentares em questão, enquanto o efeito de diferença genética aditiva entre as raças foi significativo para o RASI, indicando que, nas condições do estudo, a importância da heterose foi limitada para esta característica.

Por outro lado, para as características PL e PL305, o efeito da heterose demonstra ser de grande importância. As estimativas g e d mostram que a PL e a PL305 são favoravelmente influenciadas tanto pelo efeito genético aditivo da raça Holandesa quanto pelo efeito da heterose.

O fato de o efeito de interação epistática aditiva x aditiva (gg) ter sido significativo para a PL305 no RASI mostra que nem sempre o modelo aditivo-dominante (n^º 2) descreverá os dados da maneira mais precisa, havendo, portanto, a necessidade de sempre testar a adequação do mesmo em comparação a outros modelos.

Conclusões

Deve-se sempre investigar a ocorrência de interação genótipo-ambiente, sob risco de se incorrer em graves erros na conduta de programas de cruzamento.

Deve-se sempre testar a adequação do modelo genético de regressão aditivo-dominante, pois nem sempre este 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.

Não houve qualquer benefício em elevar a proporção de genes da raça Holandesa nas condições ambientes mais hostis. Entretanto, para as condições de manejo mais favoráveis, tal elevação deu indícios de ser interessante para o aumento da produção.

Agradecimento

À Associação Brasileira dos Criadores de Girolando, que gentilmente cedeu os dados utilizados do presente estudo e esteve sempre disponível para eventuais esclarecimentos.

Literatura Citada

Raimundo Nonato Braga Lôbo

E-mail: [e-mail]

Raimundo Martins Filho

E-mail: [e-mail]

Arlindo de Alencar Araripe Moura

E-mail: [e-mail]

Recebido em: 19/09/01

Aceito em: 11/07/02

• 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.
• CHAUDHRY, M.Z.; WILCOX, C.J.; SIMERL, N.A. Factors affecting performance of Holstein and Jersey by Sahiwal crossbred dairy cattle in Pakistan. Revista Brasileira de Genética,v.16, n.4, p.949-956, 1993.
• FREITAS, A.F.; WILCOX, C.J.; COSTA, C.N. Breed group effects on milk production of brazilian crossbred dairy cows. Journal of Dairy Science, v.81, p.2306-2311, 1998.
• FREITAS, A.F.; WILCOX, C.J.; ROMAN, R.M. Genetic parameters for milk yield and composition of crossbred dairy cattle in Brazil. Revista Brasileira de Genética,v.18, n.2, p.229-235, 1995.
• 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 da Sociedade Brasileira de Zootecnia, v.26, n.4, p.704-708, 1997.
• 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.
• 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, n.7, p.1872-1886, 1990.
• 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.
• 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.
• SAS Institute Inc. SAS/STAT. User's Guide, release6.11.ed. Cary: 1996.
• TEODORO, R.L.; VERNEQUE, R.S. Sistema de cruzamento como alternativa para o melhoramento de bovinos. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 36., 1999, Porto Alegre. Anais... Porto Alegre: Sociedade Brasileira de Zootecnia, 1999. p.127-141.

Datas de Publicação

• Publicação nesta coleção
  09 Set 2005
• Data do Fascículo
  Set 2002

Histórico

• Recebido
  19 Set 2001
• Aceito
  11 Jul 2002