Circulating leukocytes can be used as an effective model to understand the heat stress response of different cattle types and buffaloes. This investigation aimed to determine the temporal profile of HSPs (HSP40, HSP60, HSP70, and HSP90) expression in circulating peripheral blood mononuclear cells (PBMCs) of Murrah buffaloes, Holstein-Friesian (HF), and Sahiwal cows in response to sublethal heat shock at 42°C. The viability data indicated HF PBMCs to be the most affected to the heat shock, whereas Sahiwal PBMCs were least affected, indicating its better survivability during the heat stress condition. The qRT-PCR expression data showed significant increase in mRNA expression of the analyzed HSPs genes after heat stimuli to the PBMCs under in vitro condition. In each case, the HSPs were most upregulated at 2 h after the heat stress. Among the HSPs, HSP70 was relatively more expressed followed by HSP60 indicating the action of molecular chaperones to stabilize the native conformation of proteins. However, PBMCs from different cattle types and buffaloes showed difference in the extent of transcriptional response. The level of expression of HSPs throughout the time period of heat stress was highest in buffaloes, followed by HF and Sahiwal cows. The higher abundance of HSP70 mRNA at each time point after heat stress showed prolonged effect of heat stress in HF PBMCs. The data presented here provided initial evidence of transcriptional differences in PBMCs of different cattle types and buffaloes and warrant further research.

The present study aims to identify the heat responsive genes and biological pathways in heat stressed buffalo mammary epithelial cells (MECs). The primary mammary epithelial cells of riverine buffalo were exposed to thermal stress at 42°C for one hour. The cells were subsequently allowed to recover at 37°C and harvested at different time intervals (30 min to 48 h) along with control samples (un-stressed). In order to assess the impact of heat stress in buffalo MECs, several in-vitro cellular parameters (lactate dehydrogenase activity, cell proliferation assay, cellular viability, cell death and apoptosis) and transcriptional studies were conducted. The heat stress resulted in overall decrease in cell viability and cell proliferation of MECs while induction of cellular apoptosis and necrosis. The transcriptomic profile of heat stressed MECs was generated using Agilent 44 K bovine oligonucleotide array and at cutoff criteria of ≥3-or ≤3 fold change, a total of 153 genes were observed to be upregulated while 8 genes were down regulated across all time points post heat stress. The genes that were specifically up-regulated or down-regulated were identified as heat responsive genes. The upregulated genes in heat stressed MECs belonged to heat shock family viz., HSPA6, HSPB8, DNAJB2, HSPA1A. Along with HSPs, genes like BOLA, MRPL55, PFKFB3, PSMC2, ENDODD1, ARID5A, and SENP3 were also upregulated. Microarray data revealed that the heat responsive genes belonged to different functional classes viz., chaperons; immune responsive; cell proliferation and metabolism related. Gene ontology analysis revealed enrichment of several biological processes like; cellular process, metabolic process, response to stimulus, biological regulation, immune system processes and signaling. The transcriptome analysis data was further validated by RT-qPCR studies. Several HSP (HSP40, HSP60, HSP70, HSP90, and HSPB1), apoptotic (Bax and Bcl2), immune (IL6, TNFα and NF-kβ) and oxidative stress (GPX1 and DUSP1) related genes showed differential expression profile at different time points post heat stress. The transcriptional data strongly indicated the induction of survival/apoptotic mechanism in heat stressed buffalo MECs. The overrepresented pathways across all time points were; electron transport chain, cytochrome P450, apoptosis, MAPK, FAS and stress induction of HSP regulation, delta Notch signaling, apoptosis modulation by HSP70, EGFR1 signaling, cytokines and inflammatory response, oxidative stress, TNF-alpha and NF- kB signaling pathway. The study thus identified several genes from different functional classes and biological pathways that could be termed as heat responsive in buffalo MEC. The responsiveness of buffalo MECs to heat stress in the present study clearly suggested its suitability as a model to understand the modulation of buffalo mammary gland expression signature in response to environmental heat load.

Gene expression studies require appropriate normalization methods for proper evaluation of reference genes. To date, not many studies have been reported on the identification of suitable reference genes in buffaloes. The present study was undertaken to determine the panel of suitable reference genes in heat-stressed buffalo mammary epithelial cells (MECs). Briefly, MEC culture from buffalo mammary gland was exposed to 42 °C for one hour and subsequently allowed to recover at 37 °C for different time intervals (from 30 m to 48 h). Three different algorithms, geNorm, NormFinder, and BestKeeper softwares, were used to evaluate the stability of 16 potential reference genes from different functional classes. Our data identified RPL4, EEF1A1, and RPS23 genes to be the most appropriate reference genes that could be utilized for normalization of qPCR data in heat-stressed buffalo MECs.

No abstract