Using genome-wide association analysis to characterize environmental sensitivity of milk traits in dairy cattle

Genotype-by-environment interaction (GxE) has been widely reported in dairy cattle. One way to analyze GxE is to apply reaction norm models. The first derivative of a reaction norm is the environmental sensitivity (ES). In the present study we conducted a large-scale, genome-wide association analysis to identify single-nucleotide polymorphisms (SNPs) that affect general production (GP) and ES of milk traits in the German Holstein population. Sire estimates for GP and for ES were calculated from approximately 13 million daughter records by the use of linear reaction norm models. The daughters were offspring from 2297 sires. Sires were genotyped for 54k SNPs. The environment was defined as the average milk energy yield performance of the herds at the time during which the daughter observations were recorded. The sire estimates were used as observations in a genome-wide association analysis, using 1797 sires. Significant SNPs were confirmed in an independent validation set (500 sires of the same population). To separate GxE scaling and other GxE effects, the observations were log-transformed in some analyses. Results from the reaction norm model revealed GxE effects. Numerous significant SNPs were validated for both GP and ES. Many SNPs that affect GP also affect ES. We showed that ES of milk traits is a typical quantitative trait, genetically controlled by many genes with small effects and few genes with larger effect. A log-transformation of the observation resulted in a reduced number of validated SNPs for ES, pointing to genes that not only caused scaling GxE effects. The results will have implications for breeding for robustness in dairy cattle.