Genotype by environment interaction for Holstein cattle populations using autoregressive and within- and across-country multi-trait reaction norms test-day models

The progenies of international bulls in diverse climatic conditions and management levels may lead to different expressions of their genetic potential resulting in a re-ranking of these bulls. Therefore, evaluate the presence of genotype by environment interaction (G × E) within and across countries is important to guide the decision-making on alternative selection strategies. Thus, a two-step reaction norm (RN) approach was used to investigate the presence of G × E in Portuguese and Brazilian Holstein cattle. In step 1, we performed a within-country genetic evaluation using an autoregressive model to obtain precorrected phenotypes and environmental gradients (herd test-day solutions, HTD levels). In step 2, the precorrected phenotypes were considered as two distinct traits in a bi-trait RN model to estimate variance components across HTD levels, genetic correlation between HTD levels in Portugal and Brazil, and RN of the estimated breeding values. Additionally, the genetic correlation between countries using a bi-trait random regression (RR) sire model was obtained. In step 1, genetic additive variance for milk yield (MY) in Portugal was 14.1% higher than in Brazil. For somatic cell score (SCS), the genetic additive variance in Portugal was 12.7% lower than in Brazil. Although similar heritability estimates for SCS were observed in both countries, MY heritabilities were 0.31 for Portugal and 0.23 for Brazil. Genetic correlations (SD) between both countries obtained using RR sire model were 0.78 (0.051) for MY and 0.75 (0.062) for SCS. In step 2, MY genetic correlations among HTD levels within countries were higher than 0.94 for Portugal and 0.98 for Brazil. Somatic cell score genetic correlations among HTD levels ranged from 0.70 to 0.99 for Portugal and from 0.84 to 0.99 for Brazil. The average (SD) of genetic correlation estimates between Portuguese and Brazilian HTD levels were 0.74 (0.009) for MY and 0.57 (0.060) for SCS. These results suggest the presence of G × E for MY and SCS of Holstein cattle between both countries. Although there was no indication of G × E between Brazilian herd environments, the low genetic correlation for SCS indicates potential re-ranking of bulls between extreme environmental gradient in Portugal. Overall, the results of this study evidence the importance of national and international genetic evaluation systems to assist dairy farmers in the selection of the best genotypes to obtain the expected returns from investments in imported semen and to realize genetic progress in dairy populations under local environmental conditions.