用于猪亲缘关系鉴定的snp分子标记组合、应用及方法
技术领域
1.本技术涉及分子遗传学技术领域，尤其涉及一种用于猪亲缘关系鉴定的snp分子标记组合、应用及方法。


背景技术：

2.系谱信息在育种中有着非常重要的作用，系谱错误会大大降低遗传评估的准确性，影响选中选配的效果，降低群体遗传进展，给企业带来巨大的经济影响。然而，在实际的生产中，系谱错误是难免的，它受多种因素的影响，比如配种记录、产仔记录、系谱录入及整理过程中的人为错误等等。目前针对猪亲缘关系鉴定的方法基本都是基于微卫星方法建立的，如专利cn102154280a中公开了一种猪微卫星标记的复合扩增的方法及其专用引物、专利cn106244717a中公开了12个常染色体str基因座和1个性别决定基因gstram在内的13个基因座的亲权鉴定的方法。
3.snp是继rflp和str两种遗传标记之后的第三代分子标记，在生物基因组中分布最广泛，占所有已知多态性的90％以上，具有数量多、分布广、稳定性高等优点，而且由于snp属于二态的遗传标记，更适合快速、规模化筛查。
4.目前尚无利用snp标记进行杜、长、大商业化品种猪亲缘关系鉴定的报道。


技术实现要素：

5.鉴于上述内容，第一方面，本技术提供了一种用于猪亲缘关系鉴定的snp分子标记组合，其特征在于，包括：在猪的18条常染色体上的200个snp标记，所述snp标记的snp位点信息如下：
6.编号1的snp位点位于2号染色体上第9987664位；
7.编号2的snp位点位于14号染色体上第418001位；
8.编号3的snp位点位于15号染色体上第2977317位；
9.编号4的snp位点位于7号染色体上第7738134位；
10.编号5的snp位点位于10号染色体上第2775295位；
11.编号6的snp位点位于3号染色体上第11913946位；
12.编号7的snp位点位于9号染色体上第7629826位；
13.编号8的snp位点位于4号染色体上第7720753位；
14.编号9的snp位点位于16号染色体上第2857298位；
15.编号10的snp位点位于8号染色体上第1653890位；
16.编号11的snp位点位于5号染色体上第5646887位；
17.编号12的snp位点位于6号染色体上第2308923位；
18.编号13的snp位点位于17号染色体上第2842639位；
19.编号14的snp位点位于12号染色体上第1779292位；
20.编号15的snp位点位于11号染色体上第4655543位；
21.编号16的snp位点位于13号染色体上第7406226位；
22.编号17的snp位点位于18号染色体上第23529836位；
23.编号18的snp位点位于1号染色体上第1166570位；
24.编号19的snp位点位于8号染色体上第1735070位；
25.编号20的snp位点位于15号染色体上第126780202位；
26.编号21的snp位点位于13号染色体上第9797897位；
27.编号22的snp位点位于1号染色体上第277024546位；
28.编号23的snp位点位于12号染色体上第9778935位；
29.编号24的snp位点位于6号染色体上第10774932位；
30.编号25的snp位点位于9号染色体上第60564642位；
31.编号26的snp位点位于6号染色体上第11560736位；
32.编号27的snp位点位于17号染色体上第13605151位；
33.编号28的snp位点位于2号染色体上第10411684位；
34.编号29的snp位点位于6号染色体上第11229072位；
35.编号30的snp位点位于4号染色体上第15870194位；
36.编号31的snp位点位于17号染色体上第13261924位；
37.编号32的snp位点位于12号染色体上第9585887位；
38.编号33的snp位点位于8号染色体上第11917230位；
39.编号34的snp位点位于1号染色体上第3659113位；
40.编号35的snp位点位于1号染色体上第13677911位；
41.编号36的snp位点位于2号染色体上第45279469位；
42.编号37的snp位点位于16号染色体上第75677205位；
43.编号38的snp位点位于14号染色体上第1472814位；
44.编号39的snp位点位于17号染色体上第35417485位；
45.编号40的snp位点位于5号染色体上第89085559位；
46.编号41的snp位点位于6号染色体上第110566592位；
47.编号42的snp位点位于5号染色体上第26784434位；
48.编号43的snp位点位于6号染色体上第12847266位；
49.编号44的snp位点位于8号染色体上第93655842位；
50.编号45的snp位点位于14号染色体上第44862892位；
51.编号46的snp位点位于18号染色体上第24878198位；
52.编号47的snp位点位于18号染色体上第33999857位；
53.编号48的snp位点位于14号染色体上第2935893位；
54.编号49的snp位点位于10号染色体上第15659641位；
55.编号50的snp位点位于2号染色体上第126441545位；
56.编号51的snp位点位于1号染色体上第13005755位；
57.编号52的snp位点位于4号染色体上第15411794位；
58.编号53的snp位点位于1号染色体上第3067609位；
59.编号54的snp位点位于9号染色体上第151103363位；
60.编号55的snp位点位于2号染色体上第93254715位；
61.编号56的snp位点位于8号染色体上第14014245位；
62.编号57的snp位点位于7号染色体上第9853317位；
63.编号58的snp位点位于13号染色体上第16580025位；
64.编号59的snp位点位于2号染色体上第17899277位；
65.编号60的snp位点位于13号染色体上第9626564位；
66.编号61的snp位点位于1号染色体上第71814074位；
67.编号62的snp位点位于1号染色体上第195015097位；
68.编号63的snp位点位于15号染色体上第25870115位；
69.编号64的snp位点位于2号染色体上第12825813位；
70.编号65的snp位点位于14号染色体上第98924773位；
71.编号66的snp位点位于2号染色体上第120196914位；
72.编号67的snp位点位于6号染色体上第15404754位；
73.编号68的snp位点位于16号染色体上第78321074位；
74.编号69的snp位点位于15号染色体上第50703395位；
75.编号70的snp位点位于17号染色体上第5602353位；
76.编号71的snp位点位于2号染色体上第15690950位；
77.编号72的snp位点位于1号染色体上第280005928位；
78.编号73的snp位点位于6号染色体上第113590851位；
79.编号74的snp位点位于14号染色体上第16009816位；
80.编号75的snp位点位于14号染色体上第153484247位；
81.编号76的snp位点位于14号染色体上第47976495位；
82.编号77的snp位点位于12号染色体上第12928807位；
83.编号78的snp位点位于12号染色体上第31110698位；
84.编号79的snp位点位于5号染色体上第8872919位；
85.编号80的snp位点位于1号染色体上第155807226位；
86.编号81的snp位点位于6号染色体上第46463711位；
87.编号82的snp位点位于9号染色体上第11099863位；
88.编号83的snp位点位于1号染色体上第283667735位；
89.编号84的snp位点位于1号染色体上第17357393位；
90.编号85的snp位点位于13号染色体上第183593435位；
91.编号86的snp位点位于3号染色体上第115303312位；
92.编号87的snp位点位于14号染色体上第43524188位；
93.编号88的snp位点位于8号染色体上第87954260位；
94.编号89的snp位点位于8号染色体上第92363697位；
95.编号90的snp位点位于1号染色体上第69552307位；
96.编号91的snp位点位于14号染色体上第7471774位；
97.编号92的snp位点位于6号染色体上第19952359位；
98.编号93的snp位点位于13号染色体上第21156860位；
99.编号94的snp位点位于13号染色体上第188175696位；
100.编号95的snp位点位于2号染色体上第124815814位；
101.编号96的snp位点位于13号染色体上第14426507位；
102.编号97的snp位点位于11号染色体上第23153305位；
103.编号98的snp位点位于6号染色体上第43151417位；
104.编号99的snp位点位于2号染色体上第98301703位；
105.编号100的snp位点位于9号染色体上第60295990位；
106.编号101的snp位点位于2号染色体上第50604042位；
107.编号102的snp位点位于2号染色体上第152297331位；
108.编号103的snp位点位于14号染色体上第12929923位；
109.编号104的snp位点位于5号染色体上第14421184位；
110.编号105的snp位点位于12号染色体上第36700556位；
111.编号106的snp位点位于6号染色体上第127776310位；
112.编号107的snp位点位于4号染色体上第13581647位；
113.编号108的snp位点位于15号染色体上第31743469位；
114.编号109的snp位点位于8号染色体上第106292625位；
115.编号110的snp位点位于1号染色体上第65086280位；
116.编号111的snp位点位于8号染色体上第100178608位；
117.编号112的snp位点位于6号染色体上第52993684位；
118.编号113的snp位点位于1号染色体上第152542119位；
119.编号114的snp位点位于16号染色体上第74982838位；
120.编号115的snp位点位于6号染色体上第59962724位；
121.编号116的snp位点位于17号染色体上第12798122位；
122.编号117的snp位点位于12号染色体上第9053303位；
123.编号118的snp位点位于6号染色体上第157089555位；
124.编号119的snp位点位于3号染色体上第136625273位；
125.编号120的snp位点位于6号染色体上第38167864位；
126.编号121的snp位点位于15号染色体上第47966193位；
127.编号122的snp位点位于1号染色体上第11314997位；
128.编号123的snp位点位于18号染色体上第32645081位；
129.编号124的snp位点位于11号染色体上第30968353位；
130.编号125的snp位点位于3号染色体上第34875870位；
131.编号126的snp位点位于6号染色体上第10527048位；
132.编号127的snp位点位于1号染色体上第137394988位；
133.编号128的snp位点位于6号染色体上第121919072位；
134.编号129的snp位点位于15号染色体上第40337548位；
135.编号130的snp位点位于13号染色体上第58558625位；
136.编号131的snp位点位于8号染色体上第140609282位；
137.编号132的snp位点位于3号染色体上第79053118位；
138.编号133的snp位点位于1号染色体上第292925911位；
139.编号134的snp位点位于14号染色体上第68031660位；
140.编号135的snp位点位于1号染色体上第112021945位；
141.编号136的snp位点位于8号染色体上第11343109位；
142.编号137的snp位点位于15号染色体上第23522601位；
143.编号138的snp位点位于12号染色体上第47262225位；
144.编号139的snp位点位于6号染色体上第30555475位；
145.编号140的snp位点位于14号染色体上第58652121位；
146.编号141的snp位点位于15号染色体上第13833343位；
147.编号142的snp位点位于4号染色体上第105640665位；
148.编号143的snp位点位于5号染色体上第25590839位；
149.编号144的snp位点位于1号染色体上第28713935位；
150.编号145的snp位点位于7号染色体上第82090814位；
151.编号146的snp位点位于10号染色体上第14181887位；
152.编号147的snp位点位于13号染色体上第199677045位；
153.编号148的snp位点位于1号染色体上第102484718位；
154.编号149的snp位点位于3号染色体上第91245946位；
155.编号150的snp位点位于11号染色体上第77165356位；
156.编号151的snp位点位于12号染色体上第59746957位；
157.编号152的snp位点位于9号染色体上第73545907位；
158.编号153的snp位点位于10号染色体上第77737314位；
159.编号154的snp位点位于11号染色体上第18300706位；
160.编号155的snp位点位于3号染色体上第129063083位；
161.编号156的snp位点位于1号染色体上第42744126位；
162.编号157的snp位点位于12号染色体上第27939849位；
163.编号158的snp位点位于3号染色体上第27041626位；
164.编号159的snp位点位于1号染色体上第58600437位；
165.编号160的snp位点位于1号染色体上第129136174位；
166.编号161的snp位点位于7号染色体上第27082548位；
167.编号162的snp位点位于1号染色体上第173619066位；
168.编号163的snp位点位于1号染色体上第90430196位；
169.编号164的snp位点位于1号染色体上第192669302位；
170.编号165的snp位点位于4号染色体上第124706188位；
171.编号166的snp位点位于13号染色体上第179766527位；
172.编号167的snp位点位于2号染色体上第117677874位；
173.编号168的snp位点位于4号染色体上第36028125位；
174.编号169的snp位点位于3号染色体上第111456973位；
175.编号170的snp位点位于2号染色体上第146937374位；
176.编号171的snp位点位于17号染色体上第34347823位；
177.编号172的snp位点位于6号染色体上第149727715位；
178.编号173的snp位点位于15号染色体上第73085857位；
179.编号174的snp位点位于14号染色体上第39399052位；
180.编号175的snp位点位于8号染色体上第131690844位；
181.编号176的snp位点位于4号染色体上第62131813位；
182.编号177的snp位点位于17号染色体上第62012280位；
183.编号178的snp位点位于2号染色体上第44590539位；
184.编号179的snp位点位于14号染色体上第96480217位；
185.编号180的snp位点位于13号染色体上第49676934位；
186.编号181的snp位点位于9号染色体上第149185476位；
187.编号182的snp位点位于1号染色体上第226497707位；
188.编号183的snp位点位于4号染色体上第94775080位；
189.编号184的snp位点位于11号染色体上第64739915位；
190.编号185的snp位点位于3号染色体上第70024499位；
191.编号186的snp位点位于7号染色体上第122018904位；
192.编号187的snp位点位于2号染色体上第91335279位；
193.编号188的snp位点位于7号染色体上第70704925位；
194.编号189的snp位点位于9号染色体上第55233498位；
195.编号190的snp位点位于9号染色体上第120274182位；
196.编号191的snp位点位于10号染色体上第64259339位；
197.编号192的snp位点位于6号染色体上第109405647位；
198.编号193的snp位点位于1号染色体上第276846341位；
199.编号194的snp位点位于14号染色体上第150391777位；
200.编号195的snp位点位于15号染色体上第126613417位；
201.编号196的snp位点位于5号染色体上第87950981位；
202.编号197的snp位点位于16号染色体上第68067084位；
203.编号198的snp位点位于8号染色体上第83640594位；
204.编号199的snp位点位于13号染色体上第160581841位；
205.编号200的snp位点位于1号染色体上第1614730位。
206.第二方面，本技术还提供一种用于鉴定如第一方面技术方案中所述snp分子标记组合的探针。
207.第三方面，本技术还提供一种如第二方面技术方案中所述探针的试剂盒。
208.第四方面，本技术还提供一种如第一方面技术方案中所述的snp分子标记组合在猪的亲缘关系鉴定中的应用，所述猪的品种包括杜洛克猪、长白猪和大白猪。
209.第五方面，本技术还提供一种对猪群进行snp多态性分析及亲子鉴定推断的方法包括以下步骤：
210.步骤s1，制备基于靶向捕获测序的如第一方面技术方案中所述的snp分子标记组合的液相芯片；
211.步骤s2，利用待测猪基因组dna构建高通量测序文库；
212.步骤s3，将步骤s1制备的猪200snp标记的液相芯片与步骤s2构建的高通量测序文库混合，捕获猪dna高通量测序文库中包含目标snp位点的dna片段；
213.步骤s4，对步骤s3捕获的dna片段进行扩增和纯化，产物经过高通量测序后对测序结果进行参考基因组比对，获得待测猪个体的基因分型；
214.步骤s5，亲子鉴定推断：根据相反纯合子、似然比(与无关个体相比)进行过滤，排除掉大多数不可能的亲子对。
215.优选的，所述猪200snp标记的液相芯片由独立包装的探针和杂交捕获试剂构成；所述探针为单链rna探针，所述杂交捕获试剂为建库试剂盒，包括library prep kit、kit、buffer、binding buffer、wash buffer、post pcr kit。
216.本发明的有益效果：本技术通过提供一种用于猪亲缘关系鉴定的snp分子标记组合，并建立靶向捕获测序方法进行检测，来保证育种中系谱的完整和准确，进而推动猪产业的遗传育种改良的发展，经过靶向捕获测序检测snp多态性，经高深度测序，分型结果准确可靠，保证了亲子关系的精准性。
附图说明
217.此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。
218.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
219.图1为靶向捕获测序技术的原理示意图。
具体实施方式
220.下面结合说明书附图和具体实施方式，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本技术所附权利要求书所限定的范围。
221.本发明通过对3923个杜洛克个体、4440个长白个体和8620个大白个体的中芯一号芯片数据分析，筛选出18条常染色体上的200个snp位点(参考基因组sscrofa10.1)，可用于杜、长、大商业化品种猪亲缘关系鉴定的snp标记组合，其包含200个snp标记的分型及物理位置。
222.参见表1，本发明提供的用于杜、长、大商业化品种猪亲缘关系鉴定的snp分子标记组合由以下200个snp标记组成。
223.表1
224.225.226.227.228.229.230.231.232.[0233][0234]
上述200个snp标记在杜洛克、长白和大白群体中的基因频率与基因型频率如下表2所示：
[0235]
表2
[0236]
[0237]
[0238]
[0239]
[0240]
[0241]
[0242]
[0243]
[0244]
[0245][0246]
表2中，p、h、r分别代表纯合子aa基因型频率、杂合子ab基因型频率和纯合子bb基因型频率、maf代表最小等位基因频率。
[0247]
上述200个snp标记在靶向捕获测序流程中的探针序列如下表3所示：
[0248]
表3
[0249]
[0250]
[0251]
[0252]
[0253]
[0254]
[0255]
[0256]
[0257]
[0258][0259]
本发明经过靶向捕获测序检测snp多态性以及亲子推断试验验证，上述200个snp标记的按照相对严格的质控参数即call rate(》0.95)、哈温平衡(p》10-6
)、maf》0.3、ld过滤(r2》0.5)筛选出，用于确定亲子关系的置信度极高。并且实例验证用该发明的200个snp标记推断的亲子关系与实际亲本关系完全一致(实际亲本已经过猪群企业的相关记录验证是正确的)。
[0260]
本发明经过靶向捕获测序检测snp多态性(靶向捕获测序技术的原理图参见图1，包括wgs文库变性、加入捕获探针、探针与目标片段杂交、链霉亲和素磁珠结合、磁力架吸附、pcr富集目标片段以及上机测序)，经高深度测序，分型结果准确可靠，保证了亲子关系的精准性。
[0261]
靶向捕获测序技术的原理
[0262]
下面结对本发明的具体实施方式作进一步详细描述。
[0263]
实施例1：snp标记组合的选择；
[0264]
本发明对用于杜、长、大商业化品种猪群进行亲子鉴定体系的snp标记选择经过了3个过程，具体如下：
[0265]
(1)群体划分
[0266]
将杜、长、大商业化品种猪群体划分为标记筛点群体和标记验证群体，筛点群体用于筛点和验证，标记验证群体只用于标记验证。
[0267]
(2)标记筛选
[0268]
对不同群体中芯一号的基因型数据，首先按照callrate(》0.95)，哈温平衡(p》10-6
)进行过滤，剔除不符合要求的snp标记，然后筛选出maf》0.3的snp标记。最终找出所有群体符合要求的snp标记的交集。最后对这些标记进行ld过滤(r2》0.5，剔除物理位置极近的snp标记)。具体如下：
[0269]
第一步，按照callrate》0.95，哈温平衡(p》10-6
)，maf》0.3这三个条件，对1419个杜洛克个体、1472个长白个体和1052个大白个体的基因型数据进行过滤，保留了1012个三个群体的共同snp标记集合。
[0270]
第二步，采用2504个杜洛克个体、2968个长白个体和7568个大白个体的中芯一号的基因型数据，计算这1012个snp标记的哈温平衡的p值和maf，剔除哈温平衡p值《10-6
，maf《0.3的位点，保留327个snp标记。
[0271]
第三步，采用2504个杜洛克个体、2968个长白个体和7568个大白个体的中芯一号的基因型数据对这327个位点进行ld过滤(r2《0.5)，剩下202个位点。
[0272]
第四步，根据位点的物理位置，剔除了2个相距极近的标记，最终筛选出200个snp
标记。
[0273]
(3)标记验证
[0274]
采用筛选出来的200个snp标记，利用3923个杜洛克个体、4440个长白个体和8620个大白个体的基因型数据，利用sequoia软件进行系谱重建，与中芯一号芯片数据的亲子鉴定结果(真值)比对，计算错误率(eror rate)，计算公式为：
[0275]
error rate＝分配了错误亲本的个体数/总的个体数；
[0276]
计算结果参见表4。
[0277][0278][0279]
由表4可以看出，200个snp标记在杜、长、大商业化品种猪群体中进行亲子鉴定的准确率达到99.999％以上。
[0280]
实施例2：采样靶向捕获测序的技术对猪群进行snp多态性分析及亲子鉴定推断；具体包括以下爱步骤：
[0281]
步骤s1，制备基于靶向捕获测序的猪200snp标记的液相芯片；
[0282]
在此步骤中，所述猪200snp标记的液相芯片由独立包装的探针和杂交捕获试剂构成。所述猪200snp标记探针为单链rna探针，是根据筛选到的snp位点设计并合成的核苷酸序列。探针的设计原则为探针长度120bp,探针gc含量25％-75％之间、同源性比对个数≤10。根据筛选获得的snp位点设计2条覆盖所述snp位点的寡核苷酸序列，核苷酸序列5’短带有生物素基团修饰，即所述的液相芯片的探针。将合成的2条rna探针进行等摩尔质量混合，利用edta和tris-hcl混合液定容为3pmol/ml的探针混合液。
[0283]
本实施例中，所述杂交捕获试剂为康普森自主研发的建库试剂盒，包括library prep kit、kit、buffer、binding buffer、wash buffer、post pcr kit。
[0284]
步骤s2，利用待测猪基因组dna构建高通量测序文库；步骤s2具体包括以下步骤：
[0285]
步骤s2.1，提取待测猪只的基因组dna：
[0286]
本实施例选取包含父母代和子代的杜洛克、长白、大白三个猪品种群体进行检测，其中杜洛克194个个体，长白161个个体、大白220个个体，共计575个血液样本进行基因组dna的提取，具体方法如下：
[0287]
(1)向1.5ml的离心管中加入20ul蛋白酶k，之后加入200ul的血液。
[0288]
注意：
[0289]
1)冷冻抗凝血液需提前在室温(15～30℃)下放置，融化混匀。
[0290]
2)如血液体积大于或小于200ul，蛋白酶k、buffer ml和异丙醇与磁珠混合物的用量需按比例调整。
[0291]
(2)向离心管中加入200ul buffer ml，涡旋振荡5秒钟使其充分混匀后，将离心管放于56℃水浴锅中孵育15分钟，期间涡旋震荡混匀2次。
[0292]
(3)将离心管从水浴锅中取出，短暂离心后室温放置5分钟。加入彻底混匀的320ul异丙醇与magbeads混合物，涡旋震荡混匀5秒钟后将离心管放于25℃、1600rpm的恒温混匀
[0293]
仪上震荡混匀5分钟或将离心管连续颠倒混匀10分钟。
[0294]
(4)将离心管放于磁力架上静置1分钟，待magbeads完全吸附于离心管侧壁后彻底弃去溶液(保持离心管固定于磁力架上)。
[0295]
(5)将离心管从磁力架上取下，加入750ul buffer gw1(使用前请检查是否已加入无水乙醇)后涡旋点震1分钟或涡旋振荡5秒钟后放于25℃、1600rpm的恒温混匀仪上震荡混匀2分钟(震荡过程中确保magbeads处于混匀状态)。之后将离心管放于磁力架上静置1分钟，待magbeads完全吸附于离心管侧壁后轻轻颠倒磁力架将离心管盖上的杂质洗落后彻底弃去溶液(保持离心管固定于磁力架上)。
[0296]
(6)重复步骤5。
[0297]
(7)将离心管从磁力架上取下，加入750ul buffer gw2(使用前请检查是否已加入无水乙醇)后涡旋点震1分钟或涡旋振荡5秒钟后放于25℃、1600rpm的恒温混匀仪上震荡混匀2分钟(震荡过程中确保magbeads处于混匀状态)。之后将离心管放于磁力架上静置1分钟，待magbeads完全吸附于离心管侧壁后轻轻颠倒磁力架将离心管盖上的杂质洗落后彻底弃去溶液(保持离心管固定于磁力架上)。
[0298]
(8)重复步骤7。
[0299]
(9)保持离心管固定于磁力架上，用移液器进一步去除离心管管底和管盖上的溶液，之后室温放置5-10分钟，使乙醇挥发干净。
[0300]
注意：如果离心管侧壁上有液珠，可向离心管中加入750ul无水乙醇。盖盖后颠倒离心管(保持离心管固定于磁力架上)，之后彻底弃去无水乙醇。
[0301]
(10)将离心管从磁力架上取下，加入50-200ul buffer eb。涡旋震荡使磁珠完全悬浮于洗脱液中后将其放于56℃、1600rpm的恒温混匀仪上震荡洗脱10分钟，或将离心管放于56℃水浴锅中孵育10分钟，期间每隔3分钟涡旋震荡10秒钟。
[0302]
(11)将离心管放于磁力架上静置2分钟，待magbeads完全吸附于离心管侧壁后用移液器将洗脱液转移至新的离心管中-20℃保存备用。
[0303]
步骤s2.2，构建猪dna高通量测序文库：
[0304]
(1)基因组dna片段化、末端修复、3’加a，具体包括：
[0305]
在pcr管中配置总体积为30ul的反应体系，取步骤s2.1获得的猪dna50-200 ng、enzyme fragmentation mix 5ul、余量为nuclease-free water。将反应体系轻柔混匀后短暂离心将反应液收集至管底，置于pcr仪中，进行如表5所示的反应，85℃热盖。
[0306]
表5
[0307]
[0308][0309]
(2)接头连接
[0310]
将上述反应管置于冰盒上，加入adapter for mgi 5ul、5
×
ligation buffer 20ul、ligase 10ul、nuclease-free water 15ul，轻轻吹吸混匀并短暂离心将反应液收集至管底，置于pcr仪中，进行如表6所示的反应，取消热盖。
[0311]
表6
[0312]
温度(℃)时间(min)20℃15min4℃保持
[0313]
(3)连接体系纯化
[0314]
1)在上述反应体系中加入80ul的agencourtampure xp磁珠，轻柔吹打混匀，室温静置5-10min；
[0315]
2)将pcr管置于磁力架上3min使溶液澄清，移除上清；
[0316]
3)pcr管继续放置在磁力架上，向pcr管内加入200ul80％乙醇溶液，静置30s；
[0317]
4)移除上清，再向pcr管内加入200ul80％乙醇溶液，静置30s后彻底移除上清；
[0318]
5)室温静置3-5min，使残留乙醇彻底挥发；
[0319]
6)加入22ul的nuclease-free water，把pcr管从磁力架取下，轻轻吸打重悬磁珠，避免产生气泡，室温静置2min；
[0320]
7)将pcr管置于磁力架上2min使溶液澄清，用移液器吸取20ul上清液，转移到新的pcr管中(置于冰盒上)。
[0321]
(4)杂交前pcr反应：
[0322]
将上述得到的纯化产物中加入enzyme hifi pcr mix 25ul、barcode-1(20μm)2.5ul、barcode-2(20μm)2.5ul，涡旋或吹吸混匀并短暂离心将反应液收集至管底，置于pcr仪中，进行如表7所示的反应，热盖温度105℃。
[0323]
表7
[0324][0325]
(5)pcr产物纯化：
[0326]
1)在上述反应体系中加入50ul的agencourtampure xp磁珠，轻柔吹打混匀，室温静置5-10min；
[0327]
2)将pcr管置于磁力架上3min使溶液澄清，移除上清；
[0328]
3)pcr管继续放置在磁力架上，向pcr管内加入200ul80％乙醇溶液，静置30s；
[0329]
4)移除上清，再向pcr管内加入200ul80％乙醇溶液，静置30s后彻底移除上清；
[0330]
5)室温静置3min-5min，使残留乙醇彻底挥发；
[0331]
6)加入30ul的nuclease-free water，把pcr管从磁力架取下，轻轻吸打重悬磁珠，避免产生气泡，室温静置2min；
[0332]
7)将pcr管置于磁力架上2min使溶液澄清，用移液器吸取28ul上清液，转移到新的pcr管中(置于冰盒上)。
[0333]
(6)文库质检
[0334]
取1ul样品使用3.0fluorometer(qubit dsdna hs assay kit)进行文库浓度测定，记录文库浓度，如下游进行液相捕获，则需文库浓度》25ng/ul；取3ul样品进行琼脂糖凝胶电泳检测片段长度，文库长度约在300bp-500bp间，合格后进行后续实验步骤，或将文库置于-20℃冰箱中保存。
[0335]
步骤s3，将步骤s1制备的猪200snp标记液相芯片与步骤s2构建的高通量测序文库混合，捕获猪dna高通量测序文库中包含目标snp位点的dna片段。
[0336]
(1)文库制备和浓缩
[0337]
1)预先从﹣20℃保存的捕获试剂盒中，取出下表所示block试剂，放置于冰盒上溶解，溶解混匀后置于冰上或4℃待用；
[0338]
2)将样品文库与相对应的block按照下表中相对应的反应体系配制，混匀后标记为b管。(“对应”指文库index类型与block类型的对应，eg：使用了单端8bp的index类型，则参照表8所示的体系进行配制)。
[0339]
表8
[0340]
组分体积样品文库750nghyb block5μladapter-block6μl
[0341]
3)将b管放入真空浓缩离心机，打开pcr管盖，启动离心机，打开真空泵开关，开始浓缩(浓缩之前可以先用相同体积的水进行浓缩时间测试，计算样品单位体积减少的时间，确认浓缩速率，避免因浓缩时间过长导致过度干燥，造成样品损耗)；
[0342]
4)将b管浓缩至体积小于10ul，之后用nuclease-free water补足至10ul，轻轻吸打混匀，短暂离心后置于冰上待用。
[0343]
(2)探针杂交
[0344]
1)将buffer置于室温融化，混匀后置于50℃恒温混匀仪内预热，完全溶解后(无沉淀及浑浊物)每个反应取20ulbuffer置于新的200ul pcr管内，盖好管盖，置于50℃恒温混匀仪孵育待用，标记为a管；
[0345]
2)取5ul rnase block与2ul probe置于200ulpcr管内，轻轻吸打混匀，短暂离心后置于冰上或4℃待用，标记为c管；
[0346]
3)从a管吸取13ul与b管10ul加入c管中震荡混匀，短暂离心后置于冰上待用；
[0347]
4)将pcr管c置于pcr仪上，运行如表9所示的程序，热盖85℃：
[0348]
表9
[0349]
温度(℃)时间(min)80℃5min50℃hold
[0350]
(3)捕获磁珠平衡
[0351]
1)捕获磁珠(cap beads)从4℃取出，涡旋震荡重悬；
[0352]
2)取50ul磁珠置于新的pcr管内，置于磁力架上1min使溶液澄清，移除上清；
[0353]
3)从磁力架上取下pcr管，加入200ulbinding buffer轻轻吸打数次混匀，重悬磁珠；
[0354]
4)置磁力架上1min，移除上清；
[0355]
5)重复步骤3以及步骤4两次，共清洗磁珠3次；
[0356]
6)从磁力架上取下pcr管，加入180ulbinding buffer轻轻吸打6次重悬磁珠待用。
[0357]
(4)捕获目标区域dna文库
[0358]
1)保持杂交产物在pcr仪上，将上述步骤重悬后180μl cap beads加入到杂交产物中，用移液器吸打6次混匀，置于旋转混匀仪上室温结合30min；
[0359]
2)将pcr管置于磁力架上2min使溶液澄清，移除上清液；
[0360]
3)向杂交产物内加入200ul的wash buffer 1，轻轻吸打6次混匀，置于旋转混匀仪上清洗15min，然后短暂离心，将pcr管放于磁力架上2min，使溶液澄清，移除上清；
[0361]
4)加入200ul的50℃预热的wash buffer 2，轻轻吸打6次混匀，置于恒温振荡混匀仪上50℃孵育10min；
[0362]
5)短暂离心，将pcr管放于磁力架上2min，移除上清。使用buffer 2清洗2次，共计3次。最后一次彻底移除buffer 2(可用10ul移液器移除残留)；
[0363]
6)保持样品在磁力架上，向pcr管内加入200ul80％乙醇，静置30s后彻底移除乙醇溶液(可用10ul移液器移除残留)，室温晾干；
[0364]
7)向pcr管加入30ul nuclease-free water，从磁力架上取下pcr管，轻轻吸打6次重悬磁珠待用。
[0365]
步骤s4，对步骤s3捕获的dna片段进行扩增和纯化，产物经过高通量测序后对测序结果进行参考基因组比对，获得待测猪个体的基因分型。
[0366]
(1)杂交后pcr反应
[0367]
1)预先从-20℃保存的试剂盒中取出post pcr buffer、post pcr primer(25μm)试剂，放置冰上溶解，溶解混匀后置于冰上待用；取出post pcr polymerase简短离心后立即使用。
[0368]
2)捕获后需要对dna文库进行富集，根据下表配制如表10所述的反应体系：
[0369]
表10
[0370]
组分体积来自步骤s3反应结束的样品30μlpost pcr buffer18μl
post pcr primer(25μm)1μlpost pcr polymerase1μltotal volume50μl
[0371]
3)将移液器调至40ul，轻轻吸打混匀6次，然后立即置于pcr仪上，运行如表11所示的程序，热盖105℃。
[0372]
表11
[0373][0374]
(2)文库纯化
[0375]
1)pcr结束后向样品加入55ulagencourtampure xp磁珠，用移液器轻轻吸打6次混匀；
[0376]
2)室温孵育5min，把pcr管置于磁力架上3min使溶液澄清；
[0377]
3)移除上清，pcr管继续置于磁力架上，加入200ul的80％乙醇，静置30s；
[0378]
4)移除上清，再向pcr管内加入200ul的80％乙醇，静30s后彻底移除上清(可用10ul移液器移除底部残留酒精)；
[0379]
5)室温放置5min，使得残留乙醇彻底挥发；
[0380]
6)加入25ul nuclease-free water，将pcr管从磁力架拿下，轻轻吹打混匀重悬磁珠，室温放置2min；
[0381]
7)将pcr管置于磁力架上2min；
[0382]
8)用移液器吸23ul上清液转移到1.5ml离心管，标记样品信息；
[0383]
9)取1ul文库使用qubit dsdna hs assay kit进行定量，记录文库浓度，文库浓度约在1ng/ul-20ng/ul；
[0384]
10)取1ul样品使用agilent 2100bioanalyzer system(agilent dna1000kit)进行片段长度测定，文库长度约在300bp-500bp间；
[0385]
11)利用高通量测序平台进行测序，并进行参考基因组比对，获得待测猪个体的基因分型。
[0386]
步骤s5，亲子鉴定推断
[0387]
第一步：根据相反纯合子、似然比(与无关个体相比)进行过滤，排除掉大多数不可能的亲子对；
[0388]
第二步：对于每个个体(按照最早出生-最晚出生-出生日期未知的顺序)，所有前面满足过滤步骤并且出生年份更早或者未知的个体视为候选亲本，然后计算该个体与候选亲本亲子关系和全同胞关系的似然值。如果亲子关系的似然值更高，且亲子关系与全同胞关系的似然比大于给定阈值，就实现了一个亲子分配。如果存在多个符合条件的候选亲本，则分配似然比最高的候选亲本。
[0389]
该实施案例的亲子鉴定的结果文件如下：
[0390]
194个杜洛克个体的亲子关系鉴定结果参见表12，161个长白个体亲子关系鉴定结果参见表13，220个大白个体亲子关系鉴定结果参见表14。
[0391]
表12
[0392]
[0393]
[0394]
[0395]
[0396][0397]
表13
[0398]
[0399]
[0400]
[0401]
[0402]
[0403][0404]
表14
[0405]
[0406]
[0407]
[0408]
[0409]
[0410]
[0411][0412]
以上鉴定结果与系谱关系一致率100％，说明该发明专利的200个snp标记能够准确的进行杜、长、大商业化品种猪亲缘关系鉴定。
[0413]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。