Next Article in Journal
Genome-Wide Identification and Expression Profiling Analysis of the Trihelix Gene Family Under Abiotic Stresses in Medicago truncatula
Previous Article in Journal
GWAS Based on RNA-Seq SNPs and High-Throughput Phenotyping Combined with Climatic Data Highlights the Reservoir of Valuable Genetic Diversity in Regional Tomato Landraces
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 11
10.3390/genes11111388
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Identification, Evolution and Sequence Analysis of the Heat-Shock Protein Gene Family in Buffalo

by
Saif ur Rehman
1,†,
Asif Nadeem
2,†,
Maryam Javed
3,
Faiz-ul Hassan
4,
Xier Luo
1,
Ruqayya Bint Khalid
3 and
Qingyou Liu
1,*
1
State Key Laboratory for Conservation and Utilization of Subtropical Agro-Bioresources, Guangxi University, Nanning 530005, China
2
Department of Biotechnology, Virtual University of Pakistan, Lahore-54000, Pakistan
3
Institute of Biochemistry and Biotechnology, University of Veterinary and Animal Sciences, Lahore-54000, Pakistan
4
Institute of Animal and Dairy Sciences, Faculty of Animal Husbandry, University of Agriculture, Faisalabad-38040, Pakistan
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this manuscript.
Genes 2020, 11(11), 1388; https://0-doi-org.brum.beds.ac.uk/10.3390/genes11111388
Submission received: 26 October 2020 / Revised: 15 November 2020 / Accepted: 18 November 2020 / Published: 23 November 2020
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Heat-shock proteins (HSP) are conserved chaperones crucial for protein degradation, maturation, and refolding. These adenosine triphosphate dependent chaperones were classified based on their molecular mass that ranges between 10–100 kDA, including; HSP10, HSP40, HSP70, HSP90, HSPB1, HSPD, and HSPH1 family. HSPs are essential for cellular responses and imperative for protein homeostasis and survival under stress conditions. This study performed a computational analysis of the HSP protein family to better understand these proteins at the molecular level. Physiochemical properties, multiple sequence alignment, and phylogenetic analysis were performed for 64 HSP genes in the Bubalus bubalis genome. Four genes were identified as belonging to the HSP90 family, 10 to HSP70, 39 to HSP40, 8 to HSPB, one for each HSPD, HSPH1, and HSP10, respectively. The aliphatic index was higher for HSP90 and HSP70 as compared to the HSP40 family, indicating their greater thermostability. Grand Average of hydropathicity Index values indicated the hydrophilic nature of HSP90, HSP70, and HSP40. Multiple sequence alignment indicated the presence of highly conserved consensus sequences that are plausibly significant for the preservation of structural integrity of proteins. In addition, this study has expanded our current knowledge concerning the genetic diversity and phylogenetic relatedness of HSPs of buffalo with other mammalian species. The phylogenetic tree revealed that buffalo is more closely related to Capra hircus and distantly associated with Danio rerio. Our findings provide an understanding of HSPs in buffalo at the molecular level for the first time. This study highlights functionally important HSPs and indicates the need for further investigations to better understand the role and mechanism of HSPs.

1. Introduction

Heat-shock proteins (HSP) are well-known and highly conserved protein families responsible for cellular responses to various stresses. These evolutionary conserved molecular chaperones play a vital role in the survival of organisms in response to elevated ambient temperatures. Protection from cellular damage, alleviation of stresses such as energy depletion, cytokines, ischemia, and development of thermotolerance are major physiological processes controlled by HSP [1]. The pivotal role of HSP in cellular responses under stress conditions and normal physiological activities is well established [2]. HSP were first reported by the Ferruccio Ritossa in 1962 during his experiment on the Drosophila at the Genetics Institute in Pavia [3].
The HSP are highly conserved, ubiquitous, and functional proteins found in all organisms, starting from bacteria to humans. It has been suggested that HSP play an essential physiological role in normal situation as well as in cellular and systematic stresses, especially in farm animals, including cattle and buffalo [1]. Nonetheless, the most crucial function of HSP is the protection of cells from heat shock through resisting against denaturation of cellular proteins [4]. The cellular response is one of the primary pathways which help the animal to endure thermal stress. This pathway assists animals in surviving several anomalies in the functioning of cells. HSP is an important end product of cellular and tissue defense mechanism and expression of HSP gets markedly elevated during heat shock [5]. Various physiological and physical stressors are responsible for the activation of HSP proteins [6]. Elevation of HSP response can help to improve the thermotolerance of bovines [7]. Heat stress is one of the major factors affecting livestock production across the world. HSPs not only develop thermotolerance in livestock species including dairy cattle but also considered to be potential biological marker to quantify the magnitude of heat stress across different livestock species [8]. Several studies have proved the importance of HSPs in regulating various types of immunity [9]. It has been reported that the addition of HSP-derived peptides in lymphocytes of buffalo helped in increasing adaptive and innate immune effectors [9]. HSPs have also been associated with reproductive efficiency. For instance, in Holstein cows, differences in HSP40 genes have been correlated with early embryonic development in vitro. Similarly, polymorphism in HSP70 has been reported to be associated with difference in reproductive performances of Bos indicus breeds [10].
These proteins also act as molecular chaperones for various other cellular proteins, actively involved in multiple regulatory pathways, and have strong cytoprotective properties [2]. As a molecular chaperon, HSP prevent the inappropriate aggregation of proteins and regulates the final packaging, degradation, and/or repair of nascent polypeptides. HSP plays an essential role in maintaining homeostasis at cellular levels due to their specialized functions, particularly under heat stress conditions. HSP contribute significantly to the protection of cellular damage and other adverse effects of heat stress in animals. Failure of HSP to sustain their normal physiology can lead to poor or no thermotolerance rendering animals a victim of heat stress under higher ambient temperatures [11].
HSP have been classified into various families depending upon their molecular weights. HSP have been classified as HSP40, HSP60, HSP70, HSP90, HSP100 and small HSPs [12]. These HSP families are known for their distinct roles in different biological processes. For instance, various pathological conditions such as viral infections, cancer, muscular dystrophy, and neurodegeneration have been associated with HSP40 [13]. Similarly, HSP70 after attaching with protein substrate, assist in transportation as well as unfolding. Moreover, HSP70 has been involved in different pathological disorders such as infectious diseases, cancer, and neurodegenerative problems. HSP90 is known to activate certain oncogenic proteins through regulation of their conformation and stability. The inhibition of HSP90 could suppress such oncogenic signaling pathways. Due to the association of HSPs with certain diseases, these are potential targets for designing and developing therapies for these ailments [14,15].
Water buffalo is an important livestock species and second-best milk producer mammal after cattle (by contributing 5% of the total global milk production) distributed in different areas of tropics and subtropics [16,17]. Owing to diverse crucial role of HSP in cellular physiology, genomic characterization, and sequence analysis is imperative to better understand the evolutionary relationship, physiochemical properties, and chromosomal distribution in the buffalo. Our present study aimed to perform the HSP gene family analysis at the genomic level. Our findings provide insights into the understanding of the genomic architecture of the HSP family in buffalo and provide key evidence for further functional studies.

2. Materials and Methods

2.1. Identification of the HSP Genes in Buffalo

Published sequences from Bos taurus were used as queries for the genome-wide identification of HSP genes in the buffalo, so we used representative sequences of HSP10 (NP_776771.1), HSP40 (XP_005224490.1), HSP70 (NP_976067.3), HSP90 (NP_001012688.1), HSPB (NP_001020740.1), HSPD (NP_001160081.1), and HSPH1 (NP_001068770.1) from Bos taurus. We analyzed HSP protein sequences using Protein-protein basic local alignment (BLAST) search with an E value equal or less than 1.0 × e−5 of the non-redundant protein sequences in buffalo using all the parameters by default. After retrieving the required sequences, redundancy of the sequences was checked to avoid ambiguity. Chromosome locations of each HSP gene were obtained from buffalo genome resources. The corresponding gene positions files were obtained from the buffalo genome annotation file with the general feature format (GFF), which can be set as the input files of the MCScanX program, as reported previously [18]. Moreover, the exon-intron structure was analyzed by the buffalo genome annotation file using the in-house scripts. The conserved protein domains in buffalo HSP families were blasted against the conserved domain database (CDD) from NCBI.

2.2. Characterization of Physiochemical Properties of HSP Genes

To characterize the physical and chemical parameters of HSP proteins, ProtParam tool (https://web.expasy.org/protparam/) was used [19]. Protein sequences were submitted to identify the molecular weight (MW), number of amino acids, isoelectric point (pI), aliphatic index (AI), instability index (II) and grand average of hydropathicity (GRAVY). Moreover, conserved protein motifs of HSP were analyzed by using the Multiple Expectation Maximization for Motif Elicitation (MEME) programs with a maximum of 10 motifs [20]. Three-dimension protein structure for various HSP protein (HSP10, HSP40, HSP70, HSP90, HSPB1 and HSPH1) were constructed by using Phyre2 software (http://www.sbg.bio.ic.ac.uk/phyre2) with fold recognition end homology modeling for the identification of secondary structure attributes.

2.3. Phylogenetic Analysis

Results of multiple sequence alignment were further used for phylogenetic analysis using MEGA7 [21]. The evolutionary history was inferred by using the Maximum Likelihood method based on the Jones-Taylor Thornton (JTT) matrix-based model [22]. The tree with the highest log likelihood (−581.71) was selected. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model with 1000 bootstrap resampling, and then selecting the topology with superior log likelihood value. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 77 amino acid sequences. All positions containing gaps and missing data were eliminated. There was a total of 16 positions in the final dataset. Evolutionary analyses were conducted in MEGA7 [21].

2.4. Sequence Analysis

Accession no of all sequences which were used for multiple sequence alignment and subsequent analyses are presented in Tables S1–S5. Comparison of buffalo HSP was carried out with other mammalian species including cattle (Bos taurus), goat (Capra hircus), sheep (Ovis aries), horse (Equus caballus), camel (Camelus ferus), mouse (Rattus norvegicus), pig (Sus scrofa), zebra fish (Danio rerio), mouse (Mus musculus), and dog (Canis lupus familiaris). Sequence homology of buffalo HSP family was carried out along with elucidation of phylogenetic relationships using the Geneious version 2020.1 created by Biomatters (https://www.geneious.com).

2.5. Identification of Single Nucleotide Polymorphism (SNP)

Multiple sequence alignment of coding regions of all the members of two major HSP gene families (HSP90 and HSP70) in buffalo and Bos taurus were perused in DNASTAR Lasergene software to identify the number of synonymous and non-synonymous SNPs in buffalo [23]. Moreover, we also calculated the number of non-synonymous substitutions per non-synonymous site (Ka)/ the number of synonymous substitutions per synonymous site (Ks) ratio as described previously [24].

3. Results

3.1. Identification of HSP Gene Families and Their Physiochemical Properties

In the present study, we applied a comprehensive strategy to characterize the HSP in the buffalo genome. By using human and mouse sequences as queries, we detected four genes belonging to HSP90, 10 genes to HSP70, 39 genes to HSP40, 8 genes to HSPB, and one gene to each HSPD, HSPH1, and HSP10 family. Buffalo each HSP gene family tree and percentage of identity with other species is represented in Figures S1–S7. The physical and chemical properties of these HSP genes in buffalo are presented in Table 1, Table 2, Table 3 and Table 4. The functional diversity of HSP isoforms could be realized from a wide range of their MW (10 kDa to 254 kDa) and the total number of amino acids in different HSP genes that ranged from 102 (HSP10) to 2223 amino acids in DnaJ homolog subfamily C member 13 (DNAJC13).
The results of the present study revealed that all the members of the HSP90 family are unstable according to instability index. The isoelectric point indicates the acidic nature of protein members except for TRAP1, which is basic in nature. An AI greater than 65 indicates the thermostability of HSP90 protein members whereas, lower GRAVY values of HSP90 indicates its hydrophilic nature (Table 4). HSP70 members are thermostable, acidic, and hydrophilic in nature. The instability index indicates the unstable nature of HSPA14, HSPA4, and HSPA4L protein members (Table 1). Various physical and chemical parameters of HSP40 protein members are shown in Table 2. HSP40 is the largest HSP gene family, with 39 genes widely distributed over different buffalo chromosomes with MW ranging from 15 kDa to 254 kDa. The largest gene of this family is DNAJC13, which harbors 59 exons, 2243 amino acids, and MW of 254 kDa. The details about the physiochemical nature of small HSPs, including HSP10, HSPB, HSPD, and HSPH1, are given in Table 5. All members of HSPB are thermostable except HSPB7 and 8.

3.2. Phylogenetic Analysis

The evolutionary history was inferred by using the Maximum Likelihood method based on JTT matrix-based model. The bootstrap consensus tree inferred from 1000 replicates was taken to represent the evolutionary history of the taxa analyzed. Each HSP family genes clustered together in one group, revealing seven groups (Figure 1). Overall phylogenetic relationships revealed that buffalo HSP gene family is more closely related to Bos taurus, goat and sheep along with the common clade of horse and camel (Figure 2). Moreover, buffalo is distantly related to other mammals, including pig, dog, rat, mouse, and zebra fish.

3.3. Gene Structure and Motif Analysis

To gain further insight of HSP genes in buffalo, gene structures and conserved motifs were predicted. The number of exons varies in all genes; for instance, DNAJC13 has a higher number of exons i.e., 59. The exonic pattern of buffalo HSP genes are shown in Figure 3B. Ten putative motifs of heat-shock protein gene family were observed in buffalo (Table 5). MEME motif 1 and 2 (consisting of 41 and 28 amino acid) were annotated as the DnaJ domain in Pfam search, while 3 to 9 motifs were annotated as HSP70. These motifs were also confirmed from NCBI CDD. Moreover, MEME was used to predict conserved motifs present within HSP families in the buffalo (Figure 3D).
Collinearity analysis showed that HSP genes were distributed over 20 chromosomes in buffalo, while in cattle, these genes were randomly distributed over 25 chromosomes (Figure 4). Most of the genes were present on proximal or distal ends of chromosomes in buffalo.

3.4. HSP Protein Structural Configuration

Three-dimensional protein models were also predicted for various heat-shock proteins (HSP10, HSPB1, HSP90, HSP40 and HSP70) in Bos taurus to elucidate comparative structural configuration (Figure 4). It was observed that the overall structure of HSP10 in both species was unaltered with the same number (102) of amino acid residues. Similarly, in the case of HSPB1, no significant structural variation was observed among both species in 201 amino acid residues. For HSP90, considerable variation in the structure of the protein was observed in both species (Figure 5). In cattle, 380 residues were present, but in buffalo, only 339 amino acids were observed in final protein structure (Figure 5). Moreover, secondary structural elements also varied in both species as β-sheets observed in buffalo were absent in the case of Bos taurus. In cattle, protein was majorly consisting of α-helix. Similarly, the HSP40 gene also showed structural variation in buffalo as the final structure comprised of 480 amino acid residues as compared to 453 in cattle (Figure 5). However, the structure of HSP70 protein was similar in both species having amino acids 641.

3.5. SNP Analysis

Mutation analysis of 14 genes belonging to HSP70 and HSP90 families from buffalo and Bos taurus were compared to detect synonymous and non-synonymous SNPs. We identified a total of 194 SNPs in the HSP90 family genes among which only 19 were non-synonymous (Table 6). The highest number of SNPs (128) was observed in the TRP1 gene, among which 17 were non-synonymous. Moreover, the TRP1 gene also exhibited a single indel of 3 nucleotides at position 1703–1705. HSP90AA1 showed a total of 23 synonymous SNPs and only one non-synonymous SNP, while HSP90AB1 and HSP90B1 exhibited only 22 and 19 synonymous SNPs, respectively (Table S6).
The present study identified 291 SNPs in the HSP70 family genes, among which 59 were non-synonymous. The highest number of SNPs (80) was observed in the HSPA8 gene exhibiting 8 non-synonymous substitutions. The lowest number of SNPs (10) was detected in the HSPA14 gene with two non-synonymous SNPs. HSP70 gene showed highest non-synonymous to synonymous ratio by exhibiting 23 non-synonymous SNPs in a total of 62. Interestingly, the HSPA2 gene of the HSP70 family showed an indel of 15 nucleotides at position 777–791 (Table S6).

4. Discussion

HSP play crucial functions, including stress tolerance, cell survival, protein refolding, and prevention of protein denaturation. Moreover, HSP also facilitate the regulation of specific cellular processes such as signal transduction, protein synthesis, protein trafficking, and DNA replication [4,25,26]. It is well established that HSP90 plays pivotal role in controlling cell cycle, transport, folding, and degradation of proteins as well as thermotolerance and induction of signal transduction. Studies have revealed that HSP90 is present in various organisms and this HSP family is highly conserved [27]. HSPs are the molecular chaperones that help to maintain activity of cells in different stressful states. Some of these proteins are produced due to different stressful stimuli (inducible), while others are expressed constitutively. Cellular metabolism, ultraviolet radiations, alteration in protein structures, cytotoxic drugs, and heat are some of the stress stimuli that usually induce HSP expression. These proteins are highly involved in facilitating the degradation of misfolded proteins and regulating proper protein folding [28]. In this study, we analyzed a total of 64 HSP genes in buffalo. Members of HSP90 are found to be present on four chromosomes in buffalo similar to Bos taurus [29], suggesting that these protein molecules are highly conserved and share a common ancestor (Figure S3).

4.1. Physiochemical Properties of HSP Gene Families

The characterization of the physiochemical properties of proteins of various gene families is essential to identify the functions and properties of proteins encoded by these genes. Isoelectric point indicates that members of HSP90 in buffalo are acidic in nature (pI < 7.0), except for the TRAP1 gene, which was basic in nature similar to human HSP90 [30].
AI of the globular proteins is associated with their thermostability [31]. The AI, the respective volume of a protein occupied by aliphatic side chains (alanine, valine, isoleucine, and leucine), may be interpreted as a positive factor for increasing thermostability of globular proteins [31]. As the AI is associated with the thermostability of globular proteins, so a higher value of the AI reveals the high thermostability of protein. Hence, a higher AI of HSP90 and HSP70 showed that members of these proteins are more thermostable compared to HSP40 members. Similar results have also been observed in cattle genome [29]. In the present study, the HSP90 gene family of buffalo had AI greater than 65, indicating the stability of proteins at high temperatures. Similarly, the higher expression of HSP90 has been found in goat during the peak summer season, revealing the extent of thermostability of HSPs [32].
The GRAVY value is used to predict water and protein interaction. The GRAVY value for a protein or a peptide is defined as the sum of hydropathy values of all amino acids divided by the protein length. GRAVY qualifies and quantifies the hydrophobic or hydrophilic property of a protein. A positive GRAVY score indicates a globally hydrophobic protein, whereas a negative GRAVY score is related to hydrophilic protein [33]. Hence, more the negative GRAVY score, more soluble the protein is. All members of HSP40, HSP70, and HSP90 had a negative GRAVY score in the present study, suggesting these all members are soluble proteins. However, highly negative GRAVY values of HSP90 observed in the present study indicate the greater hydrophilic nature of its members. It can help to enhance the functional properties of protein binding and oligomerization [34,35]. Among heat-shock proteins, HSP90 is considered an essential and highly conserved molecular chaperone widely present in eukaryotic cells [36]. HSP90 in buffalo share high sequence similarities with other species (Figure S3). The presence of conserved motifs helps interact with other proteins and increases their chaperone function [37,38].
A total of 10 HSP70 members were observed in buffalo. Interestingly, all members of HSP70 showed acidic nature and almost similar isoelectric points indicating the possibility of functional similarity among buffalo HSP70 members. It has been suggested earlier that HSP70 has conserved biological functions [39]. GRAVY values of HSP70 protein of buffalo indicates their hydrophilic nature and high values of AI revealed their thermal stability in a wide range of temperatures. These findings confirm the consistency of its chaperoning role in protection against specific cellular and systemic stresses that can lead to protein denaturation [40]. In evolution, HSP70 is the most conserved protein, and its conserved domains and functions are considered responsible for across species conservation [11,41]. In the present study, a tree with an alignment view of HSP70 showed the presence of conserved domains. Buffalo is closely related to Bos taurus, followed by goat and sheep (Figure S2). The comparative analysis of the HSP70 protein in water buffalo revealed that most of the sequences were conserved across species. Comparative sequence analysis for bubaline HSP70.1 indicated the high similarity of genes among mammalian species, such as cattle, goat, human, mouse, and pig, respectively. Furthermore, phylogenetic analysis of bovine and bubaline HSP70.1 gene with other species revealed Indian zebu cattle grouped close to Bos taurus followed by buffalo and goat, while porcine, ruminants, and murine formed distinct cluster [42,43].
Heat-shock protein 40 (HSP40/DnaJ) is evolutionary conserved and thought to play an essential role in translocation, degradation, translation, folding, and unfolding of proteins by stimulating the ATPase activity of HSP70 protein chaperones [44,45]. All members of HSP40 have domain J, which helps in binding with HSP70. Also, 41 DnaJ domains have been reported in humans [46]. In the present study, 39 putative members of HSP40 were found. GRAVY value of HSP40 members represents the hydrophilic nature of these proteins. Furthermore, pI values of HSP40 protein members suggest that some of these proteins possess acidic nature while other proteins appear to be basic. It may explain the different functional properties of these proteins and provides insights for designing further experiments to elucidate their functional diversity. HSP40 represents a large and diverse protein family and is usually divided into three classes [47]. Type I HSP40s is known to have a glycine and phenylalanine (G/F)-rich region and a cysteine-rich, zinc finger-like region (ZFLR). Type II HSP40s only contain a G/F-rich region, while type III HSP40s have neither a G/F-rich region nor a ZFLR [48]. Comparative sequence and phylogenetic analyses revealed that the HSP40 family of buffalo is closely related to Bos taurus (Figure S1). Another study conducted on the bovine genome also indicated the presence of evolutionary conserved residues in HSP40 [29].

4.2. Sequence Analyses

Data of protein sequences can provide useful information regarding evolutionary conserved regions important for various biological functions. Multiple sequence alignment is one of the ways to identify these conserved constraints and to compile data for structural and functional analysis of proteins [49]. To investigate the protein sequence features of the heat-shock proteins, 10 motifs of HSP gene families in buffalo were predicted by the MEME tool, and the regular expression levels of the conserved motifs are listed in Table 5. These conserved motifs can play a crucial role in the structural and functional organization of proteins. Similarly, another study on Ciona savignyi, identified conserved motifs in HSP60, HSP70, and HSP90 by using MEME [50].

4.3. Phylogenetic Analysis

Molecular phylogenetic tree of heat-shock proteins (HSP40, HSP70, HSP90, HSPB, HSPD, HSP10, and HSPH1) was constructed by using the Maximum Likelihood method based on the JTT matrix-based model [22]. The bootstrap consensus tree inferred from 1000 replicates was taken to represent the evolutionary history of the taxa analyzed (Figure 1). Overall phylogenetic relationships revealed that buffalo is more closely related to goat than Bos taurus, and sheep along with horse and camel (Figure 2). Moreover, buffalo is distantly related to other mammals, including pig, dog, rat, mouse and zebra fish. A study conducted on cattle also revealed similar results as phylogenetic analysis of HSPA8 and HSPA1A indicted the high similarity between buffalo and sheep [51]. Previously, molecular characterization of HSP70.1 of goat have shown that buffalo, cattle, sheep, and goat share a common ancestor whereas human, horse and pig showed dissimilarities thus having different ancestors [52]. Moreover, Karan Fries cattle was evolutionary related to Bos taurus, Bos indicus, and buffalo, whereas distinctly related to rat and mouse [53].

4.4. Protein Structural Configuration

Three-dimensional (3-D) protein configurations of HSP in Bos taurus and buffalo was predicted. The analysis of 3-D secondary protein structures revealed the structural variations only in HSP90 and HSP40 of Bos taurus and buffalo, while other HSPs were similar or structurally related in both species. These findings confirm the results obtained in gene structure and phylogenetic analysis that HSPs are conserved proteins. Eukaryotic HSP90 proteins had a modular structure with an N-terminal ATP binding domain connected to a middle domain by a ‘linker’ of variable length. The N-terminal and middle domains together form a “split” ATPase site that is also the binding pocket for GA [54]. HSP40 can function as a molecular chaperone to bind non-native polypeptides. HSP40 pairs with HSP70 to achieve its function of protein folding, transport, and degradation functions. The J-domain of HSP40 may stimulate the ATPase activity of HSP70 and the C-terminal peptide-binding fragment of HSP40 can load the bound non-native polypeptide to HSP70 [55]. HSP70 and HSP40 are essential chaperons affecting the protein folding in the cell. Three-dimensional structural configurations of HSP40 in cattle and buffalo varied in number of amino acid residues and overall secondary structural attributes that may be associated with differential chaperon functions and stress responses in both species. Further investigations at the transcriptome level are required to corroborate these findings.

4.5. SNP Analysis

The present study revealed that HSP90AA1, HSP90AB1 and HSP90B1 are more conserved genes than TRP1 due to a higher number of non-synonymous SNPs in TRP1. Similarly, our findings indicate that HSP70.1, HSPA14, and HSPA1L are more conserved than other genes of HSP70 family. We observed a relatively higher ratio of non-synonymous to synonymous SNPs in the HSP70 gene (0.59) as compared to other genes of the HSP70 family. It indicates that HSP70 gene is under selection owing to its strong association with cellular thermotolerance [1]. Moreover, the expression dynamics of the HSP70 gene are used as a potential biomarker to reveal the level of heat stress and the ability of animals to withstand climate stress [56]. It is mainly attributed to the fact that HSP expression indicates cellular response and intensity of stress in animals [57]. At the transcription level, persuaded expression of HSP70 genes is coordinated as upstream elements at promoter region directly control the expression profiles of HSP70. Studies have shown a higher expression of HSP70 in Murrah and Tarai buffaloes [58] during the summer compared to other seasons (or thermoneutral zone). Moreover, increased expression of HSP70 in bovine lymphocytes has also been observed after heat stress [59].
Due to subsequent induction of HSP70 after heat shock and its association with physiological parameters (rectal temperature, pulse, and respiratory rates) of animals [59], serum HSP70 is considered a sensitive biomarker for heat stress management in livestock [58]. In addition to thermotolerance, studies have also reported other roles of the HSP family in animal physiology. In particular, SNPs identified in the 5′-UTR regions of HSP70 gene have been linked with different performance traits of livestock such as milk production, thermal stress and disease vulnerability [1,43]. Moreover, SNPs have also been identified in the HSP70 gene of beef cattle breeds (Angus, Brahman, and their crossbreds) with little effect on milk protein and milk fat [60]. Polymorphism (SNP) at the promoter region has been shown to alter the expression level of bovine HSP70.1 gene, and exhibited its significant association with heat tolerance and higher milk production in Frieswal cattle [43]. These findings indicate the putative effects of selective pressure (natural or artificial selection) on the HSP70 gene to ultimately favor animals with better thermotolerance, performance, and stress resilience. Our findings are in line with recent whole genome-sequencing studies on buffalo that reported evidence for the selection of genes related to social behavior, digestive physiology, fitness, and milk production [61].
HSPs are highly conserved in nature and play crucial role in stress stimulus and in several normal cellular functions. The genomic characterization of the HSP gene family and elucidation of their physiochemical properties performed in the present study would help in better understanding HSPs in buffalo. HSPs are highly conserved proteins and important from evolutionary point of view, so further investigations are warranted to study expression dynamics under various physiological conditions to elucidate their putative roles in thermotolerance and stress physiology in buffalo. In particular, non-synonymous SNPs and indels observed in the present study need to be further explored to detect their effect on protein function and subsequent physiological manifestations in buffalo. These findings might explain the better thermotolerance and fitness of buffalo compared to Bos taurus, particularly under tropical climate.

5. Conclusions

The present study revealed that the buffalo genome contains 64 genes of the HSP family in buffalo. These genes are widely dispersed across the buffalo genome. Multiple sequence alignment and phylogenetic analysis indicated the highly conserved consensus sequences and evolutionary relatedness of the HSP family. Sequence analysis revealed that the HSP90 gene family buffalo is more conserved in nature as compared to the HSP70 family that exhibited a quite higher ratio of non-synonymous substitutions. Structural variations in secondary structures of HSP40 and HSP90 were observed in buffalo and cattle. The findings of the present study provide insights into the genomic variations in the HSP gene family that could be used for a better understanding of the function of these genes and their further use for selective breeding of buffalo for better thermotolerance and stress resilience.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/11/11/1388/s1, Figure S1: Tree with alignment view of HSP40. Percentage of sequence is represented with different colors in Identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity. Figure S2: Tree with alignment view of HSP70. Percentage of sequence is represented with different colors in Identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity. Figure S3: Tree with alignment view of HSP90. Percentage of sequence is represented with different colors in Identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity. Figure S4: Tree with alignment view of HSPD. Percentage of sequence is represented with different colors in Identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity. Figure S5: Tree with alignment view of HSP10. Percentage of sequence is represented with different colors in Identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity. Figure S6: Tree with alignment view of HSPH1. Percentage of sequence is represented with different colors in Identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity. Figure S7: Tree with alignment view of HSPB. Percentage of sequence is represented with different colors in Identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity. Table S1: Accession no. of different heat-shock protein in different organisms used in phylogenetic analysis, Table S2: Accession no. of buffalo HSP90, Table S3: Accession no. of buffalo HSP70, Table S4: Accession no. of buffalo HSP40, Table S5: Accession no. of buffalo HSP family B, HSP10, HSPD, and HSPH1, Table S6: Mutations in coding region of HSP gene family in buffalo.

Author Contributions

Conceptualization, F.-u.H., A.N. and Q.L.; resources, S.u.R. and Q.L.; data curation, F.-u.H. and R.B.K; methodology, R.B.K., S.u.R., X.L. and F.-u.H.; project administration, Q.L. and A.N.; software, S.u.R.; supervision, Q.L. and X.L.; writing—original draft preparation, S.u.R., A.N., F.-u.H. and M.J.; writing—review and editing, S.u.R., A.N., M.J., F.-u.H., X.L., R.B.K. and Q.L.; funding acquisition, Q.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Chinese National Natural Science Foundation (Grant No.31860638), Guangxi Natural Science Foundation (AB16380042) and Nanning Science Foundation (2018018). None of the funding body had any influence on study design, data analysis, and interpretation, or in writing the manuscript.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Hassan, F.-U.; Nawaz, A.; Rehman, M.S.; Ali, M.A.; Dilshad, S.M.; Yang, C. Prospects of HSP70 as a genetic marker for thermo-tolerance and immuno-modulation in animals under climate change scenario. Anim. Nutr. 2019, 5, 340–350. [Google Scholar] [CrossRef]
  2. Kregel, K.C. Invited review: Heat shock proteins: Modifying factors in physiological stress responses and acquired thermotolerance. J. Appl. Physiol. 2002, 92, 2177–2186. [Google Scholar] [CrossRef] [Green Version]
  3. Ritossa, F. Discovery of the heat shock response. Cell Stress Chaperones 1996, 1, 97–98. [Google Scholar] [CrossRef] [Green Version]
  4. Feder, M.E.; Hofmann, G.E. Heat-shock proteins, molecular chaperones, and the stress response: Evolutionary and ecological physiology. Annu. Rev. Physiol. 1999, 61, 243–282. [Google Scholar] [CrossRef] [Green Version]
  5. Archana, P.R.; Aleena, J.; Pragna, P.; Vidya, M.K.; Niyas, A.P.A.; Bagath, M.; Krishnan, G.; Manimaran, A.; Beena, V.; Kurien, E.K.; et al. Role of Heat Shock Proteins in Livestock Adaptation to Heat Stress. J. Dairy Vet. Anim. Res. 2017, 5, 1–8. [Google Scholar] [CrossRef]
  6. Yang, M.; Tan, H.; Yang, Q.; Wang, F.; Yao, H.; Wei, Q.; Tanguay, R.M.; Wu, T. Association of hsp70 polymorphisms with risk of noise-induced hearing loss in Chinese automobile workers. Cell Stress Chaperones 2006, 11, 233–239. [Google Scholar] [CrossRef] [PubMed]
  7. Collier, R.J.; Collier, J.L.; Rhoads, R.P.; Baumgard, L.H. Invited Review: Genes Involved in the Bovine Heat Stress Response. J. Dairy Sci. 2008, 91, 445–454. [Google Scholar] [CrossRef] [PubMed]
  8. Mishra, S.R. Significance of molecular chaperones and micro RNAs in acquisition of thermo-tolerance in dairy cattle. Anim. Biotechnol. 2020, 1–11. [Google Scholar] [CrossRef]
  9. Mishra, A.; Hooda, O.K.; Singh, G.; Meur, S.K. Influence of induced heat stress on HSP70 in buffalo lymphocytes. J. Anim. Physiol. Anim. Nutr. 2011, 95, 540–544. [Google Scholar] [CrossRef]
  10. Cushman, R.A. Physiology and endocrinology symposium: The current status of heat shock in early embryonic survival and reproductive efficiency. J. Anim. Sci. 2013, 91, 1141–1142. [Google Scholar] [CrossRef] [Green Version]
  11. Kiang, J.G.; Tsokos, G.S. Heat shock protein 70 kDa molecular biology, biochemistry, and physiology. Pharmacol. Ther. 1998, 80, 183–201. [Google Scholar] [CrossRef]
  12. Nagarajan, N.S.; Arunraj, S.P.; Sinha, D.; Rajan, V.B.; Esthaki, V.K.; D’Silva, P. HSPIR: A manually annotated heat shock protein information resource. Bioinformatics 2012, 28, 2853–2855. [Google Scholar] [CrossRef] [Green Version]
  13. Dong, C.-W.; Zhang, Y.-B.; Zhang, Q.-Y.; Gui, J.-F. Differential expression of three Paralichthys olivaceus Hsp40 genes in responses to virus infection and heat shock. Fish Shellfish. Immunol. 2006, 21, 146–158. [Google Scholar] [CrossRef] [PubMed]
  14. Banerji, U. Heat Shock Protein 90 as a Drug Target: Some Like It Hot. Clin. Cancer Res. 2009, 15, 9–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Chen, W.; Feng, P.; Liu, T.; Jin, D. Recent Advances in Machine Learning Methods for Predicting Heat Shock Proteins. Curr. Drug Metab. 2019, 20, 224–228. [Google Scholar] [CrossRef]
  16. Mintoo, A.A.; Zhang, H.; Chen, C.; Moniruzzaman, M.; Deng, T.; Anam, M.; Huque, Q.M.E.; Guang, X.; Wang, P.; Zhong, Z.; et al. Draft genome of the river water buffalo. Ecol. Evol. 2019, 9, 3378–3388. [Google Scholar] [CrossRef]
  17. Rehman, S.U.; Shafique, L.; Yousuf, M.R.; Liu, Q.; Ahmed, J.Z.; Riaz, H. Spectrophotometric Calibration and Comparison of Different Semen Evaluation Methods in Nili-Ravi Buffalo Bulls. Pak. Vet. J. 2019, 39, 568–572. [Google Scholar] [CrossRef]
  18. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.-H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  19. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef] [Green Version]
  20. Bailey, T.L.; Williams, N.; Misleh, C.; Li, W.W. MEME: Discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34, W369–W373. [Google Scholar] [CrossRef]
  21. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Jones, D.T.; Taylor, W.R.; Thornton, J.M. The rapid generation of mutation data matrices from protein sequences. Bioinformatics 1992, 8, 275–282. [Google Scholar] [CrossRef] [PubMed]
  23. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lu, X.; Duan, A.; Liang, S.; Ma, X.; Deng, T. Genomic Identification, Evolution, and Expression Analysis of Collagen Genes Family in Water Buffalo during Lactation. Genes 2020, 11, 515. [Google Scholar] [CrossRef]
  25. Hartl, F.U. Molecular chaperones in cellular protein folding. Nature 1996, 381, 571–580. [Google Scholar] [CrossRef]
  26. Picard, D. Heat-shock protein 90, a chaperone for folding and regulation. Cell. Mol. Life Sci. 2002, 59, 1640–1648. [Google Scholar] [CrossRef]
  27. Chen, B.; Zhong, D.; Monteiro, A. Comparative genomics and evolution of the HSP90 family of genes across all kingdoms of organisms. BMC Genom. 2006, 7, 156. [Google Scholar] [CrossRef] [Green Version]
  28. Lanneau, D.; Wettstein, G.; Bonniaud, P.; Garrido, C. Heat Shock Proteins: Cell Protection through Protein Triage. Sci. World J. 2010, 10, 1543–1552. [Google Scholar] [CrossRef] [Green Version]
  29. Ajayi, O.O.; Peters, S.O.; De Donato, M.; Sowande, S.O.; Mujibi, F.D.; Morenikeji, O.B.; Thomas, B.N.; Adeleke, M.A.; Imumorin, I.G. Computational genome-wide identification of heat shock protein genes in the bovine genome. F1000Research 2018, 7, 1504. [Google Scholar] [CrossRef]
  30. Faya, N.; Penkler, D.L.; Bishop, Ö.T. Human, vector and parasite Hsp90 proteins: A comparative bioinformatics analysis. FEBS Open Bio 2015, 5, 916–927. [Google Scholar] [CrossRef] [Green Version]
  31. Ikai, A. Thermostability and Aliphatic Index of Globular Proteins. J. Biochem. 1980, 88, 1895–1898. [Google Scholar] [CrossRef] [PubMed]
  32. Yadav, V.P.; Dangi, S.S.; Chouhan, V.S.; Gupta, M.; Dangi, S.K.; Singh, G.; Maurya, V.P.; Kumar, P.; Sarkar, M. Expression analysis of NOS family and HSP genes during thermal stress in goat (Capra hircus). Int. J. Biometeorol. 2016, 60, 381–389. [Google Scholar] [CrossRef] [PubMed]
  33. Lanneau, D.; Brunet, M.; Frisan, E.; Solary, E.; Fontenay, M.; Garrido, C. Heat shock proteins: Essential proteins for apoptosis regulation. J. Cell. Mol. Med. 2008, 12, 743–761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Roy, S.; Maheshwari, N.; Chauhan, R.; Sen, N.K.; Sharma, A. Structure prediction and functional characterization of secondary metabolite proteins of Ocimum. Bioinformation 2011, 6, 315–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kadel, S.; Persico, M.; Thibodeau, J.; Lainé, C.; Bazinet, L. Use of redundancy analysis and multivariate regression models to select the significant membrane properties affecting peptide migration during electrodialysis with filtration membranes. Sep. Purif. Technol. 2019, 221, 114–125. [Google Scholar] [CrossRef]
  36. Bose, S.; Weikl, T.; Bügl, H.; Buchner, J. Chaperone Function of Hsp90-Associated Proteins. Science 1996, 274, 1715–1717. [Google Scholar] [CrossRef]
  37. Daugaard, M.; Rohde, M.; Jäättelä, M. The heat shock protein 70 family: Highly homologous proteins with overlapping and distinct functions. FEBS Lett. 2007, 581, 3702–3710. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, Q.; Tao, J.; Li, Y.; Xu, Y.; Liu, X.; Zong, S. Different Responses of Cytoplasmic and Endoplasmic Reticulum Hsp90 Genes from Eogystia hippophaecola (Lepidoptera: Cossidae) to Cold Stress. Forests 2019, 10, 1039. [Google Scholar] [CrossRef] [Green Version]
  39. Han, Y.-L.; Yang, W.-X.; Long, L.-L.; Sheng, Z.; Zhou, Y.; Zhao, Y.-Q.; Wang, Y.-F.; Zhu, J.-Q. Molecular cloning, expression pattern, and chemical analysis of heat shock protein 70 (HSP70) in the mudskipper Boleophthalmus pectinirostris: Evidence for its role in regulating spermatogenesis. Gene 2016, 575, 331–338. [Google Scholar] [CrossRef]
  40. Bukau, B.; Weissman, J.; Horwich, A.L. Molecular Chaperones and Protein Quality Control. Cell 2006, 125, 443–451. [Google Scholar] [CrossRef] [Green Version]
  41. Murphy, M.E. The HSP70 family and cancer. Carcinogenesis 2013, 34, 1181–1188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Singh, R.; Gurao, A.; Rajesh, C.; Mishra, S.K.; Rani, S.; Behl, A.; Kumar, V.; Kataria, R.S. Comparative modeling and mutual docking of structurally uncharacterized heat shock protein 70 and heat shock factor-1 proteins in water buffalo. Vet. World. 2019, 12, 2036–2045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Sodhi, M.; Mukesh, M.; Kishore, A.; Mishra, B.; Kataria, R.; Joshi, B. 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). Gene 2013, 527, 606–615. [Google Scholar] [CrossRef] [PubMed]
  44. Höhfeld, J. Regulation of the heat shock conjugate Hsc70 in the mammalian cell: The characterization of the anti-apoptotic protein BAG-1 provides novel insights. Biol. Chem. 1998, 379, 269. [Google Scholar] [PubMed]
  45. Minami, Y.; Höhfeld, J.; Ohtsuka, K.; Hartl, F.-U. Regulation of the Heat-shock Protein 70 Reaction Cycle by the Mammalian DnaJ Homolog, Hsp40. J. Biol. Chem. 1996, 271, 19617–19624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Qiu, X.-B.; Shao, Y.-M.; Miao, S.; Wang, L. The diversity of the DnaJ/Hsp40 family, the crucial partners for Hsp70 chaperones. Cell. Mol. Life Sci. 2006, 63, 2560–2570. [Google Scholar] [CrossRef]
  47. Cheetham, M.E.; Caplan, A.J. Structure, function and evolution of DnaJ: Conservation and adaptation of chaperone function. Cell Stress Chaperones 1998, 3, 28–36. [Google Scholar] [CrossRef]
  48. Sahi, C.; Craig, E.A. Network of general and specialty J protein chaperones of the yeast cytosol. Proc. Natl. Acad. Sci. USA 2007, 104, 7163–7168. [Google Scholar] [CrossRef] [Green Version]
  49. Neuwald, A.F. Gleaning structural and functional information from correlations in protein multiple sequence alignments. Curr. Opin. Struct. Biol. 2016, 38, 1–8. [Google Scholar] [CrossRef]
  50. Huang, X.; Li, S.; Gao, Y.; Zhan, A. Genome-Wide Identification, Characterization and Expression Analyses of Heat Shock Protein-Related Genes in a Highly Invasive Ascidian Ciona savignyi. Front. Physiol. 2018, 9, 1043. [Google Scholar] [CrossRef] [Green Version]
  51. Kumar, A.; Ashraf, S.; Goud, T.S.; Grewal, A.; Singh, S.; Yadav, B.; Upadhyay, R. Expression profiling of major heat shock protein genes during different seasons in cattle (Bos indicus) and buffalo (Bubalus bubalis) under tropical climatic condition. J. Therm. Biol. 2015, 51, 55–64. [Google Scholar] [CrossRef] [PubMed]
  52. Gade, N.; Mahapatra, R.K.; Sonawane, A.; Singh, V.; Doreswamy, R.; Saini, M. Molecular Characterization of Heat Shock Protein 70-1 Gene of Goat (Capra hircus). Mol. Biol. Int. 2010, 2010, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kumar, R.; Gupta, I.; Verma, A.; Kumari, R.; Verma, N. Molecular characterization and SNP identification in HSPB6 gene in Karan Fries (Bos taurus × Bos indicus) cattle. Trop. Anim. Health Prod. 2017, 49, 1059–1063. [Google Scholar] [CrossRef] [PubMed]
  54. Pearl, L.H.; Prodromou, C. Structure and mechanism of the hsp90 molecular chaperone machinery. Annu. Rev. Biochem. 2006, 75, 271–294. [Google Scholar] [CrossRef] [PubMed]
  55. Li, J.; Qian, X.; Sha, B. Heat Shock Protein 40: Structural Studies and Their Functional Implications. Protein Pept. Lett. 2009, 16, 606–612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Hansen, P.J. Physiological and cellular adaptations of zebu cattle to thermal stress. Anim. Reprod. Sci. 2004, 82, 349–360. [Google Scholar] [CrossRef] [PubMed]
  57. Neuer, A.; Spandorfer, S.D.; Giraldo, P.; Dieterle, S.; Rosenwaks, Z.; Witkin, S.S. The role of heat shock proteins in reproduction. Hum. Reprod. Updat. 2000, 6, 149–159. [Google Scholar] [CrossRef] [Green Version]
  58. Manjari, R.; Yadav, M.; Ramesh, K.; Uniyal, S.; Rastogi, S.K.; Sejian, V.; Hyder, I. HSP70 as a marker of heat and humidity stress in Tarai buffalo. Trop. Anim. Heal. Prod. 2014, 47, 111–116. [Google Scholar] [CrossRef]
  59. Liu, Y.-X.; Li, D.-Q.; Cui, Q.-W.; Shi, H.-X.; Wang, G.-L. Analysis of HSP70 mRNA level and association between linked microsatellite loci and heat tolerance traits in dairy cows. Yi Chuan Hereditas 2010, 32, 935–941. [Google Scholar] [CrossRef]
  60. Rosenkrans, C., Jr.; Banks, A.; Reiter, S.; Looper, M. Calving traits of crossbred Brahman cows are associated with Heat Shock Protein 70 genetic polymorphisms. Anim. Reprod. Sci. 2010, 119, 178–182. [Google Scholar] [CrossRef]
  61. Luo, X.; Zhou, Y.; Zhang, B.; Zhang, Y.; Wang, X.; Feng, T.; Li, Z.; Cui, K.; Wang, Z.; Luo, C.; et al. Understanding divergent domestication traits from the whole-genome sequencing of swamp- and river-buffalo populations. Natl. Sci. Rev. 2020, 7, 686–701. [Google Scholar] [CrossRef] [Green Version]
Genes 11 01388 g001 550
Figure 1. Molecular phylogenetic tree of heat-shock proteins (HSP40, HSP70, HSP90, HSPB, HSPD, HSP10, and HSPH1).
Figure 1. Molecular phylogenetic tree of heat-shock proteins (HSP40, HSP70, HSP90, HSPB, HSPD, HSP10, and HSPH1).
Genes 11 01388 g001
Genes 11 01388 g002 550
Figure 2. Phylogenetic tree with alignment view of HSP40 (sequence 1–493), HSP70 (494–1137) and HSP90 (1138–1870), HSPB (1871–2087), HSPD (2088–2662), HSPH1 (2663–3588) and HSP10KD (3589–3850). Percentage of sequence identity is represented with different colors in identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity.
Figure 2. Phylogenetic tree with alignment view of HSP40 (sequence 1–493), HSP70 (494–1137) and HSP90 (1138–1870), HSPB (1871–2087), HSPD (2088–2662), HSPH1 (2663–3588) and HSP10KD (3589–3850). Percentage of sequence identity is represented with different colors in identity bar at top of the sequences green color indicates 100% identity, green-brown shows at least 30 and under 100% while red representing below 30 identity.
Genes 11 01388 g002
Genes 11 01388 g003a 550Genes 11 01388 g003b 550
Figure 3. Phylogenetic relationships, exon motif patterns, and conserved protein motif of buffalo HSP gene family. (A) Phylogenetic tree of 64 buffalo HSP proteins. (B) Exon pattern of buffalo HSP. (C) Motif pattern of buffalo HSP proteins gene family. Ten putative motifs are indicated in different colored boxes. For details of motifs, refer to Table 5. (D) Conserved protein motif.
Figure 3. Phylogenetic relationships, exon motif patterns, and conserved protein motif of buffalo HSP gene family. (A) Phylogenetic tree of 64 buffalo HSP proteins. (B) Exon pattern of buffalo HSP. (C) Motif pattern of buffalo HSP proteins gene family. Ten putative motifs are indicated in different colored boxes. For details of motifs, refer to Table 5. (D) Conserved protein motif.
Genes 11 01388 g003aGenes 11 01388 g003b
Genes 11 01388 g004 550
Figure 4. Collinearity analysis showing distribution of HSP genes in buffalo and cattle.
Figure 4. Collinearity analysis showing distribution of HSP genes in buffalo and cattle.
Genes 11 01388 g004
Genes 11 01388 g005 550
Figure 5. Three-Dimensional protein configuration of Heat-Shock Proteins in Bos taurus and buffalo.
Figure 5. Three-Dimensional protein configuration of Heat-Shock Proteins in Bos taurus and buffalo.
Genes 11 01388 g005
Table
Table 1. Physiochemical Properties of HSP70 gene family members.
Table 1. Physiochemical Properties of HSP70 gene family members.
Sr. NoGeneChromosomeExon CountMW (kDa)A.ApIAIIIGRAVY
1HSP70.12270.276415.6785.2332.84−0.395
2HSPA1L2470.366415.8985.8332.95−0.354
3HSPA66170.966435.7883.8939.18−0.403
4HSPA816971.246505.3781.0237.59−0.454
5HSPA211169.816365.582.8135.23−0.474
6HSP701435.793236.1176.9335.1−0.473
7HSPA14141654.945095.6795.2345.33−0.085
8HSPA131551.824715.08100.1334.84−0.111
9HSPA491994.468405.1375.8945.29−0.56
10HSPA4L172294.898405.5580.4444.42−0.539
[MW (Molecular Weight in Kilo Daltons), A.A (number of amino acids), pI (Isoelectric point), AI (Aliphatic Index), II (Instability Index) and GRAVY (Grand Average of hydropathicity Index)].
Table
Table 2. Physiochemical Properties of HSP40 gene family members.
Table 2. Physiochemical Properties of HSP40 gene family members.
Sr. No.GeneChromosomeExon CountMW (kDa)A.ApIAIIIGRAVY
1DNAJA13944.883976.6571.8931.57−0.728
2DNAJA420947.974266.6374.7739.19−0.727
3DNAJA218945.764126.0668.0839.03−0.699
4DNAJB53447.174209.4269.4346.44−0.555
5DNAJB46437.823378.6663.3843.18−0.708
6DNAJB1111040.503585.9280.2537.41−0.568
7DNAJB19638.253408.7563.3549.59−0.751
8DNAJB1316836123167.6682.0336.99−0.509
9DNAJA3241252.094809.566.0245.15−0.467
10DNAJB821125.232316.464945.52−0.659
11DNAJB681331.132789.3169.1447.61−0.468
12DNAJB124941.503718.9166.3140.42−0.722
13DNAJB221032.382945.7855.1758.66−0.841
14DNAJB36126.832446.0753.1650.15−0.616
15DNAJB98325.832238.6841.6139.56−0.861
16DNAJB1471136.853298.7659.2735.91−0.874
17DNAJC514822.131984.9365.1534.45−0.476
18DNAJB74225.212177.0344.0139.69−1.102
19DNAJC5G121020.871867.9263.5540.85−0.534
20DNAJC5B15822.691995.6267.1950.49−0.474
21DNAJC189841.683588.2865.4242.81−0.892
22DNAJC21191261.675335.4553.1550.3−1.251
23DNAJC1651691.947946.8286.1840.41−0.346
24DNAJC1022491.307936.5777.3636.89−0.419
25DNAJC17111134.603046.9678.3946.35−0.906
26DNAJC3131357.705045.6478.8344.97−0.677
27DNAJC3024225.8722610.6468.7660.79−0.76
28DNAJC1151663.225598.587.7547.33−0.405
29DNAJC1141262.915438.9881.0563.77−0.833
30DNAJC731556.454946.5664.7641.63−0.722
31DNAJC45826.9323610.665.4269.46−0.856
32DNAJC144878.066998.3467.7351.48−0.614
33DNAJC2416517.201494.6978.5959.41−0.587
34DNAJC94534.303006.1170.2749.68−0.785
35DNAJC2712730.732738.5575.7133.54−0.421
36DNAJC253442.553599.2184.5462.04−0.627
37DNAJC13159254.3622436.3592.7642.37−0.234
38DNAJC124623.231985.1157.1773.75−1.024
39DNAJC1513615.8914910.188.5948.02−0.182
[MW (Molecular Weight in Kilo Daltons), A.A (number of amino acids), pI (Isoelectric point), AI (Aliphatic Index), II (Instability Index) and GRAVY (Grand Average of hydropathicity Index)].
Table
Table 3. Physiochemical Properties of HSP family B, HSP10, HSPD, and HSPH1.
Table 3. Physiochemical Properties of HSP family B, HSP10, HSPD, and HSPH1.
Sr. NoGeneChromosomeExon CountMW (kDa)A.ApIAIIIGRAVY
1HSPB124322.392015.9870.461.45−0.591
2HSPB216220.201825.0778.2446.26−0.505
3HSPB319116.731495.3396.2447.66−0.217
4HSPB618317.491645.9589.9463.95−0.121
5HSPB72318.231675.660.8459.89−0.501
6HSPB817321.671965.1359.1372.61−0.549
7HSPB93116.871578.5668.2846.32−0.38
8HSPB1161120.141766.3587.4446.75−0.16
9HSPD121360.985735.71100.9130.59−0.087
10HSP1021210.931028.8995.3923.75−0.041
11HSPH1131891.698175.3278.440.77−0.525
[MW (Molecular Weight in Kilo Daltons), A.A (number of amino acids), pI (Isoelectric point), AI (Aliphatic Index), II (Instability Index) and GRAVY (Grand Average of hydropathicity Index)].
Table
Table 4. Physiochemical Properties of HSP90 gene family members.
Table 4. Physiochemical Properties of HSP90 gene family members.
Sr. NoGeneChromosomeExon CountMW (kDa)A.ApIAIIIGRAVY
1HSP90AA1201184.747334.9379.443.03−0.751
2HSP90AB121283.257244.9681.0542.3−0.679
3HSP90B141892.438044.7676.6340.41−0.722
4TRAP1241879.537037.390.4844.82−0.336
[MW (Molecular Weight in Kilo Daltons), A.A (number of amino acids), pI (Isoelectric point), AI (Aliphatic Index), II (Instability Index) and GRAVY (Grand Average of hydropathicity Index)].
Table
Table 5. Ten differentially conserved motifs of heat-shock protein gene family observed in buffalo.
Table 5. Ten differentially conserved motifs of heat-shock protein gene family observed in buffalo.
MotifProtein SequenceLengthPfam Domain
MEME-1DYYEILGVPRGASDEEIKKAYRKLALKYHPDKNPDPGAEAE41DnaJ
MEME-2FKZIAEAYEVLSDPEKRKIYDKYGEEGL28DnaJ
MEME-3IDLGTTYSCVGVFQHGKVEIIANDQGBRTTPSYVAFTDTERLIGDAAKNQ50Hsp70
MEME-4PVTBAVITVPAYFNDSQRQATKDAGQIAGLNVLRJINEPTAAAJAYGJDK50Hsp70
MEME-5TAGDTHLGGEDFDNRLVNHFCEEFKRKHKKDISENKRALRRLRTACERAK50Hsp70
MEME-6TAGGVMTVLIKRNSTIPTKQTQTFTTYSDNQPGVLIQVYEGERAMTKDNN50Hsp70
MEME-7GIPPAPRGVPQIEVTFDIDANGILNVTAVDKSTGKENKITITNDKGRLSK50Hsp70
MEME-8DFFNGKELNKSINPDEAVAYGAAVQAAIL29Hsp70
MEME-9TLSSSTQAPLEIDSLYEGIDFYTSITRARFEELCADLFRGT41Hsp70
MEME-10LDKCQEVINWLDKNQLAEKEEFEHKQKELEQVCNPIISKLY41-
Table
Table 6. SNPs in coding region of HSP 70 and HSP90 gene family in buffalo.
Table 6. SNPs in coding region of HSP 70 and HSP90 gene family in buffalo.
GeneSynonymous (ks)Non-synonymous (ka)TotalKa/Ks Ratio
HSP70 gene family
HSPA4L186240.33
HSPA4133160.23
HSPA13249330.38
HSPA1482100.25
HSP703923620.59
HSPA8728800.11
HSP70.1171180.06
HSPA1L223250.14
HSPA6194230.21
Total232592910.25
HSP90 gene family
TRP1111171280.15
HSP90AA1231240.04
HSP90B1190190.00
HSP90AB1221230.05
Total175191940.11
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rehman, S.u.; Nadeem, A.; Javed, M.; Hassan, F.-u.; Luo, X.; Khalid, R.B.; Liu, Q. Genomic Identification, Evolution and Sequence Analysis of the Heat-Shock Protein Gene Family in Buffalo. Genes 2020, 11, 1388. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11111388

AMA Style

Rehman Su, Nadeem A, Javed M, Hassan F-u, Luo X, Khalid RB, Liu Q. Genomic Identification, Evolution and Sequence Analysis of the Heat-Shock Protein Gene Family in Buffalo. Genes. 2020; 11(11):1388. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11111388

Chicago/Turabian Style

Rehman, Saif ur, Asif Nadeem, Maryam Javed, Faiz-ul Hassan, Xier Luo, Ruqayya Bint Khalid, and Qingyou Liu. 2020. "Genomic Identification, Evolution and Sequence Analysis of the Heat-Shock Protein Gene Family in Buffalo" Genes 11, no. 11: 1388. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11111388

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop