Genetic basis of thermotolerance in 2 local dairy sheep populations in the Iberian Peninsula

Sheep milk production in the Northern Mediterranean countries heavily relies on local breeds subject to selection schemes to improve milk production. Climate change may shift the range of thermal loads on the animals and challenge their adaptation to the new thermal gradient. The objective of this study was to characterize the genetic component of thermal tolerance of 2 local breeds from the Iberian Peninsula, Latxa and Manchega, belonging to different genetic types that have evolved under different climatic environments (Oceanic for Latxa and Continentalized Mediterranean for Manchega). A total of 79,243 and 2,388,853 test day monthly records of milk, fat and protein yields from 12,882 and 277,904 ewes of Latxa and Manchega breeds, respectively, along a 12-yr period, were matched with the value of the average temperature-humidity index (THI) on the day of milk recording of the closest weather station to the flocks. These data were used to fit individual reaction norms (Legendre polynomials) describing changes in yields along the THI gradient. Genetic values for thermal tolerance were obtained from the slopes of those reaction norms under cold or heat stress and variances and covariances between yield and thermal tolerance were derived from the (co) variance matrices of the polynomial random regression coefficients. Results showed differing patterns in the 2 breeds. The Latxa breed showed clearer signs of genotype by environment interaction than did Manchega. Estimated correlations between yields under extremes of cold and heat were always above 0.8 for Manchega and around 0.4 for Latxa for all traits. Estimates of correlations between comfort and thermal stress were again over 0.8 for Manchega and lower for the higher (hot end) values (down to 0.6) than for the lower (cold end) values (down to 0.84) of the THI gradient for Latxa, indicating greater thermal stress expected from high than from low temperatures in this breed. Substantial variability in thermal tolerance under the more extreme THI values was found in both breeds. Estimated genetic correlations between yield and thermal tolerance were close to 0 and 0.2 for Latxa and Manchega in the cold and down to -0.6 and -0.3 in the heat end, although small variation was observed for fat in Latxa. Estimated realized trends from the average estimated breeding value, by year of birth, showed a positive response for yields in both breeds [around 0.1 standard deviation (SD) unit], but a detrimental correlation for thermal tolerance (down to -0.03 SD units for heat tolerance in Latxa). These results can be used to design optimal selection strategies for sustainable improvement of productivity under a rising temperature scenario associated with climate change.

Automated summary

This study aimed to characterize the genetic component of thermal tolerance in two local breeds, Latxa and Manchega, from the Iberian Peninsula. Results showed different patterns in the genetic interactions between the two breeds, with Manchega exhibiting clearer signs of genotype by environment interaction than Latxa. The findings indicate that there is a greater thermal stress expected from high temperatures in the Latxa breed compared to low temperatures, and both breeds exhibit substantial variability in thermal tolerance.