Unraveling metabolic stress response in dairy cows: genetic control of plasma biomarkers throughout lactation and the transition period.

Breeding animals able to effectively respond to stress could be a long-term, sustainable, and affordable strategy to improve resilience and welfare in livestock systems. In the present study, the concentrations of 29 plasma biomarkers were used as candidate endophenotypes for metabolic stress response in single-SNP, gene- and haplotype-based GWAS using 739 healthy lactating Italian Holstein cows and 88,271 variants. Significant genetic associations were found in all the 3 GWAS approaches for plasma γ-glutamyl transferase concentration on BTA17, for paraoxonase on BTA4, and for alkaline phosphatase and zinc on BTA2. On these chromosomes, single-SNP and gene-based chromosome-wide association studies were performed, confirming GWAS findings. The signals identified for paraoxonase, γ-glutamyl transferase, and alkaline phosphatase were in proximity of the genes coding for them. The heritability of these 4 biomarkers ranged from moderate to high (from 0.39 to 0.54). Plasma biomarkers are known to undergo large changes in concentration during metabolic stress in the transition period, with an inter-individual variability in the rate of change and recovery time. Genetics may account in part for these differences. To assess this, we studied a subset of 139 periparturient cows homozygous at 3 SNPs known to be respectively associated with concentration of plasma ceruloplasmin, paraoxonase and γ-glutamyl transferase. We compared the immune-metabolic profile measured in plasma at -7, +5 and +30 d relative to calving between groups of opposite homozygotes. A significant effect of the genotype was found on paraoxonase and γ-glutamyl transferase plasma concentration at all the 3 time points. No evidence for genotype effect was detected for ceruloplasmin. Understanding the genetic control underlying metabolic stress response may suggest new approaches to foster resilience in dairy cows.