畜牧兽医学报  2023, Vol. 54 Issue (11): 4577-4588. DOI: 10.11843/j.issn.0366-6964.2023.11.014    PDF    
永登七山羊全基因组选择信号检测分析
栗登攀1, 马克岩1, 韩金涛2, 白雅琴3, 李讨讨1, 马友记1     
1. 甘肃农业大学动物科学技术学院, 兰州 730070;
2. 永登县动物疫病预防控制中心, 兰州 730300;
3. 甘肃省畜牧技术推广总站, 兰州 730030
摘要:旨在检测永登七山羊群体基因组选择信号, 挖掘永登七山羊有价值的种质特性基因。本研究以4个绵羊群体(永登七山羊、岷县黑裘皮羊、兰州大尾羊和滩羊)共40个个体为研究对象, 利用简化基因组测序(specific-locus amplified fragment sequencing, SLAF-seq)技术检测全基因组范围内的单核苷酸多态性位点(SNPs)。基于SNPs数据集, 通过elgensoft软件进行主成分分析; 运用Treemix软件分析基因流事件; 利用群体遗传分化指数(Fst)和核苷酸多样性比值(π ratio)进行全基因组选择性清除分析, 取top 5% Fst和π ratio的交集以确定基因组受选择区域, 并对候选基因进行GO和KEGG富集分析。结果共得到1 658 596个群体SNPs; 主成分分析(PCA)发现永登七山羊能够独立分群, 基因流表明永登七山羊和兰州大尾羊存在较弱的基因交流。以永登七山羊为试验群体, 岷县黑裘皮羊、兰州大尾羊和滩羊为参考群体进行选择清除分析, 3个比较组的受选择区域分别检测出424、294、301个候选基因; GO和KEGG分析结果表明, 候选基因分别显著富集在65、79、41个GO条目及15、22、10条KEGG通路上(P<0.05)。此外, 将岷县黑裘皮羊、兰州大尾羊和滩羊3个群体的数据合并为一个数据集与永登七山羊进行比较, 共筛选到466个候选基因, 显著富集到124个GO条目及7条KEGG通路(P<0.05)。从中筛选到永登七山羊重要经济性状相关的功能基因BMP2、GRM1和ALDH1A1。研究结果表明, 在永登七山羊全基因组范围内进行选择信号检测, 鉴定到与生长发育、脂尾进化相关的候选基因, 为永登七山羊的分子遗传标记挖掘提供参考。
关键词永登七山羊    简化基因组测序    基因流    选择性清除分析    功能富集    
Detection and Analysis of Whole Genome Selection Signal of Yongdeng Qishan Sheep
LI Dengpan1, MA Keyan1, HAN Jintao2, BAI Yaqin3, LI Taotao1, MA Youji1     
1. College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730070, China;
2. Yongdeng Animal Disease Control Center, Lanzhou 730300, China;
3. Animal Husbandry Technology Extension Station of Gansu Province, Lanzhou 730030, China
Abstract: The aim of this study was to detect genomic selection signals in the Yongdeng Qishan sheep population and to mine valuable germplasm trait genes. A total of 40 individuals from 4 sheep populations (Yongdeng Qishan sheep, Minxian black fur sheep, Lanzhou fat-tailed sheep, Tan sheep) were used as study subjects to detect genome-wide single nucleotide polymorphism sites (SNPs) using specific-locus amplified fragment sequencing (SLAF-seq). Based on SNPs data, elgensoft software was used for principal component analysis (PCA), gene flow events were analyzed by Treemix software. Selective sweep analysis was performed using population genetic differentiation index (Fst) and nucleotide diversity ratio (π ratio) and the intersection of top 5% Fst and π ratio was taken to identify selected genomic regions. The GO and KEGG enrichment analysis of candidate genes was performed. A total of 1 658 596 population SNPs were obtained; Principal component analysis (PCA) found that Yongdeng Qishan sheep were able to independently group, gene flow showed that there was weak gene exchange between Yongdeng Qishan sheep and Lanzhou fat-tailed sheep. Yongdeng Qishan sheep were used as the experimental group, and Minxian black fur sheep, Lanzhou fat-tailed sheep and Tan sheep were used as the reference group for selection sweep analysis. The 424, 294 and 301 candidate genes were detected in the selected regions of the 3 comparison groups; GO and KEGG analyses showed that the candidate genes were significantly enriched in 65, 79 and 41 GO terms and 15, 22 and 10 KEGG pathways, respectively (P < 0.05). In addition, SNP datas from 3 populations of Minxian black fur sheep, Lanzhou fat-tailed sheep and Tan sheep were combined into one dataset for comparison with Yongdeng Qishan sheep and 466 candidate genes were significantly enriched to 124 GO terms and 7 KEGG pathways (P < 0.05).Functional genes BMP2, GRM1, and ALDH1A1 related to important economic traits were screened from Yongdeng Qishan sheep. The research shows that the selection signal detection in the whole genome of Yongdeng Qishan sheep can identify candidate genes related to growth and development and fat tail evolution, providing reference for molecular genetic marker mining of Yongdeng Qishan sheep.
Key words: Yongdeng Qishan sheep    specific-locus amplified fragment sequencing    gene flow    selective sweep analysis    functional enrichment    

我国地方绵羊品种资源丰富,为切实挖掘新资源并加以利用,在全国第三次畜禽资源普查中发现了永登七山羊。永登七山羊在永登县七山乡各村落均有分布,该群体具有特征一致、特性明显、遗传性能稳定等特点。研究发现,永登七山羊与兰州大尾羊、滩羊和岷县黑裘皮羊之间分化明显,独立形成聚类[1]

在动物漫长的驯化过程中,动物的行为习惯、体型外貌因受到自然选择而不断发生着巨大的变化[2-3]。使用选择信号来进行动植物的适应性进化研究越来越多[4-6]。选择信号是指动物在进化过程中受长期的自然选择和强烈的人工选择两种作用而遗留在基因组上的遗传标记,从基因组水平探讨种群间的遗传差异,有助于理解动物对选择适应的遗传机制和挖掘重要经济性状的候选基因,为种质资源的合理利用、保存提供理论依据[7-8]。Pan等[9]通过对89个中国地方绵羊品种进行全基因重测序揭示了RXFP2基因与角的独特形状相关。Li等[10]通过选择多种尾型的品种羊进行选择信号扫描,利用窗口XP-CLR和π ratio的方法,筛选出一系列(SGCZGLIS3、FGF7、PAPPA2、MAP2K3等)与脂肪沉积与代谢相关的基因。Li等[11]收集分析了中国、巴基斯坦和尼泊尔的15个代表性家山羊品种共115个样本以及4个野山羊样本的全基因组数据,筛选出调控山羊热适应、产绒和产奶的关键候选基因与SNP错义位点。这些研究结果有助于人们深刻了解动物受自然选择和人工选择的影响,并为物种起源及驯化机制提供见解。目前,对于永登七山羊的研究仅存在于遗传多样性分析上,需要进一步从基因组水平筛选鉴定可能与永登七山羊重要经济性状相关的功能基因。

本研究利用SLAF技术检测4种绵羊群体全基因组范围内的SNPs。通过主成分分析、基因流事件以及候选基因挖掘等方法,旨在从分子水平证明永登七山羊群体与其他群体的基因交流情况,并对适应性进化过程中的优势基因进行挖掘,从而为永登七山羊的保护利用提供重要的理论参考。

1 材料与方法 1.1 试验材料

本研究基于前期针对地方品种(永登七山羊)的遗传多样性分析研究[1],选择兰州大尾羊(LFT)、滩羊(TAN)、岷县黑裘皮羊(MBF)作为参考对象。如表 1所示,每个品种采集10只健康的成年母羊,颈部静脉采血3~5 mL,用5 mL EDTA抗凝管收集,置于-80 ℃冰箱保存备用。

表 1 样品信息 Table 1 Sample information
1.2 试验方法

对抗凝血样提取DNA,利用NanoDrop 2000和Qubit 2.0检测DNA的纯度和浓度。采用琼脂糖凝胶电泳检测基因组DNA的完整性。对DNA质检合格的样品采用双酶切的RAD建库方式, 将长度范围414~464 bp的酶切片段定义为SLAF标签,对SLAF标签进行处理后在Illumina HiSeqTM 2500测序仪上进行简化基因组测序。测序由北京百迈客生物科技有限公司完成。

1.3 数据处理与分析

对测序得到的原始数据按照酶切效率、测序质量和片段选择进行数据质量评估,得到各个样品的原始序列片段。其中限制性内切酶组合为Hpy166 Ⅱ+EcoR V,采用读长126 bp×2作为后续的数据评估和分析数据。利用bwa[12]将原始数据比对到绵羊参考基因组(oar_v4.0)上,使用GATK[13]和samtools[14]两种方法开发SNP,取SNP标记交集并使用snpeff软件[15]用于注释变异。

1.4 主成分分析

使用EIGENSOFT[16]软件包中的smart PCA程序基于SNP数据,在PCA二维聚类图[1]基础上重新进行数据分析,得到4个绵羊群体的三维聚类情况。

1.5 基因流分析

利用Treemix软件,使用全基因组范围内等位基因频率数据,推断多个群体分裂和混群模式。使用R的OptM判断最优m值,m值取1~10。每个m值重复5次后,使用最优m值进行Treemix分析,绘制群体最大似然树;并在最大似然树上用箭头展示渗入的迁移事件。

1.6 选择清除分析

基于SNP数据,按照次要等位基因频率(MAF:0.05)和位点完整度(INT:0.8)进行过滤,获得高一致性SNP位点,按照100 kb的窗口、10 kb的步长对染色体进行选择性清除区域检测。使用PopGenome软件包计算群体分化指数(Fst)、核苷酸多态性(π)和多样性变化倍数(θπ ratio)等群体遗传学指标。在筛选选择性清除区域时,分别以Fst和π ratio前5%分位数对应的数据作为阈值,取该区域的交集为选择性清除区域。

1.7 受选择区域筛选及基因富集分析

针对选择性清除区域进行相邻区域窗口合并,统计合并后窗口的位置信息以及窗口内基因数目,对应的基因位置。使用Enesmbl数据库的BioMart功能,利用功能分类和代谢途径信息对候选基因进行基因功能注释。利用DAVID数据库对基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)富集分析,并分别生成富集分析图表。

2 结果 2.1 SNP检测

本项目共获得709 067个SLAF标签,标签的平均测序深度为17.63x,其中,多态性SLAF标签有445 768个,共获得1 658 596个SNPs标记。根据SLAF在染色体上的分布,绘制SLAF在染色体上的分布图(图 1)。

图 1 SLAF在染色体上的分布图 Fig. 1 Distribution of SLAF on chromosomes
2.2 主成分分析

根据主成分分析推断出群体之间聚类模式的差异。由4个绵羊群体的主成分分析(图 2)可明显看出,PC1、PC2和PC3的贡献率分别为3.86%、3.43%和3.39%;LFT群个体最为聚集,MBF群体最为分散;MBF、QS和TAN群个体可独立分群。

图 2 4个绵羊群体的主成分分析 Fig. 2 Principal component analysis of 4 sheep populations
2.3 基因流分析

基因流动是指位于不同种群之间或者在种群内部产生的遗传物质的交流[17-18],既能抵消由遗传变异和自然选择引起的自我种群之间的分化,又能给整个大的群体中带来新的变异。从图 3可以看出,TAN和LFT群体间存在较强的基因交流,QS和LFT群体间发生较弱的基因交流。

图 3 基因流动 Fig. 3 Gene flow
2.4 群体间受选择区域筛选

将QS与LFT、TAN和MBF进行两两比较并将其余3个群体数据合为整体与QS进行比较。在QS与MBF之间,将Fst>0.273, π ratio>8.55的基因组交集区域定义为QS与MBF群体间选择信号(图 4A),检测到163个受选择窗口,注释到424个候选基因(图 5A);以同样的方法,在QS与LFT之间筛选到127个受选择窗口(图 4B),注释到294个候选基因(图 5B);在QS与TAN之间筛选到133个受选择窗口(图 4C),注释到301个候选基因(图 5C)。

A.QS与MBF组;B.QS与LFT组;C.QS与TAN组;D.QS与MLT组 A.QS and MBF group; B.QS and LFT group; C.QS and TAN group; D.QS and MLT group 图 4 群体间选择信号示意图 Fig. 4 Schematic of selection signatures between populations
A.QS与MBF组;B.QS与LFT组;C.QS与TAN组;D.QS与MLT组 A.QS and MBF group; B.QS and LFT group; C.QS and TAN group; D.QS and MLT group 图 5 两种方法检测的正选择基因韦恩图 Fig. 5 Venn diagram for the positively selected genes identified via two approaches

在QS与(MBF +LFT+TAN)(简称:MLT)之间,筛选到202个受选择窗口,并绘制群体间选择信号示意图(图 4D);注释到466个候选基因(图 5D)。另外,统计窗口中各个染色体的Fst和π ratio值,以及窗口内各个基因的位置,绘制QS与MLT SNPs注释表(表 2)。

表 2 QS与MLT组SNPs注释 Table 2 Annotations for SNPs in QS and MLT groups
2.5 GO富集分析

对候选基因进行GO功能富集分析。在QS与MBF组间424个候选基因中,119个候选基因显著富集于65个GO条目(P<0.05),其中包括9条细胞组分类,43条生物过程类和13条分子功能类,主要集中在细胞极性的建立或维持、肺泡发育、组蛋白H3-K4三甲基化、血管重塑等。在QS与LFT组294个候选基因中,119个候选基因显著富集于79个GO条目(P<0.05),其中包括54条生物过程类,18条分子功能类和7条细胞组分类,主要富集在mRNA稳定化、前置/后置模式规范、有丝分裂期板块聚集等。在QS与TAN组301个候选基因中,90个基因显著富集于41个GO条目(P<0.05),其中包括33条生物过程类和8条分子功能类,主要富集在烟酸和烟酰胺的代谢、金黄色葡萄球菌感染、雌性激素信号传导途径、初级胆汁酸的生物合成等。

另外,在QS与MLT组466个候选基因中,130个候选基因显著富集在124个GO条目(P<0.05),其中包括96条生物过程类,14条细胞组分类和14条分子功能类,主要富集在多细胞进程、血液循环、发育生长调节、钙离子转运的正调节等过程。富集条目、通路名称、P值以及对应基因见表 3

表 3 GO功能富集分析 Table 3 GO function enrichment analysis
2.6 KEGG富集分析

对候选基因进行KEGG通路富集分析,并挑选前10个显著的KEGG通路作富集气泡图(图 6)。结果显示, QS群体与MBF、LFT、TAN品种相比较,分别有231条、222条和214条信号通路参与主要生化代谢途径和信号转导途径。其中,QS与MBF比较组中显著富集通路有15条,其中15个候选基因(GAD1、CYP11B1、PTAFRCAMK2D、CDH1、BMP2等)可能涉及重要经济性状,主要参与类固醇激素生物合成、叶酸的生物合成、雌激素信号传导途径、金黄色葡萄球菌感染、卵巢类固醇生成、花生四烯酸代谢、精氨酸和脯氨酸的代谢等通路。QS与LFT比较组中显著富集通路有22条,涉及12个候选基因(ALDH1A1、GAD1、GRM1、PTAFRCYP11B1等),主要集中在雌性激素信号传导途径、金黄色葡萄球菌感染、视黄醇代谢、丙氨酸、天门冬氨酸和谷氨酸的代谢、苯丙氨酸的代谢、色氨酸的代谢、类固醇激素的生物合成等。QS与TAN比较组中显著富集通路有10条,涉及13个候选基因(为ALDH1A1、KMT2A、PARD6G、SETD1A、CAMK2D、GRM1等),主要集中在烟酸和烟酰胺的代谢、金黄色葡萄球菌感染、雌性激素信号传导途径、初级胆汁酸的生物合成、氮素代谢、粘着剂的交界处、赖氨酸的降解、视黄醇代谢、TGF-beta信号传导途径、轴突引导、嘧啶的代谢等。

A.QS vs MBF组;B.QS vs LFT组;C.QS vs TAN组;D.QS vs MLT组 A.QS and MBF group; B.QS and LFT group; C.QS and TAN group; D.QS and MLT group 图 6 KEGG通路分析 Fig. 6 KEGG pathway analysis

QS与MLT比较组中,参与主要生化代谢途径和信号转导途径有244条。其中PTAFRUGT2A1、KRT15等基因主要参与类固醇激素生物合成、叶酸合成、雌激素信号通路、金黄色葡萄球菌感染、卵巢类固醇生成、甘露糖型O-聚糖生物合成、SNARE在膀胱运输中的相互作用等。

值得注意的是,雌激素信号通路和金黄色葡萄球菌感染通路均显著存在于4个组别中,并从中均发现V15、K38、KRT32、KRT35、KRT36这5个基因。

3 讨论

在动物遗传学上,基于不同群体的SNP差异程度,可将各群体进行分群。本研究中,通过PCA分析,发现岷县黑裘皮羊、永登七山羊和滩羊均可与其他群体分开,除永登七山羊外,其他群体与之前一些学者的研究结果一致[19-20]

基因交流可以向其他群体中引入新的等位基因,进而影响群体遗传多样性,产生新的性状组合,发生遗传变异[21-23]。本研究发现,滩羊和兰州大尾羊之间发生过相对较强的基因交流,永登七山羊与兰州大尾羊之间发生了较弱的基因交流。其原因是滩羊、兰州大尾羊、永登七山羊均属蒙古羊,虽然永登七山羊饲养环境相对封闭,但因为和兰州大尾羊地理上接近,难免发生基因交流;也是自然选择作用的表现。因此保护这些品种生物多样性具有很大难度,可以采用等量留种的方式,尽量保持现有等位基因,防止近交衰退。

本研究取top 5% Fst和π ratio的交集用来确定基因组受选择区域,并对筛选出来的候选基因进行GO和KEGG富集分析。采用Fst和π ratio的交集来进行筛选鉴定,能够得到更加可靠的结果[24-25]。本试验中,各比较组的受选择窗口具有一定差异,且注释到的基因均随着受选择窗口的数量而变化。通过SNPs注释,发现了可能与永登七山羊重要经济性状相关功能的20个候选基因,包括CDH1、V15、K38、KRT15、KRT32、KRT35、KRT36、BMP2、GRM1、ALDH1A1、CAMK2D、MASP1、LUC7L3、UGT2A1、PTAFRPLOD1、SETD1A、PARD6G、KMT2A、SOD1。其中与绵羊生长发育有关的功能基因有BMP2、GRM1和ALDH1A1。在MBF与QS比较组中筛选出了BMP2基因,发现BMP2基因存在于Hippo信号通路、基底细胞癌、转化生长因子-β(TGF-β)和细胞因子-细胞因子受体相互作用的信号通路中。BMP2基因影响羊生长性状和脂肪酸组成[26-27],在肌肉发育中发挥重要作用[28]。与脂肪尾表型,脂质代谢相关,参与了羊尾脂肪沉积[29-31]。除此之外,BMP2可能是影响绵羊产仔数的候选基因[32]。在滩羊与永登七山羊比较组和兰州大尾羊与永登七山羊比较组均筛选出GRM1和ALDH1A1基因。发现GRM1基因存在于雌激素信号通路、磷脂酶D信号通路、逆行内源性大麻素信号、长时程增强、长时程抑制、间隙连接的信号通路,GRM1基因被认为是影响断奶后增益的重要候选基因,对绵羊生长和肉类生产性状有一定影响[33-35], 可以作为绵羊肉质性状的候选基因[36];又发现ALDH1A1基因存在于视黄醛代谢信号通路。ALDH1A1基因是绵羊脂尾进化关键基因[37],可能参与滩羊脂肪细胞分化调控,敲除滩羊尾部脂肪前体细胞中ALDH1A1基因能够抑制脂肪细胞分化过程[38]

另外,雌激素信号通路均显著存在于4个组别中。雌激素介导的oar-miR-485-5p靶向PPP1R13B调控绵羊成肌细胞增殖,参与炎症反应、新陈代谢和干细胞存活[39-40]。雌激素与雌激素受体结合后可以调节细胞核中的基因转录, 或者激活细胞质中的激酶进而发挥其作用[41],且雌激素介导的信号通路与毛囊周期性生长和发育有关[42]。V15、K38、KRT32、KRT35、KRT36均作为角蛋白,广泛存在于绵羊毛发中[43-44]。其中KRT36是影响羊毛弯曲性状的潜在候选基因[45]。在毛球下部至中部,KRT32和KRT35基因在角质层和纤维皮质层均有表达[46]KRT32在毛球中部至角质生长层皮质中广泛表达,KRT35基因和毛囊透明蛋白参与初始分化的皮质层结构的建立, 最终导致纤维弯曲[47-48]。本研究采用简化基因组测序检测永登七山羊选择信号,注释到大量候选基因。在候选基因功能显著性富集分析中,功能基因揭示了永登七山羊在尾部脂肪沉积、肌肉发育和毛发生长等方面受到了选择。

4 结论

本研究通过主成分分析和基因交流揭示4个群体的进化关系,结果表明,永登七山羊群体与其他3个群体分化明显,且与兰州大尾羊发生基因交流;此外还筛选到与绵羊生长发育有关的功能基因BMP2、GRM1和ALDH1A1等。结果为挖掘永登七山羊生长发育等相关基因提供参考,为其遗传资源保护与开发利用提供科学依据。

参考文献
[1]
马克岩, 韩金涛, 白雅琴, 等. 基于简化基因组测序的永登七山羊遗传多样性分析[J]. 畜牧兽医学报, 2023, 54(5): 1939-1950.
MA K Y, HAN J T, BAI Y Q, et al. Genetic diversity analysis of Yongdeng Qishan sheep based on specific-locus amplified fragment sequencing[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(5): 1939-1950. (in Chinese)
[2]
TONG X L, HAN M J, LU K P, et al. High-resolution silkworm pan-genome provides genetic insights into artificial selection and ecological adaptation[J]. Nat Commun, 2022, 13(1): 5619. DOI:10.1038/s41467-022-33366-x
[3]
GONG Y, LI Y F, LIU X X, et al. A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?[J]. J Anim Sci Biotechnol, 2023, 14(1): 73. DOI:10.1186/s40104-023-00860-1
[4]
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. DOI:10.1038/s41598-018-28719-w
[5]
LI R N, ZHAO Y H T, LIANG B M, et al. Genome-wide signal selection analysis revealing genes potentially related to sheep-milk-production traits[J]. Animals (Basel), 2023, 13(10): 1654.
[6]
SHI H B, LI T T, SU M C, et al. Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep[J]. BMC Genomics, 2023, 24(1): 50. DOI:10.1186/s12864-023-09146-2
[7]
薛周舣源, 宋显威, 吴林慧, 等. 畜禽选择信号检测方法及其统计学问题[J]. 畜牧兽医学报, 2018, 49(6): 1099-1107.
XUE Z Y Y, SONG X W, WU L H, et al. The identification methods of selection signatures in livestock and its statistical problems[J]. Acta Veterinaria et Zootechnica Sinica, 2018, 49(6): 1099-1107. (in Chinese)
[8]
TALEBI R, GHAFFARI M R, ZEINALABEDINI M, et al. Genetic basis of muscle-related traits in sheep: a review[J]. Anim Genet, 2022, 53(6): 723-739. DOI:10.1111/age.13266
[9]
PAN Z Y, LI S D, LIU Q Y, et al. Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization[J]. Gigascience, 2018, 7(4): giy019.
[10]
LI X, YANG J, SHEN M, et al. Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits[J]. Nat Commun, 2020, 11(1): 2815. DOI:10.1038/s41467-020-16485-1
[11]
LI Y F, GONG Y, ZHANG Z K, et al. Whole-genome sequencing reveals selection signals among Chinese, Pakistani, and Nepalese goats[J]. J Genet Genomics, 2023, 50(5): 362-365. DOI:10.1016/j.jgg.2023.01.010
[12]
LI H, DURBIN R. Fast and accurate long-read alignment with burrows-wheeler transform[J]. Bioinformatics, 2010, 26(5): 589-595. DOI:10.1093/bioinformatics/btp698
[13]
MCKENNA A, HANNA M, BANKS E, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data[J]. Genome Res, 2010, 20(9): 1297-1303. DOI:10.1101/gr.107524.110
[14]
LI H, HANDSAKER B, WYSOKER A, et al. The sequence alignment/map format and SAMtools[J]. Bioinformatics, 2009, 25(16): 2078-2079. DOI:10.1093/bioinformatics/btp352
[15]
CINGOLANI P, PLATTS A, WANG L L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3[J]. Fly (Austin), 2012, 6(2): 80-92. DOI:10.4161/fly.19695
[16]
PRICE A L, PATTERSON N J, PLENGE R M, et al. Principal components analysis corrects for stratification in genome-wide association studies[J]. Nat Genet, 2006, 38(8): 904-909. DOI:10.1038/ng1847
[17]
WEBER D, STEWART B S, GARZA J C, et al. An empirical genetic assessment of the severity of the northern elephant seal population bottleneck[J]. Curr Biol, 2000, 10(20): 1287-1290. DOI:10.1016/S0960-9822(00)00759-4
[18]
IVEY C T, HABECKER N M, BERGMANN J P, et al. Weak reproductive isolation and extensive gene flow between Mimulus glaucescens and M. guttatus in northern California[J]. Evolution, 2023, 77(5): 1245-1261. DOI:10.1093/evolut/qpad044
[19]
邱小宇, 黄勇富, 赵永聚, 等. 中国6个地方绵羊品种的微卫星分析[J]. 中国畜牧兽医, 2016, 43(4): 1024-1031.
QIU X Y, HUANG Y F, ZHAO Y J, et al. Microsatellite analysis of six Chinese indigenous sheep breeds[J]. China Animal Husbandry & Veterinary Medicine, 2016, 43(4): 1024-1031. (in Chinese)
[20]
韩旭, 刘丑生, 刘刚, 等. 中国地方绵羊线粒体DNA D-loop区遗传多样性研究[J]. 中国畜牧杂志, 2015, 51(7): 10-15.
HAN X, LIU C S, LIU G, et al. Study on genetic diversity of D-loop region of mitochondrial DNA of domestic sheep in China[J]. Chinese Journal of Animal Science, 2015, 51(7): 10-15. (in Chinese)
[21]
SUAREZ-GONZALEZ A, LEXER C, CRONK Q C B. Adaptive introgression: a plant perspective[J]. Biol Lett, 2018, 14(3): 20170688. DOI:10.1098/rsbl.2017.0688
[22]
SATOKANGAS I, NOUHAUD P, SEIFERT B, et al. Semipermeable species boundaries create opportunities for gene flow and adaptive potential[J]. Mol Ecol, 2023, 32(15): 4329-4347. DOI:10.1111/mec.16992
[23]
JANSSON E, FAUST E, BEKKEVOLD D, et al. Global, regional, and cryptic population structure in a high gene-flow transatlantic fish[J]. PLoS One, 2023, 18(3): e0283351. DOI:10.1371/journal.pone.0283351
[24]
LI X L, YUAN L F, WANG W M, et al. Whole genome re-sequencing reveals artificial and natural selection for milk traits in east Friesian sheep[J]. Front Vet Sci, 2022, 9: 1034211. DOI:10.3389/fvets.2022.1034211
[25]
LIU L L, MENG J, MA H Y, et al. Candidate genes for litter size in Xinjiang sheep identified by specific locus amplified fragment (SLAF) sequencing[J]. Anim Biotechnol, 2023, 34(7): 3053-3062. DOI:10.1080/10495398.2022.2131561
[26]
CRISPIM A C, KELLY M J, GUIMARAES S E F, et al. Multi-trait GWAS and new candidate genes annotation for growth curve parameters in Brahman cattle[J]. PLoS One, 2015, 10(10): e0139906. DOI:10.1371/journal.pone.0139906
[27]
KIM E S, ELBELTAGY A R, ABOUL-NAGA A M, et al. Multiple genomic signatures of selection in goats and sheep indigenous to a hot arid environment[J]. Heredity (Edinb), 2016, 116(3): 255-264. DOI:10.1038/hdy.2015.94
[28]
LU Z K, DU L X, LIU R Z, et al. MiR-378 and BMP-smad can influence the proliferation of sheep myoblast[J]. Gene, 2018, 674: 143-150. DOI:10.1016/j.gene.2018.06.039
[29]
KALDS P, LUO Q, SUN K, et al. Trends towards revealing the genetic architecture of sheep tail patterning: promising genes and investigatory pathways[J]. Anim Genet, 2021, 52(6): 799-812. DOI:10.1111/age.13133
[30]
SERRANITO B, CAVALAZZI M, VIDAL P, et al. Local adaptations of Mediterranean sheep and goats through an integrative approach[J]. Sci Rep, 2021, 11(1): 21363. DOI:10.1038/s41598-021-00682-z
[31]
LU Z K, LIU J B, HAN J L, et al. Association between BMP2 functional polymorphisms and sheep tail type[J]. Animals (Basel), 2020, 10(4): 739.
[32]
ZHANG Z B, LIU Q Y, DI R, et al. Single nucleotide polymorphisms in BMP2 and BMP7 and the association with litter size in Small Tail Han sheep[J]. Anim Reprod Sci, 2019, 204: 183-192. DOI:10.1016/j.anireprosci.2019.04.001
[33]
ZHANG L, LIU J S, ZHAO F P, et al. Genome-wide association studies for growth and meat production traits in sheep[J]. PLoS One, 2013, 8(6): e66569. DOI:10.1371/journal.pone.0066569
[34]
ZHU M T, ZHANG H M, YANG H, et al. Polymorphisms and association of GRM1, GNAQ and HCRTR1 genes with seasonal reproduction and litter size in three sheep breeds[J]. Reprod Domest Anim, 2022, 57(5): 532-540. DOI:10.1111/rda.14091
[35]
GOSSNER A G, HOPKINS J. Transcriptome analysis of CNS immediately before and after the detection of PrPSc in SSBP/1 sheep scrapie[J]. Vet Microbiol, 2014, 173(3-4): 201-207. DOI:10.1016/j.vetmic.2014.07.026
[36]
聂玮, 孟科, 荣轩, 等. 绵羊GRM1基因多态性及其与肉质性状的相关性[J]. 浙江农业学报, 2023, 35(4): 799-808.
NIE W, MENG K, RONG X, et al. Analysis of GRM1 gene polymorphism and its correlation with meat quality traits in sheep[J]. Acta Agriculturae Zhejiangensis, 2023, 35(4): 799-808. (in Chinese)
[37]
WANG X Y, FANG C, HE H Y, et al. Identification of key genes in sheep fat tail evolution based on RNA-seq[J]. Gene, 2021, 781: 145492. DOI:10.1016/j.gene.2021.145492
[38]
康丹菊. 绵羊脂肪沉积关键基因的筛选及ALDH1As在脂肪细胞分化中的调控研究[D]. 杨凌: 西北农林科技大学, 2018.
KANG D J. Identification of key genes related to adipose deposition in sheep and regulatory study of ALDH1As in adipocyte differentiation[D]. Yangling: Northwest A&F University, 2018. (in Chinese)
[39]
LIU S Q, LIU Z Y, WANG P, et al. Estrogen-mediated oar-miR-485-5p targets PPP1R13B to regulate myoblast proliferation in sheep[J]. Int J Biol Macromol, 2023, 236: 123987. DOI:10.1016/j.ijbiomac.2023.123987
[40]
DUAN Z, JI X, ZHU Y X, et al. Effects of sheep bone collagen peptide on liver lipid deposition in ovariectomized rats[J]. J Nutr Sci Vitaminol (Tokyo), 2022, 68(4): 320-330. DOI:10.3177/jnsv.68.320
[41]
LI Q, MILLAR R P, CLARKE I J, et al. Evidence that neurokinin b controls basal gonadotropin-releasing hormone secretion but is not critical for estrogen-positive feedback in sheep[J]. Neuroendocrinology, 2015, 101(2): 161-174. DOI:10.1159/000377702
[42]
SUN H R, MENG K, WANG Y F, et al. LncRNAs regulate the cyclic growth and development of hair follicles in Dorper sheep[J]. Front Vet Sci, 2023, 10: 118629.
[43]
WU C L, LI J Y, XU X M, et al. Effect of the FA2H Gene on cashmere fineness of Jiangnan cashmere goats based on transcriptome sequencing[J]. BMC Genomics, 2022, 23(1): 527. DOI:10.1186/s12864-022-08763-7
[44]
HE D Q, CHEN L Y, LUO F, et al. Differentially phosphorylated proteins in the crimped and straight wool of Chinese Tan sheep[J]. J Proteomics, 2021, 235: 104115. DOI:10.1016/j.jprot.2021.104115
[45]
PLOWMAN J E, HARLAND D P, RICHENA M, et al. Wool fiber curvature is correlated with abundance of K38 and specific keratin-associated proteins[J]. Proteins, 2022, 90(4): 973-981. DOI:10.1002/prot.26289
[46]
MORGAN H J, BENKETAH A, OLIVERO C, et al. Hair follicle differentiation-specific keratin expression in human basal cell carcinoma[J]. Clin Exp Dermatol, 2020, 45(4): 417-425. DOI:10.1111/ced.14113
[47]
KOWALCZYK A, CHIKINA M, CLARK N. Complementary evolution of coding and noncoding sequence underlies mammalian hairlessness[J]. Elife, 2022, 11: e76911. DOI:10.7554/eLife.76911
[48]
HUI T Y, ZHENG Y Y, YUE C, et al. Screening of cashmere fineness-related genes and their ceRNA network construction in cashmere goats[J]. Sci Rep, 2021, 11(1): 21977.

(编辑   郭云雁)