Emerging Genetic Tools to Investigate Molecular Pathways Related to Heat Stress in Chickens: A Review

2763 animals Animals : an Open Access Journal from MDPI Animals (Basel) Multidisciplinary Digital Publishing Institute (MDPI) PMC7823582 7823582 7823582 33383690 10.3390/ani11010046 Emerging Genetic Tools to Investigate Molecular Pathways Related to Heat Stress in Chickens: A Review Perini Francesco 1 Cendron Filippo 2 * Rovelli Giacomo 1 Castellini Cesare 1 Cassandro Martino 2 Lasagna Emiliano 1 1 Department of Agricultural, Food and Environmental Sciences, University of Perugia, Borgo XX Giugno, 74, 06121 Perugia (PG), Italy; francesco.perini@studenti.unipg.it (F.P.); giacomo.rovelli@studenti.unipg.it (G.R.); cesare.castellini@unipg.it (C.C.); emiliano.lasagna@unipg.it (E.L.) 2 Department of Agronomy, Food, Natural Resources, Animals and Environment, University of Padova, Viale dell’Università, 16, 35020 Legnaro (PD), Italy; martino.cassandro@unipd.it * Correspondence: filippo.cendron@unipd.it 29 12 2020 11 1 46 46 24 1 2021 © 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ( http://creativecommons.org/licenses/by/4.0/ ). Abstract Simple Summary New genomic tools have been used as an instrument in order to assess the molecular pathway involved in heat stress resistance. Local chicken breeds have a better attitude to face heat stress. This review aims to summarize studies linked to chickens, heat stress, and heat shock protein. Abstract Chicken products are the most consumed animal-sourced foods at a global level across greatly diverse cultures, traditions, and religions. The consumption of chicken meat has increased rapidly in the past few decades and chicken meat is the main animal protein source in developing countries. Heat stress is one of the environmental factors which decreases the productive performance of poultry and meat quality. Heat stress produces the over-expression of heat shock factors and heat shock proteins in chicken tissues. Heat shock proteins regulate several molecular pathways in cells in response to stress conditions, changing the homeostasis of cells and tissues. These changes can affect the physiology of the tissue and hence the production ability of chickens. Indeed, commercial chicken strains can reach a high production level, but their body metabolism, being comparatively accelerated, has poor thermoregulation. In contrast, native backyard chickens are more adapted to the environments in which they live, with a robustness that allows them to survive and reproduce constantly. In the past few years, new molecular tools have been developed, such as RNA-Seq, Single Nucleotide Polymorphisms (SNPs), and bioinformatics approaches such as Genome-Wide Association Study (GWAS). Based on these genetic tools, many studies have detected the main pathways involved in cellular response mechanisms. In this context, it is necessary to clarify all the genetic and molecular mechanisms involved in heat stress response. Hence, this paper aims to review the ability of the new generation of genetic tools to clarify the molecular pathways associated with heat stress in chickens, offering new perspectives for the use of these findings in the animal breeding field. Keywords: biodiversity, poultry, heat shock protein, climate change status released display-pdf yes is-olf no is-manuscript no is-preprint no is-journal-matter no is-scanned no is-retracted no Received 2020 Nov 9; Accepted 2020 Dec 24; Collection date 2021 Jan. 1. Introduction Poultry, especially chicken, is one of the most reared species in the world. Indeed, chickens have a pivotal economic and ecological role in the agriculture system and, additionally, are used in the households of developed countries as the main dietary source of animal protein. The consumption of chicken products has increased rapidly in the last few decades. Moreover, the most globally consumed products of animal origin are poultry products, as they do not face any religious or cultural restrictions [ 1 ]. Indeed, poultry is the fastest growing animal product, especially in developing countries. Chicken products have conquered the meat market, as they are commonly affordable and low in fat [ 2 , 3 ]. Data from 2018 showed that chicken products had dominated the market in both the egg branch and the meat industry, with 114 million tons being produced [ 4 ]. For these reasons, it is appropriate to assume that there will be an expansion of the whole poultry branch, as demand for meat and eggs is led by rising human population and urbanization. In these circumstances, the sector has to withstand unparalleled challenges [ 5 ]. In this scenario of expected growing demand [ 6 ], livestock production is likely to be undermined by climate change, competition for land and water, and food security issues [ 7 , 8 ]. In particular, climate change is involved in the reduction in pasture and water availability, the onset of new diseases, increasing numbers of drought seasons, and the occurrence of exceptional natural phenomena, which all reduce the quantity and quality of meat and eggs [ 9 , 10 , 11 , 12 ]. The alteration of environmental factors such as sunlight, temperature, and humidity affect animal metabolism and mechanisms of thermo-regulation and can cause imbalances in animal physiology [ 13 ]. Animals have proved to be adaptable in order to overcome adverse climatic conditions [ 14 ], and this permits them to survive in harsh climatic conditions [ 15 ]. Several mechanisms are the basis of this adaption ability, namely, morphological, anatomical, behavioral, physiological, biochemical, cellular, and molecular characteristics [ 16 ]. Broilers have reached high productive levels due to genetic improvements; however, their lack of effective thermoregulation and heightened metabolism entail a worse adaptation to harsh environments [ 17 ]. In contrast, native chickens generally have genetic and morphological features and rusticity that help them to face harsh conditions. Native chickens have a greater capacity to cope with heat stress [ 18 ]; however, they give a lower rate of production, since they have never been subject to breeding programs. The aim of the present review is to summarize papers concerning the identification of heat stress (HS) genetic markers in chickens using new molecular tools. In fact, both genotypic and phenotypic traits may be useful indicators of heat stress (HS) susceptibility or tolerance in different chicken breeds. Nowadays, the identification and the inclusion of appropriate biomarkers in breeding programs for HS responses in chickens could be a crucial aspect in developing thermo-tolerant breeds using marker-assisted selection [ 19 ]. 2. Climate Change and HS in Chickens: An Overview Climate change means spatial and temporal variations in the environmental climatic parameters on Earth [ 20 ]. These climatic parameters include temperature; humidity; solar radiation; precipitation; and the temperature of water, which causes glaciers to melt [ 21 ]. The Intergovernmental Panel on Climate Change’s Fifth Assessment Report shed light on a future temperature increase in the Earth’s surface ranging from 0.3 to 4.8 °C by 2100 [ 22 ]; this may affect livestock through variations in feed crops and forage, biodiversity, animal welfare, productive, and reproductive traits. Figure 1 shows climate change in the last century on the global temperature level. Figure 1 Temperature change (°C) in 2019 compared with the year 1919. Note: grey areas signify missing data (source: figure taken from https://data.giss.nasa.gov/gistemp/ and adapted for illustrative purpose only). The change that is occurring in the climate is mainly related to the constant increase in the level of greenhouse gasses (GHGs) [ 23 ]. All over the world, climate change is driving temperature rises, altered photoperiods, and reductions in precipitation, causing reduced feed quality and quantity [ 24 ]. Accordingly, climate change is one issue that the livestock sector will have to cope with in the following years. Above all, HS is the worst environmental stressor for poultry production worldwide. As early as 2003, a study estimated the annual economic loss due to HS in the US livestock industry to range from $1.69 to $2.36 billion; from these data, $128 to $165 million was related to the poultry sector [ 25 , 26 ]. Moreover, as Figure 2 shows, chicken production is mainly distributed in tropical areas, and hence, is directly affected by hot temperatures. Figure 2 These figures represent the global distribution of chicken in 2010 in intensive ( a ) and extensive ( b ) systems (source: figure taken from https://dataverse.harvard.edu/dataset.xhtml?persistentId = doi:10.7910/DVN/A7GQXG and adapted for illustrative purpose only). Exposure to HS can decrease growth rate and feed efficiency, alter immune response, damage gut microflora, and finally decrease meat quality [ 27 ]. Chickens are homoeothermic animals, able to maintain a relatively constant inner temperature, even if only in an environmental temperature of 21–28 °C. Body metabolism is the main aspect responsible for heat production, and needs heat removal to maintain the internal body temperature at 41.4–42.9 °C [ 28 ]. Indeed, hot temperatures can ruin the maintenance of this homoeothermic due to the physiological changes that affect body temperature stability. The deterioration of the regulatory mechanism of homoeothermic balance has direct and indirect consequences on chicken production. As reported in the review of Zaboli et al. [ 27 ], the negative effect of HS on meat quality is due to the combined effects of anaerobic glycolysis acceleration (metabolic acidosis) by panting [ 29 ], the increase in Reactive Oxygen Species (ROS) [ 30 ], and the corticosterone concentration in the blood [ 31 , 32 ]. HS occurs when the amount of heat produced by the animal surpasses its capacity to dissipate the extra heat to the surrounding environment [ 33 ]. Heat swap between the animal body and environment represents the major system to preserve homoeothermic balance. Heat generation, maintenance, and dissipation are ruled by biological and physical components. The intricacy of the phenomena requires a thermal index that could represent the environmental pressure on the biological aspects involved. The most common index by which to measure the impact of environment on poultry livestock is given by the Temperature Humidity Index (THI) [ 34 ]. The THI is a useful instrument with which to measure livestock productivity response as a function of climate [ 35 ]. THI takes into consideration the air temperature and relative humidity of the environment, and how they impact on diverse species [ 36 ], and it is an indicative measure of the sum of forces external to the animal that act to displace body temperature from its set point [ 37 ]. During the period of heat stress, chickens can employ some physiological heat exchange mechanisms in order to overcome the stressed status. First of all, birds change their behavior in HS conditions: less time is devoted to walking and standing [ 38 ], the consumption of feed is less, and the consumption of water is greater [ 39 ]. In an anatomical way, chickens are able to dissipate heat by radiation and conduction through an increase in blood flow to the skin surface. Moreover, partial feather loss is not uncommon in hens, mainly from the neck, back, and breast regions. Furthermore, painting is one of the visible response of poultry during exposure to heat. This specialized form of respiration dissipates heat by evaporative cooling at the surfaces of the mouth and respiratory passageways [ 40 ]. 3. The Study of Heat Adaptation: New Perspectives through Genomic Tools Over the past few decades, breeding practices have taken into consideration the estimated breeding values from a phenotypic and pedigree database, which has surely enhanced specialized chicken breeds in terms of production rate [ 41 , 42 ]. However, these approaches were limited according to the slow upgrade in chicken breeding and difficulty in splitting preferred from undesired traits [ 43 ]. To overcome these issues, Quantitative Traits Loci (QTLs) have been mapped based on both microsatellite and Single Nucleotide Polymorphism (SNP) marker technology in order to associate favorable phenotypic traits with genomic regions. Unfortunately, the QTL accuracy was low and the typical confidence range of the QTL usually covers a broad region or even a whole chromosome (CHR) [ 44 ]. Therefore, it was improbable that QTL mapping could highlight the genomic pathways related to the target traits in chicken [ 45 ]. Many QTLs have been associated with a phenotype of interest, such as body temperature ( Table S1 ). The managing of new traits, such as HS resistance, has brought about the development of new technologies to investigate genetic relations with these phenotypes. A good approach to identify the functional genes and polymorphisms correlated with HS resilience has been accomplished by virtue of high-throughput screening technologies, such as genome-wide analyses of genetic variations [ 46 , 47 ]. The Genome-Wide Association Study (GWAS) was developed to shed light on the genetic background of phenotypic traits. GWAS is based on the hypothesis of links between genetic variants and different phenotypic traits [ 48 ]. This approach is favored over other genomic tools for its high power in detecting polymorphisms and mutations and its capability to find exact genomic regions [ 49 ]. This method has been used for several domestic animals, and considers Single Nucleotide Polymorphisms (SNPs) and whole-genome sequencing (WGS) [ 44 ]. GWAS is a powerful approach with which to identify candidate genes, genomic regions, and polymorphisms related to specific phenotypic features. Its application has become increasingly widespread because it is useful for more accurate breeding programs for chickens and other animals [ 50 , 51 ]. Another high-throughput screening technology largely used is RNA sequencing (RNA-Seq). RNA molecules carry the transcribed information encoded in selected genes, which can be translated into proteins or used to check gene expression directly (i.e., miRNA) or indirectly (i.e., heterogeneous nuclear ribonucleoproteins—hnRNPs). Hence, the rate of RNAs expressed in a known circumstance shows the cell state and can disclose stress mechanisms. RNA-Seq is used in the determination of a specific trait thought the comparison of differential gene expression among different tissues [ 52 ]. Currently, RNA-Seq is the most suitable method for studying gene expression and identifying novel RNA species. Differently from the DNA microarray method (the previously used technique), RNA-Seq results in higher resolution data and a greater dynamic range of detection [ 53 ]. Moreover, RNA-Seq is a direct method with which to obtain information about sequence identity, which is important for the study of unknown genes and novel transcript isoforms. The greatest advantage of RNA-Seq, its high resolution, has promoted transcriptomic research, generating an enormous quantity of data [ 52 ]. 4. Limitation of High-Throughput Screening Technologies In spite of the general success of GWAS methodology, important limitations and failures are still present [ 54 ]. Since WGS is still expensive for mass application, genomic studies are mostly based on SNP arrays, in which only a small number of polymorphisms are drafted to be highly polymorphic and are used to represent the whole genome. Furthermore, it has been suggested that many polymorphisms associated to quantitative traits are rare variants and could be missed in low-density SNP arrays, but we expect to find them in WGS data [ 55 ]. While WGS has a greater depth of analysis, at the same time the expensive cost and the need to manage a huge amount of output data make WGS a methodology that is difficult to pursue and not within everyone’s reach. First, the identification of alleles linked to a specific traits or phenotypes requires huge sample sizes (100 to 1000) for the adequate statistical power of GWAS [ 56 ]. Hence, a lot of care is needed in all the steps of the experiment (i.e., sample collection, genotyping, data pruning, software analysis, and interpretation of results) for protecting against erroneous findings. Another difficulty of GWAS is the interpretation of individual associated variants. Even though many variants are assayed, these represent only a partial slice compared to the total genetic variation. A given gene might have a greater number of polymorphisms than those reported in the chip panel usually used for GWAS studies. The disadvantage at this stage is in representing the real quantity of SNPs in the exons of the genes (hence, the coding part of the genes) [ 57 ]. In contrast, genetic variants fall within the introns, which could affect gene expression, alternative splicing, and DNA methylation, which represent a large portion of the associated loci and may be located at various physical distances from their targets. In this context, it is really hard to evaluate how genetic variants can affect the related phenotypic trait [ 58 ]. Some limitations are also highlighted in RNA-Seq technology: when an RNA-Seq assay is running, the greater part of the transcript (around 70%) may be related to a very small number of genes [ 59 ]. Although investigating gene expression levels is relatively easy, the analysis of allelic imbalance or splice variation is more difficult. This trouble is caused by using a fractional dataset from a given gene, coming from the SNPs already known or the exon/exon junction of interest, which enhances the risk of the poor statistical power of the method. To overcome these issues, an increase in the depth of reads is required, which can bring to light the transcriptome of poorly expressed genes. Unfortunately, these methods are hard to undertake in economic terms, but there are several ways to alleviate the problem. A stratagem may be to precipitate the not-target transcripts using probes when the assay has been set on specific transcripts [ 60 ], in order to ensure that the RNA-Seq is mostly performed on targets. One of the other issues is the read length. In fact, when the reads are shorter than the target, alignment and reassembly are required. In sequencing the whole transcript in a long read, the accurate sequence would be pinpointed; in this case, the analysis of allele-specific expression will be more accurate and mutually exclusive exons in transcripts will be highlighted [ 61 ]. Lastly, the issue of RNA-Seq is the validation of the experiment; indeed, it is largely known that genes characterized as divergent expressed by an RNA-Seq experiment need to be validated by real-time PCR. This additional step increases the time and cost of the assay [ 62 ]. 5. Phenotypes against HS Poultry, like other animals, has the possibility to develop specific phenotypes advantageous for adaption to the harsh environment where they live. There are several phenotypes that mainly act for the alleviation of HS, mostly related to feather types [ 63 ]. Indeed, feathers guarantee a thermal shelter between the animal body and the environment. Plumage delays the process of heat elimination from the skin surface [ 64 ]. As an example, the Naked Neck (Na) chicken shows a better fitness under HS conditions [ 65 ]. The Na gene is a single dominant autosomal gene that allows decreasing feathers in the neck region, which permits a better dissipation of heat [ 66 ]. The Na gene reduces the plumage cover by 20% and 40% in heterozygous and homozygous Na, respectively [ 67 ]. In broilers, it is correlated with an increase in breast muscle and body weight [ 68 , 69 ], a reduction in abdominal fat [ 70 ], and a reduction in body temperature [ 71 ]. Moreover, laying Na chickens under hot temperatures showed an improvement in egg mass, deposition rate, and egg characteristics [ 72 ]. In a recent study, Galal et al. [ 73 ] compared the expression levels of HSP70 (used to assess the heat tolerance) in three Egyptian local breeds (Fayoumi, Dandarawi, and Sinai) with and without the Na gene and under normal and HS conditions. As expected, they found higher HSP70 expression levels in crossbreeds with the Na gene, suggesting that the Na gene is responsible for the up-regulation of HSP70 expression and has a positive impact on HS adaptation not only by reducing feather cover [ 73 ]. Another important feature of the plumage is the color of feathers, which can impact the ability of chickens to respond under HS conditions. A recent study demonstrated that dark chicks showed a lower expression rate of genes belonging to pathways of stress (cellular stress: SOD2 and HSPA8 ; DNA damage repair: ALKBH3 ) than paler chicks [ 74 ]. This happens because the plumage element reduced solar heat gain by 5% in both light and dark plumages. However, because the overall levels of solar heat gain were greater in dark versus light plumages, there were different fractional decreases in the heat load for light (41%) versus dark (25%) plumages [ 64 ]. Lastly, Jiang et al. highlighted the importance of density in regard to plumage. Indeed, the density of contour feathers was significantly correlated to heat tolerance under acute HS, indicating that it could be regarded as a phenotypic marker for heat tolerance in chickens [ 75 ]. 6. Physiological Biomarkers for HS Detection Recording animal performance under HS is a method to measure an animal’s capacity to face heat. The most recorded parameters are body temperature, respiration rate, and the blood level of cortisol, which could be useful as discriminant criteria in order to identify animals that are tolerant to HS conditions [ 76 , 77 ]. About 50 years ago, body temperature was suggested as the most valuable HS indicator [ 78 ]. Body temperature is assessed through Rectal Temperature (ReT). The range of physiological ReT in chickens is between 40.6 and 43 °C. In 2013, Chen et al. established a precise method to evaluate HS resilience through ReT. The correlation between ReT and HS survival time was closely related at 18 h of heat exposure. Therefore, in the choice of heat resistance, it is recommended to evaluate the ReT after 18 h HS and calculate the ΔT18 (interval between the ReT of 18 h and the initial ReT) and standard deviation (SD). Animals with less than the mean value −1/2 SD could be selected as heat tolerant [ 79 ]. Unfortunately, Van Goor et al. reported very low hereditability values for many QTLs associated with body temperature [ 80 ]. HS in chickens is regulated by a huge number of organized reactions that are difficult to measure in a few in vivo biochemical or biophysical indices over a short time period. Moreover, it has been demonstrated that chicken body weight and genotype could affect the ReT [ 81 , 82 ]. Nowadays, one of the common methods used to identify ongoing stress status in chickens is the measurement of the corticosterone level in the blood [ 83 , 84 , 85 , 86 ]. Corticosterone affects stress-induced responsiveness and can protect or destroy the coping ability of the organism [ 87 ]. High levels of corticosterone occur in chickens during HS as a defense mechanism [ 88 ]. Currently, corticosterone represents the “gold standard” of stress markers in poultry, but its reliability and validity have recently been questioned. Indeed, the assay for its identification can give a false positive due to the circadian rhythm and cross-reactivity with other glucocorticoids [ 89 ]. 7. Modulation of Inflammation and Immunity during HS In poultry, the primary production and differentiation of immune cells and antibodies occurs in immune organs such as thymus, the bursa of Fabricius, and the spleen. HS causes oxidative damage to the cell membranes in immune organs [ 90 ]; many studies have reported significant reductions in the weights of the immune organs due to HS, including the thymus, bursa, and spleen [ 91 , 92 , 93 ]. As a consequence of HS, Tang and Chen [ 94 ] detected significant reductions in B and T lymphocytes, which led to decreases in the production of antibodies, changes in cytokine secretion, and lower numbers of macrophages with a reduced phagocytic ability [ 95 , 96 ]. In 2014, Xu et al.’s study of the chicken spleen detected an increase in the mRNA levels of Tumor Necrosis Factor alpha (TNF-α) and Interleukin-4 ( IL-4 ), whereas the levels of Interferon gamma ( IFN-γ ) and Interleukin-2 ( IL-2 ) were lowered [ 97 ]. Simultaneously, HS triggers heat shock factors and proteins aimed in cellular protection [ 98 ]. It is known that fast-responding genes (i.e., pro-inflammatory genes) undergo an inhibition regulated by feedback mechanisms. For instance, Interleukin-6 ( IL-6 ), the most pyrogenic cytokine, was not increased by acute exposition to heat. This is because Heat Shock Factor 1 ( HSF1 ) down-regulates the IL-6 gene expression [ 99 ]. Han et al., in 2010, demonstrated that the response to acute HS in Ross broiler chickens is mediated by a higher activity of IL-2 in lymphocytes of the spleen [ 100 ]. Another study shows how acute exposure to HS enriched the mRNA expression of IL-4 in the spleen of HS chickens [ 101 ]. Quinteiro-Filho et al. highlighted that, in broilers exposed to chronic HS, the plasma levels of IgA and IgG were reduced, consistent with the lower robustness rate of broiler breeds regarding thermal stress [ 102 ]. Finally, Interleukin-17 ( IL-17 ) expressed the highest up-regulatory responses to acute HS. Indeed, IL-17 plays an immunological role in mucosal response and is an activator of the TLR4-dependent means of clearance of pathogenic bacteria [ 103 ]. IL-17 also contributes to the development of adaptive immunity to inflammatory agents [ 104 ]. 8. The Heat Shock Proteins (HSPs) Family HSPs are characterized by their ability to be induced by HS and, in molecular terms, by the presence of a functional heat shock element in their promoter [ 105 ]. HSPs ( Figure 3 ) act as chaperones that bind other proteins, after HS damage, with the purpose of preserving their structure, managing the proteins’ migration across membranes or organelles, or keeping the receptor availability or specific enzyme functionality under control [ 106 ]. Figure 3 Schematic mechanism of heat shock protein (HSP) in cellular response against heat stress (HS). HSF = Heat-Shock Factor; HSE = Heat-Shock Element; DP = Denatured Protein; P = Phosphate. The expression of the HSP could be physiological or induced by different inputs such as HS. HSPs might be classified in several manners; here, we report the most used classification based on molecular weights: HSP100 , HSP90 , HSP70 , HSP60 , small HSPs, and chaperonins [ 107 ]. Aside from the molecular chaperone activity, small HSPs are also involved in other cellular pathways such as stress tolerance, protein folding, cytoskeletal integrity, and cell cycle [ 108 ]. In most cases, HSPs are found in the cytoplasm, but it is also possible to find their presence in the extracellular matrix, in which they could serve as stress alarms and stimulate immune cells [ 107 ]. Table 1 briefly summarizes the most studied HSPs and HSFs, the biological processes in which they are involved, and their molecular function in chickens. Table 1 Using the Panther Classification System, the gene ontology of the most important genes related to heat stress was performed on two domains (molecular function and biological process). The current red jungle fowl ( Gallus gallus ) genome assembly, GRCg6a. Gene Name Mapped Ids Family Name Molecular Function Biological Process

HSP70

HSPA2 NC_006092.5 (GRCg6a) Heat shock 70 kDa protein 2 ATP binding; heat shock protein binding; unfolded protein binding; ATPase activity; Vesicle-mediated transport; Chaperone-mediated protein folding; Cellular response to unfolded protein

HSP90B1 NC_006088.5 (GRCg6a) Heat shock protein 90 beta family member 1 Unfolded protein binding Protein folding

HSF1 NC_006089.5 (GRCg6a) Heat shock factor protein 1 DNA-binding transcription factor activity; RNA polymerase II proximal promoter; sequence-specific DNA binding Cellular response to heat; Regulation of transcription from RNA polymerase II promoter in response to stress; Transcription by RNA polymerase II; Positive regulation of transcription by RNA polymerase II

HSF3 NC_006091.5 (GRCg6a) Heat shock factor protein 3 DNA-binding transcription factor activity; RNA polymerase II proximal promoter sequence-specific DNA binding Cellular response to heat; Regulation of transcription from RNA polymerase II promoter in response to stress; Transcription by RNA polymerase II; Positive regulation of transcription by RNA polymerase II

HSPD1

HSP60 NC_006094.5 (GRCg6a) Heat shock protein family D (Hsp60) member 1 Unfolded protein binding Protein folding

SERPINH1

HSP47 NC_006088.5 (GRCg6a) Serpin family H member 1 Endopeptidase inhibitor activity Serine-type endopeptidase activity Protease binding Proteolysis; Cellular protein metabolic process; Negative regulation of endopeptidase activity 9. HSP in Commercial Chickens It is universally known that several stress stimuli, and HS in particular, can influence HSPs’ expression in various chicken tissues. For instance, the quantity of HSP70 and ubiquitin transcripts increases when testis cells are exposed to high temperatures [ 109 ]. In the same way, a relevant increment in the HSP70 expression in female broiler chicken brains after four days of heat treatment was noticed [ 110 ]. Thermal stress caused the induction of HSP90α and HSP90ß in chicken heart, liver, and spleen, but the HSP90α and HSP90ß mRNA levels were stable in the brain [ 111 ]. The gain expression of HSP in the heart may allow protection in harsh environments. For instance, in the heart tissue of HS broiler HP60 , HSP70 , and HSP90 proteins and their relative mRNAs there is a gain after 2 h of HS, but this diminishes quickly with chronic HS [ 112 ]. Indeed, the great diversity in heat shock response among different tissues and different broiler lines is widely known. One study reports that in breast muscle, acute and chronic HS enhanced protein oxidation, but the HSP gene expression remained at physiological levels and the only trend increases were observed in the gene expression of HSP70 and HSP90 after acute HS [ 113 ]. Thermal stress up-regulates the expression of the HSP70 gene in two fast-growing broiler strains. It is interesting to note that the expression of the liver HSP70 gene in heat-stressed Ross broilers was significantly higher than that reported in Cobb [ 114 ]. Moreover, Tang et al. reported that two groups of broilers, both under heat stress and one treated with aspirin, showed different responses. The aspirin-treated group showed a greater response to HS in the kidney, with less damage to the tissue, due to the ability of aspirin to induce high expression levels of HSP47 and HSP60 [ 115 ]. Liu et al. [ 92 ] highlighted that HS stimulus is able to raise the expression of HSP27 (small HSP), HSP70 , and HSP90 mRNA in the bursa of Fabricius and the spleen of broilers; these data corroborate with the data on stressed broiler spleen from Slawinska et al.: HSP25 , HSP70 , and HSP90 were up-regulated [ 101 ]. Quite the opposite happened in the same heat condition in the thymus: HSP27 and HSP90 mRNA were significantly reduced. In addition, an HSP70 up-regulation was observed in the serum of heated layers, compared to layers under normal temperature. In some layers, HS substantially altered the gut microbiome [ 116 ]. 10. HSP in Local Chicken Breeds The genomic tools already discussed have also been applied in several studies in order to verify the HS resistance of local chickens, compared to commercial lines. The aim of this paragraph is to clarify the pathways that allow local chickens to have a greater resistance ability to HS. Indeed, Cedraz et al. reports the results of a real-time PCR that show the expression levels of several HSPs in the muscle of two local Brazilian chickens (Peloco and Caneluda breeds) and a broiler (Cobb 500) [ 117 ]. Thus, under HS the HSP70 and HSP90 genes were highly expressed in all genetic groups compared to the control conditions. At the same time, in HS conditions the HSP70 and HSP90 genes were significantly more expressed in both local chickens than in the broiler [ 117 ]. A microarray-based study highlighted the higher expression of HSP70 , HSP90AA1 , and HSP25 in the testes of a Taiwan egg-type chicken in an HS environment [ 118 ]. Moreover, the same authors compared the obtained results with the data from Taiwan meat-type chickens submitted to the same experimental design. The results of the study showed substantial changes in the expression of 28 genes in both types of chickens, and most of these genes were HSP genes. However, the testicular responses to acute HS differed between the two types of chickens: for instance, the HSP60 gene was highly expressed in HS meat-type chickens; on the contrary, the BAG3 gene was more expressed in egg-type chickens [ 119 ]. BAG3 and HSP60 genes both negatively regulate cellular apoptosis; hence, it is very interesting to see how two different chicken breeds have many genes in common, but differ for others when they have to face the same stressed situation. In a recent study by Sharma et al., the HSPß1 mRNA expression (a member of small HSP family that helps in maintaining the homeostasis of proteins by stabilizing different types of non-native proteins) was evaluated in Punjab Red and Rhode Island Red layers under HS. The results demonstrated a direct association between HS, HSPß1 mRNA expression, and serum HSPß1 concentration [ 120 ]. In 2020, Radwan described the possibility of using gene expression data to explain why some traits underwent genetic selection. In fact, the study shed light on improvements in some traits in a native Egyptian chicken line (Fayoumi) exposed to thermal stress, such as heat resistance, by selecting for the desired traits. In the study, Fayoumi chickens were reared either in a normal or a heated environment [ 121 ]. Thirty-five-week-old females from the HS group with the best egg production and the strongest eggshells were mated with their male siblings. The same mating program was used with F1 birds to obtain a second generation [ 119 ]. Radwan also evaluated the uterus mRNA level of HSP90 and ovocleidin 17 ( OC-17 ), which were demonstrated to be involved in eggshell strength [ 121 , 122 ]. The results from the real-time PCR show increased levels of HSP90 and OC-17 in the two successive generations, meaning that both eggshell strength and heat tolerance were improved thanks to selection in birds raised under conditions of HS [ 121 ]. It is noteworthy that some studies used the statistical association between SNPs and thermo-tolerance; indeed, in Zhang et al. the SNP site (C.1388 A>G) of the chicken Heat Shock Factor 3 ( HSF3 ) gene was associated with heat resistance in two chicken lines [ 123 ]. Chen et al. showed that the SNP of C.141 G>A in the HSP90β gene in chicken had an impact on the HS resistance traits, and the GG genotype represents the most suitable one under HS. The HSP90β mRNA expression was shown to be tissue-dependent by qRT-PCR. The expression of HSP90β mRNA in the heart, liver, brain, and spleen of Lingshan chickens was significantly superior compared to that of White Recessive Rock, underlining the fact that the Lingshan breed is more adapted to tropical temperatures [ 124 ]. A different study in White Recessive Rock chickens was carried out by Kong et al. and highlighted that the C.744C>G SNP in the 5’-flanking region of the Glucose-Regulated Protein 78 gene ( GRP78 ) was significantly correlated with heat tolerance parameters [ 125 ]. The gene GRP78 belongs to the HSP70 family, is a fundamental chaperon in various animals, works against apoptosis, and allows proteins and organelles to preserve their physiological functionality [ 126 ]. Finally, in Kong’s study, a qRT-PCR assay indicated a gain in the GRP78 mRNA expression in all tissues, which then decreased with chronic HS and peaked at 3 h after HS. In conclusion, the GRP78 mRNA expression depends on time and tissue [ 125 ]. In the same way, Irivboje et al. evaluated genetic polymorphisms in intron 7 and exon 8 of the HSP90AA1 gene in two exotic Nigerian chicken strains (Brown dominant and Fowl Hyline brown). Several SNPs (A7T, A160T, T223A, and C134T) were detected; however, the SNP A7T, qualified for the association analysis after the Hardy-Weinberg equilibrium test, was not significantly linked to heat resistance features [ 127 ]. An important study by Fleming et al. compared several chickens from different breeds and climates of Africa and Northern Europe for a selection signature evaluation; the aim of this work was to understand the adaptation mechanisms of the animals to their local environments. The African chickens were distinguished for their stronger resistance toward stress signaling and angiogenesis, while the Northern European chickens showed more genetic selection toward processes related to energy homeostasis. In Chromosomes 2 and 27, the populations were the most divergent in terms of selection pressure. Moreover, chromosome 27 was involved in heat tolerance in African chickens, while novel insights into unique genomic regions on chromosome 2 could be related to development and environment for Northern European chickens. In conclusion, the study shows in an excellent way how two populations subjected to divergent environments present in their genome the signs of different selection affecting different metabolic pathways [ 128 ]. A study in support of the central role of HSPs in facing HS is the recent research of Srikanth et al. Using RNA-Seq methods (in liver and heart tissues), two indigenous chicken ecotypes from Kenya were compared; the animals were sampled from the tropical climate in Mombasa (lowland) and the colder Naivasha (highland) regions in order to investigate the effects of acute (5 h, 35 °C) and chronic (3 days of 35 °C for 8 h/day) HS. Only four different differentially expressed genes (DEGs) were found in all four experimental groups and were identified as HSP70 family member 8 (HSPA8) both in acute and chronic HS and small HSP family member 7 (HSPB7) in acute HS [ 129 ]. Moreover, Fatty Acid–Binding Protein 4 ( FABP4 ) was found differentially expressed in chronic HS [ 129 ] and also in the hypothalamus of a broiler chicken [ 130 ]. The FABP4 gene is involved in the Peroxisome Proliferator-Activated Receptor ( PPAR ) signaling pathway, which is required for energy metabolism and regulating the oxidative stress-induced inflammatory response [ 129 , 131 ]. The enrichment of the PPAR pathway indicates that during adverse conditions, such as HS, PPAR promotes adaptation events [ 132 ]. An interesting study by Wells et al. reports a genome-wide SNP scan using featherless chickens; these chickens carry a single recessive mutation that causes a lack of almost all body feathers, as well as foot scales and spurs, due to a failure of skin patterning during embryogenesis [ 133 , 134 ]. The trait is potentially useful in tropical agriculture due to the ability of featherless chickens to better tolerate heat [ 135 ]. Through a cost-effective and labor-efficient SNP array mapping approach, Wells and his groups showed that a nonsense mutation in FGF20 , related to a loss of a highly conserved region, was completely associated with the featherless phenotype. This aspect was confirmed by in situ hybridization and quantitative RT-PCR assays that revealed a high FGF20 expression during the early stages of feather placode patterning. Hence, the loss of genetic function highlights the role of FGF ligand signaling in feather development and suggests FGF20 as a novel central player in the development of vertebrate skin; thus, this gene could be considered in crossbreeding in order to obtain featherless chicken lines that are less susceptible to high temperatures [ 133 ]. Recently, Te Pas et al. studied the ability of local chickens to rapidly adapt to a hotter environment. Indeed, they analyzed a transcriptome from different tissues (heart, breast muscle, and spleen) in Ethiopian lowland chicken (adapted to lowland thermal conditions) and highland chickens (non-adapted to lowland environments) in relation to the changes in temperature in lowland conditions during the day in the morning, noon, and evening [ 136 ]. The highland chickens responded rapidly to HS, and muscle tissue had more HSP regulation in highland chickens, suggesting that muscle tissue is particularly vulnerable to HS more than the other types [ 136 ]. Interestingly, through a statistical analysis of RNA-Seq data, evidence of epigenetic mechanisms was found. Indeed, epigenomic regulations on chromatin re-modeling through both DNA methylation and histone (de)acetylation mechanisms were reported. These biological mechanisms are fast regulators, much more so than mutation and selection [ 137 ]. 11. Epigenetic Mechanisms in Chicken HS Response Epigenetic mechanisms are one of the main options used by the body to better cope with adaptive developmental reprogramming [ 138 ]. In fact, the following definition clarifies how epigenetics work: “an epigenetic trait is a stably heritable phenotype resulting from changes in a chromosome without alterations in the DNA sequence.” The above-proposed description of epigenetics can involve the phenotype heritability through either mitosis or meiosis [ 139 ]. Epigenetic adaptation assumes that environmental factors—e.g., environmental temperature—can affect the physiological control system during the development phase, which might make modifications in the thermoregulatory process [ 140 ]. Genetic investigations in animals revealed that epigenetic markers may be passed down across generations, leaving a mark on the offspring phenotype [ 141 ]. Methylation is one of the most considerable epigenetic modification processes, occurring at the DNA level and playing an upstream role in the control of gene expression [ 142 ]. A selective pressure of particular stimuli can flow into the modification of the methylation level in the organism. It is known that the methylation level and target can lead to the regulation of the expression of tissue-specific genes [ 143 ]. Vinoth et al. showed the DNA methylation patterns of the HSP promoter region in poultry brain tissue in order to assess the number of methylated genes managing thermal adaptation. The authors used two chicken lines embryos (Naked Neck and Punjab Broiler-2) subjected to HS and normal temperature (named Control (C) and thermal-conditioned (TC) embryonic groups, respectively. Chickens developed from these TC and C embryos have been successively subjected to HS and normal temperature, establishing four groups: control normal (CN), control heat exposure (CHE), thermal conditioned normal (TCN), and thermal conditioned heat exposure (TCHE). In the brains of 17-day-old embryos, the HSP gene expression was up-regulated in the TC group of both Naked Neck and Punjab Broiler-2, except for HSP60 ; this implies that the embryonic HS has been stressful for both the chicken breeds. Moreover, the HSP mRNA expression in the brain tissue was lower in TCHE chicks, as a proof of the valuable role played by embryonic heat treatment. CHE chicks exhibited the highest values of all HSP genes. Comparing the expression of HSP between the CHE and TCHE groups, the authors reported an overt heat-stressed state of the CHE chicken. These differences in expression can be related to the different methylation levels in the promoter of the same HSP genes. Indeed, a higher methylation level of HSP90 α and β and HSP70 in TCHE, compared to CHE, was found [ 144 ]. These results clearly demonstrated how high levels of methylation in the promoter serves as a down-regulator of gene expression. This is in accordance with the study of Gan et al., where in chicken muscle subjected to HS an opposite association between the expression and promoter methylation of HSP70 was reported [ 145 ]. Another study, consistent with the findings of Vinoth et al., described a reduction in HSP ( 70 and 27 ) mRNA expression in different tissues (heart, liver, muscle, spleen, and bursa) in chickens exposed to HS during the late incubation period of embryos [ 146 ]. These results corroborate the thesis that a post-hatch heat stress exposure can increase the heat adaptation of chicks. 12. A Step over HSPs In the last five years, different approaches to give a more exhaustive picture of HS in chickens have been established. The subsequent studies have used different statistical analyses and different biological starting points and have obtained varied results, but have maintained the main conductive threads, with the tissue under study being the liver. The liver has become the very first organ target because it represents the metabolic center of the organism. An interesting approach has been used by Jastrebski et al. (2017), who merged the RNA-Seq gene expression approach with the study of metabolome in chronic stressed chickens [ 147 ]. In fact, one of the main functions of the liver is to maintain the homeostasis of lipids, sugars, and amino acids in the body. Generally, metabolome and transcriptome data have shown a strong response in the liver. Indeed, many pathways and metabolites had been found to be activated in order to maintain homeostasis and react to issues caused by oxidative stress. In more detail, glycogenolysis and gluconeogenesis (as also shown by Kumar et al. [ 148 ]) as well as fat deposition, glycosylation, and glutathione production were increased. Finally, the transcriptome data have shown a trend to slow down the cell cycle in order to allow time for repairing DNA damage caused by HS [ 147 ] (Jastrebski et al., 2017). A change in liver biology has also been noticed by other studies. For instance, Li et al. in 2011 found four new genes involved in the response to HS: PM20 / PM21 , ASB2 , USP45 , and TFG . In particular, the last one was involved in the activation and enrichment of two main molecular pathways, mitogen-associated protein kinase (MAPK) and nuclear factor kappa-light-chain-enhancer of activated B cells (NFKB), in the response of broilers to HS [ 149 ]. The results were also corroborated by Bertocchi et al. [ 150 ] and especially Coble et al. through RNA-Seq, also showing novel biological processes that were functionally enhanced in HS response: cellular signaling, the endocrine system, and molecular transport. Particularly, the expression patterns of the CCK , TRPC5 , DIO2 , and DIO3 genes display a direct link between feed intake, deiodinase activity, and temperature regulation in response to heat [ 151 ]. Walugembe et al. performed a GWAS in different ecotypes of chickens from Brazil, Egypt, and Sri Lanka that have to face HS in their natural environment. In particular, Sri Lanka has hot humid climatic conditions that, besides being favorable for the pathological infection of livestock, also present challenging conditions such as HS. Through selection signal methods, the authors found two regions under selection on chromosome 4 across the Sri Lanka ecotypes. This chromosomic region codes for Toll like receptor 3 ( TLR3 ) and Nuclear factor kappa B subunit 1 ( NFKB1 ). TLR3 usually activates the MAPK and NFKB pathways [ 152 ]. Additionally, Sun et al., in 2015, studied the phenomenon of HS in chickens, but used RNA-Seq to identify heat stress-responsive genes in the chicken male white leghorn hepatocellular cell line. Multiple biological processes were affected by the responsive genes, including translation, transcription, chromatin modification, DNA repair, and DNA synthesis. In addition, two signaling pathways were modulated by HS: TGFß-enriched and WNT down-regulated [ 153 ]. 13. Future Perspectives and Conclusions Due to the enormous improvement of genetic tools, data generation is becoming affordable and effortless, resulting in huge amount of information. In the last 20 years, data should have elucidated the biological system that underpins animal production, health, and welfare traits. This has shed light on the mechanisms and detection of (potential) biomarkers and improved animal breeding strategies. Several researches have demonstrated that RNA-Seq is the most suitable and promising technology that can give a definitive view of the pathways involved and the role of every single gene in the thermoregulation process. Only after the entire process is elucidated can strong biomarkers be detected. Then, the allelic variants codifying for favorable biomarkers could also be selected in specialized breeds by crossbreeding schemes, including indigenous breeds, or by biotechnological approaches (i.e., gene editing). With the threat of climate change on animal production, favorable alleles for HS reduction should be introduced in breeding plans. The climatic problem affects both intensive and extensive productive systems; the latter should be more focused on local chicken breeds, which are able to produce in harsh environments. In this context of change, a fundamental role can belong to local breeds—those not selected for commercial production, but that have genetically evolved over time thanks to the selective pressure of the environment in which they have lived. Thus, local breeds could represent an important genetic reservoir of phenotypes more adapted to HS. In this way, according to an increased demand for “ethical food,” biodiversity should have not only an ethical value but economic relevance too. Finally, this review has summarized the research regarding the use of molecular tools to investigate the HSP response against HS in chickens. The entire process of response to HS is complex and not all molecular steps are completely understood. At this stage, studies working with specific genes give only a partial view of what could happen in the tissues of chickens during HS. In particular, the role played by the different gene families in the HSP-mediated response is not yet clear. Further studies should be carried out to better clarify the mechanisms involved in HS tolerance and to understand if the HSP family (and which HSP family) could be considered a useful biomarker for detecting HS. Supplementary Materials The following are available online at https://www.mdpi.com/2076-2615/11/1/46/s1 , Table S1: All QTLs related to body temperature in chicken. Coord_A_bp = QTL starting point on a chromosome. Coord_B_bp = QTL ending point on a chromosome. All QTLs are significant. Click here for additional data file. Author Contributions Writing original draft, bibliography source, review and editing, F.P.; performed Panther analysis, writing original draft, editing and review, F.C.; editing and review, G.R.; writing original draft, editing and review, C.C.; supervision, review and editing, M.C.; writing original draft, review, editing, supervision and project administration, E.L. All authors have read and agreed to the published version of the manuscript. Funding This research received no external funding. Institutional Review Board Statement Not applicable. Informed Consent Statement Not applicable. Conflicts of Interest The authors declare no conflict of interest. Footnotes Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. References 1. OECD/FAO . OECD-FAO Agricultural Outlook 2017–2026. OECD Publishing; Paris, France: 2017. 2. Givens D.I. Milk and meat in our diet: Good or bad for health? Animals. 2010;4:1941–1952. doi: 10.1017/S1751731110001503. 3. Pereira P.M.D.C.C., Vicente A.F.D.R.B. Meat nutritional composition and nutritive role in the human diet. Meat Sci. 2013;93:586–592. doi: 10.1016/j.meatsci.2012.09.018. 4. FAOSTAT. [(accessed on 27 October 2020)]; Available online: http://www.fao.org/faostat/en/#data/QL . 5. Mottet A., Tempio G. Global poultry production: Current state and future outlook and challenges. World’s Poult. Sci. J. 2017;73:245–256. doi: 10.1017/S0043933917000071. 6. Alexandratos N., Bruinsma J. World Agriculture Towards 2030/2050: The 2012 Revision. FAO, Food and Agriculture Organization of the United Nations; Rome, Italy: 2012. ESA Working Paper No. 12-03. 7. Thornton P.K. Livestock production: Recent trends, future prospects. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2010;365:2853–2867. doi: 10.1098/rstb.2010.0134. 8. Rojas-Downing M.M., Nejadhashemi A.P., Harrigan T., Woznicki S.A. Climate change and livestock: Impacts, adaptation, and mitigation. Clim. Risk Manag. 2017;16:145–163. doi: 10.1016/j.crm.2017.02.001. 9. Thornton P.K., Ericksen P.J., Herrero M., Challinor A.J. Climate variability and vulnerability to climate change: A review. Glob. Chang. Biol. Bioenergy. 2014;20:3313–3328. doi: 10.1111/gcb.12581. 10. Bett B., Kiunga P., Gachohi J., Sindato C., Mbotha D., Robinson T., Lindahl J., Grace D. Effects of climate change on the occurrence and distribution of livestock diseases. Prev. Vet. Med. 2017;137:119–129. doi: 10.1016/j.prevetmed.2016.11.019. 11. Al-Dawood A. Towards heat stress management in small ruminants: A review. Ann. Anim. Sci. 2017;17:59–88. doi: 10.1515/aoas-2016-0068. 12. Sejian V., Bagath M., Krishnan G., Rashamol V.P., Pragna P., Devaraj C., Bhatta R. Genes for resilience to heat stress in small ruminants: A review. Small Rumin. Res. 2019;173:42–53. doi: 10.1016/j.smallrumres.2019.02.009. 13. Leinonen I., Williams A.G., Kyriazakis I. The effects of welfare-enhancing system changes on the environmental impacts of broiler and egg production. Poult. Sci. 2014;93:256–266. doi: 10.3382/ps.2013-03252. 14. Alameen A.O., Abdelatif A.M. Metabolic and endocrine responses of crossbred dairy cows in relation to pregnancy and season under tropical conditions. Am. Eurasian J. Agric. Environ. Sci. 2012;12:1065–1074. 15. McManus C.M., Paludo G.R., Louvandini H., Gugel R., Sasaki L.C.B., Paiva S.R. Heat tolerance in Brazilian sheep: Physiological and blood parameters. Trop. Anim. Health Prod. 2009;41:95–101. doi: 10.1007/s11250-008-9162-1. 16. Das R., Sailo L., Verma N., Bharti P., Saikia J. Impact of heat stress on health and performance of dairy animals: A review. Vet. World. 2016;9:260–268. doi: 10.14202/vetworld.2016.260-268. 17. Mariante A.S., Albuquerque M.S.M., Egito A.A., Mcmanus C. Present status of the conservation of livestock genetic resources in Brazil. Livest. Sci. 2009;120:204–212. doi: 10.1016/j.livsci.2008.07.007. 18. Soleimani A.F., Zulkifli I., Omar A.R., Raha A.R. Physiological responses of 3 chicken breeds to acute heat stress. Poult. Sci. 2011;90:1435–1440. doi: 10.3382/ps.2011-01381. 19. Sejian V., Bhatta R., Gaughan J.B., Dunshea F.R., Lacetera N. Adaptation of animals to heat stress. Animals. 2018;12:s431–s444. doi: 10.1017/S1751731118001945. 20. Hoffmann I. Climate change and the characterization, breeding and conservation of animal genetic resources. Anim. Genet. 2010;41:32–46. doi: 10.1111/j.1365-2052.2010.02043.x. 21. Rovelli G., Ceccobelli S., Perini F., Demir E., Mastrangelo S., Conte G., Abeni F., Marletta D., Ciampolini R., Cassandro M., et al. The genetics of phenotypic plasticity in livestock in the era of climate change: A review. Ital. J. Anim. Sci. 2020;19:997–1014. doi: 10.1080/1828051X.2020.1809540. 22. IPCC (Intergovermental Panel on Climate Change) Climate change 2013: The physical science basis. In: Stocker T.F., Qin D., Plattner G.-K., Tignor M., Allen S.K., Boschung J., Nauels A., Xia Y., Bex V., Midgley P.M., editors. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press; Cambridge, UK: New York, NY, USA: 2013. p. 1535. 23. Hertwich E.G., Wood R. The growing importance of scope 3 greenhouse gas emissions from industry. Environ. Res. Lett. 2018;13:104013. doi: 10.1088/1748-9326/aae19a. 24. Angel S.P., Bagath M., Sejian V., Krishnan G., Bhatta R. Expression patterns of candidate genes reflecting the growth performance of goats subjected to heat stress. Mol. Biol. Rep. 2018;45:2847–2856. doi: 10.1007/s11033-018-4440-0. 25. Lara L.J., Rostagno M.H. Impact of Heat Stress on Poultry Production. Animals. 2013;3:356–369. doi: 10.3390/ani3020356. 26. St-Pierre N.R., Cobanov B., Schnitkey G. Economic losses from heat stress by US livestock industries. J. Dairy Sci. 2003;86:E52–E77. doi: 10.3168/jds.S0022-0302(03)74040-5. 27. Zaboli G., Huang X., Feng X., Ahn D.U. How can heat stress affect chicken meat quality? A review. Poult. Sci. 2019;98:1551–1556. doi: 10.3382/ps/pey399. 28. Soliman A., Safwat A.M. Climate Change Impact on Immune Status and Productivity of Poultry as Well as the Quality of Meat and Egg Products. In: Ewis Omran E.S., Negm A., editors. Climate Change Impacts on Agriculture and Food Security in Egypt. Springer Water; Cham, Switzerland: 2020. pp. 481–498. 29. Wang R.H., Liang R.R., Lin H., Zhu L.X., Zhang Y.M., Mao Y.W., Dong P.C., Niu L.B., Zhang M.H., Luo X., et al. Effect of acute heat stress and slaughter processing on poultry meat quality and postmortem carbohydrate metabolism. Poult. Sci. 2017;96:738–746. doi: 10.3382/ps/pew329. 30. Klebl B.M., Ayoub A.T., Pette D. Protein oxidation, ty- rosine nitration, and inactivation of sarcoplasmic reticulum Ca2+- ATPase in low-frequency stimulated rabbit muscle. FEBS Lett. 1998;422:381–384. doi: 10.1016/S0014-5793(98)00053-2. 31. Wang X.J., Song Z.G., Jiao H.C., Lin H. Skeletal muscle fatty acids shift from oxidation to storage upon dexam- ethasone treatment in chickens. Gen. Comp. Endocrinol. 2012;179:319–330. doi: 10.1016/j.ygcen.2012.09.013. 32. Wang X.J., Wei D.I., Song Z., Jiao H.C., Lin H. Effects of fatty acid treatments on the dexamethasone-induced intramuscular lipid accumulation in chickens. PLoS ONE. 2012;7:e36663. doi: 10.1371/journal.pone.0036663. 33. Farag M.R., Alagawany M. Physiological alterations of poultry to the high environmental temperature. J. Therm. Biol. 2018;76:101–106. doi: 10.1016/j.jtherbio.2018.07.012. 34. Lin H., Jiao H.C., Buyse J., Decuypere E. Strategies for preventing heat stress in poultry. World’s Poult. Sci. J. 2006;26:71–85. doi: 10.1079/WPS200585. 35. Lallo C.H.O., Cohen J., Rankine D., Taylor M., Cambell J., Stephenson T. Characterizing heat stress on livestock using the temperature humidity index (THI)—prospects for a warmer Caribbean. Reg. Environ. Chang. 2018;18:2329–2340. doi: 10.1007/s10113-018-1359-x. 36. Hahn G.L., Gaughan J.B., Mader T.L., Eigenberg R.A. Livestock Energetics and Thermal Environment Management. American Society of Agricultural and Biological Engineers; St. Joseph, MI, USA: 2009. Thermal indices and their applications for livestock environments; pp. 113–130. 37. Dikmen S., Hansen P.J. Is the temperature-humidity index the best indicator of heat stress in lactating dairy cows in a subtropical environment? J. Dairy Sci. 2009;92:109–116. doi: 10.3168/jds.2008-1370. 38. McFarlane J.M., Curtis S.E. Multiple concurrent stressors in chicks: 3. Effects on plasma corticosterone and the heterophil: Lymphocyte ratio. Poult. Sci. 1989;68:522–527. doi: 10.3382/ps.0680522. 39. May J.D., Lott B.D. Feed and water consumption patterns of broilers at high environmental temperatures. Poult. Sci. 1992;71:331–336. doi: 10.3382/ps.0710331. 40. Etches R.J., John T.M., Gibbins A.V. Poultry Production in Hot Climates. 2nd ed. CABI; Oxfordshire, UK: 2008. Behavioural, physiological, neuroendocrine and molecular responses to heat stress; pp. 31–66. 41. Niemann H., Kues W.A. Application of transgenesis in livestock for agriculture and biomedicine. Anim. Reprod. Sci. 2003;79:291–317. doi: 10.1016/S0378-4320(03)00169-6. 42. Goddard M. Adaptation and Fitness in Animal Populations: Evolutionary and Breeding Perspectives on Genetic Resource Management. Springer; New York, NY, USA: 2009. Fitness traits in animal breeding programs; pp. 41–52. 43. Wang S.H., Zhang K., Dai Y.P. Advances in genetic engineering of domestic animals. FASE. 2016;3:1–10. doi: 10.15302/J-FASE-2016085. 44. Xu S.S., Li M.H. Recent advances in understanding genetic variants associated with economically important traits in sheep (Ovis aries) revealed by high-throughput screening technologies. FASE. 2017;4:279–288. doi: 10.15302/J-FASE-2017151. 45. Zhang H., Wang Z.P., Wang S.Z., Li H. Progress of genome wide association study in domestic animals. J. Anim. Sci. Biotechnol. 2012;3:26. doi: 10.1186/2049-1891-3-26. 46. Hoffmann T.J., Kvale M.N., Hesselson S.E., Zhan Y., Aquino C., Cao Y., Cawley S., Chung E.-, Connell S., Eshragh J., et al. Next Generation Genome-Wide Association Tool: Design and Coverage of a High-Throughput European-Optimized SNP Array. Genomics. 2011;98:79–89. doi: 10.1016/j.ygeno.2011.04.005. 47. Bush E.C., Ray F., Alvarez M.J., Realubit R., Li H., Karan C., Califano A., Sims P.A. PLATE-Seq for genome-wide regulatory network analysis of high-throughput screens. Nat. Commun. 2017;8:1–7. doi: 10.1038/s41467-017-00136-z. 48. Gebreselassie G., Berihulay H., Jiang L., Ma Y. Review on genomic regions and candidate genes associated with economically important production and reproduction traits in sheep (Ovies aries) Animals. 2020;10:33. doi: 10.3390/ani10010033. 49. Hirschhorn J.N., Daly M.J. Genome-wide association studies for common diseases and complex traits. Nat. Rev. Genet. 2005;6:95. doi: 10.1038/nrg1521. 50. Mebratie W., Reyer H., Wimmers K., Bovenhuis H., Jensen J. Genome wide association study of body weight and feed efficiency traits in a commercial broiler chicken population, a re-visitation. Sci. Rep. 2019;9:1–10. doi: 10.1038/s41598-018-37216-z. 51. Kudinov A.A., Dementieva N.V., Mitrofanova O.V., Stanishevskaya O.I., Fedorova E.S., Larkina T.A., Mishina A.I., Plemyashov K.V., Griffin D.K., Romanov M.N. Genome-wide association studies targeting the yield of extraembryonic fluid and production traits in Russian White chickens. BMC Genomics. 2019;20:270. doi: 10.1186/s12864-019-5605-5. 52. Finotello F., Di Camillo B. Measuring differential gene expression with RNA-seq: Challenges and strategies for data analysis. Brief. Funct. Genom. 2015;14:130–142. doi: 10.1093/bfgp/elu035. 53. Hrdlickova R., Toloue M., Tian B. RNA-Seq methods for transcriptome analysis. Wiley Interdiscip. Rev. RNA. 2017;8:e1364. doi: 10.1002/wrna.1364. 54. Edwards S.L., Beesley J., French J.D., Dunning M. Beyond GWASs: Illuminating the dark road from association to function. Am. J. Med. Genet. 2013;93:779–797. doi: 10.1016/j.ajhg.2013.10.012. 55. Brøndum R.F., Guldbrandtsen B., Sahana G., Lund M.S., Su G. Strategies for imputation to whole genome sequence using a single or multi-breed reference population in cattle. BMC Genom. 2014;15:728. doi: 10.1186/1471-2164-15-728. 56. Ball R.D. Genome-Wide Association Studies and Genomic Prediction. Humana Press; Totowa, NJ, USA: 2013. Designing a GWAS: Power, sample size, and data structure; pp. 37–98. 57. Young D.L., Fields S. The role of functional data in interpreting the effects of genetic variation. Mol. Biol. Cell. 2015;26:3904–3908. doi: 10.1091/mbc.E15-03-0153. 58. Shaffer J.R., Feingold E., Marazita M.L. Genome-wide association studies: Prospects and challenges for oral health. J. Dent. Res. 2012;91:637–641. doi: 10.1177/0022034512446968. 59. Shin H., Shannon C.P., Fishbane N., Ruan J., Zhou M., Balshaw R., Wilson-McManus J.E., Ng R.T., McManus B.M., Tebbutt S.J. Variation in RNA-Seq transcriptome profiles of peripheral whole blood from healthy individuals with and without globin depletion. PLoS ONE. 2014;9:e91041. doi: 10.1371/journal.pone.0091041. 60. Mercer T.R., Gerhardt D.J., Dinger M.E., Crawford J., Trapnell C., Jeddeloh J.A., Mattick J.S., Rinn J.L. Targeted RNA sequencing reveals the deep complexity of the human transcriptome. Nat. Biotechnol. 2012;30:99–104. doi: 10.1038/nbt.2024. 61. Vikman P., Fadista J., Oskolkov N. RNA sequencing: Current and prospective uses in metabolic research. J. Mol. Endocrinol. 2014;53:R93–R101. doi: 10.1530/JME-14-0170. 62. Rajkumar A.P., Qvist P., Lazarus R., Lescai F., Ju J., Nyegaard M., Mors O., Børglum A.D., Li Q., Christensen J.H. Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq. BMC Genom. 2015;16:548. doi: 10.1186/s12864-015-1767-y. 63. Bogin E.K., Avidar Y. Metabolic and biochemical changes during heat stress in relation to the feathering degree of the domestic hen. Avian Pathol. 1995;24:33–44. doi: 10.1080/03079459508419047. 64. Wolf B.O., Walsberg G.E. The role of the plumage in heat transfer processes of birds. Am. Zool. 2000;40:575–584. doi: 10.1093/icb/40.4.575. 65. Eberhart D.E., Washburn K.W. Assessing the effects of the naked neck gene on chronic heat stress resistance in two genetic populations. Poult. Sci. 1993;72:1391–1399. doi: 10.3382/ps.0721391. 66. Wasti S., Sah N., Mishra B. Impact of Heat Stress on Poultry Health and Performances, and Potential Mitigation Strategies. Animals. 2020;10:1266. doi: 10.3390/ani10081266. 67. Merat P. Potential usefulness of the Na (naked neck) gene in poultry production. World’s Poult. Sci. J. 1986;42:124–142. doi: 10.1079/WPS19860010. 68. Cahaner A., Deeb N., Gutman M. Effects of the plumage-reducing naked neck (Na) gene on the performance of fast-growing broilers at normal and high ambient temperatures. Poult. Sci. 1993;72:767–775. doi: 10.3382/ps.0720767. 69. Azoulay Y., Druyan S., Yadgary L., Hadad Y., Cahaner A. The viability and performance under hot conditions of featherless broilers versus fully feathered broilers. Poult. Sci. 2011;90:19–29. doi: 10.3382/ps.2010-01044. 70. Rajkumar U., Reddy B.L.N., Rajaravindra K.S., Niranjan M., Bhattacharya T.K., Chatterjee R.N., Panda A.K., Reddy M.R., Sharma R.P. Effect of naked neck gene on immune competence, serum biochemical and carcass traits in chickens under a tropical climate. Asian-Australas J. Anim. Sci. 2010;23:867–872. doi: 10.5713/ajas.2010.90548. 71. Yalcin S., Testik A., Ozkan S., Settar P., Celen F., Cahaner A. Performance of naked neck and normal broilers in hot, warm, and temperate climates. Poult. Sci. 1997;76:930–937. doi: 10.1093/ps/76.7.930. 72. Fathi M.M., Galal A., El-Safty S., Mahrous M. Naked neck and frizzle genes for improving chickens raised under high ambient temperature: I. Growth performance and egg production. World’s Poult. Sci. J. 2013;69:813–832. doi: 10.1017/S0043933913000834. 73. Galal A., Radwan L.M., Rezik H.H., Ayoub H. Expression levels of HSP70 and CPT-1 in three local breeds of chickens reared under normal or heat stress conditions after the introduction of the naked neck gene. J. Therm. Biol. 2019;80:113–118. doi: 10.1016/j.jtherbio.2018.12.018. 74. Diaz-Real J., Kim S.Y., Velando A. Plumage colour and the expression of stress-related genes in gull chicks. J. Avian Biol. 2017;48:1216–1225. doi: 10.1111/jav.01460. 75. Jiang R.S., Xia W.T., Chen X.Y., Geng Z.Y., Hu Z.Y. Density of contour feathers and heat tolerance in chickens. J. Appl. Anim. Res. 2010;38:253–255. doi: 10.1080/09712119.2010.10539521. 76. Carabaño M.J., Ramón M., Menéndez-Buxadera A., Molina A., Díaz C. Selecting for heat tolerance. Anim. Front. 2019;9:62–68. doi: 10.1093/af/vfy033. 77. Slawinska A., Zampiga M., Sirri F., Meluzzi A., Bertocchi M., Tavaniello S., Maiorano G. Impact of galactooligosaccharides delivered in ovo on mitigating negative effects of heat stress on performance and welfare of broilers. Poult. Sci. 2020;99:407–415. doi: 10.3382/ps/pez512. 78. Hammel H.T., Pierce J.B. Regulation of internal body temperature. Annu. Rev. Physiol. 1968;30:641–710. doi: 10.1146/annurev.ph.30.030168.003233. 79. Chen X.Y., Wei P.P., Xu S.Y., Geng Z.Y., Jiang R.S. Rectal temperature as an indicator for heat tolerance in chickens. Anim. Sci. J. 2013;84:737–739. doi: 10.1111/asj.12064. 80. Van Goor A., Bolek K.J., Ashwell C.M., Persia M.E., Rothschild M.F., Schmidt C.J., Lamont S.J. Identification of quantitative trait loci for body temperature, body weight, breast yield, and digestibility in an advanced intercross line of chickens under heat stress. Genet. Sel. Evol. 2015;47:96. doi: 10.1186/s12711-015-0176-7. 81. Ait-Boulahsen A., Garlich J.D., Edens F.W. Effect of fasting and acute heat stress on body temperature, blood acid-base and electrolyte status in chickens. Comp. Biochem. Physiol. Part A Physiol. 1989;94:683–687. doi: 10.1016/0300-9629(89)90617-8. 82. Berrong S.L., Washburn K.W. Effects of Genetic Variation on Total Plasma Protein, Body Weight Gains, and Body Temperature Responses to Heat Stress. Poult. Sci. 1998;77:379–385. doi: 10.1093/ps/77.3.379. 83. Edens F.W., Siegel H.S. Modification of corticosterone and glucose responses by sympatholytic agents in young chickens during acute heat exposure. Poult. Sci. 1976;55:1704–1712. doi: 10.3382/ps.0551704. 84. Nathan D.B., Heller E.D., Perek M. The effect of short heat stress upon leucocyte count, plasma corticosterone level, plasma and leucocyte ascorbic acid content. Br. Poult. Sci. 1976;17:481–485. doi: 10.1080/00071667608416303. 85. Lin H., Decuypere E., Buyse J. Oxidative stress induced by corticosterone administration in broiler chickens (Gallus gallus domesticus): 1. Chronic exposure. Comp. Biochem. Physiol. Part B Biochem. 2004;139:737–744. doi: 10.1016/j.cbpc.2004.09.013. 86. Ramiah S.K., Awad E.A., Mookiah S., Idrus Z. Effects of zinc oxide nanoparticles on growth performance and concentrations of malondialdehyde, zinc in tissues, and corticosterone in broiler chickens under heat stress conditions. Poult. Sci. 2019;98:3828–3838. doi: 10.3382/ps/pez093. 87. Ferrante V., Mugnai C., Ferrari L., Marelli S.P., Spagnoli E., Lolli S. Stress and reactivity in three Italian chicken breeds. Ital. J. Anim. Sci. 2016;15:303–309. doi: 10.1080/1828051X.2016.1185978. 88. Marelli S.P., Terova G., Cozzi M.C., Lasagna E., Sarti F.M., Cavalchini L.G. Gene expression of hepatic glucocorticoid receptor NR3C1 and correlation with plasmatic corticosterone in Italian chickens. Anim. Biotechnol. 2010;21:140–148. doi: 10.1080/10495391003608621. 89. Greene E.S., Rajaei-Sharifabadi H., Dridi S. Feather HSP70: A novel non-invasive molecular marker for monitoring stress induced by heat exposure in broilers. Poult. Sci. 2019;98:3400–3404. doi: 10.3382/ps/pez120. 90. Xu D., Tian Y. Selenium and Polysaccharides of Atractylodes macrocephala Koidz Play Different Roles in Improving the Immune Response Induced by Heat Stress in Chickens. Biol. Trace Elem. Res. 2015;168:235–241. doi: 10.1007/s12011-015-0351-2. 91. Niu Z.Y., Liu F.Z., Yan Q.L., Li W.C. Effects of different levels of vitamin E on growth performance and immune responses of broilers under heat stress. Poult. Sci. 2009;88:2101–2107. doi: 10.3382/ps.2009-00220. 92. Liu L.L., He J.H., Xie H.B., Yang Y.S., Li J.C., Zou Y. Resveratrol induces antioxidant and heat shock protein mRNA expression in response to heat stress in black-boned chickens. Poult. Sci. 2014;93:54–62. doi: 10.3382/ps.2013-03423. 93. Jahanian R., Rasouli E. Dietary chromium methionine supplementation could alleviate immunosuppressive effects of heat stress in broiler chicks. J. Anim. Sci. 2015;93:3355–3363. doi: 10.2527/jas.2014-8807. 94. Tang J., Chen Z. The protective effect of γ-aminobutyric acid on the development of immune function in chickens under heat stress. J. Anim. Physiol. Anim. Nutr. 2016;100:768–777. doi: 10.1111/jpn.12385. 95. Ohtsu H., Yamazaki M., Abe H., Murakami H., Toyomizu M. Heat stress modulates cytokine gene expression in the spleen of broiler chickens. Poult. Sci. J. 2015;52:282–287. doi: 10.2141/jpsa.0150062. 96. Varasteh S., Braber S., Akbari P., Garssen J., Fink-Gremmels J. Differences in susceptibility to heat stress along the chicken intestine and the protective effects of galacto-oligosaccharides. PLoS ONE. 2015;10:e0138975. doi: 10.1371/journal.pone.0138975. 97. Xu D., Li W., Huang Y., He J., Tian Y. The effect of selenium and polysaccharide of Atractylodes macrocephala Koidz. (PAMK) on immune response in chicken spleen under heat stress. Biol. Trace Elem. Res. 2014;160:232–237. doi: 10.1007/s12011-014-0056-y. 98. Santoro M.G. Biochemical Pharmacology. Elsevier Inc.; Amsterdam, The Netherlands: 2000. Heat shock factors and the control of the stress response; pp. 55–63. 99. Inouye S., Izu H., Takaki E., Suzuki H., Shirai M., Yokota Y., Ichikawa H., Fujimoto M., Nakai A. Impaired IgG production in mice deficient for heat shock transcription factor 1. J. Biol. Chem. 2004;279:38701–38709. doi: 10.1074/jbc.M405986200. 100. Han A.Y., Zhang M.H., Zuo X.L., Zheng S.S., Zhao C.F., Feng J.H., Cheng C. Effect of acute heat stress on calcium concentration, proliferation, cell cycle, and interleukin-2 production in splenic lymphocytes from broiler chickens. Poult. Sci. 2010;89:2063–2070. doi: 10.3382/ps.2010-00715. 101. Slawinska A., Mendes S., Dunislawska A., Siwek M., Zampiga M., Sirri F., Meluzzi A., Tavaniello S., Maiorano G. Avian model to mitigate gut-derived immune response and oxidative stress during heat. BioSystems. 2019;178:10–15. doi: 10.1016/j.biosystems.2019.01.007. 102. Quinteiro-Filho W.M., Calefi A.S., Cruz D.S.G., Aloia T.P.A., Zager A., Astolfi-Ferreira C.S., Piantino Ferreira J.A., Sharif S., Palermo-Neto J. Heat stress decreases expression of the cytokines, avian β-defensins 4 and 6 and Toll-like receptor 2 in broiler chickens infected with Salmonella Enteritidis. Vet. Immunol. Immunopathol. 2017;186:19–28. doi: 10.1016/j.vetimm.2017.02.006. 103. Jin W., Dong C. IL-17 cytokines in immunity and inflammation. Emerging Microbes and Infections 2. Emerg. Microbes Infect. 2013;2:1–5. doi: 10.1038/emi.2013.58. 104. Van Maele L., Carnoy C., Cayet D., Songhet P., Dumoutier L., Ferrero I., Janot L., Erard F., Bertout J., Leger H., et al. TLR5 Signaling Stimulates the Innate Production of IL-17 and IL-22 by CD3 neg CD127 + Immune Cells in Spleen and Mucosa. J. Immunol. 2010;185:1177–1185. doi: 10.4049/jimmunol.1000115. 105. Sharp F.R., Massa S.M., Swanson R.A. Heat-shock protein protection. Trends Neurosci. 1999;22:97–99. doi: 10.1016/S0166-2236(98)01392-7. 106. Kampinga H.H., Craig E.A. The HSP70 chaperone machinery: J proteins as drivers of functional specificity. Nat. Rev. Mol. 2010;11:579–592. doi: 10.1038/nrm2941. 107. De Maio A., Vazquez D. Extracellular heat shock proteins: A new location, a new function. Shock. 2015;40:239–246. doi: 10.1097/SHK.0b013e3182a185ab. 108. Bakthisaran R., Tangirala R., Rao C.M. Small heat shock proteins: Role in cellular functions and pathology. Biochim. Biophys. Acta Proteins Proteom. 2015;1854:291–319. doi: 10.1016/j.bbapap.2014.12.019. 109. Mezquita B., Mezquita C., Mezquita J. Marked differences between avian and mammalian testicular cells in the heat shock induction and polyadenylation of Hsp70 and ubiquitin transcripts. FEBS Lett. 1998;436:382–386. doi: 10.1016/S0014-5793(98)01172-7. 110. Zulkifli I., Liew P.K., Israf D.A., Omar A.R., Hair-Bejo M. Effects of early age feed restriction and heat conditioning on heterophil/lymphocyte ratios, heat shock protein 70 expression and body temperature of heat-stressed broiler chickens. J. Therm. Biol. 2003;28:217–222. doi: 10.1016/S0306-4565(02)00058-X. 111. Mahmoud K.Z., Edens F.W., Eisen E.J., Havenstein G.B. The effect of dietary phosphorus on heat shock protein mRNAs during acute heat stress in male broiler chickens (Gallus gallus) Comp. Biochem. Physiol. C Toxicol. Pharmacol. 2004;137:11–18. doi: 10.1016/j.cca.2003.10.013. 112. Yu J., Bao E., Yan J., Lei L. Expression and localization of Hsps in the heart and blood vessel of heat-stressed broilers. Cell Stress Chaperones. 2008;13:327–335. doi: 10.1007/s12192-008-0031-7. 113. Xie J., Tang L., Lu L., Zhang L., Xi L., Liu H.C., Odle J., Luo X. Differential expression of heat shock transcription factors and heat shock proteins after acute and chronic heat stress in laying chickens (Gallus gallus) PLoS ONE. 2014;9:e102204. doi: 10.1371/journal.pone.0102204. 114. Roushdy E.M., Zaglool A.W., El-Tarabany M.S. Effects of chronic thermal stress on growth performance, carcass traits, antioxidant indices and the expression of HSP70, growth hormone and superoxide dismutase genes in two broiler strains. J. Therm. Biol. 2018;74:337–343. doi: 10.1016/j.jtherbio.2018.04.009. 115. Tang S., Zhou S., Yin B., Xu J., Di L., Zhang J., Bao E. Heat stress-induced renal damage in poultry and the protective effects of HSP60 and HSP47. Cell Stress Chaperones. 2018;23:1033–1040. doi: 10.1007/s12192-018-0912-3. 116. Zhu L., Liao R., Wu N., Zhu G., Yang C. Heat stress mediates changes in fecal microbiome and functional pathways of laying hens. Appl. Microbiol. Biotechnol. 2019;103:461–472. doi: 10.1007/s00253-018-9465-8. 117. Cedraz H., Gromboni J.G.G., Garcia A.A.P., Farias Filho R.V., Souza T.M., De Oliveira E.R., De Oliveira E.B., Do Nascimento C.S., Meneghetti C., Wenceslau A.A. Heat stress induces expression of HSP genes in genetically divergent chickens. PLoS ONE. 2017;12:e0186083. doi: 10.1371/journal.pone.0186083. 118. Wang S.H., Cheng C.Y., Tang P.C., Chen C.F., Chen H.H., Lee Y.P., Huang S.Y. Differential gene expressions in testes of L2 strain Taiwan country chicken in response to acute heat stress. Theriogenology. 2013;79:374–382. doi: 10.1016/j.theriogenology.2012.10.010. 119. Wang S.H., Cheng C.Y., Tang P.C., Chen C.F., Chen H.H., Lee Y.P., Huang S.Y. Acute heat stress induces differential gene expressions in the testes of a broiler-type strain of Taiwan country chickens. PLoS ONE. 2015;10:e0125816. doi: 10.1371/journal.pone.0125816. 120. Sharma A., Kumar B.V.S., Dash S., Singh S., Verma R. Heat shock protein B1 expression is associated with age at sexual maturity in Rhode Island Red and Punjab Red layers under heat stress. Int. J. Biometeorol. 2020;64:1133–1143. doi: 10.1007/s00484-020-01887-z. 121. Radwan L.M. Genetic improvement of egg laying traits in Fayoumi chickens bred under conditions of heat stress through selection and gene expression studies. J. Therm. Biol. 2020;89:102546. doi: 10.1016/j.jtherbio.2020.102546. 122. Zhang Q., Liu L., Zhu F., Ning Z., Hincke M.T., Yang N., Hou Z. Integrating de novo transcriptome assembly and cloning to obtain chicken ovocleidin-17 full-length cDNA. PLoS ONE. 2014;9:e93452. doi: 10.1371/journal.pone.0093452. 123. Zhang W.W., Kong L.N., Zhang D.X., Ji C.L., Zhang X.Q., Luo Q.B. Effect of the C.–1 388 A> G polymorphism in chicken heat shock transcription factor 3 gene on heat tolerance. J. Integr. Agric. 2015;14:1808–1815. doi: 10.1016/S2095-3119(14)60943-6. 124. Chen Z.Y., Gan J.K., Xiao X., Jiang L.Y., Zhang X.Q., Luo Q.B. The association of SNPs in Hsp90β gene 5′ flanking region with thermo tolerance traits and tissue mRNA expression in two chicken breeds. Mol. Biol. Rep. 2013;40:5295–5306. doi: 10.1007/s11033-013-2630-3. 125. Kong L.N., Zhang D.X., Ji C.L., Zhang X.Q., Luo Q.B. Association analysis between SNPs in the 5’-flanking region of the chicken GRP78 gene, thermotolerance parameters, and tissue mRNA expression. Genet. Mol. Res. 2015;14:6110–6123. doi: 10.4238/2015.June.8.9. 126. Endo S., Hiramatsu N., Hayakawa K., Okamura M., Kasay A., Tagawa Y., Sawada N., Yao J., Kitamura M. Geranylgeranylacetone, an inducer of the 70-kDa heat shock protein (HSP70), elicits unfolded protein response and coordinates cellular fate independently of HSP70. Mol. Pharmacol. 2007;72:1337–1348. doi: 10.1124/mol.107.039164. 127. Irivboje Y.I., Sanni M.T., Fafiolu A.O., Olowofeso O., Ikeobi C.O.N. Genetic polymorphisms in part of intron 7 and exon 8 of HSP90AA1 gene and its association with heat tolerance traits in two exotic layer chicken strains. Trop. Anim. Health Prod. 2020;52:969–977. doi: 10.1007/s11250-019-02079-7. 128. Fleming D.S., Weigend S., Simianer H., Weigend A., Rothschild M., Schmidt C., Ashwell C., Persia M., Reecy J., Lamont S.J. Genomic comparison of indigenous African and northern European chickens reveals putative mechanisms of stress tolerance related to environmental selection pressure. G3-Genes Genom. Genet. 2017;7:1525–1537. doi: 10.1534/g3.117.041228. 129. Srikanth K., Kumar H., Park W., Byun M., Lim D., Kemp S., te Pas M.F.W., Kim J.M., Park J.E. Cardiac and Skeletal Muscle Transcriptome Response to Heat Stress in Kenyan Chicken Ecotypes Adapted to Low and High Altitudes Reveal Differences in Thermal Tolerance and Stress Response. Front. Genet. 2019;10:993. doi: 10.3389/fgene.2019.00993. 130. Sun H., Jiang R., Xu S., Zhang Z., Xu G., Zheng J., Qu L. Transcriptome responses to heat stress in hypothalamus of a meat-type chicken. J. Anim. Sci. Biotechnol. 2015;6:6. doi: 10.1186/s40104-015-0003-6. 131. Kim J.S., Lee Y.H., Chang Y.U., Yi H.K. PPARγ regulates inflammatory reaction by inhibiting the MAPK/NF-κB pathway in C2C12 skeletal muscle cells. J. Physiol. Biochem. 2017;73:49–57. doi: 10.1007/s13105-016-0523-3. 132. Kumar H., Iskender A.U., Srikanth K., Kim H., Zhunushov A.T., Chooq H., Jang G.W., Lim Y., Song K.D., Park J.E. Transcriptome of chicken liver tissues reveals the candidate genes and pathways responsible for adaptation into two different climatic conditions. Animals. 2019;9:1076. doi: 10.3390/ani9121076. 133. Wells K.L., Hadad Y., Ben-Avraham D., Hillel J., Cahaner A., Headon D.J. Genome-wide SNP scan of pooled DNA reveals nonsense mutation in FGF20 in the scaleless line of featherless chickens. BMC Genom. 2012;13:257. doi: 10.1186/1471-2164-13-257. 134. Widelitz R.B., Jiang T.X., Lu J., Chuong C.M. Beta-catenin in epithelial morphogenesis: Conversion of part of avian foot scales into feather buds with a mutated beta-catenin. Dev. Biol. 2000;219:98–114. doi: 10.1006/dbio.1999.9580. 135. Deeb N., Cahaner A. The effects of naked neck genotypes, ambient temperature, and feeding status and their interactions on body temperature and performance of broilers. Poult. Sci. 1999;78:1341–1346. doi: 10.1093/ps/78.10.1341. 136. Te Pas M.F.W., Park W., Srikanth K., Kemp S., Kim J.M., Lim D., Park J.E. Transcriptomic profiles of muscle, heart, and spleen in reaction to circadian heat stress in Ethiopian highland and lowland male chicken. Cell Stress Chaperones. 2019;24:175–194. doi: 10.1007/s12192-018-0954-6. 137. Berger S. The complex language of chromatin regulation during transcription. Nature. 2007;447:407–412. doi: 10.1038/nature05915. 138. Duncan E.J., Gluckman P.D., Dearden P.K. Epigenetics, plasticity, and evolution: How do we link epigenetic change to phenotype? J. Exp. Zool. B Mol. Dev. Evol. 2014;322:208–220. doi: 10.1002/jez.b.22571. 139. Berger S.L., Kouzarides T., Shiekhattar R., Shilatifard A. An operational definition of epigenetics. Genes Dev. 2009;23:781–783. doi: 10.1101/gad.1787609. 140. Tzschentke B., Basta D. Early development of neuronal hypothalamic thermosensitivity in birds: Influence of epigenetic temperature adaptation. Comp. Biochem. Physiol. Part A: Mol. Integr. Physiol. 2002;131:825–832. doi: 10.1016/S1095-6433(02)00020-X. 141. Jablonka E., Raz G. Transgenerational Epigenetic Inheritance: Prevalence, Mechanisms, and Implications for the Study of Heredity and Evolution. Q. Rev. Biol. 2009;84:131–176. doi: 10.1086/598822. 142. Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16:6–21. doi: 10.1101/gad.947102. 143. Yang C., Zhang M., Niu W., Yang R., Zhang Y., Qiu Z., Sun B., Zhao Z. Analysis of DNA methylation in various swine tissues. PLoS ONE. 2011;6:e16229. doi: 10.1371/journal.pone.0016229. 144. Vinoth A., Thirunalasundari T., Shanmugam M., Uthrakumar A., Suji S., Rajkumar U. Evaluation of DNA methylation and mRNA expression of heat shock proteins in thermal manipulated chicken. Cell Stress Chaperones. 2018;23:235–252. doi: 10.1007/s12192-017-0837-2. 145. Gan J.K., Zhang D.X., He D.L., Zhang X.Q., Chen Z.Y., Luo Q.B. Promoter methylation negatively correlated with mRNA expression but not tissue differential expression after heat stress. Genet. Mol. Res. 2013;12:809–819. doi: 10.4238/2013.March.15.1. 146. Rajkumar U., Vinoth A., Shanmugam M., Rajaravindra K.S., Rama Rao S.V. Effect of Embryonic Thermal Exposure on Heat Shock Proteins (HSPs) Gene Expression and Serum T3 Concentration in Two Broiler Populations. Anim. Biotechnol. 2015;26:260–267. doi: 10.1080/10495398.2015.1022183. 147. Jastrebski S.F., Lamont S.J., Schmidt C.J. Chicken hepatic response to chronic heat stress using integrated transcriptome and metabolome analysis. PLoS ONE. 2017;12 doi: 10.1371/journal.pone.0181900. 148. Kumar H., Choo H., Iskender A.U., Srikanth K., Kim H., Zhunushov A.T., Jang G.W., Lim Y., Song K.D., Park J.E. RNA seq analyses of chicken reveals biological pathways involved in acclimation into different geographical locations. Sci. Rep. 2020;10:19288. doi: 10.1038/s41598-020-76234-8. 149. Li C., Wang X., Wang G., Li N., Wu C. Expression analysis of global gene response to chronic heat exposure in broiler chickens (Gallus gallus) reveals new reactive genes. Poult. Sci. 2011;90:1028–1036. doi: 10.3382/ps.2010-01144. 150. Bertocchi M., Zampiga M., Luise D., Vitali M., Sirri F., Slawinska A., Tavaniello S., Palumbo O., Archetti I., Maiorano G., et al. In ovo injection of a galacto-oligosaccharide prebiotic in broiler chickens submitted to heat-stress: Impact on transcriptomic profile and plasma immune parameters. Animals. 2019;9:1067. doi: 10.3390/ani9121067. 151. Coble D.J., Fleming D., Persia M.E., Ashwell C.M., Rothschild M.F., Schmidt C.J., Lamont S.J. RNA-seq analysis of broiler liver transcriptome reveals novel responses to high ambient temperature. BMC Genom. 2014;15:1084. doi: 10.1186/1471-2164-15-1084. 152. Walugembe M., Bertolini F., Dematawewa C.M.B., Reis M.P., Elbeltagy A.R., Schmidt C.J., Lamont S.J., Rothschild M.F. Detection of selection signatures among Brazilian, Sri Lankan, and Egyptian chicken populations under different environmental conditions. Front. Genet. 2019;10:737. doi: 10.3389/fgene.2018.00737. 153. Sun L., Lamont S.J., Cooksey A.M., McCarthy F., Tudor C.O., Vijay-Shanker K., DeRita R.M., Rothschild M., Ashwell C., Persia M.E., et al. Transcriptome response to heat stress in a chicken hepatocellular carcinoma cell line. Cell Stress Chaperones. 2015;20:939–950. doi: 10.1007/s12192-015-0621-0. Associated Data Supplementary Materials Click here for additional data file.