一种与二花脸母猪产仔性状相关的SNP标记

摘要

本发明公开了与二花脸母猪产仔性状相关的SNP标记、其检测方法及应用。所述SNP标记位于猪8号染色体上的UCHL1基因和LIMCH1基因的基因之间的间隔区，所述SNP标记的位点为国际猪基因组10.2版本参考序列猪8号染色体上g.33985796核苷酸位点，且具有T/C多态性，所述SNP标记与二花脸母猪窝产总仔数极显著相关。一种用于检测权利要求1所述的SNP标记的引物对，上游引物为：SEQ ID NO：2，下游引物为：SEQ ID NO：3。本发明提供的SNP标记与二花脸母猪的产仔性能相关，因此，可以通过鉴定该SNP标记来筛选高产的二花脸母猪品系，所得的二花脸母猪高产品系具有重要的经济效益与社会价值。

权利要求

1.一种用于检测与二花脸母猪产仔性状相关的SNP标记的引物对在筛选高产仔二花脸母猪品系中的应用，所述的引物对的上游引物为：SEQ ID NO：2，下游引物为：SEQ ID NO：3，所述SNP标记位于猪8号染色体上的UCHL1基因和LIMCH1基因的基因之间的间隔区，所述SNP标记的位点为国际猪基因组10.2版本参考序列猪8号染色体上g. 33985796核苷酸位点，且具有T/C多态性，所述SNP标记与二花脸母猪窝产总仔数极显著相关。

2.一种筛选高产仔二花脸母猪品系的方法，其特征在于包括检测二花脸母猪国际猪基因组10.2版本参考序列猪8号染色体上g. 33985796核苷酸位点的基因型，该SNP位点具有TT基因型的二花脸个体母猪窝总产仔数显著高于具有CC基因型的二花脸个体母猪窝总产仔数，选育该位点的TT型个体作为种猪。

说明书

一种与二花脸母猪产仔性状相关的SNP标记

技术领域

[0001] 本发明属于分子生物学技术领域，涉及与二花脸母猪产仔性状相关的SNP标记、其检测方法及应用。

背景技术

[0002] 猪产仔数是重要的经济性状，产仔数的提高将会大大地增加商品猪肉的供应量，给现代养猪生产带来巨大的经济效益。随着我国经济持续快速的发展，人们生活水平逐渐提高，对猪肉的需求量也越来越大。2014年12月底我国能繁母猪存栏数4300万头，如果母猪平均每胎多产仔1头，按每年每头母猪产2胎计算，则每年可为整个行业多提供约9000万头猪。因此，人们也越来越关注如何提高猪的产仔数性能。然而，产仔数为复杂多基因控制的经济性状，遗传力低，通过传统的选育方法提高猪产仔数收效甚微，所以开发和利用新的分子选育标记来提高猪的繁殖性能倍受重视。

[0003] 二花脸是位于我国太湖流域的国家级畜禽遗传资源保护品种，是我国劳动人民千百年来选育的产仔性能极其优秀的地方猪种资源。根据发明人所在单位对二花脸猪群最大基地“焦溪二花脸猪专业合作社”内177头二花脸约1000窝产仔数据分析来看，二花脸群体内产仔性能已出现显著分离，尤其是初产总产仔数和产活仔数变异系数分别达到22.14％和26.97％；经产也分别为19.87％和22.14％。然而目前导致二花脸猪品种内产仔数变异的遗传机理尚不清楚。

[0004] 多年以来虽然国内已有大量的研究机构利用二花脸来鉴别二花脸猪种高繁殖力的遗传机理，但限于研究方法、手段、材料及繁殖性状本身复杂性等众多因素的制约，二花脸猪种高繁殖力遗传机制并没有得到充分的揭示和有效的利用。因此，为保护好我国人民千百年来选育的二花脸品种的高产仔性能，我们急需鉴别出影响二花脸猪品种内产仔数变异的基因和标记，通过分子标记选择技术加快恢复和提升二花脸猪种的产仔性能，培育性能相对稳定的高产群体、巩固太湖流域这些地方猪种的高产优势。

[0005] 发明人利用一个177头纯种二花脸母猪的群体，开展了全基因组QTL定位研究，并在8号染色体上定位到了一个品种内群体间遗传分化显著的QTL。从国际猪QTL数据库网站(http://www.animalgenome.org/cgi‐bin/QTLdb/SS/index)可知，目前在猪除10、11号染色体及性染色体上没有定位到影响总产仔数的QTL，其它常染色体上都已定位到影响总产仔数的QTL，但大部分是利用微卫星标记定位的QTL，置信区间多在10‐20cM，无法确定真正的主效基因及其关键变异位点，因此难以直接应用于种猪选育改良。

发明内容

[0006] 本发明的目的在于针对现有技术不足，产仔数的低遗传力，提供与母猪窝总产仔数相关的SNP标记。

[0007] 本发明的另一个目的在于提供用于检测上述SNP标记的引物和检测方法。

[0008] 本发明的另一个目的在于提供上述SNP标记的用途。

[0009] 一种与二花脸母猪产仔性状相关的SNP标记，所述SNP标记位于猪8号染色体上的UCHL1基因和LIMCH1基因的基因之间的间隔区，所述SNP标记的位点为国际猪基因组10.2版本参考序列猪8号染色体上g.33985796核苷酸位点，且具有T/C多态性，所述SNP标记与二花脸母猪窝产总仔数极显著相关。该SNP位点具有TT基因型的二花脸个体母猪窝总产仔数显著高于具有CC基因型的二花脸个体母猪窝总产仔数。

[0010] 一种基于本发明所述的SNP开发分子标记的方法，以含有本发明所述的SNP标记的核苷酸序列为基础序列，设计引物对，以二花脸母猪基因组DNA为模板进行PCR扩增，使本发明所述的SNP标记转化为分子标记。

[0011] 其中，所述的引物对序列优选为上游引物：SEQ ID NO：2，下游引物：SEQ ID NO：3；所述的分子标记序列优选如SEQ ID NO：1所示，所述的SNP位点位于第806位，存在T/C多态性。

[0012] 按照本发明所述的方法得到的分子标记。

[0013] 所述的分子标记序列优选如SEQ ID NO：1所示，所述的SNP位点位于第806位，存在T/C多态性。

[0014] 一种用于检测权利要求1所述的SNP标记的引物对，上游引物为：SEQ ID NO：2，下游引物为：SEQ ID NO：3。

[0015] 一种检测本发明所述的SNP标记的方法，包含PCR扩增二花脸母猪基因组中含有本发明所述的SNP标记的一段序列，对扩增产物进行测序，判读该位点的T/C多态性。

[0016] 本发明所述的检测本发明所述的SNP标记的方法优选包括以下步骤：

[0017] (1)取一头二花脸母猪的耳组织样品并提取总DNA；

[0018] (2)用已提取的二花脸母猪基因组DNA为模板，使用本发明所述的引物进行PCR扩增；

[0019] (3)扩增产物进行测序，分析测序结果，判读在SEQ ID NO：1第806位的T/C多态性。

[0020] 其中，步骤(2)所述的PCR扩增反应体系为：DNA模板2.5μL、SEQ ID NO：2和SEQ ID NO：3所示的引物各1.25μL、PCR Mix试剂25μL、双蒸水20μL；其中所述DNA模板浓度为30ng/μL，所述引物的浓度为10mol/L，所述PCR Mix试剂为南京欧科生物技术有限公司的P394961L型号试剂；PCR扩增的反应程序为：预变性96℃2min；变性96℃20s；退火62℃30s，延伸72℃45s，35个循环；延伸72℃10min。

[0021] 本发明所述的SNP标记、所述的分子标记、所述的引物在筛选高产二花脸母猪品系中的应用。

[0022] 一种筛选高产二花脸母猪品系的方法，包括检测二花脸母猪g.33985796核苷酸位点的基因型，选育g.33985796核苷酸位点的TT型个体作为种猪。

[0023] 有益效果：

[0024] 本发明提供的SNP标记与二花脸母猪的产仔性能相关，因此，可以通过鉴定该SNP标记来筛选高产的二花脸母猪品系，所得的二花脸母猪高产品系具有重要的经济效益与社会价值。

附图说明

[0025] 图1为高产与相对低产二花脸群体间染色体上SNP标记Fst值的分布情况。其中，猪的18条常染色体及X染色体信息标识于X轴。

[0026] 图2为使用本发明的引物扩增g.33985796的电泳图。

[0027] 图3为g.33985796的突变位点不同基因型的DNA测序结果峰图。

具体实施方式

[0028] 以下实施例用于说明本发明，但不用来限制本发明的范围。在不背离本发明精神和本质的情况下，对本发明方法、步骤或条件所作的修改或替换，均属于本发明的范围。

[0029] 实施例1

[0030] 1、实验动物来源

[0031] 江苏常州焦溪二花脸猪专业合作社。

[0032] 计算177头二花脸母猪的育种值，计算模型为

[0033] Y(窝产总仔数)＝parity(胎次)+farm(场)+year(年度)+season(季节)+age(母猪分娩年龄)+sire(与配公猪)+permanent effect(母猪的永久效应)+additive effect(个体的加性效应)+e(残差)，

[0034] 其中包括固定效应‐胎次、母猪分娩的场/年/季，协变量‐母猪分娩的年龄，随机[0035] 效应‐与配公猪，永久效应‐母猪，个体加性遗传值。选择育种值排在前10％与后10％的个体。

[0036] 2、提取基因组DNA

[0037] 采集36头母猪的耳组织样品，放置于装有70％酒精的离心管内，‐20℃冰箱保[0038] 存备用。

[0039] 使用传统酚/氯仿法提取耳组织基因组DNA，所需试剂包括：

[0040] 裂解液实验室配备

[0041] 蛋白酶K(德国MERCK生物科技有限公司)

[0042] Tris饱和酚(北京索莱宝生物科技有限公司)

[0043] Tris饱和酚：氯仿：异戊醇(25：24：1)(北京索莱宝生物科技有限公司)[0044] 氯仿(江苏永华精细化学品有限公司)

[0045] 无水乙醇(广东光华科技股份有限公司)

[0046] 3M乙酸钠(北京索莱宝生物科技有限公司)

[0047] 具体步骤如下所述：

[0048] (1)取黄豆大小组织样，尽量剪碎放入2ml离心管中；

[0049] (2)加入裂解液(自己配备)800μL，及蛋白酶K 30μL(0mg/ml)；

[0050] (3)样品置于55℃恒温箱中孵育过夜，至管中无组织块为止；

[0051] (4)加入Tris饱和酚800μL，轻微混匀10min，4℃12000r/min离心12min；

[0052] (5)取650μL上清加Tris饱和酚：氯仿：异戊醇(25：24：1)800μL，混摇10min，4℃12000r/min离心12min；

[0053] (6)取550μL上清，加氯仿800μL，混摇10min，4℃12000r/min离心12min；以下步骤换1.5ml的离心管

[0054] (7)取450μL上清，加无水乙醇800μL，3M乙酸钠40μL，混摇6min，4℃1000r/min离心8min；

[0055] (8)弃上清留下DNA沉淀团，加入1000μL 70％乙醇(自己配备)，混摇5min，4℃1000r/min离心5min，弃上清(如需要可重复一次)；

[0056] (9)将离心管放入通风橱，吹干至管内无小滴；

[0057] (10)样品加100μL超纯水，轻微吹打至DNA溶解，经过Nanodrop‐100分光光度计检测质量与浓度后将浓度同一稀释到50ng/μL于‐20℃下保存备用。

[0058] 3、猪全基因组6万个(60K)SNP基因型检测

[0059] 上述个体的DNA在Illumina Beadstation平台上根据公司标准流程进行猪全基因组60K SNP(Illumina，美国)基因型判定。利用PLINK(1.9)对所有样本60K芯片数据进行质量控制，剔除检出率低于0.95、家系孟德尔错误率高于0.05的个体；最小等位基因频率小于0.05的SNP标记。

[0060] 4、高产和相对低产群体Fst值的计算

[0061] 使用Powermarker软件包对对分型个体的60K SNP标记型数据进行处理，计算得到两个群体的每一个SNP位点遗传分化系数Fst值。结果显示，Fst值最高点位于猪8号染色体上，该SNP标记极有可能与总产仔数性状相关(图1)。

[0062] 实施例2

[0063] 本实施例为实施例1中得到的SNP位点g.33985796T/C在二花脸母猪群体内的验证。

[0064] 1、提取二花脸母猪基因组DNA

[0065] 采集具有准确窝总产仔数记录的195头纯种二花脸母猪的耳组织样品，放置于装有

[0066] 70％酒精的离心管内，‐20℃冰箱保存备用。利用上述方法提取耳组织基因组DNA，经过质量、浓度检测后将浓度稀释到30ng/μL于‐20℃下保存备用。

[0067] 2、目的片段PCR扩增与测序

[0068] 用已提取的DNA为模板，根据所设计的引物，进行PCR扩增：取DNA模板2.5μL、SEQ ID NO：2和SEQ ID NO：3所示的引物各1.25μL、PCR Mix试剂25μL、双蒸水20μL；设置PCR扩增体系：预变性96℃2min；变性96℃20s；退火62℃30s；延伸72℃45s；35个循环；然后延伸10min。

[0069] PCR产物在1.2％琼脂糖凝胶中电泳检测，扩增的目的片段大小为996bp，电泳图见图2，将剩余的扩增产物进行测序，测序结果用DNAman软件与GenBank中猪的相关基因片段序列比对、分析，判读g.33985796T/C的基因型，然后利用SAS软件进行基因型对表型的影响效应分析。分析模型为Yijklm＝u+Gj+Bk+Pl+eijklm

[0070] 其中：Yijkm为猪的产仔数；Gj代表第j个SNP的基因型固定效应；Bk是第k个批次固定效应；Pl是胎次的随机效应，不同胎次的产仔数记录作为重复数据处理；eijklm为残差。

[0071] 显著性的P值经过10000次的随机抽样校正。

[0072] 表1给出了g.33985796T/C突变位点在纯种二花脸群体中对窝总产仔数的影响效应。由表1可知，纯种二花脸猪中，g.33985796T/C位点的TT基因型个体与CC型个体相比较：窝总产仔数平均增加1.87头。由此可见，在二花脸猪种中，继代选育g.33985796T/C位点的TT型个体，均可逐步提高二花脸母猪的窝总产仔数，达到提高二花脸母猪繁殖性能的目的。

[0073] 表1 g.33985796T/C SNP位点与二花脸母猪窝总产仔数的关联分析

[0074]