Estimation of genomic inbreeding coefficient, detection of ROH Islands and Related Genes in different Egyptian sheep breeds adapted to different environment

Document Type : Research Paper

Authors

1 Department of Animal Sciences, Faculty of Agriculture and Natural Resources, Arak University, Arak, Iran

2 Department of Animal Sciences, Faculty of Agriculture, University of Ilam. Ilam, Iran.

Abstract

In this study in order to estimate inbreeding coefficient and identify the ROH Islands associated with the genes under selection, the 50k BeadChip genotyped data of 206 sheep from 3 different Egyptian breeds were used. After quality control, 48361 SNPs in 204 sheep were remained for the future analysis. Inbreeding coefficient was calculated using four methods including, genomic relationship matrix (FGRM), excess of homozygosity (FHOM), correlation between uniting gametes (FUNI) using the GCTA 1.0 software and run of homozygosity (FROH) using the PLINK 1.9 software. One percent of SNP with the highest frequency in ROH were considered as ROH Islands. The lowest rate of inbreeding according to (FGRM, FHOM, and FUNI) was related to Wahati and the highest was related to Saidi breed. The highest amount of FROH (0.043) was observed in Saidi and the lowest amount (0.018) was observed in Barki breed. Average length of ROH ranged from 45.02 to 205.87 Mb, while the average number of ROH ranged from 8.14 to 14.07. The highest number ROH was observed on chromosome 2, while the lowest was on chromosome 26. Average inbreeding coefficient from FROH in Barki, Saidi and Wahati were estimated 0.018, 0.043 and 0.027 respectively. A total of 62 ROH Islands with length: 24.60 Kb to 13 Mb were identified, which covering less than 1% of the sheep genome. The ROH Islands was not distributed across the genome uniform and varied among breeds, but some common were identified. Bioinformatics analysis demonstrated that some of these genomic regions overlapped with reported genes that directly or indirectly influenced traits for adaptation and immune system (SEMA3D, CSF2, ITPR1), development of the skeletal muscle (UGGT1, ITGA2), and reproduction (ABHD16B). The results of this study revealed that, the selection processes in different sheep breeds for economic traits during several years, has led to the formation of many ROH islands in sheep genome, therefore scanning these regions at the genome level can be an alternative strategy to identify genes and associated loci with economic traits.

Keywords

Main Subjects


Extended Abstract

Introduction

SNP genotyping technology is being used to study animal breed diversity and population structure. Understanding this diversity may lead to inferences about the interactions between different breeds that share geographic regions, including further information about the genetic evolution of animals and their historical origins in ancient parts of the world. Using genotyping technologies can help to identify animals with superior adaptive traits, such as the ability to withstand environmental stressors common in harsh environments, which will ensure the long-term viability of future breed improvement strategies. It can provide a genetic selection model for breeding programs to produce elite and well-adapted breeds. A run of homozygosity (ROH) is a consecutive tract of homozygous genotypes in an individual that indicates it has inherited the same ancestral haplotype from both parents. ROH one of the most methods was used to detecting the genomic inbreeding. The locations of ROHs which are under positive selection, or laboring favorable allele in population, tend to be fixed in the genome and formation of ROH Island during long times. Genomic regions enriched with ROH may be indicative of selection sweeps and are known as ROH islands. As detecting the ROH Islands, the genomic regions contain economic traits could be detectable. Three local sheep breeds (Saidi, Wahati, and Barki) were studied. The animals were exercised under natural heat stress. The heat tolerance index of the animals was calculated to identify animals with high and low heat tolerance based on their response to meteorological and physiological parameters.

 

Materials and methods

In this study in order to estimate inbreeding coefficient and identify the ROH Islands associated with the genes under selection, the 50k BeadChip genotyped data of 206 sheep from 3 different Egyptian breeds were used. Genomic variation in these breeds was assessed using 64,756 single nucleotide polymorphic markers (SNPs). All the animals were non-lactating, non-pregnant females ranging in parity from 2 to 4. A) The Barki desert sheep, a native of the CZWD known for its adaptation to desert environments, is a fat-tailed sheep with open, coarse wool, a white body, and a colored head. For centuries, they have been subject to natural selection in desert areas known for high temperatures, intense sunlight, frequent drought, and extensive grazing. B) Wahati sheep have a narrow, medium-sized body and a medium-sized head with a straight profile. Most rams are unhorned and have white fleece and small ears. The fat tail has a broad base that ends in a cylindrical part. The animals have a white body with a brown head or an all-white and C) Saidi sheep, the native Upper Egyptian breed, are bred in small flocks; they are known for their tolerance to heat stress and the prevailing hot, dry environment. The Saidi sheep are considered the oldest Egyptian sheep breed and are characterized by good fertility but high mortality of the young animals

After quality control, 48361 SNPs in 204 sheep were remained for the future analysis. inbreeding coefficient was calculated using four methods including, genomic relationship matrix (FGRM), excess of homozygosity (FHOM), correlation between uniting gametes (FUNI) using the GCTA 1.0 software and run of homozygosity (FROH) using the PLINK 1.9 software. One percent of SNP with the highest frequency in ROH were considered as ROH Islands.

 

Results and discussion

The lowest rate of inbreeding according to (FGRM, FHOM, and FUNI) was related to Wahati and the highest was related to Saidi breed. The highest amount of FROH (0.043) was observed in Saidi and the lowest amount (0.018) was observed in Barki breed. Average length of ROH ranged from 45.02 to 205.87 Mb, while the average number of ROH ranged from 8.14 to 14.07. The highest number ROH was observed on chromosome 2, while the lowest was on chromosome 26. Average inbreeding coefficient from FROH in Barki, Saidi and Wahati were estimated 0.018, 0.043 and 0.027 respectively. A total of 62 ROH Islands with length: 24.60 Kb to 13 Mb were identified, which covering less than 1% of the sheep genome. The ROH Islands was not distributed across the genome uniform and varied among breeds, but some common were identified. Bioinformatics analysis demonstrated that some of these genomic regions overlapped with reported genes that directly or indirectly influenced traits for adaptation and  immune system (SEMA3D, CSF2, ITPR1), development of the skeletal muscle (UGGT1, ITGA2), and reproduction (ABHD16B).

 

Conclusion

The results of this study revealed that, the selection processes in different sheep breeds for economic traits during several years, has led to the formation of many ROH islands in sheep genome, therefore scanning these regions at the genome level can be an alternative strategy to identify genes and associated loci with economic traits. Also, our findings contribute to the understanding of genetic diversity and population demography, and help design and implement breeding and conservation strategies for study sheep.

محمّدی؛ حسین، رأفت؛ عباس، مرادی شهر بابک؛ حسین، شجاع؛ جلیل و مرادی، محمد حسین (1397). برآورد ضریب همخونی ژنومی و اندازه مؤثر جمعیت در گوسفندان زندی با استفاده از تراشه متراکم نشانگری. مجله پژوهش و سازندگی، 119 (2)، 142-129.
موسی نژاد خبیصی؛ مژده، اسمعیلی زاده؛ علی و اسدی فوزی، مسعود (1401). بررسی میزان همخونی ژنومی در گوسفندان بومی ایران با استفاده از نشانگرهای متراکم (600k SNP). مجله پژوهش­های تولیدات دامی، 13 (35)، 167-158.
نصرتی؛ مریم و محمّد آبادی، محمد رضا (1401). پویش ژنومی جزایر همخونی و ژن­های مرتبط با آن در جمعیت گوسفندان دنیا. مجله پژوهش­های تولیدات دامی، 13 (36)، 146-137.
 
REFERENCES
Addo, S., Klingel, S., Thaller, G., & Hinrichs, D. (2021). Genetic diversity and the application of runs of homozygosity-based methods for inbreeding estimation in German White-headed Mutton sheep. PLoS One, 16, e0250608.
Aboul-Naga, A.M., Alsamman, A.M., El Allali, A., Elshafie, M.H., Abdelal, E.S., Abdelkhalek, T.M., & Abdelsabour, T.H. (2022). Genome-wide analysis identified candidate variants and genes associated with heat stress adaptation in Egyptian sheep breeds. Frontiers in Genetics, 13, 898522.
Almamun, H. A., Clark, S. A., Kwan, P. & Gondro, C. (2015). Genome-wide linkage disequilibrium and genetic diversity in five populations of Australian domestic sheep. Genetic Selection Evolution, 47, 90.
Curik, I., Ferenčaković, M. & Sölkner, J. (2014). Inbreeding and runs of homozygosity: a possible solution to an old problem. Livestock Science, 166, 26-34.
Elshazly, A., & Youngs, C.R. (2019). Feasibility of utilizing advanced reproductive technologies for sheep breeding in Egypt. Part 1. Genetic and nutritional resources. Egyptian Journal of Sheep and Goats Sciences, 14, 39–52.
Eydivandi. S., Roudbar, M. A., Karimi, M. O. & Sahana, G. (2021). Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia. Scientific Reports, 11(1), 2834.
Gomez-Raya, L., Rodríguez, C., Barragán, C. & Silió, L. (2015). Genomic inbreeding coefficients based on the distribution of the length of runs of homozygosity in a closed line of Iberian pigs. Genetic Selection Evolution, 47,81.
Liu, J., Shi, L., Li, Y., Chen, L., Garrick, D., Wang, L. & Zhao, F. (2021). Estimates of genomic inbreeding and identification of candidate regions that differ between Chinese indigenous sheep breeds. Journal of Animal Science and Biotechnology, 12(1), 95.
Machová, K., Marina, H., Arranz, J. J., Pelayo, R., Rychtářová, J., Milerski, M., Vostrý, L. & Suárez-Vega, A. (2023). Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms. Animal, 17(1), 100690.
Mastrangelo, S., M. Tolone, M.T. Sardina, G. Sottile, A.M. Sutera, R. Di Gerlando & Portolano, B. (2017). Genome-wide scan for runs of homozygosity identifies potential candidate genes associated with local adaptation in Valle del Belice sheep. Genetics Selection Evolution, 49, e84.
Mastrangelo, S., E. Ciani, M. T., Sardina, G. Sottile, F. Pilla & Portolano, B. (2018). Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds. Animal Genetics, 49, 71-81.
Mohammadi, H., Rafat, A., Moradi Shahrebabak, H., Shodja, J. & Moradi, M. H. (2018). Estimation of genomic inbreeding coefficient and effective population size in Zandi sheep breed using density SNP markers (50K SNPChip). Animal Sciences, 119,129-142. (In Persian)
Moosanezhad Khabisi, M., Esmailizadeh, A. & Asadi Fozi, M. (2022). Evaluation of genomic inbreeding rate in Iranian native sheep using dense SNP markers (600K). Research on Animal Production, 13(35),158-167. (In Persian)
Moioli, B., Pilla, F. & Ciani, E. (2015). Signatures of selection identify loci associated with fat tail in sheep. Journal of Animal Science, 93, 4660–4669.
Nosrati, M. & Mohammad Abadi, M. R. (2022). Genome-wide scans of ROH islands and related genes in the world's sheep populations. Research on Animal Production, 13 (36), 137-146 (In Persian).
Pasandideh, M., Gholizadeh, M. & Rahimi Mianji, G. (2020). Estimation of effective population size and genomic inbreeding coefficients in Baluchi sheep using genome wide single nucleotide polymorphisms (SNPs). Iranian Journal of Applied Animal Sciences (IJAS), 10(2), 287-295.
Peripolli, E., Munari, D. P., Silva, M. V. G. B., Lima, A. L. F., Irgang, R. & Baldi, F. (2017). Runs of homozygosity: current knowledge and applications in livestock. Animal Genetics, (3), 255-271.
Purfield, D. C., S. McParland, E. Wall & Berry, D. P. (2017). The distribution of runs of homozygosity and selection signatures in six commercial meat sheep breeds. PLoS One, 12, e176780.
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R. & Bender, D. (2007). PLINK: a toolset for whole-genome association and population-based linkage analysis. The American Journal of Human Genetics, 81, 559-575.
Srikanth, K., Lee, S. H., Chung, K. Y., Park, J. E., Jang, G. W., Park, M. R., Kim, N. Y., Kim, T. H. & Lim, D. (2020). A Gene-set enrichment and protein-protein interaction network-based gwas with regulatory snps identifies candidate genes and pathways associated with carcass traits in hanwoo cattle. Genes (Basel), 11(3), 316.
Tsartsianidou, V., Sánchez-Molano, E., Kapsona, V. V., Basdagianni, Z., Chatziplis, D., Arsenos, G., Triantafyllidis, A. & Banos, G. (2021). A comprehensive genome-wide scan detects genomic regions related to local adaptation and climate resilience in Mediterranean domestic sheep. Genetic Selection Evolution, 53(1), 90.
VanRaden, P. M. (2008). Efficient methods to compute genomic predictions. Journal of Dairy Science, 91(11), 4414–4423.
Wann A. (2012). A role for IFT88/the primary cilium in the inflammatory response to interleukin-1. Cilianal, 1:60-61.
Xu, Z. Y., Yang, H., Xiong, Y. Z., Deng, C. Y. & Li, F. E. (2010). Identification of three novel SNPs and association with carcass traits in porcine TNNI1 and TNNI2. Molecular Biology Reports, 37, 3609–3613.
Yang, L., Qin, Y. & Jian, C. (2021). Screening for core genes related to pathogenesis of Alzheimer’s disease. Frontiers Cell Development Biology, 9, 668738.
Yurchenko, A. A., Daetwyler, H. D. & Yudin, N. (2018). Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation. Scientific Reports, 8, 12984.
Zhao, F., Deng, T., Shi, L., Wang, W., Zhang, Q., Du, L. & Wang L. (2020). Genomic scan for selection signature reveals fat deposition in Chinese indigenous sheep with extreme tail types. Animals (Basel), 10(5), 773.