Changes in m6A RNA methylation of goat lung following PPRV infection

Peste des petits ruminants (PPR) is an acute, highly contagious viral disease of goats and sheep, caused by the Peste des petits ruminants virus (PPRV). Earlier studies suggest the involvement of diverse regulatory mechanisms in PPRV infection. Methylation at N6 of Adenosine called m6A is a type RNA modification that influences various physiological and pathological phenomena. As the lung tissue represents the primary target organ of PPRV, the present study explored the m6A changes and their functional significance in PPRV disease pathogenesis. m6A-seq analysis revealed 1289 m6A peaks to be significantly altered in PPRV infected lung in comparison to normal lung, out of which 975 m6A peaks were hypomethylated and 314 peaks were hypermethylated. Importantly, hypomethylated genes were enriched in Interleukin-4 and Interleukin-13 signaling and various processes associated with extracellular matrix organization. Further, of the 843 differentially m6A-containing cellular transcripts, 282 transcripts were also found to be differentially expressed. Functional analysis revealed that these 282 transcripts are significantly enriched in signaling by Interleukins, extracellular matrix organization, cytokine signaling in the immune system, signaling by receptor tyrosine kinases, and Toll-like Receptor Cascades. We also found m6A reader HNRNPC and the core component of methyltransferase complex METTL14 to be highly upregulated than the m6A readers - HNRNPA2B1 and YTHDF1 at the transcriptome level. These findings suggest that alteration in the m6A landscape following PPRV is implicated in diverse processes including Interleukin signaling.

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

Protective Efficacy of Multiple Epitope-Based Vaccine against Hyalomma anatolicum, Vector of Theileria annulata and Crimean-Congo Hemorrhagic Fever Virus

Hyalomma anatolicum is the principal vector for Theileria annulata, T. equi, and T. Lestoquardi in animals and the Crimean-Congo hemorrhagic fever virus in humans. Due to the gradual loss of efficacy of the available acaricides against field tick populations, the development of phytoacaricides and vaccines has been considered the two most critical components of the integrated tick management strategies. In the present study, in order to induce both cellular and humoral immune responses in the host against H. anatolicum, two multi-epitopic peptides (MEPs), i.e., VT1 and VT2, were designed. The immune-stimulating potential of the constructs was determined by in silicoinvestigation on allergenicity (non-allergen, antigenic (0.46 and 1.0046)), physicochemical properties (instability index 27.18 and 35.46), as well as the interaction of constructs with TLRs by docking and molecular dynamics analysis. The immunization efficacy of the MEPs mixed with 8% Montanide^TM gel 01 PR against H. anatolicum larvae was determined as 93.3% and 96.9% in VT1- and VT2-immunized rabbits, respectively. Against adults, the efficacy was 89.9% and 86.4% in VT1- and VT2-immunized rabbits, respectively. A significant (p < 0.001) reduction in the anti-inflammatory cytokine (IL-4) and significantly higher IgG response was observed in a VT1-immunized group of rabbits as compared with the response observed in the control group. However, in the case of the VT2-immunized rabbits, an elevated anti-VT2 IgG and pro-inflammatory cytokine (IL-2) (>30 fold) along with a decreased level of anti-inflammatory cytokine IL-4 (0.75 times) was noted. The efficacy of MEP and its potential immune stimulatory responses indicate that it might be useful for tick management.

N6-Methyladenosine RNA Modification in Host Cells Regulates Peste des Petits Ruminants Virus Replication

N6-methyladenosine (m6A) modification is a major RNA epigenetic regulatory mechanism. The dynamics of m6A levels in viral genomic RNA and their mRNAs have been shown to have either pro- or antiviral functions, and therefore, m6A modifications influence virus-host interactions. Currently, no reports are available on the effect of m6A modifications in the genome of Peste des petits ruminants virus (PPRV). In the present study, we took PPRV as a model for nonsegmented negative-sense single-stranded RNA viruses and elucidate the role of m6A modification on viral replication. We detected m6A-modified sites in the mRNA of the virus and host cells, as well as the PPRV RNA genome. Further, it was found that the level of m6A modification in host cells alters the viral gene expression. Knockdown of the METTL3 and FTO genes (encoding the m6A RNA modification writer and eraser proteins, respectively) results in alterations of the levels of m6A RNA modifications in the host cells. Experiments using these genetically modified clones of host cells infected with PPRV revealed that both higher and lower m6A RNA modification in the host cells negatively affect PPRV replication. We found that m6A-modified viral transcripts had better stability and translation efficiency compared to the unmodified mRNA. Altogether, from these data, we conclude that the m6A modification of RNA regulates PPRV replication. These findings contribute toward a way forward for developing novel antiviral strategies against PPRV by modulating the dynamics of host m6A RNA modification. IMPORTANCE Peste des petits ruminants virus (PPRV) causes a severe disease in sheep and goats. PPRV infection is a major problem, causing significant economic losses to small ruminant farmers in regions of endemicity. N6-methyladenosine (m6A) is an important RNA modification involved in various functions, including virus-host interactions. In the present study, we used stable clones of Vero cells, having knocked down the genes encoding proteins involved in dynamic changes of the levels of m6A modification. We also used small-molecule compounds that interfere with m6A methylation. This resulted in a platform of host cells with various degrees of m6A RNA modification. The host cells with these different microenvironments were useful for studying the effect of m6A RNA modification on the expression of viral genes and viral replication. The results pinpoint the level of m6A modifications that facilitate the maximum replication of PPRV. These findings will be useful in increasing the virus titers in cultured cells needed for the economical development of the vaccine. Furthermore, the findings have guiding significance for the development of novel antiviral strategies for limiting PPRV replication in infected animals.

Differential expression of long non-coding RNAs under Peste des petits ruminants virus (PPRV) infection in goats

Peste des petits ruminants (PPR) characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia, is an acute, highly contagious viral disease of sheep and goats. The role of long non-coding RNAs (lncRNAs) in PPRV infection has not been explored to date. In this study, the transcriptome profiles of virulent Peste des petits ruminants virus (PPRV) infected goat tissues – lung and spleen were analyzed to identify the role of lncRNAs in PPRV infection. A total of 13,928 lncRNA transcripts were identified, out of which 170 were known lncRNAs. Intergenic lncRNAs (7625) formed the major chunk of the novel lncRNA transcripts. Differential expression analysis revealed that 15 lncRNAs (11 downregulated and 4 upregulated) in the PPRV infected spleen samples and 16 lncRNAs (13 downregulated and 3 upregulated) in PPRV infected lung samples were differentially expressed as compared to control. The differentially expressed lncRNAs (DElncRNAs) possibly regulate various immunological processes related to natural killer cell activation, antigen processing and presentation, and B cell activity, by regulating the expression of mRNAs through the cis- or trans-regulatory mechanism. Functional enrichment analysis of differentially expressed mRNAs (DEmRNAs) revealed enrichment of immune pathways and biological processes in concordance with the pathways in which correlated lncRNA-neighboring genes were enriched. The results suggest that a coordinated immune response is raised in both lung and spleen tissues of the goat through mRNA-lncRNA crosstalk.

Canine Parvovirus and Its Non-Structural Gene 1 as Oncolytic Agents: Mechanism of Action and Induction of Anti-Tumor Immune Response

The exploration into the strategies for the prevention and treatment of cancer is far from complete. Apart from humans, cancer has gained considerable importance in animals because of increased awareness towards animal health and welfare. Current cancer treatment regimens are less specific towards tumor cells and end up harming normal healthy cells. Thus, a highly specific therapeutic strategy with minimal side effects is the need of the hour. Oncolytic viral gene therapy is one such specific approach to target cancer cells without affecting the normal cells of the body. Canine parvovirus (CPV) is an oncolytic virus that specifically targets and kills cancer cells by causing DNA damage, caspase activation, and mitochondrial damage. Non-structural gene 1 (NS1) of CPV, involved in viral DNA replication is a key mediator of cytotoxicity of CPV and can selectively cause tumor cell lysis. In this review, we discuss the oncolytic properties of Canine Parvovirus (CPV or CPV2), the structure of the NS1 protein, the mechanism of oncolytic action as well as role in inducing an antitumor immune response in different tumor models.

Identification of oxidant susceptible proteins in Salmonella Typhimurium

The human gut pathogen, Salmonella Typhimurium (S. Typhimurium) not only survives but also replicates inside the phagocytic cells. Bacterial proteins are the primary targets of phagocyte generated oxidants. Because of the different amino acid composition, some proteins are more prone to oxidation than others. Many oxidant induced modifications to amino acids have been described. Introduction of carbonyl group is one of such modifications, which takes place quite early following exposure of proteins to oxidants and is quite stable. Therefore, carbonyl groups can be exploited to identify oxidant susceptible proteins. Hypochlorous acid (HOCl) is one of the most potent oxidants produced by phagocytes. Incubation of S. Typhimurium with 3 mM HOCl resulted in more than 150 folds loss of bacterial viability. Proteins extracted from HOCl exposed S. Typhimurium cells showed about 60 folds (p < 0.001) more carbonyl levels as compared to unexposed cells. Similarly, 2, 4-Dinitrophenylhydrazine (2, 4-DNPH) derivatized proteins of HOCl treated S. Typhimurium cultures reacted strongly with anti-DNP antibodies as compared to buffer treated counterpart. Next, we have derivatized carbonyl groups on the proteins with biotin hydrazide. The derivatized proteins were then isolated by avidin affinity chromatography. Mass spectrometry based analysis revealed the presence of 204 proteins.

Cross-talk between energy metabolism and epigenetics during temperature stress response in C2C12 myoblasts

Objective: Environmental stress induces disturbances in cell energy metabolism and may cause epigenetic modifications. This study aimed to understand the possible impact of temperature stress (35 °C, 39 °C and 41 °C, compared to control 37 °C) on energy metabolism and epigenetic modifications, such as DNA methylation and histone H4 acetylation, as well as its effects on the expression of genes responsible for epigenetic changes, in mouse skeletal myoblasts (C2C12 cells). Methods: The results showed significantly reduced maximal respiration and spare respiratory capacity under heat stress (39 °C and 41 °C), suggesting that mitochondrial functions were compromised under these conditions. The glycolytic capacity and glycolysis markedly increased following low-temperature stress (35 °C). The results suggested that, under cold stress, cells prefer glycolysis as a rapid compensatory mechanism to meet energy requirements for adaptive thermogenic response. Results: Epigenetic changes (histone H4 acetylation and global DNA methylation) were observed under both heat and cold stress. Among the genes coding for DNA methyltransferases, the Dnmt3a was significantly increased under high-temperature conditions (39 °C and 41 °C), while Dnmt1 expression was significantly increased at low temperature (35 °C), indicating that under these conditions the cells preferred maintenance of methylation to de novo methylation activity. An expression pattern similar to Dnmt3a was observed for Gcn5, encoding for a histone acetyltransferase. The study revealed that temperature stress induced changes in the metabolic profiles, as well as epigenetic modifications, including the dynamics of the key enzymes. Conclusion: The results indicated the existence of crosstalk mechanisms between energy metabolism and epigenetics during cell stress response.

DNA methylation analysis of porcine mammary epithelial cells reveals differentially methylated loci associated with immune response against Escherichia coli challenge

Background

Epigenetic changes such as cytosine (CpG) DNA methylations regulate gene expression patterns in response to environmental cues including infections. Microbial infections induce DNA methylations that play a potential role in modulating host-immune response. In the present study, we sought to determine DNA methylation changes induced by the mastitis causing Escherichia coli (E. coli) in porcine mammary epithelial cells (PMEC). Two time points (3 h and 24 h) were selected based on specific transcriptomic changes during the early and late immune responses, respectively.

Results

DNA methylation analysis revealed 561 and 898 significant (P < 0.01) differentially methylated CpG sites at 3 h and 24 h after E. coli challenge in PMEC respectively. These CpG sites mapped to genes that have functional roles in innate and adaptive immune responses. Significantly, hypomethylated CpG sites were found in the promoter regions of immune response genes such as SDF4, SRXN1, CSF1 and CXCL14. The quantitative transcript estimation indicated higher expression associated with the DNA CpG methylation observed in these immune response genes. Further, E. coli challenge significantly reduced the expression levels of DNMT3a, a subtype of de novo DNA methylation enzyme, in PMEC indicating the probable reason for the hypomethylation observed in the immune response genes.

Conclusions

Our study revealed E. coli infection induced DNA methylation loci in the porcine genome. The differentially methylated CpGs were identified in the regulatory regions of genes that play important role in immune response. These results will help to understand epigenetic mechanisms for immune regulation during coliform mastitis in pigs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5976-7) contains supplementary material, which is available to authorized users.

Contrasting Gene Expression Profiles of Monocytes and Lymphocytes From Peste-Des-Petits-Ruminants Virus Infected Goats

In this study, transcriptome analysis of PPRV infected PBMC subsets-T helper cells, T cytotoxic cells, monocytes, and B lymphocytes was done to delineate their role in host response. PPRV was found to infect lymphocytes and not monocytes. The established receptor for PPRV-SLAM was found downregulated in lymphocytes and non-differentially expressed in monocytes. A profound deviation in the global gene expression profile with a large number of unique upregulated genes (851) and downregulated genes (605) was observed in monocytes in comparison to lymphocytes. ISGs-ISG15, Mx1, Mx2, RSAD2, IFIT3, and IFIT5 that play a role in antiviral response and the genes for viral sensors-MDA5, LGP2, and RIG1, were found to be upregulated in lymphocytes and downregulated in monocytes. The transcription factors-IRF-7 and STAT-1 that regulate expression of most of the ISGs were found activated in lymphocytes and not in monocytes. Interferon signaling pathway and RIG1 like receptor signaling pathway were found activated in lymphocytes and not in monocytes. This contrast in gene expression profiles and signaling pathways indicated the predominant role of lymphocytes in generating the antiviral response against PPRV in goats, thus, giving us new insights into host response to PPRV.

Molecular epidemiological analysis of wild animal rabies isolates from India

Aim

This study was conducted to know the genetic variability of rabies viruses (RVs) from wild animals in India.

Materials and Methods

A total of 20 rabies suspected brain samples of wild animals from different states of India were included in the study. The samples were subjected for direct fluorescent antibody test (dFAT), reverse transcription polymerase chain reaction (RT-PCR), and quantitative reverse transcriptase real-time PCR (RT-qPCR). The phylogenetic analysis of partial nucleoprotein gene sequences was performed.

Results

Of 20 samples, 11, 10, and 12 cases were found positive by dFAT, RT-PCR, and RT-qPCR, respectively. Phylogenetic analysis showed that all Indian wild RVs isolates belonged to classical genotype 1 of Lyssavirus and were closely related to Arctic/Arctic-like single cluster indicating the possibility of a spillover of rabies among different species.

Conclusion

The results indicated the circulation of similar RVs in sylvatic and urban cycles in India. However, understanding the role of wild animals as reservoir host needs to be studied in India.

VARV B22R homologue as phylogenetic marker gene for Capripoxvirus classification and divergence time dating

Sheeppox disease is associated with significant losses in sheep production world over. The sheep pox virus, the goatpox virus, and the lumpy skin disease virus cannot be distinguished by conventional serological tests. Identification of these pathogens needs molecular methods. In this study, seven genes viz. EEV maturation protein-F12L, Virion protein-D3R, RNA polymerase subunit-A5R, Virion core protein-A10L, EEV glycoprotein-A33R, VARV B22R homologue, and Kelch like protein-A55R that cover the start, middle, and end of the genome were selected. These genes were amplified from Roumanian-Fanar vaccine strain and Jaipur virulent strain, cloned, and sequenced. On analysis with the available database sequences, VARV B22R homologue was identified as a marker for phylogenetic reconstruction for classifying the sheeppox viruses of the ungulates. Further, divergence time dating with VARV B22R gene accurately predicted the sheeppox disease outbreak involving Jaipur virulent strain.

Molecular epidemiology of rabies virus circulating in domestic animals in India

Rabies is a neglected viral zoonotic disease affecting humans, domestic and wild animals and is endemic in most parts of the India. Dog mediated rabies is more predominant than other forms of rabies and molecular epidemiology is poorly understood in both reservoir and susceptible hosts. In the present study, a total of 140 rabies suspected brain samples from different species of animals from different geographical regions of India were used. The samples were parallelly tested by direct fluorescent antibody test, reverse transcriptase PCR and real-time PCR. Thirty positive samples were subjected for partial nucleoprotein gene sequencing and phylogenetic analysis. On sequence and phylogenetic analysis, it was observed that all Indian rabies viruses belonged to classical rabies virus of genotype 1 of Lyssavirus and formed two distinct groups. The majority of isolates were in group-1 and are closely related to arctic/arctic like lineage, whereas group-II isolated are closely related to cosmopolitan lineage. These results indicated there is simultaneous existence of two distinct lineages of rabies viruses in Indian subcontinent. Further whole genome studies are needed for better understanding of molecular epidemiology of rabies virus circulating in animals for control and prevention of rabies in India.

Effect of nicotinic acetylcholine receptor alpha 1 (nAChRα1) peptides on rabies virus infection in neuronal cells

Rabies virus (RABV) is neurotropic and causes acute progressive encephalitis. Herein, we report the interaction of nAChRα1-subunit peptides with RABV and the effect of these peptides on RABV infection in cultured neuronal cells. Peptide sequences derived from torpedo, bovine, human and rats were synthesized and studied for their interactions with RABV using virus capture ELISA and peptide immunofluorescence. The results showed specific binding of the nAChRα1-subunit peptides to the RABV. In the virus adsorption assay, these peptides were found to inhibit the attachment of the RABV to the neuronal cells. The nAChRα1-subunit peptides inhibited the RABV infection and reduced viral gene expression in the cultured neuroblastoma (N2A) cells. Torpedo peptide sequence (T-32) had highest antiviral effect (IC50=14±3.01μM) compared to the other peptides studied. The results of the study indicated that nAChRα1-subunit peptides may act as receptor decoy molecules and inhibit the binding of virus to the native host cell receptors and hence may reduce viral infection.

Differential effect of thermal stress on HSP70 expression, nitric oxide production and cell proliferation among native and crossbred dairy cattle

In a tropical country like India, thermal stress is one of the major factors which significantly affects the productivity of dairy cattle. The present study was aimed to identify the effect of heat and cold stress on cell viability, mitogen stimulation indices, nitric oxide production and HSP70 expression in Sahiwal and Holstein crossbred (Frieswal) population in India. The results indicated that the Sahiwal breed can better withstand the effect of heat and cold stress significantly (P<0.05) when compared to the crossbred cattle due to the higher survivability of the Peripheral Blood Mononuclear Cells (PBMCs) and Phytohemagglutinin (PHA-P) mitogen based stimulation indices. The study also revealed the significant differences (P<0.05) in the level of nitric oxide (µM) production amongst the pre and post thermal stressed samples of Sahiwal and Frieswal crossbred samples. Further, the expression of HSP70 was significantly (P<0.05) higher in Sahiwal compared to Frieswal immediately after heat/cold shock to 6h of recovery as indirect ELISA analysis showed gradual rise in the Hsp70 protein concentration (ng/ml) immediately after heat and cold stress (0h) and reached the peak at 6h of recovery. Western blot and immune fluorescent assay results were also corroborated with the findings of indirect ELISA. In Sahiwal cattle the mRNA expression of HSP70 and its protein concentration were higher (P<0.05) during peak summer (44°C) and winter (10°C) as compared to Frieswal cattle. This investigation supports the earlier information on the higher adaptability of indigenous cattle breeds to hot and humid conditions compared to the crossbreds of temperate cattle breeds.

Understanding the mechanisms of ATPase beta family genes for cellular thermotolerance in crossbred bulls

Na+/K+-ATPase is an integral membrane protein composed of a large catalytic subunit (alpha), a smaller glycoprotein subunit (beta), and gamma subunit. The beta subunit is essential for ion recognition as well as maintenance of the membrane integrity. Present study was aimed to analyze the expression pattern of ATPase beta subunit genes (ATPase B1, ATPase B2, and ATPase B3) among the crossbred bulls under different ambient temperatures (20-44 °C). The present study was also aimed to look into the relationship of HSP70 with the ATPase beta family genes. Our results demonstrated that among beta family genes, transcript abundance of ATPase B1 and ATPase B2 is significantly (P < 0.05) higher during the thermal stress. Pearson correlation coefficient analysis revealed that the expression of ATPase Β1, ATPase B2, and ATPase B3 is highly correlated (P < 0.01) with HSP70, representing that the change in the expression pattern of these genes is positive and synergistic. These may provide a foundation for understanding the mechanisms of ATPase beta family genes for cellular thermotolerance in cattle.

Effect of heat stress on the expression profile of Hsp90 among Sahiwal (Bos indicus) and Frieswal (Bos indicus × Bos taurus) breed of cattle: a comparative study

We evaluated the effect of thermal challenge on the expression profile of heat shock protein 90 (Hsp90) among Sahiwal (Bos indicus) and Frieswal (Bos indicus × Bos taurus) breeds of cattle. The present investigation was focused on the comparative studies on Hsp90 expression among Frieswal and Sahiwal under in vitro and environmental heat stress. Measured immediately after the in vitro heat shock to the peripheral blood mononuclear cells (PBMCs), the relative expression of Hsp90 mRNA was significantly (P<0.05) higher in Sahiwal compared to those in Frieswal. In later intervals of time, the differences in the expression levels between the two breeds become negligible coming down towards the basal level. A similar pattern was observed in the protein concentration showing significantly (P<0.05) higher levels in Sahiwal compared to those in Frieswal. The second sets of experiments were undertaken during summer months (March to May) when temperature peaked from 37 to 45 °C. During these months, Frieswal cows consistently recorded higher rectal temperatures than the Sahiwal breed. Further during this peak summer stress, Sahiwal showed significantly higher levels of mRNA transcripts as well as protein concentration compared to the Frieswal breed. Our findings also interestingly showed that, the cell viability of PBMC are significantly higher among the Sahiwal than Frieswal. Taken together, the experiments of both induced in vitro and environmental stress conditions indicate that, Sahiwal may express higher levels of Hsp90 then Frieswal to regulate their body temperature and increase cell survivality under heat stressed conditions.

Promoter variants at AP2 box region of Hsp70.1 affect thermal stress response and milk production traits in Frieswal cross bred cattle

Heat shock proteins (Hsp) are known to play major role in protection of cells from thermal stress. Nucleotide polymorphisms within the promoter of Hsp affect degree of expression and inducibility of Hsp mRNA. The present study aimed to investigate the effect of polymorphism within promoter region on the cellular expression of Hsp70.1 mRNA and association of identified polymorphisms with the physiological parameters during summer stress and milk production traits in dairy cattle. Two hundred Frieswal cows were genotyped using double PCR-RFLP to identify deletion of cytosine within the Hsp70.1 promoter AP2 box at base position 895. Homozygous wild type genotypes (CC) were found in lower frequency (39.29, n=78) than heterozygous cytosine deletion mutant genotypes (C-) (60.71, n=122). In the observed physiological parameters (rectal temperature, respiration rate and heat tolerance coefficient), cows that were homozygous wild types had better significant (P<0.05) summer tolerance than the heterozygous deletion genotypes. Cytosine deletion mutation in the promoter region negatively affected (P<0.01) the expression of Hsp70.1 mRNA in peripheral bovine mononuclear cells (PBMC) subjected to in vitro heat stress. Further association of observed polymorphism with the milk production traits was significant as the heterozygous cytosine deletion cows had lower total milk yield, peak yield, yield at 300 days, protein% (P<0.01) and fat% (P<0.05) than the native wild type promoter cows. The results from the present study suggest that the promoter region of bovine hsp70.1 gene is polymorphic and may be useful in selection of dairy cows for relatively better thermotolerance and higher milk production.