Acessibilidade / Reportar erro

Heterogeneidade de variâncias na avaliação genética de búfalas no Brasil

Heterogeneity of variances on genetic evaluation of buffaloes in Brazil

Resumos

Registros de produção de leite de 754 búfalas da raça Murrah foram utilizados com o objetivo de avaliar o efeito da heterogeneidade de variâncias na avaliação genética. Os componentes de covariância foram estimados pelo método da máxima verossimilhança restrita utilizando-se quatro modelos bicaracterísticos, considerando, como efeitos fixos, estação de parto e rebanho-ano de parto, e idade da vaca como covariável (efeito linear e quadrático). Os quatro modelos utilizados foram: modelo aditivo; modelo de repetibilidade; modelo aditivo com inclusão interação reprodutor x rebanho-ano; modelo de repetibilidade com inclusão da interação reprodutor x rebanho-ano. Os rebanhos foram classificados em duas classes de desvio-padrão fenotípico para produção de leite e análises bicaracterísticas foram realizadas considerando cada classe de desvio-padrão como característica diferente. Foi conduzida também uma análise unicaracterística desconsiderando as classes de desvio-padrão fenotípico, incluindo o efeito da interação reprodutor x rebanho-ano. As estimativas de componentes de variância genética aditiva foram maiores na classe de alto desvio-padrão, comparadas às de baixo desvio-padrão. A maioria dos animais selecionados nos arquivos sem estratificação foi selecionada para alto desvio-padrão. Apesar do aumento nas variâncias aditivas e do erro nas de classes de alto desvio-padrão, suas herdabilidades foram menores, com exceção do modelo 2, cujo herdabilidade foi maior para a classe de alto desvio-padrão. Quando rebanhos são classificados em alto e baixo desvio-padrão fenotípico e a produção de leite nas diferentes classes é avaliada em modelo multicaracterística, a avaliação genética considera a heterogeneidade de variâncias entre rebanhos.

estratificação dos dados; melhoramento genético de búfalos


Milk yield records of 754 Murrah female buffaloes were used to evaluate the effects of heterogeneity of variance among herds on genetic evaluation. The restricted maximum likelihood method was used to estimate the (co)variance components using four bi-trait models, considering season and herd-year of birth as fixed effects and age of the cow as covariable (linear and quadratic effects). The following models were used: additive; repeatability; additive with sire x herd-year interaction; and repeatability with sire x herd-year interaction. The herds were classified in two classes of phenotipic standard deviation for milk production and bi-traits analyses were carried out considering each class of standard deviation as a different characteristic. A single trait analysis was also carried out, disregarding phenotypic standard deviation classes, including sire x herd-year interaction effect. The estimates of additive genetic variance components were higher in the high standard deviation class than those of low standard deviation. Most of the animals selected from files without stratification was selected for high standard deviation. Despite of the increase in additive variances and the error in high standard deviation classes, their heritability were lower, except for model 2, whose heritability was higher for the class with high standard deviation. When herds are classified into high and low phenotypic standard deviation and milk production in the different classes is evaluated in a model trait, genetic evaluation takes into account the heterogeneity of variances among herds.

buffalo breeding; data stratification


MELHORAMENTO, GENÉTICA E REPRODUÇÃO

Heterogeneidade de variâncias na avaliação genética de búfalas no Brasil

Heterogeneity of variances on genetic evaluation of buffaloes in Brazil

Antonia Kécya França MoitaI,* * Correspondências devem ser enviadas para: kecyamoita@yahoo.com.br ; Paulo Sávio LopesII; Robledo de Almeida TorresII; Ricardo Frederico EuclydesII; Humberto TonhatiIII; Ary Ferreira de FreitasIV

IDoutoranda em Genética e Melhoramento/UFV

IIDepartamento de Zootecnia/UFV

IIIDepartamento de Zootecnia/FCAV/UNESP, Jaboticabal, SP

IVEmbrapa Gado de Leite

RESUMO

Registros de produção de leite de 754 búfalas da raça Murrah foram utilizados com o objetivo de avaliar o efeito da heterogeneidade de variâncias na avaliação genética. Os componentes de covariância foram estimados pelo método da máxima verossimilhança restrita utilizando-se quatro modelos bicaracterísticos, considerando, como efeitos fixos, estação de parto e rebanho-ano de parto, e idade da vaca como covariável (efeito linear e quadrático). Os quatro modelos utilizados foram: modelo aditivo; modelo de repetibilidade; modelo aditivo com inclusão interação reprodutor x rebanho-ano; modelo de repetibilidade com inclusão da interação reprodutor x rebanho-ano. Os rebanhos foram classificados em duas classes de desvio-padrão fenotípico para produção de leite e análises bicaracterísticas foram realizadas considerando cada classe de desvio-padrão como característica diferente. Foi conduzida também uma análise unicaracterística desconsiderando as classes de desvio-padrão fenotípico, incluindo o efeito da interação reprodutor x rebanho-ano. As estimativas de componentes de variância genética aditiva foram maiores na classe de alto desvio-padrão, comparadas às de baixo desvio-padrão. A maioria dos animais selecionados nos arquivos sem estratificação foi selecionada para alto desvio-padrão. Apesar do aumento nas variâncias aditivas e do erro nas de classes de alto desvio-padrão, suas herdabilidades foram menores, com exceção do modelo 2, cujo herdabilidade foi maior para a classe de alto desvio-padrão. Quando rebanhos são classificados em alto e baixo desvio-padrão fenotípico e a produção de leite nas diferentes classes é avaliada em modelo multicaracterística, a avaliação genética considera a heterogeneidade de variâncias entre rebanhos.

Palavras-chave: estratificação dos dados, melhoramento genético de búfalos

ABSTRACT

Milk yield records of 754 Murrah female buffaloes were used to evaluate the effects of heterogeneity of variance among herds on genetic evaluation. The restricted maximum likelihood method was used to estimate the (co)variance components using four bi-trait models, considering season and herd-year of birth as fixed effects and age of the cow as covariable (linear and quadratic effects). The following models were used: additive; repeatability; additive with sire x herd-year interaction; and repeatability with sire x herd-year interaction. The herds were classified in two classes of phenotipic standard deviation for milk production and bi-traits analyses were carried out considering each class of standard deviation as a different characteristic. A single trait analysis was also carried out, disregarding phenotypic standard deviation classes, including sire x herd-year interaction effect. The estimates of additive genetic variance components were higher in the high standard deviation class than those of low standard deviation. Most of the animals selected from files without stratification was selected for high standard deviation. Despite of the increase in additive variances and the error in high standard deviation classes, their heritability were lower, except for model 2, whose heritability was higher for the class with high standard deviation. When herds are classified into high and low phenotypic standard deviation and milk production in the different classes is evaluated in a model trait, genetic evaluation takes into account the heterogeneity of variances among herds.

Key words: buffalo breeding, data stratification

Introdução

De acordo com Tonhati (1997), a quase totalidade dos criatórios de búfalos no Brasil adota regime de criação extensivo, caracterizado pela falta de controle zootécnico, sanitário e nutricional, desconhecendo o desempenho dos búfalos quanto ao ganho em peso, produção e qualidade do leite e características reprodutivas.

Estimativas de componentes de variâncias têm se mostrado heterogêneas de acordo com diferentes níveis de produção e classes de desvio-padrão fenotípico dos rebanhos (Araújo et al., 2001; Torres et al., 1999). Esta propriedade das variâncias proporciona desvios, acarretando riscos de serem selecionadas maiores proporções de animais dos rebanhos de maior variabilidade fenotípica e não necessariamente de maior valor genético, podendo, assim, reduzir o progresso genético esperado (Hill et al., 1983).

Gianola (1986) sugeriu uma análise multicaracterística, na qual os arquivos dos rebanhos com diferentes variâncias poderiam ser considerados como características diferentes. Entretanto, esse método requer estimativas de compo-nentes de variância para cada rebanho.

Uma alternativa para solucionar os problemas associados à heterogeneidade de variâncias seria a estratificação dos dados baseada em algum critério, e posterior análise de características múltiplas (Gianola, 1986), tendo como princípio a proposta apresentada por Falconer (1952), de considerar a expressão de um genótipo em diferentes ambientes, como características distintas.

De acordo com Martins (2002), de forma geral, três situações de ocorrência de heterogeneidade de variância podem ser discriminadas. A primeira situação é aquela em que há apenas heterogeneidade de variância residual, indicando que a fonte de heterogeneidade não afeta as diferenças entre genótipos, porém altera a magnitude dos resíduos. A segunda é aquela em que há apenas heterogeneidade de variância genética, indicando que a fonte de heterogeneidade de variância afeta a magnitude das diferenças entre os genótipos, mas não afeta a magnitude dos resíduos. A terceira é aquela em que há heterogeneidade de variância genética e residual, indicando que a fonte de heterogeneidade afeta tanto a magnitude das diferenças entre genótipos como dos resíduos.

Com este estudo, objetivou-se avaliar o efeito da heterogeneidade de variâncias na avaliação genética de búfalas leiteiras no Brasil.

Material e Métodos

Os dados são provenientes de rebanhos bubalinos, coletados pelo Controle Leiteiro do Departamento de Zootecnia da Universidade Estadual de São Paulo - UNESP, câmpus de Jaboticabal.

Os dados foram editados de forma a eliminar registros de lactações incompletas e encerradas por causas anormais de secagem, lactações inferiores a 64 dias e registros após 305 dias de lactação. A idade das búfalas ao parto variou de 24 a 266 meses.

Os meses de parto das búfalas foram agrupados em duas épocas: época 1, de abril a setembro, que corresponde aos meses de menor incidência de chuvas; e época 2, de outubro a março, em que há maior precipitação pluviométrica, à exceção de um rebanho situado na região Nordeste, onde se considerou o inverso. Posteriormente, os efeitos fixos de rebanho e ano foram agrupados em subclasses, sendo exigidas, no mínimo, quatro observações por subclasse de rebanho-ano.

Feitas as edições necessárias, o arquivo de dados continha 1.774 lactações, oriundas de 754 búfalas da raça Murrah, filhas de 39 reprodutores, as quais pariram no período de 1987 a 2005, distribuídas em treze rebanhos, doze no Estado de São Paulo e um no Estado do Ceará. A duração da lactação foi de 64 a 305 dias e a idade ao parto variou de 24 a 185 meses. A conexidade dos dados foi feita para avaliar reprodutor x rebanho em função das filhas, para isso, foi exigido que cada reprodutor tivesse, no mínimo, duas filhas distribuídas em mais de um rebanho.

As estimativas dos componentes de variância e parâmetros genéticos e os valores genéticos preditos foram obtidos utilizando o programa MTDFREML (Multiple Trait Derivative Free Restricted Maximum Likelihood), descrito por Boldman et al. (1995), utilizando modelo animal. Esse programa utiliza o algoritmo simplex para localizar o mínimo de -2logeL, sendo L a função de verossimilhança; os componentes de covariância que minimizam a função -2logeL são estimativas de verossimilhança. Utilizou-se como critério de convergência a variância dos valores simples, igual a 10-9. Para conferir a convergência repetiu-se a análise e observou-se o valor do -2logeL, de modo a observar se os valores foram iguais até a quarta casa decimal.

Com o intuito de investigar a heterogeneidade de variâncias, utilizou-se o desvio-padrão fenotípico dos rebanhos para produção de leite até 305 dias de lactação, como critério para estratificar os dados em duas classes, com número aproximadamente igual de observações. Os rebanhos foram classificados em duas classes de desvio-padrão fenotípico para produção de leite e efetivaram-se análises bicaracterísticas, considerando-se cada classe de desvio-padrão como característica diferente. Foi conduzida, também, uma análise unicaracterística, na qual foram desconsideradas as classes de desvio-padrão fenotípico, incluindo o efeito da interação reprodutor x rebanho-ano. A classe de baixo desvio-padrão fenotípico corresponde aos dados cujo desvio-padrão da produção de leite do rebanho foi menor ou igual a 600 kg, e a de alto desvio-padrão fenotípico, maior que 600 kg.

O arquivo de pedigree que originou a matriz de numeradores dos coeficientes de parentesco (NRM), utilizada em todas as análises, continha 864 animais diferentes, 1.776 elementos não-zero na NRM e nenhum animal endogâmico, consequentemente, o coeficiente médio de endogamia foi igual a zero.

As análises foram conduzidas utilizando-se quatro modelos bcaracterísticos, nos quais foram considerados como efeitos fixos estação de parto e rebanho-ano de parto e idade da vaca como covariável (efeito linear e quadrático). Os efeitos aleatórios dos quatro modelos são descritos a seguir: modelo bicaracterístico aditivo (efeito aditivo e erro); modelo bicaracterístico de repetibilidade (modelo 1 mais efeito de ambiente permanente); modelo bicaracterístico aditivo com a interação reprodutor x rebanho-ano (modelo 1 incluindo a interação reprodutor x rebanho-ano). Modelo bicaracterístico de repetibilidade com a interação reprodutor x rebanho-ano (modelo 2 incluindo a interação reprodutor x rebanho-ano).

A estimativa da estatística do teste da razão de verossimilhança (LR) foi comparada com o valor obtido por meio da distribuição de qui-quadrado, com 1 grau de liberdade. A estimativa foi obtida pela expressão:

Na qual LRij = estatística do teste da razão de verossimilhança para modelos sequencialmente reduzidos; Li = máximo da verossimilhança para o modelo completo i; e Lj = máximo da verossimilhança para o modelo reduzido j.

A hipótese de nulidade a ser testada implicou que as funções de verossimilhança dos modelos completo e reduzido não diferiam entre si, ou seja,

A regra decisória empregada foi: se LRij > , o teste era significativo e o modelo completo fornecia maior valor da função de verossimilhança em relação ao modelo reduzido; se LRij < , o teste era não-significativo e o modelo completo não fornecia maior valor da função de verossimilhança em relação ao modelo reduzido.

Os valores genéticos aditivos preditos dos reprodutores para produção de leite até 305 dias de lactação, incluindo ou não no modelo o termo da interação reprodutor x rebanho-ano, foram organizados em arquivos, objetivando verificar possíveis alterações no ordenamento das predições dos valores genéticos aditivos dos reprodutores, quando o termo da interação reprodutor x rebanho-ano foi incluído no modelo por meio da obtenção da correlação de Spearman e de Pearson.

Resultados e Discussão

Os rebanhos que foram classificados como de baixo desvio-padrão apresentaram baixa produção, com exceção do rebanho 1011, que possui baixo desvio-padrão e alta produção, assim como os que apresentaram alto desvio-padrão apresentaram alta produção.

Araújo (2000) e Torres et al. (1999) trabalharam com bovinos da raça Holandesa e também observaram que as médias de produção de leite do rebanho aumentaram da classe de baixo desvio-padrão fenotípico para a classe de alto desvio-padrão fenotípico do rebanho.

As inclusões do efeito de ambiente permanente (modelo 2 vs modelo 1), da interação reprodutor x rebanho-ano no modelo de aditivo (modelo 3 vs 1) e da interação reprodutor x rebanho-ano no modelo de repetibilidade (modelo 4 vs modelo 2) não foram significativas.

As estimativas de componentes de variância genética aditiva para produção de leite foram maiores na classe de alto desvio-padrão, comparadas às de baixo desvio-padrão. Entre os modelos, a maior variância aditiva foi verificada no modelo 1, para a classe de alto desvio-padrão. Enquanto as estimativas de componentes de variâncias residuais foram maiores no modelo 2, para a classe de alto desvio-padrão.

Torres et al. (1999), Costa (1999) e Short et al. (1990), em pesquisa com vacas da raça Holandesa, também observaram maiores valores nas estimativas de componentes de variância genética aditiva e residual nas classes de maior desvio-padrão fenotípico. Ao passo que Boldman & Freeman (1990) obtiveram os mesmos resultados trabalhando com médias de produção em vacas.

Houve diminuição na estimativa do componente de variância do efeito de ambiente permanente quando o efeito da interação reprodutor x rebanho-ano foi acrescentado.

A inclusão do efeito de ambiente permanente no modelo diminuiu o valor da estimativa de herdabilidade de 0,50 para 0,13 na classe de baixo desvio-padrão, e de 0,44 para 0,40 nas classes de alto desvio-padrão nos modelos 1 e 2. A inclusão da interação reprodutor x rebanho-ano no modelo aditivo não modificou os valores estimados para herdabilidade nas classes de baixo desvio-padrão e reduziu de 0,44 para 0,43 nas classes de alto desvio-padrão. Entretanto, a inclusão da interação reprodutor x rebanho-ano no modelo de repetibilidade, causou aumento na estimativa de herdabilidade de 0,40 para 0,41 na classe de alto desvio-padrão e aumento de 0,13 para 0,49 na classe de baixo desvio-padrão fenotípico.

Com exceção do modelo 2, as estimativas de herdabilidade nas classes de baixo desvios-padrão apresentaram-se maiores que nas de alto desvio-padrão. Entretanto, Torres et al. (1999), Short et al. (1990) e Dong & Mao (1990) verificaram maiores estimativas de herdabilidade para produção de leite em bovinos da raça Holandesa nas classes de alto desvio-padrão fenotípico. Mesmo resultado foi obtido por Hill et al. (1983), trabalhando com médias de produção.

A menor amplitude entre os valores genéticos preditos para os animais ocorreu na classe de baixo desvio-padrão do modelo 2, e a maior na de alto desvio-padrão do modelo 1 (Tabela 5). Todavia, não houve grandes diferenças entre os quatro modelos. O modelo 4 apresentou os menores valores para média (alto desvio-padrão), ao passo que o modelo 2 apresentou o menor valor de desvio-padrão (baixo desvio-padrão).

As correlações entre as classes de baixo e alto desvio-padrão fenotípico foram de 0,92; 0,99, 0,94 e 0,89; respectivamente para os modelos M1, M2, M3 e M4.

As correlações de Pearson e Spearman obtidas entre os dados do arquivo sem estratificação (análise geral) e os dados das classes de alto desvio-padrão foram maiores em todos os modelos quando comparados com os das classes de baixo desvio-padrão. Isso indica que os animas selecionados nos arquivos sem estratificação provavelmente seriam, na maioria, da classe de alto desvio-padrão (Tabela 6).

As menores correlações de Pearson e Spearman foram observadas no modelo 4, entre o arquivo sem estratificação e a classe de baixo desvio-padrão, e entre as classes de baixo e alto desvio-padrão fenotípico. Resultados diferentes foram obtidos por Costa (1999), que encontrou correlações genéticas próximas à unidade para produção de leite em bovinos da raça Holandesa entre as classes de alto e baixo desvio-padrão.

De acordo com Martins (2002), quando as correlações genéticas entre as classes da fonte de heterogeneidade de variância não forem iguais à unidade, conforme verificado no presente trabalho, uma consequência adicional ocorreria porque a avaliação genética tenderia a indicar como superiores os genótipos cuja capacidade adaptativa fosse mais estável, o que não significa, obrigatoriamente, que tais genótipos sejam superiores em todos os ambientes.

Houve grande alteração no ordenamento dos animais com base nos valores genéticos preditos, quando se utilizaram as classes de baixo desvio-padrão nos quatro modelos estudados. Nos 20 melhores animais, 45% nas classes de baixo desvio-padrão e que 80% foram comuns nas classes alto desvio-padrão (Tabela 7) e variou de 0,45 a 0,90 quando se compararam baixo e alto desvios-padrão fenotípico dos modelos 2 e 3, respectivamente.

Ao comparar os 20 melhores animais do modelo com a interação reprodutor x rebanho-ano, sem estratificação ao modelo 1, observa-se que 60% deles são comuns aos de baixo desvio-padrão e 85% deles são comuns aos de alto desvio-padrão (Tabela 8). Esses resultados confirmam os valores obtidos nas correlações de Pearson e de Spearman (Tabela 6).

Conclusões

Ao selecionar animais sem considerar a heterogeneidade de variância, há tendência em favorecer aqueles pertencentes a rebanhos de alta variância, em detrimento aos de menor variabilidade fenotípica. A estratificação dos rebanhos de búfalos para a característica produção de leite em classes de alto e baixo desvio-padrão fenotípico corrige para a heterogeneidade de variância. Quando realizada a estratificação do rebanho pelo nível de produção, é recomendável na avaliação genética o uso de um modelo mais simples (apenas efeito aditivo). Caso não seja feita a estratificação dos rebanhos, é necessário utilizar um modelo mais completo (efeitos aditivo, permanente e da interação reprodutor x rebanho-ano).

Agradecimentos

Aos criadores de búfalos, pela perseverança na dedicação ao melhoramento da espécie. À Associação Brasileira de Criadores de Búfalos. Ao professor Humberto Tonhati, pela cessão dos dados sobre búfalos. Ao Conselho Nacional de Desenvolvimento Científico e Tecnológico - (CNPq) e à Universidade Federal de Viçosa, pela oportunidade de estudar, aprender e desenvolver esta pesquisa.

Recebido em 23/8/2007 e aprovado em 3/7/2009.

  • ARAÚJO, C.V.; TORRES, R.A.; COSTA, C.N. Et Al. Efeito da interação reprodutor x rebanho na acurácia da avaliação genética de reprodutores da raça holandesa no Brasil. Revista Brasileira de Zootecnia, v.30, n.3, p.1000-1006, 2001 (supl. 1).
  • BOLDMAN, K.G.; KRIESE, L.A.; VAN VLECK, L.D. et al. A manual for use of MTDFREML: a set of programs to obtain estimates of variances and covariances (DRAFT). Lincoln: Department of Agriculture/Agriculture Research Service, 1995. 86p.
  • BOLDMAN, K.G.; FREEMAN, A.E. Adjustment for heterogeneity of variance by herd production level in dairy cow and sire evaluation. Journal of Dairy Science, v.73, n.2, p.503-512, 1990.
  • COSTA, C.N. Interação genótipo x ambiente em gado de leite. In: SIMPÓSIO INTERNACIONAL DE GENÉTICA E MELHORAMENTO ANIMAL, 1., 1999, Viçosa, MG. Anais... Viçosa, MG: UFV, 1999. p.161-191.
  • DONG, M.C.; MAO, I.L. Heterogeneity of (co)variance and heritability in different levels of intraherd milk production variance and of herd average. Journal of Dairy Science, v.73, p.843-851, 1990.
  • FALCONER, D.S. The problem of environment and selection. American Nature, v.86, n.830, p.293-298, 1952.
  • GIANOLA, D. Selection criteria and estimation of parameters when the variance is heterogeneous. Theoretical Applied Genetic, v.72, n.5, p.671-677, 1986.
  • HILL, W.G.; EDWARDS, M.R.; AHMED, M.K. A. et al. Heritability of milk yield and composition at different levels and variability of production. Animal Production, v.36, n.1, p.59-68, 1983.
  • MARTINS, E.N. Avaliação genética e heterogeneidade de variância. In: Reunião Anual da Sociedade Brasileira de Zootecnia, 34., 2002, Recife. Anais... Recife: Sociedade Brasileira de Zootecnia, 2002. p.205-214.
  • SHORT, T.H.; BLAKE, R.W.; QUAAS, R.L. et al. Heterogeneous within-herd variance: 1 - Genetic parameters for first and second lactation milk yield of grade Holstein cows. Journal of Dairy Science, v.73, n.11, p.3312-3320, 1990.
  • TONHATI, H. Melhoramento genético de bubalinos para carne e leite. In: OLIVEIRA, G.J.C; ALMEIDA, A.M.L.; SOUZA FILHO, U.A.S. (Eds.) O búfalo no Brasil Cruz das almas: UFBA – Escola de Agronomia, 1997. p.101-113.
  • TORRES, R.A.; BERGMANN, J.A.G.; COSTA, C.N. et al. Ajustamento de variâncias para a produção de leite entre rebanhos da raça holandesa no Brasil. Revista Brasileira de Zootecnia, v.28, n.2, p.295-303, 1999.
  • *
    Correspondências devem ser enviadas para:
  • Datas de Publicação

    • Publicação nesta coleção
      31 Ago 2010
    • Data do Fascículo
      Jul 2010

    Histórico

    • Aceito
      03 Jul 2009
    • Recebido
      23 Ago 2007
    Sociedade Brasileira de Zootecnia Universidade Federal de Viçosa / Departamento de Zootecnia, 36570-900 Viçosa MG Brazil, Tel.: +55 31 3612-4602, +55 31 3612-4612 - Viçosa - MG - Brazil
    E-mail: rbz@sbz.org.br