Genomic analyses and biological validation of candidate genes for rectal temperature as an indicator of heat stress in Holstein cattle

Genetic background of hematological parameters in Holstein cattle based on genome-wide association and RNA sequencing analyses

Hematological parameters refer to the assessment of changes in the number and distribution of blood cells, including leukocytes (LES), erythrocytes (ERS), and platelets (PLS), which are essential for the early diagnosis of hematological system disorders and other systemic diseases in livestock. In this context, the primary objectives of this study were to investigate the genomic background of 19 hematological parameters in Holstein cattle, focusing on LES, ERS, and PLS blood components. Genetic and phenotypic (co)variances of hematological parameters were calculated based on the average information restricted maximum likelihood method and 1,610 genotyped individuals and 5,499 hematological parameter records from 4,543 cows. Furthermore, we assessed the genetic relationship between these hematological parameters and other economically important traits in dairy cattle breeding programs. We also carried out genome-wide association studies and candidate gene analyses. Blood samples from 21 primiparous cows were used to identify candidate genes further through RNA sequencing (RNA-seq) analyses. Hematological parameters generally exhibited low-to-moderate heritabilities ranging from 0.01 to 0.29, with genetic correlations between them ranging from -0.88 ± 0.09 (between mononuclear cell ratio and lymphocyte cell ratio) to 0.99 ± 0.01 (between white blood cell count and granulocyte cell count). Furthermore, low-to-moderate approximate genetic correlations between hematological parameters with one longevity, 4 fertility, and 5 health traits were observed. One hundred ninety-nine significant SNP located primarily on the Bos taurus autosomes (BTA) BTA4, BTA6, and BTA8 were associated with 16 hematological parameters. Based on the RNA-seq analyses, 6,687 genes were significantly downregulated and 4,119 genes were upregulated when comparing 2 groups of cows with high and low phenotypic values. By integrating genome-wide association studies (GWAS), RNA-seq, and previously published results, the main candidate genes associated with hematological parameters in Holstein cattle were ACRBP, ADAMTS3, CANT1, CCM2L, CNN3, CPLANE1, GPAT3, GRIP2, PLAGL2, RTL6, SOX4, WDFY3, and ZNF614. Hematological parameters are heritable and moderately to highly genetically correlated among themselves. The large number of candidate genes identified based on GWAS and RNA-seq indicate the polygenic nature and complex genetic determinism of hematological parameters in Holstein cattle.

Genetic parameters for novel climatic resilience indicators derived from automatically-recorded vaginal temperature in lactating sows under heat stress conditions

BACKGROUND

Longitudinal records of automatically-recorded vaginal temperature (TV) could be a key source of data for deriving novel indicators of climatic resilience (CR) for breeding more resilient pigs, especially during lactation when sows are at an increased risk of suffering from heat stress (HS). Therefore, we derived 15 CR indicators based on the variability in TV in lactating sows and estimated their genetic parameters. We also investigated their genetic relationship with sows' key reproductive traits.

RESULTS

The heritability estimates of the CR traits ranged from 0.000 ± 0.000 for slope for decreased rate of TV (SlopeDe) to 0.291 ± 0.047 for sum of TV values below the HS threshold (HSUB). Moderate to high genetic correlations (from 0.508 ± 0.056 to 0.998 ± 0.137) and Spearman rank correlations (from 0.431 to 1.000) between genomic estimated breeding values (GEBV) were observed for five CR indicators, i.e. HS duration (HSD), the normalized median multiplied by normalized variance (Nor_medvar), the highest TV value of each measurement day for each individual (MaxTv), and the sum of the TV values above (HSUA) and below (HSUB) the HS threshold. These five CR indicators were lowly to moderately genetically correlated with shoulder skin surface temperature (from 0.139 ± 0.008 to 0.478 ± 0.048) and respiration rate (from 0.079 ± 0.011 to 0.502 ± 0.098). The genetic correlations between these five selected CR indicators and sow reproductive performance traits ranged from - 0.733 to - 0.175 for total number of piglets born alive, from - 0.733 to - 0.175 for total number of piglets born, and from - 0.434 to - 0.169 for number of pigs weaned. The individuals with the highest GEBV (most climate-sensitive) had higher mean skin surface temperature, respiration rate (RR), panting score (PS), and hair density, but had lower mean body condition scores compared to those with the lowest GEBV (most climate-resilient).

CONCLUSIONS

Most of the CR indicators evaluated are heritable with substantial additive genetic variance. Five of them, i.e. HSD, MaxTv, HSUA, HSUB, and Nor_medvar share similar underlying genetic mechanisms. In addition, individuals with higher CR indicators are more likely to exhibit better HS-related physiological responses, higher body condition scores, and improved reproductive performance under hot conditions. These findings highlight the potential benefits of genetically selecting more heat-tolerant individuals based on CR indicators.

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.

Single-step genome-wide association analyses for selected infrared-predicted cheese-making traits in Walloon Holstein cows

This study aimed to perform genome-wide association study to identify genomic regions associated with milk production and cheese-making properties (CMP) in Walloon Holstein cows. The studied traits were milk yield, fat percentage, protein percentage, casein percentage (CNP), calcium content, somatic cell score (SCS), coagulation time, curd firmness after 30 min from rennet addition, and titratable acidity. The used data have been collected from 2014 to 2020 on 78,073 first-parity (485,218 test-day records), 48,766 second-parity (284,942 test-day records), and 21,948 third-parity (105,112 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA) of 6,617 animals (1,712 males), were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of ∼216 KB) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for positional candidate genes. Heritability estimates for the studied traits ranged from 0.10 (SCS) to 0.53 (CNP), 0.10 (SCS) to 0.50 (CNP), and 0.12 (SCS) to 0.49 (CNP) in the first, second, and third parity, respectively. Genome-wide association analyses identified 6 genomic regions (BTA1, BTA14 [4 regions], and BTA20) associated with the considered traits. Genes including the SLC37A1 (BTA1), SHARPIN, MROH1, DGAT1, FAM83H, TIGD5, MROH6, NAPRT, ADGRB1, GML, LYPD2, JRK (BTA14), and TRIO (BTA20) were identified as positional candidate genes for the studied CMP. The findings of this study help to unravel the genomic background of a cow's ability for cheese production and can be used for the future implementation and use of genomic evaluation to improve the cheese-making traits in Walloon Holstein cows.

Key parameters of physiological responses to acute heat stress in two commercial layers determined by Fisher discriminant analyses

Given the escalating global warming and the intense nature of modern poultry production, layers are becoming increasingly susceptible to heat stress. This stress disrupts the physiological processes of layers, which leads to reduced productivity and welfare. To address this issue, it is crucial to first evaluate the stress response systematically. However, such evaluations are still lacking in this field. The objective of this study was to accurately monitor the impact of thermal stress and identify common and key indicators that would support decision-making to maintain layer welfare and productivity under stress. We constructed two heat stress models to reflect moderate (32 °C) to severe (36 °C) stress effects and obtained a comprehensive profile of blood physiological parameters associated with the layers' responses to heat stress. We found that genetic differences had limited influence on their physiological responses to heat stress after 32 °C heat challenges. Using 8 selected and significantly changed parameters, layers' physiological status under heat stress could be accurately determined (judgmental accuracy of 98%). As ambient temperature increased to 36 °C, birds suffered more severe challenges that parameters changed in larger percentages. Additionally, breed variations of the physiological responses became apparent, a Fisher discriminant function based on 5 selected parameters could distinguish heat stress effects at 32 °C or 36 °C with 80% accuracy. The results obtained from this study provide two discriminant models for assessing heat stress and shed lights on developing effective and widely applicable heat stress mitigation strategies targeting these indicators.

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.

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.

Invited review: From heat stress to disease-Immune response and candidate genes involved in cattle thermotolerance

Heat stress implies unfavorable effects on primary and functional traits in dairy cattle and, in consequence, on the profitability of the whole production system. The increasing number of days with extreme hot temperatures suggests that it is imperative to detect the heat stress status of animals based on adequate measures. However, confirming the heat stress status of an individual is still challenging, and, in consequence, the identification of novel heat stress biomarkers, including molecular biomarkers, remains a very relevant issue. Currently, it is known that heat stress seems to have unfavorable effects on immune system mechanisms, but this information is of limited use in the context of heat stress phenotyping. In addition, there is a lack of knowledge addressing the molecular mechanisms linking the relevant genes to the observed phenotype. In this review, we explored the potential molecular mechanisms explaining how heat stress affects the immune system and, therefore, increases the occurrence of immune-related diseases in cattle. In this regard, 2 relatively opposite hypotheses are under focus: the immunosuppressive action of cortisol, and the proinflammatory effect of heat stress. In both hypotheses, the modulation of the immune response during heat stress is highlighted. Moreover, it is possible to link candidate genes to these potential mechanisms. In this context, immune markers are very valuable indicators for the detection of heat stress in dairy cattle, broadening the portfolio of potential biomarkers for heat stress.

Genetic Markers Associated with Milk Production and Thermotolerance in Holstein Dairy Cows Managed in a Heat-Stressed Environment

Dairy production in Holstein cows in a semiarid environment is challenging due to heat stress. Under such conditions, genetic selection for heat tolerance appears to be a useful strategy. The objective was to validate molecular markers associated with milk production and thermotolerance traits in Holstein cows managed in a hot and humid environment. Lactating cows (n = 300) exposed to a heat stress environment were genotyped using a medium-density array including 53,218 SNPs. A genome-wide association study (GWAS) detected six SNPs associated with total milk yield (MY305) that surpassed multiple testing (p < 1.14 × 10^-6). These SNPs were further validated in 216 Holstein cows from two independent populations that were genotyped using the TaqMan bi-allelic discrimination method and qPCR. In these cows, only the SNPs rs8193046, rs43410971, and rs382039214, within the genes TLR4, GRM8, and SMAD3, respectively, were associated (p < 0.05) with MY305, rectal temperature (RT), and respiratory rate. Interestingly, these variables improved as the number of favorable genotypes of the SNPs increased from 0 to 3. In addition, a regression analysis detected RT as a significant predictor (R² = 0.362) for MY305 in cows with >1 favorable genotype, suggesting this close relationship was influenced by genetic markers. In conclusion, SNPs in the genes TLR4, GRM8, and SMAD3 appear to be involved in the molecular mechanism that regulates milk production in cows under heat-stressed conditions. These SNPs are proposed as thermotolerance genetic markers for a selection program to improve the milk performance of lactating Holstein cows managed in a semiarid environment.

Weighted single-step GWAS and RNA sequencing reveals key candidate genes associated with physiological indicators of heat stress in Holstein cattle

Background

The study of molecular processes regulating heat stress response in dairy cattle is paramount for developing mitigation strategies to improve heat tolerance and animal welfare. Therefore, we aimed to identify quantitative trait loci (QTL) regions associated with three physiological indicators of heat stress response in Holstein cattle, including rectal temperature (RT), respiration rate score (RS), and drooling score (DS). We estimated genetic parameters for all three traits. Subsequently, a weighted single-step genome-wide association study (WssGWAS) was performed based on 3200 genotypes, 151,486 phenotypic records, and 38,101 animals in the pedigree file. The candidate genes located within the identified QTL regions were further investigated through RNA sequencing (RNA-seq) analyses of blood samples for four cows collected in April (non-heat stress group) and four cows collected in July (heat stress group).

Results

The heritability estimates for RT, RS, and DS were 0.06, 0.04, and 0.03, respectively. Fourteen, 19, and 20 genomic regions explained 2.94%, 3.74%, and 4.01% of the total additive genetic variance of RT, RS, and DS, respectively. Most of these genomic regions are located in the Bos taurus autosome (BTA) BTA3, BTA6, BTA8, BTA12, BTA14, BTA21, and BTA24. No genomic regions overlapped between the three indicators of heat stress, indicating the polygenic nature of heat tolerance and the complementary mechanisms involved in heat stress response. For the RNA-seq analyses, 2627 genes were significantly upregulated and 369 downregulated in the heat stress group in comparison to the control group. When integrating the WssGWAS, RNA-seq results, and existing literature, the key candidate genes associated with physiological indicators of heat stress in Holstein cattle are: PMAIP1, SBK1, TMEM33, GATB, CHORDC1, RTN4IP1, and BTBD7.

Conclusions

Physiological indicators of heat stress are heritable and can be improved through direct selection. Fifty-three QTL regions associated with heat stress indicators confirm the polygenic nature and complex genetic determinism of heat tolerance in dairy cattle. The identified candidate genes will contribute for optimizing genomic evaluation models by assigning higher weights to genetic markers located in these regions as well as to the design of SNP panels containing polymorphisms located within these candidate genes.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00748-6.

Genome-wide associations for heat stress response suggest potential candidate genes underlying milk fatty acid composition in dairy cattle

Contents of milk fatty acids (FA) display remarkable alterations along climatic gradients. Detecting candidate genes underlying such alterations might be beneficial for the exploration of climate sensitivity in dairy cattle. Consequently, we aimed on the definition of FA heat stress indicators, considering FA breeding values in response to temperature-humidity index (THI) alterations. Indicators were used in GWAS, in ongoing gene annotations and for the estimation of chromosome-wide variance components. The phenotypic data set consisted of 39,600 test-day milk FA records from 5,757 first-lactation Holstein dairy cows kept in 16 large-scale German cooperator herds. The FA traits were C18:0, polyunsaturated fatty acids (PUFA), saturated fatty acids (SFA), and unsaturated fatty acids (UFA). After genotype quality control, 40,523 SNP markers from 3,266 cows and 930 sires were considered. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI, which were allocated to 10 different THI classes. The same FA from 3 stages of lactation were considered as different, but genetically correlated traits. Consequently, a 3-trait reaction norm model was used to estimate genetic parameters and breeding values for FA along THI classes, considering either pedigree (A) or genomic (G) relationship matrices. De-regressed proofs and genomic estimated breeding values at the intermediate THI class 5 and at the extreme THI class 10 were used as pseudophenotypes in ongoing genomic analyses for thermoneutral (TNC) and heat stress conditions (HSC), respectively. The differences in de-regressed proofs and in genomic estimated breeding values from both THI classes were pseudophenotypes for heat stress response (HSR). Genetic correlations between the same FA under TNC and HSC were smallest in the first lactation stage and ranged from 0.20 for PUFA to 0.87 for SFA when modeling with the A matrix, and from 0.35 for UFA to 0.86 for SFA when modeling with the G matrix. In the first lactation stage, larger additive genetic variances under HSC compared with TNC indicate climate sensitivity for C18:0, PUFA, and UFA. Climate sensitivity was also reflected by pronounced chromosome-wide genetic variances for HSR of PUFA and UFA in the first stage of lactation. For all FA under TNC, HSC, and HSR, quite large genetic variance proportions were explained by BTA14. In GWAS, 30 SNP (within or close to 38 potential candidate genes) overlapped for HSR of the different FA. One unique potential candidate gene (AMFR) was detected for HSR of PUFA, 15 for HSR of SFA (ADGRB1, DENND3, DUSP16, EFR3A, EMP1, ENSBTAG00000003838, EPS8, MGP, PIK3C2G, STYK1, TMEM71, GSG1, SMARCE1, CCDC57, and FASN) and 3 for HSR of UFA (ENSBTAG00000048091, PAEP, and EPPK1). The identified unique genes play key roles in milk FA synthesis and are associated with disease resistance in dairy cattle. The results suggest consideration of FA in combination with climatic responses when inferring genetic mechanisms of heat stress in dairy cows.

Genome-wide association study of a thermo-tolerance indicator in pregnant ewes exposed to an artificial heat-stressed environment

Environmental heat stress negatively influences sheep production in warm semi-arid regions. An animal's ability to tolerate warm weather is difficult to measure naturally due to environmental variability and genetic variation between animals. In this study we developed a thermo-tolerance indicator (TTI) to define heat stress tolerance in pregnant sheep in a controlled environment. Next, we performed a genome-wide association study (GWAS) to identify genomic regions and target genes associated with thermo-tolerance in sheep. Pregnant Columbia-Rambouillet crossbred ewes (n = 127) were heat-stressed inside a climate-controlled chamber for 57 days by increasing the temperature-humidity index to ≥30. Rectal temperature (RT) and feed intake (FI) data were collected daily and used for the predictive TTI analysis. After the tenth day of heat stress, the regression analyses revealed that FI was stable; however, when the ewe's RT exceeded 39.8 °C their FI was less than thermo-tolerant ewes. This average predicted temperature was used to classify each ewe as heat stress tolerant (≤39.8 °C) and non-heat stress tolerant (>39.8 °C). A GWAS analysis was performed and genomic regions were compared between heat stress tolerant and non-tolerant ewes. The single-marker genomic analysis detected 16 single nucleotide polymorphisms (SNP) associated with heat stress tolerance (P < 0.0001), whereas the multi-marker Bayesian analysis identified 8 overlapped 1-Mb chromosomal regions accounting for 11.39% of the genetic variation associated with tolerance to heat stress. Four intragenic SNP showed a remarkable contribution to thermo-tolerance, and these markers were within the genes FBXO11 (rs407804467), PHC3 (rs414179061), TSHR (rs418575898) and STAT1 (rs417581105). In conclusion, genomic regions harboring four intragenic SNP were associated with heat stress tolerance, and these candidate genes are proposed to influence heat tolerance in pregnant ewes subjected to an artificially induced warm climate. Moreover, these genetic markers could be suitable for use in further genetic selection programs in sheep managed in semi-arid regions.

Signature selection analysis reveals candidate genes associated with production traits in Iranian sheep breeds

BACKGROUND

Sheep were among the first animals to be domesticated. They are raised all over the world and produce a major scale of animal-based protein for human consumption and play an important role in agricultural economy. Iran is one of the important locations for sheep genetic resources in the world. Here, we compared the Illumina Ovine SNP50 BeadChip data of three Iranian local breeds (Moghani, Afshari and Gezel), as a population that does not undergone artificial breeding programs as yet, and five other sheep breeds namely East Friesian white, East Friesian brown, Lacaune, DorsetHorn and Texel to detect genetic mechanisms underlying economical traits and daptation to harsh environments in sheep.

RESULTS

To identify genomic regions that have been targeted by positive selection, we used fixation index (Fst) and nucleotide diversity (Pi) statistics. Further analysis indicated candidate genes involved in different important traits such as; wool production included crimp of wool (PTPN3, NBEA and KRTAP20-2 genes), fiber diameter (PIK3R4 gene), hair follicle development (LHX2 gene), the growth and development of fiber (COL17A1 gene)), adaptation to hot arid environments (CORIN gene), adaptive in deficit water status (CPQ gene), heat stress (PLCB4, FAM107B, NBEA, PIK3C2B and USP43 genes) in sheep.

CONCLUSIONS

We detected several candidate genes related to wool production traits and adaptation to hot arid environments in sheep that can be applicable for inbreeding goals. Our findings not only include the results of previous researches, but also identify a number of novel candidate genes related to studied traits. However, more works will be essential to acknowledge phenotype- genotype relationships of the identified genes in our study.

Automated Monitoring of Cattle Heat Stress and Its Mitigation

Climate change related global warming is likely to continue, despite all mitigation measures taken by humans, due to the lag effect of long-term anthropogenic activities. Warming of the atmosphere can impact worldwide cattle production directly by compromising health, welfare and productivity, and indirectly by reducing the quality and quantity of animal feed. Under warm thermal conditions, cattle adjust their physiological and behavioural responses as an integral part of thermoregulation to maintain internal body temperature within a safe range. However, a greater intensity and duration of heat exposure can exceed thermoregulatory capacity leading to an increase in internal body temperature beyond the normal limit that ultimately evokes different animal responses to heat. In cattle, response to heat stress can be visually observed as elevated respiration rate or panting, but continuous visual monitoring is labour intensive, time consuming and subjective. Therefore, different weather-based indices have been developed such as the temperature humidity index (THI) and heat load index (HLI) which are commonly used weather-based indices for monitoring cattle heat stress at commercial level. However, the thermal comfort level of cattle based on weather-based indices has limited use at a microclimatic and individual animal level. Varying sensor-based approaches have shown promise to shift the focus of heat stress management to the individual level. Monitoring individual animal response and mitigation strategies for isolated heat-susceptible cattle could save on heat management costs whilst improving animal welfare and productivity. Here we review the technologies that enable automatic, continuous, and real-time cattle heat stress monitoring and mitigation under commercial conditions. Future platforms for autonomous monitoring and mitigation of heat stress in cattle are likely to be based on minimally-invasive smart technologies either singly, or in an integrated system, enabling real-time solutions to animal responses under various production systems and environmental conditions.