Heritability of rectal temperature and genetic correlations with production and reproduction traits in dairy 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.

Genetic analysis of health traits and their associations with longevity, fertility, production, and conformation traits in Holstein cattle

Health traits have high economic values in dairy cattle breeding, which can cause considerable financial loss through involuntary culling. In this study, fourteen health traits were analysed, including five composite health traits: reproductive disorders, udder health (UH), digestive disorders, metabolic disorders, locomotory diseases (LD), and nine independent health traits: gestation disorders and peripartum disorders, irregular estrus cycle and sterility, metritis (ME), mastitis (MA), abomasal displacement (AD), enteritis (EN), and ketosis, claw diseases (CD), laminitis complex. This study analysed variance components for health traits through both single and bivariate repeatability animal models. All health traits showed low heritability, ranging from 0.001 to 0.025. Most of the health traits in five categories showed negative genetic correlations, ranging from -0.012 (CD and EN) to -0.634 (ME and EN). Strong positive genetic correlations appeared within the same category, ranging from 0.469 (EN and AD) to 0.994 (UH and MA, LD and CD). Furthermore, approximate genetic correlations were evaluated between health traits and routinely collected traits (longevity, fertility, production, and conformation). In general, the low to moderate approximate genetic correlations were estimated between health traits and routinely collected traits. The estimated correlations between health traits and longevity, fertility, production, and conformation traits could provide an indirect reference for disease-resistance breeding in Holstein cattle.

Genetic analysis of phenotypic indicators for heat tolerance in crossbred dairy cattle

Climate change-induced rise in global temperatures has intensified heat stress on dairy cattle and is contributing to the generally observed low milk productivity. Selective breeding aimed at enhancing animals' ability to withstand rising temperatures while maintaining optimal performance is crucial for ensuring future access to dairy products. However, phenotypic indicators of heat tolerance are yet to be effectively factored into the objectives of most selective breeding programs. This study investigated the response of milk production to changing heat load as an indication of heat tolerance and the influence of calving season on this response in multibreed dairy cattle performing in three agroecological zones Kenya. First-parity 7-day average milk yield (65 261 milk records) of 1 739 cows were analyzed. Based on routinely recorded weather data that were accessible online, the Temperature-Humidity Index (THI) was calculated and used as a measure of heat load. THI measurements used represented averages of the same 7-day periods corresponding to each 7-day average milk record. Random regression models, including reaction norm functions, were fitted to derive two resilience indicators: slope of the reaction norm (Slope) and its absolute value (Absolute), reflecting changes in milk yield in response to the varying heat loads (THI 50 and THI 80). The genetic parameters of these indicators were estimated, and their associations with average test-day milk yield were examined. There were no substantial differences in the pattern of milk yield response to heat load between cows calving in dry and wet seasons. Animals with ≤50% Bos taurus genes were the most thermotolerant at extremely high heat load levels. Animals performing in semi-arid environments exhibited the highest heat tolerance capacity. Heritability estimates for these indicators ranged from 0.06 to 0.33 and were mostly significantly different from zero (P < 0.05). Slope at THI 80 had high (0.64-0.71) negative correlations with average daily milk yield, revealing that high-producing cows are more vulnerable to heat stress and vice versa. A high (0.63-0.74) positive correlation was observed between Absolute and average milk yield at THI 80. This implied that low milk-producing cows have a more stable milk production under heat-stress conditions and vice versa. The study demonstrated that the slope of the reaction norms and its absolute value can effectively measure the resilience of crossbred dairy cattle to varying heat load conditions. The implications of these findings are valuable in improving the heat tolerance of livestock species through genetic selection.

A comprehensive meta-analysis of genetic parameters for resilience and productivity indicator traits in Holstein cattle

Selection for resilience indicator (RIND) traits in Holstein cattle is becoming an important breeding objective as the worldwide population is expected to be exposed to increased environmental stressors due to both climate change and changing industry standards. However, genetic correlations between RIND and productivity indicator (PIND) traits, which are already being selected for and have the most economic value, are often unfavorable. As a result, it is necessary to fully understand these genetic relationships when incorporating novel traits into selection indices, so that informed decisions can be made to fully optimize selection for both groups of traits. In the past 2 decades, there have been many estimates of RIND traits published in the literature, albeit in small populations. To provide valuable pooled summary estimates, a random-effects meta-analysis was conducted for heritability and genetic correlation estimates for PIND and RIND traits in worldwide Holstein cattle. In total, 926 heritability estimates for 9 PIND and 27 RIND traits, along with 362 estimates of genetic correlation (PIND × RIND traits) were collected. Resilience indicator traits were grouped into the following subgroups: Metabolic Diseases, Hoof Health, Udder Health, Fertility, Heat Tolerance, Longevity, and Other. Pooled estimates of heritability for PIND traits ranged from 0.201 ± 0.05 (energy-corrected milk) to 0.377 ± 0.06 (protein content), while pooled estimates of heritability for RIND traits ranged from 0.032 ± 0.02 (incidence of lameness, incidence of milk fever) to 0.497 ± 0.05 (measures of body weight). Pooled estimates of genetic correlations ranged from -0.360 ± 0.25 (protein content vs. milk acetone concentration) to 0.535 ± 0.72 (measures of fat-to-protein ratio vs. milk acetone concentration). Additionally, out of 243 potential genetic correlations between PIND and RIND traits that could have been reported, only 40 had enough published estimates to implement the meta-analysis model. Our results confirmed that the interactions between PIND and RIND traits are complex, and all relationships should be evaluated when incorporating novel traits into selection indices. This study provides a valuable reference for breeders looking to incorporate RIND traits for Holstein cattle into selection indices.

Optimal sampling interval for characterisation of the circadian rhythm of body temperature in homeothermic animals using periodogram and cosinor analysis

Core body temperature (T c) is a critical aspect of homeostasis in birds and mammals and is increasingly used as a biomarker of the fitness of an animal to its environment. Periodogram and cosinor analysis can be used to estimate the characteristics of the circadian rhythm of T c from data obtained on loggers that have limited memory capacity and battery life. The sampling interval can be manipulated to maximise the recording period, but the impact of sampling interval on the output of periodogram or cosinor analysis is unknown. Some basic guidelines are available from signal analysis theory, but those guidelines have never been tested on T c data. We obtained data at 1-, 5- or 10-min intervals from nine avian or mammalian species, and re-sampled those data to simulate logging at up to 240-min intervals. The period of the rhythm was first analysed using the Lomb-Scargle periodogram, and the mesor, amplitude, acrophase and adjusted coefficient of determination (R 2) from the original and the re-sampled data were obtained using cosinor analysis. Sampling intervals longer than 60 min did not affect the average mesor, amplitude, acrophase or adjusted R 2, but did impact the estimation of the period of the rhythm. In most species, the period was not detectable when intervals longer than 120 min were used. In all individual profiles, a 30-min sampling interval modified the values of the mesor and amplitude by less than 0.1°C, and the adjusted R 2 by less than 0.1. At a 30-min interval, the acrophase was accurate to within 15 min for all species except mice. The adjusted R 2 increased as sampling frequency decreased. In most cases, a 30-min sampling interval provides a reliable estimate of the circadian T c rhythm using periodogram and cosinor analysis. Our findings will help biologists to select sampling intervals to fit their research goals.

Vaginal temperature of lactating cows during heat waves or normal summer day and effect of additional daily cooling treatments as heat load mitigation strategy

The vaginal temperature (VT) of lactating Holstein cows was monitored in not heat wave (NHW) and in heat wave (HW) summer days. Temperature humidity index (THI) was monitored and assigned to four classes of heat load (HL): THI < 68 null; 68 < THI < 74 low; 74 < THI < 80 moderate; and THI > 80 high.Five daily treatments consisting of continuous forced ventilation and sprayed water (1' on follow by 5' off) were assumed as control cooling protocol (CC) and compared with two experimental cooling protocols (EC) applied in the feed bunk and based on the CC plus two additional cooling treatments which lasted a total of 90' (EC90) or 150' (EC150) in the day.Sixty lactating cows were enrolled in two summer trials carried out in NHW or HW. In each trial, 10 cows were cooled by CC, 10 by EC90 and 10 by EC150. Twenty additional cows were monitored in a fall trail to have reference value of THI and VT under thermoneutral conditions (TN). Each trial lasted 72 h, and measurements of VT were carried out by intra-vaginal data loggers.The 33% of observed THI was within the high class of HL during HW, whereas THI never exceeded the upper threshold of moderate or low class of HL in NHW and TN, respectively.Multiparous and high yielding cows were more sensitive to HL, and the increased daily cooling treatments reduced heat load during hot conditions. However, during heat waves a certain degree of hyperthermia occurred even with intensive cooling management based on seven daily treatments.

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.

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.

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.

Impact of heat stress on dairy cattle and selection strategies for thermotolerance: a review

Climate change is a problem that causes many environmental issues that impact the productivity of livestock species. One of the major issues associated with climate change is an increase of the frequency of hot days and heat waves, which increases the risk of heat stress for livestock species. Dairy cattle have been identified as being susceptible to heat stress due to their high metabolic heat load. Studies have shown heat stress impacts several biological processes that can result in large economic consequences. When heat stress occurs, dairy cattle employ several physiological and cellular mechanisms in order to dissipate heat and protect cells from damage. These mechanisms require an increase and diversion in energy toward protection and away from other biological processes. Therefore, in turn heat stress in dairy cattle can lead numerous issues including reductions in milk production and reproduction as well as increased risk for disease and mortality. This indicates a need to select dairy cattle that would be thermotolerant. Various selection strategies to confer thermotolerance have been discussed in the literature, including selecting for reduced milk production, crossbreeding with thermotolerant breeds, selecting based on physiological traits and most recently selecting for enhanced immune response. This review discusses the various issues associated with heat stress in dairy cattle and the pros and cons to the various selection strategies that have been proposed to select for thermotolerance in dairy cattle.

Differences in body temperature regulation during heat stress and seasonal depression in milk yield between Holstein, Brown Swiss, and crossbred cows

It is not clear whether cattle that are genetically superior in regulation of body temperature during heat stress are also better able to sustain milk production during hot conditions. Objectives were to evaluate differences in body temperature regulation during heat stress between Holstein, Brown Swiss, and crossbred cows under semi-tropical conditions and test whether the seasonal depression in milk yield was greater for genetic groups less able to regulate body temperature. For the first objective, conducted during heat stress, vaginal temperature was measured at 15-min intervals for 5 d in 133 pregnant lactating cows. Vaginal temperatures were affected by time and interaction between genetic group and time. Vaginal temperatures were higher for Holsteins for most times of the day. Moreover, the maximum daily vaginal temperature was higher for Holstein (39.8 ± 0.1°C) than for Brown Swiss (39.3 ± 0.2°C) or crossbreds (39.2 ± 0.1°C). For the second objective, 6,179 lactation records from 2,976 cows were analyzed to determine effects of genetic group and season of calving (cool season = Oct to March; warm season = April to Sept) on 305-d milk yield. Milk yield was affected by genetic group and season but not by the interaction of genetic group and season. The difference in average 305-d milk yield between cows calving in cool versus hot weather was 310 kg (4% decrease) for Holstein, 480 kg (7% decrease) for Brown Swiss, and 420 kg (6% decrease) for crossbreds. In conclusion, Brown Swiss and crossbreds regulated body temperature during heat stress better than Holsteins but these breeds were not more resistant to heat stress with respect to milk yield. Thus, genetic differences in thermotolerance are likely to exist that are independent of regulation of body temperature.

Colony stimulating factor 2 protects the preimplantation bovine embryo from heat shock

Heat stress can have severe deleterious effects on embryo development and survival. The present study evaluated whether CSF2 can protect the developmental competence of the bovine embryo following exposure to a heat shock of 41°C at the zygote and morula stages. In the first experiment, putative zygotes and 2-cell embryos were assigned to receive either 10 ng/ml CSF2 or vehicle, and then cultured for 15 h at either 38.5°C or 41°C and then at 38.5°C until day 7.5. Heat shock reduced blastocyst development for embryos treated with vehicle but not for embryos cultured with CSF2. In the second experiment, day 5 embryos (morula) were treated with CSF2 or vehicle and then cultured for 15 h at either 38.5°C or 41°C and then at 38.5°C until day 7.5. Temperature treatment did not affect development to the blastocyst stage and there was no effect of CSF2 treatment or the interaction. Results indicate that CSF2 can reduce the deleterious effects of heat shock at the zygote or two-cell stage when the embryo is transcriptionally inactive.

Phenotypic and genetic parameters of circadian rhythms from core body temperature profiles and their relationships with beef steers' production efficiency profiles during successive winter feeding periods

This 2-year study evaluated differences in circadian parameters obtained from measures of core body temperatures using telemetric reticulo-rumen and rectal devices during two winter feeding regimes in western Canada. The study also estimated phenotypic correlations and genetic parameters associated with circadian parameters and other production traits in each feeding regime. Each year, 80 weaned steer calves (initial age: 209 ± 11 days; BW: 264 ± 20 kg) from the same cohort were tested over two successive regimes, Fall-Winter (FW) and Winter-Spring (WS) at Lanigan, Saskatchewan, Canada. The steers received forage-based rations in both regimes where the individual feed intake was measured with automatic feeding units. During the trial, the reticulo-rumen (RTMP) and rectal (RCT) temperatures were simultaneously measured every 5 min using telemetric devices. These were used to calculate the circadian parameters (Midline Estimating Statistic Of Rhythms, amplitude, and acrophase/peak time) for both temperature measures. Growth and efficiency performance traits were also determined for all steers. Each steer was assigned into inefficient, neutral, and efficient classes based on the SD of the residual feed intake (RFI), residual gain (RG), and residual intake and gain (RIG) within each year and feeding regime. Higher (p < 0.0003) RTMP and rectal temperature MESORs were observed in the Fall-Winter compared to the Winter-Spring regime. While the two test regimes were different (p < 0.05) for the majority of the RTMP or RCT temperature parameters, they did not differ (p > 0.10) with the production efficiency profiles. The heritability estimates were higher in FW (0.78 ± 0.18 vs. 0.56 ± 0.26) than WS (0.50 ± 0.18 vs. 0.47 ± 0.22) for the rumen and rectal MESORs, respectively. There were positive genetic correlations between the two regimes for the RTMP (0.69 ± 0.21) and RCT (0.32 ± 0.59). There was a negative correlation (p < 0.001) between body temperature and ambient temperature. The high heritability estimates and genetic correlations for rumen and rectal temperature parameters demonstrate their potential as beef genetic improvement tools of economic traits associated with the parameters. However, there are limited practical implications of using only the core-body temperature as a proxy for production efficiency traits for beef steers during winter.

Genes and models for estimating genetic parameters for heat tolerance in dairy cattle

Dairy cattle are highly susceptible to heat stress. Heat stress causes a decline in milk yield, reduced dry matter intake, reduced fertility rates, and alteration of physiological traits (e.g., respiration rate, rectal temperature, heart rates, pulse rates, panting score, sweating rates, and drooling score) and other biomarkers (oxidative heat stress biomarkers and stress response genes). Considering the significant effect of global warming on dairy cattle farming, coupled with the aim to reduce income losses of dairy cattle farmers and improve production under hot environment, there is a need to develop heat tolerant dairy cattle that can grow, reproduce and produce milk reasonably under the changing global climate and increasing temperature. The identification of heat tolerant dairy cattle is an alternative strategy for breeding thermotolerant dairy cattle for changing climatic conditions. This review synthesizes information pertaining to quantitative genetic models that have been applied to estimate genetic parameters for heat tolerance and relationship between measures of heat tolerance and production and reproductive performance traits in dairy cattle. Moreover, the review identified the genes that have been shown to influence heat tolerance in dairy cattle and evaluated the possibility of using them in genomic selection programmes. Combining genomics information with environmental, physiological, and production parameters information is a crucial strategy to understand the mechanisms of heat tolerance while breeding heat tolerant dairy cattle adapted to future climatic conditions. Thus, selection for thermotolerant dairy cattle is feasible.

Growth Dynamics on the Skin and the Coat in Normandian and Simmental Cattle During Adaptation to the Farming Technology in Southwestern Part of Bulgaria

The aim of the present study was to identify the adaptability of Normandy cows raised in a loose system to the climatic and technological conditions and to compare their adaptability with already acclimatized Simmental cows in one farm in Southwestern Bulgaria. Total of 20 cows at same age (II – IV lactation) and same body weight (630-660 kg) from both breeds were classified into 2 groups (each having 10 cows) and subjected to skin thickness measurement, fibers total weight determination and fiber categorization during the winter and summer season. The results have shown that during the winter season, the Normandy breed had significantly higher skin thickness at the elbow compared to the Simmental breed, whereas no differences were observed in the skin thickness between both breeds either at the neck or at the middle of the last rib. During the summer season, the differences on the skin thickness compared to the Simmental breed were observed at the neck (p<0.05), whilst at the elbow and at the middle of the last rib no differences were recorded (p>0.05). During the winter seasons, the coat of Normandy cows contained more soft fibers in comparison to the Simmentals cows. Furthermore, the changes in the observed parameters, influenced by climatic conditions showed similar pattern in both breeds. It can be concluded that the two breeds have emphasized their genetic potential and have a good adaptability to the temperate continental climate of the region.

Metabolome profiling of plasma reveals different metabolic responses to acute cold challenge between Inner-Mongolia Sanhe and Holstein cattle

Low-temperature conditions influence cattle productivity and survivability. Understanding the metabolic regulations of specific cattle breeds and identifying potential biomarkers related to cold challenges are important for cattle management and optimization of genetic improvement programs. In this study, 28 Inner-Mongolia Sanhe and 22 Holstein heifers were exposed to -25°C for 1 h to evaluate the differences in metabolic mechanisms of thermoregulation. In response to this acute cold challenge, altered rectal temperature was only observed in Holstein cattle. Further metabolome analyses showed a greater baseline of glycolytic activity and mobilization of AA in Sanhe cattle during normal conditions. Both breeds responded to the acute cold challenge by altering their metabolism of volatile fatty acids and AA for gluconeogenesis, which resulted in increased glucose levels. Furthermore, Sanhe cattle mobilized the citric acid cycle activity, and creatine and creatine phosphate metabolism to supply energy, whereas Holstein cattle used greater AA metabolism for this purpose. Altogether, we found that propionate and methanol are potential biomarkers of acute cold challenge response in cattle. Our findings provide novel insights into the biological mechanisms of acute cold response and climatic resilience, and will be used as the basis when developing breeding tools for genetically selecting for improved cold adaptation in cattle.

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.

Effectiveness of the Australian breeding value for heat tolerance at discriminating responses of lactating Holstein cows to heat stress

Heat stress has negative consequences for milk production and reproduction of dairy cattle. These adverse effects are likely to increase because of climate change and anticipated increases in milk yield. Some of the variation among cows in ability to resist effects of heat stress is genetic. The current objective of this observational study was to assess the effectiveness of the Australian breeding value for heat tolerance (ABVHT) based on the decline in milk yield with heat stress for predicting cow differences in effects of heat stress on regulation of body temperature, milk production, and reproductive function. Genomic breeding values for heat tolerance were calculated for 12,487 cows from a single California dairy farm. Rectal temperature in the afternoon (1100-2045 h) was measured on a subset of 626 lactating cows with ABVHT ≥102 (heat tolerant) or <102 (heat sensitive). Rectal temperature was 0.12°C lower for heat-tolerant cows than heat-sensitive cows. Vaginal temperatures were measured every 15 min for 5 d in 118 cows with ABVHT ≥108 (extreme heat tolerant) or <97 (extreme heat sensitive). Vaginal temperature was 0.07°C lower for extreme heat-tolerant cows than extreme heat-sensitive cows. Lactation records for 4,703 cows with ABVHT were used to evaluate seasonal variation in first 90-d milk yield, fat percent, and protein percent for each ABVHT quartile. Overall, cows with higher ABVHT had lower milk yield, fat percentage, and protein percentage and higher first service pregnancy rate. There was no summer depression in production or reproduction or interactions between season and ABVHT quartile. We observed that ABVHT can successfully identify heat-tolerant cows that maintain lower body temperatures during heat stress. The lack of a pronounced seasonality in milk production or reproduction precluded evaluation of whether ABVHT is related to the magnitude of effect of heat stress on those traits.

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.

Genetic Modeling and Genomic Analyses of Yearling Temperament in American Angus Cattle and Its Relationship With Productive Efficiency and Resilience Traits

Cattle temperament has been considered by farmers as a key breeding goal due to its relevance for cattlemen's safety, animal welfare, resilience, and longevity and its association with many economically important traits (e.g., production and meat quality). The definition of proper statistical models, accurate variance component estimates, and knowledge on the genetic background of the indicator trait evaluated are of great importance for accurately predicting the genetic merit of breeding animals. Therefore, 266,029 American Angus cattle with yearling temperament records (1-6 score) were used to evaluate statistical models and estimate variance components; investigate the association of sex and farm management with temperament; assess the weighted correlation of estimated breeding values for temperament and productive, reproductive efficiency and resilience traits; and perform a weighted single-step genome-wide association analysis using 69,559 animals genotyped for 54,609 single-nucleotide polymorphisms. Sex and extrinsic factors were significantly associated with temperament, including conception type, age of dam, birth season, and additional animal-human interactions. Similar results were observed among models including only the direct additive genetic effect and when adding other maternal effects. Estimated heritability of temperament was equal to 0.39 on the liability scale. Favorable genetic correlations were observed between temperament and other relevant traits, including growth, feed efficiency, meat quality, and reproductive traits. The highest approximated genetic correlations were observed between temperament and growth traits (weaning weight, 0.28; yearling weight, 0.28). Altogether, we identified 11 genomic regions, located across nine chromosomes including BTAX, explaining 3.33% of the total additive genetic variance. The candidate genes identified were enriched in pathways related to vision, which could be associated with reception of stimulus and/or cognitive abilities. This study encompasses large and diverse phenotypic, genomic, and pedigree datasets of US Angus cattle. Yearling temperament is a highly heritable and polygenic trait that can be improved through genetic selection. Direct selection for temperament is not expected to result in unfavorable responses on other relevant traits due to the favorable or low genetic correlations observed. In summary, this study contributes to a better understanding of the impact of maternal effects, extrinsic factors, and various genomic regions associated with yearling temperament in North American Angus cattle.

Associations of Wearable Sensor Measures With Feed Intake, Production Traits, Lactation, and Environmental Parameters Impacting Feed Efficiency in Dairy Cattle

Feed efficiency is an important trait to dairy production because of its impact on sustainability and profitability. Measuring individual cow feed intake on commercial farms would be unfeasibly costly at present. Thus, developing cheap and portable indicators of feed intake would be highly beneficial for genetic selection and precision feeding management tools. Given the growing use of automated sensors on dairy farms, the objective of this study was to determine the relationship between measurements recorded from multiple wearable sensors and feed intake. A total of three different wearable sensors were evaluated for their association with dry mater intake (DMI). The sensors measured activity (sensors = 3), rumination (sensors = 1), ear temperature (sensors = 1), rumen pH (sensors = 1) and rumen temperature (sensors = 1). A range of 56–340 cows with assorted sensors from 24 to 313 days in milk (DIM) were modeled to evaluate associations with DIM, parity, and contemporary group (CG; comprised of pen and study cohort). Models extending upon these variables included known energy sinks (i.e., milk production, milk fat/protein and metabolic body weight), to characterize the association of sensors measures and DMI. Statistically significant (i.e., P &lt; 0.05) regression coefficients for individual sensor measures with DMI ranged from 9.01E-07 to −3.45 kg DMI/day. When integrating all measures from a single sensor in a model, estimated regression coefficients ranged 8.83E-07 to −3.48 kg DMI/day. Significant associations were also identified for milk production traits, parity, DIM and CG. Associations tended to be highest for timepoints around the time of feeding and when multiple measurements within a sensor were integrated in a single model. The findings of this study indicate sensor measures are associated with feed intake and other energy sink traits and variables impacting feed efficiency. This information would be helpful to improve feed and feeding efficiency on commercial farms as proxy measurements for feed intake.

Effects of heat stress on the endometrial epidermal growth factor profile and fertility in dairy cows

The endometrial epidermal growth factor (EGF) profile is an indicator of uterine function and fertility in cattle. The present study aimed to investigate the effects of heat stress on the endometrial EGF profile and fertility in lactating Holstein cows. The endometrial EGF profiles of 365 cows in the Hokkaido and Kyushu regions were examined between June and September (heat stress period, n = 211) and between October and January (control period, n = 154). EGF profiles were investigated using uterine endometrial tissues obtained by biopsy 3 days after estrus (Day 3). The proportion of cows with an altered EGF profile was higher between June and September than between October and January (41.2 vs. 16.2%, P < 0.05). The effects of rectal temperature on Days 0 and 3 on the endometrial EGF profile were also assessed in cows (n = 79) between June and September in the Kyushu region. A single embryo was transferred to cow on Day 7 to evaluate fertility (n = 67). Regardless of the rectal temperature on Day 3, the proportion of cows with an altered EGF profile was higher (64.1 vs. 30.0%, P < 0.05) and the pregnancy rate after embryo transfer (ET) was lower (26.7 vs. 51.4%, P < 0.05) in cows with a rectal temperature ≥ 39.5°C on Day 0 than in cows with a rectal temperature < 39.5°C on Day 0. The present results indicate that alterations in the endometrial EGF profile induced by an elevated body temperature on Day 0 contributed to reductions in fertility in lactating dairy cows during the heat stress period.

The Genetics of Thermoregulation in Pigs: A Review

Heat stress (HS) affects pig performance, health and welfare, resulting in a financial burden to the pig industry. Pigs have a limited number of functional sweat glands and their thermoregulatory mechanisms used to maintain body temperature, are challenged by HS to maintain body temperature. The genetic selection of genotypes tolerant to HS is a promising long-term (adaptation) option that could be combined with other measures at the production system level. This review summarizes the current knowledge on the genetics of thermoregulation in pigs. It also discusses the different phenotypes that can be used in genetic studies, as well as the variability in thermoregulation between pig breeds and the inheritance of traits related to thermoregulation. This review also considers on-going challenges to face for improving heat tolerance in pigs.

Whole-Genome comparative analysis reveals genetic mechanisms of disease resistance and heat tolerance of tropical Bos indicus cattle breeds

The Bos indicus cattle breeds have been naturally selected over thousands of years for disease resistance and thermo-tolerance. However, a genetic mechanism of these specific inherited characteristics needs to be discovered. Hence, in this study, the whole-genome comparative analysis of Bos indicus cattle breeds of Kangayam, Tharparkar, Sahiwal, Red Sindhi, and Hariana of the Indian subcontinent was conducted. The genetic variants identification analysis revealed a total of 15,58,51,012 SNPs and 1,00,62,805 InDels in the mapped reads across all Bos indicus cattle breeds. The functional annotation of 17,252 genes that comprised both, SNPs and InDels, of high functional impact on proteins, has been carried out. The functional annotation results revealed the pathways that were involved in the innate immune response including toll-like receptors, a retinoic acid-inducible gene I like receptors, NOD-like receptors, Jak-STAT signaling pathways, and the non-synonymous variants in the candidate immune genes. Further, we also identified several pathways involved in heat shock response, hair and skin properties, oxidative stress response, osmotic stress response, thermal sweating, feed intake, metabolism, and the non-synonymous variants in the candidate thermo-tolerant genes. These pathways and genes were directly or indirectly contributing to the disease resistance and thermo-tolerance adaptations of Bos indicus cattle breeds.

Body temperature is a repeatable trait in a free-ranging passerine bird

Body temperature (Tb) affects animal function through its influence on rates of biochemical and biophysical reactions, the molecular structures of proteins and tissues, and, ultimately, organismal performance. Despite its importance in driving physiological processes, there are few data on how much variation in Tb exists within populations of organisms, and whether this variation consistently differs among individuals over time (i.e. repeatability of a trait). Here, using thermal radio-frequency identification implants, we quantified the repeatability of Tb, both in the context of a fixed average environment (∼21°C) and across ambient temperatures (6-31°C), in a free-living population of tree swallows (Tachycineta bicolor, n=16). By experimentally trimming the ventral plumage of a subset of female swallows (n=8), we also asked whether the repeatability of Tb is influenced by the capacity to dissipate body heat. We found that both female and male tree swallow Tb was repeatable at 21°C (R=0.89-92), but female Tb was less repeatable than male Tb across ambient temperature (Rfemale=0.10, Rmale=0.58), which may be due to differences in parental investment. Trimmed birds had on average lower Tb than control birds (by ∼0.5°C), but the repeatability of female Tb did not differ as a function of heat dissipation capacity. This suggests that trimmed individuals adjusted their Tb to account for the effects of heat loss on Tb. Our study provides a first critical step toward understanding whether Tb is responsive to natural selection, and for predicting how animal populations will respond to climatic warming.

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.

Methylome Patterns of Cattle Adaptation to Heat Stress

Heat stress has a detrimental impact on cattle health, welfare and productivity by affecting gene expression, metabolism and immune response, but little is known on the epigenetic mechanisms mediating the effect of temperature at the cellular and organism level. In this study, we investigated genome-wide DNA methylation in blood samples collected from 5 bulls of the heat stress resilient Nellore breed and 5 bulls of the Angus that are more heat stress susceptible, exposed to the sun and high temperature-high humidity during the summer season of the Brazilian South-East region. The methylomes were analyzed during and after the exposure by Reduced Representation Bisulfite Sequencing, which provided genome-wide single-base resolution methylation profiles. Significant methylation changes between stressful and recovery periods were observed in 819 genes. Among these, 351 were only seen in Angus, 366 were specific to Nellore, and 102 showed significant changes in methylation patterns in both breeds. KEGG and Gene Ontology (GO) enrichment analyses showed that responses were breed-specific. Interestingly, in Nellore significant genes and pathways were mainly involved in stress responses and cellular defense and were under methylated during heat stress, whereas in Angus the response was less focused. These preliminary results suggest that heat challenge induces changes in methylation patterns in specific loci, which should be further scrutinized to assess their role in heat tolerance.

Climate-Resilient Dairy Cattle Production: Applications of Genomic Tools and Statistical Models

The current changing climate trend poses a threat to the productive efficacy and welfare of livestock across the globe. This review is an attempt to synthesize information pertaining to the applications of various genomic tools and statistical models that are available to identify climate-resilient dairy cows. The different functional and economical traits which govern milk production play a significant role in determining the cost of milk production. Thus, identification of these traits may revolutionize the breeding programs to develop climate-resilient dairy cattle. Moreover, the genotype–environment interaction also influences the performance of dairy cattle especially during a challenging situation. The recent advancement in molecular biology has led to the development of a few biotechnological tools and statistical models like next-generation sequencing (NGS), microarray technology, whole transcriptome analysis, and genome-wide association studies (GWAS) which can be used to quantify the molecular mechanisms which govern the climate resilience capacity of dairy cows. Among these, the most preferred option for researchers around the globe was GWAS as this approach jointly takes into account all the genotype, phenotype, and pedigree information of farm animals. Furthermore, selection signatures can also help to demarcate functionally important regions in the genome which can be used to detect potential loci and candidate genes that have undergone positive selection in complex milk production traits of dairy cattle. These identified biomarkers can be incorporated in the existing breeding policies using genomic selection to develop climate-resilient dairy cattle.

Genetic parameters for rectal temperature, respiration rate, and drooling score in Holstein cattle and their relationships with various fertility, production, body conformation, and health traits

Genetic selection for improved climatic resilience is paramount to increase the long-term sustainability of high-producing dairy cattle, especially in face of climate change. Various physiological indicators, such as rectal temperature (RT), respiration rate score (RR), and drooling score (DS), can be used to genetically identify animals with more effective coping mechanisms in response to heat stress events. In this study, we investigated genetic parameters for RT, RR (score from 1-3), and DS (score from 1-3). Furthermore, we assessed the genetic relationship among these indicators and other economically important traits for the dairy cattle industry. After data editing, 59,265 (RT), 30,290 (RR), and 30,421 (DS) records from 13,592 lactating Holstein cows were used for the analyses. Variance components were estimated based on a multiple-trait repeatability animal model. The heritability ± standard error estimate for RT, RR, and DS was 0.06 ± 0.01, 0.04 ± 0.01, and 0.02 ± 0.01, respectively, whereas their repeatability was 0.19, 0.14, and 0.14, respectively. Moderate genetic correlations of RR with RT and DS (0.26 ± 0.11 and 0.25 ± 0.16) and nonsignificant correlation between RT and DS (-0.11 ± 0.14) were observed. Furthermore, the approximate genetic correlations between RT, RR, and DS with 12 production, 29 conformation, 5 fertility and reproduction, 5 health, and 9 longevity-indicator traits were assessed. In general, the approximate genetic correlations calculated were low to moderate. In summary, 3 physiological indicators of heat stress response were measured in a large number of animals and shown to be lowly heritable. There is a value in developing a selection index including all the 3 indicators to improve heat tolerance in dairy cattle. All the unfavorable genetic relationships observed between heat tolerance and other economically important traits can be accounted for in a selection index to enable improved climatic resilience while also maintaining or increasing productivity in Holstein cattle.

Assessment of Performance and Some Welfare Indicators of Cows in Vietnamese Smallholder Dairy Farms

Smallholder dairy farms (SDFs) are distributed widely across lowland and highland regions in Vietnam, but data on the productivity and welfare status of these cows remains limited. This cross-sectional study was conducted to describe and compare the productivity and welfare status of SDF cows across contrasting regions. It was conducted in autumn 2017 on 32 SDFs randomly selected from four typical but contrasting dairy regions (eight SDFs per region); a south lowland, a south highland, a north lowland, and a north highland region. Each farm was visited over a 24-h period (an afternoon followed by a morning milking and adjacent husbandry activities) to collect data of individual lactating cows (n = 345) and dry cows (n = 123), which included: milk yield and concentrations, body weight (BW), body condition score (BCS, 5-point scale, 5 = very fat), inseminations per conception, and level of heat stress experienced (panting score, 4.5-point scale, 0 = no stress). The high level of heat stress (96% of lactating cows were moderate to highly heat-stressed in the afternoon), low energy corrected milk yield (15.7 kg/cow/d), low percentage of lactating cows (37.3% herd), low BW (498 and 521 kg in lactating and dry cows, respectively), and low BCS of lactating cows (2.8) were the most important productivity and welfare concerns determined and these were most serious in the south lowland. By contrast, cows in the north lowland, a relatively hot but new dairying region, performed similarly to those in the south highland; a region historically considered to be one of the most suitable for dairy cows in Vietnam due to its cool environment. This indicates the potential to mitigate heat stress through new husbandry strategies. Cows in the north highland had the highest BW (535 and 569 kg in lactating and dry cows, respectively) and the highest energy corrected milk yield (19.2 kg/cow/d). Cows in all regions were heat-stressed during the daytime, although less so in the highlands compared to the lowlands. Opportunities for research into improving the productivity and welfare of Vietnamese SDF cows are discussed.

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

Heat stress is a major cause of welfare issues and economic losses to the worldwide dairy cattle industry. Genetic selection for heat tolerance has a great potential to positively affect the dairy industry, as the gains are permanent and cumulative over generations. Rectal temperature (RT) is hypothesized to be a good indicator trait of heat tolerance. Therefore, this study investigated the genetic architecture of RT by estimating genetic parameters, performing genome-wide association studies, and biologically validating potential candidate genes identified to be related to RT in Holstein cattle. A total of 33,013 RT records from 7,598 cows were used in this study. In addition, 1,114 cows were genotyped using the Illumina 150K Bovine BeadChip (Illumina, San Diego, CA). Rectal temperature measurements taken in the morning (AMRT) and in the afternoon (PMRT) are moderately heritable traits, with estimates of 0.09 ± 0.02 and 0.04 ± 0.01, respectively. These 2 traits are also highly genetically correlated (r = 0.90 ± 0.08). A total of 10 SNPs (located on BTA3, BTA4, BTA8, BTA13, BTA14, and BTA29) were found to be significantly associated with AMRT and PMRT. Subsequently, gene expression analyses were performed to validate the key functional genes identified (SPAG17, FAM107B, TSNARE1, RALYL, and PHRF1). This was done through in vitro exposure of peripheral blood mononuclear cells (PBMC) to different temperatures (37°C, 39°C, and 42°C). The relative mRNA expression of 2 genes, FAM107B and PHRF1, significantly changed between the control and heat stressed PBMC. In summary, RT is heritable, and enough genetic variability exists to enable genetic improvement of heat tolerance in Holstein cattle. Important genomic regions were identified and biologically validated; FAM107B and PHRF1 are the main candidate genes identified to influence heat stress response in dairy cattle.

A study on stress response and fertility parameters in phenotypically thermotolerant and thermosensitive dairy cows during summer heat stress

It is well documented that heat stress (HS) causes subfertility in dairy cows. However, during the last ten years we have been observing that, under high temperature-humidity index (THI ≥ 75), despite the overall reduced fertility, some cows conceive at the first artificial insemination (AI). Here, we examined distinctive features of cows with conserved fertility under severe HS. From the databases of three herds, 167 lactating Holstein cows were selected; group TT cows (n = 57) conceived in the previous summer (THI ≥ 75) at the 1st AI, and group TS (n = 110) failed to conceive at the same period after at least 2 consecutive AIs. The animals calved in spring, and in August, blood samples were collected during a hot day (THI ≥ 81) for the determination of cortisol and HSP70 concentrations. In one farm, the validity of fertility data of the previous year was re-examined. In 28 cows from group TT and in 39 cows from group TS, the conception rate was examined during July and August. In 6 cows from each group (TT and TS) the oestrous cycles were synchronized, ovulation was induced with GnRH (THI = 80), and the concentration of the pre-ovulatory LH surge was determined in 9 blood samples. The progesterone concentration in the ensuing cycle was determined in blood samples collected every other day. Overall, cortisol and HSP70 were significantly lower in TT group compared to TS. More (p < .05) animals from group TT conceived at the first AI compared with those from group TS. The induced pre-ovulatory LH surge peaked at higher level (p < .002) in group TT than in group TS, while no difference was recorded among groups either in mean progesterone concentrations or in the duration of the ensuing oestrous cycle. These results are highly suggestive that thermotolerance in some dairy cows is an inherent characteristic, warranting further genetic investigation.

A comprehensive evaluation of microchips to measure temperature in dairy calves

Elevated temperature is often an indicator of an immune response and used in the diagnosis of illness in dairy calves; however, measuring rectal temperature is labor intensive and often not measured daily on the farm. The objective of this study was to measure body temperature using a microchip and determine an appropriate implant site that would passively read body temperature in dairy calves. First, the precision of the temperature microchips and the rectal thermometer were tested ex vivo. Then, Holstein bull calves (n = 12) at 14 ± 12 d (mean ± SD) of age were implanted with microchips subcutaneously by the scutulum of the ear, subcutaneously in the upper scapula (SCAP), and intramuscularly in the trapezius muscle of the neck. One week after implantation, a temperature reading was taken for every microchip implant site using a radiofrequency ID (RFID) reader, as well as rectally and in the tympanic membrane using a digital thermometer every 60 min for 24 h in each calf (hereafter, the hourly study). Additionally, microchip readings and rectal temperatures were taken daily at 0800 h from 8 wk of age (n = 9; 57 ± 12 d of age) until 2 wk postweaning for a subset of the bull calves used in the hourly study (hereafter, the daily study). In the ex vivo trial, the microchip readings were very highly correlated with the rectal thermometer (r = 0.96), and the average coefficient of variation between microchip readings was very low (0.12 ± 0.03%). The relationships between the microchip readings within ear, SCAP, and neck and rectal and tympanic temperatures were analyzed using Pearson correlations and Bland-Altman plots. The ear and neck readings were strongly correlated for the hourly study [individual animal correlation; median (Q1, Q3), r = 0.78 (0.73, 0.84)] and for the daily study [r = 0.79 (0.73, 0.89)] across calves. However, rectal temperature was not significantly correlated with ear, SCAP, neck, or tympanic temperature for the daily and hourly studies. Results suggest that temperature microchips measure temperature appropriately, but temperature is dependent upon the implant site in calves, and temperature measured at ear, SCAP, and neck implant sites cannot be used to estimate rectal temperature. Future research should determine thresholds for fever that are specific to implant sites in calves.

Prospects for gene introgression or gene editing as a strategy for reduction of the impact of heat stress on production and reproduction in cattle

In cattle, genetic variation exists in regulation of body temperature and stabilization of cellular function during heat stress. There are opportunities to reduce the impact of heat stress on cattle production by identifying the causative mutations responsible for genetic variation in thermotolerance and transferring specific alleles that confer thermotolerance to breeds not adapted to hot climates. An example of a mutation conferring superior ability to regulate body temperature is the group of frame-sift mutations in the prolactin receptor gene (PRLR) that lead to a truncated receptor and development of cattle with a short, sleek hair coat. Slick mutations in PRLR have been found in several extant breeds derived from criollo cattle. The slick mutation in Senepol cattle has been introgressed into dairy cattle in Puerto Rico, Florida and New Zealand. An example of a mutation that confers cellular protection against elevated body temperature is a deletion mutation in the promoter region of a heat shock protein 70 gene called HSPA1L. Inheritance of the mutation results in amplification of the transcriptional response of HSPA1L to heat shock and increased cell survival. The case of PRLR provides a promising example of the efficacy of the genetic approach outlined in this paper. Identification of other mutations conferring thermotolerance at the whole-animal or cellular level will lead to additional opportunities for using genetic solutions to reduce the impact of heat stress.

Symposium review: Genomic selection for reducing environmental impact and adapting to climate change

The world has been warming as greenhouse gases accumulate. Worldwide from 1880 to 2012, the average surface temperature has increased by about 0.85°C and by 0.12°C per decade since 1951. The world's cattle population is a contributor to atmospheric methane, a potent greenhouse gas, in addition to suffering from high temperatures combined with humidity. This makes research into reducing the global footprint of dairy cows of importance on a long-term horizon, while improving tolerance to heat could alleviate the effects of rising temperatures. In December 2017, genomic estimated breeding values for heat tolerance in dairy cattle were released for the first time in Australia. Currently, heat tolerance is not included in the Balanced Performance Index (Australia's national selection index), and the correlation between heat tolerance breeding values and Balanced Performance Index is -0.20, so over time, heat tolerance has worsened due to lack of selection pressure. However, in contrast, sizable reductions in greenhouse gas emissions have been achieved as a favorable response to selecting for increased productivity, longevity, and efficiency, with opportunities for even greater gains through selecting for cow emissions directly. Internationally considerable research effort has been made to develop breeding values focused on reducing methane emissions using individual cow phenotypes. This requires (1) definition of breeding objectives and selection criteria and (2) assembling a sufficiently large data set for genomic prediction. Selecting for heat tolerance and reduced emissions directly may improve resilience to changing environments while reducing environmental impact.

The Influence of Temperament on Body Temperature Response to Handling in Angus Cattle

Simple Summary

Understanding animal responses to stressful stimuli is a fundamental aspect to evaluating animal welfare. Stress-induced hyperthermia (SIH) is a term used to describe a short-term increase in body temperature that occurs in response to stressful stimuli. Recently there has been increasing interest in SIH as a physiological measure of psychological stress in livestock species. Previously, studies have suggested that cattle with more excitable temperaments exhibit an increased stress response. This study evaluated the influence of temperament on SIH, during a standardized handling procedure in Bos taurus cattle. In this study, body temperature increased, regardless of sex or temperament traits, characterizing SIH. Nevertheless, both flight speed (FS) and crush score (CS) were associated with an elevated rectal temperature (T_REC) 30 min prior to the handling procedure, and this continued from the start of handling (T0) to 10 min post-handling (T10). The results from this study suggest that temperament may be related to variation in SIH in cattle during handling. Understanding the variation in behavioral and physiological response to stressful events may enable the development of new measures for genetic selection in cattle.

Abstract

Previous studies have indicated that cattle with more excitable temperaments exhibit an increased stress response. The objective of this study was to investigate the relationship between temperament traits, handling, and stress-induced hyperthermia (SIH) in beef cattle. Rectal temperatures (T_REC, °C) of 60 purebred Angus cattle (30 heifers, 30 steers; 235.2 ± 5.11 kg) were recorded at 20 s intervals from 30 min prior to handling until two hours post handling. All cattle were exposed to a standardized handling procedure consisting of (i) being restrained in a weighing box for 30 s; (ii) being held within the crush for 30 s; and then (iii) being restrained in a head bail for 60 s. Cattle temperaments were evaluated via three traits: (1) agitometer score (AG); (2) crush score (CS); and (3) flight speed (FS) during the handling procedure. Agitometer scores and FS measures were used to describe an AG category (AG_CAT) and an FS category (FS_CAT) that were used to classify animals into three temperament categories: 1, calm; 2, intermediate; and 3, temperamental. Pearson’s correlation coefficients were used to evaluate the associations between (i) AG, CS, FS, and T_REC 30 min prior to entry into the weighing box (T-30) and then at 1 min intervals between time of entry into the weighing box (T0) until 10 min post-weighing (T10); and (ii) the relationship between AG, CS, and FS. The relationship between T_REC and temperament traits over the 2.5 h were modeled by using a first-order autoregressive repeated measures model. Flight speed had strong to moderate associations with T_REC at T-30 (r ≥ 0.37; p ≤ 0.006) and between T0 and T10 (r ≥ 0.36; p ≤ 0.01). There were moderate associations amongst T_REC between T0 and T10 and CS (r ≥ 0.31; p ≤ 0.01). A weak relationship existed with CS (r = 0.16; p = 0.16). There were no associations between AG and T_REC at T-30 (r ≥ −0.15; p = 0.84) or between T0 and T10 (r ≤ 0.04; p ≥ 0.4). Rectal temperature, irrespective of sex and temperament traits, was influenced by time (p < 0.0001), and maximum T_REC (39.3 ± 0.04 °C) occurred between 4 and 5.7 min after entry into the weighing box. In addition, CS (p = 0.007) influenced T_REC in these cattle. There were also time × temperament trait × sex interactions with the CS (p = 0.0003) and FS_CAT (p = 0.043) categories; however, time × temperament trait interactions were not statistically significant. Results from this study suggest that cattle with excitable temperaments, as evaluated by FS and CS, have a greater increase in T_REC. In addition, these results suggest that a relationship exists between basal T_REC and FS and CS. Together, these results highlight that temperament, as assessed by FS and CS, influences both basal T_REC and the peak temperature recorded following handling but does not influence the magnitude of change in T_REC post handling.

Genetic Selection for Thermotolerance in Ruminants

Variations in climatic variables (temperature, humidity and solar radiation) negatively impact livestock growth, reproduction, and production. Heat stress, for instance, is a source of huge financial loss to livestock production globally. There have been significant advances in physical modifications of animal environment and nutritional interventions as tools of heat stress mitigation. Unfortunately, these are short-term solutions and may be unsustainable, costly, and not applicable to all production systems. Accordingly, there is a need for innovative, practical, and sustainable approaches to overcome the challenges posed by global warming and climate change-induced heat stress. This review highlights attempts to genetically select and breed ruminants for thermotolerance and thereby sustain production in the face of changing climates. One effective way is to incorporate sustainable heat abatement strategies in ruminant production. Improved knowledge of the physiology of ruminant acclimation to harsh environments, the opportunities and tools available for selecting and breeding thermotolerant ruminants, and the matching of animals to appropriate environments should help to minimise the effect of heat stress on sustainable animal genetic resource growth, production, and reproduction to ensure protein food security.

Genetic parameters for hair characteristics and core body temperature in a multibreed Brahman-Angus herd1

Thermal stress in hot humid conditions limits cattle production. The objectives for this study were to estimate genetic parameters for hair characteristics and core body temperature under low and high temperature humidity index (THI) conditions. Hair samples were collected and measured for length and diameter. Core body temperature was measured as vaginal temperature every 15 min over a 5-d period using an iButton temperature measuring device implanted in a blank CIDR in 336 heifers from the University of Florida multibreed herd (ranging from 100% Angus to 100% Brahman). Restricted maximum likelihood procedures were used to estimate heritabilities from multiple bivariate animal models using the WOMBAT program. Estimates of heritability for hair diameter, undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions were 0.50, 0.67, 0.42, 0.32, and 0.26, respectively. The genetic parameters estimated in this study indicate a large, exploitable genetic variance which can be selected upon to improve tolerance in cattle. Breed effects for differing compositions of Brahman and Angus were also estimated. As Brahman breed composition increased by 25% undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions decreased by 1.32 mm, 2.94 mm, 0.11 °C, and 0.14 °C, respectively. Under both low and high THI conditions, cattle with 25% Brahman breed composition or greater maintained a significantly lower body temperature than the 100% Angus breed group. The incorporation of Brahman germplasm is recommended for herds that often experience heat stress conditions in order to increase resilience to heat stress.

Identification of genetic variants in HSF1 gene and their association with heat tolerance in Murrah buffaloes

The study was carried out to identify Single Nucleotide Polymorphism (SNPs) and their association with thermostolerance traits in 200 murrah Buffalos. Genomic DNA was extracted from frozen/thawed blood samples collected in Beckton-Dickinson vacutainer containing 0.5% (10 μl/ml of blood) anticoagulant EDTA, using phenol-chloroform extraction method. Further the samples were processed for checked quality of DNA on 0.8% agarose gel electrophoresis while its quantification was done using Biospec-nano spectrophotometer method. Custom sequencing results revealed 7 SNPs at the position of A15732G located in intron 4, C17061T in intron 6, C17202T and A17226G in intron 7, G17454A in intron 8, C17605T in exon 9 and T18421C in intron 9 of HSF1 gene sequence. Association analysis showed that the thermal tolerance trait in Murrah buffaloes was significantly affected by three SNP locus, viz. A15732G, C17061T and T18421C. The association among the different genotype of this SNP locus with thermo-tolerance was analyzed using Generalized Linear Model procedure in Statistical Analysis System. Animals of GG genotype at locus A15732 G, TT genotype at locus T18421C and C17061T locus showed lower respiration rate and least HTC was observed in animals belonging to GG genotype at locus A15732G. At linkage, disequilibrium and haplotype construction were analysed using SHEsis software. Haplotypes (49) were constructed, and out of these seven haplotypes (>3 sample size) were considered for association studies. The individuals with Hap6 (ACCAGCC) haplotype combination had lower respiration rate (RR) than other haplotype combinations and these individuals may have better thermal adaptability in comparison to others animals.

Multi-breed genome-wide association studies across countries for electronically recorded behavior traits in local dual-purpose cows

Basic bovine behavior is a crucial parameter influencing cattle domestication. In addition, behavior has an impact on cattle productivity, welfare and adaptation. The aim of the present study was to infer quantitative genetic and genomic mechanisms contributing to natural dual-purpose cow behavior in grazing systems. In this regard, we genotyped five dual-purpose breeds for a dense SNP marker panel from four different European countries. All cows from the across-country study were equipped with the same electronic recording devices. In this regard, we analyzed 97,049 longitudinal sensor behavior observations from 319 local dual-purpose cows for rumination, feeding, basic activity, high active, not active and ear temperature. According to the specific sensor behaviors and following a welfare protocol, we computed two different welfare indices. For genomic breed characterizations and multi-breed genome-wide association studies, sensor traits and test-day production records were merged with 35,826 SNP markers per cow. For the estimation of variance components, we used the pedigree relationship matrix and a combined similarity matrix that simultaneously included both pedigree and genotypes. Heritabilities for feeding, high active and not active were in a moderate range from 0.16 to 0.20. Estimates were very similar from both relationship matrix-modeling approaches and had quite small standard errors. Heritabilities for the remaining sensor traits (feeding, basic activity, ear temperature) and welfare indices were lower than 0.09. Five significant SNPs on chromosomes 11, 17, 27 and 29 were associated with rumination, and two different SNPs significantly influenced the sensor traits “not active” (chromosome 13) and “feeding” (chromosome 23). Gene annotation analyses inferred 22 potential candidate genes with a false discovery rate lower than 20%, mostly associated with rumination (13 genes) and feeding (8 genes). Mendelian randomization based on genomic variants (i.e., the instrumental variables) was used to infer causal inference between an exposure and an outcome. Significant regression coefficients among behavior traits indicate that all specific behavioral mechanisms contribute to similar physiological processes. The regression coefficients of rumination and feeding on milk yield were 0.10 kg/% and 0.12 kg/%, respectively, indicating their positive influence on dual-purpose cow productivity. Genomically, an improved welfare behavior of grazing cattle, i.e., a higher score for welfare indices, was significantly associated with increased fat and protein percentages.

Reproductive physiology of the heat-stressed dairy cow: implications for fertility and assisted reproduction

Heat stress causes a large decline in pregnancy success per insemination during warm times of the year. Improvements in fertility are possible by exploiting knowledge about how heat stress affects the reproductive process. The oocyte can be damaged by heat stress at the earliest stages of folliculogenesis and remains sensitive to heat stress in the peri-ovulatory period. Changes in oocyte quality due to heat stress are the result of altered patterns of folliculogenesis and, possibly, direct effects of elevated body temperature on the oocyte. While adverse effects of elevated temperature on the oocyte have been observed in vitro, local cooling of the ovary and protective effects of follicular fluid may limit these actions in vivo. Heat stress can also compromise fertilization rate. The first seven days of embryonic development are very susceptible to disruption by heat stress. During these seven days, the embryo undergoes a rapid change in sensitivity to heat stress from being very sensitive (2- to 4-cell stage) to largely resistant (by the morulae stage). Direct actions of elevated temperature on the embryo are likely to be an important mechanism for reduction in embryonic survival caused by heat stress. An effective way to avoid effects of heat stress on the oocyte, fertilization, and early embryo is to bypass the effects through embryo transfer because embryos are typically transferred into females after acquisition of thermal resistance. There may be some opportunity to mitigate effects of heat stress by feeding antioxidants or regulating the endocrine environment of the cow but neither approach has been reduced to practice. The best long-term solution to the problem of heat stress may be to increase genetic resistance of cows to heat stress. Thermotolerance genes exist within dairy breeds and additional genes can be introgressed from other breeds by traditional means or gene editing.

Whole Genome Mapping Reveals Novel Genes and Pathways Involved in Milk Production Under Heat Stress in US Holstein Cows

Heat stress represents a major environmental factor that negatively affects the health and performance of dairy cows, causing huge economic losses to the dairy industry. Identifying and selecting animals that are thermotolerant is an attractive alternative for reducing the negative effects of heat stress on dairy cattle performance. As such, the objectives of the present study were to estimate genetic components of milk yield, fat yield, and protein yield considering heat stress and to perform whole-genome scans and a subsequent gene-set analysis for identifying candidate genes and functional gene-sets implicated in milk production under heat stress conditions. Data consisted of about 254k test-day records from 17,522 Holstein cows. Multi-trait repeatability test day models with random regressions on a function of temperature-humidity index (THI) values were used for genetic analyses. The models included herd-test-day and DIM classes as fixed effects, and general and thermotolerance additive genetic and permanent environmental as random effects. Notably, thermotolerance additive genetic variances for all milk traits increased across parities suggesting that cows become more sensitive to heat stress as they age. In addition, our study revealed negative genetic correlations between general and thermotolerance additive effects, ranging between −0.18 to −0.68 indicating that high producing cows are more susceptible to heat stress. The association analysis identified at least three different genomic regions on BTA5, BTA14, and BTA15 strongly associated with milk production under heat stress conditions. These regions harbor candidate genes, such as HSF1, MAPK8IP1, and CDKN1B that are directly involved in the cellular response to heat stress. Moreover, the gene-set analysis revealed several functional terms related to heat shock proteins, apoptosis, immune response, and oxidative stress, among others. Overall, the genes and pathways identified in this study provide a better understanding of the genetic architecture underlying dairy cow performance under heat stress conditions. Our findings point out novel opportunities for improving thermotolerance in dairy cattle through marker-assisted breeding.