Inclusão da epistasia em modelo de avaliação genética de bovinos de corte compostos

The effect of the inclusion of epistasis on the genetic evaluation model of beef cattle composite

ZOOTECNIA E TECNOLOGIA E INSPEÇÃO DE PRODUTOS DE ORIGEM ANIMAL

Inclusão da epistasia em modelo de avaliação genética de bovinos de corte compostos

The effect of the inclusion of epistasis on the genetic evaluation model of beef cattle composite

R.S. Bueno^I; R.A. Torres^II; J.B.S. Ferraz^I; P.S. Lopes^II; J.P. Eler^I; M. Almeida e Silva^III; R.F. Euclydes^II;E.C. Mattos^I

^IFaculdade de Zootecnia e Engenharia de Alimentos - Universidade de São Paulo. Av. Duque de Caxias Norte, nº 225 - Campus da USP. 13635-900 - Pirassununga, SP

^IIUniversidade Federal de Viçosa - Viçosa, MG

^IIIUniversidade Federal de Minas Gerais - Belo Horizonte, MG

RESUMO

Dados de bovinos compostos foram analisados para avaliar o efeito da epistasia nos modelos de avaliação genética. As características analisadas foram os pesos aos 205 (P205) e 390 dias (P390) e perímetro escrotal aos 390 dias (PE390). As análises foram realizadas pela metodologia de máxima verossimilhança considerando-se dois modelos: o modelo 1 incluiu como covariáveis os efeitos aditivos diretos e maternos, e os não aditivos das heterozigoses para os efeitos diretos e para o materno total, e o modelo 2 considerou também o efeito direto de epistasia. Para comparação dos modelos, foram utilizados o critério de informação de Akaike (AIC) e o critério de informação Bayesiano de Schwartz (BIC), e o teste de razão de verossimilhança. A inclusão da epistasia no modelo de avaliação genética pouco alterou as estimativas de componentes de (co)variâncias genéticas aditivas e, consequentemente, as herdabilidades. O teste de verossimilhança e o critério de Akaike sugeriram que o modelo 2, que inclui a epistasia, apresentou maior aderência aos dados para todas as características analisadas. O critério BIC indicou este modelo como o melhor apenas para P205. Para análise genética dessa população, o modelo que considerou o efeito de epistasia foi o mais adequado.

Palavras-chave: cruzados, efeitos não aditivos, perímetro escrotal, pesos

ABSTRACT

Composite bovine data was analyzed with the objective of evaluating the effect of the epistasis parameter in the models of genetic evaluation. The analyzed characteristics were weight at 205 (W205) and 390 days (P390), and scrotal circumference at 390 days (SC390). The analysis were done by the maximum likelihood method, considering two models: model 1, which included as covariates the direct and maternal additive effects, and non-additive of the heterozygosis for the direct and total maternal, and model 2, which also considered the direct epistasis direct. The Akaike Information Criteria (AIC) and the Bayesiano of Schwartz Information Criteria (BIC) were used for the comparison of the models and the test of ratio of likelihood. The inclusion of the epistasis effects on the model of the genetic evaluation did not alter much the estimation of the genetic additive (co)variances components and, consequently the heritability. However, it was significantly superior by the likelihood ratio test for the studied characteristics. Through the BIC, model 2 was more adequate only for W205. For the genetic analysis of that population the model that considers the epitasis is the more adequate.

Keywords: crossbreed, non-additive genetic effects, scrotal circumference, weights

INTRODUÇÃO

Aumento considerável na produção de bovinos de corte por meio do melhoramento genético tem sido alcançado por seleção, cruzamento e pela combinação de ambos, entretanto os modelos usados em avaliações genéticas e programas de melhoramento de animais cruzados e compostos consideram apenas efeitos aditivos e de dominância. Justificativa para o uso de modelos aditivo-dominante a fim de medir o desempenho esperado em populações de animais cruzados baseia-se na pressuposição de que a heterose é principalmente atribuída aos efeitos de dominância.

Estudos desenvolvidos por Gregory et al. (1991) sugerem que a retenção de heterose é linearmente proporcional à heterozigose. Dickerson (1973) propôs a utilização do termo "perdas por recombinação", que seriam proporcionais aos desvios de linearidade entre a heterose e a heterozigose.

A relação de proporcionalidade linear entre heterose e heterozigose também foi observada por Arthur et al. (1999) e Fries et al. (2000). No entanto, esses autores recomendaram que a adição dos efeitos de epistasia aos programas de avaliação genética pode melhorar as estimativas dos valores genéticos dos dados multirraciais ou de cruzamentos e que os efeitos epistáticos são importantes na constituição do genótipo animal.

Apesar de estarem fortemente associados aos componentes aditivos e de dominância, o que dificulta sua estimação, existem evidências de efeitos epistáticos sobre características de produção em cruzamentos de Bos indicus e Bos taurus. Este fato justifica a sistemática adição e testes de componentes epistáticos aos modelos operacionais, para populações multirraciais ou oriundas de cruzamentos avançados.

O objetivo deste trabalho foi avaliar a inclusão do efeito de epistasia em modelos de avaliação genética para características de peso e perímetro escrotal em uma população de bovinos de corte compostos.

MATERIAL E MÉTODOS

Foram analisados dados de peso ao desmame (P205) e aos 390 dias (P390) e do perímetro escrotal aos 390 dias (CE390) de bovinos compostos Montana Tropical^®, relativos a animais nascidos no período de 1994 a 2004. Os animais foram agrupados pela composição em tipos biológicos, N (Bos indicus, zebuínos), A (Bos taurus, taurinos adaptados às condições tropicais), B (Bos taurus, taurinos de origem britânica) e C (Bos taurus, taurinos de origem continental), segundo suas semelhanças de tipo, função, fisiologia, aspectos de crescimento e reprodução.

Os dados foram analisados, por meio da metodologia da máxima verossimilhança, utilizando-se o procedimento MIXED do software SAS (SAS, 1999). Dois modelos que diferiam quanto aos efeitos genéticos foram empregados: o modelo 1 incluiu como covariáveis os efeitos aditivos diretos e maternos, e os não aditivos das heterozigoses para os efeitos diretos e para o materno total, e o modelo 2, que considerou todos os efeitos do modelo acrescido do efeito de epistasia direto. Os efeitos ambientais incluídos nos modelos foram grupo contemporâneo e classe de idade da mãe ao parto. Os grupos de contemporâneos foram definidos como animais de mesmo sexo, nascidos na mesma safra e fazenda, e criados no mesmo grupo de manejo.

Os efeitos de heterose direta e materna foram obtidos a partir das heterozigoses diretas, calculadas segundo a fórmula, , sendo r a porcentagem de contribuição do tipo biológico baseado no sistema NABC para pai i e mãe j.

As heterozigoses diretas foram consideradas: NxA = porcentagem de cruzamento de raças do tipo biológico N e A existente na composição racial do animal; NxB = porcentagem de cruzamento de raças do tipo biológico N e B existente na composição racial do animal; NxC = porcentagem de cruzamento de raças do tipo biológico N e C existente na composição racial do animal; AxB = porcentagem de cruzamento de raças do tipo biológico A e B existente na composição racial do animal; AxC = porcentagem de cruzamento de raças do tipo biológico A e C existente na composição racial do animal; BxC = porcentagem do cruzamento de raças do tipo biológico B e C existente na composição racial do animal. A heterozigose materna total (MT) foi expressa como a soma das heterozigoses parciais, como obtida acima, porém foi calculado um valor total para a mãe do animal.

Quanto ao efeito de epistasia, foi incluso o coeficiente proposto por Dickerson (1973), segundo Koch (1985), em que é a proporção de cada raça no individuo, e heterozigose do indivíduo.

A comparação dos modelos foi feita pelo critério de informação de Akaike (AIC), ou seja, critério de informação Bayesiano de Schwartz (BIC), e o teste de razão de verossimilhança foi aplicado aos modelos.

A estimativa da estatística do teste da razão de verossimilhança (LR) (Rao, 1973) foi comparada com o valor obtido pela distribuição de qui-quadrado (χ²), com 1 grau de liberdade e nível de significância fixo de 5%. Para AIC e BIC, o valor para comparação foi obtido por: AIC=-2lnL+2p e BIC=-2lnL+pln(N-r), em que: p refere-se ao número de parâmetros do modelo; N, ao número total de informações; e r, ao posto da matriz X, matriz de incidência dos efeitos fixos. Menores valores de AIC e BIC indicam modelos mais adequados.

Os dados foram submetidos ao pré-ajuste para os efeitos não aditivos estimados, conforme recomendações de Mourão (2005). Estimativas de componentes de (co)variâncias e de parâmetros genéticos foram obtidas pelo procedimento de máxima verossimilhança restrita usando um algoritmo livre de derivadas, sob um modelo animal que considerou como efeitos aleatórios, aditivo direto e materno e residual e para P205, também o efeito permanente de meio materno, utilizando-se o programa MTDFREML (Boldman et al., 1995), que foram comparadas com as análises univariadas para as características em estudo realizadas por Mourão (2005).

RESULTADOS E DISCUSSÃO

A estatística descritiva dos dados do P205, P390 e PE390 é apresentada na Tab. 1. Médias, desvio-padrão e coeficientes de variação das características são similares às reportadas na literatura (Eler et al., 2000, Mourão, 2007).

Os valores do logaritmo natural da função de verossimilhança (LnL) para os modelos, em cada característica, os valores da razão de verossimilhança (LR) e dos critérios de informação de Akaike (AIC) e Bayesiano (BIC) são apresentados na Tab. 2.

A inclusão do efeito de epistasia no modelo provocou reduções significativas para P205 (P<0,01), P390 e CE390 (P<0,05) no logaritmo natural da função de verossimilhança, conforme teste da razão de verossimilhança que indica superioridade do modelo 2 em relação ao 1. Demeke et al. (2003b), ao testarem modelos diferentes quanto à quantificação dos efeitos não aditivos, também verificaram que os modelos que consideraram os efeitos epistáticos apresentaram maior aderência aos dados em relação aos modelos de dominância. Entre os modelos epistáticos, o modelo de Dickerson apresentou melhor R² em relação aos demais.

Kahi et al. (2000) compararam quatro modelos genéticos: modelo dominante, modelo que consideraram os efeitos de recombinação (ou epistasia), Dickerson e Kinghorn, e um modelo que considerou os efeitos de dominância e interações aditiva x aditiva, na estimação de efeitos aditivos e não aditivos de uma população de bovinos compostos do Kenya. O modelo de Dickerson foi o mais adequado para analisar a referida população.

Para todas as características em estudo, a covariância residual estimada pelos modelos foi menor no modelo 2 (informação não apresentada).

Segundo os critérios de AIC e BIC, um modelo é mais apropriado quanto menor são essas estatísticas. Assim, pelo AIC, o modelo que proporcionou melhor aderência aos dados foi o que incluiu o efeito de epistasia (modelo 2), para todas as características. No entanto, pelo BIC, o modelo 2 apresentou maior aderência somente para P205. Vale ressaltar que ambos, AIC e BIC, penalizam modelos com maior número de parâmetros. Entretanto, para BIC, esta penalidade é mais rigorosa e favorece modelos mais parcimoniosos (Nunez-Antón e Zimmerman, 2000).

Os componentes de (co)variâncias e parâmetros genéticos estimados quando os dados foram ajustados para os efeitos não aditivos dos modelos 1 (heterozigoses diretas e heterozigose materna total) e 2 (heterozigoses diretas, heterozigose materna total e epistasia direta) estão na Tab. 3.

Ao se comparar os modelos, as estimativas de componentes de (co)variâncias tiveram comportamento similar, e as magnitudes não diferiram.

Para P205 e P390, os componentes de variâncias genéticas, aditiva direta e materna foram menores quando estimadas pelo modelo 2, e a covariância genética direta e materna foi negativa para todos os modelos e características em estudo.

As diferenças observadas para os componentes de (co)variância foram pequenas, constatadas nas estimativas de herdabilidade que diferiram em 0,01 para P205 e P390, e não diferiram para PE390 entre os modelos.

A estimativas de herdabilidade direta obtidas para P205 foram 0,26 e 0,25, e materna 0,23 e 0,22, respectivamente, para os modelos 1 e 2. O componente de variância para o efeito permanente de meio foi 0,07 para ambos os modelos. Para P390, as herdabilidades diretas foram 0,22 e 0,21, e maternas de 0,12 e 0,11, para os modelos 1 e 2, respectivamente.

As estimativas de herdabilidades direta encontradas neste estudo para características de pesos foram maiores que as encontradas por Demeke et al. (2003a) para uma população multirracial na Etiópia, os quais obtiveram estimativas de 0,09 e 0,12, para os pesos a desmama e aos 12 meses. Gregory et al. (1995) relataram maiores herdabilidades, com estimativas de 0,21 e 0,38 para peso aos 200 dias; de 0,31 e 0,43 para peso aos 368 dias, e de 0,24 e 0,40 para peso aos 410 dias, em animais compostos e cruzados, respectivamente.

Estimativas de herdabilidades semelhantes foram obtidas por Mourão (2005), ao trabalhar com a mesma população de bovinos compostos.

Em geral, as estimativas encontradas neste estudo estão dentro da variação apresentada na literatura para as características.

CONCLUSÕES

O uso de modelos mais simples para avaliação de bovinos compostos pode erroneamente identificar a contribuição dos efeitos aditivos e não aditivos como efeito de heterose. Modelos que incluam os efeitos de epistasia são mais adequados para análise de populações multirraciais.

AGRADECIMENTOS

À Agropecuária CFM Ltda., pelo fornecimento dos dados, e ao CNPq, pelo apoio financeiro.

Recebido em 20 de julho de 2010

Aceito em 9 de junho de 2011

E-mail: [e-mail]

• ARTHUR, P.F.; HEARSHAW, H.; STEPHENSON, P.D. Direct and maternal additive and heterosis effect from crossing Bos indicus and Bos taurus cattle: Cow and calf performance in two environments. Livest. Prod. Sci., v.57, p.231-241, 1999.
• 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. 125p.
• DEMEKE, S.; NESER, F.W.C.; SCHOEMAN, S.J. Early growth performance of Bos taurus × Bos indicus cattle crosses in Ethiopia: Evaluation of different crossbreeding models. J. Anim. Breed. Gen, v.120, p.39-50, 2003a.
• DEMEKE, S.; NESER, F.W.C.; SCHOEMAN, S.J. Variance components and genetic parameters for early growth traits in a mixed population of purebred Bos indicus and crossbred cattle. Livest. Prod. Sci., v.84, p.11-21, 2003b.
• DICKERSON, G.E. Inbreeding and heterosis in animals. In: Proceedings of the Animal Breeding and Genetics Symposium. Am. Soc. Anim. Sci., p.54-77, 1973.
• ELER, J.P.; FERRAZ, J.B.S.; GOLDEN, B.L. et al Influência da interação touro x rebanho na estimação da correlação entre efeitos genéticos direto e materno em bovinos da raça Nelore. Rev. Bras. Zootec., v.29, p.1642-1648, 2000.
• FRIES, L.A.; JOHNSTON, D.J.; HEARNSHAW, H. et al Evidence of epistatic effects on weaning weight in crossbreed beef cattle. Asian-Aust. J. Anim. Sci., v.13, suppl. B, p.242, 2000.
• GREGORY, K.E.; LUNSTRA, D.D.; CUNDIFF, L.V. et al Breed effects and heterosis in advanced generations of composite populations for puberty and scrotal traits of beef cattle. J. Anim. Sci., v.69, p.2795-2807, 1991.
• GREGORY, K.E.; CUNDIFF, L.V.; KOCH, R.M. Genetic and phenotypic (co)variances for growth and carcass traits of purebred and composite populations of beef cattle. J. Anim. Scie., v.73, p.1920-1926, 1995.
• KAHI, A.K.; THORPE, W.; NITTER, G. et al Crossbreeding for dairy production in the lowland tropics of Kenya: I. Estimation of individual crossbreeding effects on milk production reproductive traits and on cow live weight. Livest. Prod. Sci., v.63, p.39-54, 2000.
• MOURÃO, G.B.; FERRAZ, J.B.S.; ELER, J.P. et al Genetic parameters for growth traits of a Brazilian Bos taurus x Bos indicus beef composite. Genet. Mol. Res., v.6, p.1190-1200, 2007.
• MOURAO, G.B. Estimação de efeitos genéticos aditivos diretos e maternos e não aditivos e predição do desempenho de pesos, perímetro escrotal e musculosidade em uma população de bovinos de corte compostos (Bos taurus x Bos indicus) 2005. 109f. Tese (Doutorado em Zootecnia) - Faculdade de Zootecnia e Engenharia de Alimentos/ Universidade de São Paulo, Pirassununga, SP, 2005.
• NUNEZ-ANTÓN, V.N.; ZIMMERMAN, D.L. Modelling nonstationary longitudinal data. Biometrics., v.56, p.699-705, 2000.
• RAO, C.R. Linear statistical inference and its aplications 2.ed. New York: John Wiley & Sons. 1973. p.417-420.
• STATISTICAL ANALYSES SYSTEM - SAS. SAS/STAT user'guide, version 8.0 ed. Cary: 1999. v.1, 943p.

Datas de Publicação

• Publicação nesta coleção
  06 Set 2011
• Data do Fascículo
  Ago 2011

Histórico

• Aceito
  09 Jun 2011
• Recebido
  20 Jul 2010