一种基于mSNP技术检测白萝卜种子纯度的混样检测方法与流程

一种基于msnp技术检测白萝卜种子纯度的混样检测方法
技术领域
1.本发明属于种子纯度检测领域，具体涉及一种基于msnp技术检测白萝卜种子纯度的混样检测方法。


背景技术：

2.白萝卜是十字花科萝卜属一年生作物，染色体数目18对，营养成分丰富，在世界各地均有种植，开发价值较大。中国是白萝卜的发源地，是我国重要的蔬菜之一，已有2700年以上的栽培历史，具有丰富的种质资源。
3.随着市场种子产业的兴起，种子产品日渐丰富，种子的生产经营单位也越来越多，市场上种子纯度问题普遍存在，严重干扰市场秩序。种子纯度为本品种的种子数占供检本作物样品种子数的百分率。为对种子进行有效区分，保证种子纯度，在萝卜杂交种的新品种推广中一般要试种一季，通过田间种植的方法确定杂交种的纯度之后才能在市场上推广应用，这个过程不仅增加了田间管理和储藏成本，而且延误了种子上市的时机。为了加快优质种质的上市，促进萝卜产业的发展，对快速准确的的种子纯度鉴定方法有迫切的需求。
4.种子纯度鉴定经过形态学、生物化学水平，已经进入了dna水平的研究水平。dna分子标记信息量大、遗传性稳定、不易受内外环境影响、可以直接反映dna在分子水平上的遗传变异、对所需的dna样品量较少、操作方便，是可用于应用的理想的检测方法。得益于dna分子标记的诸多优点，近年来在作物领域如水稻、小麦等物种的种子纯度鉴定中应用已经比较广泛。常用的分子标记rapd、ssr、rflp等；rflp是最早出现的第一代分子标记，优点是遗传信息丰富、稳定性和重复性较好；缺点是费时、对dna样本的质量要求高、需要用到放射性探针等。ssr标记的优点是重复性高、多态丰富、共显性、操作简单；缺点是引物设计及筛选工作量较大，成本较高。rapd标记优点是操作灵活、所需的dna样品少、操作方便；缺点是对实验条件依赖度高，影响结果的可靠性。


技术实现要素：

5.本发明的目的在于克服现有技术的缺陷，提供一种基于msnp技术检测白萝卜种子纯度的混样检测方法，其高效、准确、成本低。
6.为了实现上述目的，本发明采取的技术方案如下：
7.技术方案一：
8.一种于msnp技术的白萝卜种子纯度检测的引物组，其特征在于，其包括引物对1f/r、引物对2f/r、引物对3f/r、引物对4f/r、引物对5f/r、引物对6f/r、引物对7f/r、引物对8f/r、引物对9f/r、引物对10f/r、引物对11f/r、引物对12f/r、引物对13f/r、引物对14f/r、引物对15f/r、引物对16f/r、引物对17f/r、引物对18f/r、引物对19f/r、引物对20f/r；其中每条引物对均由正向引物和反应引物组成；
9.所述引物对1f/r中，f引物的序列seqi dno.1所示，r引物的序列如seqi dno.2所示；
10.所述引物对2f/r中，f引物的序列seqi dno.3所示，r引物的序列如seqi dno.4所示；
11.所述引物对3f/r中，f引物的序列seqi dno.5所示，r引物的序列如seqi dno.6所示；
12.所述引物对4f/r中，f引物的序列seqi dno.7所示，r引物的序列如seqi dno.8所示；
13.所述引物对5f/r中，f引物的序列seqi dno.9所示，r引物的序列如seqi dno.10所示；
14.所述引物对6f/r中，f引物的序列seqi dno.11所示，r引物的序列如seqi dno.12所示；
15.所述引物对7f/r中，f引物的序列seqi dno.13所示，r引物的序列如seqi dno.14所示；
16.所述引物对8f/r中，f引物的序列seqi dno.15所示，r引物的序列如seqi dno.16所示；
17.所述引物对9f/r中，f引物的序列seqi dno.17所示，r引物的序列如seqi dno.18所示；
18.所述引物对10f/r中，f引物的序列seqi dno.19所示，r引物的序列如seqi dno.20所示；
19.所述引物对11f/r中，f引物的序列seqi dno.21所示，r引物的序列如seqi dno.22所示；
20.所述引物对12f/r中，f引物的序列seqi dno.23所示，r引物的序列如seqi dno.24所示；
21.所述引物对13f/r中，f引物的序列seqi dno.25所示，r引物的序列如seqi dno.26所示；
22.所述引物对14f/r中，f引物的序列seqi dno.27所示，r引物的序列如seqi dno.28所示；
23.所述引物对15f/r中，f引物的序列seqi dno.29所示，r引物的序列如seqi dno.30所示；
24.所述引物对16f/r中，f引物的序列seqi dno.31所示，r引物的序列如seqi dno.32所示；
25.所述引物对17f/r中，f引物的序列seqi dno.33所示，r引物的序列如seqi dno.34所示；
26.所述引物对18f/r中，f引物的序列seqi dno.35所示，r引物的序列如seqi dno.36所示；
27.所述引物对19f/r中，f引物的序列seqi dno.37所示，r引物的序列如seqi dno.38所示；
28.所述引物对20f/r中，f引物的序列seqi dno.39所示，r引物的序列如seqi dno.40所示。
29.进一步的，引物对1f/r～20f/r由msnp技术获得。
30.技术方案二：
31.一种根据上述的引物组检测白萝卜种子纯度的混样检测方法，其包括如下步骤：
32.步骤1、选材：选取1个或多个白萝卜品种；每份份白萝卜样品至少采用96粒种子；
33.步骤2、对白萝卜基因组dna进行准确定量；
34.步骤3、将所述的引物组中进行引物合成，每个引物对中的正向引物和反向引物合成时，均合成10条带有不同目标标签的引物；然后按照指定的标签组合进行引物的混合制备引物混合液；
35.步骤4、以白萝卜基因组dna为模板，用引物混合液分别对白萝卜的基因组dna进行一轮pcr扩增，获得目标区域；
36.步骤5、将所获得的pcr扩增产物，进行等量混合；
37.步骤6、混合后的产物进行片段筛选；
38.步骤7、对筛选后所得体系中的单链dna进行消化；
39.步骤8、对消化后的产物进行纯化；
40.步骤9、在步骤8中获得的体系中配置二轮pcr体系；
41.步骤10、对二轮pcr产物进行纯化，完成测序文库的制备；
42.步骤11、将测序文库等质量混合后上机测序，获得测序数据；
43.步骤12、对获得的测试数据，再次根据标签组合将样本拆分开；
44.步骤13、识别测试样目标位点的基因型结果，通过位点基因型情况判定种子的纯度。
45.进一步的，引物对中的正向引物和反向引物合成时，均合成10条带有不同目标标签的引物；
46.每个引物对的标签组合中，所述正向引物的标签序列与所述反向引物的标签序列不同。
47.进一步的，步骤3中，每个引物对中，合成的10条带有标签的正向引物中的标签序列分别如seqi dno.45～54所示；
48.每个引物对中，合成的10条带有标签的反向引物中的标签序列分别如seqi dno.55～64所示。
49.进一步的，每对引物对合成的10条带有标签的引物对中，正向引物和反向引物的标签组合方式选自表1中的任意一种或几种。
50.表1：标签组合方式
51.52.[0053][0054]
进一步的，步骤3中，引物对1f/r～20f/r中的f引物还包括f端通用引物，所述f端通用引物的序列如seqi dno.41所示；引物对1f/r～20f/r中的r引物还包括r端通用引物，所述r端通用引物的序列如seqi dno.42所示；
[0055]
步骤9中、二轮pcr所用的frimer f的序列如seqi dno.43所示，
[0056]
所用的primer r的序列如seqi dno.44所示。
[0057]
更进一步的，当白萝卜品种为多个时，所述primer r的序列还包括用于区别白萝卜品种的条形码序列。
[0058]
进一步的，步骤2中，每个引物对的标签组合中，所述正向引物的标签序列与所述反向引物的标签序列不同。
[0059]
进一步的，在步骤1中，采用白萝卜种子基因组dna提取试剂盒对白萝卜种子的基因组dna进行提取。
[0060]
进一步的，在步骤4中，所述一轮pcr扩增体系：引物混合液8μl；dna用量100ng；3*t酶15μl；加水补足45μl；
[0061]
所述一轮pcr扩增程序：95℃3min；(95℃30s，60℃4min，72℃30s)28个循环；72℃4min。
[0062]
进一步的，在步骤6中，混合后的产物进行片段筛选，具体操作如下：
[0063]
步骤6.1、加入一轮pcr体积0.4倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，取上清转移至新的管中；
[0064]
步骤6.2、加入一轮pcr体积0.6倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0065]
步骤6.3、添加一轮pcr体积0.9倍的磁珠悬浮液，重悬磁珠，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0066]
步骤6.4、加入100μl的体积浓度为80％的乙醇，用磁力架反复在不同的两面吸附磁珠，使磁珠得到充分的洗涤，用磁力架吸附2min，去除上清，室温放置至乙醇挥发干净。
[0067]
在步骤7中，对筛选后所得体系中的单链dna进行消化，具体操作步骤为：
[0068]
步骤7.1、向所获产物中加入20μl水，混匀磁珠；
[0069]
步骤7.2、吸附磁珠，转移16μl上清至新的ep管中；
[0070]
步骤7.3、向体系中加入exo i 2μl，10
×
reaction buffer 2μl；
[0071]
步骤7.4、消化体系的消化程序为：37℃30min；85℃15min；
[0072]
在步骤8中，对消化后的产物进行纯化的具体操作步骤为：
[0073]
步骤8.1、加入0.9倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0074]
步骤8.2、添加等pcr体积的磁珠重悬液，重悬磁珠，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0075]
步骤8.3、加入100μl的体积浓度为80％的乙醇，用磁力架反复在不同的两面吸附磁珠，使磁珠得到充分的洗涤，用磁力架吸附2min，去除上清，室温放置至乙醇挥发干净。
[0076]
进一步的，在步骤9中，所述二轮pcr体系：3*t酶10μl；primer f；primer r；h2o18μl；
[0077]
所述二轮pcr程序：95℃3min；(95℃15s，58℃15s，72℃30s)12个循环；72℃4min；
[0078]
进一步的，在步骤10中，对二轮pcr产物进行纯化，利用的是0.80倍的磁珠，具体操作如下：
[0079]
步骤10.1、加入0.8倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0080]
步骤10.2、添加等pcr体积的磁珠重悬液，重悬磁珠，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0081]
步骤10.3、加入100μl的80％乙醇，用磁力架反复在不同的两面吸附磁珠，使磁珠得到充分的洗涤，用磁力架吸附2min，去除上清，室温放置至乙醇挥发干净；
[0082]
步骤10.4、加入23μl elution buffer，充分悬浮磁珠，室温静置2min以洗脱dna，将磁珠用磁铁吸附，所得到的上清dna溶液吸至一新的管中，获得测序文库；所述elution buffer为10mm tris
‑
hcl，ph 8.0
‑
8.5。
[0083]
与现有技术相比，本发明所取得的有益效果如下：
[0084]
msnp技术：本发明采用msnp技术，一个扩增子内对应多个snp，最大限度的利用了每个扩增子片段所获得的的信息，可以在扩增子不变的情况下检测到尽量多的snp变异。且这些snp之间可构成单倍型，提高了变异的检测效率。这样不仅可以采用msnp位点内和位点
间的变异，而且还可以采用单倍型和snp两种方式进行检测，使得遗传变异的检测更加精细，提高了标记鉴定的准确度和灵敏度。本专利中采用了20对引物对，实际检测到的变异信息有100多个，有助于筛选出更多的符合需求的多态性位点用于纯度鉴定分析。与常规snp检测相比，减少了引物对数的使用，降低了成本。
[0085]
成本低：采用msnp技术，扩增子不变的情况下可以检测到更多的snp变异；采用一轮扩增后至少96个测试样产物混样检测的方案，减少二轮扩增工作量的同时降低了测序成本。
[0086]
高效：采用混样检测的方案，一轮扩增后至少96个检测样的产物混样进行后续的扩增和测序，可在一天时间内最多可完成1440粒种子的检测。与其他检测方法相比，操作简单，不用田间种植、不需要工作量大、检测难度大的实验操作等，可以快速进行种子纯度的检测工作。
[0087]
准确：本方法是利用测序技术直接读取单核苷酸多态，通过程序直接进行纯度结果的判读，结果直观，可靠，避免受种子发育期、及结果判读时主观因素的影响。
附图说明
[0088]
图1是本发明混样检测方法的流程图。
具体实施方式
[0089]
实施例1：特异性引物的获得方法
[0090]
本发明特异性引物的获得方法，具体为：
[0091]
利用白萝卜的全基因组重测序数据，采用bwa
‑
mem(http://bio
‑
bwa.sourceforge.net/)回贴到白萝卜参考基因组上，利用gatk(https://software.broadinstitute.org/gatk/)进行单核苷酸变异鉴定。
[0092]
鉴定出的单核苷酸变异位点集，筛选最小等位变异频率>0.02、杂合率＜15％、缺失率＜20％，将单核苷酸变异位点进行合并，筛选单核苷酸位点数位于6
‑
10的区段，即msnp(多聚单苷酸多态性)位点，相较于传统的snp(单核苷酸多态性)位点，msnp位点可最大限度的利用每个引物对所获得的信息，即在引物对不变的情况下检测到尽可能多的snp位点，且同一引物对内的所有snp还可以组合成单倍型，使其多态性更高。例如一个snp存在a/t两种变异，其可区分的多态为aa、at、tt共计3种，若检测的为msnp位点，假如一个引物对内有3个snp(a/t、g/t、c/a)，则可以存在8种(agc、aga、ata、atc、tgc、tga、tta、ttc)多态。这使得遗传变异的检测更加精细，同时提高了标记鉴定的准确度和灵敏度。本专利筛选多态性比较高的区段共计37个。
[0093]
对37个目标区段进行引物设计，并对引物特异性进行筛选，共获得染色体特异引物20对，共计131个单核苷酸变异位点。从检测成本和实用角度综合考虑，同样按白萝卜样本间5
‑
10％的多样性原则，我们最终选取20对引物混合进行白萝卜种子纯度的检测，共计可检测130个左右snp位点。
[0094]
在本发明中，共计获得20组特异性引物对，即引物对1f/r～20f/r，引物对1f/r～20f/r的基因序列如seqi dno.1～40所示。
[0095]
实施例2：白萝卜种子纯度检测所用引物组
[0096]
白萝卜种子纯度检测所用的引物组，包括引物对1f/r～20f/r，在引物对1f/r～20f/r中，不仅包括如seqi dno.1～40所示特异性引物序列，还包括通用引物序列，引物对1f/r～20f/r中的f引物的f端通用引物序列如seqi dno.41所示；引物对1f/r～20f/r中的r引物的r端通用引物的序列如seqi dno.42所示；
[0097]
实施例3：引物混合液的获得方法：
[0098]
获得特异性引物后，设计特异的标签序列，再重新进行引物合成，此时引物合成时会添加上特异的标签序列，本实施例采用96组特异的标签组合，具体为：
[0099]
依据特定标签组合情况，将每条特异性引物合成10条带有不同目标标签的引物，带有目标标签的引物序列形式如表2所示，将带有目标标签的所有f(正向引物)和所有r(反向引物)，每条引物取10μl，定容至10ml；每条引物的浓度0.1μm，本实施例共计制备出含有96组特异标签组合的引物对，即引物混合液；
[0100]
表2引物组
[0101]
引物对f正向引物(5’～3’)r反向引物(5’～3’)1f/rfffyyygtgtaaagacattgacaagcagacarrryyytgtttacaagaataacaatagttgtcca2f/rfffyyycatccttgtacttcaccaatccttgrrryyycatagaaacccaacggtaatgatgt3f/rfffyyyacgaagttcacttagttttctctgtrrryyygtgtacagagataaacttgtccgat4f/rfffyyyatcagtttggtcctaatcgtgatccrrryyycatttactaggaccagcagagatct5f/rfffyyyggaacatatttgaacagctgcagaarrryyyttgaaagcaacaatcgtagcaatgc6f/rfffyyyaagatgagggatgtcttggaatgatrrryyyaatgcaagcaatccgaactccattt7f/rfffyyycttggagaatcaggtttgtttgctarrryyyttcttctgctctttttccttcaacc8f/rfffyyycaaaacgccacaaaagaaagatcagrrryyyctctttggtgcgtagaacatacttc9f/rfffyyycttcatcgccatgtttcctctttttrrryyycaaggaagctagctttagtgattcg10f/rfffyyycaaagaacgtcaagtttcaacacacrrryyycaatcgtgatgaagagatgagcatc11f/rfffyyygatcatatggacggtagggaaagaarrryyytacaacaaggacaatcgaacaaaca12f/rfffyyytgttccaccgtaagtacatagagaarrryyyggaacgatgcgatggaagat13f/rfffyyyaggattcagatgttgtagctattgarrryyyccttgctcacaacttccctca14f/rfffyyyactggaaatgaagaaactatgtctctrrryyytagtgtcttctcgaatgtttcccat15f/rfffyyyttcttttgcttcttgagatcaggaarrryyyttcttatcctcccttaaactcagcg16f/rfffyyycattaactggtcccgtctctgtctcrrryyyaccaaaacaaaagatagtaatcgaaagga17f/rfffyyyggtagatttgaatgttgaggtctgcrrryyyacccaacctaaagagcttgtcttaa18f/rfffyyyagtaaagcatgattccatcctaaagcrrryyytctggactagaagatgttgaggttg19f/rfffyyyaatgatgagaagggagatggagatgrrryyyttcagatgtttacagaggcacaaag20f/rfffyyyacacataatttcagttttcagacagarrryyyttgaagacttactatgatctcggcc
[0102]“fff”为f端通用引物序列，所述f端通用引物序列为aacgacatggctacgatccgactt，如seqi dno.41所示；
[0103]“rrr”为r端通用引物序列，所述r端通用引物序列为ctaagaccgcttggcctccgactt，如seqi dno.42所示。
[0104]
其中“yyy”为标签序列，
[0105]
每个引物对中，合成的10条带有标签的正向引物中的标签序列分别为ccttc、accga、atgtg、aatgc、ttcgg、aaggt、cccat、atgga、acgat、ctctg，即分别如seqi dno.45～54
所示；
[0106]
每个引物对中，合成的10条带有标签的反向引物中的标签序列分别为atccg、tatcg、actcg、taacc、cttac、tccta、acact、tacgt、tcacg、acgca，即分别如seqi dno.55～64所示。
[0107]
每对引物对合成的10条带有标签的引物对中，正向引物和反向引物的标签组合方式选自表2中的任意一种或几种。
[0108]
实施例4：纯种鉴定
[0109]
步骤1、材料的选取，选取2份白萝卜(每份白萝卜分别选取96粒种子)，分别编号为cc2
‑
22和cc2
‑
18，提取白萝卜种子的基因组dna，采用石家庄博瑞迪生物技术有限公司生产的白萝卜种子基因组dna提取试剂盒对提取后的dna进行准确定量。
[0110]
步骤2、以白萝卜种子基因组dna为模板，用引物混合液进行一轮pcr扩增，获得目标区域；
[0111]
所述一轮pcr扩增体系：实施例3获得的引物混合液8μl；dna用量100ng；3*t酶15μl；加水补足45μl。
[0112]
所述一轮pcr扩增程序：95℃3min；(95℃30s，60℃4min，72℃30s)28个循环；72℃4min。
[0113]
步骤3、将所获得的pcr扩增产物，进行等量混合，只有扩增产物为带有不同组合的特定标签才可进行混合，混合时直接按照体系1:1混合即可。等量混合后，对混合后的产物进行产物纯化即片段筛选，具体步骤为；
[0114]
步骤3.1、加入一轮pcr体积0.4倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，取上清转移至新的管中；
[0115]
步骤3.2、加入一轮pcr体积0.6倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0116]
步骤3.3添加一轮pcr体积0.9倍的磁珠悬浮液，重悬磁珠，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0117]
步骤3.4加入100μl的体积浓度为80％的乙醇，用磁力架反复在不同的两面吸附磁珠，使磁珠得到充分的洗涤。用磁力架吸附2min，去除上清，室温放置至乙醇挥发干净；
[0118]
所述磁珠为：诺唯赞磁珠
[0119]
步骤4、对筛选后所得体系中的单链dna进行消化；
[0120]
在步骤3中获得的含磁珠的体系中，进行如下操作：
[0121]
步骤4.1、向所获产物中加入20μl水，混匀磁珠；
[0122]
步骤4.2、吸附磁珠，转移16μl上清至新的ep管中；
[0123]
步骤4.3、向上述体系中加入exo i 2μl，10*reaction buffer 2μl。
[0124]
步骤4.4、上述体系的消化程序为：37℃30min；85℃15min。
[0125]
步骤5、对消化后的产物进行纯化：
[0126]
步骤5.1、加入0.9倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0127]
步骤5.2、添加等pcr体积的磁珠重悬液，重悬磁珠，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0128]
步骤5.3、加入100μl的体积浓度为80％的乙醇，用磁力架反复在不同的两面吸附磁珠，使磁珠得到充分的洗涤。用磁力架吸附2min，去除上清，室温放置至乙醇挥发干净。
[0129]
步骤6、在步骤5获得的体系中配置二轮pcr体系：
[0130]
在步骤5中获得的含磁珠的体系中，配置二轮pcr体系，进行二轮pcr扩增；
[0131]
所述二轮pcr体系：3*t酶10μl；primer f；primer r；h2o 18μl
[0132]
所述二轮pcr程序：95℃3min；(95℃15s，58℃15s，72℃30s)12个循环；72℃4min。
[0133]
所述primer f的序列如seqi dno.43所示，为gaacgacatggctacgatccgactt；所述primer r的序列如seqi dno.44所示，为tgtgagccaaggagttgttgtcttcctaagaccgcttggcctccgactt；
[0134]
由于本实施例采用的是2份白萝卜品种，为了对样本进行区分，因此，primer r的序列除包括如seqi dno.44所示的序列外，还包括唯一条形码序列barcode；
[0135]
带条形码的primer r的序列为：
[0136]
tgtgagccaaggagttgxxxxxxxxxxttgtcttcctaagaccgcttggcctccgactt；
[0137]
其中“xxxxxxxx”为用于识别样品的唯一条形码barcode，以区分样本。
[0138]
本实施例中2份白萝卜品种样本的barcode序列分别为cggctaaa、tccccgtg；即分别如seqi dno.65
‑
66所示。
[0139]
步骤7、利用0.80倍的磁珠对二轮pcr扩增产物进行纯化，完成测序文库的制备；
[0140]
步骤7.1加入0.8倍的磁珠，用移液器上下吹打混匀，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0141]
步骤7.2添加等pcr体积的磁珠重悬液，重悬磁珠，静置2min，用磁力架吸附，至溶液澄清，去除上清；
[0142]
步骤7.3加入100μl的体积浓度为80％的乙醇，用磁力架反复在不同的两面吸附磁珠，使磁珠得到充分的洗涤。用磁力架吸附2min，去除上清，室温放置至乙醇挥发干净。
[0143]
步骤7.4加入23μl elution buffer，充分悬浮磁珠，室温静置2min以洗脱dna。将磁珠用磁铁吸附，所得到的上清dna溶液吸至一新的管中，获得测序文库(elution buffer为10mm tris
‑
hcl，ph 8.0
‑
8.5)；
[0144]
步骤8、将测序文库等质量混合后上机测序，获得测序数据。
[0145]
步骤9、识别目标dna的基因型结果，从96粒种子中分别选取检测合成的种子用于后续纯度的判断依据，通过多态性高的位点基因型情况判定这批种子的纯度，本实施例的检测结果见表3。
[0146]
表3检测结果
[0147]
样品编号自交种子数种子纯度结论cc2
‑
22495.74％同常规鉴定结果cc2
‑
18297.85％同常规鉴定结果
[0148]
种子纯度程序判读原则：
[0149]
首先对每个位点是杂交/异交还是自交进行判断，判断标准：计算一批种子在每个位点基因型种类所占的比例。如果这个位点的基因型只有一种且纯合，这个位点舍弃，因为无法判断是自交还是杂交。如果一种杂合基因型的占比超过90％，判断为杂交；一个样本在这个位点的基因型是纯合就判断是自交；如果是其他杂合基因型就判断为异交。
[0150]
然后对样本进行自交、异交和杂交的判定。最后统计这批种子的整体情况。
[0151]
纯度计算，用一批种子纯度百分率表示：
[0152][0153]
以上所述实施方式仅为本发明的优选实施例，而并非本发明可行实施的穷举。对于本领域一般技术人员而言，在不背离本发明原理和精神的前提下对其所作出的任何显而易见的改动，都应当被认为包含在本发明的权利要求保护范围之内。