全基因组选择信号揭示绒山羊绒毛性状相关的候选基因

doi:10.11843/j.issn.0366-6964.2020.12.009

摘要/Abstract

摘要： 旨在对辽宁绒山羊和内蒙古绒山羊的基因组选择信号进行筛选。本研究利用Illumina 50K的山羊芯片，对内蒙古绒山羊、辽宁绒山羊和黄淮山羊进行基因型检测，质控后获得50 010个共同SNPs位点。通过群体分化系数Fst和XP-EHH，以黄淮山羊为参考群体分别对内蒙古绒山羊和辽宁绒山羊进行选择信号检测，Fst和XP-EHH值的top 5%作为阈值。结果，利用Fst在绒山羊基因组中筛选到501个候选基因，利用XP-EHH在绒山羊基因组中筛选到145个候选基因。其中有12个SNPs在绒山羊基因组中受到较强选择。通过基因注释筛选到21个候选基因，包括EXOC4、ASIC2、PCDH9、RHBDD1、IRS1和PDE10A等。这些基因主要富集在PI3K-Akt信号通路、蛋白质消化吸收和松弛素信号通路等。研究结果发现，绒山羊在驯化过程中其基因组很多与毛囊相关的基因受到了强烈的选择。这些发现也有助于更好地了解绒山羊的选择进展，为中国绒山羊品种的种质资源保护和利用提供重要参考。

关键词: 内蒙古绒山羊, 辽宁绒山羊, Fst, XP-EHH

Abstract: The purpose of this study was to screen the genome-wide selection signals of Liaoning cashmere goat and Inner Mongolia cashmere goat. Based on the Illumina Caprine 50K SNP chip, Liaoning cashmere goat, Inner Mongolia cashmere goat and Huanghuai goat were genotyped, and the common 50 010 SNPs were obtained after quality control. Using population differentiation coefficients Fst and XP-EHH, Huanghuai goats were used as reference populations to detect selection signals for Inner Mongolia cashmere goats and Liaoning cashmere goats. The top 5% of Fst and XP-EHH values was used as the threshold. The result showed that 501 candidate genes were selected by Fst and 145 candidate genes were selected by XP-EHH in cashmere goats. Based on these SNPs, we detected 12 SNPs under strong selection by Fst and XP-EHH. After gene annotation, 21 candidate genes were identified, such as EXOC4, ASIC2, PCDH9, RHBDD1, IRS1 and PDE10A. These genes were found to be mainly enriched in PI3K-Akt signaling pathway, protein digestion and absorption and relaxin signaling pathways. The study indicated that genes associated with cashmere traits were subjected to strong artificial selection pressure during domestication of the cashmere goats. These findings also help us better understand the selection progress of cashmere goats and provide a new theoretical basis for the protection and utilization of germplasm resources of Chinese cashmere goat breeds.

Key words: Inner Mongolia cashmere goat, Liaoning cashmere goat, Fst, XP-EHH

中图分类号:

S827.2

金美林, 陆健, 费晓娟, 卢曾奎, 权凯, 储明星, 狄冉, 王慧华, 魏彩虹. 全基因组选择信号揭示绒山羊绒毛性状相关的候选基因[J]. 畜牧兽医学报, 2020, 51(12): 2991-3000.

JIN Meilin, LU Jian, FEI Xiaojuan, LU Zengkui, QUAN Kai, CHU Mingxing, DI Ran, WANG Huihua, WEI Caihong. Genome-wide Selection Signals Reveal Candidate Genes Associated with the Cashmere Traits of Cashmere Goats[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(12): 2991-3000.

参考文献 49

[1]	DAI B,ZHANG M,YUAN J L,et al.Integrative analysis of methylation and transcriptional profiles to reveal the genetic stability of cashmere traits in the Tβ4 overexpression of cashmere goats[J].Animals (Basel),2019,9(12):1002.
[2]	GAO Y,WANG X L,YAN H L,et al.Comparative transcriptome analysis of fetal skin reveals key genes related to hair follicle morphogenesis in cashmere goats[J].PLoS One,2016,11(3):e0151118.
[3]	LIU Y,WANG L L,LI X Y,et al.High-throughput sequencing of hair follicle development-related micrornas in cashmere goat at various fetal periods[J].Saudi J Biol Sci,2018,25(7):1494-1508.
[4]	SU R,FAN Y X,QIAO X,et al.Transcriptomic analysis reveals critical genes for the hair follicle of Inner Mongolia cashmere goat from catagen to telogen[J].PLoS One,2018,13(10):e0204404.
[5]	LIU B,GAO F Q,GUO J,et al.A microarray-based analysis reveals that a short photoperiod promotes hair growth in the Arbas cashmere goat[J].PLoS One,2016,11(1):e0147124.
[6]	国家畜禽遗传资源委员会组.中国畜禽遗传资源志·羊志[M].北京:中国农业出版社,2011.National Livestock and Poultry Genetic Resources Committee.Animal Genetic Resources in China·Sheep and Goats[M].Beijing:China Agricultural Press,2011.(in Chinese)
[7]	周占琴.中国绒山羊业发展现状、前景与对策[J].中国畜牧杂志,2008,44(4):42-45.ZHOU Z Q.The present situation and future and countermeasures of the fine-wool purpose goat development in China[J].Chinese Journal of Animal Science,2008,44(4):42-45.(in Chinese)
[8]	GUO J Z,TAO H X,LI P F,et al.Whole-genome sequencing reveals selection signatures associated with important traits in six goat breeds[J].Sci Rep,2018,8(1):10405.
[9]	WEIR B S,COCKERHAM C C.Estimating f-statistics for the analysis of population structure[J].Evolution,1984,38(6):1358-1370.
[10]	马云龙,张勤,丁向东.利用高密度SNP检测不同猪品种间X染色体选择信号[J].遗传,2012,34(10):1251-1260.MA Y L,ZHANG Q,DING X D.Detecting selection signatures on X chromosome in pig through high density SNPs[J]. Hereditas, 2012,34(10):1251-1260.(in Chinese)
[11]	ISLAM R,LI Y F,LIU X X,et al.Genome-wide runs of homozygosity,effective population size,and detection of positive selection signatures in six Chinese goat breeds[J].Genes,2019,10(11):938.
[12]	GUO J Z,TAO H X,LI P F,et al.Whole-genome sequencing reveals selection signatures associated with important traits in six goat breeds[J].Sci Rep,2018,8(1):10405.
[13]	LI X K,SU R,WAN W T,et al.Identification of selection signals by large-scale whole-genome resequencing of cashmere goats[J].Sci Rep,2017,7(1):15142.
[14]	KIJAS J W,ORTIZ J S,MCCULLOCH R,et al.Genetic diversity and investigation of polledness in divergent goat populations using 52088 SNPs[J].Anim Genet,2013,44(3):325-335.
[15]	PURCELL S,NEALE B,TODD-BROWN K,et al.PLINK:a tool set for whole-genome association and population-based linkage analyses[J].Am J Hum Genet,2007,81(3):559-575.
[16]	BARENDSE W,HARRISON B E,BUNCH R J,et al.Genome wide signatures of positive selection:the comparison of independent samples and the identification of regions associated to traits[J].BMC Genomics,2009,10:178.
[17]	刘雪雪,阿地力江·卡德尔,董坤哲,等.德保矮马X染色体选择信号筛选[J].畜牧兽医学报,2015,46(12):2161-2168.LIU X X,ADLIJIANG H,DONG K Z,et al.Detecting selection signatures on X chromosome of Debao pony[J].Acta Veterinaria et Zootechnica Sinica,2015,46(12):2161-2168.(in Chinese)
[18]	SABETI P C,REICH D E,HIGGINS J M,et al.Detecting recent positive selection in the human genome from haplotype structure[J].Nature,2002,419(6909):832-837.
[19]	李长青,尹俊,张燕军,等.内蒙古绒山羊与辽宁绒山羊皮肤毛囊周期性变化的比较研究[J].畜牧兽医学报,2005, 36(7):674-679.LI C Q,YIN J,ZHANG Y J,et al.Comparative study on skin and hair follicles cycling between Inner Mongolia and Liaoning cashmere goats[J].Acta Veterinaria et Zootechnica Sinica,2005,36(7):674-679.(in Chinese)
[20]	STEFAN B.Means for stimulation and activation of hair growth by il-15:US,2008069791A1[P]. 2008-03-20.
[21]	BAI W L,YIN R H,YIN R L,et al.IGF1 mRNA splicing variants in liaoning cashmere goat:identification,characterization,and transcriptional patterns in skin and visceral organs[J].Anim Biotechnol,2013,24(2):81-93.
[22]	PLIKUS M V,MAYER J A,DE LA CRUZ D,et al.Cyclic dermal BMP signalling regulates stem cell activation during hair regeneration[J].Nature,2008,451(7176):340-344.
[23]	LARAMIE J M,WILK J B,WILLIAMSON S L,et al.Polymorphisms near EXOC4 and LRGUK on chromosome 7q32 are associated with Type 2 Diabetes and fasting glucose;the NHLBI family heart study[J].BMC Med Genet,2008,9:46.
[24]	KIKUCHI S,NINOMIYA T,KAWAMATA T,et al.Expression of ASIC2 in ciliated cells and stereociliated cells[J].Cell Tissue Res,2008,333(2):217-224.
[25]	CHIARELLI N,CARINI G,ZOPPI N,et al.Transcriptome-wide expression profiling in skin fibroblasts of patients with joint hypermobility syndrome/ehlers-danlos syndrome hypermobility type[J].PLoS One,2016,11(8):e0161347.
[26]	ZHANG M M,MIAO F,HUANG R,et al.RHBDD1 promotes colorectal cancer metastasis through the Wnt signaling pathway and its downstream target ZEB1[J].J Exp Clin Cancer Res,2018,37(1):22.
[27]	PELEGRINELLI F F F,THIRONE A C P,GASPARETTI A L,et al.Early steps of insulin action in the skin of intact rats[J].J Invest Dermatol,2001,117(4):971-976.
[28]	DUAN C H,XU J H,SUN C M,et al.Effects of melatonin implantation on cashmere yield,fibre characteristics,duration of cashmere growth as well as growth and reproductive performance of Inner Mongolian cashmere goats[J].J An Sci Biotechnol, 2015,6(1):22.
[29]	王林枫,杨改青,杨耀胜,等.光照和褪黑激素对非生绒期绒山羊激素分泌和绒毛生长的影响[J].中国畜牧兽医, 2008, 35(12):29-32.WANG L F,YANG G Q,YANG Y S,et al.Effects of photoperiod and melatonin on endocrine and cashmere growth in cashmere goats in telogen[J].China Animal Husbandry & Veterinary Medicine,2008,35(12):29-32.(in Chinese)
[30]	GE W,WANG S H,SUN B,et al.Melatonin promotes cashmere goat (Capra hircus) secondary hair follicle growth:a view from integrated analysis of long non-coding and coding RNAs[J].Cell Cycle,2018,17(10):1255-1267.
[31]	SPIWOKS-BECKER I,WOLLOSCHECK T,RICKES O,et al.Phosphodiesterase 10A in the rat pineal gland:localization,daily and seasonal regulation of expression and influence on signal transduction[J].Neuroendocrinology,2011,94(2):113-123.
[32]	DESTEFANO G M,KURBAN M,ANYANE-YEBOA K,et al.Mutations in the cholesterol transporter gene ABCA5 are associated with excessive hair overgrowth[J].PLoS Genet,2014,10(5):e1004333.
[33]	FANTAUZZO K A,KURBAN M,LEVY B,et al.Trps1 and its target gene Sox9 regulate epithelial proliferation in the developing hair follicle and are associated with hypertrichosis[J].PLoS Genet,2012,8(11):e1003002.
[34]	CASTELA M,LINAY F,ROY E,et al.Igf1r signalling acts on the anagen-to-catagen transition in the hair cycle[J].Exp Dermatol,2017,26(9):785-791.
[35]	RAKOWIECKI S,EPSTEIN D J.Divergent roles for Wnt/β-catenin signaling in epithelial maintenance and breakdown during semicircular canal formation[J].Development,2013,140(8):1730-1739.
[36]	GENTILE P,GARCOVICH S.Advances in regenerative stem cell therapy in androgenic alopecia and hair loss:wnt pathway,growth-factor,and mesenchymal stem cell signaling impact analysis on cell growth and hair follicle development[J].Cells,2019,8(5):466.
[37]	RISHIKAYSH P,DEV K,DIAZ D,et al.Signaling involved in hair follicle morphogenesis and development[J].Int J Mol Sci,2014, 15(1):1647-1670.
[38]	HE N,DONG Z G,TAI D P,et al.The role of Sox9 in maintaining the characteristics and pluripotency of Arbas cashmere goat hair follicle stem cells[J].Cytotechnology,2018,70(4):1155-1165.
[39]	VELTRI A,LANG C,LIEN W H.Concise review:wnt signaling pathways in skin development and epidermal stem cells[J].Stem Cells,2018,36(1):22-35.
[40]	KIMURA R,WATANABE C,KAWAGUCHI A,et al.Common polymorphisms in WNT10A affect tooth morphology as well as hair shape[J].Hum Mol Genet,2015,24(9):2673-2680.
[41]	SULAYMAN A,TIAN K C,HUANG X X,et al.Genome-wide identification and characterization of long non-coding RNAs expressed during sheep fetal and postnatal hair follicle development[J].Sci Rep,2019,9(1):8501.
[42]	REALE S,YAHYAOUI M H,FOLCH J M,et al.Genetic polymorphism of the K-casein (CSN3) gene in goats reared in Southern Italy[J].Ital J Anim Sci,2005,4(1):97-101.
[43]	LEGRAND J M D,ROY E,ELLIS J J,et al.STAT5 activation in the dermal papilla is important for hair follicle growth phase induction[J].J Invest Dermatol,2016,136(9):1781-1791.
[44]	LEGRAND J M D.Mesenchymal-epidermal interactions in hair follicle cycling and regeneration[D].Brisbane:The University of Queensland,2015.
[45]	SHIN H Y,HENNIGHAUSEN L,YOO K H.STAT5-driven enhancers tightly control temporal expression of mammary-specific genes[J].J Mammary Gland Biol Neoplasia,2019,24(1):61-71.
[46]	CATAISSON C,OHMAN R,PATEL G,et al.Inducible cutaneous inflammation reveals a protumorigenic role for keratinocyte CXCR2 in skin carcinogenesis[J].Cancer Res,2009,69(1):319-328.
[47]	YU Y C,SU Y J,OPALENIK S R,et al.Short tail with skin lesion phenotype occurs in transgenic mice with keratin-14 promoter-directed expression of mutant CXCR2[J].J Leukoc Biol,2008,84(2):406-419.
[48]	KARANTH S,ZINKHAN E K,HILL J T,et al.FOXN3 regulates hepatic glucose utilization[J].Cell Rep,2016,15(12):2745-2755.
[49]	SCHIAVO G,BERTOLINI F,UTZERI V J,et al.Taking advantage from phenotype variability in a local animal genetic resource:identification of genomic regions associated with the hairless phenotype in Casertana pigs[J].Anim Genet,2018,49(4):321-325.