Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats

BACKGROUND

Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate.

RESULTS

To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds.

CONCLUSIONS

The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Unraveling inbreeding patterns and selection signals in Alpine Grey cattle

Inbreeding plays a crucial role in livestock breeding, influencing genetic diversity and phenotypic traits. Genomic data have helped address limitations posed by incomplete pedigrees, providing deeper insights into breed genetic diversity. This study assesses inbreeding levels via pedigree and genomic approaches and analyzes old and recent inbreeding using runs of homozygosity (ROH), and selection signals in Alpine Grey cattle. Pedigree data from 165 575 individuals, analyzed with INBUPGF90 software, computed inbreeding coefficients. Genomic-based coefficients derived from PLINK v1.9. or DetectRUNS R package analyses of 1 180 individuals' genotypes. Common single nucleotide polymorphisms within ROH pinpointed genomic regions, aggregating into "ROH islands" indicative of selection pressure. Overlaps with USCS Genome Browser unveiled gene presence. Moderate correlations (0.20-0.54) existed between pedigree and genomic coefficients, with most genomic estimators having higher (>0.8) correlation values. Inbreeding averaged 0.04 in < 8 Mb ROH segments, and 0.03 in > 16 Mb segments; > 90% of ROHs were < 8 Mb, indicating ancient inbreeding prevalence. Recent inbreeding proved less detrimental than in cosmopolitan breeds. Two major ROH islands on chromosomes 6 and 7 harbored genes linked to immune response, disease resistance (PYURF, HERC3), and fertility (EIF4EBP3, SRA1). This study underscores the need for detailed inbreeding analyses to understand genetic characteristics and historical changes in local breeds like Alpine Grey cattle. Genomic insights, especially from ROH, facilitated overcoming pedigree limitations, illuminating breed genetic diversity. Our findings reveal ancient inbreeding's enduring genetic impact and ROH islands potential for selective sweeps, elucidating traits in Alpine Grey cattle.

Selection signatures and landscape genomics analysis to reveal climate adaptation of goat breeds

Goats have achieved global prominence as essential livestock since their initial domestication, primarily owing to their remarkable adaptability to diverse environmental and production systems. Differential selection pressures influenced by climate have led to variations in their physical attributes, leaving genetic imprints within the genomes of goat breeds raised in diverse agroecological settings. In light of this, our study pursued a comprehensive analysis, merging environmental data with single nucleotide polymorphism (SNP) variations, to unearth indications of selection shaped by climate-mediated forces in goats. Through the examination of 43,300 SNPs from 51 indigenous goat breeds adapting to different climatic conditions using four analytical methods: latent factor mixed models (LFMM), F-statistics (Fst), Extended haplotype homozygosity across populations (XPEHH), and spatial analysis method (SAM), A total of 74 genes were revealed to display clear signs of selection, which are believed to be influenced by climatic conditions. Among these genes, 32 were consistently identified by at least two of the applied methods, and three genes (DENND1A, PLCB1, and ITPR2) were confirmed by all four approaches. Moreover, our investigation yielded 148 Gene Ontology (GO) terms based on these 74 genes, underlining pivotal biological pathways crucial for environmental adaptation. These pathways encompass functions like vascular smooth muscle contraction, cellular response to heat, GTPase regulator activity, rhythmic processes, and responses to temperature stimuli. Of significance, GO terms about endocrine regulation and energy metabolic responses, key for local adaptation were also uncovered, including biological processes, such as cell differentiation, regulation of peptide hormone secretion, and lipid metabolism. These findings contribute to our knowledge of the genetic structure of climate-triggered adaptation across the goat genome and have practical implications for marker-assisted breeding in goats.

Combined Metabolomics and Biochemical Analyses of Serum and Milk Revealed Parity-Related Metabolic Differences in Sanhe Dairy Cattle

The production performance of dairy cattle is closely related to their metabolic state. This study aims to provide a comprehensive understanding of the production performance and metabolic features of Sanhe dairy cattle across different parities, with a specific focus on evaluating variations in milk traits and metabolites in both milk and serum. Sanhe dairy cattle from parities 1 to 4 (S1, n = 10; S2, n = 9; S3, n = 10; and S4, n = 10) at mid-lactation were maintained under the same feeding and management conditions. The milk traits, hydrolyzed milk amino acid levels, serum biochemical parameters, and serum free amino acid levels of the Sanhe dairy cattle were determined. Multiparous Sanhe dairy cattle (S2, S3, and S4) had a greater milk protein content, lower milk lactose content, and lower solids-not-fat content than primiparous Sanhe dairy cattle (S1). Moreover, S1 had a higher ratio of essential to total amino acids (EAAs/TAAs) in both the serum and milk. The serum biochemical results showed the lower glucose and total protein levels in S1 cattle were associated with milk quality. Furthermore, ultra-high-resolution high-performance liquid chromatography with tandem MS analysis (UPLC-MS/MS) identified 86 and 105 differential metabolites in the serum and milk, respectively, and these were mainly involved in amino acid, carbohydrate, and lipid metabolism. S1 and S2/S3/S4 had significantly different metabolic patterns in the serum and milk, and more vitamin B-related metabolites were significantly higher identified in S1 than in multiparous cattle. Among 36 shared differential metabolites in the serum and milk, 10 and 7 metabolites were significantly and strongly correlated with differential physiological indices, respectively. The differential metabolites identified were enriched in key metabolic pathways, illustrating the metabolic characteristics of the serum and milk from Sanhe dairy cattle of different parities. L-phenylalanine, dehydroepiandrosterone, and linoleic acid in the milk and N-acetylornithine in the serum could be used as potential marker metabolites to distinguish between Sanhe dairy cattle with parities of 1-4. In addition, a metabolic map of the serum and milk from the three aspects of carbohydrates, amino acids, and lipids was created for the further analysis and exploration of their relationships. These results reveal significant variations in milk traits and metabolites across different parities of Sanhe dairy cattle, highlighting the influence of parity on the metabolic profiles and production performance. Tailored nutritional strategies based on parity-specific metabolic profiles are recommended to optimize milk production and quality in Sanhe cattle.

Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle

The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from -0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075-35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.

Genomic insights into Yak (Bos grunniens) adaptations for nutrient assimilation in high-altitudes

High-altitude environments present formidable challenges for survival and reproduction, with organisms facing limited oxygen availability and scarce nutrient resources. The yak (Bos grunniens), indigenous to the Tibetan Plateau, has notably adapted to these extreme conditions. This study delves into the genomic basis of the yak's adaptation, focusing on the positive selection acting on genes involved in nutrient assimilation pathways. Employing techniques in comparative genomics and molecular evolutionary analyses, we selected genes in the yak that show signs of positive selection associated with nutrient metabolism, absorption, and transport. Our findings reveal specific genetic adaptations related to nutrient metabolism in harsh climatic conditions. Notably, genes involved in energy metabolism, oxygen transport, and thermoregulation exhibited signs of positive selection, suggesting their crucial role in the yak's successful colonization of high-altitude regions. The study also sheds light on the yak's immune system adaptations, emphasizing genes involved in response to various stresses prevalent at elevated altitudes. Insights into the yak's genomic makeup provide valuable information for understanding the broader implications of high-altitude adaptations in mammalian evolution. They may contribute to efforts in enhancing livestock resilience to environmental challenges.

Conservation Mitonuclear Replacement: Facilitated mitochondrial adaptation for a changing world

Most species will not be able to migrate fast enough to cope with climate change, nor evolve quickly enough with current levels of genetic variation. Exacerbating the problem are anthropogenic influences on adaptive potential, including the prevention of gene flow through habitat fragmentation and the erosion of genetic diversity in small, bottlenecked populations. Facilitated adaptation, or assisted evolution, offers a way to augment adaptive genetic variation via artificial selection, induced hybridization, or genetic engineering. One key source of genetic variation, particularly for climatic adaptation, are the core metabolic genes encoded by the mitochondrial genome. These genes influence environmental tolerance to heat, drought, and hypoxia, but must interact intimately and co-evolve with a suite of important nuclear genes. These coadapted mitonuclear genes form some of the important reproductive barriers between species. Mitochondrial genomes can and do introgress between species in an adaptive manner, and they may co-introgress with nuclear genes important for maintaining mitonuclear compatibility. Managers should consider the relevance of mitonuclear genetic variability in conservation decision-making, including as a tool for facilitating adaptation. I propose a novel technique dubbed Conservation Mitonuclear Replacement (CmNR), which entails replacing the core metabolic machinery of a threatened species-the mitochondrial genome and key nuclear loci-with those from a closely related species or a divergent population, which may be better-adapted to climatic changes or carry a lower genetic load. The most feasible route to CmNR is to combine CRISPR-based nuclear genetic editing with mitochondrial replacement and assisted reproductive technologies. This method preserves much of an organism's phenotype and could allow populations to persist in the wild when no other suitable conservation options exist. The technique could be particularly important on mountaintops, where rising temperatures threaten an alarming number of species with almost certain extinction in the next century.

Longitudinal genomic analyses of automatically-recorded vaginal temperature in lactating sows under heat stress conditions based on random regression models

BACKGROUND

Automatic and continuous recording of vaginal temperature (TV) using wearable sensors causes minimal disruptions to animal behavior and can generate data that enable the evaluation of temporal body temperature variation under heat stress (HS) conditions. However, the genetic basis of TV in lactating sows from a longitudinal perspective is still unknown. The objectives of this study were to define statistical models and estimate genetic parameters for TV in lactating sows using random regression models, and identify genomic regions and candidate genes associated with HS indicators derived from automatically-recorded TV.

RESULTS

Heritability estimates for TV ranged from 0.14 to 0.20 over time (throughout the day and measurement period) and from 0.09 to 0.18 along environmental gradients (EG, - 3.5 to 2.2, which correspond to dew point values from 14.87 to 28.19 ˚C). Repeatability estimates of TV over time and along EG ranged from 0.57 to 0.66 and from 0.54 to 0.77, respectively. TV measured from 12h00 to 16h00 had moderately high estimates of heritability (0.20) and repeatability (0.64), indicating that this period might be the most suitable for recording TV for genetic selection purposes. Significant genotype-by-environment interactions (GxE) were observed and the moderately high estimates of genetic correlations between pairs of extreme EG indicate potential re-ranking of selection candidates across EG. Two important genomic regions on chromosomes 10 (59.370-59.998 Mb) and16 (21.548-21.966 Mb) were identified. These regions harbor the genes CDC123, CAMK1d, SEC61A2, and NUDT5 that are associated with immunity, protein transport, and energy metabolism. Across the four time-periods, respectively 12, 13, 16, and 10 associated genomic regions across 14 chromosomes were identified for TV. For the three EG classes, respectively 18, 15, and 14 associated genomic windows were identified for TV, respectively. Each time-period and EG class had uniquely enriched genes with identified specific biological functions, including regulation of the nervous system, metabolism and hormone production.

CONCLUSIONS

TV is a heritable trait with substantial additive genetic variation and represents a promising indicator trait to select pigs for improved heat tolerance. Moderate GxE for TV exist, indicating potential re-ranking of selection candidates across EG. TV is a highly polygenic trait regulated by a complex interplay of physiological, cellular and behavioral mechanisms.

Genome-wide landscape of runs of homozygosity and differentiation across Egyptian goat breeds

Understanding the genomic features of livestock is essential for successful breeding programs and conservation. This information is scarce for local goat breeds in Egypt. In the current study, genomic regions with selection signatures were identified as well as runs of homozygosity (ROH), genomic inbreeding coefficients (FROH) and fixation index (FST) were detected in Egyptian Nubian, Damascus, Barki and Boer goat breeds. A total of 46,268 SNP markers and 337 animals were available for the genomic analyses. On average, 145.44, 42.02, 87.90 and 126.95 ROHs were detected per individual in the autosomal genome of the respective breeds. The mean accumulative ROH lengths ranged from 46.5 Mb in Damascus to 360 Mb in Egyptian Nubian. The short ROH segments (< 2 Mb) were most frequent in all breeds, while the longest ROH segments (> 16 Mb) were exclusively found in the Egyptian Nubian. The highest average FROH was observed in Egyptian Nubian (~ 0.12) followed by Boer (~ 0.11), while the lowest FROH was found in Damascus (~ 0.05) and Barki breed (~ 0.03). The estimated mean FST was 0.14 (Egyptian Nubian and Boer), 0.077 (Egyptian Nubian and Barki), 0.075 (Egyptian Nubian and Damascus), 0.071 (Barki and Boer), 0.064 (Damascus and Boer), and 0.015 (Damascus and Barki), for each pair of breeds. Interestingly, multiple SNPs that accounted for high FST values were observed on chromosome 6 in regions harboring ALPK1 and KCNIP4. Genomic regions overlapping both FST and ROH harbor genes related to immunity (IL4R, PHF23, GABARAP, GPS2, and CD68), reproduction (SPATA2L, TNFSF12, TMEM95, and RNF17), embryonic development (TCF25 and SOX15) and adaptation (MC1R, KDR, and KIT), suggesting potential genetic adaptations to local environmental conditions. Our results contribute to the understanding of the genetic architecture of different goat breeds and may provide valuable information for effective preservation and breeding programs of local goat breeds in Egypt.

Candidate genes associated with heat stress and breeding strategies to relieve its effects in dairy cattle: a deeper insight into the genetic architecture and immune response to heat stress

The need for food products of animal origin is increasing worldwide. Satisfying these needs in a way that has minimal impact on the environment requires cutting-edge technologies and techniques to enhance the genetic quality of cattle. Heat stress (HS), in particular, is affecting dairy cattle with increasing frequency and severity. As future climatic challenges become more evident, identifying dairy cows that are more tolerant to HS will be important for breeding dairy herds that are better adapted to future environmental conditions and for supporting the sustainability of dairy farming. While research into the genetics of HS in the context of the effect of global warming on dairy cattle is gaining momentum, the specific genomic regions involved in heat tolerance are still not well documented. Advances in omics information, QTL mapping, transcriptome profiling and genome-wide association studies (GWAS) have identified genomic regions and variants associated with tolerance to HS. Such studies could provide deeper insights into the genetic basis for response to HS and make an important contribution to future breeding for heat tolerance, which will help to offset the adverse effects of HS in dairy cattle. Overall, there is a great interest in identifying candidate genes and the proportion of genetic variation associated with heat tolerance in dairy cattle, and this area of research is currently very active worldwide. This review provides comprehensive information pertaining to some of the notable recent studies on the genetic architecture of HS in dairy cattle, with particular emphasis on the identified candidate genes associated with heat tolerance in dairy cattle. Since effective breeding programs require optimal knowledge of the impaired immunity and associated health complications caused by HS, the underlying mechanisms by which HS modulates the immune response and renders animals susceptible to various health disorders are explained. In addition, future breeding strategies to relieve HS in dairy cattle and improve their welfare while maintaining milk production are discussed.

Copy Number Variation Regions Differing in Segregation Patterns Span Different Sets of Genes

Copy number variations regions (CNVRs) can be classified either as segregating, when found in both parents, and offspring, or non-segregating. A total of 65 segregating and 31 non-segregating CNVRs identified in at least 10 individuals within a dense pedigree of the Gochu Asturcelta pig breed was subjected to enrichment and functional annotation analyses to ascertain their functional independence and importance. Enrichment analyses allowed us to annotate 1018 and 351 candidate genes within the bounds of the segregating and non-segregating CNVRs, respectively. The information retrieved suggested that the candidate genes spanned by segregating and non-segregating CNVRs were functionally independent. Functional annotation analyses allowed us to identify nine different significantly enriched functional annotation clusters (ACs) in segregating CNVR candidate genes mainly involved in immunity and regulation of the cell cycle. Up to five significantly enriched ACs, mainly involved in reproduction and meat quality, were identified in non-segregating CNVRs. The current analysis fits with previous reports suggesting that segregating CNVRs would explain performance at the population level, whereas non-segregating CNVRs could explain between-individuals differences in performance.

Association of Phenotypic Markers of Heat Tolerance with Australian Genomic Estimated Breeding Values and Dairy Cattle Selection Indices

Dairy cattle predicted by genomic breeding values to be heat tolerant are known to have less milk production decline and lower core body temperature increases in response to elevated temperatures. In a study conducted at the University of Melbourne's Dookie Robotic Dairy Farm during summer, we identified the 20 most heat-susceptible and heat-tolerant cows in a herd of 150 Holstein Friesian lactating cows based on their phenotypic responses (changes in respiration rate, surface body temperature, panting score, and milk production). Hair samples were collected from the tip of the cows' tails following standard genotyping protocols. The results indicated variation in feed saved and HT genomic estimated breeding values (GEBVs) (p ≤ 0.05) across age, indicating a potential for their selection. As expected, the thermotolerant group had higher GEBVs for HT and feed saved but lower values for milk production. In general, younger cows had superior GEBVs for the Balanced Performance Index (BPI) and Australian Selection Index (ASI), whilst older cows were superior in fertility, feed saved (FS), and HT. This study demonstrated highly significant (p ≤ 0.001) negative correlations (-0.28 to -0.74) between HT and GEBVs for current Australian dairy cattle selection indices (BPI, ASI, HWI) and significant (p ≤ 0.05) positive correlations between HT and GEBVs for traits like FS (0.45) and fertility (0.25). Genomic selection for HT will help improve cow efficiency and sustainability of dairy production under hot summer conditions. However, a more extensive study involving more lactating cows across multiple farms is recommended to confirm the associations between the phenotypic predictors of HT and GEBVs.

Heat stress effects on milk yield traits and metabolites and mitigation strategies for dairy cattle breeds reared in tropical and sub-tropical countries

Heat stress is an important problem for dairy industry in many parts of the world owing to its adverse effects on productivity and profitability. Heat stress in dairy cattle is caused by an increase in core body temperature, which affects the fat production in the mammary gland. It reduces milk yield, dry matter intake, and alters the milk composition, such as fat, protein, lactose, and solids-not-fats percentages among others. Understanding the biological mechanisms of climatic adaptation, identifying and exploring signatures of selection, genomic diversity and identification of candidate genes for heat tolerance within indicine and taurine dairy breeds is an important progression toward breeding better dairy cattle adapted to changing climatic conditions of the tropics. Identifying breeds that are heat tolerant and their use in genetic improvement programs is crucial for improving dairy cattle productivity and profitability in the tropics. Genetic improvement for heat tolerance requires availability of genetic parameters, but these genetic parameters are currently missing in many tropical countries. In this article, we reviewed the HS effects on dairy cattle with regard to (1) physiological parameters; (2) milk yield and composition traits; and (3) milk and blood metabolites for dairy cattle reared in tropical countries. In addition, mitigation strategies such as physical modification of environment, nutritional, and genetic development of heat tolerant dairy cattle to prevent the adverse effects of HS on dairy cattle are discussed. In tropical climates, a more and cost-effective strategy to overcome HS effects is to genetically select more adaptable and heat tolerant breeds, use of crossbred animals for milk production, i.e., crosses between indicine breeds such as Gir, white fulani, N'Dama, Sahiwal or Boran to taurine breeds such as Holstein-Friesian, Jersey or Brown Swiss. The results of this review will contribute to policy formulations with regard to strategies for mitigating the effects of HS on dairy cattle in tropical countries.

Signatures of selection in indigenous Chinese cattle genomes reveal adaptive genes and genetic variations to cold climate

Cold climate shapes the genome of animals and drives them to carry sufficient genetic variations to adapt to changes in temperature. However, limited information is available about the genome-wide pattern of adaptations to cold environments in cattle. In the present study, we used 777K SNP genotyping (15 Chinese cattle breeds, 198 individuals) and whole genome resequencing data (54 cattle breeds of the world, 432 individuals) to disentangle divergent selection signatures, especially between the cold-adapted (annual average temperature of habitat, 6.24 °C to 10.3 °C) and heat-adapted (20.2 °C to 24.73 °C) Chinese native cattle breeds. Genomic analyses revealed a set of candidate genes (e.g., UQCR11, DNAJC18, EGR1, and STING1) were functionally associated with thermogenesis and energy metabolism. We also characterized the adaptive loci of cattle exposed to cold temperatures. Our study finds new candidate genes and provides new insights into adaptations to cold climates in cattle.

Invited review: integration of technologies and systems for precision animal agriculture-a case study on precision dairy farming

Precision livestock farming (PLF) offers a strategic solution to enhance the management capacity of large animal groups, while simultaneously improving profitability, efficiency, and minimizing environmental impacts associated with livestock production systems. Additionally, PLF contributes to optimizing the ability to manage and monitor animal welfare while providing solutions to global grand challenges posed by the growing demand for animal products and ensuring global food security. By enabling a return to the "per animal" approach by harnessing technological advancements, PLF enables cost-effective, individualized care for animals through enhanced monitoring and control capabilities within complex farming systems. Meeting the nutritional requirements of a global population exponentially approaching ten billion people will likely require the density of animal proteins for decades to come. The development and application of digital technologies are critical to facilitate the responsible and sustainable intensification of livestock production over the next several decades to maximize the potential benefits of PLF. Real-time continuous monitoring of each animal is expected to enable more precise and accurate tracking and management of health and well-being. Importantly, the digitalization of agriculture is expected to provide collateral benefits of ensuring auditability in value chains while assuaging concerns associated with labor shortages. Despite notable advances in PLF technology adoption, a number of critical concerns currently limit the viability of these state-of-the-art technologies. The potential benefits of PLF for livestock management systems which are enabled by autonomous continuous monitoring and environmental control can be rapidly enhanced through an Internet of Things approach to monitoring and (where appropriate) closed-loop management. In this paper, we analyze the multilayered network of sensors, actuators, communication, networking, and analytics currently used in PLF, focusing on dairy farming as an illustrative example. We explore the current state-of-the-art, identify key shortcomings, and propose potential solutions to bridge the gap between technology and animal agriculture. Additionally, we examine the potential implications of advancements in communication, robotics, and artificial intelligence on the health, security, and welfare of animals.

The profile of HSPA1A gene expression and its association with heat tolerance in crossbred cattle and the tropically adapted dwarf Vechur and Kasaragod

Certain livestock breeds are adapted to hot and humid environments, and these breeds have genetics that could be useful in a changing climate. The expression of several genes has been identified as a useful biomarker for heat stress. In this study, the responses to heat exposure of heat-tolerant Vechur and Kasaragod cattle found in Kerala state in India (also known as dwarf Bos taurus indicus) were compared to crossbred cattle (crosses of Bos t. taurus with Bos t. indicus). At various time points during heat exposure, rectal temperature and the expression of HSPA1A were determined, and the relationship between them was characterized. We characterized HSPA1A mRNA in Vechur cattle and performed molecular clock analysis. The expression of HSPA1A between the lineages and at different temperature humidity index (THI) was significant. There were significant differences between the expression profiles of HSPA1A in Kasaragod and crossbred (p < 0.01) and Vechur and crossbred (p < 0.01) cattle, but no significant difference in expression was observed between Vechur and Kasaragod cattle. The genetic distance between Vechur, B. grunniens, B. t. taurus, and B. t. indicus was 0.0233, 0.0059, and 0.007, respectively. The genetic distance between Vechur and the Indian dwarf breed Malnad Gidda was 0.0081. A molecular clock analysis revealed divergent adaptive evolution of Vechur cattle to B. t. taurus, with adaptations to the high temperatures and humidity that are prevalent in their breeding tract in Kerala, India. These results could also prove useful in selecting heat-tolerant animals using HSPA1A as a marker.

Assessing Genetic Diversity and Searching for Selection Signatures by Comparison between the Indigenous Livni and Duroc Breeds in Local Livestock of the Central Region of Russia

Indigenous pig breeds are mainly associated with the adaptive capacity that is necessary to respond adequately to climate change, food security, and livelihood needs, and natural resources conservation. Livni pigs are an indigenous fat-type breed farmed in a single farm in the Orel region and located in the Central European part of the Russian Federation. To determine the genomic regions and genes that are affected by artificial selection, we conducted the comparative study of two pig breeds with different breeding histories and breeding objectives, i.e., the native fat-type Livni and meat-type Duroc breeds using the Porcine GGP HD BeadChip, which contains ~80,000 SNPs. To check the Livni pigs for possible admixture, the Landrace and the Large White breeds were included into the study of genetic diversity as these breeds participated in the formation of the Livni pigs. We observed the highest level of genetic diversity in Livni pigs compared to commercial breeds (UHE = 0.409 vs. 0.319–0.359, p < 0.001; AR = 1.995 vs. 1.894–1.964, p < 0.001). A slight excess of heterozygotes was found in all of the breeds. We identified 291 candidate genes, which were localized within the regions under putative selection, including 22 and 228 genes, which were specific for Livni and Duroc breeds, respectively, and 41 genes common for both breeds. A detailed analysis of the molecular functions identified the genes, which were related to the formation of meat and fat traits, and adaptation to environmental stress, including extreme temperatures, which were different between breeds. Our research results are useful for conservation and sustainable breeding of Livni breed, which shows a high level of genetic diversity. This makes Livni one of the valuable national pig genetic resources.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The f_dM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

Improving Genomic Selection for Heat Tolerance in Dairy Cattle: Current Opportunities and Future Directions

Heat tolerance is the ability of an animal to maintain production and reproduction levels under hot and humid conditions and is now a trait of economic relevance in dairy systems worldwide because of an escalating warming climate. The Australian dairy population is one of the excellent study models for enhancing our understanding of the biology of heat tolerance because they are predominantly kept outdoors on pastures where they experience direct effects of weather elements (e.g., solar radiation). In this article, we focus on evidence from recent studies in Australia that leveraged large a dataset [∼40,000 animals with phenotypes and 15 million whole-genome sequence variants] to elucidate the genetic basis of thermal stress as a critical part of the strategy to breed cattle adapted to warmer environments. Genotype-by-environment interaction (i.e., G × E) due to temperature and humidity variation is increasing, meaning animals are becoming less adapted (i.e., more sensitive) to changing environments. There are opportunities to reverse this trend and accelerate adaptation to warming climate by 1) selecting robust or heat-resilient animals and 2) including resilience indicators in breeding goals. Candidate causal variants related to the nervous system and metabolic functions are relevant for heat tolerance and, therefore, key for improving this trait. This could include adding these variants in the custom SNP panels used for routine genomic evaluations or as the basis to design specific agonist or antagonist compounds for lowering core body temperature under heat stress conditions. Indeed, it was encouraging to see that adding prioritized functionally relevant variants into the 50k SNP panel (i.e., the industry panel used for genomic evaluation in Australia) increased the prediction accuracy of heat tolerance by up to 10% units. This gain in accuracy is critical because genetic improvement has a linear relationship with prediction accuracy. Overall, while this article used data mainly from Australia, this could benefit other countries that aim to develop breeding values for heat tolerance, considering that the warming climate is becoming a topical issue worldwide.

Genome-Wide Estimates of Runs of Homozygosity, Heterozygosity, and Genetic Load in Two Chinese Indigenous Goat Breeds

Runs of homozygosity (ROH) and heterozygosity (ROHet) are windows into population demographic history and adaptive evolution. Numerous studies have shown that deleterious mutations are enriched in the ROH of humans, pigs, cattle, and chickens. However, the relationship of deleterious variants to ROH and the pattern of ROHet in goats have been largely understudied. Here, 240 Guangfeng and Ganxi goats from Jiangxi Province, China, were genotyped using the Illumina GoatSNP50 BeadChip and genome-wide ROH, ROHet, and genetic load analyses were performed in the context of 32 global goat breeds. The classes with the highest percentage of ROH and ROHet were 0.5–2 Mb and 0.5–1 Mb, respectively. The results of inbreeding coefficients (based on SNP and ROH) and ROHet measurements showed that Guangfeng goats had higher genetic variability than most Chinese goats, while Ganxi goats had a high degree of inbreeding, even exceeding that of commercial goat breeds. Next, the predicted damaging homozygotes were more enriched in long ROHs, especially in Guangfeng goats. Therefore, we suggest that information on damaging alleles should also be incorporated into the design of breeding and conservation programs. A list of genes related to fecundity, growth, and environmental adaptation were identified in the ROH hotspots of two Jiangxi goats. A sense-related ROH hotspot (chromosome 12: 50.55–50.81 Mb) was shared across global goat breeds and may have undergone selection prior to goat domestication. Furthermore, an identical ROHet hotspot (chromosome 1: 132.21–132.54 Mb) containing two genes associated with embryonic development (STAG1 and PCCB) was detected in domestic goat breeds worldwide. Tajima’s D and BetaScan2 statistics indicated that this region may be caused by long-term balancing selection. These findings not only provide guidance for the design of conservation strategies for Jiangxi goat breeds but also enrich our understanding of the adaptive evolution of goats.

Interplay Between Stress and Reproduction: Novel Epigenetic Markers in Response to Shearing Patterns in Australian Merino Sheep (Ovis aries)

In this study, we determined the effect(s) of early shearing on Australian Merino ewes (Ovis aries) and their lambs. To test this research question, we used a suite of field and laboratory methods including GPS collars, wool cortisol, and epigenetic change between ewes and lambs identified using Illumina NovaSeq RRBS. Once shorn ewes (n = 24) were kept on their full fleece throughout the entire gestation period, whereas twice (early) shorn ewes (n = 24) had their wool shorn pre-joining. Top-knot wool sample was taken from ewes during pre-joining, day 50 (mid-gestation), and day 90 (late gestation) for laboratory analysis. Ewes were pregnancy scanned at mid-gestation to determine whether they were early or late parturition (this confirmation is provided by the pregnancy scanner based on fetus size). Top-knot wool sample was also taken from the lambs at weaning for hormone and wool quality testing. Ear tissue was taken from ewes at day 50 (mid-gestation) and from lambs at lamb marking for DNA analysis. Results showed that twice or early shorn ewes grazed 10% higher and maintained stronger body condition than once shorn ewes. Wool cortisol levels were also significantly lower in the early shorn ewes between mid- and late gestation. Lambs bred from twice shorn ewes had on average better visual wool quality parameters in terms of micron, spin finesses, and curvature. For the DNA methylation results, when comparing a group of once sheared with twice sheared ewes, we have discovered one locus (Chr20:50404014) that was significantly differentially methylated [False Discovery Rate (FDR) = 0.005]. This locus is upstream of a protein-coding gene (ENSOARG00000002778.1), which shows similarities to the forkhead box C1 (FOXC1) mRNA using BLAST searches. To further our understanding of the potential interaction between pregnancy status and shearing frequency of the ewes, we performed further differential methylation analysis using a combination of shearing treatment and pregnancy scanning status. The comparisons (1) late pregnancy vs. early pregnancy for ewes with one shearing treatment and (2) late pregnancy vs. early pregnancy for sheep with two shearing treatments were carried out to identify associations between loci and pregnancy duration for sheep with either one or two shearing events. We discovered that 36 gene loci were significantly modulated either between different shearing treatments or late vs. early pregnancy status of ewes. This result suggests that maternal pregnancy and nutritional status during gestation influence DNA methylation. We further investigated DNA methylation in lambs and identified 16 annotated gene loci that showed epigenetic modulation as a result of being born from an early or late stage pregnancy. From the genomics data, we pointed out that ewes go through epigenetic modifications during gestation, and there is a degree of intra-individual variation in the reproductive performance of ewes, which could be due to combination of intrinsic (genetic and physiological) and extrinsic (management and climatic) factors. Collectively, this research provides novel dataset combining physiological, molecular epigenetics, and digital tracking indices that advances our understanding of how Merino ewes respond to shearing frequency, and this information could guide further research on Merino sheep breeding and welfare.

Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines

Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.

Limiting factors for milk production in dairy cows: perspectives from physiology and nutrition

Milk production in dairy cows increases worldwide since many decades. With rising milk yields, however, potential limiting factors are increasingly discussed. Particularly, the availability of glucose and amino acids is crucial to maintain milk production as well as animal health. Limitations arise from feed sources, the rumen and digestive tract, tissue mobilization, intermediary metabolism and transport, and the uptake of circulating nutrients by the lactating mammary gland. The limiting character can change depending on the stage of lactation. Although physiological boundaries are prevalent throughout the gestation-lactation cycle, limitations are aggravated during the early lactation period when high milk production is accompanied by low feed intake and high mobilization of body reserves. The knowledge about physiological constraints may help to improve animal health and make milk production more sustainably. The scope of this review is to address contemporary factors related to production limits in dairy cows from a physiological perspective. Besides acknowledged physiological constraints, selected environmental and management-related factors affecting animal performance and physiology will be discussed. Potential solutions and strategies to overcome or to alleviate these constraints can only be presented briefly. Instead, they are thought to address existing shortcomings and to identify possibilities for optimization. Despite a scientific-based view on physiological limits, we should keep in mind that only healthy animals could use their genetic capacity and produce high amounts of milk.

Genomic Diversity Profiling and Breed-Specific Evolutionary Signatures of Selection in Arunachali Yak

Arunachali yak, the only registered yak breed of India, is crucial for the economic sustainability of pastoralist Monpa community. This study intended to determine the genomic diversity and to identify signatures of selection in the breed. Previously available double digest restriction-site associated DNA (ddRAD) sequencing data of Arunachali yak animals was processed and 99,919 SNPs were considered for further analysis. The genomic diversity profiled based on nucleotide diversity, π (π = 0.041 in 200 bp windows), effective population size, Ne (Ne = 83) and Runs of homozygosity (ROH) (predominance of shorter length ROHs) was found to be optimum. Subsequently, 207 regions were identified to be under selective sweeps through de-correlated composite of multiple signals (DCMS) statistic which combined three individual test statistics viz. π, Tajima’s D and |iHS| in non-overlapping 100 kb windows. Mapping of these regions revealed 611 protein-coding genes including KIT, KITLG, CDH12, FGG, FGA, FGB, PDGFRA, PEAR1, STXBP3, olfactory receptor genes (OR5K3, OR5H6 and OR1E1) and taste receptor genes (TAS2R1, TAS2R3 and TAS2R4). Functional annotation highlighted that biological processes like platelet aggregation and sensory perception were the most overrepresented and the associated regions could be considered as breed-specific signatures of selection in Arunachali yak. These findings point towards evolutionary role of natural selection in environmental adaptation of Arunachali yak population and provide useful insights for pursuing genome-wide association studies in future.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.

Genetic Diversity and Population Structure for Resistance and Susceptibility to Mastitis in Braunvieh Cattle

Mastitis is a disease that causes significant economic losses, since resistance to mastitis is a difficult trait to be improved due to its multifactorial occurrence. Therefore, our objective was to characterize a Mexican Braunvieh cattle population for genetic resistance and susceptibility to mastitis. We used 66 SNP markers for 45 candidate genes in 150 animals. The average heterozygosity was 0.445 ± 0.076, a value higher than those reported for some European breeds. The inbreeding coefficient was slightly negative for resistance to subclinical (−0.058 ± 0.055) and clinical (−0.034 ± 0.076) mastitis, possibly due to low selection for the immunological candidate genes that influence these traits. The genotypic profiles for the candidate loci per K-means group were obtained, as well as the group distribution through the graphics of the principal component analysis. The genotypic profiles showed high genetic diversity among groups. Resistance to clinical mastitis had the lowest presence of the heterozygous genotypes. Although the percentage of highly inbred animals (>50%) is up to 13.3%, there are highly heterozygous groups in terms of the studied traits, a favorable indicator of the presence of genetic diversity. The results of this study constitute evidence of the genetic potential of the Mexican Braunvieh population to improve mastitis-related traits.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

Breeding Strategies for Weather Resilience in Small Ruminants in Atlantic and Mediterranean Climates

Many efforts are being made to cope with negative consequences of climate change (CC) on livestock. Among them, selective breeding of resilient animals to CC is presented as an opportunity to maintain high levels of performance regardless of variation in weather. In the present work, we proposed a set of breeding strategies to improve weather resilience in dairy goats raised in north-western European Atlantic conditions and dairy sheep raised in Mediterranean conditions while improving production efficiency at the same time. Breeding strategies differed in the selection emphasis placed on resilience traits, ranging from 0 to 40% in the index. Simulations were carried out mimicking real breeding programs including: milk yield, length of productive life, age at first kidding and mastitis incidence in dairy goats and milk, fat and protein yields, and fertility for dairy sheep. Considering the particular climatic conditions in the two regions, the predicted future climate scenarios, and genetic correlations among breeding objectives, resilience was defined as stability to weather changes for dairy goats and as the ability to improve performance under heat stress for dairy sheep. A strategy giving a selection weight of 10 and 20% for goat and sheep resilience, respectively, resulted in the best overall genetic response in terms of both, production and resilience ability. Not considering resilience in breeding programs could lead to a major production loss in future climate scenarios, whereas putting too much emphasis on resilience would result in a limited progress in milk production.

Whole-Genome Resequencing Points to Candidate DNA Loci Affecting Body Temperature under Cold Stress in Siberian Cattle Populations

Despite the economic importance of creating cold resilient cattle breeds, our knowledge of the genetic basis of adaptation to cold environments in cattle is still scarce compared to information on other economically important traits. Herein, using whole-genome resequencing of animals showing contrasting phenotypes on temperature maintenance under acute cold stress combined with the existing SNP (single nucleotide polymorphism) functional annotations, we report chromosomal regions and candidate SNPs controlling body temperature in the Siberian cattle populations. The SNP ranking procedure based on regional FST calculations, functional annotations, and the allele frequency difference between cold-tolerant and cold-sensitive groups of animals pointed to multiple candidate genes. Among these, GRIA4, COX17, MAATS1, UPK1B, IFNGR1, DDX23, PPT1, THBS1, CCL5, ATF1, PLA1A, PRKAG1, and NR1I2 were previously related to thermal adaptations in cattle. Other genes, for example KMT2D and SNRPA1, are known to be related to thermogenesis in mice and cold adaptation in common carp, respectively. This work could be useful for cattle breeding strategies in countries with harsh climates, including the Russian Federation.

Investigating the Short-Term Effects of Cold Stress on Metabolite Responses and Metabolic Pathways in Inner-Mongolia Sanhe Cattle

Simple Summary

Cold stress is a major environmental stressor affecting cattle performance in temperate regions, which causes impaired welfare and economic losses to cattle producers. The identification of biological mechanisms associated with cold stress response is paramount for developing effective mitigation strategies, such as genomic selection. In this study, we assessed the short-term effects of hyper-cold stress on metabolite responses and metabolic pathways in the serum of Inner-Mongolia Sanhe cattle. Moreover, 19 differential metabolites were found, mainly involved in amino acid metabolism. A further integration of metabolome results and gene expression highlighted the regulation of metabolic changes and related pathways in severe cold exposure, such as “aminoacyl-tRNA biosynthesis” and “valine, leucine, and isoleucine degradation”. In summary, we presented new insights on the short-term effects of severe cold stress as well as metabolites and metabolic pathways associated with cold stress response in Inner-Mongolia Sanhe cattle.

Abstract

Inner-Mongolia Sanhe cattle are well-adapted to low-temperature conditions, but the metabolic mechanisms underlying their climatic resilience are still unknown. Based on the ¹H Nuclear Magnetic Resonance platform, 41 metabolites were identified and quantified in the serum of 10 heifers under thermal neutrality (5 °C), and subsequent exposure to hyper-cold temperature (−32 °C) for 3 h. Subsequently, 28 metabolites were pre-filtrated, and they provided better performance in multivariate analysis than that of using 41 metabolites. This indicated the need for pre-filtering of the metabolome data in a paired experimental design. In response to the cold exposure challenge, 19 metabolites associated with cold stress response were identified, mainly enriched in “aminoacyl-tRNA biosynthesis” and “valine, leucine, and isoleucine degradation”. A further integration of metabolome and gene expression highlighted the functional roles of the DLD (dihydrolipoamide dehydrogenase), WARS (tryptophanyl-tRNA synthetase), and RARS (arginyl-tRNA synthetase) genes in metabolic pathways of valine and leucine. Furthermore, the essential regulations of SLC30A6 (solute carrier family 30 (zinc transporter), member 6) in metabolic transportation for propionate, acetate, valine, and leucine under severe cold exposure were observed. Our findings presented a comprehensive characterization of the serum metabolome of Inner-Mongolia Sanhe cattle, and contributed to a better understanding of the crucial roles of regulations in metabolites and metabolic pathways during cold stress events in cattle.

New loci and neuronal pathways for resilience to heat stress in cattle

While understanding the genetic basis of heat tolerance is crucial in the context of global warming's effect on humans, livestock, and wildlife, the specific genetic variants and biological features that confer thermotolerance in animals are still not well characterized. We used dairy cows as a model to study heat tolerance because they are lactating, and therefore often prone to thermal stress. The data comprised almost 0.5 million milk records (milk, fat, and proteins) of 29,107 Australian Holsteins, each having around 15 million imputed sequence variants. Dairy animals often reduce their milk production when temperature and humidity rise; thus, the phenotypes used to measure an individual's heat tolerance were defined as the rate of milk production decline (slope traits) with a rising temperature-humidity index. With these slope traits, we performed a genome-wide association study (GWAS) using different approaches, including conditional analyses, to correct for the relationship between heat tolerance and level of milk production. The results revealed multiple novel loci for heat tolerance, including 61 potential functional variants at sites highly conserved across 100 vertebrate species. Moreover, it was interesting that specific candidate variants and genes are related to the neuronal system (ITPR1, ITPR2, and GRIA4) and neuroactive ligand-receptor interaction functions for heat tolerance (NPFFR2, CALCR, and GHR), providing a novel insight that can help to develop genetic and management approaches to combat heat stress.