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

The HSP90AA1 gene is involved in heat stress responses and its functional genetic polymorphisms are associated with heat tolerance in Holstein cows

As the stress-inducible isoform of the heat-shock protein 90 (HSP90), the HSP90AA1 gene encodes HSP90α and plays an important role in heat stress (HS) response. Therefore, this study aimed to investigate the role of the HSP90AA1 gene in cellular responses during HS and to identify functional SNPs associated with thermotolerance in Holstein cattle. For the in vitro validation experiment of acute HS, cells from the Madin-Darby bovine kidney cell line were exposed to 42°C for 1 h, and various parameters were assessed, including cell apoptosis, cell autophagy, and the cellular functions of HSP90α by using its inhibitor 17-allylamino-17-demethoxygeldanamycin (17-AAG). Furthermore, the polymorphisms identified in the HSP90AA1 gene and their functions related to HS were validated in vitro. Acute HS exposure induced cell apoptosis, cell autophagy, and upregulated expression of the HSP90AA1 gene. Inhibition of HSP90α by 17-AAG treatment had a significant effect on the expression of the HSP90α protein and increased cell apoptosis. However, autophagy decreased in comparison to the control treatment when cells were exposed to 42°C for 1 h. Five SNPs identified in the HSP90AA1 gene were significantly associated with rectal temperature and respiration score in Holstein cows, in which the rs109256957 SNP is located in the 3' untranslated region (3' UTR). Furthermore, we demonstrated that the 3' UTR of HSP90AA1 is a direct target of bta-miR-1224 by cell transfection with exogenous microRNA (miRNA) mimic and inhibitor. The luciferase assays revealed that the SNP rs109256957 affects the regulation of bta-miR-1224 binding activity and alters the expression of the HSP90AA1 gene. Heat stress-induced HSP90AA1 expression maintains cell survival by inhibiting cell apoptosis and increasing cell autophagy. The rs109256957 located in the 3' UTR region is a functional variation and it affects the HSP90AA1 expression by altering its binding activity with bta-miR-1224, thereby associating with the physiological parameters of Holstein cows.

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.

Whole genome discovery of regulatory genes responsible for the response of chicken to heat stress

Long noncoding RNAs (lncRNAs) are functional bridges connecting the genome with phenotypes by interacting with DNA, mRNA, and proteins. Using publically available acute heat stress (AHS)-related RNA-seq data, we discovered novel lncRNAs and tested their association with AHS along with ~ 8800 known lncRNAs and ~ 28,000 mRNA transcripts. Our pipeline discovered a total of 145 potentially novel-lncRNAs. One of them (Fishcomb_p-value = 0.06) along with another novel transcript (annotated as protein-coding; Fishcomb_p-value = 0.03) were identified as significantly associated with AHS. We found five known-lncRNAs and 134 mRNAs transcripts that were significantly associated with AHS. Four novel lncRNAs interact cis-regulated with 12 mRNA transcripts and are targeted by 11 miRNAs. Also six meta-lncRNAs associate with 134 meta-mRNAs through trans-acting co-expression, each targeted by 15 and 216 miRNAs, respectively. Three of the known-lncRNAs significantly co-expressed with almost 97 of the significant mRNAs (Pearson correlation p-value < 0.05). We report the mentioned three known-lncRNAs (ENSGALT00000099876, ENSGALT00000107573, and ENSGALT00000106323) as the most, significantly regulatory elements of AHS in chicken. It can be concluded that in order to alleviate the adverse effects of AHS on chicken, the manipulation of the three regulatory lncRNAs could lead to a more desirable result than the manipulation of the most significant mRNAs.

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.

Short heat shock has a long-term effect on mesenchymal stem cells' transcriptome

The adverse effects of heat stress (HS) on physiological systems are well documented, yet the underlying molecular mechanisms behind it remain poorly understood. To address this knowledge gap, we conducted a comprehensive investigation into the impact of HS on mesenchymal stem cells (MSCs), focusing on their morphology, phenotype, proliferative capacity, and fate determination. Our in-depth analysis of the MSCs' transcriptome revealed a significant influence of HS on the transcriptional landscape. Notably, even after a short period of stress, we observed a persistent alteration in cell identity, potentially mediated by the activation of bivalent genes. Furthermore, by comparing the differentially expressed genes following short HS with their transcriptional state after recovery, we identified the transient upregulation of MLL and other histone modifiers, providing a potential mechanistic explanation for the stable activation of bivalent genes. This could be used to predict and modify the long-term effect of HS on cell identity.

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).

Extraction and Elevation of Cell-Free DNA under Mastitis and Heat Stress in Dairy Cattle

In this study, four methods (phenol-chloroform protocol, sodium iodide kit, QIAamp DNA Blood Mini Kit, and TIANamp Micro DNA Kit) were used to extract cell-free DNA (cfDNA) from cattle blood, and the yield and purity of cfDNA varied in four different methods from 0.36 to 0.84 ng/mL for yield and 0.67 to 1.80 (A260/A280) for purity. Compared with other methods, the TIANamp Micro DNA kit performed better in both cfDNA amount and purity (p < 0.05); furthermore, blood cfDNA levels were significantly increased in Holstein dairy cows under the influence of heat stress (p < 0.01) and mastitis (p < 0.0001), which showed a potential power to discriminate mastitis (AUC = 0.99, 95% CI = 0.97 to 1.00) or heat stress (AUC = 0.86, 95% CI = 0.73 to 0.98) in cows. In brief, we established a complete experimental system for the extraction of cfDNA from cattle blood based on the high-yielding method of the TIANamp Micro DNA Kit and showed the effect of mastitis and heat stress on cfDNA levels in cattle blood for the first time. Our findings suggested that cfDNA in cattle blood may be a useful marker to measure mastitis and heat stress in dairy cattle.

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.

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.

Min pig skeletal muscle response to cold stress

The increased sensitivity of pigs to ambient temperature is due to today’s intensive farming. Frequent climate disasters increase the pressure on healthy pig farming. Min pigs are an indigenous pig breed in China with desirable cold resistance characteristics, and hence are ideal for obtaining cold-resistant pig breeds. Therefore, it is important to discover the molecular mechanisms that are activated in response to cold stress in the Min pig. Here, we conducted a transcriptomic analysis of the skeletal muscle of Min pigs under chronic low-temperature acclimation (group A) and acute short cold stress (group B). Cold exposure caused more genes to be upregulated. Totals of 125 and 96 differentially expressed genes (DEGs) were generated from groups A and B. Sixteen common upregulated DEGs were screened; these were concentrated in oxidative stress (SRXN1, MAFF), immune and inflammatory responses (ITPKC, AREG, MMP25, FOSL1), the nervous system (RETREG1, GADD45A, RCAN1), lipid metabolism (LRP11, LIPG, ITGA5, AMPD2), solute transport (SLC19A2, SLC28A1, SLCO4A1), and fertility (HBEGF). There were 102 and 73 genes that were specifically differentially expressed in groups A and B, respectively. The altered mRNAs were enriched in immune, endocrine, and cancer pathways. There were 186 and 91 differentially expressed lncRNAs generated from groups A and B. Analysis of the target genes suggested that they may be involved in regulating the MAPK signaling pathway for resistance to cold. The results of this study provide a comprehensive overview of cold exposure–induced transcriptional patterns in skeletal muscle of the Min pig. These results can guide future molecular studies of cold stress response in pigs for improving cold tolerance as a goal in breeding programs.

Transcriptome Analyses Reveal Essential Roles of Alternative Splicing Regulation in Heat-Stressed Holstein Cows

Heat stress (HS) severely impacts the productivity and welfare of dairy cows. Investigating the biological mechanisms underlying HS response is crucial for developing effective mitigation and breeding strategies. Therefore, we evaluated the changes in milk yield, physiological indicators, blood biochemical parameters, and alternative splicing (AS) patterns of lactating Holstein cows during thermoneutral (TN, N = 19) and heat-stress (HS, N = 17) conditions. There was a significant (p < 0.05) decline in milk yield as physiological indicators increased after exposure to natural HS conditions. The levels of eight out of 13 biochemical parameters of HS were also significantly altered in the presence of HS (p < 0.05). These results demonstrate that HS negatively influences various biological processes of Holstein cows. Furthermore, we investigated AS events based on the RNA-seq data from blood samples. With HS, five common types of AS events were generally increased by 6.7−38.9%. A total of 3470 AS events corresponding to 3143 unique genes were differentially alternatively spliced (DSGs) (p-adjusted < 0.05) between TN and HS groups. The functional annotation results show that the majority of DSGs are involved in mRNA splicing and spliceosomal complex, followed by enrichment in immune and metabolic processes. Eighty-seven out of 645 differentially expressed genes (DEGs) (fold change ≥ 1.5 and false discovery rate < 0.05) overlapped with DSGs. Further analyses showed that 20 of these genes were significantly enriched for the RNA splicing, RNA binding, and RNA transport. Among them, two genes (RBM25 and LUC7L3) had strong interrelation and co-expression pattern with other genes and were identified as candidate genes potentially associated with HS responses in dairy cows. In summary, AS plays a crucial role in changing the transcriptome diversity of heat-stress-related genes in multiple biological pathways and provides a different regulation mechanism from DEGs.

Effect of acute heat shock on stress gene expression and DNA methylation in zebu (Bos indicus) and crossbred (Bos indicus × Bos taurus) dairy cattle

Environmental temperature is one of the major factors to affect health and productivity of dairy cattle. Gene expression networks within the cells and tissues coordinate stress response, metabolism, and milk production in dairy cattle. Epigenetic DNA methylations were found to mediate the effect of environment by regulating gene expression patterns. In the present study, we compared three Indian native zebu cattle, Bos indicus (Sahiwal, Tharparkar, and Hariana) and one crossbred Bos indicus × Bos taurus (Vrindavani) for stress gene expression and differences in the DNA methylation patterns. The results indicated acute heat shock to cultured PBMC affected their proliferation, stress gene expression, and DNA methylation. Interestingly, expressions of HSP70, HSP90, and STIP1 were found more pronounced in zebu cattle than the crossbred cattle. However, no significant changes were observed in global DNA methylation due to acute heat shock, even though variations were observed in the expression patterns of DNA methyltransferases (DNMT1, DNMT3a) and demethylases (TET1, TET2, and TET3) genes. The treatment 5-AzaC (5-azacitidine) that inhibit DNA methylation in proliferating PBMC caused significant increase in heat shock-induced HSP70 and STIP1 expression indicating that hypomethylation facilitated stress gene expression. Further targeted analysis DNA methylation in the promoter regions revealed no significant differences for HSP70, HSP90, and STIP1. However, there was a significant hypomethylation for BDNF in both zebu and crossbred cattle. Similarly, NR3C1 promoter region showed hypomethylation alone in crossbred cattle. Overall, the results indicated that tropically adapted zebu cattle had comparatively higher expression of stress genes than the crossbred cattle. Furthermore, DNA methylation may play a role in regulating expression of certain genes involved in stress response pathways.

A Post-GWAS Functional Analysis Confirming Effects of Three BTA13 Genes CACNB2, SLC39A12, and ZEB1 on Dairy Cattle Reproduction

In our previous GWAS of Chinese and Nordic dairy cattle, genes CACNB2, SLC39A12, and ZEB1 locating on BTA 13 were suggested as candidate genes for reproduction. In this study, validation of these associations was performed in an independent population with records of nine reproductive traits. More importantly, functions of these genes in the reproductive process were verified by employing the expression data of ovarian follicles. The potential variants within the three genes were firstly detected in 68 Chinese Holstein bulls, and then screened in 1,588 Chinese Holstein cows using the KASP (Kompetitive allele-specific PCR) method. There were nine variants with polymorphisms in CACNB2, five in SLC39A12, and four in ZEB1, respectively, of which one SNP was in the upstream regulatory region, two in exon region, four in downstream regulatory region, and 11 SNPs in intronic regions. Amongst the 18 variants, g.33267056T/G in CACNB2 explained the largest phenotypic variance for age at first calving (0.011%), interval from first to last insemination (0.004%), and calving ease (0.002%), while g.32751518G/A in SLC39A12 contributed the most to stillbirth in heifers (0.038%). Two haplotype blocks were constructed for CACNB2 while one each for SLC39A12 and ZEB1, which were significantly associated with five reproductive traits, including age at the first service, age at the first calving, calving ease in heifers and cows, and the interval from calving to the first insemination. We then studied the profile of gene expression in granulosa cells isolated from four developmental stages of ovarian follicles from eight dairy cows. All three genes were differentially expressed between ovarian follicles with different sizes (p < 0.05), indicating their potential roles in the reproductive process of dairy cows. This study successfully demonstrated the associations of three BTA 13 genes CACNB2, SLC39A12, and ZEB1 with reproduction and further examined their expression levels in ovarian follicles directly. These findings can be beneficial for the ongoing genomic selection program for reproductive traits which have long been considered as traits that are difficult to achieve genetic improvement due to the lack of efficient genetic markers.

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

Simple Summary

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

Abstract

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