SMARTER-database: a tool to integrate SNP array datasets for sheep and goat breeds.

Underutilized sheep and goat breeds can adapt to challenging environments due to their genetics. Integrating publicly available genomic datasets with new data will facilitate genetic diversity analyses; however, this process is complicated by data discrepancies, such as outdated assembly versions or different data formats. Here, we present the SMARTER-database, a collection of tools and scripts to standardize genomic data and metadata, mainly from SNP chip arrays on global small ruminant populations, with a focus on reproducibility. SMARTER-database harmonizes genotypes for about 12,000 sheep and 6,000 goats to a uniform coding and assembly version. Users can access the genotype data via File Transfer Protocol and interact with the metadata through a web interface or using their custom scripts, enabling efficient filtering and selection of samples. These tools will empower researchers to focus on the crucial aspects of adaptation and contribute to livestock sustainability, leveraging the rich dataset provided by the SMARTER-database. Availability and implementation: The code is available as open-source software under the MIT license at https://github.com/cnr-ibba/SMARTER-database.

African Local Pig Genetic Resources in the Context of Climate Change Adaptation.

Africa is home to a wide diversity of locally adapted pig breeds whose genetic architecture offers important insights into livestock adaptation to climate change. However, the majority of these inherent traits have not been fully highlighted. This review presents an overview of the current state of African pig genetic resources, providing highlights on their population and production statistics, production system, population diversity indices, and genomic evidence underlying their evolutionary potential. The study results reveal an incomplete characterization of local pig genotypes across the continent. The characterized population, however, demonstrates moderate to high levels of genetic diversity, enough to support breeding and conservation programs. Owing to low genetic differentiation and limited evidence of distinct population structures, it appears that most local pig populations are strains within larger breeds. Genomic evidence has shown a higher number of selection signatures associated with various economically important traits, thus making them potential candidates for climate change adaptation. The reportedly early evidence of hybridization with wild suid groups further suggests untapped insights into disease resistance and resilience traits that need to be illuminated using higher-density markers. Nevertheless, gene introgression from commercial breeds is prevalent across Africa; thus, efforts to realize and utilize these traits must increase before they are permanently depleted.

Assessing the population structure and genetic variability of Kenyan native goats under extensive production system.

Indigenous goats are important to many livelihoods. Despite this, they are subjected to indiscriminate crossbreeding. This affects their genetic variability which is needed to survive in current regime of climate change. The study assessed population structure and genetic diversity of Galla and Small East African goats (SEA) using pedigree information. A total of 7384 animals, 5222 of the Galla and 2162 of the SEA breeds, born between the years 1983 and 2022, were utilized. Individuals with known parents were defined as reference population. From the results, the maximum generation traced for Galla and SEA populations was 14.6 and 14.5, respectively. However, only 6 and 5 generations for Galla and SEA were complete. Pedigree completeness increased with the increasing number of generations. The average generation interval (GI) for Galla and SEA was 3.84 ± 0.04 and 4.4 ± 0.13 years. The average increase in the rate of inbreeding per generation for Galla and SEA was 0.04 and 0.05, with the effective ancestors and founders (fa/fe) ratio being same (1.00) for both breeds. Fifty percent (50%) of the genetic variability in the populations was contributed by 3 and 1 ancestor for Galla SEA, respectively. The effective population size (Ne) was 5.19 and 4.77 for Galla and SEA. Therefore, the current breeding programs should be changed to avoid future genetic bottlenecks in this population. These findings offer an opportunity to enhance the current genetic status and management of Kenyan native goats and other regions with similar production systems.

Insights into the genomic homogeneity of Moroccan indigenous sheep breeds though the lens of runs of homozygosity.

Numerous studies have indicated that Morocco's indigenous sheep breeds are genetically homogenous, posing a risk to their survival in the challenging harsh climate conditions where they predominantly inhabit. To understand the genetic behind genetic homogeneity through the lens of runs of homozygosity (ROH), we analyzed the whole genome sequences of five indigenous sheep breeds (Beni Guil, Ouled Djellal, D'man, Sardi, Timahdite and Admixed).The results from principal component, admixture, Fst, and neighbour joining tree analyses consistently showed a homogenous genetic structure. This structure was characterized by an average length of 1.83 Mb for runs of homozygosity (ROH) segments, with a limited number of long ROH segments (24-48 Mb and > 48 Mb). The most common ROH segments were those ranging from 1-6 Mb. The most significant regions of homozygosity (ROH Islands) were mostly observed in two chromosomes, namely Chr1 and Chr5. Specifically, ROH Islands were exclusively discovered in the Ouled Djellal breed on Chr1, whereas Chr5 exhibited ROH Islands in all breeds. The analysis of ROH Island and iHS technique was employed to detect signatures of selection on Chr1 and Chr5. The results indicate that Chr5 had a high level of homogeneity, with the same genes being discovered across all breeds. In contrast, Chr1 displays some genetic variances between breeds. Genes identified on Chr5 included SLC39A1, IL23A, CAST, IL5, IL13, and IL4 which are responsible for immune response while genes identified on Chr1 include SOD1, SLAMF9, RTP4, CLDN1, and PRKAA2. ROH segment profile and effective population sizes patterns suggests that the genetic uniformity of studied breeds is the outcome of events that transpired between 250 and 300 generations ago. This research not only contributes to the understanding of ROH distribution across breeds but helps design and implement native sheep breeding and conservation strategies in Morocco. Future research, incorporating a broader sample size and utilizing the pangenome for reference, is recommended to further elucidate these breeds' genomic landscapes and adaptive mechanisms.

Buzzing with Intelligence: Current Issues in Apiculture and the Role of Artificial Intelligence (AI) to Tackle It.

Honeybees (Apis mellifera L.) are important for agriculture and ecosystems; however, they are threatened by the changing climate. In order to adapt and respond to emerging difficulties, beekeepers require the ability to continuously monitor their beehives. To carry out this, the utilization of advanced machine learning techniques proves to be an exceptional tool. This review provides a comprehensive analysis of the available research on the different applications of artificial intelligence (AI) in beekeeping that are relevant to climate change. Presented studies have shown that AI can be used in various scientific aspects of beekeeping and can work with several data types (e.g., sound, sensor readings, images) to investigate, model, predict, and help make decisions in apiaries. Research articles related to various aspects of apiculture, e.g., managing hives, maintaining their health, detecting pests and diseases, and climate and habitat management, were analyzed. It was found that several environmental, behavioral, and physical attributes needed to be monitored in real-time to be able to understand and fully predict the state of the hives. Finally, it could be concluded that even if there is not yet a full-scale monitoring method for apiculture, the already available approaches (even with their identified shortcomings) can help maintain sustainability in the changing apiculture.

A Review on Indigenous Goats of East Africa: A Case for Conservation and Management.

Indigenous goats are important in the livelihoods of rural households in East African countries. This is due to their ability to produce and reproduce in different environments and climatic conditions. Even though these indigenous goats are important, there is little available information on the genetic characterization of these breeds in Africa and at the international level. This paper reviews the status of indigenous goats, highlighting their production systems, phenotypic and genetic characteristics, and genetic diversity, and proposes potential ways for sustainable improvement and conservation in East African countries. Most households use traditional production systems with various uncharacterized goat breeds and ecotypes, which are hence named after the tribe or locality in which they are found. Most of these goats are classified as small East African breeds, with significant variability in morphological features. Some of the challenges to goat production in this region are indiscriminate crossbreeding, lack of pedigree records, parasites and disease incidences, low-quality pastures, and low levels of management. There is a need for a collaborative approach amongst the actors in goat breeding value chains as well as integrating modern genomic tools into breeding programs to enhance selection. This will ensure the resilience and sustainability of these unique indigenous goat populations in East Africa.

Climate change and dairy farming sustainability; a causal loop paradox and its mitigation scenario.

It is arguable at this time whether climate change is a cause or effect of the disruption in dairy farming. Climate change drastically affects the productive performance of livestock, including milk and meat production, and this could be attributed to the deviation of energy resources towards adaptive mechanisms. However, livestock farming also contributes substantially to the existing greenhouse gas pool, which is the causal of the climate change. We gathered relevant information from the recent publication and reviewed it to elaborate on sustainable dairy farming management in a changing climatic scenario, and efforts are needed to gather this material to develop methods that could help to overcome the adversities associated with livestock industries. We summarize the intervention points to reverse these adversities, such as application of genetic technology, nutrition intervention, utilization of chemical inhibitors, immunization, and application of metagenomics, which may help to sustain farm animal production in the changing climate scenario.

White and other fur colourations and hybridization in golden jackals (Canis aureus) in the Carpathian basin.

The golden jackal (Canis aureus) is a reoccurring species in the centre of the Carpathian basin, in Hungary. In total, 31 golden jackal tissue samples were collected, from 8 white-coated, 2 black-coated and one mottled animal across Hungary. Sequences and fragment length polymorphisms were studied for white colour (MC1R), and for black coat colouration (CBD103). In each white animal, the most widespread mutation causing white fur colour in dogs in homozygous form was detected. Three animals were found to carry the mutation in heterozygous form. The two black golden jackals were heterozygous for the 3 bp deletion in CBD103 that mutation for black coat colouration in dogs, and one of them also carried the mutation causing white fur. None of the white animals showed signs of hybridization, but both the black and the mottled coloured individuals were found to be hybrids based on genetic testing. Kinship was found three times, twice between white animals, and once between a white animal and an agouti animal carrying the mutation of white coat. Our results confirm the findings that golden jackal-dog hybrids may occur without human intervention, and the detected mutation causing white fur colour in golden jackals could possibly be due to an early hybridization event.

Hungarian indigenous Tsigai, a promising breed for excellent heat tolerance and immunity.

The adverse effects of climate change on sheep farming have become more noticeable in recent decades. Extensive efforts have been made to untangle the complex relationship between heat tolerance, animal health, and productivity, also to identify a resilient and economically suitable breed for selection that can be resilient to future climate change conditions. Using quantitative real-time polymerase chain reaction (qRT-PCR), we observed the seasonal variations in the expression of several important genes related to heat stress and immunity (HSP70, IL10, TLR2, TLR4, and TLR8) in three of the most widely kept sheep breeds in Hungary: The indigenous Tsigai, Hungarian Merino, and White Dorper. We found that the seasonal stressor affected the relative gene expression of all genes in this study. Notably, The Hungarian indigenous Tsigai was the most robust breed adapted to the Hungarian continental (hot summer, cold winter) environment, with excellent thermotolerance and immunity. Furthermore, despite suffering from heat stress in the summer, Hungarian Merino maintained their robust immune system well throughout the year.

Effect of ad libitum feeding of Holstein Friesian calves on immunological parameters and molecular stress on a transcriptional level.

Our aim was to compare the health and performance of ad libitum (ADLIB) and restrictedly fed Holstein Friesian heifer calves. Calves were selected to ADLIB (n = 13) and control (n = 13) groups randomly. The period of ADLIB feeding lasted for 3 weeks after colostrum supplementation. The calves in the control group received the same milk replacer, which was supplied according to the restrained feeding schedule of the farm. There was no difference between the two groups in weight, weight gain or movement activity, furthermore in the values of glucose, albumin, total protein, BHB, AST, oxidant and antioxidant status incl. dROM, PAT and OSI. The IL8 gene had higher levels (non significant, p > 0.05) of expression in the ADLIB group during the first 20 days of life, which indicates that ADLIB feeding might potentiate a stronger immune response to environmental stress. The IGF1 gene showed increased expression in the ADLIB fed group at almost all time points, however the difference was already detected on the first day of the study, indicating the importance of individual differences even within the same breed. During the first 10 days INS expression was higher in the restricted group, followed by a shift by day 20 and after, when the ADLIB group showed a higher relative expression level. The observed values describe a trend that, although not significant (p > 0.05), would seem to indicate that ADLIB feeding might potentiate a stronger immune response to environmental stress.

Genetic polymorphism in European and African sheep breeds reared in Hungary based on 48 SNPs associated with resistance to gastrointestinal parasite infection using KASP-PCR technique.

This pilot study used an alternative and economically efficient technique, the Kompetitive Allele-Specific Polymerase Chain Reaction (KASP-PCR) to examine 48 SNPs from 11 parasite-resistance genes found on 8 chromosomes in 110 animals from five sheep breeds reared in Hungary; Hungarian Tsigai, White Dorper, Dorper, Ile de France, and Hungarian Merino. Allele and genotype frequencies, fixation index, observed heterozygosity, expected heterozygosity, F statistic, and their relationship with the Hardy-Weinberg equilibrium (WHE) and the polymorphic information content (PIC) were determined, followed by principal component analysis (PCA). As much as 32 SNPs out of the 48 initially studied were successfully genotyped. A total of 9 SNPs, 4 SNPs in TLR5, 1 SNP in TLR8, and 4 SNPs in TLR2 genes, were polymorphic. The variable genotype and allele frequency of the TLRs gene indicated genetic variability among the studied sheep breeds, with the Hungarian Merino exhibiting the most polymorphisms, while Dorper was the population with the most SNPs departing from the HWE. According to the PIC value, the rs430457884-TLR2, rs55631273-TLR2, and rs416833129-TLR5 were found to be informative in detecting polymorphisms among individuals within the populations, whereas the rs429546187-TLR5 and rs424975389-TLR5 were found to have a significant influence in clustering the population studied. This study reported a moderate level of genetic variability and that a low to moderate within-breed diversity was maintained in the studied populations.

Kompetitive Allele Specific PCR Genotyping of 89 SNPs in Romanian Spotted and Romanian Brown Cattle Breeds and Their Association with Clinical Mastitis.

Mastitis is the most common production disease in the dairy sector worldwide, its incidence being associated with both cows' exposure to bacteria and the cows' genetic make-up for resistance to pathogens. The objective of our study was to analyse 89 missense SNPs belonging to six genes (CXCR2, CXCL8, TLR4, BRCA1, LTF, BOLA-DRB3), which were found to be associated with genetic resistance or susceptibility to mastitis. A total of 298 cattle (250 Romanian Spotted and 48 Romanian Brown) were genotyped by Kompetitive Allele Specific PCR (KASP) and a chi-squared test was used for genetic association studies with clinical mastitis. A total of 35 SNPs (39.3%) among the selected 89 SNPs were successfully genotyped, of which 31 markers were monomorphic. The polymorphic markers were found in two genes: TLR4 (rs460053411) and BOLA-DRB3 (rs42309897, rs208816121, rs110124025). The polymorphic SNPs with MAF > 5% and call rates > 95% were used for the association study. The results showed that rs110124025 in the BOLA-DRB3 gene was significantly associated with mastitis prevalence (p ≤ 0.05) in both investigated breeds. Current results show that the SNP rs110124025 in the BOLA-DRB3 gene can be used as a candidate genetic marker in selection for mastitis resistance in Romanian dairy cattle.

A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds.

BACKGROUND: To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. RESULTS: The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. CONCLUSIONS: To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

Honey bee (Apis mellifera) wing images: a tool for identification and conservation.

BACKGROUND: The honey bee (Apis mellifera) is an ecologically and economically important species that provides pollination services to natural and agricultural systems. The biodiversity of the honey bee in parts of its native range is endangered by migratory beekeeping and commercial breeding. In consequence, some honey bee populations that are well adapted to the local environment are threatened with extinction. A crucial step for the protection of honey bee biodiversity is reliable differentiation between native and nonnative bees. One of the methods that can be used for this is the geometric morphometrics of wings. This method is fast, is low cost, and does not require expensive equipment. Therefore, it can be easily used by both scientists and beekeepers. However, wing geometric morphometrics is challenging due to the lack of reference data that can be reliably used for comparisons between different geographic regions. FINDINGS: Here, we provide an unprecedented collection of 26,481 honey bee wing images representing 1,725 samples from 13 European countries. The wing images are accompanied by the coordinates of 19 landmarks and the geographic coordinates of the sampling locations. We present an R script that describes the workflow for analyzing the data and identifying an unknown sample. We compared the data with available reference samples for lineage and found general agreement with them. CONCLUSIONS: The extensive collection of wing images available on the Zenodo website can be used to identify the geographic origin of unknown samples and therefore assist in the monitoring and conservation of honey bee biodiversity in Europe.

Assessing the Genomics Structure of Dorper and White Dorper Variants, and Dorper Populations in South Africa and Hungary.

Dorper sheep was developed for meat production in arid and semi-arid regions under extensive production systems in South Africa. Two variants with distinct head and neck colors were bred during their development process. White Dorper have a white coat while Dorper have a black head and neck. Both variants have grown in popularity around the world. Therefore, understanding the genomic architecture between South African Dorpers and Dorper populations adapted to other climatic regions, as well as genomic differences between Dorper and White Dorper variants is vital for their molecular management. Using the ovine 50K SNP chip, this study compared the genetic architecture of Dorper variants between populations from South Africa and Hungary. The Dorper populations in both countries had high genetic diversity levels, although Dorper in Hungary showed high levels of inbreeding. White Dorpers from both countries were genetically closely related, while Dorpers were distantly related according to principal component analysis and neighbor-joining tree. Additionally, whereas all groups displayed unique selection signatures for local adaptation, Dorpers from Hungary had a similar linkage disequilibrium decay. Environmental differences and color may have influenced the genetic differentiation between the Dorpers. For their molecular management and prospective genomic selection, it is crucial to understand the Dorper sheep's genomic architecture, and the results of this study can be interpreted as a step in this direction.

Application potential of chicken DNA chip in domestic pigeon species - Preliminary results.

Introducing the SNP technology to pigeon breeding will enhance the competitiveness of a sector that produces one of the healthiest and best quality meats. The present study aimed to test the applicability of the Illumina Chicken_50K_CobbCons array on 24 domestic pigeon individuals from the Mirthys hybrids and Racing pigeon breeds. A total of 53,313 SNPs were genotyped. Principal component analysis shows a significant overlap between the two groups. The chip performed poorly in this data set, with a call rate per sample of 0.474 (49%). The low call rate was likely due to an increase in the evolutionary distance. A total of 356 SNPs were retained after a relatively strict quality control. We have demonstrated that it is technically feasible to use a chicken microarray chip on pigeon samples. Presumably, with a larger sample size and by assigning phenotypic data, efficiency would be improved, allowing more thorough analyses, such as genome-wide association studies.

An On-Station Analysis of Factors Affecting Growth Traits of Pure Red Maasai and Dorper Sheep Breeds under an Extensive Production System.

This comparative study aimed to evaluate genetic and environmental factors' effects on the growth traits of lambs in Dorper and Red Maasai (RedM) sheep breeds. The data analyzed contained the following measurements: birth weight (Bwt), weaning weight (Wwt), yearling weight (Ywt), birth type (single or twins), and details on each lamb's dam (dam ID and age) and sire. Except for the RedM, whose birth weight decreased with time, both breeds generally showed an increase in other growth weights across the study period, with yearly variations affecting both breeds. Additionally, analysis by a linear mixed model with restricted maximum likelihood (REML) showed that only breed as a factor significantly (p < 0.05) influenced birth weight; breed, sex, and birth type all significantly (p < 0.05) influenced weaning weight, whereas season, sex, and dam age significantly (p < 0.05) influenced yearling weight. The RedM breed outperformed the Dorper breed in post-weaning growth rate, demonstrating early resilience and adaptation to local environmental factors independent of maternal influence. Breed and non-genetic factors play a vital role in the growth of lambs, and the results of this study offer an opportunity for improved farm management under an extensive production system and selection for the conservation of the indigenous Red Maasai breed.

A review on the potential effects of environmental and economic factors on sheep genetic diversity: Consequences of climate change.

Climate change has a significant effect on the productivity of livestock including milk, meat, and reproduction. This could be attributed to the internal diversion of energy resources towards adaptive mechanisms. Among the climate change variables, thermal stress seems to be the major limiting factor in animal agriculture. A better understanding of the effects of climate change-influenced ecological factors on the genetic diversity of livestock species is warranted. Sheep is an ideal livestock species to be used in investigating environmental adaptation due to its wide range of agroecological habitats, genetic and phenotypic variability. There is a heavy reliance on sheep genetic diversity for future animal protein security, but the implications of climate change on their genetic diversity receive less attention. Here, the potential environmental factors influencing natural selection in sheep populations are presented. We argue that prolonged exposure to these factors plays a major role in influencing the development of adaptation traits in indigenous sheep breeds, consequently leading to the alteration of genetic diversity at specific loci. The factors discussed include hot temperatures (heat stress), insufficient water, low quantity and quality of forage, and prevalence of parasites, pests, and diseases. In addition, genetic diversity, some signatures of selection for adaptation and economic angles of selection are also briefly discussed. A better understanding of environmental factors influencing the genetic diversity of sheep populations will inform breeding and management programs and may offer an opportunity for greater production efficiency with low input costs.

Genetic diversity and substructuring of the Hungarian merino sheep breed using microsatellite markers.

The Hungarian Merino sheep breed (Ovis aries) is the most significant animal resource of the Hungarian sheep sector which, unfortunately, has gone through a huge reduction in number during the last decades and became endangered in 2014. A modern molecular genetic survey is now becoming more than necessary in order to characterize the within-breed genetic diversity and structure. For that reason, six Hungarian Merino flocks were genotyped in 16 microsatellite markers. In total, 144 different alleles were found and the mean values of observed and expected heterozygosity were 0.714 and 0.705, respectively, suggesting a noticeable genetic variability of the breed. The genetic differentiation of the Hungarian flocks was generally low, as reflected by the estimated total FST value (0.036), the extended pattern of admixture in Structure analysis, as well as, by the noticeable level of genetic clustering in UPGMA and FCA analyses. However, two out of the six studied flocks tended to be genetically more distant. The outcome of our study could be a starting point for a planned breeding strategy of the Hungarian Merino breed, by keeping the within-flock genetic variability in priority, as well as, by preserving the potential genetic uniqueness with close monitoring of the inbreeding.

Validation of SNP markers for thermotolerance adaptation in Ovis aries adapted to different climatic regions using KASP-PCR technique.

A study on 51 SNPs belonging to 29 genes related to heat stress was carried out in 720 sheep from 17 different breeds adapted to different climates from Hungary, Bosnia and Herzegovina, Morocco and Romania, using Kompetitive Allele-Specific Polymerase Chain Reaction. Genotype frequency and the Hardy-Weinberg equilibrium were calculated, followed by a clustering using the Principal Component Analysis. We analyzed the polymorphisms in the following genes analyzed: HSPA12A, HSP90AA1, IL33, DIO2, BTNL2, CSN2, ABCG1, CSN1S1, GHR, HSPA8, STAT3, and HCRT. We emphasized on HSPA12A and HSPA8 genes as they were successfully genotyped in all studied flocks in which genotype frequency patterns were identified. Contrary to previous findings, the A allele for HSPA8 SNP was not observed in the heat tolerant breeds, being found exclusively in cold-tolerant breeds. The principal component analysis could not clearly differentiate the breeds, while plot concentration was slightly varied among the three groups, with HSP90AA1 and IL33 SNPs' loading values significantly contributing to PC1 and PC2. We confirmed previous works that the HSPA12A, HSPA8, HSP90AA1 and IL33 SNPs are potential candidate markers for thermotolerance adaptation in sheep. This research contributes to the genetic variability of SNPs for thermotolerance adaptability in sheep.