Objetivou-se com o presente estudo foram caracterizar e definir ambientes homogêneos de produção de bovinos de corte compostos no Brasil com relação às variáveis climáticas e geográficas, utilizando técnicas exploratórias multivariadas. Verificar a presença de interação genótipo-ambiente (GxE) nas características peso ao nascimento (PN), peso a desmama (PD), ganho de peso da desmama ao sobreano (GP), perímetro escrotal (PE) e musculosidade (MUS). Pela análise de agrupamento não-hierárquico foram agrupadas as regiões similares com relação às variáveis ambientais. Foram formados seis grupos de fazendas. A inclusão do efeito de interação touro-grupo foi avaliada em análises uni-característica. Comparou-se um modelo com o efeito de interação touro-grupo com outro sem esse efeito. Incluir o efeito de interação touro-GEO no modelo de avaliação genética do PN, PD e PE não resultou melhor ajuste aos dados, no entanto não deve ser descartada a hipótese de se considerar outros tipos de efeitos de GxE. Foram estimados parâmetros genéticos por meio de análises multi-característica, considerando-se a mesma característica como diferente em cada grupo de fazendas. Foi verificada heterogeneidade de variância para todas as características. Os coeficientes de herdabilidade nos grupos de fazendas para PN, PD, GP, PE e MUS variaram de 0,15 a 0,25; 0,16 a 0,25; 0,10 a 0,20; 0,17 a 0,31 e 0,17 a 0,24, respectivamente. As correlações genéticas variaram de 0,19 a 0,90 para PN, -0,02 a 0,92 para PD, 0,31 a 0,93 para GP, 0,64 a 0,89 para PE e de 0,18 a 0,80 para MUS nos grupos fazendas. As diferentes estimativas de herdabilidade obtidas entre grupos de fazendas implicam resposta à seleção diferenciada conforme o ambiente em que os animais são criados e selecionados. Pelas correlações genéticas entre as características nas diversas regiões, constatou-se GxE, indicando que os melhores reprodutores para uma determinada região não são sempre os mesmos para as demais. Um modelo hierárquico de norma de reação sob abordagem Bayesiana também foi utilizado para estimação dos componentes de variância, parâmetros genéticos e verificação da existência de GxE. Os gradientes ambientais baseados nas soluções para o efeito de grupo de contemporâneos para PN, PD, GP e PE foram -6,45 a +4,75 kg, -65 a +65 kg, -72 a +112 kg e -6.5 a +5.5 cm, respectivamente. As estimativas de herdabilidade foram crescentes no gradiente ambiental, PN (0,04 a 0,55), PD (0,39 a 0,47), GP (0,01 a 0,43) e PE (0,21 a 0,23). A correlação entre o nível e a inclinação da norma de reação para PN e GP foi de alta magnitude, indicando que os animais de maior valor genético médio foram os que apresentaram maior resposta à melhoria das condições ambientais, caracterizando o efeito de escala da GxE. Para PD e PE, a correlação entre intercepto e inclinação foi baixa implicando reclassificação dos animais em ambientes diferentes. O modelo hierárquico de normas de reação foi útil para descrever alterações nos componentes de variância decorrentes do ambiente e para descrever a presença de GxE nas características estudadas de bovinos compostos. Existe variação genética com respeito à sensibilidade dos animais, o que possibilita a seleção de genótipos mais plásticos ou mais robustos.

The objectives of this study were to characterize and define homogenous production environments of composite beef cattle in Brazil in terms of climatic and geographic variables using multivariate exploratory techniques; to evaluate the presence of genotype by environment interaction (GxE) for birth weight (BW), weaning weight (WW), postweaning gain (PWG), scrotal circumference (SC) and muscling. Nonhierarchical cluster analysis was used to group farms located in regions with similar environmental variables into clusters. Six clusters of farms were formed. The effect of sire-cluster interaction was tested by single-trait analysis. The inclusion of sire-cluster interaction in the genetic evaluation model may not result in better fit to the data for BW, WW and SC. Genetic parameters were estimated by multiple-trait analysis considering the same trait to be different in each cluster. The heritability coefficient in the clusters for BW, WW, PWG, SC and muscling ranged from 0.15 to 0.25; 0.16 to 0.25; 0.10 to 0.20; 0.17 to 0.31 and 0.17 to 0.24, respectively. The genetic correlations ranged from 0.19 to 0.90 for BW, -0.02 to 0.92 for WW, 0.31 to 0.93 for PWG, 0.64 a 0.89 for SC and 0.18 to 0.80 for muscling in the clusters of farms. The different heritability estimates between groups of farms indicates that the response to selection varies with the environment in which animals are selected. The low genetic correlations between traits in the different regions demonstrated the presence of GxE, indicating that the best sires in a certain region are not the same for the other regions. A reaction norm hierarchical model using Bayesian approach was also used for estimation of variance components, genetic parameters and to verify the existence of GxE. Environmental gradients based in solutions for the effect of contemporary groups for BW, WW, PWG and SC were -6.45 to +4.75 kg, -65 kg to +65, -72 to +112 kg and -6.5 to +5.5 cm, respectively. Heritability estimates were increasing in the environmental gradient, BW (0.04 to 0.55), WW (0.39 to 0.47), PWG (0.01 to 0.43) and SC (0.21 to 0.23). The correlation between the level and slope of reaction norm for BW and PWG was of high magnitude, indicating that animals of higher average breeding value were the ones which presented a best response to environmental improvement, characterizing a scale effect on GxE. For WW and SC, the correlation between intercept and slope was low implying reranking of animals in different environments. The reaction norm hierarchical model has been useful to describe changes in the variance components due to the environment and to describe the presence of GxE traits in composite beef cattle. There is genetic variation with respect to the sensitivity of the animals, which enables the selection of genotypes most plastics or more robust.