Molecular Characterization and Seasonal Variation in Expression of HSP70.1 Gene in Gangatiri Cattle and Its Comparison with Buffalo

Under tropical climate heat stress is a major challenge for livestock production. HSP70.1 is a ubiquitously expressed protein maintaining cellular machinery through proper folding of denatured proteins and prevents cellular apoptosis and protect cell from heat stress. Therefore, present investigation was undertaken to explore genetic variability in HSP70.1 gene in Gangatiri cattle, its comparison with buffalo sequences and differential expression in different season. The allelic variant was identified by sequencing amplified PCR product of HSP70.1 gene by primer walking. Season-wise total RNA samples was prepared for differential expression study. Brilliant SYBR Green QPCR technique was used to study the expression kinetics of this gene. DNA sequencing by primer walking identified four allelic variants in Gangatiri cattle. Sequence alignment study revealed four, six and one substitutions in the 5' untranslated region (5'UTR), coding and 3' untranslated region ((3'UTR) of HSP70.1 gene, respectively. Comparative analysis of HSP70.1 gene revealed that Cattle has shorter 5'UTR and 3' UTR than the buffalo. In Gangatiri cattle, summer season has significantly higher (P ≤ 0.05) expression of HSP70.1 than the spring and winter. The relative expression of HSP70.1 was increased by more than six folds in summer and nearly 1.5 folds higher in winter in comparison to the spring season. Therefore, HSP70.1 may be considered to have a critical role in the development of thermal tolerance in Gangatiri cattle.

Comparative miRNA signatures among Sahiwal and Frieswal cattle breeds during summer stress

MicroRNAs (miRNAs) are known to take part in different biological mechanisms, including biotic as well as abiotic cellular stresses. The present investigation was aimed to identify comparative expression profile of differentially expressed miRNAs among Sahiwal (Bos indicus) and Frieswal (Bos indicus × Bos taurus) cattle breeds during summer stress. Stress responses in animals were characterized by recording various physiological parameters, biochemical assays and expression profiling of heat shock protein 70 (Hsp70) during elevated environmental temperature. Ion Torrent-based deep sequencing as well as CLC-genomic analysis identified 322 and 420 Bos taurus annotated miRNAs among Sahiwal and Frieswal, respectively. A total 69 common miRNAs were identified to be differentially expressed during summer among the breeds. Out of the 69, a total 14 differentially expressed miRNAs viz. bta-mir 6536-2, bta-mir-2898, bta-mir-let-7b, bta-mir-425, bta-mir-2332, bta-mir-2478, bta-mir-150, bta-mir142, bta-mir-16a, bta-mir-2311, bta-mir-1839, bta-mir-1248-1, bta-mir-103-2 and bta-mir-181b were randomly selected for qRT-PCR-based validation. bta-mir-2898, bta-mir-6536-1, bta-mir-let-7b, bta-mir-2478, bta-mir-150, bta-mir-16a, bta-mir-2311, bta-mir-1032-b and bta-mir-181-b were significantly (p < 0.01) upregulated during summer among Frieswal in comparison to Sahiwal while, bta-mir 6536-2, bta-mir-2332, bta-mir142, bta-mir-1839 and bta-mir-1248-1 was significantly (p < 0.01) expressed at higher level in Sahiwal in contrast to Frieswal correlation coefficient analysis revealed that bta-mir(s)-150, 16a and 181b are negatively correlated (p < 0.05) with Hsp70 expression. Thus, this study identified that miRNA expression during summer stress can vary between the breeds which may reflect their differential post-transcriptional regulation.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-020-02608-4.

Novel polymorphisms in UTR and coding region of inducible heat shock protein 70.1 gene in tropically adapted Indian zebu cattle (Bos indicus) and riverine buffalo (Bubalus bubalis)

Due to evolutionary divergence, cattle (taurine, and indicine) and buffalo are speculated to have different responses to heat stress condition. Variation in candidate genes associated with a heat-shock response may provide an insight into the dissimilarity and suggest targets for intervention. The present work was undertaken to characterize one of the inducible heat shock protein genes promoter and coding regions in diverse breeds of Indian zebu cattle and buffaloes. The genomic DNA from a panel of 117 unrelated animals representing 14 diversified native cattle breeds and 6 buffalo breeds were utilized to determine the complete sequence and gene diversity of HSP70.1 gene. The coding region of HSP70.1 gene in Indian zebu cattle, Bos taurus and buffalo was similar in length (1,926 bp) encoding a HSP70 protein of 641 amino acids with a calculated molecular weight (Mw) of 70.26 kDa. However buffalo had a longer 5' and 3' untranslated region (UTR) of 204 and 293 nucleotides respectively, in comparison to Indian zebu cattle and Bos taurus wherein length of 5' and 3'-UTR was 172 and 286 nucleotides, respectively. The increased length of buffalo HSP70.1 gene compared to indicine and taurine gene was due to two insertions each in 5' and 3'-UTR. Comparative sequence analysis of cattle (taurine and indicine) and buffalo HSP70.1 gene revealed a total of 54 gene variations (50 SNPs and 4 INDELs) among the three species in the HSP70.1 gene. The minor allele frequencies of these nucleotide variations varied from 0.03 to 0.5 with an average of 0.26. Among the 14 B. indicus cattle breeds studied, a total of 19 polymorphic sites were identified: 4 in the 5'-UTR and 15 in the coding region (of these 2 were non-synonymous). Analysis among buffalo breeds revealed 15 SNPs throughout the gene: 6 at the 5' flanking region and 9 in the coding region. In bubaline 5'-UTR, 2 additional putative transcription factor binding sites (Elk-1 and C-Re1) were identified, other than three common sites (CP2, HSE and Pax-4) observed across all the analyzed animals. No polymorphism was found within the 3'-UTR of Indian cattle or buffalo as it was found to be monomorphic. The promoter sequences generated in 117 individuals showed a rich array of sequence elements known to be involved in transcription regulation. A total of 11 nucleotide changes were observed in the promoter sequence across the analyzed species, 3 of these changes were located within the potential transcription factor binding domains. We also identified 4 microsatellite markers within the buffalo HSP70.1 gene and 3 microsatellites within bovine HSP70.1. The present study identified several distinct changes across indicine, taurine and bubaline HSP70.1 genes that could further be evaluated as molecular markers for thermotolerance.