Efeito da interação genótipo × ambiente sobre o peso ao nascimento, aos 205 e aos 550 dias de idade de bovinos da raça Nelore na Região Sul do Brasil

Genotype and environment interaction effect on weights at birth, 205 and 550 days of age of Nellore cattle in the South Region of Brazil

MELHORAMENTO, GENÉTICA E REPRODUÇÃO

Efeito da interação genótipo × ambiente sobre o peso ao nascimento, aos 205 e aos 550 dias de idade de bovinos da raça Nelore na Região Sul do Brasil

Genotype and environment interaction effect on weights at birth, 205 and 550 days of age of Nellore cattle in the South Region of Brazil

Jader Silva Lopes^I; Paulo Roberto Nogara Rorato^II; Tomás Weber^III; Arione Augusti Boligon^IV; Juliana Grigoletto Comin^V; Mariana de Almeida Dornelles^V

^ICurso de Graduação em Zootecnia - Departamento de Zootecnia - UFSM. Bolsista PIBIC-CNPq

^IIDepartamento de Zootecnia da UFSM

^IIICurso de Pós-Graduação em Zootecnia - Departamento de Zootecnia - UFSM. Bolsista CAPES

^IVCurso de Pós-Graduação em Genética e Melhoramento Animal - Departamento de Zootecnia - FCAV/UNESP - Jaboticabal, SP

^VCurso de Graduação em Zootecnia - Departamento de Zootecnia - UFSM

RESUMO

Com o objetivo de avaliar o efeito da interação genótipo × ambiente sobre os pesos ao nascimento (PN), à desmama (P205) e ao sobreano (P550), foram analisados registros de peso de 10.874 animais da raça Nelore, filhos de 425 touros e de 7.629 vacas, criados em 46 fazendas nos três estados do Sul do Brasil no período de 1976 a 2001. Os componentes de (co)variância utilizados para estimar os parâmetros genéticos foram obtidos pelo REML, sob modelo animal, em análises uni e bivariadas. Os modelos consideraram aleatórios os efeitos genéticos aditivos direto (PN, P205 e P550) e materno (PN e P205) e o efeito residual e, como fixos, os efeitos de grupo de contemporâneos (PN, P205 e P550), regime alimentar (P205 e P550), condição de criação (P205 e P550), além da idade da vaca ao parto (PN e P205) como covariável, com os efeitos lineares e quadráticos. As herdabilidades diretas para cada estado variaram de 0,38 ± 0,04 a 0,56 ± 0,08; de 0,22 ± 0,03 a 0,47 ± 0,05 e de 0,35 ± 0,03 a 0,51 ± 0,05, respectivamente, para PN, P205 e P550; as herdabilidades maternas variaram de 0,19 ± 0,03 a 0,31 ± 0,03 e de 0,16 ± 0,08 a 0,48 ± 0,07, respectivamente, para PN e P205. As correlações genéticas estimadas entre os estados foram 0,75 ± 0,03; 0,25 ± 0,05; e 0,64 ± 0,07 para PN; 0,79 ± 0,08; -0,04 ± 0,07; e 0,41 ± 0,07 para P205; e 0,53 ± 0,02; 0,08 ± 0,15; e 0,76 ± 0,09 para P550, respectivamente, para os estados do Rio Grande do Sul e Santa Catarina, Rio Grande do Sul e Paraná e Santa Catarina e Paraná. De modo geral, as correlações de Spearman entre as DEP preditas nos estados foram baixas e indicam alterações na ordem de classificação dos reprodutores, como conseqüência da interação genótipo × ambiente.

Palavras-chave: avaliação genética, bovinos de corte, classificação, correlação genética, herdabilidade

ABSTRACT

This study aimed to evaluate the genotype × environment interaction effect on birth (PN), weaning (P205) and yearling (P550) weights of Nellore cattle raised in 46 farms in three states (Rio Grande do Sul-RS, Santa Catarina-SC and Paraná-PR) of the South Region of Brazil, from 1976 to 2001. Records on 10,874 animals, sired by 425 bulls and 7,629 dams were used to estimate (co)variance components by REML. Univariate and bivariate models included the contemporaneous group, feeding type (P205 and P550), management type (P205 and P550), age of cow at calving (PN and P205) covariate (linear and quadratic components) as fixed effects and direct additive genetic (PN, P205 and P550), maternal (PN and P205) and residual as random effects. Direct heritability estimates of PN, P205 and P550 within each State ranged from 0.38 ± 0.04 to 0.56 ± 0.08; from 0.22 ± 0.03 to 0.47 ± 0.05 and from 0.35 ± 0,03 to 0.51 ± 0,05, respectively. Maternal heritability estimates for the States ranged from 0.19 ± 0.03 to 0.31 ± 0.03 (PN) and from 0.16 ± 0.08 to 0.48 ± 0.07 (P205). Genetic correlations for PN, P205 and P550 between states were 0.75 ± 0.03, 0.25 ± 0.05 and 0.64 ± 0.07; 0.79 ± 0.08, -0.04 ± 0.07 and 0.41 ± 0.07 and 0.53 ± 0.02, 0.08 ± 0.15 and 0.76 ± 0.09 respectively for the RS-SC, RS-PR and SC-PR pairs. Overall, the Spearman correlations between breeding values (DEP) for PN, P205 and P550 predicted in different States were low. Changes in sire ranks between states indicate the presence of genotype× environment interaction in growth traits of Nellore cattle.

Key Words: beef cattle, genetic correlation, genetic evaluation, heritability, sire ranks

Introdução

O Brasil é um país de proporções continentais, com sistemas heterogêneos de exploração de bovinos determinados, em grande parte, pelas diferenças climáticas e econômicas e pela disponibilidade de recursos naturais relacionados à produção animal. Esta diversidade de ambientes proporciona oportunidades diferentes de expressão para o mesmo genótipo, o que dificulta a identificação de indivíduos geneticamente superiores para a reprodução.

Há basicamente duas proposições na realização da seleção: uma sugere que os animais devem ser selecionados nos melhores ambientes para que possam expressar seu potencial genético (Hammond, 1947); e a outra sugere que o conjunto de genes responsável pela expressão de determinada característica pode variar de acordo com o ambiente a que o genótipo é exposto. Desse modo, a mesma característica em ambientes diferentes pode ser considerada características diferentes (Falconer, 1952).

Considerando que a variação fenotípica é parcialmente controlada pela herança genética, uma vez que a expressão dos genes é afetada pelo meio, e que, com a difusão da inseminação artificial, a utilização de touros se expandiu atingindo diversos ambientes, é necessário que a magnitude da influência do ambiente sobre a produção seja esclarecida. Neste contexto, os estudos relacionados ao efeito da interação genótipo × ambiente (IGA) sobre o desempenho dos animais assumem maior importância.

Há na literatura diversas definições para IGA, todavia, todas comparam o desempenho de dois ou mais genótipos em pelo menos dois ambientes diferentes (Haldane, 1946; McBride, 1958; Dumlop, 1962). Considerando a teoria de Haldane (1946) e a classificação (rank) dos genótipos em ambientes diferentes, a significância estatística do efeito da interação e os tipos diferentes de interação relatados por Van Vleck et al. (1962; 1963), Pani (1971) propôs a classificação das interações em quatro tipos. No primeiro, apesar das diferenças entre os ambientes, o comportamento dos genótipos é similar e, neste caso, não há interação; no segundo, há uma pequena inversão na ordem de classificação (rank) dos genótipos nos ambientes diferentes, uma vez que neste caso há interação, porém não-significativa; no terceiro, não há inversão na ordem de classificação (rank) dos genótipos, embora o desempenho desses genótipos apresente grandes diferenças nos ambientes diversos. Neste caso, a interação é considerada significativa; finalmente, no quarto tipo, observa-se importante inversão na ordem de classificação (rank) dos genótipos nos ambientes diferentes e a interação é considerada significativa.

Na atividade produtiva, é importante avaliar a magnitude do efeito desta interação e suas conseqüências econômicas. As avaliações genéticas de bovinos de corte abrangem rebanhos submetidos a grandes variações climáticas e de manejo alimentar e sanitário. Por isso, os rebanhos podem diferir muito na média de produção e nas variâncias fenotípicas e genéticas para as principais características de importância econômica.

Neste estudo, avaliou-se a ocorrência de interação genótipo × ambiente e seus efeitos sobre as características de desempenho de bovinos da raça Nelore criados em diferentes estados do Sul do Brasil, por meio das estimativas dos coeficientes de herdabilidade, das correlações genéticas e dos valores genéticos preditos.

Material e Métodos

Para estudar o efeito da IGA sobre as características peso ao nascimento (PN), peso à desmama ajustado para 205 dias de idade (P205) e peso ao sobreano ajustado para 550 dias de idade (P550), foram utilizados registros de desempenho de 10.874 animais nascidos nos estados do Rio Grande do Sul (1.499), Santa Catarina (2.332) e Paraná (7.043) no período de 1976 e 2001. Os animais eram filhos de 425 touros e de 7.629 vacas e foram distribuídos em 46 fazendas participantes do serviço de Controle de Desenvolvimento Ponderal da Associação Brasileira de Criadores de Zebu (ABCZ).

Foram definidos 669 grupos de contemporâneos (GC), constituídos de animais nascidos no mesmo estado, rebanho e ano, em mesma estação, e pertencentes ao mesmo sexo (Tabela 1). O ano foi dividido em duas estações de nascimento, de acordo com as especificidades de clima e o regime de chuvas na Região Sul (abril a setembro e outubro a março). O sistema de criação foi constituído de dois grupos distintos (animais em aleitamento aos 205 dias e animais desmamados aos 205 dias). O regime alimentar foi dividido em três grupos diferentes: o de animais criados extensivamente (59% dos dados), composto de animais mantidos exclusivamente em pastejo sobre gramíneas recebendo apenas sal proteinado; o de animais em semiconfinamento (30% dos dados), no qual os animais receberam suplementação com volumoso ou concentrado ou tinham acesso a pastagens consorciadas de gramíneas com leguminosas; e de animais em confinamento (11%), no qual os animais recebiam toda a alimentação diária (rações, fenos e/ou silagens) no comedouro.

Foram excluídos do arquivo os registros de touros com menos de cinco filhos, as fazendas com menos de 20 animais, os grupos contemporâneos com menos de cinco representantes, as vacas com idade ao parto superior a 15 e inferior a dois anos e os touros com filhos em apenas um dos estados.

O efeito da interação genótipo × ambiente foi avaliado considerando a mesma característica expressa em ambientes diferentes como características distintas determinadas por grupos diferentes de genes, conforme proposto por Falconer (1952). Primeiramente, assumiu-se que a diferença entre as magnitudes das herdabilidades estimadas nos diferentes estados e que as correlações genéticas de baixa magnitude entre os estados, para uma mesma característica, são indícios de IGA.

O arquivo de trabalho foi dividido em três, cada um com as observações dos filhos de um mesmo grupo de reprodutores em dois dos estados, da seguinte forma: Rio Grande do Sul - Santa Catarina; Rio Grande do Sul - Paraná; e Santa Catarina - Paraná. As correlações genéticas para a mesma característica nos estados diferentes foram estimadas por meio de análises bivariadas. Para estimar as herdabilidades, cada um dos três arquivos foi dividido em dois, de acordo com os estados, e foram realizadas análises univariadas.

Para as análises univariadas, adotou-se um modelo animal que considerou como aleatórios os efeitos genéticos aditivos diretos (para PN, P205 e P550) e maternos (para PN e P205) e o residual; e, como fixos, os efeitos de grupo de contemporâneos (para PN, P205 e P550), regime alimentar (para P205 e P550) e condição de criação (para P205 e P550) e a idade da vaca ao parto como covariável, com os efeitos lineares e quadráticos (para PN e P205). Para as análises univariadas e bivariadas realizadas com o arquivo completo, foi considerado, além dos efeitos descritos, o efeito de estado da federação. O efeito de ambiente permanente da vaca não foi considerado no modelo, uma vez que a média de partos por vaca foi inferior a um parto e meio por vaca (1,45 parto/vaca). O modelo geral de análise utilizado pode ser descrito como:

em que: y = vetor das observações para as características PN, P205 ou P550 em cada estado (análises univariadas) e PN, P205 e P550 da progênie de um mesmo reprodutor em dois estados (análises bivariadas); X = matriz de incidência associada aos efeitos fixos (grupo de contemporâneos ao nascimento, estado, condição de criação e regime alimentar); b = vetor de soluções para os efeitos fixos; Z₁ = matriz de incidência associada aos efeitos genéticos aditivos diretos; Z₂ = matriz de incidência associada aos efeitos genéticos aditivos maternos; a = vetor de efeitos genéticos aditivos diretos do animal; m = vetor de efeitos genéticos aditivos maternos; e = vetor dos resíduos.

Os componentes de (co)variância utilizados para estimar os parâmetros e os valores genéticos foram obtidos pelo Método da Máxima Verossimilhança Restrita livre de derivada utilizando-se o aplicativo MTDFREML (Boldman et al., 2001).Aplicou-se um algoritmo simplex para obtenção de -2 Log da função de verossimilhança. Os valores dos componentes de (co)variância que minimizam a função -2 Log da função de verossimilhança são as estimativas de máxima verossimilhança desses componentes. O critério de convergência considerado foi 10^-6, repetindo-se as análises do quanto necessário para que a diferença entre as estimativas permanecesse menor que 10^-6, sempre utilizando os valores das análises anteriores.

As pressuposições assumidas em relação aos vetores y, a, m e e nas análises bivariadas foram:

em que G_a = matriz de covariâncias genéticas aditivas; G_m = a matriz de covariâncias genéticas maternas; R_e = matriz de covariância residual; A = matriz do numerador de relações genético-aditivas; I = matriz identidade; e = produto direto entre matrizes. Assumiu-se ainda que os vetores a , m e e são não-correlacionados entre si. Para as análises univariadas, as pressuposições assumidas foram:

em que: A = matriz de parentesco, I = matriz identidade, R = matriz de variâncias e covariâncias residuais, = componente de variância genética aditiva direta, = componente de variância genética aditiva materna, = componente de variância residual.

Outra forma de avaliar o efeito da IGA foi por meio da mudança na ordem de classificação dos reprodutores. Os reprodutores foram classificados (rank) com base na média de produção de suas progênies (DEP) nos três arquivos, agrupando os estados dois a dois. A similaridade entre a classificação para os reprodutores nos dois estados foram avaliadas por meio da correlação de Spearman (SAS, 2001) considerando três pressões de seleção distintas, 10, 20 e 30%.

Resultados e Discussão

A maior média para PN ocorreu no Paraná (30,44 ± 2,43 kg) e foi 2% maior que a verificada em Santa Catarina (29,12 ± 2,88 kg). A média obtida no Paraná foi superior à observada por Cubas et al. (2001) para esta mesma raça no mesmo estado (28,50 ± 0,38 kg). A maior média para P205 ocorreu no Rio Grande do Sul (177,66 ± 30,68 kg) e foi 14% superior à observada em Santa Catarina (151,97 ± 33,11 kg). As médias observadas nos três estados foram próximas às encontradas na raça Nelore por Eler et al. (2000) para o Sudeste e Centro-Oeste brasileiro e por Toral et al. (2004) para o Mato Grosso do Sul; todavia, são superiores às relatadas por Cubas et al. (2001) para o PR (141,30 ± 1,47kg) e inferiores às relatadas por Ferreira et al. (2001) para o estado de Minas Gerais (189,38 ± 28,25 e 183,62 ± 32,76), também com a mesma raça. O maior valor para P550 foi observado no Rio Grande do Sul (331,97 ± 25,78 kg) e foi 9% maior que o menor valor observado em Santa Catarina (303,35 ± 24,41 kg). Os menores pesos para as três características foram observados em Santa Catarina, apesar de o PN no Rio Grande do Sul ter sido praticamente o mesmo que o verificado em Santa Catarina. O P205 e o P550 foram, respectivamente, 4 e 2% maiores no Rio Grande do Sul em relação ao Paraná. Os menores coeficientes de variação (CV) para as três características ocorreram no Paraná e os maiores em Santa Catarina. Entretanto, apesar de o coeficiente de variação para PN ter sido 20% menor no Paraná em relação ao Rio Grande do Sul, essa diferença foi de apenas 1 e 2%, respectivamente, para P205 e P550.

Tabela 2

As estimativas de herdabilidade direta obtidas por análises univariadas (Tabela 3) para PN variaram de 0,38 ± 0,04 (completo) a 0,56 ± 0,08 (Rio Grande do Sul) e estão no intervalo (0,16 a 0,66) observado por Toral et al. (2004) em três microrregiões do Mato Grosso do Sul e por Machado et al. (1999) e Eler et al. (2000) para a raça Nelore. Para P205, as estimativas de herdabilidade variaram de 0,22 ± 0,03 (Paraná) a 0,47 ± 0,05 (Rio Grande do Sul) e foram próximas aos valores relatados porEler et al. (2000), Ferreira et al. (2001) eToral et al. (2004), de 0,16 a 0,59; entretanto, são superiores aos descritos por Ribeiro et al. (2001), 0,16, e Silveira et al. (2004), 0,17, para a mesma raça. Para P550, as estimativas de herdabilidades variaram de 0,35 ± 0,03 (Paraná) a 0,51 ± 0,05 (Rio Grande do Sul) e corroboram parcialmente os valores relatados por Toral et al. (2004), de 0,35 a 0,45. No entanto, foram superiores às obtidas por Ferreira et al. (2001), 0,17, e Ribeiro et al. (2001), 0,76, para a raça Nelore. As maiores estimativas de herdabilidade para as três características estudadas foram observadas no Rio Grande do Sul e comprovam efeito menor do ambiente sobre o desempenho nas condições de criação adotadas, no entanto, o menor número de observações coletadas em menor número de propriedades com manejo mais homogêneo neste estado pode ter reduzido a variância de ambiente e superestimado esses valores.

As diferenças entre as magnitudes das estimativas de herdabilidade de um estado para outro foram 7, 12 e 7% para PN; 13, 53 e 22% para P205 e 16, 29 e 11% para P550, respectivamente, para os arquivos Rio Grande do Sul - Santa Catarina, Rio Grande do Sul - Paraná e Santa Catarina - Paraná. A maior diferença foi observada sempre para o arquivo Rio Grande do Sul - Paraná, onde os ambientes são mais contrastantes. Essas altas estimativas indicam que grande parte da variação fenotípica na população estudada está relacionada aos efeitos aditivos dos genes, o que possibilita obter ganho genético considerável por meio da seleção. Entretanto, a diferença na magnitude dessas estimativas sugere que a resposta fenotípica das progênies de um mesmo reprodutor é diferente nos três estados.

A habilidade materna expressa a capacidade da vaca em criar suas progênies, por isso, assume grande importância nos programas de seleção. As estimativas de herdabilidade materna obtidas por análises univariadas (Tabela 3) para PN variaram de 0,19 ± 0,03 (Paraná) a 0,31 ± 0,03 (completo) e são superiores aos valores relatados por Eler et al. (2000), de 0,03 a 0,09. Para P205, as estimativas variaram de 0,16 ± 0,06 (Santa Catarina) a 0,48 ± 0,07 (Rio Grande do Sul) e corroboram, em parte, as obtidas por Ribeiro et al. (2001), de 0,16. Entretanto, foram superiores às obtidas por Eler et al. (2000), Ferreira et al. (2001) e Silveira et al. (2004), que variaram de 0,06 a 0,15. A diferença entre a maior e a menor estimativa para PN entre os estados foi de seis pontos percentuais, o que significa que a diferença de ambiente entre os estados não exerceu grande influência sobre o efeito materno no período pré-natal; todavia, esta diferença para P205 foi de 32 pontos percentuais e comprova que houve maior influência no período pós-natal. As diferenças entre as magnitudes destas estimativas de um estado para outro foram 16, 4 e 24% para PN e de 47, 60 e 11% para P205, respectivamente, para os arquivos RS-SC, RS-PR e SC-PR. Esses valores sugerem que a habilidade materna é afetada pelo ambiente e que as filhas de um mesmo reprodutor tendem a apresentar desempenhos diferentes, em maior magnitude para P205, nos diferentes estados.

As correlações genéticas (Tabela 3) entre os desempenhos das progênies de um mesmo reprodutor nos diferentes estados variaram de 0,25 ± 0,05 (Rio Grande do Sul e Paraná) a 0,75 ± 0,03 (Rio Grande do Sul e Santa Catarina) para PN; de -0,04 ± 0,07 (Rio Grande do Sul e Paraná) a 0,79 ± 0,08 (Rio Grande do Sul e Santa Catarina) para P205 e de 0,08 ± 0,15 (Rio Grande do Sul e Paraná) a 0,76 ± 0,09 (Santa Catarina e Paraná) para P550. As menores correlações genéticas foram observadas entre Rio Grande do Sul e Paraná, estados com condições de ambiente mais contrastantes, especialmente se comparadas às regiões oeste e noroeste do Paraná com o Rio Grande do Sul. Essas estimativas diferem, em magnitude, daquelas relatadas por Toral et al. (2004), em torno de 0,32, para a raça Nelore em diferentes microrregiões do estado de Mato Grosso do Sul e por Fridrich et al.(2003) para a raça Tabapuã entre as regiões Sul e Nordeste (-0,86). Entretanto, nesses trabalhos, as correlações genéticas inferiores a 0,80, de acordo com o proposto por Robertson (1959), evidenciam a ocorrência da IGA. Correlação genética baixa significa que os genes relacionados ao crescimento em um estado não são os mesmos genes responsáveis pela expressão desta característica no outro estado. Esses resultados indicam que a inclusão da região de nascimento no modelo de análise para avaliações genéticas aumenta a acurácia da predição do valor genético e sugerem que avaliações genéticas por estado talvez sejam mais eficientes para identificar animais geneticamente superiores.

As correlações de Spearman entre as classificações (rank) dos reprodutores com base nas DEP preditas nos diferentes estados, considerando três intensidades de seleção distintas, variaram de 0,11 (Rio Grande do Sul e Paraná) a 0,75 (Santa Catarina e Paraná) para PN; de 0,15 (Rio Grande do Sul e Paraná) a 0,41 (Santa Catarina e Paraná) para P205; e de 0,12 (Rio Grande do Sul e Paraná) a 0,51 (Santa Catarina e Paraná) para P550 (Tabela 4). Esses valores indicam que as diferenças de ambiente nos diferentes estados promoveram mudanças na ordem relativa de classificação dos touros avaliados geneticamente. O mesmo foi observado por Toral et al. (2004) para diferentes regiões do Mato Grosso do Sul e por Alencar et al. (2005) para a raça Canchim no estado de São Paulo. Assim, é evidente que um reprodutor mais adequado para um estado não o é, necessariamente, para o outro. Considerando a intensidade de seleção de 10% (Tabela 4), o valor máximo de coincidência entre a classificação de um reprodutor em dois estados foi de apenas 35%, portanto, houve 65% de probabilidade de os filhos desses reprodutores apresentarem desempenhos diferenciados. O valor mínimo para a correlação de ordem ocorreu sempre para o conjunto de dados Rio Grande do Sul e Paraná para as três intensidades de seleção (10, 20 e 30%) consideradas.

O sucesso de um programa de melhoramento genético animal depende da proporção do ganho genético transferido de uma geração para outra. Os resultados deste trabalho sugerem que, quando o processo de mudança na média da população é incrementado por reprodutores e/ou sêmen procedentes de outros criatórios, deve ser dispensada atenção especial ao ambiente de origem deste material genético, pois, quanto maior a semelhança entre os ambientes, melhor a resposta fenotípica do material genético imigrante. Portanto, caso a interação genótipo × ambiente não seja considerada, pode ocorrer a escolha de reprodutores inadequados, o que prejudica o progresso genético na população em seleção.

Conclusões

Houve interação genótipo × ambiente nas populações estudadas, portanto, é necessário considerar esse efeito nas avaliações genéticas para características de desempenho da Região Sul do Brasil. Considerando a importância do tema e o volume relativamente pequeno de dados utilizados neste estudo, recomenda-se que outras pesquisas sejam realizadas com amostras mais representativas.

Literatura Citada

Recebido em 22/8/2006 e aprovado em 9/8/2007

Correspondências devem ser enviadas para [e-mail]

• ALENCAR, M.M.; MASCIOLI, A.S.; FREITAS, A.R. Evidências de interação genótipo × ambiente sobre características de crescimento em bovinos de corte. Revista Brasileira de Zootecnia, v.34, n.2, p.489-495, 2005.
• BOLDMAN, K.H.; KRIESE, L.A.; Van VLECK, L.D. et al. A manual for use for MTDFREML. A set of programs to obtain estimates of variance and covariances (Draft) Lincoln: Departament of Agriculture/Agricultural Research Service, 2001. 120p.
• CUBAS, A.C.; PEROTTO, D.; ABRAHÃO, J.J.P. Desempenho até a desmama de bezerros Nelore e cruzados com Nelore. Revista Brasileira de Zootecnia, v.30, n.3 p.694-701, 2001.
• DUMLOP, A.A. Interaction between heredity and environment in the Australian Merino. I. Strain x location interactions in wool traits. Australian Journal of Agricultural Research, v.13, p.503-531, 1962.
• ELER, J.P.; FERRAZ, J.B.S.; GOLDEN, B.L. et al. Influência da interação touro × rebanho na estimação da correlação entre efeitos genéticos direto e materno em bovinos da raça Nelore. Revista Brasileira de Zootecnia, v.29, n.6, p.1642-1648, 2000.
• FERREIRA, V.C.P.; PENA, V.M.; BERGMANN, J.A.G. et al. Interação genótipo-ambiente em algumas características produtivas de gado de corte no Brasil.Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.53, n.3, p.385-392, 2001.
• FALCONER, D.S. The problem of environment and selection. American Naturalist, v.86, p.293-298, 1952.
• FRIDRICH, A.B.; SILVA, M.A.; FRIDRICH, D. et al. Interação genótipo × ambiente em algumas características de produção de bovinos de corte da raça Tabapuã. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 40., 2003, Santa Maria. Anais... Santa Maria: Sociedade Brasileira de Zootecnia, 2003.
• HALDANE, J.B.S. The interaction of nature and nurture. Annals of Eugenics, v.13, p.197-205, 1946.
• HAMMOND, J. Animal breeding in relation to nutrition and environment conditions. Journal of Animal Science, v.22, n.2, p.195-213, 1947.
• MCBRIDE, C. The environment and animal breeding problems. Animal Breeding Abstract, v.26, n.4, p.349-358, 1958.
• MACHADO, P.F.A.; AQUINO, L.H.; GONÇALVES, T.M. Estimativas de parâmetros genéticos e critérios de seleção em características ponderais de bovinos Nelore. Ciência e Agrotecnia, v.23, n.1, p.197-204, 1999.
• PANI, S.N. Genetic x environment interaction in sire evaluation Missouri: 1971. 138p.
• RIBEIRO, M.N.; PIMENTA FILHO, E.C.; MARTINS, G.A. et al. Herdabilidade para efeitos direto e materno de características de crescimento de bovinos Nelore no estado da Paraíba. Revista Brasileira de Zootecnia, v.30, n.4, p.1224-1227, 2001.
• ROBERTSON, A. The sampling variance of genetic correlation coefficient. Biomethic, v.15, p.469-485, 1959.
• STATISTICAL ANALYSIS SYSTEMS -SAS. User's guide: Stat, Version 8 ed. Cary: 2001. (CD-ROM).
• SILVEIRA, J.C.; McMANUS, C.; MACIOLI, A.S. et al. Fatores ambientais e parâmetros genéticos para características produtivas e reprodutivas em um rebanho Nelore no estado do Mato Grosso do Sul. Revista Brasileira de Zootecnia, v.33, n.6, p.1432-1444, 2004.
• TORAL, F.L.B.; SILVA, L.O.C.; MARTINS, E.N. et al. Interação genótipo × ambiente em características de crescimento de bovinos da raça Nelore no Mato Grosso do Sul. Revista Brasileira de Zootecnia, v.33, n.6, p.1445-1455, 2004.
• Van VLECK, L.D. Genotype and environment in sire evaluation. Journal of Dairy Science, v.46, p.983-989, 1963.
• Van VLECK, L.D.; KING, S.C.; DOLITTLE, D.P. Sources of variations in the Cornell controls at two locations. Poultry Science, v.42, p.1114-1125, 1962.

Datas de Publicação

• Publicação nesta coleção
  21 Jan 2008
• Data do Fascículo
  Jan 2008

Histórico

• Recebido
  22 Ago 2006
• Aceito
  09 Ago 2007