Determine genetic variations in heat shock factor gene family (HSFs) and study their effect on the functional and structural characterization of protein in Tali goat.

In this study, the effect of genetic variations of four heat shock transcription factor genes (HSF1, HSF2, HSF4, and HSF5) on the 3 D protein structure and function were studied. We defined the breed-specific genetic variations of pooled DNA of Tali goat that differed from the goat reference sequence (CHI2.0). Disordered regions of HSF proteins were predicted using PONDR. Post-translation changes were studied by several predicted online servers. Then, the structure of the order region of proteins was anticipated by using the Swiss model. Tali goat HSF genes contain a total number of 181, 679, 91, and 301 SNPs for HSF1, 2, 4, and 5, respectively. Also, 5 and 3 variants were identified as nsSNPs in the coding region of HSF4 and HSF5, respectively. (r.145A/S), (r.322P/Y), (r.379T/C) in HSF4 and (r.300Q/P), (r.573E/Q) in HSF5 obtained the tolerant and high confidence (SIFT score) for nsSNPs. More than half of these proteins are predicted to be disordered (56, 50, 52, and 50%, respectively for HSF1, 2, 4, and 5). Phosphorylation, acetylation, glycosylation, and Sumoylation sites of HSFs were compared between Tali goat and reference goat. Three residues S145, S263, and S322 of HSF4 in Tali goat were phosphorylation sites, and in HSF5, the reference goat has a phosphorylation site in S593.