Scientia Agricultura Sinica ›› 2014, Vol. 47 ›› Issue (6): 1190-1199.doi: 10.3864/j.issn.0578-1752.2014.06.015

• ANIMAL SCIENCE·VETERINARY SCIENCERE·SOURCE INSECT • Previous Articles     Next Articles

Genome-wide Detection of Selection Signature on Sunite Sheep

1、Institute of Animal Science and Veterinary Medicine, Chinese Academy of Agricultural Sciences, Beijing 100193;
2、College of Animal Science and Technology, Sichuan Agricultural University, Ya’an 625014, Sichuan   

  1. 1、Institute of Animal Science and Veterinary Medicine, Chinese Academy of Agricultural Sciences, Beijing 100193;
    2、College of Animal Science and Technology, Sichuan Agricultural University, Ya’an 625014, Sichuan
  • Received:2013-07-12 Online:2014-03-15 Published:2013-12-20

Abstract: 【Objective】Selection signatures are the selective footprints across the genome because of the effects of selection in the process of species under natural and artificial selection. It could not only reflect the effect of selection in new breed cultivation but also used as a method for QTL mapping which is correlated with economically important traits. Sunite sheep is one of the excellent indigenous sheep breeds distributed in Inner Mongolia through long time artificial selection-, which could adapt to the harsh Gobi natural environment condition. Selection signature detection could be used to search candidate genes due to positive selection and reveal genetic mechanism of economically important traits. Moreover, it could provide genetic evidences for traits undergoing positive selection for a long time in the breed formation process of Sunite sheep. 【Method】 The integrated haplotype score (iHS) method was used to detect genomic selection signatures in the Sunite sheep population based on the Illumina Ovine SNP50K BeadChip data. After quality control, 42616 SNP markers were retained for linkage disequilibrium (LD) analysis and haplotype construction. According to ancestral allelic information, 30537 SNP makers were left to calculate the iHS values. All the iHS values within one window with 500kb length, which was split the whole genome into non-overlapping segments, were averaged. After significant test, the genomic regions with selection signals were annotated. To investigate the biological function of candidate genes, gene ontology enrichment analysis was carried out. 【Result】 The LD decay map of Sunite sheep was constructed. Linkage disequilibrium analysis suggested that LD decreased with marker distances increase, but few pair markers with long distance had high LD levels. There were 204 genomic regions with selection signatures harboring 845 candidate genes were detected. For example, RXFP2 confers the absence of horns in sheep and ASIP could regulate a series of alleles for black and white coat color. HTR4 and SOX10 involved on nervous system development, SOX10 was closely related to neural crest development in embryonic period. E2F2 had an important effect in skeletal development which could activate transcription factor-2 in skeletal growth control. PLA2G6 could affect growth of skeleton and muscle. RPL7 promoted synthesis of ribosomal protein and POL was linked with Ovine Pulmonary Carcinoma. MATR3 could regulate innate and adaptive immune response. Furthermore, candidate genes also included ZWINT, PPP1R1B, GPR98, LUC7L3, CAPZA1 and MYT1L connected with the development of nervous system and disease traits. Bioinformatics analysis found 24 GO items in biological process, 4 GO items in molecular function and 2 GO items in cellular component, respectively. These GO items were mainly related to protein synthesis, macromolecular substances metabolic degradation, targeted transportation of protein, molecular activity, ribosomal components and ribosomal subunits. 【Conclusion】 The first LD map and genome-wide selection signatures map on Chinese sheep breeds were constructed. Many candidate genes related with economically important traits were found through genome-wide selection signatures detection. Some of them have very important significance in the process of sheep domestication. This study provided important reference for further understanding the process of artificial selection on sheep and also provided a theoretical basis for the breeding of Chinese sheep.

Key words: whole genome , linkage disequilibrium , selection signature , Sunite sheep

[1]国家畜禽遗传资源委员会,中国畜禽遗传资源志, 羊志, 北京: 中国农业出版社: 2011: 40-42.

China National Commission of Animal Genetic Resources. Animal Genetic Resources in China. Sheep and Goat. Beijing: Chinese Agricultural Press, 2011: 40-42. (in Chinese)

[2]Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 1989, 123(3):585-595.

[3]Fay J C, Wu C I. Hitchhiking under positive Darwinian selection. Genetics, 2000, 155(3):1405-1413.

[4]Lewontin R C, Krakauer J. Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics, 1973, 74(1):175-195.

[5]Sabeti P C, Reich D E, Higgins J M, Levine H Z, Richter D J, Schaffner S F, Gabriel S B, Platko J V, Patterson N J, McDonald G J, Ackerman H C, Campbell S J, Altshuler D, Cooper R, Kwiatkowski D, Ward R, Lander E S. Detecting recent positive selection in the human genome from haplotype structure. Nature, 2002, 419(6909):832-837.

[6]Sabeti P C, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie X, Byrne E H, McCarroll S A, Gaudet R, Schaffner S F, Lander E S, Frazer K A, Ballinger D G, Cox D R, Stewart J. Genome-wide detection and characterization of positive selection in human populations. Nature, 2007, 449(7164):913-918.

[7]Voight B F, Kudaravalli S, Wen X, Pritchard J K. A map of recent positive selection in the human genome. PLoS Biology, 2006, 4(3):e72.

[8]Qanbari S, Gianola D, Hayes B, Schenkel F, Miller S, Moore S, Thaller G, Simianer H. Application of site and haplotype-frequency based approaches for detecting selection signatures in cattle. BMC Genomics, 2011, 12:318.

[9]Browning B L, Browning S R. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. American Journal of Human Genetics, 2009, 84(2):210-223.

[10]Lewontin R C. The Interaction of Selection and Linkage. I. General Considerations; Heterotic Models. Genetics, 1964, 49(1):49-67.

[11]Hill W G. Estimation of linkage disequilibrium in randomly mating populations. Heredity, 1974, 33(2):229-239.

[12]Zhao H, Nettleton D, Dekkers J C. Evaluation of linkage disequilibrium measures between multi-allelic markers as predictors of linkage disequilibrium between single nucleotide polymorphisms. Genetical Research, 2007, 89(1):1-6.

[13]Kijas J W, Lenstra J A, Hayes B, Boitard S, Porto Neto L R, San Cristobal M, Servin B, McCulloch R, Whan V, Gietzen K, Paiva S, Barendse W, Ciani E, Raadsma H, McEwan J, Dalrymle B, other members of the International Sheep Genomics Consortium. Genome-wide analysis of the world's sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biology, 2012, 10(2):e1001258.

[14]Biswas S, Akey J M. Genomic insights into positive selection. Trends in Genetics : TIG, 2006, 22(8):437-446.

[15]Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics, 2009, 10:48.

[16]Farnir F, Coppieters W, Arranz J J, Berzi P, Cambisano N, Grisart B, Karim L, Marcq F, Moreau L, Mni M, Nezer C, Simon P, Vanmanshoven P, Wagenaar D, Georges M. Extensive genome-wide linkage disequilibrium in cattle. Genome Research, 2000, 10(2): 220-227.

[17]Smith E M, Wang X, Littrell J, Eckert J, Cole R, Kissebah A H, Olivier M. Comparison of linkage disequilibrium patterns between the HapMap CEPH samples and a family-based cohort of Northern European descent. Genomics, 2006, 88(4):407-414.

[18]Pullen T J, Sylow L, Sun G, Halestrap A P, Richter E A, Rutter G A. Overexpression of monocarboxylate transporter-1 (SLC16A1) in mouse pancreatic beta-cells leads to relative hyperinsulinism during exercise. Diabetes, 2012, 61(7):1719-1725.

[19]Laayouni H, Montanucci L, Sikora M, Mele M, Dall'Olio GM, Lorente-Galdos B, McGee KM, Graffelman J, Awadalla P, Bosch E, Comas D, Navarro A, Calafell F, Casals F, Bertranpetit J. Similarity in recombination rate estimates highly correlates with genetic differentiation in humans. PloS One, 2011, 6(3):e17913.

[20]Amaral A J, Ferretti L, Megens H J, Crooijmans R P, Nie H, Ramos-Onsins S E, Perez-Enciso M, Schook L B, Groenen M A. Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA. PloS One, 2011, 6(4):e14782.

[21]Giuffra E, Tornsten A, Marklund S, Bongcam-Rudloff E, Chardon P, Kijas J M, Anderson S I, Archibald A L, Andersson L. A large duplication associated with dominant white color in pigs originated by homologous recombination between LINE elements flanking KIT. Mammalian Genome: Official Journal of the International Mammalian Genome Society, 2002, 13(10):569-577.

[22]Stella A, Ajmone-Marsan P, Lazzari B, Boettcher P. Identification of selection signatures in cattle breeds selected for dairy production. Genetics, 2010, 185(4):1451-1461.

[23]Johnston S E, McEwan J C, Pickering N K, Kijas J W, Beraldi D, Pilkington J G, Pemberton J M, Slate J. Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population. Molecular Ecology, 2011, 20(12):2555-2566.

[24]Yuan F P, Li X, Lin J, Schwabe C, Bullesbach E E, Rao C V, Lei Z M. The role of RXFP2 in mediating androgen-induced inguinoscrotal testis descent in LH receptor knockout mice. Reproduction, 2010, 139(4):759-769.

[25]Gautier M, Naves M. Footprints of selection in the ancestral admixture of a New World Creole cattle breed. Molecular Ecology, 2011, 20(15):3128-3143.

[26]Gautier M, Flori L, Riebler A, Jaffrezic F, Laloe D, Gut I, Moazami-Goudarzi K, Foulley JL. A whole genome Bayesian scan for adaptive genetic divergence in West African cattle. BMC Genomics, 2009, 10:550.

[27]Mollaaghababa R, Pavan W J. The importance of having your SOX on: role of SOX10 in the development of neural crest-derived melanocytes and glia. Oncogene, 2003, 22(20):3024-3034.

[28]Montero J A, Giron B, Arrechedera H, Cheng Y C, Scotting P, Chimal-Monroy J, Garcia-Porrero J A, Hurle J M. Expression of Sox8, Sox9 and Sox10 in the developing valves and autonomic nerves of the embryonic heart. Mechanisms of Development, 2002, 118(1-2): 199-202.

[29]Liu J, Zhang L, Xu L, Ren H, Lu J, Zhang X, Zhang S, Zhou X, Wei C, Zhao F, Du L. Analysis of copy number variations in the sheep genome using 50K SNP BeadChip array. BMC Genomics, 2013, 14:229.

[30]Ebermann I, Wiesen M H, Zrenner E, Lopez I, Pigeon R, Kohl S, Lowenheim H, Koenekoop R K, Bolz H J. GPR98 mutations cause Usher syndrome type 2 in males. Journal of Medical Genetics, 2009, 46(4):277-280.

[31]Gao G, Dudley S C, Jr. RBM25/LUC7L3 function in cardiac sodium channel splicing regulation of human heart failure. Trends in Cardiovascular Medicine, 2013, 23(1):5-8.

[32]Qanbari S, Pimentel E C, Tetens J, Thaller G, Lichtner P, Sharifi A R, Simianer H. A genome-wide scan for signatures of recent selection in Holstein cattle. Animal Genetics, 2010, 41(4):377-389.

[33]Luvalle P, Ma Q, Beier F. The role of activating transcription factor-2 in skeletal growth control. The Journal of Bone and Joint Surgery American Volume, 2003, 85-A (Suppl. 2):133-136.

[34]Zhang L, Liu J, Zhao F, Ren H, Xu L, Lu J, Zhang S, Zhang X, Wei C, Lu G Zheng Y, Du L. Genome-wide association studies for growth and meat production traits in sheep. PloS One, 2013, 8(6):e66569.

[35]王宇, 刘淑英, 韩敏, 李建云,内源性绵羊肺腺瘤病毒NM株pol基因的克隆与序列分析,中国兽医杂志, 2007: 8-10.

Wang Y, Liu S, Han M, Li J. Cloning and sequence analysis of pol gene of the Inner Mongolia strain of endogenous jaagsiekte sheep retrovirus. Chinese Journal of Veterinary Medicine, 2007, 43:8-10. (in Chinese)

[36]Palmarini M, Sharp J M, Lee C, Fan H. In vitro infection of ovine cell lines by Jaagsiekte sheep retrovirus. Journal of Virology, 1999, 73(12):10070-10078.

[37]Wool I G, Chan Y L, Gluck A. Structure and evolution of mammalian ribosomal proteins. Biochemistry and Cell Biology = Biochimie et Biologie Cellulaire, 1995, 73(11-12):933-947.

[38]Lee Y J, Jeong S H, Hong S C, Cho B I, Ha W S, Park S T, Choi S K, Jung E J, Ju Y T, Jeong C Y, Kim J W, Lee C W, Yoo J, Ko G H. Prognostic value of CAPZA1 overexpression in gastric cancer. International Journal of Oncology, 2013, 42(5):1569-1577.

[39]Stevens S J, van Ravenswaaij-Arts C M, Janssen J W, Klein Wassink-Ruiter J S, van Essen A J, Dijkhuizen T, van Rheenen J, Heuts-Vijgen R, Stegmann A P, Smeets E E, Engelen J J. MYT1L is a candidate gene for intellectual disability in patients with 2p25.3 (2pter) deletions. American Journal of Medical Genetics Part A, 2011, 155A(11):2739-2745.

[40]Roxburgh R H, Marquis-Nicholson R, Ashton F, George A M, Lea R A, Eccles D, Mossman S, Bird T, van Gassen K L, Kamsteeg E J, Love D R. The p.Ala510Val mutation in the SPG7 (paraplegin) gene is the most common mutation causing adult onset neurogenetic disease in patients of British ancestry. Journal of Neurology, 2013, 260(5):1286-1294.
[1] TU YunJie,JI GaiGe,ZHANG Ming,LIU YiFan,JU XiaoJun,SHAN YanJu,ZOU JianMin,LI Hua,CHEN ZhiWu,SHU JingTing. Screening of Wnt3a SNPs and Its Association Analysis with Skin Feather Follicle Density Traits in Chicken [J]. Scientia Agricultura Sinica, 2022, 55(23): 4769-4780.
[2] JunYi GAI,JianBo HE. Major Characteristics, Often-Raised Queries and Potential Usefulness of the Restricted Two-Stage Multi-Locus Genome-Wide Association Analysis [J]. Scientia Agricultura Sinica, 2020, 53(9): 1699-1703.
[3] DongHua LI,YaWei FU,ChenXi ZHANG,YanFang CAO,WenTing LI,ZhuanJian LI,XiangTao KANG,GuiRong SUN. Genome-Wide Identification and Characterization of Transposable Elements in Xichuan Black-Bone Chicken [J]. Scientia Agricultura Sinica, 2020, 53(7): 1491-1500.
[4] WANG HaiGang,WEN QiFen,MU ZhiXin,QIAO ZhiJun. Population Structure and Association Analysis of Main Agronomic Traits of Shanxi Core Collection in Foxtail Millet [J]. Scientia Agricultura Sinica, 2019, 52(22): 4088-4099.
[5] WANG Jia, ZENG ZhaoQiong, LIANG JianQiu, YU XiaoBo, WU HaiYing, ZHANG MingRong. Development New Molecular Markers for Quantitative Trait Locus (QTL) Analysis of the Seed Protein Content Based on Whole Genome Re-Sequencing in Soybean [J]. Scientia Agricultura Sinica, 2019, 52(16): 2743-2757.
[6] WANG YuLong, HUANG BingYan, WANG SiYu, DU Pei, QI FeiYan, FANG YuanJin, SUN ZiQi, ZHENG Zheng, DONG WenZhao, ZHANG XinYou. Development and Characterization of Whole Genome SSR in Tetraploid Wild Peanut (Arachis monticola) [J]. Scientia Agricultura Sinica, 2019, 52(15): 2567-2580.
[7] REN YiYing, CUI Cui, WANG Qian, TANG ZhangLin, XU XinFu, LIN Na, YIN JiaMing, LI JiaNa, ZHOU QingYuan. Genome-Wide Association Analysis of Silique Density on Racemes and Its Component Traits in Brassica napus L. [J]. Scientia Agricultura Sinica, 2018, 51(6): 1020-1033.
[8] GAO BaoZhen, LIU Bo, LI ShiKai, LIANG JianLi, CHENG Feng, WANG XiaoWu, WU Jian. Genome-Wide Association Studies for Flowering Time in Brassica rapa [J]. Scientia Agricultura Sinica, 2017, 50(17): 3375-3385.
[9] FAN Hu, ZHAO Tuan-Jie, DING Yan-Lai, XING Guang-南, GAI Jun-Yi. Genetic Analysis of the Characteristics and Geographic Differentiation of Chinese Wild Soybean Population [J]. Scientia Agricultura Sinica, 2012, 45(3): 414-425.
[10] CHU Ming-xing,ZHANG Bao-yun,WANG Ping-qing,FANG Li,DI Ran,MA Yue-hui,LI Kui
. Polymorphic and Linkage Analysis of Microsatellite OarJL36 and FecB Gene in Sheep
[J]. Scientia Agricultura Sinica, 2009, 42(6): 2133-2141 .
[11]

. Identification of QTL for Length and Diameter of Porcine Limb Bones in a White Duroc× Erhualian F2 Resource Population
[J]. Scientia Agricultura Sinica, 2009, 42(4): 1378-1385 .
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!