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

Impact of Heat Stress on Milk Yield, Milk Fat-to-Protein Ratio, and Conception Rate in Thai-Holstein Dairy Cattle: A Phenotypic and Genetic Perspective

Heat stress severely affects dairy cattle production and reproduction performances in tropical regions. Genetic selection to maintain adequate yield and reproductive performance while enhancing their ability to withstand heat is essential for improving the genetics of dairy cows. Therefore, in this study, we aimed to estimate genetic parameters affecting production and reproduction performances under heat stress conditions in dairy cattle and to investigate the threshold point of heat stress for milk yield (MY), milk fat-to-protein ratio (FPR), and conception rate (CR) in Thai-Holstein dairy cattle. The data included 168,124 records related to MY and milk FPR and 21,278 records of CR in Thai-Holstein dairy cattle, covering the period from 1990 to 2007. A multiple-trait threshold-linear random regression model based on a Bayesian approach via Gibbs sampling was used to estimate variance components, genetic parameters (heritability values, and genetic correlations), and decline rates for each studied trait. The threshold point of heat stress was identified as a temperature and humidity index (THI) of 76. At THI76, a decline was observed in the MY, milk FPR, and CR of Thai dairy cattle. The heritability estimates for MY, milk FPR and CR were 0.347 ± 0.032, 0.293 ± 0.021, and 0.032 ± 0.001, respectively. The genetic correlation between MY and milk FPR and MY and CR were -0.24 and -0.53, respectively, whereas those between milk FPR and heat tolerance as well as between CR and heat tolerance were -0.48 and -0.49, respectively. In addition, the decline rates in MY, milk FPR, and CR were found to be associated with a high percentage of Holstein genetics. In conclusion, the results obtained in this study reveal that the simultaneous consideration of the MY, milk FPR, CR, and heat tolerance traits of Thai-Holstein dairy cattle is possible. In addition, developing a genetic model that incorporates THI is essential for sustainably addressing heat stress problems.

Functional transcriptome analysis revealed major changes in pathways affecting systems biology of Tharparkar cattle under seasonal heat stress

Heat stress significantly disturbs the production, reproduction, and systems biology of dairy cattle. A complex interaction among biological systems helps to combat and overcome heat stress. Indicine cattle breed Tharparkar has been well known for its thermal adaptability. Therefore, present investigation considered RNA-seq technology to explore the functional transcriptomics of Tharparkar cattle with the help of samples collected in spring and summer season. Among differentially expressed genes, about 3280 genes were highly dysregulated, in which 1207 gene were upregulated and 2073 genes were downregulated (|log2fold change|≥ 1 and p ≤ 0.05). Upregulated genes were related to insulin activation, interferons, and potassium ion transport. In contrast, downregulated genes were related to RNA processing, translation, and ubiquitination. Functional annotation revealed that the pathways associated with nervous system (NPFFR1, ROBO3) and metal ion transport (KCNG2, ATP1A2) were highly activated while mRNA processing and translation (EIF4A, EIF4B) and protein processing pathway (VPS4B, PEX13) were highly downregulated. Protein-protein interactions identified hub genes such as ATP13A3, IFNGR2, UBXN7, EIF4A2, SLC12A8 found to play an important role in immune, ubiquitination, translation and transport function. Co-expression network includes LYZ, PNRC1, SQSTM1, EIF4AB and DDX17 genes which are involved in lysosomal activity, tumor inhibition, ubiquitination, and translation initiation. Chemokine signaling pathway associated with immune response was highly upregulated in cluster analysis. The findings of this study provide insights into transcriptome expression and regulation which may better explain complex thermal resilience mechanism of Tharparkar cattle in heat stress under natural conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04018-2.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes

BACKGROUND

Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected.

RESULTS

A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3.

CONCLUSIONS

The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Naturally segregating genetic variants contribute to thermal tolerance in a Drosophila melanogaster model system

Thermal tolerance is a fundamental physiological complex trait for survival in many species. For example, everyday tasks such as foraging, finding a mate, and avoiding predation are highly dependent on how well an organism can tolerate extreme temperatures. Understanding the general architecture of the natural variants within the genes that control this trait is of high importance if we want to better comprehend thermal physiology. Here, we take a multipronged approach to further dissect the genetic architecture that controls thermal tolerance in natural populations using the Drosophila Synthetic Population Resource as a model system. First, we used quantitative genetics and Quantitative Trait Loci mapping to identify major effect regions within the genome that influences thermal tolerance, then integrated RNA-sequencing to identify differences in gene expression, and lastly, we used the RNAi system to (1) alter tissue-specific gene expression and (2) functionally validate our findings. This powerful integration of approaches not only allows for the identification of the genetic basis of thermal tolerance but also the physiology of thermal tolerance in a natural population, which ultimately elucidates thermal tolerance through a fitness-associated lens.

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.

Genetic factors influencing milk and fat yields in tropically adapted dairy cattle: insights from quantitative trait loci analysis and gene associations

OBJECTIVE

The objective of this study was to identify genes associated with 305-day milk yield (MY) and fat yield (FY) that also influence the adaptability of the Thai multibreed dairy cattle population to tropical conditions.

METHODS

A total of 75,776 imputed and actual single nucleotide polymorphisms (SNPs) from 2,661 animals were used to identify genomic regions associated with MY and FY using the single-step genomic best linear unbiased predictions. Fixed effects included herd-yearseason, breed regression, heterosis regression and calving age regression effects. Random effects were animal additive genetic and residual. Individual SNPs with a p-value smaller than 0.05 were selected for gene mapping, function analysis, and quantitative trait loci (QTL) annotation analysis.

RESULTS

A substantial number of QTLs associated with MY (9,334) and FY (8,977) were identified by integrating SNP genotypes and QTL annotations. Notably, we discovered 17 annotated QTLs within the health and exterior QTL classes, corresponding to nine unique genes. Among these genes, Rho GTPase activating protein 15 (ARHGAP15) and catenin alpha 2 (CTNNA2) have previously been linked to physiological traits associated with tropical adaptation in various cattle breeds. Interestingly, these two genes also showed signs of positive selection, indicating their potential role in conferring tolerance to trypanosomiasis, a prevalent tropical disease.

CONCLUSION

Our findings provide valuable insights into the genetic basis of MY and FY in the Thai multibreed dairy cattle population, shedding light on the underlying mechanisms of tropical adaptation. The identified genes represent promising targets for future breeding strategies aimed at improving milk and fat production while ensuring resilience to tropical challenges. This study significantly contributes to our understanding of the genetic factors influencing milk production and adaptability in dairy cattle, facilitating the development of sustainable genetic selection strategies and breeding programs in tropical environments.

Skin transcriptomic analysis reveals candidate genes and pathways associated with thermotolerance in hair sheep

The skin plays an important role in thermoregulation. Identification of genes on the skin that contribute to increased heat tolerance can be used to select animals with the best performance in warm environments. Our objective was to identify candidate genes associated with the heat stress response in the skin of Santa Ines sheep. A group of 80 sheep assessed for thermotolerance was kept in a climatic chamber for 8 days at a stress level temperature of 36 °C (10 am to 04 pm) and a maintenance temperature of 28 °C (04 pm to 10 am). Two divergent groups, with seven animals each, were formed after ranking them by thermotolerance using rectal temperature. From skin biopsy samples, total RNA was extracted, quantified, and used for RNA-seq analysis. 15,989 genes were expressed in sheep skin samples, of which 4 genes were differentially expressed (DE; FDR < 0.05) and 11 DE (FDR 0.05-0.177) between the two divergent groups. These genes are involved in cellular protection against stress (HSPA1A and HSPA6), ribosome assembly (28S, 18S, and 5S ribosomal RNA), and immune response (IGHG4, GNLY, CXCL1, CAPN14, and SAA-4). The candidate genes and main pathways related to heat tolerance in Santa Ines sheep require further investigation to understand their response to heat stress in different climatic conditions and under solar radiation. It is essential to verify whether these genes and pathways are present in different breeds and to understand the relationship between heat stress and other genes identified in this study.

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.

MicroRNA breed and parent-of-origin effects provide insights into biological pathways differentiating cattle subspecies in fetal liver

Introduction: MicroRNAs (miRNAs) play a crucial role in regulating gene expression during key developmental processes, including fetal development. Brahman (Bos taurus indicus) and Angus (Bos taurus taurus) cattle breeds represent two major cattle subspecies with strikingly different phenotypes. Methods: We analyzed miRNA expression in liver samples of purebred and reciprocal crosses of Angus and Brahman to investigate breed and parent-of-origin effects at the onset of accelerated fetal growth. Results: We identified eight novel miRNAs in fetal liver samples and 14 differentially expressed miRNAs (DEMs) between purebred samples. Correlation of gene expression modules and miRNAs by breed and parent-of-origin effects revealed an enrichment of genes associated with breed-specific differences in traits such as heat tolerance (Brahman) and fat deposition (Angus). We demonstrate that genes predicted to be targets of DEMs were more likely to be differentially expressed than non-targets (p-value < 0.05). We identified several miRNAs (bta-miR-187, bta-miR-216b, bta-miR-2284c, bta-miR-2285c, bta-miR-2285cp, bta-miR-2419-3p, bta-miR-2419-5p, and bta-miR-11984) that showed similar correlation patterns as bta-miR-2355-3p, which has been associated with the glutamatergic synapse pathway, a key facilitator of heat tolerance. Furthermore, we report Angus-breed-specific miRNAs (bta-miR-2313-5p, btamiR-490, bta-miR-2316, and bta-miR-11990) that may be involved in fat deposition. Finally, we showed that the DEMs identified in fetal liver are involved in Rap1, MAPK, and Ras signalling pathways, which are important for fetal development, muscle development and metabolic traits such as fat metabolism. Conclusion: Our work sheds light on the miRNA expression patterns that contribute to gene expression differences driving phenotypic differences in indicine and taurine cattle.

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.

Naturally segregating genetic variants contribute to thermal tolerance in a D. melanogaster model system

Thermal tolerance is a fundamental physiological complex trait for survival in many species. For example, everyday tasks such as foraging, finding a mate, and avoiding predation, are highly dependent on how well an organism can tolerate extreme temperatures. Understanding the general architecture of the natural variants of the genes that control this trait is of high importance if we want to better comprehend how this trait evolves in natural populations. Here, we take a multipronged approach to further dissect the genetic architecture that controls thermal tolerance in natural populations using the Drosophila Synthetic Population Resource (DSPR) as a model system. First, we used quantitative genetics and Quantitative Trait Loci (QTL) mapping to identify major effect regions within the genome that influences thermal tolerance, then integrated RNA-sequencing to identify differences in gene expression, and lastly, we used the RNAi system to 1) alter tissue-specific gene expression and 2) functionally validate our findings. This powerful integration of approaches not only allows for the identification of the genetic basis of thermal tolerance but also the physiology of thermal tolerance in a natural population, which ultimately elucidates thermal tolerance through a fitness-associated lens.

Heat Stress: A Serious Disruptor of the Reproductive Physiology of Dairy Cows

Global warming is a significant threat to the sustainability and profitability of the dairy sector, not only in tropical or subtropical regions but also in temperate zones where extreme summer temperatures have become a new and challenging reality. Prolonged exposure of dairy cows to high temperatures compromises animal welfare, increases morbidity, and suppresses fertility, resulting in devastating economic losses for farmers. To counteract the deleterious effects of heat stress, cattl e employ various adaptive thermoregulatory mechanisms including molecular, endocrine, physiological, and behavioral responses. These adaptations involve the immediate secretion of heat shock proteins and cortisol, followed by a complex network of disrupted secretion of metabolic and reproductive hormones such as prolactin, ghrelin, ovarian steroid, and pituitary gonadotrophins. While the strategic heat stress mitigation measures can restore milk production through modifications of the microclimate and nutritional interventions, the summer fertility records remain at low levels compared to those of the thermoneutral periods of the year. This is because sustainment of high fertility is a multifaceted process that requires appropriate energy balance, undisrupted mode of various hormones secretion to sustain the maturation and fertilizing competence of the oocyte, the normal development of the early embryo and unhampered maternal-embryo crosstalk. In this review, we summarize the major molecular and endocrine responses to elevated temperatures in dairy cows, as well as the impacts on maturing oocytes and early embryos, and discuss the consequences that heat stress brings about in dairy cattle fertility.

RNA-sequencing Reveals Differentially Expressed Genes of Laying Hens Fed Baihu Decoction Under Heat Shock

Egg production, an important economic trait in the poultry industry, is sensitive to heat stress. The hypothalamus is a crucial center for thermoregulation by detecting temperature changes and regulating the autonomic nervous system in poultry. Baihu decoction (BH), which contains four ingredients (Rhizoma Anemarrhenae, Gypsum Fibrosum, Radix Glycyrrhizae, and Semen Oryzae Nonglutinosae), is a traditional Chinese medicinal formula for clearing heat. Our study aimed to investigate the changes in gene transcription levels in the hypothalamus of laying hens treated with heat stress with and without BH using RNA sequencing. A total of 223 differentially expressed genes (DEGs) were identified in the heat-treated group compared with the control group and 613 DEGs were identified in the BH group compared with the heat-treated group. Heat shock led to significant changes in the expression of multiple genes involved in the "neuroactive ligand-receptor interaction" pathway. Moreover, feeding BH led to significant upregulation in the expression of eight genes encoding heat shock proteins (HSPs), which were highlighted as candidates to control the "protein processing in the endoplasmic reticulum (ER)" pathway. These results provide the novel insight that BH responds to heat stress by participating in regulation of the ER signaling pathway and HSPs expression.

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.

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.

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.

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

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

Runs of Homozygosity Analysis Reveals Genomic Diversity and Population Structure of an Indigenous Cattle Breed in Southwest China

In aiming to achieve sustainable development goals in the livestock industry, it is becoming increasingly necessary and important for the effective conservation of genetic resources. There are some indigenous cattle breeds in Sichuan, southwest China, for which, however, the genetic diversity and population structures still remain unknown because of the unavailability of systematic breeding programs and pedigree information. Xieka cattle are an indigenous breed locally distributed in southeastern Sichuan and have a long-term evolutionary adaptation to local environments and climates. In this study, we obtained 796,828 single nucleotide polymorphisms (SNPs) through sequencing the genomes of 30 Xieka cattle and used them for analyzing the genetic diversity and runs of homozygosity (ROH). The mean nucleotide diversity was 0.28 and 72% of SNPs were found to be in the heterozygous states. A total of 4377 ROH were detected with even distribution among all autosomes, and 74% of them were lower than 1 Mb in length. Meanwhile, only five ROH were found longer than 5 Mb. We further determined 19 significant genomic regions that were obviously enriched by ROH, in which 35 positional candidate genes were found. Some of these genes have been previously reported to be significantly associated with various production traits in cattle, such as meat quality, carcass performances, and diseases. In conclusion, the relatively high degree of genetic diversity of Xieka cattle was revealed using the genomic information, and the proposed candidate genes will help us optimize the breeding programs regarding this indigenous breed.

Genomic analyses of claw disorders in Holstein cows: Genetic parameters, trait associations, and genome-wide associations considering interactions of SNP and heat stress

The aim of the present study was an in-depth genomic analysis to understand the genomic mechanisms of the 3 claw disorders dermatitis digitalis (DD), interdigital hyperplasia (HYP), and sole ulcer (SU). In this regard, we estimated genetic parameters based on genomic relationship matrices, performed genome-wide association studies, annotated potential candidate genes, and inferred genetic associations with breeding goal traits considering the most important chromosomal segments. As a further novelty of this study, we inferred possible SNP × heat stress interactions for claw disorders. The study consisted of 17,264 first-lactation Holstein Friesian cows kept in 50 large-scale contract herds. The disease prevalence was 15.96, 2.36, and 8.20% for DD, HYP, and SU, respectively. The remaining breeding goal traits consisted of type traits of the feet and leg composite, female fertility, health traits, and 305-d production traits. The final genotype data set included 44,474 SNPs from the 17,264 genotyped cows. Heritabilities for DD, HYP, and SU were estimated in linear and threshold models considering the genomic relationship matrix (G matrix). Genetic correlations with breeding goal traits based on G were estimated in a series of bivariate linear models, which were verified via SNP effect correlations for specific chromosome segments (i.e., segments harboring potential candidate genes for DD, HYP, and SU). Genome-wide association studies were performed for all traits in a case-control design by applying a single SNP linear mixed model. Furthermore, for DD, HYP, and SU, we modeled SNP × heat stress interactions in genome-wide association studies. Single nucleotide polymorphism-based heritabilities were 0.04 and 0.08 for DD, 0.03 and 0.10 for SU, and 0.03 and 0.23 for HYP from linear and threshold models, respectively. The genetic correlations between DD, HYP, and SU with conformation traits from the feet and leg composite were positive throughout, indicating the value of indirect selection on conformation traits to improve claw health. Genetic correlations between DD, SU, and HYP with other breeding goal traits indicated impaired female fertility, impaired udder health status, and productivity decline of diseased cows. Genetic correlations among DD, SU, and HYP were moderate to large, indicating that different claw disorders have similar genetic mechanisms. Nevertheless, we identified disease-specific potential candidate genes, and genetic associations based on the surrounding SNPs partly differed from the genetic correlations. Especially for candidate genes contributing to 2 traits simultaneously, correlations based on SNP effects from the respective chromosome segment were close to 1 or to -1. In this regard, we annotated the candidate genes KRT33A and KRT33B for HYP and DD, KIF27 for HYP and calving to first insemination, and MAN1A1 for SU and the production traits. For SNP × heat stress interactions, we identified significant SNPs on BTA 2, 4, 5, 7, 8, 9, 13, 22, 25, and 28, and we annotated the potential candidate genes FSIP2, CLCN1, ADGRV1, DOP1A, THBD, and RHOBTB1. Results indicate gene-specific mechanisms of the claw disorders only in specific environments.

Direct genetic effects, maternal genetic effects, and maternal genetic sensitivity on prenatal heat stress for calf diseases and corresponding genomic loci in German Holsteins

The aim of this study was to infer the effects of heat stress (HS) of dams during late gestation on direct and maternal genetic parameters for pneumonia (PNEU, 112,563 observations), diarrhea (DIAR, 176,904 observations), and omphalitis (OMPH, 176,872 observations) in Holstein calves kept in large-scale co-operator herds. The genotype dataset included 41,135 SNPs from 19,247 male and female cattle. Temperature-humidity indices (THI) during the last 8 wk of pregnancy were calculated, using the climate data from the nearest public weather station for each herd. Heat load effects were considered for average weekly THI larger than 60. Phenotypically, regression coefficients of calf diseases on prenatal THI during the last 8 wk of gestation were estimated in 8 consecutive runs. The strongest detrimental effects of prenatal HS on PNEU and DIAR were identified for the last week of pregnancy (wk 1). Thus, only wk 1 was considered in ongoing genetic and genomic analyses. In an advanced model considering prenatal HS, random regression coefficients on THI in wk 1 nested within maternal genetic effects (maternal slope effects for heat load) were considered as parameters to infer maternal sensitivity in response to prenatal THI alterations. Direct heritabilities from the advanced model ranged from 0.10 (THI 60) to 0.08 (THI 74) for PNEU and were close to 0.16 for DIAR. Maternal heritabilities for PNEU increased from 0.03 to 0.10 along the THI gradient. For DIAR, the maternal heritability was largest (0.07) at the minimum THI (THI = 60) and decreased to 0.05 at THI 74. Genetic correlations smaller than 0.80 for PNEU and DIAR recorded at THI 60 with corresponding diseases at THI 74 indicated genotype by climate interactions for maternal genetic effects. Genome-wide associations studies were performed using de-regressed proofs of genotyped sires for direct genetic, maternal genetic, and maternal slope effects. Thirty suggestive and 2 significant SNPs were identified from the GWAS. Forty-three genes located close to the suggestive SNPs (±100 kb) were annotated as potential candidate genes. Three biological processes were inferred on the basis of the these genes, addressing the negative regulation of the viral life cycle, innate immune response, and protein ubiquitination. Hence, the genetics of prenatal heat stress mechanisms are associated with immune physiology and disease resistance mechanisms.

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.

Identification of important genomic footprints using eight different selection signature statistics in domestic cattle breeds

In the present study, the population genomic data of different cattle breeds were explored to decipher the genomic regions affected due to selective events and reflected in the productive, reproductive, thermo-tolerance, and health-related traits. To find out these genomic deviations due to selective sweeps, we used eight different statistical tools (Tajima's D, Fu & Li's D*, CLR, ROH, iHS, F_ST, FLK, and hapFLK) on seven indigenous and five exotic cattle breeds. We further performed composite analysis by comparing their covariance matrix. Several candidate genes were found to be related to milk production (ADARB, WDR70, and CA8), reproductive (PARN, FAM134B2, and ZBTB20), and health-related traits (SP110, CXCL2, CLXCL3, CXCL5, IRF8, and MYOM1). The outcome of this investigation provides a basis for detecting selective sweeps that explain the genetic variation of traits. They may possess functional importance for multiple cattle breeds in different subcontinents. However, further studies are required to improve the findings using high-density arrays or whole-genome sequencing with higher resolution and greater sample sizes.

Functionally prioritised whole-genome sequence variants improve the accuracy of genomic prediction for heat tolerance

BACKGROUND

Heat tolerance is a trait of economic importance in the context of warm climates and the effects of global warming on livestock production, reproduction, health, and well-being. This study investigated the improvement in prediction accuracy for heat tolerance when selected sets of sequence variants from a large genome-wide association study (GWAS) were combined with a standard 50k single nucleotide polymorphism (SNP) panel used by the dairy industry.

METHODS

Over 40,000 dairy cattle with genotype and phenotype data were analysed. The phenotypes used to measure an individual's heat tolerance were defined as the rate of decline in milk production traits with rising temperature and humidity. We used Holstein and Jersey cows to select sequence variants linked to heat tolerance. The prioritised sequence variants were the most significant SNPs passing a GWAS p-value threshold selected based on sliding 100-kb windows along each chromosome. We used a bull reference set to develop the genomic prediction equations, which were then validated in an independent set of Holstein, Jersey, and crossbred cows. Prediction analyses were performed using the BayesR, BayesRC, and GBLUP methods.

RESULTS

The accuracy of genomic prediction for heat tolerance improved by up to 0.07, 0.05, and 0.10 units in Holstein, Jersey, and crossbred cows, respectively, when sets of selected sequence markers from Holstein cows were added to the 50k SNP panel. However, in some scenarios, the prediction accuracy decreased unexpectedly with the largest drop of - 0.10 units for the heat tolerance fat yield trait observed in Jersey cows when 50k plus pre-selected SNPs from Holstein cows were used. Using pre-selected SNPs discovered on a combined set of Holstein and Jersey cows generally improved the accuracy, especially in the Jersey validation. In addition, combining Holstein and Jersey bulls in the reference set generally improved prediction accuracy in most scenarios compared to using only Holstein bulls as the reference set.

CONCLUSIONS

Informative sequence markers can be prioritised to improve the genomic prediction of heat tolerance in different breeds. In addition to providing biological insight, these variants could also have a direct application for developing customized SNP arrays or can be used via imputation in current industry SNP panels.

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.

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

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