本发明涉及猪多组学整合精准育种方法，具体涉及用以提高猪的繁殖效率的多组学整合精准育种方法。
背景技术：
我国养猪业历史悠久，猪种资源丰富，分布广泛，有着巨大的利用潜力和广阔的发展前景。种猪生产是养猪业生产的基础，只有通过科学地饲养管理，种猪才能繁殖出数量多、质量高的仔猪，才能提高养猪业的经济效益。
随着我国经济持续快速的发展，人们的生活水平逐渐提高，对猪肉的需求量也越来越高，因此，人们也越来越关注如何提高猪的繁殖性能。繁殖性能是猪的重要经济性状，包括产仔数、初生重和初生窝重、断乳窝重、初产日龄和产仔间隔等。母猪繁殖力的选育对于提高养猪生产整体效益起到至关重要作用，给养猪业带来巨大经济效益，是影响养猪业生产效率的重要因素之一。
一方面，目前获得高繁殖效率猪的方法主要集中在基因组编辑。在转基因技术和体细胞克隆技术的基础上发展起来的基因组编辑技术，可以通过对目标基因进行“编辑”，包括对特定dna片段的敲除、特异突变引入和定点dna片段转入等，从而实现对基因组的精确修饰。我国在基因组编辑猪相关方面研究进展迅速，并获得了一系列农用、医用工程猪。科研人员通过基因操作将有利基因引入特定的位点，从而提高外源基因的稳定性和效率，改变动物的遗传特定，如提高动物的生产性能，最终育成满足人们需要的高产、优质、抗病动物新品种。此外，也可以对于一些负调控动物生长的基因进行修饰或敲除，培育出高产的转基因动物。
然而，人们在对猪肉产量要求提高的同时，越来越注重猪的培育方法对猪品质的作用，对人类健康的影响。由于科学界难以在短时间内证实转基因动物对人类的健康完全无副作用，以及广大群众对转基因动物的认知不足，人们对食用非转基因动物存在强烈的渴求，同时造成了非转基因食用动物价格高于转基因动物。
另一方面，除去基因组编辑技术，其他育种理论的实现同样需要对猪基因组的了解。然而，目前对猪基因组的了解和认识远远不足，对影响猪繁殖效率(包括猪产仔数和平均胎重等)的基因掌握尚不全面，有必要进一步进行研究，满足猪培育方法的需求。
因此，亟需建立一种全新的高繁殖效率猪的培育方法，以满足市场需求，提高经济效益，实现我国猪育种理论和技术水平的提升，实现我国猪种业的跨越式发展。
技术实现要素：
本发明人进行了锐意研究，针对上述猪育种方法中存在的问题，通过在转录水平上对候选基因的定量，确定了hsd17b1基因对猪产仔性状的影响；同时通过全基因组关联分析的方法，丰富了影响猪繁殖效率的snp位点；将转录水平和基因水平上影响特征相结合，提供一种用以提高猪的繁殖效率的多组学整合精准育种方法，该方法未采用基因组编辑手段，提高广大群众接受度，同时克服经典遗传育种理论只考虑dna序列变异体对性状遗传变异影响的弊端，综合多层次遗传变异效应，从而完成本发明。
本发明的目的在于提供以下方面：
(1)猪多组学整合精准育种方法，其中，所述方法包括测定hsd17b1基因表达量、以及hsd17b1基因上snp位点基因型的步骤。
在一种优选的实施方式中，通过测定猪hsd17b1基因mrna的表达量，判定猪产仔性状：hsd17b1基因mrna的表达量越低，猪具有更高的产仔数和/或窝重。
在一种优选的实施方式中，所述育种方法中通过获得猪卵巢组织或其他组织细胞的总rna，经反转录、定量pcr确定猪hsd17b1基因mrna的表达量；
优选地，定量pcr步骤中，包括用于扩增hsd17b1基因mrna反转录产物的扩增引物p3和p4，所述扩增引物p3和p4的核苷酸序列分别包括seqidno.2和seqidno.3所示序列。
在一种优选的实施方式中，所述snp位点位于猪hsd17b1基因第323407415位，该snp位点的两种等位基因为a或g；
snp位点对应的三种基因型分别为aa、gg和ag，aa基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为a的纯合体；gg基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为g的纯合体；ag基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为a与g的杂合体；
通过测定猪hsd17b1基因第323407415位snp位点对应的三种基因型，判定猪产仔性状：基因型为ag和gg的猪相较于基因型为aa的猪有更高的产仔数和/或窝重。
在一种优选的实施方式中，用以检测所述snp位点的方法中包括一组扩增引物p1和p2、以及延伸引物，
所述扩增引物p1和p2的核苷酸序列分别包括seqidno.4和seqidno.5所示序列。
在一种优选的实施方式中，所述延伸引物的核苷酸序列分别包括seqidno.6所示序列。
在一种优选的实施方式中，snp位点处基因型的测定方法，所述测定方法包括：
方法1，检测猪hsd17b1基因第323407415位脱氧核糖核苷酸为a或g、或a和g，以确定待测猪的基因型是aa、gg或ag；和/或
方法2，通过sequenommassarray方法直接确定基因型。
在一种优选的实施方式中，所述育种方法包括以下步骤：
步骤1，选择hsd17b1基因mrna低表达量、且hsd17b1基因第323407415位为gg或ag基因型的猪作为亲本杂交，获得f1代仔猪；
步骤2，对f1代仔猪进行hsd17b1基因mrna定量分析、以及第323407415位基因型确定，选育hsd17b1基因mrna的低表达量、且第323407415位为gg或ag基因型的仔猪，成年后继续杂交得f2代仔猪；
步骤3，对f2代仔猪进行hsd17b1基因mrna定量分析、以及第323407415位基因型确定，选育hsd17b1基因mrna的低表达量、且第323407415位为gg或ag基因型的仔猪，成年后继续杂交得f3代仔猪。
重复上述步骤2或3，最终得到hsd17b1低表达量、基因型gg的高产仔数和/或窝重的种猪。
在一种优选的实施方式中，所述育种猪为大白猪和/或梅山猪。
(2)猪多组学整合精准育种方法在培育/鉴定高产仔数和窝重的种猪方面的用途。
根据本发明提供的猪多组学整合精准育种方法，具有以下有益效果：
1、本发明提供的基于hsd17b1基因和/或snp位点的育种方法，在自然杂交的基础上进行人工选择，不依赖于基因组编辑，可更好满足市场对非转基因肉类的需求；同时将转录层面和基因层面与猪产仔性状相关的因素相结合，实现多组学整合精准育种；
2、本发明提供了用于检测snp位点和基因分型的试剂盒，其使用方便，操作简单，成本低；
3、本发明提出的多组学整合精准育种方法为我国猪种业的发展提供了新的思路，可加快各系列种猪的遗传进展。
附图说明
图1是实施例1中梅山猪与大白猪hsd17b1基因表达量结果图。
具体实施方式
下面通过对本发明进行详细说明，本发明的特点和优点将随着这些说明而变得更为清楚、明确。在这里专用的词“示例性”意为“用作例子、实施例或说明性”。这里作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。
本发明人长期致力于种猪培育工作，经过大量研究和实践发现：dna水平变异对猪群体和个体表型差异的作用被高估，基因组表达调控机制不完全清楚；多层次组学对表型的效应被低估，猪种个性化遗传机制未受到关注。猪品种间、个体间巨大的表型遗传差异不仅取决其基因组结构和dna顺序的变异，更主要取决于基因组多层次表达调控。为了实现基因组信息的高效应用，有必要开展猪多组学变异研究，整合多组学信息、建立种群个性化多组学整合精准育种理论和技术体系，实现多组学选择。实现多组学整合精准育种理论创新和应用，将显著提升我国猪育种理论和技术水平。
基于上述考量，本发明在多层次如转录水平上和基因水平上对影响猪生产效率的因素进行了分析和验证。
为克服克服现有技术的缺陷，本发明一方面，提供了与猪生产效率即产仔性状相关的基因，所述基因为hsd17b1基因。
hsd17b1基因是hsd17b家族的一个成员，主要在雌激素生物合成的最后一步发挥作用，在动物生长、发育和生殖过程中具有重要的调控作用，能够将雌激素转化为更有生物活性的雌二醇。目前对hsd17b1与一些生殖疾病的发生、发展进行了关联，但尚未有明确的结论。有研究者在乳腺癌、子宫内膜癌中都检测到hsd17b1基因的表达；不同的研究者还分别对hsd17b1基因的多态性或者变异与乳腺癌、子宫内膜癌的相关性进行研究，结果表明，hsd17b1基因的多态性或者变异都与研究的癌症无显著性相关。
同时了解到，猪生产效率指标如产仔数、胎重和窝重等与猪的子宫内膜或乳腺的生长发育相关，进而可能与hsd17b1基因的表达相关：猪妊娠期，绒毛膜上皮和子宫内膜上皮形成皱褶，皱褶折叠水平影响着胎盘的相对大小，胎盘的相对大小是影响胎儿存活的主要因素之一，胎儿的存活率又与猪初生重具有重要的关系，初生重大可以刺激母猪多产奶，而hsd17b1基因在乳腺癌、子宫内膜癌中均有表达，因此，本发明人认为与猪产奶、胎盘大小有间接关系的hsd17b1基因为影响猪生产效率的产仔性状候选基因。然而目前，尚未有hsd17b1基因与猪产仔性状关联的报道。
本发明人经过深入研究，发现了猪hsd17b1基因与生产效率即产仔性状的相关性，所述hsd17b1基因具有seqidno.1所示的核苷酸序列；通过测定猪hsd17b1基因mrna的表达量，判定猪产仔性状；hsd17b1基因mrna的表达量越低，猪具有更高的产仔数和/或窝重。
本发明另一方面，提供了一个与猪生产效率即产仔性状相关的snp位点。优选地，选择大白猪或梅山猪进行产仔性状相关的基因和snp位点的鉴定。大白猪和梅山猪在我国的猪种培育中占较高比重，选择上述两种猪，便于后续对我国猪种的培育。
本发明运用全基因组关联分析(gwas)的方法寻找与猪产仔性状相关的snp位点。gwas是利用全基因组范围内筛选出高密度的分子标记对所研究的群体进行扫描，分析扫描得出的分子标记数据与表型性状之间关联关系的方法，该方法也是目前进行相关snp位点筛选最为广泛使用的方法。然而，运用全基因组关联分析的方法寻找与猪产仔性状相关的snp位点，得到的结果庞杂无序，大多数的基因变异与目标性状并不相关。在此庞杂无序的结果中选择相关snp位点是相当困难的。因而，有必要预先选择特定的基因，结合gwas方法确定位于特定基因上的相关snp位点。
本发明在众多筛选的基因中选择特定的基因，进而确定该基因上snp位点，极大提高了对于相关snp位点确定的成功率。基于对hsd17b1基因作为与猪产仔性状相关基因的考量，所述特定的基因是hsd17b1基因。
本发明人经过深入研究，首次在猪hsd17b1基因中确定了产仔性状snp位点。所述snp位点位于hsd17b1基因第323407415位(rs323407415)，该snp位点的两种等位基因为a或g。优选地，核苷酸a的基因频率为0.159，核苷酸g的基因频率为0.841。
该snp位点对应的三种基因型分别为aa、gg和ag，aa基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为a的纯合体；gg基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为g的纯合体；ag基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为a与g的杂合体。经过大量研究发现，基因型为ag和gg的猪相较于基因型为aa的猪有更高的产仔数和窝重。
基于多组学整合精准育种的考量，本发明提供了一种多组学整合精准育种方法，所述方法包括测定hsd17b1基因表达量、以及hsd17b1基因第323407415位snp位点上基因型的步骤。
在一种优选的实施方式中，测定猪hsd17b1基因mrna的表达量的方法包括：获得猪卵巢组织或其他组织细胞的总rna，经反转录、定量pcr确定猪hsd17b1基因mrna的表达量。mrna的表达量越低，猪有更高的产仔数和/或窝重。
其中，定量pcr扩增程序为：预变性95℃，5min；变性95℃，5s；延伸60℃，30s；重复40个循环；95℃，15s；60℃，1min；95℃，15s；60℃，15s。
在进一步优选的实施方式中，本发明提供了一种用于测定hsd17b1基因mrna的表达量的pcr扩增引物p3和p4，p3和p4的碱基序列分别如seqidno.2和seqidno.3所示。
在进一步优选的实施方式中，本发明提供了一种用于测定hsd17b1基因mrna的表达量的试剂盒，所述试剂盒包括扩增引物p3和p4，p3和p4的碱基序列分别如seqidno.2和seqidno.3所示。
在一种优选的实施方式中，本发明中利用sequenommassarray检测系统确定该snp位点的基因分型。sequenommassarray检测系统的基本原理为基质辅助激光解吸电离飞行时间质谱(maldi-tofms)技术，首选pcr扩增目标序列，然后加入snp序列特异延伸引物，在snp位点上，延伸1个碱基。将制备的样品分析物与芯片基质共结晶后在质谱仪的真空管经瞬时纳秒(10-9s)强激光激发，核酸分子解吸附并转变为亚稳态离子(多为单电荷离子)，电厂中离子飞行时间与离子质量成反比，通过飞行时间检测器检测核酸分子在真空管中的飞行时间而获得样品分析物的精确分子量，从而检测出snp位点信息。
由上述知，snp位点基因分型测定过程中，涉及猪基因组的dna的提取，pcr扩增反应和单碱基延伸反应等步骤。
在进一步优选的实施方式中，所述pcr扩增反应中扩增引物p1和p2的核苷酸序列分别如seqidno.4和seqidno.5所示。
在进一步优选的实施方式中，所述单碱基延伸反应中延伸引物的核苷酸序列如seqidno.6所示。
在另一种优选的实施方式中，提供了一种与猪产仔性状相关的snp位点的检测方法，所述检测方法包括直接测序或者通过试剂盒测定。
所述试剂盒包括(1)snp位点的pcr延伸引物；优选还包括(2)pcr反应试剂：包括pcr缓冲液、dntp、pfudna聚合酶，snapshot混合液；(3)pcr扩增引物p1和p2；(4)pcr产物纯化试剂：包括核酸外切酶、虾-碱性磷酸酶，以及纯化用缓冲液；(5)snp位点对照标本：包括纯合突变阳性对照标本、杂合突变阳性对照标本、阴性对照标本。
所述snp纯合突变阳性对照标本为含有猪rs323407415位点人工合成为g碱基的质粒。
所述snp杂合突变阳性对照标本为含有猪rs323407415位点人工合成a和g碱基的两种质粒的混合物。
所述snp阴性对照标本为人工合成的带有随机碱基的质粒。使用三种对照标本对反应体系进行质量控制，可保证检出结果的一致性。
本发明中提供的检测试剂盒的使用方法为：从猪耳或其它组织细胞中提取dna样品；配制pcr反应体系，进行pcr扩增，纯化pcr产物；pcr产物与snp位点的延伸引物同时进行延伸反应；对延伸产物进行毛细管电泳分析；snp位点分析，进而还可得到基因分型信息。上述snp位点分析可采用但不限于abi3730xl自动测序仪。
在另一种优选的实施方式中，还提供了一种与猪产仔性状相关的snp位点基因分型的检测方法，该方法包括以下步骤：
(1)提取猪的基因组dna；
(2)以待测猪的基因组dna为模板，利用扩增引物和延伸引物进行sequenommassarray检测，确定猪hsd17b1基因第323407415位的基因型。
其中，所述pcr扩增反应中扩增引物(p1和p2)如seqidno.4和seqidno.5所示；延伸引物如seqidno.6所示。
本发明中，具体地，多组学整合精准育种方法包括以下步骤：
步骤1，选择hsd17b1基因mrna低表达量、且hsd17b1基因第323407415位为gg或ag基因型的猪作为亲本(f0)杂交，获得第一代(f1代)仔猪；
步骤2，对f1代仔猪进行hsd17b1基因mrna定量分析、以及第323407415位基因型确定，选育hsd17b1基因mrna的低表达量、且第323407415位为gg或ag基因型的仔猪，成年后继续杂交得f2代仔猪；
步骤3，对f2代仔猪进行hsd17b1基因mrna定量分析、以及第323407415位基因型确定，选育hsd17b1基因mrna的低表达量、且第323407415位为gg或ag基因型的仔猪，成年后继续杂交得f3代仔猪。
重复上述步骤2或3，最终得到hsd17b1低表达量、基因型gg的高产仔数和/或窝重的种猪。
上述方法遵从自然繁殖，通过较少人工干预，获得高产仔数/窝重的种猪；同时，所述育种方法为非转基因方法，群众的接受度更高，可有效促进养猪业发展，提升我国猪种业的核心竞争力。
本发明中，还提供了一种鉴定或辅助鉴定猪产仔性状的方法，至少包括以下两步骤：
步骤(1)，测定猪hsd17b1基因mrna的表达量，mrna的表达量低的猪，具有较高的产仔数和/或窝重；和
步骤(2)，检测猪hsd17b1基因第323407415位脱氧核糖核苷酸为a、或g、或a和g，以确定待测猪的基因型是aa、gg或ag，或者直接通过sequenommassarray等方法确定基因型，根据基因型确定猪产仔数和/或窝重：基因型为ag和gg的猪相较于基因型为aa的猪有更高的产仔数和/或窝重，
其中，aa基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为a的纯合体；gg基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为g的纯合体；ag基因型为猪hsd17b1基因第323407415位脱氧核糖核苷酸为a与g的杂合体。
本发明中，根据上述hsd17b1基因、pcr定量扩增引物、表达量测试方法、表达量检测试剂盒、snp位点、snp位点延伸引物、基因扩增引物、snp位点测试方法、snp位点检测试剂盒、snp位点基因分型，还可提供基于上述公开在鉴定和/或选育猪种方面的应用，所述猪种具有较高的生产效率，所述生产效率包括产仔数和窝重。
结合上述基因层面的分析和转录层面的分析，利用多为变异组学技术手段系统阐述遗传变异对选择和驯化的效应，实现了目标群体的多维组学选择，提高可猪育种的准确性，可精准获得高生产效率的种猪。
实施例
实施例1hsd17b1基因表达量检测
1、实验样品采集
用于检测hsd17b1基因表达的大白猪和梅山猪，来源于天津市武清区中国农业科学院北京畜牧兽医研究所猪场。采集怀孕期49天母猪的卵巢，样品均储存在液氮中，提取总rna。每个品种均采集三头个体作为生物学重复。
2、hsd17b1基因表达量检测
2.1卵巢组织总rna提取
rna提取过程中所用玻璃器皿及镊子需经过200℃干烤4h以灭活rnase，eppendorf管及移液器枪头均需使用rnase灭活产品，整个实验过程中实验人员需佩戴口罩及乳胶手套。
(1)样品采集：将收集的卵巢组织置于1.5ml的eppendorf离心管中，按照rna提取说明书进行操作；
(2)用研磨杵冰上研磨组织，加入1mltrizol裂解液,剧烈震荡30s，室温静置10min，充分裂解组织细胞；
(3)加入0.2ml氯仿(氯仿:trizol体积比为1:5)，剧烈振荡15s，室温静置10min，在4℃的情况下，12000g离心15min，重复该步骤直至抽提完全；
(4)吸取上层水相(约450μl)到新的eppendorf管中，加入等体积异丙醇，轻轻混匀然后室冰上静止20min，在4℃的情况下，12000g离心10min；
(5)离心后观察到离心管底有少量白色胶状沉淀，弃上清，然后加入1ml75％乙醇(用depc水配置)洗涤沉淀，在4℃的情况下，12000g离心10min，弃去乙醇，于超净工作台风干5-10min至沉淀呈半透明状；
(6)加入depc水溶解rna沉淀，视沉淀量多少决定加入depc水量；
(7)紫外分光光度计测定样品rna浓度。
2.2总rna质量检测方法
取所提的样本rna1μl，利用agilent2100bioanalyzer和琼脂糖凝胶电泳的方法检测rna的质量。取rin值(rnaintegritynumber，rna分子完整数)大于8的样品为建库样品。同时采用2％的琼脂糖凝胶电泳检测rna的完整性，选取完整性高的样品进行后续mrna定量。
2.3mrna的反转录
用紫外分光光度计测定rna浓度后，取2μgrna加入2μloligo(dt)，用depc水补足体积至13μl，混匀后72℃水浴5min，然后冰上放置3-5min，离心后加入如下物质，混匀后，42℃水浴60min，-20℃保存。
表1反转录反应体系
2.4定量pcr
real-timequantitativepcr反应按照takarasybrpremixextaq(takara，drr041s)试剂盒说明书进行。每个样品设三个重复，反应在abiprism7500实时荧光定量pcr系统中进行，β-actin为内参。pcr反应体系见表2，反应条件见表3，引物序列见4。
表2rt-qpcr反应体系
表3rt-qpcr扩增程序
表4mrnaqpcr引物
2.5结果分析
梅山猪原产于我国上海嘉定、江苏太仓等地区，是我国猪种中繁殖力最强、产仔数最多的品种之一，比大白猪等现代欧洲猪种同窝产仔数高30％～40％，也是迄今为止世界上繁殖力最强的猪种之一。将梅山猪与大白猪hsd17b1基因表达量进行比较，可有效证明hsd17b1基因表达量对产仔性状的影响。检测结果见图1。
由图1可知，高产仔数的梅山猪中hsd17b1的表达量显著低于大白猪。hsd17b1是hsd17β家族中短链脱氢还原酶的典型代表，在雌激素合成过程中有着非常重要的功能。hsd17b1基因可能通过影响激素的合成及代谢，从而影响猪产仔性能。
实施例2snp分子标记的获得
1、实验样品采集
来自山东某种猪场的375头大白猪，用于hsd17b1基因snp分型。逐个采集耳组织样，放于75％的酒精中并-20℃保存，待提dna。
2、hsd17b1基因snp分型
2.1猪基因组dna的提取与检测
(1)剪取适量的猪耳组织，剪碎将其放进1.5ml的axgen管中。
(2)取50mlbd离心管，将终浓度为0.4mg/ml的蛋白酶k和裂解缓冲液进行混合均匀，向装有猪耳朵组织的1.5ml的离心管中加入0.5ml裂解液进行裂解。
(3)离心管平行均匀的放置于恒温杂交炉的摇板上(严封管盖，以免液体外渗)。55℃放置6个小时以上(消化过程中样品的充分混合是非常重要的。判断依据是无明显肉眼可见的猪耳组织，消化后的混合液成乳状)；
(4)待样品充分裂解后，取出离心管，在每管中加入0.3ml饱和氯化钠溶液，进行颠倒充分混合6-8次，然后放于冰上，冰浴15分钟；
(5)冰浴过后，12000rpm室温离心15分钟，小心缓慢的将上清转移至新的1.5mlaxgen离心管中(小心避免沉淀随上清一起倒出，保持倒的手法要一致，保持每管倒得上清在量上保持相同)；
(6)加入0.7ml的异丙醇(加入异丙醇的量随倒出的上清的量的变化而变化，二者等体积)于各管，颠倒直至溶液中有絮状沉淀出现(若无絮状沉淀，可将溶液-20℃冰箱放置2小时或4℃冰箱放置过夜)；
(7)12000rpm室温离心15分钟，去掉上清液体(此过程中留意观察离心管底部的白色dna沉淀，切勿将其与上清一起倒掉)；
(8)在各离心管中加入0.5ml70％乙醇，轻轻颠倒以充分冲洗上述沉淀出来的dna；
(9)10000rmp离心30s，用200μl微量移液枪将离心管中的乙醇吸掉，将沉淀的dna留在管中(当心勿将管中的dna沉淀吸走，此步骤所用的枪头在各离心管间操作时可不用更换)；
(10)自然风干dna10分钟；
(11)移液枪取0.1ml的te缓冲液于各管中，将上述dna沉淀重新溶解，置其于55℃放置2小时，定时摇晃数次，以使dna充分溶解；
(12)待dna充分溶解后，用紫外分光光度计测定所提浓度(基因分型需要标化样本dna浓度和od值，dna浓度15-20ng/μl，体积为30μl，a260/230介于1.5～2.3)，并琼脂糖凝胶电泳检测所提质量(可见单一条带)。
(13)将dna溶液放置4℃过夜，次日取1μl进行pcr(若数星期内使用所提dna，可在4℃保存；如若长时间不使用，需放置于-20℃环境下)。
2.2snp芯片基因型判定及基因型数据的质控
将375个大白猪的dna样品，snp的相关信息交给北京康普森生物有限公司进行sequenommassarray法核酸质谱测序。
以样本dna为模板，加入pcr扩增引物(p1和p2)对以及单碱基延伸引物进行pcr扩增反应，反应产物纯化后与质谱芯片共结晶，用飞行质谱法检测并分型，检测目的基因的snp位点。
2.3数据整理与分析
1)表型数据分析
利用sas9.2统计分析软件，对猪初产产仔测定值和经产(2-8次)产仔测定值进行描述性统计分析。
2)单倍型分析
单核苷酸多态性的基因型和等位基因频率的计算是使用popgene3.2实现的。利用sas9.2软件中的glm过程分析snp和产仔数性状的关联分析，固定效应模型如下：
y＝μ+gi+mk+e
y为繁殖性状的表型记录；μ为性状的总平均值；gi为基因型效应；mk为生产月份效应；e为随机误差。数据以概率值和平均值±标准差表示，p值＜0.05表示具有显著性差异。
2.4结果分析
通过上述sequenommassarray，测得大白猪hsd17b1基因snp分型结果，如表5所示。
表5
基因型检测结果表明：5头猪的基因型为aa基因型，19头猪的基因型为ag基因型，67头猪的基因型为gg基因型。猪hsd17b1基因在检测猪群中的基因型频率结果：aa基因型频率为0.055，ag基因型频率为0.209，gg基因型频率为0.736，gg和ag的基因型频率高于aa，gg和ag基因型为优势基因。
利用sas9.2软件对基因型和繁殖性能进行统计分析，并做样本间多重比较。结果如表6所示。
表6
从表6可以看出，hsd17b1不同基因型间初产产仔数和平均胎重差异极显著。aa型的产仔数低于ag和gg型，ag型的初产平均胎重低于aa和gg型，但窝重高于aa和gg型。经产数据中，ag和gg型的产仔数、窝重和平均胎重均高于aa型。
因此，snp(rs323407415)位点明显对猪的产仔性状有影响，在实际的猪育种中，ag和gg型基因型猪具有更高的繁殖效率。
综上所述，可以通过确定猪hsd17b1基因rs323407415位点的核苷酸来确定猪个体为aa基因型还是ag和gg基因型，从而辅助鉴定初产/经产产仔数和平均胎重：ag和gg型基因型猪具有更高的繁殖效率。
aa基因型为猪hsd17b1基因第323407415位的核苷酸为a的纯合体。
ag基因型为猪hsd17b1基因第323407415位的核苷酸为a与g的杂合体；
gg基因型为猪hsd17b1基因第323407415位的核苷酸为g的纯合体。
实施例3检测相关snp位点的试剂盒
本发明中snp位点检测试剂盒包括以下组分：
(1)，dntp(2.5mm)；
(2)，5×pcr缓冲液；
(3)，fastpfudna聚合酶(250u)；
(4)，snapshot混合液；
(5)，mgso4溶液(50mm)；
(6)，6×dnaloadingbuffer；
(7)，pcr扩增引物(p1和p2)：其核苷酸序列为seqidno.4和seqidno.5所示；
(8)，延伸引物：其核苷酸序列为seqidno.6所示；
(9)，核酸外切酶(250u)；
(10)，虾-碱性磷酸酶(100u)；
(11)，纯化用缓冲液：；
(12)，ddntp；
(13)，snp位点纯合突变的阳性对照标本：含有猪rs323407415位点人工合成为g碱基的质粒；
(14)，snp位点杂合突变的阳性对照标本：含有猪rs323407415位点人工合成a和g碱基的两种质粒的混合物；
(15)，snp位点阴性对照标本：人工合成的带有随机碱基的质粒。
通过本发明的snp位点检测试剂盒及其使用方法，对猪hsd17b1基因第323407415位snp位点进行检测。
1、dna提取
剪取适量的猪耳组织，提取全基因组dna。
2、pcr扩增反应
pcr反应试剂(20μl)包括：1μl正向引物组(10μm)，1μl反向引物组(10μm)，4μl5×pcr缓冲液，2μldntp(2.5mm)，0.2μlfastpfudna聚合酶(250u)，h2o补充至20μl。待测dna样品(200ng)、snp纯合突变阳性对照标本(100ng)、snp杂合突变阳性对照标本(100ng)、snp阴性对照标本(100ng)与上述pcr反应试剂构成pcr反应体系。
3、pcr产物鉴定与纯化
pcr扩增结束后，琼脂糖凝胶电泳进行pcr产物鉴定。pcr扩增产物用虾-碱性磷酸酶和核酸外切酶进行纯化处理，37℃酶切60min，80℃15min灭活虾-碱性磷酸酶和核酸外切酶。纯化体系包括：2μl终浓度虾-碱性磷酸酶(0.6u)，5μl终浓度核酸外切酶(1.2u)，50ngpcr扩增产物，和2μl纯化用缓冲液。
4、延伸反应
对扩增产物进行回收，之后在含有0.1mmddntp的体系中使用延伸引物进行pcr扩增。
5、毛细管电泳及snp位点分析
使用abi3730xl测序仪对延伸引物测序。使用genemapper系列软件对结果峰图进行snp位点上脱氧核苷酸的判定。
以上结合具体实施方式和范例性实例对本发明进行了详细说明，不过这些说明并不能理解为对本发明的限制。本领域技术人员理解，在不偏离本发明精神和范围的情况下，可以对本发明技术方案及其实施方式进行多种等价替换、修饰或改进，这些均落入本发明的范围内。本发明的保护范围以所附权利要求为准。
sequencelisting
<110>中国农业科学院北京畜牧兽医研究所
<120>猪多组学整合精准育种方法
<130>2010
<160>6
<170>patentinversion3.5
<210>1
<211>2213
<212>dna
<213>猪（susscrofa）
<400>1
gcctcacaggagttccccaggcctgggtgaaggaactcatgctggtgggatcaggcatgg60
aaggtttggaggggcctgagctgccttcctccctcatctcaaatttgggatcaaagtact120
gttttaatacaacttaataaactgatgaatcacctggacctgtgtgtcggaatgtcagtg180
ccatgacacctgggaggactcaaggctatcgcatagacactggtgaggccacttctccct240
gcccccttgcctcagctcccccattttaaaaagtaacctttggcctttgaacggagactg300
ccctggagttcccttcgtggcacagcggtaatgaaccagactagtatccatgaggatgtg360
tattcaatccctggcctcgctcagtgggttaaggatctggcgctgccatgagctgtggtg420
taggtcacagactccttggatcctgagttgttgtgtagagtggcagctacagctctgatt480
cgacccctagcctgagaacttccatatgcccctaaaaagacaataagtaaacgaacagac540
acactgagaatgcccttccgaagcatttgctttcacaccctcagcaacaaaccttggcaa600
ctaggattgtccctgtgtgtcagaggcttgctaggggaagatcagctctggggcactcag660
tgcccacacctcctgtttgcaaagcttcacccaccctctgaccagcaactaaaaggggat720
gttccaggcagtggcctcatctctgccttgacctagagaaaattggggagagcaggccag780
caacccacaggagggctcagccttgagattgtgctgaagacataataaggggtggggacc840
ccagcgggcaggggcgggtcacagcagatgataactatcccggagaccagagcccagcat900
tcctcgcctttccaaactgcacaacccaatggaccgcaccgtggtgctcatcaccggctg960
ttcctccggcatcggcctgcacctggccctgcgtctggcatctgacccatctcagagctt1020
caaagtgtacgccaccctaagggacttgacggcacaaggcccgctgttggaggccgccca1080
ggcccgagggtgccccccgggctccctggagaccttgcagctggatgtgagggatgcaga1140
ttcaatagcggatgcccgggcacacgtgaccgagggacgcgtggatgtgctggtgtgcaa1200
cgcaggccggggcctggtggggcccctggaggcacacagggagggcgccgtggactccgt1260
gctggatgtgaacctggccggcaccgtgcggatgatgcaggccttcctgccagacatgaa1320
gcggcgccgctcgggtcgcatattggtgaccgggagcctgggaggcctgatagggcttcc1380
cttcaacgctgtttactgcgccagcaagtttgcaatcgaaggcttatgcgagagtctggc1440
ggttgtgctgcagagcttcggggtccacgtgagcgtcatcgagtgcggcccggtgcgtac1500
tgccttccccgagaagctggaggacggcctgggcgggatactggatcgcgcggacgccga1560
gacccgtgacctcttctcccgctacctgagccacttcgagcagacgttccttgaggcggc1620
gcaggacccgaaggaggtggtcgaggtcttcctccaggctctgcgcgcgccgcgcccggc1680
tctgcgctacttcaccaccgagcacttccatcccctgatcaagctgcgcttctccgaccc1740
cagcagcagcagctacgtcgcagccgagcaccaaagagtgttcggcgaccaggccaccga1800
gggcttcgagggcactgattgcggcgaggcggaggcgggagctggagacctggggccgtc1860
tgagctcggcgctcctcttgccaccccgcaataagggcttggtctgccgctgtctcccgc1920
actctcctttgcggcccggaaggtgtgcggtttagggatggaaaggggctcctgactgct1980
ctgcacttgcgcctgcagagggtgagccgtcaggaaaggtcgtgttagcccgggttgtac2040
ctactgaaaccgaagcccagaaagggaaggagctttgcccaaggtcactcctcgtcctcc2100
acagggcagggagagtggtcctcgatccccaaggggtgtctgcagcaggaggcgctagtc2160
ggacttgagtgcggggagagttggcttgtgacttggcttccaccagtcatggg2213
<210>2
<211>20
<212>dna
<213>猪hsd17b1基因mrna的表达量pcr扩增引物（susscrofa）
<400>2
tccaatggagaagatggcgg20
<210>3
<211>20
<212>dna
<213>猪hsd17b1基因mrna的表达量pcr扩增引物（susscrofa）
<400>3
ggagatgtccggttcagagc20
<210>4
<211>29
<212>dna
<213>猪hsd17b1基因snp位点扩增引物（susscrofa）
<400>4
acgttggatgacgggcttcgcttgcttct29
<210>5
<211>30
<212>dna
<213>猪hsd17b1基因snp位点扩增引物（susscrofa）
<400>5
acgttggatgaaacagcgttgaagggaagc30
<210>6
<211>15
<212>dna
<213>猪hsd17b1基因snp位点延伸引物（susscrofa）
<400>6
cttgcttctgctgct15