Genetic diversity and thermotolerance in Holstein cows: Pathway analysis and marker development using whole‐genome sequencing

Heat stress causes extensive losses in the dairy sector, due to negative effects on milk production and reproduction. Cows have evolved a series of protective mechanisms, (physiological, biochemical, behavioural) to cope with the thermostressing environments, which have allowed the preservation of productive and reproductive potential of specific animals during summer; these animals are considered thermotolerant and could be used to design programs of selective breeding. These programs, targeting the generations of a population of heat-resistant animals, would increase the frequency of the desired phenotypes, tackling the financial losses on one hand and reducing the carbon footprints of the dairy sector on the other. The development of genomics techniques has enabled genome wide variant calling, to detect SNPs associated with the desired phenotypes. In this study, we used a comparative genomics approach to detect genetic variation associated with thermotolerance and to design molecular markers for characterizing the animals as tolerant/sensitive. A total of 40 cows from each group were split in four sequencing pools and a whole-genome sequencing approach was used. Results and conclusion: Genome-wide genetic variation between groups was characterized and enrichment analysis revealed specific pathways which participate in the adaptive mechanisms of thermotolerance, implicated into systemic and cellular responses, including the immune system functionality, Heat Stress and Unfolded Protein Response. The markers made a promising set of results, as specific SNPs in five genes encoding for Heat Shock Proteins were significantly associated with thermotolerance.