The pattern of runs of homozygosity, heterozygosity and genomic inbreeding in the population of sensitive and resistant sheep to footrot

Document Type : Research Paper

Authors

1 Department of Animal Science, Faculty of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University, Sari-Iran

2 Department of Animal Science, Faculty of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University, SARI-IRAN

Abstract

In this study, the genomic inbreeding coefficient based on runs of homozygosity (FROH) in the genomic regions of sheep sensitive and resistant to footrot was estimated. Swiss Alpine sheep data genotyped with Ovine SNP 600K BeadChip in a case-control study (154 healthy controls versus 76 severe cases) were used. FROH estimates ranged from 0.03 to 0.13, with a mean of 0.09. The highest FROH percentage was related to chromosome 2 (85.5%) and 1 (84.3%). The lowest FROH percentage was related to chromosome 8 (0.65%). 68959 runs of homozygosity (ROH) and 731026 runs of heterozygosity (ROHet) were detected in the healthy population and, 34569 ROH and 364125 ROHet in the case population. Using ROH, the highest percentage of observations of the healthy population was related to the 0-6 Mbp class (98.79%), and the lowest percentage of observations related to the class more than 48 Mbp (0.01%). The highest percentage of observations of the case population was related to Mbp 0-6 class (98.81%), and the lowest percentage of observations was related to 24-48 Mbp class. For ROHet, the highest percentage of observations of the healthy population was related to 0-6 Mbp class (99.99%), and the lowest percentage of observations was related to 6-12 Mbp class; for the case population, the lowest percentage was related to 6-12 Mbp and Mbp 24-48 classes. Totally, the healthy population had more number of homozygosity and heterozygosity pattern as well as the average length of heterozygosity pattern and a more ancient inbreeding than the case population.

Keywords

Main Subjects


Extended Abstract

Introduction

  The present study aimed to estimate the genomic inbreeding coefficient based on runs of homozygosity (FROH), to identify runs of homozygosity (ROH), and heterozygosity (ROHet), and to map selection signatures in the relevant genomic regions in Swiss white Alpine sheep.

 

Materials and Methods

  Swiss White Alpine sheep were genotyped using 600K single nucleotide polymorphism (SNP) arrays, and a case-control study was conducted in the form of 154 healthy controls versus 76 severe cases. For the quality control, SNPs that could not be assigned a position on the ARS-UI_Ramb_v2.0, and duplicate SNPs were discarded from further analysis. Individuals with missing genotypes larger than 0.05 were excluded. Also, SNPs with a call rate less than 0.05, a minor allele frequency (MAF) less than 0.05, and a significant deviation from Hardy-Weinberg equilibrium (HWE) (P < 1×10-6) were eliminated. 229 sheep and, a total of 461195 SNPs and 459100 SNPs, respectively, in the healthy and case groups met the desired criteria and were subsequently included in the final analysis.

 

Results and discussion

  FROH estimates ranged from 0.03 to 0.13, with an average of 0.09 for the entire genome. Chromosome 2 (85.5% by 65 counts) and 1 (84.3% by 129 counts) had the highest FROH percentage and chromosome 8 (0.65% by 1 count) had the lowest FROH percentage. In the healthy population, 68959 ROH and 731026 ROHet were observed, while in the case population, 34569 ROH and 364125 ROHet were detected. Based on the run lengths for the ROH and ROHet, five classes have been considered for each population (0-6, 6-12,12-24, 24-48, and >48 Mbp). Using ROH, the highest percentage of observations for healthy population was detected for the 0-6 Mbp class (98.79%,) and the lowest percentage of observations was for the class >48 Mbp, (0.01%). On the other hand, the highest percentage of observations in the case population was for the class 0-6 Mbp (98.81%), and the lowest percentage of observations was for the class 24-48 Mbp (0.01%). In the 0-6 Mbp class for the healthy population, there were 68128 ROH with an average size of 0.671 Mbp. For the case population, there were 34157 ROH with an average size of 0.680 Mbp. On the other hand, in the >48 Mbp class, the average size for the 7 ROH in the healthy population was 97.212 Mbp, and in the case population was 103.804 Mbp. For both healthy and severely case populations, chromosomes 1 and 26 had the highest and lowest number of ROH, respectively. The highest and lowest number of ROHet for both healthy and case populations were on chromosomes 1 and 24, respectively. There were 17 and 10 representative regions in ROH islands with frequencies >70% in the healthy and case populations, respectively. By ROHet, 6 and 19 regions presented frequencies of >60% in the healthy and the case populations, respectively. Using ROH analysis, the genes identified for healthy animals were distributed on chromosomes 1, 2, 8, 13, 15, 18, and 22, and for case animals on chromosomes 2, 4, 6, 8, 13, and 18. According to ROHet analysis, 527 genes had a known function and were distributed on chromosome 24 for both healthy and diseased populations. In total, 553 candidate genes were identified using ROH and ROHet analyses within the selection signatures. ROH islands had genes associated with reproductive traits (FSIP2), production traits (POPDC3), immune system (RAB39A), adaptation (PTPN9, LDB1), milk production (SIN3A), mastitis (MAN2C1) and disease resistance (COMMD4). ROHet islands had genes related to the immune system (MAPK8IP, ADCY9, IL32), inflammatory response (TNFRSF12A), adaptation (CLCN7, E4F1, CLDN9, HMOX2, ZNF598, PAQR4), body size (TRAF7, SEC14L5), body weight and body structure traits (GLIS2, VASN, TFAP4), litter size (NUDT16L1, ANKS3, ZNF500), milk traits (AMDHD2), calf birth weight (ZNF75A), lactation continuity (EEF2KMT, ALG1, SEC14L5, RBFOX1, NAGPA) ), growth traits (SOX8, SSTR5, IGFALS, NPW) and meat quality (RPL3L, SLC9A3R2). Using the ROH method, some of these genes identified in the healthy population were also identified in the case population such as FSIP2 on chromosome 2, POPDC3 on chromosome 8, and PTPN9, SIN3A, MAN2C1, and COMMD4 on chromosome 18. By the ROHet method, some of the genes identified in the healthy population in this study were consistent with the candidate genes reported in previous studies. including TEKT4, MGRN1, ROGDI, GLYR1, UBN1, PPL, UNKL, CCDC154, IFT140, DNAJA3, NMRAL1, HCFC1R1, THOC6, CLDN9, CLDN6, TNFRSF12A, BICDL2 and MMP25, MAPK8IP3, ADCY9, IL32, TNFRSF12A, CLCN7, E4F1, CLDN9, ZNF598, HMOX2, PAQR4, PTX4, TRAF7, SEC14L5, AMDHD2, ZNF75A, NAA60, NLRC3, TFAP4, GLIS2, VASN, NUDT16L1, ANKS3, ZNF500, EEF2KMT, ALG1, SEC14L5, RBFOX1, NAGPA, ADRA1D, SOX8 , SSTR5, IGFALS, NPW, RPL3L, SLC9A3R2 and RNF151, which were all located on chromosome 24. In the case population, GNG1 was detected. The results of the gene ontology (GO) showed that the most enriched terms identified in cellular components, molecular function, and KEGG pathway were related to the Apicolateral plasma membrane, Metal ion binding, and Fanconi anemia pathway, respectively.

 

Conclusion

  The results of the present study revealed that the healthy population of Swiss White Alpine sheep exhibited ancient inbreeding, a greater number of ROH and ROHet, and a higher average length of ROHet compared to the case population. The number of islands detected using the ROH was higher in the healthy population compared to the case population, whereas the ROHet method identified more islands in the case population than in the healthy population. According to the comparison made between two groups of case and healthy sheep, this research provides a new insight in understanding the mechanism of resistance to footrot in sheep.

Author Contributions

Conceptualization, M.G.; methodology, M.G. and F.E.; software, M.G. and F.E.; formal analysis, M.G. and F.E.; writing—original draft preparation, F.E.; writing—review and editing, M.G. and A.F.; supervision, M.G. and A.F. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The data that support the findings of this study can be retrieved at https://www.animalgenome.org/repository/pub/BERN2017.0821/.

 

Acknowledgements

We used genotyping data retrieved at public database (https://www.animalgenome). The authors are grateful for access to the data.

Ethical considerations

The study was approved by the Ethics Committee of the University of Sari Agricultural Sciences and Natural Resources University (Ethical code: SANRU.1403.03). The authors avoided data fabrication, falsification, plagiarism, and misconduct.

REFERENCES

Andreini, C., Banci, L., Bertini, I. & Rosato, A. (2006). Counting the zinc-proteins encoded in the human genome. Journal of Proteome Research, 5(1), 196-201.
Biscarini, F., Cozzi, P., Gaspa, G. & Marras, G. (2018). detectRUNS: Detect runs of homozygosity and runs of heterozygosity in diploid genomes.
Biscarini, F., Mastrangelo S., Catillo, G., Senczuk, G., & Ciampolini, R. (2020). Insights into genetic diversity, runs of homozygosity and heterozygosity-rich regions in Maremmana semi-feral cattle using pedigree and genomic data. Animals (Basel), 10, 2285.
Bjelland, D.W., Weigel, K.A., Vukasinovic N., & Nkrumah. J.D. (2013). Evaluation of inbreeding depression in Holstein cattle using whole genome SNP markers and alternative measures of genomic inbreeding. Journal of Dairy Science, 96, 4697-4706.
Burren, A., Flury, C., Aeschlimann, C., Hagger, C. & Rieder, S. (2012). Population structure and genetic diversity of Swiss sheep breeds. Agrarforschung Schweiz, 19(3), 140-147.
Carracelas, B., Navajas, E.A., Vera, B. & Ciappesoni, G. (2022). Genome-Wide Association Study of Parasite Resistance to Gastrointestinal Nematodes in Corriedale Sheep. Genes, 13(9), 1548.
Ceballos, F. C., Joshi, P. K., Clark, D. W., Ramsay, M., & Wilson, J. F. (2018). Runs of Homozygosity: Windows into Population History and Trait Architecture. Nature Review Genetics, 19, 220–234.
Curik, I., Ferenčaković, M., & Sölkner, J. (2014). Inbreeding and runs of homozygosity: A possible solution to an old problem. Livestock Sciience, 166, 26–34,
Dhungyel, O., Schiller, N., Eppleston, J., Lehmann, D., Nilon, P., Ewers, A. & Whittington, R. (2013). Outbreak-specific monovalent/bivalent vaccination to control and eradicate virulent ovine footrot. Vaccine, 31(13), 1701-1706.
Dzomba, E.F., Chimonyo, M., Pierneef, R. and Muchadeyi, F.C. (2021). Runs of homozygosity analysis of South African sheep breeds from various production systems investigated using OvineSNP50k data. BMC genomics, 22(1), 1-14.
Hamazaki, Y., Itoh, M., Sasaki, H., Furuse, M. & Tsukita S. (2002). Multi-PDZ domain protein 1 (MUPP1) is concentrated at tight junctions through its possible interaction with claudin-1 and junctional adhesion molecule. Journal of Biological Chemistry, 277, 455–61.
Hildersley, K.A., McNeilly, T.N., Gillan, V., Otto, T.D., Löser, S., Gerbe, F., Jay, P., Maizels, R.M., Devaney, E. & Britton, C. (2021). Tuft cells increase following ovine intestinal parasite infections and define evolutionarily conserved and divergent responses. Frontiers in Immunology, 12, 781108.
Huang, D.W., Sherman, B.T., & Lempicki, R.A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols, 4(1), 44-57.
Lanaspa, M.A., Almeida, N.E., Andres-Hernando, A., Rivard, C.J., Capasso, J.M. & Berl, T. (2007). The tight junction protein, MUPP1, is up-regulated by hypertonicity and is important in the osmotic stress response in kidney cells. Proceedings of the National Academy of Sciences, 104(34), 13672-13677.
Manunza, A., Diaz, J.R., Sayre, B.L., Cozzi, P., Bobbo, T., Deniskova, T., Dotsev, A., Zinovieva, N. & Stella, A. (2023). Discovering novel clues of natural selection on four worldwide goat breeds. Scientific Reports, 13(1), 2110.
Marras, G., Gaspa, G., Sorbolini, S., Dimauro, C., Ajmone-Marsan, P., Valentini, A., Williams, J. L., & Macciotta, N. P. P. (2015). Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy. Animal Genetics, 46, 110–121.
Marras, G., Wood, B. J., Makanjuola, B., Malchiodi, F., Peeters, K., Van As, P., Baes, C. F., Biscarini, F., & Turkeys, H. (2018). Characterization of runs of homozygosity and heterozygosity-rich regions in a commercial Turkey (Meleagris gallopavo) population. Proceedings of the World Congress on Genetics Applied to Livestock Production, 11, 763.
McQuillan, R., Leutenegger, A.L., Abdel-Rahman, R., Franklin, C.S., Pericic, M., & Barac-Lauc, L. (2008). Runs of homozygosity in European populations. Am. J. Hum. Gene, 83, 359–372.
Mucha, S., Bunger, L., & Conington. J. (2015). Genome-wide association study of footrot in Texel sheep. Genetics Selection Evolution, 47, 35.
Niggeler, A., Tetens, J., Stäuble, A., Steiner A., & Drögemüller C. (2017). A genome-wide significant association on chromosome 2 for footrot resistance/susceptibility in Swiss White Alpine sheep. Animal Genetics, 48(6):712-715.
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M.A., Bender, D., et al. (2007). PLINK: a tool set for whole genome association and population-based linkage analyses. American Journal of Humman Genetics, 81, 559–575.
Saif-ur-Rehman, M., Hassan, F.U., Reecy, J. & Deng, T. (2023). Whole-genome SNP markers reveal runs of homozygosity in indigenous cattle breeds of Pakistan. Animal Biotechnology, 34(4), 1384-1396.
Santos, W., Schettini, G., Fonseca, M. G., Pereira, G. L., Chardulo, L. A., Neto, O., et al. (2021). Fine‐scale Estimation of Inbreeding Rates, Runs of Homozygosity and Genome-wide Heterozygosity Levels in the Mangalarga Marchador Horse Breed. Journal of Animal Breeding and Genetics, 138, 161–173.
Santos, W. B., Pereira, C. B., Maiorano, A. M., Arce, C. D. S., Baldassini, W. A., Pereira, G. L., ... & Curi, R. A. (2023). Genomic inbreeding estimation, runs of homozygosity, and heterozygosity-enriched regions uncover signals of selection in the Quarter Horse racing line. Journal of Animal Breeding and Genetics.
Signer-Hasler, H., Burren, A., Ammann, P., Drögemüller, C. & Flury, C. (2019). Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds. Animal genetics, 50(5), 512-525.
Smedley, D., Haider, S., Ballester, B., Holland, R., London, D., Thorisson, G., & Kasprzyk, A. (2009). BioMart – biological queries made easy. BMC Genomics, 10, 22.
Solkner, J., Ferencakovic, M., Gredler, B., & Curik, I. (2010). Genomic metrics of individual autozygosity, applied to a cattle population, in: Proceedings of the 61st Annual Meeting of the European Association of Animal Production, Heraklion, Greece.
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. Genetics Selection Evolution, 53(1), 1-17.
Uimari, P., & Tapio, M. (2011). Extent of linkage disequilibrium and effective population size in Finnish Landrace and Finnish Yorkshire pig breeds. Journal of animal science, 89(3), 609-614.
Williams, J. L., Hall, S. J. G., Del Corvo, M., Ballingall, K. T., Colli, L., Ajmone Marsan, P., et al. (2016). Inbreeding and Purging at the Genomic Level: The Chillingham Cattle Reveal Extensive, Non-random SNP Heterozygosity. Animal Genetics, 47, 19–27.