一种家猪家系划分及亲子鉴定试剂盒及方法

no.11所示，反向引物序列如seq id no.12所示；s0101微卫星座正向引物序列如seq id no.13所示，反向引物序列如seq id no.14所示；s0225微卫星座正向引物序列如seq id no.15所示，反向引物序列如seq id no.16所示；s0386微卫星座正向引物序列如seq id no.17所示，反向引物序列seq id no.18所示；s0090微卫星座正向引物序列如seq id no.19所示，反向引物序列如seq id no.20所示；s0355微卫星座正向引物序列如seq id no.21所示，反向引物序列如seq id no.22所示；s0026微卫星座正向引物序列如seq id no.23所示，反向引物序列如seq id no.24所示；sw24微卫星座正向引物序列如seq id no.25，反向引物序列如seq id no.26所示；sw2156微卫星座正向引物序列如seq id no.27所示，反向引物序列如seq id no.28所示。
8.进一步地位点s0155、sw963、s0386荧光标记为5`
‑
fam，位点sw240、sw72、sw122、sw24、s0026荧光标记为5`
‑
hex，位点s0227、s0101、s0090、s0355荧光标记为5`tamra；位点s0155、sw240、sw72、s0227、sw963、s0101、s0225、s0386、s0090、s0355、sw24的pcr退火温度为55℃，位点sw122的pcr退火温度为58.5℃，位点sw2156退火温度为60℃。
9.进一步地，所述试剂盒为25μl反应体系，包括1μl上游引物，1μl下游引物，1μl mix dntr，2.5μl包括mgcl2taq buffer，0.5μl taq酶，19μl ddh2o。
10.进一步地，所述14个ssr位点可应用于安庆六白猪。
11.本发明还提供一种家猪家系划分及亲子鉴定方法，以上述试剂盒合成ssr引物进行pcr从而进行家系划分及亲子鉴定。
12.进一步地，鉴定方法具体步骤如下所述：(1)样品采集；(2)dna提取：采用北京艾德莱生物科技公司组织基因组dna快速提取试剂盒进行猪dna提取；(3)ssr引物序列合成；(4)pcr扩增；(5)pcr产物毛细管电泳检测；(6)数据处理及结果分析。
13.进一步地，步骤(6)数据处理及结果分析具体步骤为：使用popgene3.2软件对微卫星位点多态性进行等位基因数目、有效等位基因数、观测纯和度、期望纯和度、观测杂合度、期望杂合度、香农指数、hard
‑
weinberg平衡指标分析；使用pic_cale0.6软件进行多态信息含量计算；利用diploiddate程序对个体间的nei氏遗传距离进行计算并获得，以非加权组平均法(unweighted pair
‑
group method with arithmetic mean，upgma)构建的亲缘聚类图。
14.进一步地，所述家猪为安庆六白猪；所述dna采集位点位于猪耳。
15.进一步地，步骤(4)pcr过程为：、95℃预变性3分钟；30秒94℃变性，30秒60℃退火，30秒72℃延伸，10个循环；30秒94℃变性，30秒55℃退火，30秒72℃延伸，35个循环；8分钟72℃修复延伸。
16.本发明具有如下优点：
17.1.针对安庆六白猪家系划分和亲子鉴定方面的空白，特异性地选择14个ssr位点，通过合成引物进行pcr，获取ssr标记信息的家系划分与生产记录相符，因此，可以通过分子鉴定快速获得安庆六白猪的家系判断。选择的14个微卫星位点的ne
‑
1p、ne
‑
2p、ne
‑
pp累积排除率均在99％以上，满足生产实践中不同个体间亲子鉴定要求。
18.2.选择的14个位点分属不同染色体，彼此间连锁可能性较小；14个微卫星位点在安庆六白猪群体中等位基因型分布具有特异性。
19.3.本方法dna采集位点为猪耳，可以以最小伤害获得猪的dna，有利于生产实践。根
据本方法基于14个ssr位点合成引物对进行pcr，可以提高家猪尤其是安庆六白猪育种准确性，加快其育种进程，具有良好的应用前景和经济效益。
附图说明
20.图1荧光标记5`
‑
fam下s0155位点65号样本电泳结果
21.图2荧光标记5`
‑
hex下sw122位点8号样本电泳结果
22.图3荧光标记5`
‑
tamra下s0090位点95号样本电泳结果
具体实施方式
23.下面结合具体实施例对本发明进行具体阐述。
24.实施例1：
25.1材料方法：
26.1.1样品采集
27.本发明共采集98头份纯种安庆六白猪耳组织样本，由安徽省花亭湖绿色食品开发有限公司提供(安徽省安庆市太湖县晋熙镇程岭村)。耳样保存于装有500μl 75％酒精的1.5ml离心管中，带回实验室后于
‑
20℃条件下保存。
28.1.2 dna提取
29.安庆六白猪组织dna采用艾德莱生物科技公司(北京)组织基因组dna快速提取试剂盒进行提取，提取后的dna样品使用紫外分光光度计进行浓度及纯度检查，质量合格的dna样品存放在
‑
80℃环境下备用。
30.1.3微卫星位点选择及引物设计
31.本发明选择的14个位点信息如表1所示，根据位点信息合成对应位点引物序列，位点信息及引物序列详见表1。
32.表1微卫星位点及引物信息
33.[0034][0035]
1.4 pcr扩增体系及反应条件
[0036]
pcr反应采用25μl反应体系，包括1μl上游引物，1μl下游引物，1μl dntr(mix)，2.5μl taq buffer(with mgcl2)，0.5μl taq酶，19μl ddh2o；pcr程序设定为：95℃预变性(3min)，10循环(94℃变性(30sec)，60℃退火(30sec)，72℃延伸(30sec))，35循环(94℃变性(30sec)，55℃退火(30sec)，72℃延伸(30sec))，72℃修复延伸(8min).
[0037]
1.5 pcr产物毛细管电泳检测
[0038]
使用abi3730xl基因测序仪对pcr产物进行毛细管电泳，使用genemapper软件读取分析毛细管电泳结果，判断各样本在每个微卫星位点上的等位基因大小及基因型。
[0039]
1.6数据处理及分析
[0040]
整理各微卫星位点在安庆六白猪个体中基因型。使用popgene3.2软件对微卫星位点多态性进行以下信息分析：等位基因数目(number of alleles)，有效等位基因数(effective number of alleles)，观测纯和度(observed homozygote)，期望纯和度(expected homozygote)，观测杂合度(observed heterozygosity)、期望杂合度(expected heterozygosity)、香农指数(shannon's information index)、hard
‑
weinberg平衡等指标(hard
‑
weinberg equilibrium test)，使用pic_cale0.6软件进行多态信息含量计算(polymorphism information content)。本发明将每个样本视为一个群体，利用diploiddate程序对个体间的nei氏遗传距离进行计算，并获得以非加权组平均法(unweighted pair
‑
group method with arithmetic mean，upgma)构建的亲缘聚类图。使用cervus 3.07软件计算14个微卫星位点无效等位基因频率(nullallele frequency)及非排除率(combined exclusion probability)，排除率计算方法为1
‑
非排除率。
[0041]
2结果与分析
[0042]
2.1安庆六白猪14个微卫星位点遗传多样性分析
[0043]
本发明共选用14个微卫星位点对98个安庆六白猪个体进行遗传多样性分析。部分
微卫星位点pcr产物毛细管电泳结果如图1、2、3所示：如65号样本在s0155位点电泳结果只存在一个峰(161bp)，说明为纯合子；8号样本sw122位点电泳结果存在两个峰(111/118bp)，为杂合子；95号样本在s0090位点存在两个峰为杂合子。
[0044]
各微卫星位点在安庆六白猪群体中等位基因及基因型频率平见表2(按片段大小从低到高排序分别命名为等位基因a至等位基因l)。
[0045]
表2 14个微卫星位点在安庆六白猪群体中等位基因及基因频率统计
[0046][0047][0048]
14个位点共检测出119个等位基因，平均每个位点等位基因数为8.5，其中s0255位点等位基因数最多(n＝12)，片段范围为167bp
‑
191bp，sw72位点等位基因数最少(n＝5)，片段范围为99bp
‑
113bp。在14个位点中，sw72位点b等位基因(100bp)在安庆六白猪群体中等位基因频率最高为0.4323；s0255位点l等位基因(191bp)在安庆六白猪群体中等位基因频率最低为0.0052。各位点等位基因频率大于0.1的等位基因型数量为2
‑
5个。
[0049]
安庆六白猪14个微卫星位点遗传多样性分析见表3。各位点有效等位基因数处于3.1927
‑
7.1078，平均有效等位基因数为4.8161，其中s0026位点有效等位基因数最少，sw2156位点有效等位基因数最多。各位点观测纯合度处于0.1684
‑
0.6421之间，其中sw240位点观测纯合度值最小，sw963位点观测纯合度值最大。各位点期望纯合度处于0.3579
‑
0.8316，其中sw963位点期望纯合度值最小，sw240位点期望纯合度值最大。各位点观测杂合度范围在0.1358
‑
0.3097之间，sw2156位点观测杂合度值最小，s0026位点观测杂合度值最
大。各位点期望杂合度范围在0.6903
‑
0.8626之间，s0026位点期望杂合度值最低，s0255位点期望杂合度值最高。
[0050]
表3 14个微卫星位点在安庆六白猪群体中遗传信息分析
[0051][0052][0053]
注：na(观察等位基因数)；ne(有效等位基因数)；hom(obs)观测纯合度；hom(exp)期望纯合度；het(obs)观测杂合度；het(exp)期望杂合度；i：香农指数；pic多态信息含量；h
‑
w：hard
‑
weinberg平衡检测，**表示位点未达到h
‑
w动态平衡，无**标记表示位点已达到h
‑
w动态平衡。
[0054]
各位点香农指数处于1.3263
‑
2.1314之间，sw72位点的香农指数值最低，sw2156位点的香农指数值最高。各位点多态信息含量范围在0.6444
‑
0.8424之间，s0026位点的多态信息含量值最低，s0255位点的多态信息含量最高。各位点经过卡方检验以研究其是否处于hard
‑
weinberg动态平衡，经计算s0155，s0227，sw963，sw122，s0255，s0386，s0355，sw2156不处于hard
‑
weinberg动态平衡；sw240，sw72，s0101，s0090，sw24，s0026已处于hard
‑
weinberg动态平衡。
[0055]
2.2基于ssr标记的家系划分及亲缘鉴定
[0056]
根据nei氏遗传距离运算结果构建98头安庆六白猪的亲缘关系树，可发现本实施例样本可分为7个类群，与该保种场保种群的家系组成记录相符。部分样本间的nei氏遗传距离、遗传一致度系数见表4，结果可知遗传一致度系数与遗传距离系数之间呈负相关。
[0057]
表4 1
‑
12样本nei氏遗传距离及遗传一致系数
[0058][0059][0060]
注释：表上方为遗传一致度系数，表下方为遗传距离系数。
[0061]
14个微卫星位点非排除率见表5。当一个候选亲本和子代的基因型已知(ne
‑
1p)，各位点排除率在0.276
‑
0.419，14个位点累计排除率为0.9995；一个已知亲本、性别不同的另一个候选亲本和子代的基因型已知(ne
‑
2p)，各位点排除率在0.448
‑
0.713，14个位点累计排除率为0.9999；候选父母双方和子代基因型已知(ne
‑
pp)，各位点排除率在0.626
‑
0.887，14个位点累计排除率为0.9999。
[0062]
本发明选取的14个微卫星位点在安庆六白猪群体中共观测到119个等位基因，每个位点平均等位基因数为8.500
±
2.1031；平均有效等位基因数为4.8162
±
1.3217。安庆六白猪微卫星位点平均有效等位基因数高于安徽省其他地方猪种(圩猪、皖南黑猪、皖南花猪)，也高于部分省外地方猪种如枫泾猪、撒坝猪、大河猪、巴马香猪。检测各位点的ne/na比值，发现s0355位点ne/na比值最高，为0.7138；sw24位点比值最低，为0.3614，说明s0355位点各等位基因在安庆六白猪群体中分布较为均匀。此外，s0227位点在安庆六白猪群体中的等位基因数为7个，在泰国地方猪中只检测出4个等位基因，安庆六白猪其他微卫星位点与墨西哥无毛猪、古巴克里奥尔猪、伊比利亚猪等猪种相比，等位基因数也有所差异。这表明14个位点等位基因在安庆六白猪群体中分布存在特异性。
[0063]
表5 14个微卫星位点排除概率
[0064]
locine
‑
1pne
‑
2pne
‑
ppne
‑
ine
‑
sif(null)s01550.4190.60.790.9290.6250.0115sw2400.5340.70.8740.960.6620.0019sw720.2780.4480.6260.8580.564
‑
0.0080
s02270.3440.5220.7080.8960.5970.0741sw9630.4310.610.80.9320.6290.3810sw1220.4560.6320.8140.940.64
‑
0.0060s01010.4080.5880.7820.9230.6180.0055s02550.5520.7130.8790.9630.6690.2350s03860.3430.5240.7190.8960.5920.2846s00900.3790.560.750.9130.6110.0353s03550.4820.6550.8310.9470.650.1759sw240.3510.5370.7440.9010.5890.0410s00260.2760.4530.6430.8590.5570.0173sw21560.5620.7220.8870.9290.6250.1713累计排除率0.99950.99990.99990.99990.9999 [0065]
有效等位基因数、多态信息含量、杂合度是研究群体遗传变异的重要指标，hardy
‑
weinberg动态平衡则可反映群体的迁移，选择等情况。杂合度可以用来度量群体的变异程度，杂合度值越低，说明该群体遗传变异少，多样性低。安庆六白猪14个微卫星位点平均期望杂合度值为0.7824
±
0.0573，在国内地方猪中处于中等水平。多态信息含量可作为衡量位点多态性的指标之一，若多态信息含量值高，则代表位点在群体中杂合子的比例较高，其遗传潜力也越大。14个微卫星位点在安庆六白猪群体中多态信息含量值都高于0.5，其平均多态信息含量值为0.74
±
0.06。当某遗传位点达到h
‑
w动态平衡时，其等位基因频率及基因型频率在群体中可稳定遗传。本发明中，14个微卫星位点在安庆六白猪群体中有s0155等八个位点未达到h
‑
w动态平衡，sw240等六个位点达到h
‑
w动态平衡。本发明共选择14个微卫星位点进行多态性检测，且各位点分属不同染色体，彼此间连锁可能性较小。根据微卫星结果绘制的遗传距离聚类图可将98头安庆六白猪区分为7个家系，与其保种群构建之初的引种记录一致，说明本发明的可信度较高。本发明所选择的14个微卫星位点的ne
‑
1p、ne
‑
2p、ne
‑
pp累积排除率均在99％以上，满足生产实践中不同个体间亲子鉴定要求。
[0066]
可以知道，上述实施例仅为了说明发明原理而采用的示例性实施方式，然而本发明不仅限于此，本领域技术人员在不脱离本发明实质情况下，可以做出各种改进和变更，这些改进和变更也属于本发明的保护范围。