Abstract
Background
Natural and artificial selection tend to cause variability that contributes to shape the genome of livestock in a way that differentiates them among the animal kingdom. The particular aim here is to identify positive selection signatures with whole genome pooled-sequence data of Pakistani Teddy goat.
Methods and results
Paired-end alignment of 635,357,043 reads of Teddy goat with (ARS1) reference genome assembly was carried out. Pooled-Heterozygosity (Hp) and Tajima’s D (TD) are applied for validation and getting better hits of selection signals, while pairwise FST statistics is conducted on Teddy vs. Bezoar (wild goat ancestor) for genomic differentiation, moreover annotation of regions under positive selection was also performed. Hp score with − ZHp > 5 detected six windows having highest hits on Chr. 29, 9, 25, 15 and 14 that harbor HRASLS5, LACE1 and AXIN1 genes which are candidate for embryonic development, lactation and body height. Secondly, − ZTD value of > 3.3 showed 4 windows with very strong hits on Chr.5 & 9 which harbor STIM1 and ADM genes related to body mass and weight. Lastly, − ZFST < − 5 generated four strong signals on Chr.5 & 12 harbor LOC102183233 gene. Other significant selection signatures encompass genes associated with wool production, prolificacy and coat colors traits in this breed.
Conclusions
In brief, this study identified the genes under selection in Pakistani Teddy goat that will be helpful to refining the marker-assisted breeding policies and converging required production traits within and across other goat breeds and to explore full genetic potential of this valued species of livestock.
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11033-021-06726-x/MediaObjects/11033_2021_6726_Fig1_HTML.jpg)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11033-021-06726-x/MediaObjects/11033_2021_6726_Fig2_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11033-021-06726-x/MediaObjects/11033_2021_6726_Fig3_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11033-021-06726-x/MediaObjects/11033_2021_6726_Fig4_HTML.png)
Similar content being viewed by others
Data availability
Relevant data is available in the manuscript including supplementary files and supporting information. This article has already been submitted to a preprint platform distributed under a Creative Commons CC BY license having DOI number https://doi.org/10.20944/preprints202012.0421.v1.
Code availability
FST, TD scripts are available publicly, while \(Hp\) statistics was in-house written ruby script. Manhattan and other supplementary plots were plotted using Bioconductor qqMan package using R.
References
Naderi S et al (2008) The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals. Proc Natl Acad Sci USA 105:17659–17664. https://doi.org/10.1073/pnas.0804782105
Frantz LA, Bradley DG, Larson G, Orlando L (2020) Animal domestication in the era of ancient genomics. Nat Rev Genet 21:449–460. https://doi.org/10.1038/s41576-020-0225-0
Ejaz Wasti MSH, Asif M, Shoaib Malik M, Husnain Z, Khatoon S, Shah A, Qureshi A, Farooq O, Saddiq S, Mazhar N, Shujaat Malik Awan MAH, Arslan M, Syed Muhammad Azeem H, Anwar F, Faisal Shamim M (2018–2019) Pakistan economic survey (ed. by Economic Adviser’s Wing FD)
Kuthu Z, Javed K, Awan K, Ahmad N, Ahad A (2017) A study on phenotypic and genetic correlations between birth weight and other growth traits in Teddy goat. J Dairy Vet Anim Res. https://doi.org/10.15406/jdvar.2017.05.00158
Henkel J et al (2019) Selection signatures in goats reveal copy number variants underlying breed-defining coat color phenotypes. PLoS Genet 15:e1008536. https://doi.org/10.1371/journal.pgen.1008536
Saif R, Henkel J, Jagannathan V, Drögemüller C, Flury C, Leeb T (2020) The LCORL locus is under selection in large-sized Pakistani goat breeds. Genes 11:168. https://doi.org/10.3390/genes11020168
TIANGEN BIOTECH (BEIJING) CO., LTD. https://www.tiangen.com/en/?product1/t1/4.html.
Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. https://doi.org/10.1093/bioinformatics/btu170
Li H et al (2009) The sequence alignment/map format and SAMtools. Bioinformatics 25:2078–2079. https://doi.org/10.1093/bioinformatics/btp352
Rubin C-J et al (2012) Strong signatures of selection in the domestic pig genome. Proc Natl Acad Sci USA 109:19529–19536. https://doi.org/10.1073/pnas.1217149109
Kofler R, Pandey RV, Schlötterer C (2011) PoPoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq). Bioinformatics 27:3435–3436. https://doi.org/10.1093/bioinformatics/btr589
Korneliussen TS, Moltke I, Albrechtsen A, Nielsen R (2013) Calculation of Tajima’s D and other neutrality test statistics from low depth next-generation sequencing data. BMC Bioinform 14:1–14. https://doi.org/10.1186/1471-2105-14-289
Wang X et al (2016) Whole-genome sequencing of eight goat populations for the detection of selection signatures underlying production and adaptive traits. Sci Rep 6:38932. https://doi.org/10.1038/srep38932
Guo J et al (2018) Whole-genome sequencing reveals selection signatures associated with important traits in six goat breeds. Sci Rep 8:1–11. https://doi.org/10.1038/s41598-018-28719-w
Zhou C et al (2019) Genome-wide association study for milk protein composition traits in a Chinese Holstein population using a single-step approach. Front Genet 10:72. https://doi.org/10.3389/fgene.2019.00072
Shaffer L, Young T, Guess F, Bensmail H, León R (2008) Using R software for reliability data analysis. Int J Reliabil Appl 91:53–70
Raza SHA et al (2020) Bioinformatics analysis and genetic polymorphisms in genomic region of the bovine SH2B2 gene and their associations with molecular breeding for body size traits in qinchuan beef cattle. Biosci Rep. https://doi.org/10.1042/BSR20192113
Tao L, He X, Pan L, Wang J, Gan S, Chu M (2020) Genome-wide association study of body weight and conformation traits in neonatal sheep. Anim Genet 51:336–340. https://doi.org/10.1111/age.12904
Zonaed Siddiki A, Miah G, Islam M, Kumkum M, Rumi MH, Baten A, Hossain MA (2020) Goat genomic resources: the search for genes associated with its economic traits. Int J Genomics. https://doi.org/10.1155/2020/5940205
El-Halawany N, Zhou X, Al-Tohamy AF, El-Sayd YA, Shawky A, Michal JJ, Jiang Z (2016) Genome-wide screening of candidate genes for improving fertility in Egyptian native Rahmani sheep. Anim Genet. https://doi.org/10.1111/age.12437
Yurchenko AA et al (2019) High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia. BMC Genomics 20:294. https://doi.org/10.1186/s12864-019-5537-0
Nayeri S et al (2019) Genome-wide association analysis for β-hydroxybutyrate concentration in milk in Holstein dairy cattle. BMC Genet 20:58. https://doi.org/10.1186/s12863-019-0761-9
Zhang S et al (2019) Comparative transcriptomic and proteomic analyses identify key genes associated with milk fat traits in Chinese Holstein cows. Front Genet 10:672. https://doi.org/10.3389/fgene.2019.00672
Carty CL, Johnson NA, Hutter CM, Reiner AP, Peters U, Tang H, Kooperberg C (2012) Genome-wide association study of body height in African Americans: The women’s health initiative SNP health association resource (share). Hum Mol Genet 21:711–720. https://doi.org/10.1093/hmg/ddr489
Xu S-S, Li M-H (2017) Recent advances in understanding genetic variants associated with economically important traits in sheep (Ovis aries) revealed by high-throughput screening technologies. Front Agric Sci Eng 4:279–288
Nonneman D, Shackelford S, King D, Wheeler T, Wiedmann R, Snelling W, Rohrer G (2013) Genome-wide association of meat quality traits and tenderness in swine. J Anim Sci 91:4043–4050. https://doi.org/10.2527/jas.2013-6255
Santana M et al (2014) Genome-wide association study for feedlot average daily gain in Nellore cattle (Bos indicus). J Anim Breed Genet 131:210–216. https://doi.org/10.1111/jbg.12084
Jeong H et al (2015) Exploring evidence of positive selection reveals genetic basis of meat quality traits in Berkshire pigs through whole genome sequencing. BMC Genet 16:1–9. https://doi.org/10.1186/s12863-015-0265-1
Zhang L et al (2020) Genomic variants associated with the number and diameter of muscle fibers in pigs as revealed by a genome-wide association study. Animal 14:475–481. https://doi.org/10.3389/fgene.2019.00672
Gao Y et al (2017) CNV discovery for milk composition traits in dairy cattle using whole genome resequencing. BMC Genomics 18:265. https://doi.org/10.1186/s12864-017-3636-3
Zheng Z et al (2020) The origin of domestication genes in goats. Sci Adv 6:eaaz5216. https://doi.org/10.1126/sciadv.aaz5216
Acknowledgements
Authors are obliged to Prof. Dr. Tosso Leeb & Dr. Vidhya Jagannathan for their continuous support and valued time for accomplishing this work. Thanks to Dr. Mushtaq Ahmad Gondal and Mr. Muhammad Sidique for helping in sample collection.
Funding
Sequencing of Pakistani goats were funded by Swiss National Science Foundation (31003A_172964), when RS was postdoc fellow at University of Bern under Swiss Government Excellence Scholarship/Hans Sigrist Foundation and now further analyzing this already generated publicly available data.
Author information
Authors and Affiliations
Contributions
RS envisaged the idea, involved in critical thinking, analysis of data, editing, proofread and correspondence with journal. JH was involved in understandings the software and Linux features, analysis of data and provision of in-house \(Hp\) script. TM, AE and FA helped in data analysis and initial write-up. SZ helped in understanding the statistical methods and proof reading of the manuscript.
Corresponding author
Ethics declarations
Conflict of interest
There are no competing interests among the authors.
Ethical approval
All Pakistani animals were sampled with the consent of their owners by obeying the local regulations laid down by “Ethical Review Committee” UVAS-Pakistan, while sampling of Swiss animals was permitted by the “Cantonal Committee for Animal Experiments” (Canton of Bern; permit 75/16).
Consent to participation
Authors informed to the farmers and took their consent before collecting samples of their animals.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Saif, R., Henkel, J., Mahmood, T. et al. Detection of whole genome selection signatures of Pakistani Teddy goat. Mol Biol Rep 48, 7273–7280 (2021). https://doi.org/10.1007/s11033-021-06726-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11033-021-06726-x