Introduction

Climate change is a complex phenomenon with disastrous consequences in almost all aspects of human societies, including in the livestock sector. Human activities would have contributed to a 1.0 °C increase in global temperature above pre-industrial levels. If current trends continue, global warming will likely approach 1.5 °C between 2030 and 20521. Furthermore, according to reports from several international organizations involved in climate change research, practically all sectors of the economy and all regimes around the world will suffer from the negative effects of climate change, although to varying degrees, such as the irreversible loss of some natural ecosystems2,3,4,5. Particular attention should be paid to temperature variations as the most significant stressor on the efficiency of livestock production systems, affecting animal development, growth, production outputs and reproduction efficiency.

Compared to other ruminants, sheep tend to perform better in harsh environments, nevertheless, heat stress was shown to negatively impact physiology6,7,8, production parameters (milk production quality and quantity9,10, growth rates11), reproduction (male12 and female fertility13), sheep health14, and sheep welfare15, lowering farms economic returns and efficiency. Van Wettere et al.16 reported that when temperatures drop below 12 °C or increase over 31 °C, thermoregulatory mechanisms are severely disrupted and sheep's ability to maintain homeothermy is disturbed, with the impacts on performance and welfare being severe.

Heat stress is cytotoxic as it changes biological molecules, disrupts cell activities, affects metabolic responses, causes oxidative cell damage, and activates apoptosis and necrosis pathways17,18,19. In ewes, heat stress reduces estradiol concentration and aromatase activity, negatively influencing estrous occurrence and duration20,21,22. For each additional day of high-temperature exposure (> 32 °C) in a week prior to estrous onset, 2.7% of the fertilization rate and 3.5% of the lambing rate reduction were observed16. Heat stress also leads to higher circulating progesterone levels during the luteal phase, implying increased progesterone synthesis and decreased clearance, which impacts the timing and development of pre-ovulatory follicle growth13. Research with Malpura sheep housed in a climatic chamber for 6 h at 42 °C and 54% relative humidity found a significant drop in hemoglobin and packed cell volume levels as more water is transported through the thermoregulation circulatory system. Furthermore, higher plasma cortisol levels and increased cholesterol catabolism were involved in improving gluconeogenesis and providing more energy to heat-stressed rams. Because of the lower Gonadotropin releasing hormone (GnRH) release from the hypothalamus, plasma triiodothyronine (T3), thyroxine (T4) and testosterone levels were lowered, as were the sexual behaviour and semen volume, with less progressive sperm motility and sperm concentration23,24.

A critical component of adaptation to climate change is represented by the animals' genetic capacity to survive under stressful climatic conditions. Since it results in a permanent and cumulative transformation, the selection of thermoresistant animals could represent an efficient method of increasing livestock productivity during periods of high environmental temperature. Identifying and utilizing thermo-tolerant genotypes in sheep is critical due to the changing climate scenario, having the potential to significantly influence sheep productivity.

Thermotolerant genes in sheep have been studied by employing various genomic tools, for example, in heat stress protein genes; HSP90 and HSP70 genes polymorphism were investigated in Indian sheep using PCR–RFLP method25, FGF2, GNAI13 and PLCB1 melanogenesis candidate genes were studied in Egyptian Barki sheep breed, generated by 50K SNPs Beadchips26, and HSP5 and HSP40 in Brazilian sheep using the 50K SNP Chip27. Given that genetics has a complex role in affecting an individual’s capacity to respond to a stressful situation, more research into mechanisms and the development of more appropriate tools is required.

This study aimed to investigate the polymorphism of 51 SNPs in 29 genes involved in thermotolerance throughout the use of Kompetitive Allele Specific PCR (KASP-PCR) technique in 17 sheep breeds originating from different climatic conditions. The main goal was to identify and describe polymorphisms related to climatic adaptation that could be used to develop future thermal resilience in sheep through genomic selection.

Results

Allele and genotype frequency

A total of 601 animals (83.47%) from the total 720 animals were successfully genotyped and 32 SNPs (62.74%) among the initial set of 51 SNPs were successfully genotyped in this study; 17 of them were found to be polymorphic (33.33%) (see Table S1 of the Supplementary data), which were HSPA12A, HSP90AA1, IL33, DIO2, BTNL2, CSN2, ABCG1, CSN1S1, GHR, HSPA8, STAT3, and HCRT. Allelic and genotypic frequencies presented in Tables S2 and S3 of the Supplementary data, were different from one population to another. Four SNPs were successfully genotyped in all breeds: rs161504783-HSP12A, rs588145625-HSPA8, rs588498137-STAT3 and rs602521720-HCRT.

The heterozygote TC for rs161504783-HSPA12A was dominant, except for Hungarian Racka, Transylvanian Merino, Hungarian Merino, Botosani Karakul and Sardi. The T and C allele frequencies were almost equally frequent in most breeds, except in Sardi (T allele = 0.107 and C alleles = 0.828) and Botosani Karakul (T allele = 0.735 and C alleles = 0.265).

The heterozygote GA for SNP rs588145625-HSPA8 was absent in the heat tolerance breeds, except the Transylvanian Merino, and in breeds with high cold tolerance (Hungarian Racka, Babolna Tetra, Hungarian Tsigai, Romanian Racka, Pramenka and Turcana). The homozygote GG was present in all heat-tolerant breeds and some cold-tolerant breeds (Suffolk, Ile de France and Hungarian Merino). The G allele was dominant in all breeds, ranging from 0.760 to 1, and the A allele has only appeared in cold-tolerant breeds and Transylvanian Merino with frequency ranging from 0 to 0.308.

No patterns have been found for allelic and genotypic frequencies for rs588498137-STAT3 and rs602521720-HCRT SNPs.

The GG genotype and G allele were dominant in all breeds for SNP rs588498137-STAT3, with G allelic frequency varying from 0.750 to 1. Similarly, the CC genotype and C allele for SNP rs602521720-HCRT was dominant in all breeds, with C allelic frequency varying from 0.546 to 1.

The Hardy–Weinberg equilibrium (HWE) was investigated with Chi-square (x2) test using allelic frequencies, observed and expected genotypes, and P value of the polymorphic genes are summarized in Table S2. Botosani Karakul breed from Romania was found to be the breed with most SNPs deviated from HWE (P ≤ 0.05); rs416259751-IL33, rs411181557-DIO2, rs414917134-BTNL2 and rs420611298-ABCG1.

Genetic diversity and interrelationship between SNPs

SNPs data were used to perform Principal Component Analysis (PCA) to highlight differences between either the breeds (Fig. 1a) or the climatic characteristics (Fig. 1b).

Figure 1
figure 1

Score biplot of Principal Component Analysis of the 17 SNPs and 601 animals. (a) Individuals of different breeds are differently colored; (AWAS Hungarian Awassi, BG Béni Guil, DM D’Man, HUMER Hungarian Merino, HURAC Hungarian Racka, HUTSI Hungarian Tsigai, IDF Ile de France, KAR Botosani Karakul, PRA Pramenka, RORAC Romanian Racka, ROTSI Romanian Tsigai, SAR Sardi, SUF Suffolk, TET Babolna Tetra, TIM Timahdite, TRANSMER Transylvanian Merino and TUR Turcana). (b) Breeds grouped by climatic characteristics; COLD: cold-tolerant breeds (Suffolk, Babolna Tetra, Ile de France, Hungarian Tsigai, Hungarian Racka, Hungarian Merino, Pramenka, Romanian Racka and Turcana), HOT: heat tolerant breeds originated from Morocco (Béni Guil, D’Man, Timahdite and Sardi), HOT EU: Heat tolerant breeds reared in Europe (Hungarian Awassi, Botosani Karakul, Transylvanian Merino and Romanian Tsigai).

As displayed in Fig. 1, the PC1 and PC2 account for 11.21% and 9.98% of the total variation in the 17 breeds, respectively. Both PC1 and PC2 were unable to clearly separate neither the breeds (Fig. 1a), nor the climatic regions (Fig. 1b) even after we separated the heat tolerant breeds reared in the EU and in Morocco. All breeds were mostly concentrated in − 2.50 < PC1 < 1.12; cold tolerant breeds were outspread in PC1 score of − 5.593 to − 5.013 and PC2 score of − 4.149 to 4.020, heat-tolerant breeds were outspread in PC1 score of − 0.452 to 0.460 and PC2 score of − 0.157 to 1.183, while heat tolerant breeds reared in Europe had PC1 and PC2 ranging from − 0.881 to 0.686 and − 3.987 to 2.931, respectively. The highest contributions to the principal component were by rs397514117-HSP90AA1 (c) and rs397514272-HSP90AA1 (e), whereas rs410259751-IL33 (g) also had a positive contribution (Fig. 2). The loading value and score value of PCA are available on Tables S4 and S5 of the Supplementary data.

Figure 2
figure 2

Loading biplot of principal component analysis of the 17 SNPs for 17 sheep breeds.

Discussion

The genetic background of thermotolerance adaptation in sheep represents a great interest for current research, especially when faced with the imminent increase of the earth's temperature because of climate change. Different breeds from different regions were used in this study, with the aim to see the different genetic backgrounds of thermotolerance adaptation according to geographical origins.

The application of KASP-PCR assay in this study did not show a high assay success rate as in previous studies carried out on goats28 and dairy cattle29. As much as 62.74% of SNPs were successfully genotyped in this study. These relatively high (37.26%) unsuccessful genotypes could be the consequence of sample damage during transportation from the laboratory to the subcontracting laboratory. Furthermore, quality control results in lower data for the analysis.

Only four SNPs were successfully genotyped in all investigated populations, and the focus was given to SNPs HSPA12A and HSPA8 as they are the only ones that showed allelic and genotypic frequency patterns with the climatic characteristics. HSPs are a wide family of chaperone proteins categorized according to their molecular size and amino acid sequence similarity. The HSPA12A and HSPA8 are members of the HSP70 family with a molecular weight of 70 kDa30. It is the biggest, most numerous and most conserved protein family throughout evolution, as well as the most extensively studied protein family across a wide range of species31,32,33. It is found on chromosome 15 and is made up of nine exons in sheep. This gene is widely studied for thermal adaptability. It controls cellular survival to heat stress with a highly dynamic nature, and these proteins are responsible for maintaining the organism's equilibrium and acclimating to heat stress34,35.

HSPA12A affects aspects like cellular senescence and how cells respond to heat stress36. It was found to be more active in ruminants during the summer, which helps them adapt to a hostile environment. This gene was found to be highly associated with heat tolerance in indicine cattle (Bos taurus indicus) and also has an important role in the water holding capacity in beef breeds, which might also act as a mechanism in surviving heat stress through the increased adrenaline and pH change in muscles37,38. According to the genotype and allele frequency results, the heterozygote TC for rs161504783-HSPA12A was dominant in most of the breeds, except Hungarian Racka, Hungarian Merino, Transylvanian Merino, Botosani Karakul and Sardi. Béni Guil, D’Man and Timahdite, breeds from a hot region, were not in HWE for HSPA12A.

In HSP70, exon 1 is non-coding, while the following eight exons combine to form the HSPA8 protein, with 650 amino acids of 71 kDa. HSPA8 aids in the day-to-day cell functions of protein folding and unfolding, polypeptide aggregation prevention, disassembly of large protein complexes and protein translocation across cellular compartments35,39. Given that increasing HSPA8 levels have been discovered to be positively correlated with heat tolerance, this gene has been employed as a candidate for heat resistance in many livestock species40,41,42.

In this study, the homozygote GG was carried by all heat-tolerant breeds and some cold-tolerant breeds (Suffolk, Ile de France and Romanian Tsigai), while the heterozygote GA for SNP rs588145625-HSPA8 was only found in cold-tolerant breeds (Hungarian Racka, Babolna Tetra, Hungarian Tsigai, Romanian Racka, Transylvanian Merino, Pramenka and Turcana). The A allele was absent in heat tolerant breeds, which is opposed to the observation made by25 in a study of gene expression in Indian sheep, which discovered that animals with A allele have better hot climate adaptability, compared to animals with the G allele. The AA genotype is more adaptable to a hot environment and has lower HSPA8 gene expression, than animals with the AG genotype. Similarly, a study by41 showed that the GG genotype has the least ability to survive heat stress in Awassi and Arabi sheep.

In this study, Botosani Karakul has shown to be the breed with the largest number of SNPs deviating from HWE. Lower heterozygosity was observed in 10 SNPs, and 4 of them deviated from HWE. One possible reason for this deviation of HWE was because Botosani Karakul is one of the Romanian crossbreeds that has undergone extensive genetic mixing from its German and Austria lines, and other Romanian breeds since its introduction from Russia at the beginning of the nineteenth century43,44.

PCA was unable to differentiate each breed in a clear and distinct manner; however, a relative clustering was observed based on the three different categories of climatic regions (Fig. 1b). This could be attributed to the fact that heat-tolerant breeds that are kept in the EU have become acclimated to the subtropical environment, in addition to the fact that genetic ad-mixture has occurred as a result of the widespread use of reproductive technologies, which has led to a less distinct genetic divergence between breeds, with previous studies reporting high levels of ad-mixture between different sheep breeds reared in Hungary and Romania, with the aim to increase the productivity45,46,47,48,49. We acknowledge that the cold breeds used in this study were not from year-around cold regions (e.g. Iceland, Finland, Norway), which was one of our limitations due to our inability to obtain samples from these regions. However, we strongly believe that the cold tolerance breeds from temperate climates utilized in this study are sufficiently contrasting with the Moroccan heat tolerant breeds.

From the loading biplot, HSP90AA1 and IL33 SNPs significantly contributed to PC1 and PC2. HSP90AA1 has been confirmed to be associated with thermal stress susceptibility of sheep in previous studies50,51, while IL33 was found to affect the sheep immunity and resistance to gastrointestinal intestinal nematode infection52,53, which is also associated because heat stress promotes to immune suppression and increases animal vulnerability to illnesses54. Our PCA findings validate these two SNPs as potential candidates for heat adaptability across the sheep breeds investigated in the current study, although more research is needed in order to clarify the relationship between the investigated SNPs and heat resistance.

Conclusion

Based on our 17 SNP polymorphism analyses performed on 601 animals, we found that the KASP-PCR method represents a feasible method for investigating polymorphisms in different sheep breeds. Furthermore, based on allele and genotype frequency, we validated that HSPA12A and HSPA8 SNPs are potential candidate markers for thermotolerance adaptation in sheep, whereas principal component analysis confirmed that HSP90AA1 and IL33 SNPs were the primary potential candidates. The results contribute to an increase in knowledge regarding the genetic variability of SNPs for thermotolerance adaptation in sheep. However, more studies are needed in order to clarify the relationship between the studied SNPs and heat resistance in sheep.

Method

Genomic DNA extraction

Samples were collected from 720 sheep belonging to17 breeds adapted to different climatic conditions, originating from 4 countries (Table 1). The majority of the breeds studied are indigenous to the country of origin, with some exotic breeds being included to determine if their acclimatization has contributed to adaptation. The breed characteristic (hot or cold tolerant) was determined based on the origin and development history of each breed in that particular country, as well as the environmental conditions under which the sample was collected (Table 2). All research activities were conducted in compliance with the European Union's Animal Experimentation Directive (Directive 2010/63/EU). FAO/IAEA55 recommended method from hair follicles and56 method from the blood was used for DNA genomic isolation. The DNA was kept at − 20 °C until analysis. The concentration of DNA was determined using a NanoDrop Spectrophotometer (Thermo Scientific, Waltham, MA, USA). All samples were diluted to a uniform concentration, and genotyping was performed using the equivalent of 50 ng of DNA per sample.

Table 1 Samples origin and breed characteristics.
Table 2 Sampling region climatological details.

Selection of SNPs

Based on the findings of previous genome-wide association studies (GWAS) and marker-assisted selection studies across the sheep genome25,41,57,58,59, a SNP panel of 51 SNPs from 29 genes related to heat stress was selected, which included loci distributed on 18 chromosomes (Table 3). The Ovis SNP data were obtained from the Single Nucleotide Polymorphism Database (dbSNP) from the National Center for Biotechnology Information (NCBI) or Ensembl a few years ago, some of the RS reference IDs might have changed up to the present time. Due to the frequently changing genotyping database and varying reference sequences from those available at the European Variation Archive (EVA), only 10 SNPs out of 17 SNPs are deposited on EVA.

Table 3 Selected SNPs used for studying thermotolerance adaptation genes in various sheep breeds.

Genotyping and quality control

The bi-allelic discrimination of the selected 51 SNPs was performed using Kompetitive Allele Specific PCR (KASPTM, LGC Genomics, Teddington, Middlesex, UK). SNP Viewer software version 1.99 (Hoddesdon, UK) was used to visualize the results. All genotype data were exported for statistical analysis. Only SNPs that appeared in at least 50% of the breeds were considered. Data quality control of genotyped data included discarding animals with a call rate of less than 50% and the SNPs with call rates < 50%. This led to discrepancies in either the number of animals per breed or the number of SNPs per animal.

Data analysis

The raw allele calls obtained from LGC Genomics were analyzed using LGC Genomics' KlusterCaller program. Gene diversity, allele and genotype frequencies, and their accordance with or deviation from the Hardy–Weinberg equilibrium were determined by POPGENE software version 1.3160

The Principal Component Analysis was done using FactoMineR61 and ggplot262 packages from the R Program63 to visualize the genetic divergences between sheep breeds that were divided according to their climatic characteristics; cold-tolerant breeds (Babolna Tetra, Hungarian Merino, Hungarian Racka, Hungarian Tsigai, Ile de France, Pramenka, Romanian Racka, Suffolk and Turcana), heat tolerance breeds originated from Morocco (Béni Guil, D’Man, Timahdite and Sardi), and heat tolerant breeds reared in Europe (Hungarian Awassi, Botosani Karakul, Transylvanian Merino, and Romanian Tsigai).

Ethical approval

The authors confirm that the experiment complied with the European Union's Directive on Animal Experimentation (Directive 2010/63/EU) and ARRIVE guidelines. All animals in the experiment underwent standard procedures without experiencing any harm or discomfort, and all procedures were carried out in compliance with applicable guidelines and regulations. The study was approved by the Scientific and Ethics Committee of Centre for Agricultural Genomics and Biotechnology, University of Debrecen (Ethics statement No. 07).