14次 摘要：为了解二毛期时滩羊串子花型裘皮与其它类型裘皮中蛋白质的差异,本试验采用iTRAQ技术及LC-MS/MS蛋白质组学研究方法,对串子花型、软大花型、绿豆丝型及其它不规则型花穗裘皮蛋白质进行鉴定和筛选,并运用Proteome Discoverer l.4软件进行定量分析,结合数 --> 　　摘要：为了解二毛期时滩羊串子花型裘皮与其它类型裘皮中蛋白质的差异,本试验采用iTRAQ技术及LC-MS/MS蛋白质组学研究方法,对串子花型、软大花型、绿豆丝型及其它不规则型花穗裘皮蛋白质进行鉴定和筛选,并运用Proteome Discoverer l.4软件进行定量分析,结合数据库搜索,鉴定出具有显著表达差异的蛋白,同时应用生物学技术对其进行GO和Pathway分析。结果显示4类花穗型裘皮共检测出2 886个蛋白,其中有135个、142个、113个差异蛋白分别存在于软大花型与串子花型、绿豆丝型与串子花型、其它不规则型与串子花型3个对比组中。对有表达差异的蛋白进行分析,发现膜联蛋白与血管内皮生长因子可能与软大花型毛股形成相关,KAP3和KAP6可能与滩羊串子花型毛股结构相关。研究发现滩羊不同二毛裘皮蛋白水平上的差异,可为选育优良的滩羊串子花型二毛裘皮提供理论基础。　　关键词：滩羊; iTRAQ; 裘皮; 差异蛋白;　　滩羊是一种珍贵而著名的产裘皮绵羊品种，是在特殊的地理和气候条件环境下形成其特有的以轻、暖、花穗美丽而著称的二毛裘皮性状(赵金宇,2013)。花穗一词概括出了滩羊二毛皮的综合性状，其概念主要包括花形和毛型比例两方面内容。其美观程度由花形决定，保暖性能则取决于毛型。故以其毛股花形为区分标准，滩羊二毛皮花穗可分为平波形和不规则形两类。　　平波形的特点是其毛股上部具有弧度均匀的波形弯曲，且各个波纹在同一平面上排列、延伸，其弯曲规律、整齐可形成花案；凡不属于平波状的均归为不规则花形，将其统称为“不规则花穗”。以毛型比例及毛股的粗细为区分依据，平波形花穗可分为3类，其中毛股中无髓毛含量(50%)、毛股粗细适中，下部结构无显著增大或稍有增大，且绒毛长度适中的称为串子花型；而以毛股中含绒毛60%以上且毛股下部显著增大为特点的称为软大花型；第3种绿豆丝型花穗毛股中含30%～40%的无髓毛，较前两种少，绒毛较短，毛股较细，结构紧密(张幼麟等,1965,宁夏农林科技,12:33-34)。羊皮肤是滩羊毛股花穗的重要载体，覆盖于体表并直接与外界接触，起重要的保护作用。　　国内外目前对绒山羊、细毛羊和超细毛羊等羊皮肤蛋白组学方面的研究较多，而对于滩羊皮肤蛋白组学方面的相关报道却是屈指可数。众所周知，滩羊毛以两型毛为主，但其形成机理目前尚不甚明了，故本试验以滩羊二毛期皮肤为研究对象，利用i TRAQ研究其不同花穗型皮肤的差异表达蛋白，以期从蛋白质层面找出滩羊二毛期出现不同花穗型裘皮的原因，为串子花型滩羊二毛裘皮的选育提供科学依据。　　1结果与分析　　1.1 SDS-PAGE电泳　　4种花穗型裘皮的8份样品经蛋白质定量(BCA法,上样量为20μg)后，用SDS-PAGE电泳进行分离，最终获得的蛋白质条带(图1)，其平行度较好，条带较均一。　　1.2蛋白鉴定及其整体分布　　4种花穗裘皮样品经i TRAQ技术中的强阳离子交换柱分离，获得8组肽，经MS仪器分析，结果经Mascot查库，之后以Peptide FDR≤0.01为参照筛选过滤。结果有12 988个肽段(Peptide)被鉴定出，得到共计2 886个蛋白质组(Protein group)。以其相对分子质量、肽段数量及肽段序列覆盖率为依据对鉴定出的所有蛋内质进行统计，结果显示利用i TRAQ技术鉴定出的不同花穗裘皮样品中的蛋白，其中分子量在20～80 k D(A)之间的蛋白质数量较多，包含5个以内(B)的肽段的蛋白数量居多，蛋白覆盖率则大多集中于0%～20%之间(C)(图2)。　　图1 样品SDS-PAGE电泳　　　　Figure 1 Sample SDS-PAGE electrophoresis　　注:A1,A2:串子花型;B1,B2:软大花型;C1,C2:绿豆丝型;D1,D2:其它不规则型　　1.3各组显著性差异蛋白定量分析　　分别统计软大花型(B)与串子花型(A)、绿豆丝型(C)与串子花型(A)、其它不规则型(D)与串子花型(A)3个比对组得到的差异蛋白数量，对其进行分析，并定义其中ratio>1.3且p<0.05的蛋白为上调的蛋白，符合ratio<0.77且p<0.05的蛋白为下调的蛋白，分析各比较组的差异蛋白具体统计结果(表1)。　　依照各比较组筛选出的差异蛋白的数目，除未命名蛋白外，列出各比较组已鉴定出的存在显著差异的蛋白，得到其统计结果(表2)。　　统计分析3组经对比得到的差异蛋白，发现共有56个差异蛋白共同存在于B VS A与C VS A组中，其中上调蛋白21个，下调蛋白35个；而有46个相同的差异蛋白同时存在于B VS A与D VS A两组中，其中上调蛋白数19个，下调蛋白数27个；有59个相同的差异蛋白存在于C VS A与D VS A两组中，属于上调蛋白的22个，属于下调蛋白的37个；并且有36个相同的差异蛋白同时存在于B VS A、C VS A、D VS A组中，其中上调蛋白数为13个，下调蛋白数为23个，得到统计结果(图3)。所列出的差异蛋白主要是促甲状腺激素释放激素降解酶、细胞周期蛋白依赖性激酶5、非特异性丝氨酸/苏氨酸蛋白激酶等蛋白酶类(表2)。　　1.4 GO分析　　3个比较组得到的差异蛋经过GO(Gene ontology)分析，结果在B VS A组的差异蛋白中发现其所涉及到的分子功能共有11种，生物学过程与细胞组成分别为22种和16种(图4A)；在C VS A组中发现的差异蛋白共涉及到分子功能12种，生物学过程22种，细胞组成15种(图4B);D VS A组中涉及到的分子功能共7种，生物学过程共22种，细胞组成共11种(图4C)。本试验检测出的3个比较组的差异蛋白主要涉及到催化、转运蛋白、抗氧化、分子功能调节剂等众多重要分子功能；涉及到生物的生长、生殖、代谢及免疫系统等一系列重要生物学过程；涉及到大分子复合物、细胞连接、细胞器、细胞外基质等细胞组成。　　图2 i TRAQ鉴定蛋白的总概况　　　　Figure 2 i TRAQ to identify the total protein profiles　　注:A:鉴定蛋白的质量分布;B:蛋白所含肽段的数量分布;C:肽段的序列覆盖度　　图3 各比对组差异蛋白分析　　　　Figure 3 Analysis of each match differentially expressed proteins　　1.5 KEGG分析　　对三比较组得到的差异蛋白进行KEGG生物通路富集分析，按p值分析各组的KEGG通路，由左至右p值增大，以p<0.05与p<0.01对其进行分类(图5),B VS A组中的差异蛋白涉及碳代谢、VEGF(Vascular endothelial growth factor)信号通路、丙酮酸代谢、花生四烯酸代谢、NOD样受体信号通路、FcεRI信号通路等共计14条通路(图5A)。C VS A组中比对得到的差异蛋白涉及丁酸代谢、维生素B6代谢、各种氨基酸代谢与降解及脂肪酸的代谢与降解等共计13条具体代谢通路(图5B)。在D VS A比对组中得到的差异蛋白涉及丙酸酯代谢、氨基酸代谢及降解、脂肪酸代谢及降解、芳香化合物的降解等共计21条代谢通路(图5C)。　　表1 差异蛋白统计　　　　表2 各比对组差异表达蛋白　　　　　　　　2讨论　　2.1 KAP3　　羊毛纤维几乎全部由蛋白质组成，主要组成成分是角蛋白与角蛋白相关蛋白(keratin-associated proteins,KAPs)，而KAPs则通过与角蛋白的相互作用影响毛纤维的机械性能(Plowman et al.,2015)。以毛纤维中所含半胱氨酸或甘氨酸与酪氨酸残基的比例为依据，人为地将KAPs又分为高硫蛋白(high-sulphur proteins,HSPs)、超高硫蛋白(ultrahigh-sulphur proteins,UHSPs)和高甘氨酸-酪氨酸(high-glycine-tyrosine proteins,HGTPs)蛋白(Rogers,2006)。有研究发现不同KAPs基因与羊毛绒品质具有一定的关联性(高建军等,2014)，这一观点与之前研究发现的细毛羊毛纤维中KAPs中的HSPs基因与毛纤维细度有着极大的关联(刘桂芬等,2007)相吻合。KAP3属于HSPs关联蛋白家族的一种，其中还包括KAP1和KAP2两种蛋白(Plowman et al.,2009)。已有相关研究发现KAP3.2基因与藏系绵羊毛长和产毛量性状存在相关性，并可作为这两个性状的分子标记(王志有等,2011)。有研究者利用RT-PCR技术研究藏绵羊KAP3.2基因的结构与功能，结果发现该基因在皮肤组织的表达存在特异性，对羊毛性状具有重要影响(杨珂伟等,2014)，而且KAP3.2基因在新吉毛羊超细型中的表达量极显著地高于其在新吉细型细毛羊中的表达量，可能影响毛纤维机制表达的代谢(刘春洁等,2015)。研究发现KAP3.3在羊毛纤维皮质层表达(Plowman et al.,2009)，进一步研究发现KAP3.3基因以组织特异性方式在羊皮肤毛囊中高度表达(Zhao et al.,2016)。本试验中KAP3是在C VS A对比组中发现的差异蛋白，且为下调蛋白，较绿豆丝型裘皮而言，串子花型裘皮中的KAP3高度表达，故而KAP3可能在滩羊串子花型毛股结构的形成过程中起重要作用。　　2.2 KAP6　　H GTPs家族包括KAP6、KAP7和KAP8等蛋白，在毛纤维的主要结构及机械性能方面扮演着重要的角色(Liu et al.,2014)。KAP6在羊毛纤维中的表达量会对毛纤维直径产生一定的影响(Cockett et al.,2001)。已有研究显示绵羊羊毛纤维直径的变化受到了KAP6.1基因变化的影响，并且在KAP6.1基因的57 bp片段缺失后，粗羊毛纤维直径的相关特性会发生极大变化(Zhou et al.,2015)。有研究者用PCR-SSCP方法对绒山羊KAP6.3基因编码区和3'-非翻译区进行多态性扫描，并分析其变异与绒山羊绒毛性状的关系，结果发现绒毛的长度、细度、产绒量均受KAP6.3不同SNPs位点的影响(刘海英,2007)。之后有人用同样的方法研究了高原型藏山羊产绒性状与KAP6.1和KAP6.2基因位点间的关系，结果显示KAP6.1主要影响产绒量性状，而KAP6.2基因与藏山羊的绒毛长度、细度等存在更大地相关性(王杰等,2010)，并于2016年，有研究者发现了羊KAP6.4和KAP6.5这两个新的基因，进一步完善了KAP6基因家族(Zhou et al.,2016)。在本试验中KAP6同样是在C VS A对比组中发现的差异蛋白，其ratio<0.77，故较KAP6在绿豆丝型裘皮中的表达而言，KAP6在串子花型裘皮中属于高度表达，从KAP6在滩羊不同花穗型裘皮的表达量可推测，KAP6在滩羊串子花型裘皮中高度表达与其串子花型毛股形成有着紧密的联系。　　图4 差异蛋白GO注释　　　　Figure 4 Differential protein GO annotation　　注:A:软大花型/串子花型;B:绿豆丝型/串子花型;C:其它不规则型/串子花型　　图5 差异蛋白KEGG富集分析　　　　Figure 5 Differential protein KEGG enrichment analysis　　注:A:软大花型/串子花型;B:绿豆丝型/串子花型;C:其它不规则型/串子花型　　2.3膜联蛋白　　膜联蛋白是一类磷脂结合蛋白，它的表达受到钙离子的调控。作为膜联蛋白多基因家族成员之一膜联蛋白A2则广泛表达于多细胞生物体内，并在多种类型的细胞中均有存在，起初发现于Rous肉瘤病毒转化的鸡胚胎成纤维细胞中(Hastie et al.,2008)。研究表明，膜联蛋白A2与多种炎症机制的反应有关，而且参与细胞粘附及胞吐等过程，并且在宫颈组织中的过度表达与宫颈癌的发生发展关系密切(赵丽萍和陆晓媛,2010)。有研究表明，膜联蛋白A2与人类多种心血管、肿瘤及病毒感染疾病有关，当膜联蛋白A2的表达量较高时，纤溶酶的表达量增多，加快细胞外基质的降解，对血管再生有一定的积极作用(聂纪芹等,2012;Yi et al.,2016)，而且李璟波等(2015)还发现膜联蛋白A2在非小细胞肺癌的诊断中具有一定的意义。本试验中，膜联蛋白相对而言在软大花型中的表达量较高，为上调蛋白，通过生物信息分析，发现其主要涉及到代谢过程的调节、营养水平反应的调节、对细胞外环境刺激的应答反应等生物学过程，并且与细胞的多种组成成分密不可分，涉及到阴阳离子结合、钙依赖性磷脂结合等系列分子功能。因此可推测其与软大花型毛股的形成具有一定的关联。　　2.4 VEGF信号通路　　血管内皮生长因子(VEGF)是主要作用于血管内皮细胞的血管生成因子，与血管正常的生长、发育以及伤口愈合有着密切的关系，与血管肿瘤生长和转移等过程有关，是一种促进血管生成活性的功能性蛋白，目前发现的VEGF主要有胎盘生长因子与VEGF-A、VEGF-B、VEGF-C、VEGF-D等5个成员(Nash et al.,2006;Saito et al.,2007)。VEGF信号通路具有启动并调节胚胎血管的发育的作用(Aguago et al.,2002;姜怀志等,2010)。VEGF特异性作用于内皮细胞调控着血管新生，影响着血管的通透性，并与其它因子协同作用，促进血管内皮细胞增殖、迁移，诱导血管生长、形成及其相关蛋白酶的合成(Jalali et al.,2001;Wagatsuma,2007;姜怀志,2012)。VEGF对皮肤毛囊血管新生具有重要意义，调控着毛囊的周期性变化(姜怀志,2012)。已有研究发现VEGF可以促进毛囊发育,并影响毛囊直径(Araújo et al.,2011)。无独有偶，有研究者测定了辽宁绒山羊常年长绒和季节长绒两个品系皮肤毛囊中VEGF的相对表达量，发现辽宁绒山羊常年长绒型品系常年长绒可能与VEGF基因m R-NA的相对表达量常年维持在一个较高的水平有关(姚纪元等,2010)。之后又有人发现在胎儿期，绒山羊体内的VEGF与皮肤上微血管密度呈显著强的正相关(谷博等,2013)，这说明VEGF可增加毛囊的微血管密度，以利于毛囊发育。并且早期有研究已表明VEGF能维持毛囊超微结构的完整性，促进毛囊生长(Bruno et al.,2009)。本试验中对B VS A组的差异蛋白进行KEGG通路分析发现其涉及VEGF信号通路，在软大花型裘皮中涉及VEGF信号通路的差异蛋白有较高表达量，属于上调蛋白，故而可推测VEGF信号通路所涉及的差异蛋白在滩羊二毛期花穗形成过程中的高表达有利于产生更多的绒毛，形成软大花型毛股。　　3材料与方法　　3.1样品采集　　4种滩羊二毛期的不同花穗型裘皮样品均采集于宁夏盐池的滩羊选育场(图6)，取样部位为羊体侧肩胛骨后缘1 cm2处，每种花穗型裘皮的滩羊采2份样品，共8份，采集的样品经生理盐水处理，立即置于装有干冰的泡沫盒中，带至实验室，保存于超低温冰箱(-80℃)备用，得到样品信息(表3)。　　3.2实验仪器与试剂　　E PS601 600V电泳仪、AKTA Purifier 100纯化仪均由GE Healthcare提供；Q-Exactive质谱仪、纳升级液相色谱仪(型号:Easy n LC1000)均购于Thermo Finnigan；低温高速离心机，型号为Eppendorf 5430R;Eppendorf Concentrator plus真空离心浓缩仪；UV-2100型可见紫外分光光度计(尤尼柯WFZ)。　　图6 滩羊二毛期不同毛股花穗型　　　　Figure 6 Tan sheep er-mao different hairs spike type　　注:A:串子花型;B:软大花型;C:绿豆丝型;D:其它形态　　表3 样品信息　　　　SDS(161-0302,Bio-Rad),Urea(161-0731,Bio-rad),Trisbase(161-0719,Bio-rad),DTT(161-0404,Bio-rad),IAA(Bio-rad,163-2109),KH2PO4(10017618,国药)，NH4HCO3(Sigma,A6141);i TRAQ Reagent-8plex Multiplex Kit(AB,SCIEX),Acetonitrile,ACN(I59223-0123,Merck);C18Cartridge(66872-U,Sigma);5×上样缓冲液；SDT Buffer;UA Buffer;Dissolution Buffer(AB,SCIEX)。　　3.3样品的制备　　将不同类型的皮肤样品分别置于研钵中，加入适量液氮研磨至粉状，加入STD Buffer(150 mmol/L Tris-HCl,4%SDS,1 mmol/L DTT,p H=8.0)溶液600μL，匀浆，然后沸水浴5 min，用超声波破碎(80 W,每次超声8 s,间隔10 s,持续5 min)，置于沸水浴5 min,12 000 g离心持续15 min，提取上清液，进行BCA蛋白质定量。　　3.4 SDS-P AGE电泳　　在上述BCA法定量后的8份样品中各取20μg，加入5×上样缓冲液(体积比5∶1)，经5 min沸水水浴，10 min 14 000 g的高速离心后，提取上清液，用于进行12.5%的SDS-PAGE电泳(恒流:14 m A,持续90 min)，用考马斯亮蓝染色。　　3.5胰蛋白酶酶解　　各取200μg的裘皮样品，用胰蛋白酶酶解，具体步骤详见杨佐青等(2017)，得到酶解肽。　　3.6 i TRAQ标记　　参照OD280肽段定量的结果，分别取约80μg的4种不同毛股花穗型裘皮样品，共8份，依据i TRAQ试剂盒说明书进行标记，将标记后的各组肽段混合，得到样品标记信息(表4)。　　3.7毛细管高效液相色谱分析　　用Easy n LC液相系统分离上述标记的样品，具体方法参照“应用i TRAQ技术对滩羊裘皮差异表达蛋白的研究”一文(杨佐青等,2017)。　　表4 样品标记　　　　3.8质谱鉴定　　用Q-Exactive质谱仪对上述经色谱分离后的样品进行MS分析，具体分析模式参照杨佐青等(2017)。　　3.9质谱数据分析　　用Mascot2.2和Proteome Discoverer1.4(Thermo)软件对经质谱鉴定产生的原始数据进行查库(uniprot Ovis aries蛋白质库)鉴定和定量分析，以FDR≤0.01为蛋白鉴定的过滤参数，肽段报告离子峰强度值选用Proteome Discoverer软件进行定量分析。　　3.10差异蛋白分析　　以A1和A2两份样品的113和114两组通道平均值为参数，其余各通道标签与参数的比值即为各组蛋白质的i TRAQ比值。将鉴定出的蛋白质进行T检验与显著性分析(p value)确定出存在显著差异的蛋白。对得到的差异蛋白进行分析，当某一差异蛋白的丰度比在1.3以上或0.77以下，且其p<0.05时，则将该蛋白定义为不同样品间的显著差异蛋白。　　3.11生物信息学分析　　根据相似性原理，即序列相似的蛋白其功能也可能相似，将各比较组得到的差异蛋白与GO数据库以及KEGG pathway数据库进行对比，并对其进行GO注释和KEGG通路分析，最终将具有同源性的蛋白的注释信息分别转嫁到各自所对应的差异蛋白上，以研究不同花穗型皮肤样中存在表达差异的蛋白所起的生物学作用及其可能参与的机体代谢过程。　　参考文献　　[]Aguago A.,Kantarjian H.M.,Estey E.H.,Giles F.J.,Verstovsek S.,Manshouri T.,Gidel C.,O'Brien S.,Keating M.J.,and Albitar M.,2002,Plasma vascular endothelial growth factor levels have prognostic significance in patients with acute myeloid leukemia but not in patients with myelodysplastic syndromes,Cancer,95(9):1923-1930　　[]Araújo V.R.,Silva G.M.,Duarte A.B.,Magalhúes D.M.,Almeida A.P.,Gon觔alves R.F.,Bruno J.B.,Silva T.F.,Campello C.C.,Rodrigues A.P.,and Figueiredo J.R.,2011,Vascular endothelial growth factor-A165(VEGF-A165)stimulates the in vitro development and oocyte competence of goat preantral follicles,Cell Tissue Res.,346(2):273-281　　[]Bruno J.B.,Celestino J.J.,Lima-Verde I.B.,Lima L.F.,Matos M.H.T.,Araújo V.R.,Saraiva M.V.A.,Martins F.S.,Name K.P.O.,Campello C.C.,Báo S.N.,Silva J.R.V.,and Figueiredo J.R.,2009,Expression of vascular endothelial growth factor(VEGF)receptor in goat ovaries and improvement of in vitro caprine preantral follicle survival and growth with VEGF,Reproduction,Fertility and Development,21(5):679-687　　[]Cockett N.E.,Shay T.L.,and Smit M.,2001,Analysis of the sheep genome,Physiological Genomics,7(2):69-78　　[]Gao J.J.,Xie T.T.,Li S.W.,and Wang H.J.,2014,Prokaryotic expression and bioinformation analysis of wool follicle KAP7 of Hetian sheep,Zhongguo Xumu Zazhi(Chinese Journal of Animal Science),50(21):63-68(高建军,谢婷婷,李树伟,王海军,2014,和田羊毛囊KAP7基因的原核表达及其生物信息学分析,中国畜牧杂志,50(21):63-68)　　[]Gu B.,Sun L.M.,Chang Q.,and Jiang H.Z.,2013,Study on VEGF in the expression of skin hair follicle and the relationship with microvessel density by liaoning cashmere goat in fetal period,Zhongguo Xumu Shouyi(China Animal Husbandry&Veterinary Medicine),40(6):158-161(谷博,孙丽敏,常青,姜怀志,2013,血管内皮生长因子基因在辽宁绒山羊胎儿期皮肤毛囊发育中表达及其与微血管密度关系的研究,中国畜牧兽医,40(6):158-161)　　[]Hastie C.,Masters J.R.,Moss S.E.,and Naaby-Hansen S.,2008,Interferon reduces cell surface expression of Annexin2 and suppresses the invasive capacity of prostate cancer cells,Biol.Chem.,283(18):12595-12603　　[]Jalali S.,del Pozo M.A.,Chen K.,Miao H.,Li Y.S.,Schwartz M.A.,Shyy J.Y.J.,and Chien S.,2001,Integrin-mediated mechanotransduction requires its dynamic interaction with specific extracellular matrix(ECM)ligands,Proc.Nat.Acad.Sci.,98(3):1042-1046　　[]Jiang H.Z.,Chen Y.,and Chang Q.,2010,Effect of vascular endothelial growth factor(VEGF)on regulation of angiogenesis around skin follicle of mammal,Zhongguo Xumu Shouyi(China Animal Husbandry&Veterinary Medicine),37(5):47-49(姜怀志,陈洋,常青,2010,血管内皮生长因子在哺乳动物皮肤毛囊周围血管新生过程中的调控作用,中国畜牧兽医,37(5):47-49)　　[]Jiang H.Z.,2012,Hair follicle structure traits and developmental mechanisms of Chinese cashmere goat,Jilin Nongye Daxue Xuebao(Journal of Jilin Agricultural University),34(5):473-482(姜怀志,2012,中国绒山羊的毛囊结构特性与发育机制,吉林农业大学学报,34(5):473-482)　　[]Li J.B.,Yu C.H.,Guo N.N.,Tang J.,Liu Y.,and Zhang J.H.,2015,Research of enzyme-linked protein A2 and cathepsin B expression in non-small cell lung cancer,Linchuang Feike Zazhi(Journal of Clinical Pulmonary Medicine),20(5):815-817(李璟波,于长海,郭楠楠,唐健,刘颖,张建华,2015,膜联蛋白A2及组织蛋白酶B在肺鳞癌和腺癌中表达的研究,临床肺科杂志,20(5):815-817)　　[]Liu C.J.,Fu X.F.,Tian K.C.,Huang X.X.,Xu X.M.,Zhan Z.H.,and Yu L.J.,2015,Screening and expression of keratin ESTin Xinji fine wool sheep,Zhongguo Xumu Shouyi(China Animal Husbandry&Veterinary Medicine),42(11):2873-2879(刘春洁,付雪峰,田可川,黄锡霞,徐新明,詹振宏,于丽娟,2015,新吉细毛羊角蛋白的EST筛选及表达研究,中国畜牧兽医,42(11):2873-2879)　　[]Liu G.F.,Tian K.C.,Zhang E.P.,Huang X.X.,and Zhang Y.H.,2007,Candidate gene analysis of high quality merino sheep,Yichuan(Genetics),29(1):70-74(刘桂芬,田可川,张恩平,黄锡霞,张艳华,2007,优质细毛羊羊毛细度的候选基因分析,遗传,29(1):70-74)　　[]Liu H.Y.,2007,Influence of nucleotide variations in KAP8,KA-P6.3 and FGF5 genes on fibre traits and change law of fibre quality,Dissertation for Ph.D.,China Agricultural University,Supervisor:Jia Z.H.,pp.61-79(刘海英,2007,KAP8,KAP6.3和FGF5基因变异对绒毛性状的影响及绒毛品质变化规律的研究,博士学位论文,中国农业大学,导师:贾志海,pp.61-79)　　[]Liu Y.X.,Shi G.Q.,Wang H.X.,Wan P.C.,Tang H.,Yang H.,and Guan F.,2014,Polymorphisms of KAP6,KAP7,and KAP8 genes in four Chinese sheep breeds,Gene.Mol.Res.,13(2):3438-3445　　[]Nash A.D.,Baca M.,Wright C.,and Scotney P.D.,2006,The biology of vascular endothelial growth factor-B(VEGF-B),Pulm.Pharmacol.Ther.,19(1):61-69　　[]Nie J.Q.,Wang G.Z.,and Liu C.Q.,2012,Prokaryotic expression and purification of the human annexinA2 protein and preparation of its polyclonal antibody,Zhongguo Mianyixue Zazhi(Chinese Journal of Immunology),28(1):67-70(聂纪芹,汪龚泽,刘朝奇,2012,annexinA2基因片段的克隆,表达,纯化及多克隆抗体的制备,中国免疫学杂志,28(1):67-70)　　[]Plowman J.E.,Debchoudhury S.,Bryson W.G.,Clerens S.,and Dyer J.M.,2009,Protein expression in orthocortical and paracortical cells of merino wool fibers,J.Agric.Food Chem.,57(6):2174-2180　　[]Plowman J.E.,Harland D.P.,Ganeshan S.,Woods J.L.,Van Shaijik B.,Deb-Choudhury S.,Thomas A.,Clerens S.,and Scobie D.R.,2015,The proteomics of wool fibre morphogenesis,J.Structural Biol.,191(3):341-351　　[]Rogers G.E.,2006,Biology of the wool follicle:an excursion into a unique tissue interaction system waiting to be rediscovered,Exp.Dermatol.,15(12):931-949　　[]Saito A.,Sugawara A.,Uruno A.,Kudo M.,Kagechika H.,Sato Y.,Owada Y.,Kondo H.,Sato M.,Kurabayashi M.,Lmaizumi M.,Tsuchiya S.,and Ito S.,2007,All-trans retinoic acid induces in vitro angiogenesis via retinoic acid receptor:Possible involvement of paracrine effects of endogenous vascular endothelial growth factor signaling,Endocrinol.,148(3):1412-1423　　[]Wagatsuma A.,2007,Endogenous expression of angiogenesis-related factors in response to muscle injury,Mol.Cell Biochem.,298(1):151-159　　[]Wang J.,Dai Y.J.,Wang Y.,Zi X.D.,Ouyang X.,and Liu L.S.,2010,KAP 6.1 and KAP 6.2 loci and its relation to down traits of tibetan goats,Xinan Minzu Daxue Xuebao(Journal of Southwest University for Nationalities(Natural Science Edition)),36(6):962-966(王杰,代怡倩,王永,字向东,欧阳熙,刘鲁蜀,2010,高原型藏山羊KAP6.1和KAP6.2位点与产绒性状的关系研究,西南民族大学学报(自然科学版),36(6):962-966)　　[]Wang Z.Y.,Chen Y.L.,Xu Q.L.,Ma Z.R.,and Qi Q.Q.,2011,Polymorphism of KAP3.2 gene and its effect on partial economic traits in Tibetan sheep,Xumu Shouyi Xuebao(Chinese Journal of Animal and Veterinary Sciences),42(2):284-288(王志有,陈玉林,徐秋良,马忠仁,祁全青,2011,藏系绵羊KAP3.2基因多态性及其对部分经济性状的影响,畜牧兽医学报,42(2):284-288)　　[]Yang K.W.,Hu L.,Luo B.,Liu S.,Zheng Y.C.,and Zi X.D.,2014,Cloning and tissue expression of KAP3.2 gene in tibetan sheep,Xinan Minzu Daxue Xuebao(Journal of Southwest University for Nationalities(Natural Science Edition)),40(6):809-813(杨珂伟,胡亮,罗斌,刘霜,郑玉才,字向东,2014,藏绵羊KAP3.2基因的cDNA克隆,序列分析及组织表达的研究,西南民族大学学报(自然科学版),40(6):809-813)　　[]Yang Z.Q.,Tao J.Z.,Li Y.K.,Ding W.,Liu Y.B.,and Wei B.D.,2017,Study on differentially expressed proteins of Tan sheep skin using i TRAQ technology,Zhongguo Xumu Shouyi(China Animal Husbandry&Veterinary Medicine),44(9):2682-2691(杨佐青,陶金忠,李颖康,丁伟,刘永斌,魏炳栋,2017,应用i TRAQ技术对滩羊裘皮差异表达蛋白的研究,中国畜牧兽医,44(9):2682-2691)　　[]Yao J.Y.,Bao H.X.,Luan W.M.,and Jiang H.Z.,2010,Study on expression of VEGF gene in hair follicles of Liaoning cashmere goats,Zhongguo Xumu Shouyi(China Animal Husbandry&Veterinary Medicine),37(12):139-141(姚纪元,包红喜,栾维民,姜怀志,2010,辽宁绒山羊皮肤毛囊血管内皮生长因子基因的表达研究,中国畜牧兽医,37(12):139-141)　　[]Yi J.,Zhu Y.,Jia Y.,Jiang H.D.,Zheng X.,Liu D.J.,Gao S.X.,Sun M.J.,Hu B.,Jiao B.H.,Wang L.H.,and Wang K.H.,2016,The annexin a2 promotes development in arthritis through neovascularization by amplification hedgehog pathway,PLo S One,11(3):e0150363　　[]Zhao J.Y.,2013,Study on the growth-development and productivity of Tan sheep in Ningxia,Thesis for M.S.,Northwest A&F University,Supervisor:Luo J.,pp.1-6(赵金宇,2013,宁夏滩羊早期生长发育规律和生产性能研究,硕士学位论文,西北农林科技大学,导师:罗军,pp.1-6)　　[]Zhao Z.C.,Liu G.B.,Li X.Y.,Huang J.,Xiao Y.J.,Du X.Y.,and Yu M.,2016,Characterization of the promoter regions of two sheep keratin associated protein genes for hair cortexspecific expression,PLoS One,11(4):e0153936　　[]Zhao L.P.,and Lu X.Y.,2010,Expression of annexin A2 and Bcl-2in cervical cancer and its significance,Xuzhou Yixueyuan Xuebao(Acta Academiae Medicinae Xuzhou),30(8):507-510(赵丽萍,陆晓媛,2010,Annexin A2,Bcl-2在宫颈癌中的表达及意义,徐州医学院学报,30(8):507-510)　　[]Zhou H.,Gong H.,Li S.,Luo Y.,and Hickford J.G.H.,2015,A57 bp deletion in the ovine KAP6-1 gene affects wool fibre diameter,Journal of Animal Breeding and Genetics,13(4):301-307　　[]Zhou H.T.,Gong H.,Wang J.Q.,Dyer J.M.,Luo Y.Z.,and Hickford J.G.H.,2016,Identification of four new gene members of the KAP6 gene family in sheep,Sci.Rep.,6:24074