Unraveling candidate genes related to heat tolerance and immune response traits in some native sheep using whole genome sequencing data

• Published in: Small ruminant research Vol. 225; p. 107018
• Main Authors: Saadatabadi, Leila Mohammadipour, Mohammadabadi, Mohammadreza, Nanaei, Hojjat Asadollahpour, Ghanatsaman, Zeinab Amiri, Stavetska, Ruslana Volodymyrivna, Kalashnyk, Oleksandr, Kochuk-Yashchenko, Oleksandr Anatoliiovych, Kucher, Dmytro Mykolaiovych
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.08.2023
• Subjects: Adaptation; Candidate gene; Immune; Iranian; Sheep; Breed n; SRA; FST; Immune; Θπ; GCTA; BAM; PCA; GATK; SNP; NCBI; Sheep; Iranian; Candidate gene; AMSL; Adaptation
• ISSN: 0921-4488; 1879-0941
• DOI: 10.1016/j.smallrumres.2023.107018

In the sheep breeding industry, heat stress caused by the general increase in global temperature has become a significant issue, having both direct and indirect effects on animal health and productivity. The identification of effective pathways for responding to heat stress and immune responses, as well as the discovery of related genes, can improve the production and sustainability of genetic diversity in this area. In this study, we used nucleotide diversity (θπ) and FST statistical measurements to analyze the genomic data of the native Iranian sheep in order to find potential genes related to heat adaptation and immune response. We also compared the whole genome sequencing data of 29 indigenous Iranian sheep (the Afshari, Ghezel, Makuei, Moghani, Shal, Zel, Karakul, Grey-Shiraz, Baluchi, and Kermani breeds) with a number of other sheep breeds from Asia (n = 28), Europe (n = 28), and Africa (n = 25) to evaluate the genetic structure of the Iranian sheep population. The results from the sheep population genetic analysis showed a clear separation between different populations that corresponds well with their geographic origins. Iranian sheep were further divided into a southern and a northern group, which coincided with the most prominent climatic division in Iran. Our search to identify potential genomic regions under selection showed several candidate genes involved in (1) response to heat stress (SIK2, FER, ATP1A1, CDK5RAP3, and TLR4), (2) immune response in hot and dry environments (CD109, CR2, EOMES, and MARCHF1), (3) response to drought stress and adaptation to desert areas (ZFP1, PLCB1, and PDGFD), and also (4) response to heat stress by controlling digestive metabolism (HTR4, TRHDE, and ALDH1A3). The findings of this study may aid in our understanding of the molecular processes by which tiny ruminants adapt to hot, dry environments. In addition to the results of earlier studies, our findings also revealed a number of novel candidate genes related to heat adaptation. These genes will be valuable for future studies to choose livestock breeds that live in desert regions. •SIK2, FER, ATP1A1, CDK5RAP3, and TLR4 genes were identified as potential genes for the response to heat stress.•Candidate genes associated with the immune response in hot and dry environments included CD109, CR2, EOMES, and MARCHF1.•As potential genes for the response to drought stress and fitness to arid condition, ZFP1, PLCB1, and PDGFD were found.•HTR4, TRHDE, and ALDH1A3 genes were identified as candidates responding to heat stress by controlling digestive metabolism.