Characterization of genes and pathways that respond to heat stress in Holstein calves through transcriptome analysis

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.

Skin transcriptome analysis in Brangus cattle under heat stress

Heat stress is a major factor that negatively affects animal welfare and production systems. Livestock should adapt to tropical and subtropical areas and to meet this, composite breeds have been developed. This work aimed to evaluate gene expression profiles in the skin of Brangus cattle under heat stress using a case-control design, and to correlate this with skin histological characteristics. Two groups of bulls were set using rectal temperature as a criterion to define stress conditions: stressed (N = 5) and non-stressed (N = 5) groups. Skin transcriptomics was performed and correlations between breed composition, phenotypic and skin histological traits were evaluated. Results showed 4309 differentially expressed genes (P < 0.01), 2113 downregulated and 2196 upregulated. Enrichment and ontology analyses revealed 132 GO terms and 67 pathways (P < 0.01), including thermogenesis, glycolysis, gluconeogenesis, mitochondrial activity, antioxidant and immune response, and apoptosis. The identity of the terms and pathways indicated the diversity of mechanisms directed to relieve the animals' suffering, acting from simple passive mechanisms (conduction, convection and radiation) to more complex active ones (behavioural changes, evaporation, vasodilation and wheezing). Furthermore, significant differences between phenotypic and skin histological traits and correlations between pairs of traits suggested a direction towards heat dissipation processes. In this sense, number of vessels was positively correlated with number of sweat glands (P < 0.001) and both were positively correlated with zebuine genetic content (P < 0.05 and P < 0.01, respectively), gland size was positively correlated with epidermal thickness and negatively with hair length (P < 0.05), and epidermal thickness was negatively correlated with gland-epidermis distance (P < 0.0005). These results support the notion that response to heat stress is physiologically complex, producing significant changes in the expression of genes involved in several biological pathways, while the animal's ability to face it depends greatly on their skin features.

Genomic insights into the population history and adaptive traits of Latin American Criollo cattle

Criollo cattle, the descendants of animals brought by Iberian colonists to the Americas, have been the subject of natural and human-mediated selection in novel tropical agroecological zones for centuries. Consequently, these breeds have evolved distinct characteristics such as resistance to diseases and exceptional heat tolerance. In addition to European taurine (Bos taurus) ancestry, it has been proposed that gene flow from African taurine and Asian indicine (Bos indicus) cattle has shaped the ancestry of Criollo cattle. In this study, we analysed Criollo breeds from Colombia and Venezuela using whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) array data to examine population structure and admixture at high resolution. Analysis of genetic structure and ancestry components provided evidence for African taurine and Asian indicine admixture in Criollo cattle. In addition, using WGS data, we detected selection signatures associated with a myriad of adaptive traits, revealing genes linked to thermotolerance, reproduction, fertility, immunity and distinct coat and skin coloration traits. This study underscores the remarkable adaptability of Criollo cattle and highlights the genetic richness and potential of these breeds in the face of climate change, habitat flux and disease challenges. Further research is warranted to leverage these findings for more effective and sustainable cattle breeding programmes.

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.

Climate Stressors and Physiological Dysregulations: Mechanistic Connections to Pathologies

This review delves into the complex relationship between environmental factors, their mechanistic cellular and molecular effects, and their significant impact on human health. Climate change is fueled by industrialization and the emission of greenhouse gases and leads to a range of effects, such as the redistribution of disease vectors, higher risks of disease transmission, and shifts in disease patterns. Rising temperatures pose risks to both food supplies and respiratory health. The hypothesis addressed is that environmental stressors including a spectrum of chemical and pathogen exposures as well as physical and psychological influences collectively impact genetics, metabolism, and cellular functions affecting physical and mental health. The objective is to report the mechanistic associations linking environment and health. As environmental stressors intensify, a surge in health conditions, spanning from allergies to neurodegenerative diseases, becomes evident; however, linkage to genetic-altered proteomics is more hidden. Investigations positing that environmental stressors cause mitochondrial dysfunction, metabolic syndrome, and oxidative stress, which affect missense variants and neuro- and immuno-disorders, are reported. These disruptions to homeostasis with dyslipidemia and misfolded and aggregated proteins increase susceptibility to cancers, infections, and autoimmune diseases. Proposed interventions, such as vitamin B supplements and antioxidants, target oxidative stress and may aid mitochondrial respiration and immune balance. The mechanistic interconnections of environmental stressors and disruptions in health need to be unraveled to develop strategies to protect public health.

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.

Heat stress adaptation in cows - Physiological responses and underlying molecular mechanisms

Heat stress is a key abiotic stressor for dairy production in the tropics which is further compounded by the ongoing climate change. Heat stress not only adversely impacts the production and welfare of dairy cows but severely impacts the economics of dairying due to production losses and increased cost of rearing. Over the years, selection has ensured development of high producing breeds, however, the thermotolerance ability of animals has been largely overlooked. In the past decade, the ill effects of climate change have made it pertinent to rethink the selection strategies to opt for climate resilient breeds, to ensure optimum production and reproduction. This has led to renewed interest in evaluation of the impacts of heat stress on cows and the underlying mechanisms that results in their acclimatization and adaptation to varied thermal ambience. The understanding of heat stress and associated responses at various level of animal is crucial to device amelioration strategies to secure optimum production and welfare of cows. With this review, an effort has been made to provide an overview on temperature humidity index as an important indicator of heat stress, general effect of heat stress in dairy cows, and impact of heat stress and subsequent response at physiological, haematological, molecular and genetic level of dairy cows.

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.

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.

Genome-wide expression analysis reveals different heat shock responses in indigenous (Bos indicus) and crossbred (Bos indicus X Bos taurus) cattle

Environmental heat stress in dairy cattle leads to poor health, reduced milk production and decreased reproductive efficiency. Multiple genes interact and coordinate the response to overcome the impact of heat stress. The present study identified heat shock regulated genes in the peripheral blood mononuclear cells (PBMC). Genome-wide expression patterns for cellular stress response were compared between two genetically distinct groups of cattle viz., Hariana (B. indicus) and Vrindavani (B. indicus X B. taurus). In addition to major heat shock response genes, oxidative stress and immune response genes were also found to be affected by heat stress. Heat shock proteins such as HSPH1, HSPB8, FKB4, DNAJ4 and SERPINH1 were up-regulated at higher fold change in Vrindavani compared to Hariana cattle. The oxidative stress response genes (HMOX1, BNIP3, RHOB and VEGFA) and immune response genes (FSOB, GADD45B and JUN) were up-regulated in Vrindavani whereas the same were down-regulated in Hariana cattle. The enrichment analysis of dysregulated genes revealed the biological functions and signaling pathways that were affected by heat stress. Overall, these results show distinct cellular responses to heat stress in two different genetic groups of cattle. This also highlight the long-term adaptation of B. indicus (Hariana) to tropical climate as compared to the crossbred (Vrindavani) with mixed genetic makeup (B. indicus X B. taurus).

Molecular regulation, breed differences and genes involved in stress control in farm animals

Stress is a state of disturbed homeostasis evoking a multiplicity of somatic and mental adaptive reactions resulting from any of the 5 freedoms of animals being violated. Many environmental forces disrupt homeostasis in farm animals, such as extreme temperatures, poor nutrition, noise, hunger, and thirst. During stressful situations, neuronal circuits in the limbic system and prefrontal cortex are activated, which lead to the release of adrenalin and noradrenalin. The hormones released during stress are needed for adaptation to acute stress and are regulated by many genes. This review examined molecular regulation, breed differences, and genes involved in stress control in farm animals. Major molecular regulation of stress, such as oxidative, cytosolic heat shock, unfolded protein, and hypoxic responses, were discussed. The responses of various poultry, ruminant, and pig breeds to different stress types were also discussed. Gene expressions and polymorphisms in the neuroendocrine and neurotransmitter pathways were also elucidated. The information obtained from this review will help farmers mitigate stress in farm animals through appropriate breed and gene-assisted selection. Also, information obtained from this review will add to the field of stress genetics since stress is a serious welfare issue in farm animals.

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.

Metabolomic and transcriptomic study to understand changes in metabolic and immune responses in steers under heat stress

Heat stress (HS) damages livestock by adversely affecting physiological and immunological functions. However, fundamental understanding of the metabolic and immunological mechanisms in animals under HS remains elusive, particularly in steers. To understand the changes on metabolic and immune responses in steers under HS condition, we performed RNA-sequencing and proton nuclear magnetic resonance spectroscopy-based metabolomics on HS-free (THI value: 64.92 ± 0.56) and HS-exposed (THI value: 79.13 ± 0.56) Jersey steer (n = 8, body weight: 559.67 ± 32.72 kg). This study clarifies the metabolic changes in 3 biofluids (rumen fluid, serum, and urine) and the immune responses observed in the peripheral blood mononuclear cells of HS-exposed steers. This integrated approach allowed the discovery of HS-sensitive metabolic and immunological pathways. The metabolomic analysis indicated that HS-exposed steers showed potential HS biomarkers such as isocitrate, formate, creatine, and riboflavin (P < 0.05). Among them, there were several integrative metabolic pathways between rumen fluid and serum. Furthermore, HS altered mRNA expression and immune-related signaling pathways. A meta-analysis revealed that HS decreased riboflavin metabolism and the expression of glyoxylate and dicarboxylate metabolism-related genes. Moreover, metabolic pathways, such as the hypoxia-inducible factor-1 signaling pathway, were downregulated in immune cells by HS (P < 0.05). These findings, along with the datasets of pathways and phenotypic differences as potential biomarkers in steers, can support more in-depth research to elucidate the inter-related metabolic and immunological pathways. This would help suggest new strategies to ameliorate the effects of HS, including disease susceptibility and metabolic disorders, in Jersey steers.

Identification of Novel mRNA Isoforms Associated with Acute Heat Stress Response Using RNA Sequencing Data in Sprague Dawley Rats

The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This study identified the differentially expressed mRNA isoforms which narrowed down the most reliable DEG markers and molecular pathways that underlie the mechanisms of thermoregulation. This experiment was performed on Sprague Dawley rats housed at 22 °C (control group; CT), and three acute heat-stressed groups housed at 42 °C for 30 min (H30), 60 min (H60), and 120 min (H120). Earlier, we demonstrated that acute heat stress increased the rectal temperature of rats, caused abnormal changes in the blood biochemical parameters, as well as induced dramatic changes in the expression levels of genes through epigenetics and post-transcriptional regulation. Transcriptomic analysis using RNA-Sequencing (RNA-Seq) data obtained previously from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5764, and 4988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3549 and 120 in the liver, as well as 3078 and 220 in adrenal glands, respectively. About 35 genes were involved in the heat stress response, out of which, Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3, and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between the liver and adrenal glands. In conclusion, this study would enhance our understanding of the complex underlying mechanisms of acute heat stress, and the identified mRNA isoforms and genes can be used as potential candidates for thermotolerance selection in mammals.

Molecular responses of a key Antarctic species to sedimentation due to rapid climate change

Rapid regional warming causing glacial retreat and melting of ice caps in Antarctica leads benthic filter-feeders to be exposed to periods of food shortage and high respiratory impairment as a consequence of seasonal sediment discharge in the West Antarctic Peninsula coastal areas. The molecular physiological response and its fine-tuning allow species to survive acute environmental stress and are thus a prerequisite to longer-term adaptation to changing environments. Under experimental conditions, we analyzed here the metabolic response to changes in suspended sediment concentrations, through transcriptome sequencing and enzymatic measurements in a highly abundant Antarctic ascidian. We found that the mechanisms underlying short-term response to sedimentation in Cnemidocarpa verrucosa sp. A involved apoptosis, immune defense, and general metabolic depression. These mechanisms may be understood as an adaptive protection against sedimentation caused by glacial retreat. This process can strongly contribute to the structuring of future benthic filter-feeder communities in the face of climate change.

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.

Epigenomics of rats' liver and its cross-species functional annotation reveals key regulatory genes underlying short term heat-stress response

Molecular responses to heat stress are multifaceted and under a complex cellular post-transcriptional control. This study explores the epigenetic and transcriptional alterations induced by heat stress (42 °C for 120 min) in the liver of rats, by integrating ATAC-seq, RNA-Seq, and WGBS information. Out of 2586 differential ATAC-seq peaks induced by heat stress, 36 up-regulated and 22 down-regulated transcript factors (TFs) are predicted, such as Cebpα, Foxa2, Foxo4, Nfya and Sp3. Furthermore, 150,189 differentially methylated regions represent 2571 differentially expressed genes (DEGs). By integrating all data, 45 DEGs are concluded as potential heat stress response markers in rats. To comprehensively annotate and narrow down predicted markers, they are integrated with GWAS results of heat stress parameters in cows, and PheWAS data in humans. Besides better understanding of heat stress responses in mammals, INSR, MAPK8, RHPN2 and BTBD7 are proposed as candidate markers for heat stress in mammals.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Time Series Ovarian Transcriptome Analyses of the Porcine Estrous Cycle Reveals Gene Expression Changes during Steroid Metabolism and Corpus Luteum Development

Simple Summary

The estrous cycle, which is divided into follicular and luteal phases based on ovulation, is influenced by reproductive hormones which affect reproduction and cause changes in the reproductive system of the pig. As the main reproductive organ, the ovary is involved in ovulation and changes in the corpus luteum. We aimed to identify dynamic changes in gene expression through differentially expressed gene profiling and to provide a comprehensive understanding of the molecular mechanisms that occur in the pig ovary during the estrous cycle. The transcriptome analysis revealed that the dynamic change in gene expression was more activated in the luteal phase than in the follicular phase. Functional analysis revealed that the metestrus and diestrus periods are important in preparation for pregnancy or the next estrous cycle after ovulation.

Abstract

The porcine estrous cycle is influenced by reproductive hormones, which affect porcine reproduction and result in physiological changes in the reproductive organs. The ovary is involved in ovulation, luteinization, corpus luteum development, and luteolysis. Here, we aimed to provide a comprehensive understanding of the gene expression patterns in porcine ovarian transcriptomes during the estrous cycle through differentially expressed genes profiling and description of molecular mechanisms. The transcriptomes of porcine ovary were obtained during the estrous cycle at three-day intervals from day 0 to day 18 using RNA-seq. At seven time points of the estrous cycle, 4414 DEG were identified; these were classified into three clusters according to their expression patterns. During the late metestrus and diestrus periods, the expression in cluster 1 increased rapidly, and steroid biosynthesis was significant in the pathway. Cluster 2 gene expression patterns represented the cytokine–cytokine receptor interaction in significant pathways. In cluster 3, the hedgehog signaling pathway was selected as the significant pathway. Our study exhibited dynamic gene expression changes with these three different patterns of cluster 1, 2, and 3. The results helped identify the functions and related significant genes especially during the late metestrus and diestrus periods in the estrous cycle.

Transcriptome analysis and potential mechanisms of bovine oocytes under seasonal heat stress

Heat stress is the major factor affecting cattle fertility but molecular mechanisms of deleterious impacts of elevated temperature on oocyte are still not well known. Therefore, the aim of this study is to gain a better understanding of the underlying molecular mechanism of how heat stress affects GV-stage and MII-stage oocytes and discover hub genes to heat resistance for cow oocytes. In this study, we used the bioinformatics approach to discover the differentially expressed genes between GV-stage and MII-stage oocytes, which were collected during spring and summer. When GV-stage oocytes were compared to MII-stage oocytes collected in July (Jul DEGs group) a total of 1068 genes were found as differentially expressed as a result of heat stress. Also, HSPA8, COPS5, POLR2L, PSMC6, and TPI1 were identified as the common top ranked genes for the Jul DEGs group. The highest connected hub gene for the Jul DEGs group was determined as HSPA8. Our results showed that different heat response mechanisms might be activated to protect oocytes from elevated temperatures in cattle. The identified genes and their associated pathways might play an important role in the response to heat stress that affects the oocytes in cattle.

The Role of microRNAs in the Mammary Gland Development, Health, and Function of Cattle, Goats, and Sheep

Milk is an integral and therefore complex structural element of mammalian nutrition. Therefore, it is simple to conclude that lactation, the process of producing milk, is as complex as the mammary gland, the organ responsible for this biochemical activity. Nutrition, genetics, epigenetics, disease pathogens, climatic conditions, and other environmental variables all impact breast productivity. In the last decade, the number of studies devoted to epigenetics has increased dramatically. Reports are increasingly describing the direct participation of microRNAs (miRNAs), small noncoding RNAs that regulate gene expression post-transcriptionally, in the regulation of mammary gland development and function. This paper presents a summary of the current state of knowledge about the roles of miRNAs in mammary gland development, health, and functions, particularly during lactation. The significance of miRNAs in signaling pathways, cellular proliferation, and the lipid metabolism in agricultural ruminants, which are crucial in light of their role in the nutrition of humans as consumers of dairy products, is discussed.

The Effects of Heat Shock Protein 70 Addition in the Culture Medium on the Development and Quality of In Vitro Produced Heat Shocked Bovine Embryos

The aims of the present study were to examine the effects of HSP70 addition in the in vitro culture medium of day 3 embryos on their developmental competence and quality. Bovine oocytes (n = 1442) were in vitro matured, inseminated and cultured for the first two days according to standardized methods. The presumptive zygotes were randomly allocated in three experimental groups: Control, C (embryos cultured at 39 °C throughout the culture period), group C41 (temperature was raised to 41 °C from the 48th to 72nd h post insemination (p.i.) and then it returned at 39 °C for the remaining culture period), and group H41 (the temperature modification was the same as in C41 and during heat exposure, HSP70 was added in the culture medium). Cleavage and embryo yield were assessed 48 h p.i. and on days 7, 8, 9, respectively and gene expression in day 7 blastocysts was assessed by RT-PCR. Blastocyst yield was the highest in group C39; and higher in group H41 compared to group C41. From the gene expression analyses, altered expression of 11 genes was detected among groups. The analysis of the orchestrated patterns of gene expression differed between groups. The results of this study confirm the devastating effects of heat stress on embryo development and provide evidence that HSP70 addition at the critical stages can partly counterbalance, without neutralizing, the negative effects of the heat insult on embryos, acting mainly through mechanisms related to energy deployment.

Heat stress modulates polymorphonuclear cell response in early pregnancy cows: I. interferon pathway and oxidative stress

One of the major causes of early pregnancy loss is heat stress. In ruminants, interferon tau (IFNT) is the embryo signal to the mother. Once the interferon signaling pathway is activated, it drives gene expression for interferon-stimulated genes (ISGs) and alters neutrophils responses. The aim of the present study was to evaluate interferon (IFN) pathway, ISGs and gene expression in polymorphonuclear leukocytes (PMN) and oxidative stress in dairy cows under heat stress. Pregnant cows had their estrous cycle synchronized and randomly assigned to a comfort or heat stress group. Blood samples were collected at artificial insemination (AI) and on Days 10, 14 and 18 following AI. Pregnant cows were pregnancy checked by ultrasound on Day 30 and confirmed on Day 60 post-AI. Results are presented as mean ± SEM. The corpus luteum (CL) diameter was not different between groups of pregnant cows; concentration of progesterone of pregnant cows on Day 18 following AI was greater in comfort group compared to heat stressed group. Comfort pregnant cows had higher expression of all analyzed genes from interferon pathway, except for IFNAR1, on both Days 14 and 18. Conversely, heat stressed cows did not show altered expression of IFNT pathway genes and ISGs between Days 10, 14, and 18 after AI. The oxidative stress, determined as malondialdehyde (MDA) levels, was greater in heat stress group on Days 10, 14 and 18, independent of pregnancy status. Heat stress negatively influences expression of ISGs, IFN pathway gene expression in neutrophils, and oxidative stress. Our data suggest that lower conception rates in cows under heat stress are multifactorial, with the association of interferon pathway activation and the unbalanced oxidative stress being main contributing factors.

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.

Impacts of heat stress on immune responses and oxidative stress in farm animals and nutritional strategies for amelioration

Heat stress is one of the greatest challenges for the global livestock industries as increased environmental temperature and humidity compromises animal production during summer leading to devastating economic consequences. Over the last 30 years, significant developments have been achieved in cooling and provision of shade and shelter to mitigate heat stress reducing some of the losses associated with heat stress in farm animals. However, the recent increase in the incidence of heat waves which are also becoming more severe and lasting longer, due to climate change, further accentuates the problem of heat stress. Economic losses associated with heat stress are both direct due to loss in production and animal life, and indirect due to poorer quality products as a result of poor animal health and welfare. Animal health is affected due to impaired immune responses and increased reactive oxygen species production and/or deficiency of antioxidants during heat stress leading to an imbalance between oxidant and antioxidants and resultant oxidative stress. Research over the last 20 years has achieved partial success in understanding the intricacies of heat stress impacts on oxidative stress and immune responses and developing interventions to ameliorate impacts of heat stress, improving immune responses and farm animal health. This paper reviews the body of knowledge on heat stress impacts on immune response in farm animals. The impacts of heat stress on both cell-mediated and humoral immune responses have been discussed identifying the shift in immune response from cell-mediated towards humoral response, thereby weakening the immune status of the animal. Both species and breed differences have been identified as influencing how heat stress impacts the immune status of farm animals. In addition, crosstalk signaling between the immune system and oxidative stress has been considered and the role of antioxidants as potential nutritional strategies to mitigate heat stress has been discussed.

Dynamic changes in blood immune cell composition and function in Holstein and Jersey steers in response to heat stress

Heat stress has detrimental effects on livestock via diverse immune and physiological changes; heat-stressed animals are rendered susceptible to diverse diseases. However, there is relatively little information available regarding the altered immune responses of domestic animals in heat stress environments, particularly in cattle steers. This study aimed to determine the changes in the immune responses of Holstein and Jersey steers under heat stress. We assessed blood immune cells and their functions in the steers of two breeds under normal and heat stress conditions and found that immune cell proportions and functions were altered in response to different environmental conditions. Heat stress notably reduced the proportions of CD21+MHCII+ B cell populations in both breeds. We also observed breed-specific differences. Under heat stress, in Holstein steers, the expression of myeloperoxidase was reduced in the polymorphonuclear cells, whereas heat stress reduced the WC1+ γδ T cell populations in Jersey steers. Breed-specific changes were also detected based on gene expression. In response to heat stress, the expression of IL-10 and IL-17A increased in Holstein steers alone, whereas that of IL-6 increased in Jersey steers. Moreover, the mRNA expression pattern of heat shock protein genes such as Hsp70 and Hsp90 was significantly increased in only Holstein steers. Collectively, these results indicate that altered blood immunological profiles may provide a potential explanation for the enhanced susceptibility of heat-stressed steers to disease. The findings of this study provide important information that will contribute to developing new strategies to alleviate the detrimental effects of heat stress on steers.

Identification of key Genes and Pathways Associated With Thermal Stress in Peripheral Blood Mononuclear Cells of Holstein Dairy Cattle

The objectives of the present study were to identify key genes and biological pathways associated with thermal stress in Chinese Holstein dairy cattle. Hence, we constructed a cell-model, applied various molecular biology experimental techniques and bioinformatics analysis. A total of 55 candidate genes were screened from published literature and the IPA database to examine its regulation under cold (25°C) or heat (42°C) stress in PBMCs. We identified 29 (3 up-regulated and 26 down-regulated) and 41 (15 up-regulated and 26 down-regulated) significantly differentially expressed genes (DEGs) (fold change ≥ 1.2-fold and P < 0.05) after cold and heat stress treatments, respectively. Furthermore, bioinformatics analyses confirmed that major biological processes and pathways associated with thermal stress include protein folding and refolding, protein phosphorylation, transcription factor binding, immune effector process, negative regulation of cell proliferation, autophagy, apoptosis, protein processing in endoplasmic reticulum, estrogen signaling pathway, pathways related to cancer, PI3K- Akt signaling pathway, and MAPK signaling pathway. Based on validation at the cellular and individual levels, the mRNA expression of the HIF1A gene showed upregulation during cold stress and the EIF2A, HSPA1A, HSP90AA1, and HSF1 genes showed downregulation after heat exposure. The RT-qPCR and western blot results revealed that the HIF1A after cold stress and the EIF2A, HSPA1A, HSP90AA1, and HSF1 after heat stress had consistent trend changes at the cellular transcription and translation levels, suggesting as key genes associated with thermal stress response in Holstein dairy cattle. The cellular model established in this study with PBMCs provides a suitable platform to improve our understanding of thermal stress in dairy cattle. Moreover, this study provides an opportunity to develop simultaneously both high-yielding and thermotolerant Chinese Holstein cattle through marker-assisted selection.

Heat stress and immune response phenotype affect DNA methylation in blood mononuclear cells from Holstein dairy cows

Heat stress negatively affects health and production in cows. Examining the cellular response to heat stress could reveal underlying protective molecular mechanisms associated with superior resilience and ultimately enable selection for more resilient cattle. This type of investigation is increasingly important as future predictions for the patterns of heat waves point to increases in frequency, severity, and duration. Cows identified as high immune responders based on High Immune Response technology (HIR) have lower disease occurrence compared to their average and low immune responder herd-mates. In this study, our goal was to identify epigenetic differences between high and low immune responder cows in response to heat stress. We examined genome-wide DNA methylation of blood mononuclear cells (BMCs) isolated from high and low cows, before and after in vitro heat stress. We identified differential methylation of promoter regions associated with a variety of biological processes including immune function, stress response, apoptosis, and cell signalling. The specific differentially methylated promoter regions differed between samples from high and low cows, and results revealed pathways associated with cellular protection during heat stress.

The effects of heat stress on the behaviour of dairy cows – a review

Heat stress in livestock is a function of macro- and microclimatic factors, their duration and intensity, the environments where they occur and the biological characteristics of the animal. Due to intense metabolic processes, high-producing dairy cows are highly vulnerable to the effects of heat stress. Disturbances in their thermoregulatory capability are reflected by behavioural, physiological and production changes. Expression of thermoregulatory behaviour such as reduction of activity and feed intake, searching for a cooler places or disturbances in reproductive behaviours may be very important indicators of animal welfare. Especially maintain of standing or lying position in dairy cattle may be a valuable marker of the negative environmental impact. Highly mechanized farms with large numbers of animals have the informatic system can detect alterations automatically, while small family farms cannot afford these type of equipments. Therefore, observing and analysing behavioural changes to achieve a greater understanding of heat stress issue may be a key factor for developing the effective strategies to minimize the effects of heat stress in cattle. The aim of this review is to present the state of knowledge, over the last years, regarding behavioural changes in dairy cows (Bos Taurus) exposed to heat stress conditions and discuss some herd management strategies provided mitigation of the overheat consequences.

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.

Changes in Blood Metabolites and Immune Cells in Holstein and Jersey Dairy Cows by Heat Stress

Simple Summary

As global temperatures rise, thermal stress can be a major problem affecting cows. If they are subjected to heat stress, they are likely to exhibit abnormal metabolic reactions and affect their immune system. However, the relationship between metabolism and immunity during thermal stress and these crosstalk mechanisms remain unclear. Therefore, the aim of this study was to understand the changes in blood immune cell response with the physiological metabolism changes of Holstein and Jersey cows through the biochemistry and flow cytometry branches under thermal stress conditions. We found that various blood metabolites were reduced in both Holsteins and Jerseys by heat stress conditions. There were breed-specific variations in the immune cell population in Holstein and Jersey cows under different environmental conditions. The main findings of this study provide information on the metabolism and immunity changes of two types of cow under heat stress, broadening the potential relationship of these changes.

Abstract

Owing to increasing global temperatures, heat stress is a major problem affecting dairy cows, and abnormal metabolic responses during heat stress likely influence dairy cow immunity. However, the mechanism of this crosstalk between metabolism and immunity during heat stress remains unclear. We used two representative dairy cow breeds, Holstein and Jersey, with distinct heat-resistance characteristics. To understand metabolic and immune responses to seasonal changes, normal environmental and high-heat environmental conditions, we assessed blood metabolites and immune cell populations. In biochemistry analysis from sera, we found that variety blood metabolites were decreased in both Holstein and Jersey cows by heat stress. We assessed changes in immune cell populations in peripheral blood mononuclear cells (PBMCs) using flow cytometry. There were breed-specific differences in immune-cell population changes. Heat stress only increased the proportion of B cells (CD4–CD21+) and heat stress tended to decrease the proportion of monocytes (CD11b+CD172a+) in Holstein cows. Our findings expand the understanding of the common and specific changes in metabolism and immune response of two dairy cow breeds under heat stress conditions.

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

Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (F ST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.

Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia

The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.

Common and Differential Dynamics of the Function of Peripheral Blood Mononuclear Cells between Holstein and Jersey Cows in Heat-Stress Environment

Simple Summary

Seasonal change, particularly changing to hot and humid season, has a negative effect on dairy cows in various ways, including productivity, reproduction, metabolism, and immunity. In high-temperature and humid weather, dairy cows are vulnerable to diseases by weakened immune system. However, the cause of this has not been fully described. Therefore, this study aims to understand changes of specific gene expression and immune pathways based on transcriptome analysis from peripheral blood mononuclear cells of Holstein and Jersey dairy cows between normal and heat-stress environmental conditions. We observed that the two breeds of dairy cow have common and different immune shifts according to the changes of temperature and humidity condition. Overall, the findings of this study improve the understanding of the underlying mechanisms by which seasonal changes affect dairy cow immunity.

Abstract

Heat stress has been reported to affect the immunity of dairy cows. However, the mechanisms through which this occurs are not fully understood. Two breeds of dairy cow, Holstein and Jersey, have distinct characteristics, including productivity, heat resistance, and disease in high-temperature environments. The objective of this study is to understand the dynamics of the immune response of two breeds of dairy cow to environmental change. Ribonucleic acid sequencing (RNA-seq) results were analyzed to characterize the gene expression change of peripheral blood mononuclear cells (PBMCs) in Holstein and Jersey cows between moderate temperature-humidity index (THI) and high THI environmental conditions. Many of the differentially expressed genes (DEGs) identified are associated with critical immunological functions, particularly phagocytosis, chemokines, and cytokine response. Among the DEGs, CXCL3 and IL1A were the top down-regulated genes in both breeds of dairy cow, and many DEGs were related to antimicrobial immunity. Functional analysis revealed that cytokine and chemokine response-associated pathways in both Holstein and Jersey PBMCs were the most important pathways affected by the THI environmental condition. However, there were also breed-specific genes and pathways that altered according to THI environmental condition. Collectively, there were both common and breed-specific altered genes and pathways in Holstein and Jersey cows. The findings of this study expand our understanding of the dynamics of immunity in different breeds of dairy cow between moderate THI and high THI environmental conditions.

Inheritance patterns of leukocyte gene expression under heat stress in F1 hybrid cattle and their parents

Crossbreeding capitalizes on heterosis effects and results in increased performance of crossbred animals. Dominance hypothesis and overdominance hypothesis are 2 common models proposed to explain heterosis. Differential gene expression between parents and hybrids is hypothesized to be responsible for heterosis. This study aimed to investigate the heat tolerance and inheritance patterns of leukocyte transcriptomics in F₁ hybrid cattle (Angus males × Droughtmaster females) and their parents Red Angus (AN) and Droughtmaster (DR) under heat stress. According to the respiratory rate and heat tolerance coefficient index, DR was better adapted to heat stress than AN. The physiological responses to heat stress of F₁ hybrids were similar to AN. We identified 802 differentially expressed genes in leukocytes between AN and DR under heat stress using mRNA sequencing. Compared with AN, upregulated genes in DR were enriched in biological processes of response to stress, external and chemical stimulus, and cytokine, cell surface receptor signaling pathway, and cardiovascular system development. In contrast, upregulated genes in AN were enriched in B cell activation and regulation of B cell activation. Gene expression levels can be inherited additively or nonadditively and are classified into additive (35%), dominance (44%), and overdominance and underdominance (18%) modes in F₁ hybrids and their parents. Inheritance patterns of gene expression showed that 97% (249/255) of the dominant genes were classified as paternal AN dominant in hybrids. The paternal imprinted PEG10 gene and its regulatory transcription factor MYC showed an AN dominant expression pattern. The MYC interacted with most AN dominant genes. These transcriptomic analyses revealed that DR and AN had specific cellular and humoral immunity and cardiovascular systems development function under heat stress. Inheritance pattern analyses from gene expression partly explained phenotypic differences between parents and F₁ hybrids. The paternal imprinted PEG10 gene interaction with transcription factor MYC may contribute to explaining paternal dominant gene expression in hybrids.

Effects of chronic heat stress on mRNA and miRNA expressions in dairy cows

Heat stress has a negative impact on dairy cow productivity. In order to reveal the mechanisms of heat-stress response, the mRNA and miRNA expression profiles of five cows under chronic heat-stress and thermoneutral conditions were assayed in blood by high-throughput sequencing technology. A total of 540 mRNAs and 9 miRNAs were expressed differently under heat-stress and thermoneutral conditions. Functional analyses revealed that MAPK signaling pathway, cellular senescence, circadian entrainment, aldosterone synthesis and secretion, and pathways in cancer were enriched for differently expressed mRNAs; meanwhile cGMP-PKG signaling pathway, thermogenesis, and protein digestion and absorption were enriched for differently expressed miRNAs. In addition, GADD45G, TGFB2, and GNG11 may play roles in the heat stress, and bta-miR-423-5p might be one of the regulators of heat-stress response in cows as potential mediators of chronic heat-stress response. In conclusion, the present study described the mRNA and miRNA expression patterns in blood extracted from cows during the transition from heat-stress to thermoneutral conditions. The results provide new data that could help in identifying mechanisms that mediate cows' response to chronic heat stress.

Genome-Wide Transcriptome and Metabolome Analyses Provide Novel Insights and Suggest a Sex-Specific Response to Heat Stress in Pigs

Heat stress (HS) negatively impacts pig production and swine health. Therefore, to understand the genetic and metabolic responses of pigs to HS, we used RNA-Seq and high resolution magic angle spinning (HR-MAS) NMR analyses to compare the transcriptomes and metabolomes of Duroc pigs (n = 6, 3 barrows and 3 gilts) exposed to heat stress (33 °C and 60% RH) with a control group (25 °C and 60% RH). HS resulted in the differential expression of 552 (236 up, 316 down) and 879 (540 up, 339 down) genes and significant enrichment of 30 and 31 plasma metabolites in female and male pigs, respectively. Apoptosis, response to heat, Toll-like receptor signaling and oxidative stress were enriched among the up-regulated genes, while negative regulation of the immune response, ATP synthesis and the ribosomal pathway were enriched among down-regulated genes. Twelve and ten metabolic pathways were found to be enriched (among them, four metabolic pathways, including arginine and proline metabolism, and three metabolic pathways, including pantothenate and CoA biosynthesis), overlapping between the transcriptome and metabolome analyses in the female and male group respectively. The limited overlap between pathways enriched with differentially expressed genes and enriched plasma metabolites between the sexes suggests a sex-specific response to HS in pigs.

Genome-wide detection of signatures of selection in three Valdostana cattle populations

The Valdostana is a local dual purpose cattle breed developed in Italy. Three populations are recognized within this breed, based on coat colour, production level, morphology and temperament: Valdostana Red Pied (VPR), Valdostana Black Pied (VPN) and Valdostana Chestnut (VCA). Here, we investigated putative genomic regions under selection among these three populations using the Bovine 50K SNP array by combining three different statistical methods based either on allele frequencies (F_ST ) or extended haplotype homozygosity (iHS and Rsb). In total, 8, 5 and 8 chromosomes harbouring 13, 13 and 16 genomic regions potentially under selection were identified by at least two approaches in VPR, VPN and VCA, respectively. Most of these candidate regions were population-specific but we found one common genomic region spanning 2.38 Mb on BTA06 which either overlaps or is located close to runs of homozygosity islands detected in the three populations. This region included inter alia two well-known genes: KDR, a well-established coat colour gene, and CLOCK, which plays a central role in positive regulation of inflammatory response and in the regulation of the mammalian circadian rhythm. The other candidate regions identified harboured genes associated mainly with milk and meat traits as well as genes involved in immune response/inflammation or associated with behavioural traits. This last category of genes was mainly identified in VCA, which is selected for fighting ability. Overall, our results provide, for the first time, a glimpse into regions of the genome targeted by selection in Valdostana cattle. Finally, this study illustrates the relevance of using multiple complementary approaches to identify genomic regions putatively under selection in livestock.

Heat Stress Impacts Immune Status in Cows Across the Life Cycle

Heat stress has a myriad of effects on dairy cattle across the life cycle. Whereas, the most commonly recognized impacts are associated with production responses, emerging evidence indicates that heat stress profoundly alters the immune response of calves and cows, from the prenatal stage through lactation. For example, in utero heat stress reduces passive immune transfer regardless of colostrum source, relative to normothermic conditions in late gestation. Dry cows exposed to heat stress have lower immunoglobulin responses to ovalbumin vaccination, but this effect dissipates with cooling following parturition. Conversely, cows under heat stress when dry exhibit carryover effects on the innate arm of the immune system in early lactation. In this paper we review the effects of heat stress throughout the life cycle of the dairy cow, with particular emphasis on the impact of heat stress during late gestation on the cow and the developing fetus, both before and after parturition. In addition, the impact of altered immune status under heat stress on other physiological systems, especially those supporting milk production, are considered. Finally, management interventions to prevent and reverse the effect of heat stress are presented.

Clinical Ketosis-Associated Alteration of Gene Expression in Holstein Cows

Ketosis is one of the most prevalent transition metabolic disorders in dairy cows, and has been intrinsically influenced by both genetic and nutritional factors. However, altered gene expression with respective to dairy cow ketosis has not been addressed yet, especially at the genome-wide level. In this study, we recruited nine Holsteins diagnosed with clinical ketosis and ten healthy controls, for which whole blood samples were collected at both prepartum and postpartum. Four groups of blood samples were defined: from cows with ketosis at prepartum (PCK, N = 9) and postpartum (CK, N = 9), respectively, and controls at prepartum (PHC, N = 10) and postpartum (HC, N = 10). RNA-Seq approach was used for investigating gene expression, by which a total of 27,233 genes were quantified with four billion high-quality reads. Subsequently, we revealed 75 and four differentially expressed genes (DEGs) between sick and control cows at postpartum and prepartum, respectively, which indicated that sick and control cows had similar gene expression patterns at prepartum. Meanwhile, there were 95 DEGs between postpartum and prepartum for sick cows, which showed depressed changes of gene expression during this transition period in comparison with healthy cows (428 DEGs). Functional analyses revealed the associated DEGs with ketosis were mainly involved in biological stress response, ion homeostasis, AA metabolism, energy signaling, and disease related pathways. Finally, we proposed that the expression level of STX1A would be potentially used as a new biomarker because it was the only gene that was highly expressed in sick cows at both prepartum and postpartum. These results could significantly help us to understand the underlying molecular mechanisms for incidence and progression of ketosis in dairy cows.

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.

Cardiac and Skeletal Muscle Transcriptome Response to Heat Stress in Kenyan Chicken Ecotypes Adapted to Low and High Altitudes Reveal Differences in Thermal Tolerance and Stress Response

Heat stress (HS) negatively affects chicken performance. Agricultural expansion will happen in regions that experience high ambient temperatures, where fast-growing commercial chickens are vulnerable. Indigenous chickens of such regions, due to generations of exposure to environmental challenges, might have higher thermal tolerance. In this study, two indigenous chicken ecotypes, from the hot and humid Mombasa (lowland) and the colder Naivasha (highland) regions, were used to investigate the effects of acute (5 h, 35°C) and chronic (3 days of 35°C for 8 h/day) HS on the cardiac and skeletal muscle, through RNA sequencing. The rectal temperature gain and the number of differentially expressed genes (DEGs) [False Discovery Rate (FDR) < 0.05] were two times higher in the acute stage than in the chronic stage in both ecotypes, suggesting that cyclic exposure to HS can lead to adaptation. A tissue- and stage-specific difference in response to HS was observed, with peroxisome proliferator-activated-receptor (PPAR) signaling and mitogen-activate protein kinase (MAPK) signaling pathways, enriched in heart and skeletal muscle, respectively, and the p53 pathway enriched only in the acute stage in both tissues. The acute and chronic stage DEGs were integrated by a region-specific gene coexpression network (GCN), and genes with the highest number of connections (hub genes) were identified. The hub genes in the lowland network were CCNB2, Crb2, CHST9, SESN1, and NR4A3, while COMMD4, TTC32, H1F0, ACYP1, and RPS28 were the hub genes in the highland network. Pathway analysis of genes in the GCN showed that p53 and PPAR signaling pathways were enriched in both low and highland networks, while MAPK signaling and protein processing in endoplasmic reticulum were enriched only in the gene network of highland chickens. This shows that to dissipate the accumulated heat, to reduce heat induced apoptosis, and to promote DNA damage repair, the ecotypes activated or suppressed different genes, indicating the differences in thermal tolerance and HS response mechanisms between the ecotypes. This study provides information on the HS response of chickens, adapted to two different agro climatic environments, extending our understanding of the mechanisms of HS response and the effect of adaptation in counteracting HS.

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.

Serial gene co-expression network approach to mine biological meanings from integrated transcriptomes of the porcine endometrium during estrous cycle

The estrous cycle is a complex process regulated by several hormones. To understand the dynamic changes in gene expression that takes place in the swine endometrium during the estrous cycle relative to the day of estrus onset, we performed RNA-sequencing analysis on days 0, 3, 6, 9, 12, 15, and 18, resulting in the identification of 4495 differentially expressed genes (DEGs; Q ≤ 0.05 and |log₂FC| ≥ 1) at various phases in the estrous cycle. These DEGs were integrated into multiple gene co-expression networks based on different fold changes and correlation coefficient (R²) thresholds and a suitable network, which included 899 genes (|log2FC| ≥ 2 and R² ≥ 0.99), was identified for downstream analyses based on the biological relevance of the Gene Ontology (GO) terms enriched. The genes in this network were partitioned into 6 clusters based on the expression pattern. Several GO terms including cell cycle, apoptosis, hormone signaling, and lipid biosynthetic process were found to be enriched. Furthermore, we found 15 significant KEGG pathways, including cell adhesion molecules, cytokine-cytokine receptor signaling, steroid biosynthesis, and estrogen signaling pathways. We identified several genes and GO terms to be stage-specific. Moreover, the identified genes and pathways extend our understanding of porcine endometrial regulation during estrous cycle and will serve as a good resource for future studies.

Invited review: Advances and challenges in application of feedomics to improve dairy cow production and health

Dairy cattle science has evolved greatly over the past century, contributing significantly to the improvement in milk production achieved today. However, a new approach is needed to meet the increasing demand for milk production and address the increased concerns about animal health and welfare. It is now easy to collect and access large and complex data sets consisting of molecular, physiological, and metabolic data as well as animal-level data (such as behavior). This provides new opportunities to better understand the mechanisms regulating cow performance. The recently proposed concept of feedomics could help achieve this goal by increasing our understanding of interactions between the different components or levels and their impact on animal production. Feedomics is an emerging field that integrates a range of omics technologies (e.g., genomics, epigenomics, transcriptomics, proteomics, metabolomics, metagenomics, and metatranscriptomics) to provide these insights. In this way, we can identify the best strategies to improve overall animal productivity, product quality, welfare, and health. This approach can help research communities elucidate the complex interactions among nutrition, environment, management, animal genetics, metabolism, physiology, and the symbiotic microbiota. In this review, we summarize the outcomes of the most recent research on omics in dairy cows and highlight how an integrated feedomics approach could be applied in the future to improve dairy cow production and health. Specifically, we focus on 2 topics: (1) improving milk yield and milk quality, and (2) understanding metabolic physiology in transition dairy cows, which are 2 important challenges faced by the dairy industry worldwide.