本发明属于羊品种鉴定技术领域，尤其涉及一种豫西脂尾羊snp标记及其筛选方法和应用。
背景技术：
豫西脂尾羊(yuxifat-tailedsheep)是河南省一个地方优良品种，1985年被编入《河南省畜禽品种志》一书。多年来，由于环境影响小，豫西脂尾羊品质的变化不明显，其优良性状耐粗饲、抗病力强、耐炎热和长膘快的特点得以保存。它适应性强，生长发育快，性成熟早，羔羊成活率高，肉质好，屠宰率高达57％。其肉质鲜美、细嫩，营养丰富，风味独到，色正味佳，温而不热，油而不腻，用豫西脂尾羊肉熬制的羊肉汤呈乳白色，香、鲜、浓，香而不膻，爽而不黏，是豫西特色名吃。市场开发前景广阔。
了解本土肉用绵羊品种的遗传多样性及其遗传关系，对于品种起源进化研究、种质资源保护以及科学开发利用均具有重要的指导意义。随着杂交育种的推广，一方面，高代杂交群体的个体表型与亲本的高度相似，无法单纯地利用豫西脂尾羊的毛色、角型、体型等表型性状特征进行准确区分；另一方面，生产制作出的动物产品已经不能再表现出所属品种的表型性状特征，因此也就无法利用表型性状对其进行鉴别。因此，仅利用表型性状特征进行品种资源的鉴定不够准确、全面和科学。传统的标记鉴定方法诸如细胞学标记、生物化学标记以及免疫学标记，虽然已经在品种鉴定中进行应用，但由于这些标记的多态性低和信息量小，且是对基因的间接反映，受到基因与环境互作的影响比较大，在生产实践中有很大的局限性。
dna分子遗传标记问世为克服这一难题提供了有效的解决途径。dna分子遗传标记具有不受表型、时空的限制、多态性好等特点，因此常常被用来进行品种鉴定的工作。其中，最常用的dna分子标记是微卫星标记和snp(singlenucleotidepolymorphism)标记。
微卫星标记被首先用于进行个体鉴定，其缺点是由于不同实验室的条件不同而难以实现标准化和规模化。
单核苷酸多态性(singlenucleotidepolymorphism，snp)是指在基因组dna水平上发生的单个核苷酸变异所引起的dna序列多态性。随着山羊与绵羊的全基因组测序工作的完成，羊基因组的信息基本趋于完善。目前已知羊的基因组大约有26亿对碱基，其中可能发生变异的位点多达46401347(ncbidbsnp142)，平均每60个碱基就可能出现一个snp。研究表明健康个体之间基因组碱基差异大约为0.1％，即两个随机个体中每1200-1490个碱基就会有一个差异碱基。实验室现有羊基因组信息的个体数量约为600只，snp位点信息约有40-50m。这些snps不仅是造成品种间差异的重要因素，而且大量研究还证实，有些snp与羊的多种经济性状都有关系。随着snp分型技术的发展，snp标记的发现和定位越来越多。相对于以限制性片段长度多态性为代表的第一代遗传标记和微卫星多态性为代表的第二代遗传标记，snp具有分布广泛、数量多等特点，成为了第三代遗传标记，更加适合于基因性状及品种鉴定的研究。
技术实现要素：
针对现有技术的不足，本发明提供一种豫西脂尾羊snp标记及其筛选方法和应用。本发明基于snp位点对豫西脂尾羊品种的真实性进行鉴定，克服了表型性状特征进行品种鉴定的缺陷性。
本发明的技术方案如下：一种豫西脂尾羊snp标记，所述豫西脂尾羊的特异snp位点如表1所示：
表1
突变类型中“/”左侧表示其他绵羊在该snp位点的碱基，“/”右侧表示豫西脂尾羊在该snp位点的碱基，snp位点位置指斜线“/”左右两侧第一个碱基的位置。
本发明第二方面，还提供所述豫西肥尾羊snp标记的筛选方法，所述豫西肥尾羊snp标记的筛选方法如下：从世界范围内选出多个绵羊个体涵盖多个绵羊品种作为总样本，总样本中豫西肥尾羊个数占1～3％，并获得总样本中各绵阳个体10x的全基因组snp数据，然后从总样本中取出40％～70％个绵羊个体包括5～15只豫西肥尾羊作为实验样本，从实验样本的allele频率数据中筛选出仅在豫西肥尾羊中高频出现的snp位点；将总样本中剩余的豫西肥尾羊和其他品种绵羊做验证实验。
筛选时若位点满足在豫西肥尾羊中alt频率与其他各种绵羊在该位点的alt频率的差值在0.7及以上，则该位点即被判断为豫西肥尾羊的特有的高频snp位点，在此基础上，将样本中豫西肥尾羊作为一个群体，其他品种的绵阳作为另一个群体，计算豫西肥尾羊与其他品种绵阳之间的fst值，仅保留fst在0.9以上的snp位点，最后筛选出149个snp位点作为杜泊羊的特征snp位点。
本发明第三方面，还提供基于所述snp标记的豫西肥尾羊品种或其畜产品的鉴定方法，所述鉴定方法利用焦磷酸测序法、sanger测序法、基因芯片法、或飞行时间质谱法进行鉴定snp分型分析，鉴定供试品种是否为豫西肥尾羊。
进一步地，利用焦磷酸测序法进行snp分型分析，具体操作如下：
步骤(1)采用对应所述snp位点的pcr引物对(-f和-r)及gotaqmastermix(除了模板、引物和水的其它组份)，分别以供试品种和对照品种的基因组dna为模板进行pcr扩增；
步骤(2)应用pcr引物(-f和-r)，将供试品种和豫西脂尾羊对照品种的pcr扩增产物以及其他原料置于焦磷酸测序仪上进行焦磷酸测序；
步骤(3)测序结束后使用“snp”模式对30个豫西脂尾羊snp位点进行分型分析，比较供试品种和对照品种在所述snp位点上的分型情况，只要有一个snp位点分型相同即为同一品种，否则，为不同品种。
步骤(1)中pcr扩增反应体系为：总体积为50ul，其中：mastermix(除了模板、引物和水的其它组份，提供反应环境)25ul，10umol/lpcr引物各1ul，dna模板2ul，用灭菌去离子水补齐至50ul。
进一步地，利用基因芯片法进行snp分型分析，具体操作如下：
步骤a、取实验动物，耳静脉取血，从血样中提取基因组dna。
步骤b、取步骤a得到的基因组dna，与固有30个探针的核酸芯片进行杂交；
步骤c、完成步骤b后，将核酸芯片中的各个点进行末端延伸，从而获知基因组dna中149个snp位点的分型；如果有超过40个snp位点分型相同即为同一品种，否则为不同品种。
本发明的特点如下：本发明首先通过全基因组范围内的snp数据筛选得到杜泊羊特异的snp位点，为保证所得位点的特异性，本发明实施例中用以分析比较的绵羊个体达800多只，且样本品种覆盖全世界主要大洲超过106个品种；得到的杜泊羊的特异性位点准确性高，基于所述snp位点可以快速进行杜泊羊个体及其肉质来源的品种鉴定，鉴定时绵羊供试品种有100个及以上的snp位点分型与豫西肥尾羊相同即为同一品种，否则为不同品种。基于所述snp位点和出现概率，确定了有高准确度的判断依据并经过了严格的检验，用以进行杜泊羊品种鉴定，克服了表型性状特征进行品种鉴定的缺陷性。
作为优选，判断检验样本是否为豫西肥尾羊时，检验样本在所述豫西肥尾羊snp标记位点中有高于100个alt阳性位点即可判定样本个体为豫西肥尾羊。
本发明第四方面，还提供所述豫西肥尾羊snp标记在生物遗传多样性或物种适应性进化研究中的应用。
与现有技术相比，本发明具有以下有益效果：本发明首先筛选得到豫西脂尾羊可信度极高的特异的snp位点，并基于所述snp位点和出现概率，确定了有高准确度的判断依据并经过了严格的检验，用以进行豫西脂尾羊品种鉴定；
本发明的品种真实性鉴定方法与常规测定方法的对比，其操作简便、无位点间的干扰，适于高通量操作；本发明以焦磷酸测序技术为基础，具有很高的准确性；快速高效，5h之内即可完成品种鉴定工作；结果判定简单：鉴定结果可以直接从可见光信号峰值读出。
具体实施方式
下面结合具体实施例对本发明做进一步说明，但本发明并不局限于以下技术方案。
以下实施例中涉及的snp位点筛选方法如下：从世界范围内选出885个绵羊个体涵盖106个绵羊品种，并获得其10x的全基因组snp数据。其中的豫西脂尾羊有25只，然后从432个绵羊个体包括15只豫西脂尾羊的allele频率数据中，筛选出仅在豫西脂尾羊中高频出现的snp位点，筛选策略如下：若位点满足在豫西脂尾羊中alt频率与其他各种绵羊在该位点的alt频率的差值在0.7及以上，则该位点即被判断为豫西脂尾羊的特有的高频snp位点，基于此种方法筛选出286个符合要求的snp位点。在此基础上，将总样本中25只豫西脂尾羊作为一个群体，其他个品种的绵阳作为另一个群体，计算二者之间的fst值，我们仅保留fst在0.9以上的点，故最后筛选出149个snp位点作为豫西脂尾羊的特征snp。验证：将另外的2只豫西脂尾羊与453只其他品种绵羊做验证实验，结果为：2只豫西脂尾羊的alt阳性位点均高于100，其他品种绵羊的alt阳性位点均少于20。另依据筛选出的位点在453只绵羊数据集中做pca分析，结果显示豫西脂尾羊群体被明显区分。
实施例1焦磷酸法测定豫西脂尾羊品种
(1)试剂：promega公司生产的gotaqmastermix溶液；上海生工合成的特异性扩增引物和测序引物；biotage公司生产的sepharosebead。
(2)扩增反应体系及扩增程序：扩增反应的总体积为50μl，其各种成分分别为：2×gotaqmastermix25μl，10μmol/l引物各1μl，对供试品种和对照品种dna模板2μl(黄瓜种子dna)，用灭菌去离子水补齐至50μl；反应程序：95℃预变性10min，94℃变性30s，50℃退火30s，72℃延伸45s，循环50次，最后72℃延伸7min，4℃保存。
(3)测序反应体系：50μlpcr产物中加入47μlbindingbuffer和3μlsepharosebeads，1300rpm涡旋混匀15min，经vacuumprepworkstation单链分离，释放到预先加入38.8μlannealingbuffer和1.2μl测序引物的psq96测序反应板中，80℃金属加热块上放置2min后冷却至室温；将酶(dna聚合酶、atp硫酸化酶、荧光素酶和三磷酸腺苷双磷酸酶)、底物和a、t、c、g成分加入试剂仓，即可上机测序。
(4)测序时间8～10min，测序结束后使用“snp”模式对供试品种和对照品种的30个绵羊snp位点进行分型分析，比较供试品种和对照品种在149个snp位点上的分型情况。对照品种和供试品种1为豫西脂尾羊，供试品种2为豫西脂尾羊杂种羊。结果表明，供试品种1与对照品种在87个snp分型结果上相同，所以供试品种与对照品种为同一品种，即是豫西脂尾羊。供试品种2与对照品种在24个snp分型结果上相同，所以供试品种2是混有豫西脂尾羊血统的杂种羊。
焦磷酸法测序时：测序引物与pcr扩增的dna模板相结合。然后将其与dna聚合酶、atp硫酸化酶、荧光素酶和三磷酸腺苷双磷酸酶，以及底物aps和荧光素一起孵育；四种dntp(datp，dttp，dctp，dgtp)之一被加入反应体系，如与模扳配对(a—t，c—g)，此dntp与引物的末端形成共价键，dntp的焦磷酸基团(ppi)释放出来。而且释放出来的ppi的量与和模板结合的dntp的量成正比；atp硫酸化酶在三磷酸腺苷双磷酸酶存在的情况下催化ppi形成atp，atp驱动荧光素酶介导的荧光素向氧化荧光素的转化，氧化荧光素发出与atp量成正比的可见光信号，光信号由ccd摄像机检测，并由可见光信号峰表示。每个光信号的峰高与反应中掺入的核苷酸数目成正比；然后加入下一种dntp。最终待测序列顺序，即可从反应光强的信号峰中读出。
实施例2pcr扩增法测定豫西脂尾羊品种
随机选择24只阿勒泰羊，24只巴音布鲁克羊和29只豫西脂尾羊作为样品个体，采集这些羊的血液，使用肝素纳抗凝，采取天根血液基因组试剂盒提取基因组dna，0.8％琼脂糖凝胶电泳检测，核酸蛋白定量仪测定其浓度和纯度，稀释至100ng/ul，保存备用。采用表2所列引物以待测样本基因组dna为模板进行pcr扩增，其中pcr扩增反应体系为：加入0.4μl的taqdnapolymerase(5u/μl)，5μl的10×pcrbuffer，4μl的mgcl2(25mmol/μl)，1μl的forwardprimer(10pmol/μl)，1μl的reverseprimer(10pmol/μl)，4μl的dntp(2.5mmol/μl)，2μl的模板dna和32.6μl的ddh2o，共50μl。
采用琼脂糖凝胶dna回收试剂盒进行pcr产物的回收与纯化。收集产物进行测序，之后对测序结果进行进行序列匹配，找到snp位置并分析，确定snp类型。比较待测样本和对照品种在149个snp位点上的分型情况，最终29份血液样本在149个位点上至少有110个以上的匹配，最多达到147个匹配；其余48份样本全无位点匹配，对于豫西脂尾羊的鉴定的准确性均达到了100％，其中全部匹配为豫西脂尾羊，否则，不是豫西脂尾羊。
实施例3基因芯片法测定豫西脂尾羊品种
步骤1、取实验动物，耳静脉取血，从血样中提取基因组dna。
步骤2、取步骤1得到的基因组dna，与固有149个探针(30个探针分别为序列表2的序列1至序列30所示的单链dna分子)的核酸芯片进行杂交。
步骤3、完成步骤2后，将核酸芯片中的各个点进行末端延伸，从而获知基因组dna中149个snp位点的基因型。
步骤4、取步骤1得到的基因组dna，进行全基因组测序，获知基因组dna中149个snp位点的基因型。结果表明，步骤3的结果和步骤4的结果完全一致。豫西脂尾羊的48头实验动物的149个snp位点中能够至少匹配上99-136个位点。非豫西脂尾羊的16头实验动物的149个snp位点匹配个数都低于20。