本发明涉及分子育种技术领域，更具体地，涉及一种基于叶绿体基因组序列的火炬松个体鉴定方法。
背景技术：
火炬松是一种原产北美东南部的高大乔木，目前在我国多地引种种植。其具有生长快、适应能力强、干型通直无节，侧枝较细、自然整枝能力强、材质好、松脂产量高等特点，已经被许多国家和地区广泛引种栽培。
火炬松遗传改良离不开种质资源，建立育种群体是保护种质资源的重要形式，是火炬松遗传改良的重要内容。火炬松遗传改良遵循的是多世代轮回选择的模式。在育种进程中，为了短期内获得高的遗传增益，必须提高选择强度；而随着世代和育种轮回的推进，在选择压力下，群体的遗传基础不可避免地逐步变窄，遗传多样性水平逐步下降，选择强度越大，遗传多样性水平下降越快，最终导致群体内亲缘关系普遍存在，近郊衰退严重，损害长期增益。育种工作者通过科学的群体管理，既获得较高的短期增益，又能满足长期育种的需要。群体管理的核心内容就是个体间遗传关系的管理和调控，需要从交配设计、选择方式、选择强度等方面着手进行。
目前，在火炬松高世代实际育种中，为了全面了解杂交组合的效应和避免群体内亲缘关系过快增加，采用完全谱系交配设计和人工控制授粉，这样不仅杂交组合多，而且耗时耗力、投入大。为了解决高世代育种面临的问题，可以通过混合花粉授粉的方法获得半同胞子代，根据子代测定的结果筛选出优良单株，继而利用分子标记鉴定拟选优树的父本来源，保证入选优树亲本来源清晰，从而控制育种群体的亲缘关系。
目前普遍利用基因组分子标记鉴定后代的父本来源，但是松树核基因组不仅体量大，而且重复序列多，在有性繁殖过程中基因发生重组，针对核基因组进行的亲本鉴定有很大难度。松树叶绿体基因组相对核基因组有许多优势。首先，叶绿体基因组包含适量的遗传信息，为比较研究提供了一个相对适中的数据基础；其次，叶绿体dna的核苷酸变异率适中，在应用上很有价值；除此以外，叶绿体基因组dna序列较小，便于测序；这为鉴定优树的亲本提供了十分有利的条件。火炬松基因组大小约为22gb，基因组中包含大量重复序列。重复序列的存在导致pcr扩增的特异性明显降低，绝大部份的分子标记扩增都需要设计巢式pcr引物进行扩增，特异性低，实验复杂。
叶绿体是绿色植物细胞中半自主的细胞器，有独立的基因组，基因组主要为环状双链分子。以往的研究发现，松属树种的叶绿体基因组属于父系遗传，但是在结构上极其高度保守。一直以来，叶绿体基因组序列主要是用于分类、进化和基因工程方面的研究，只限于物种水平的研究。
技术实现要素：
本发明的目的是为了克服现有技术的不足，提供一种基于叶绿体基因组序列的火炬松个体鉴定方法。
本发明的第一个发明目的是提供一种鉴定火炬松个体的种质或父系祖先的方法。
本发明的第二个发明目的是提供火炬松叶绿体基因组的分子标记在鉴定火炬松个体的父本来源中的应用。
本发明的第三个发明目的是提供一个鉴定火炬松个体的种质或父系祖先的分子标记组合。
本发明的第四个发明目的是提供所述分子标记组合在鉴定火炬松个体的种质或父系祖先中的应用。
本发明的第五个发明目的是提供一个鉴定火炬松个体的种质或父系祖先的引物组合。
本发明的第六个发明目的是提供所述引物组合在鉴定火炬松的种质或父系祖先、或检测所述分子标记组合，或鉴定火炬松的种质或父系祖先试剂盒中的应用。
本发明的第七个发明目的是提供一个鉴定火炬松的种质或父系祖先的试剂盒。
为了实现上述目的，本发明是通过以下技术方案予以实现的：
一种鉴定火炬松个体的种质或父系祖先的方法，检测待测样本的叶绿体基因组的分子标记。
优选地，所述种质或父系祖先为“中国森林植物种质资源信息系统”中“种质中文名”为“火炬松武32”、“火炬松p043”、“火炬松p040”、“火炬松p051”、“火炬松w26”、“火炬松w16”、“火炬松w11”、“火炬松q6”、“火炬松q13”、“火炬松n3”、“火炬松g16”、“火炬松g10”、“火炬松g01”、“火炬松a279”、“火炬松a270”、“火炬松a262”、“火炬松a259”、“火炬松a201”、“火炬松24”、“火炬松18”或“火炬松13”中的一种或几种，所述分子标记为叶绿体基因组上的21个分子标记中的一个或几个，所述21个分子标记为genbank登录号nc_021440.1的序列上的第29448、29449、37238、37239、37240、39801、97585、100523、100730、100732、101076、101085、101103、101112、101121、101130、101139、101157、101166、104944和105020个碱基。
“中国森林植物种质资源信息系统”的网址为：www.fgr.cn。
优选的，所述21个分子标记为genbank登录号nc_021440.1的序列上的第29448个碱基的t/g突变、第29449个碱基的t/g突变、第37238个碱基的t/a突变、第37239个碱基的t/a突变、第37240个碱基的t/a突变、第39801个碱基的a/g突变、第97585个碱基的t/g突变、第100523个碱基的g/t突变、第100730个碱基的a/g突变、第100732个碱基的g/t突变、第101076个碱基的t/g突变、第101085个碱基的t/g突变、第101103个碱基的t/g突变、第101112个碱基的t/g突变、第101121个碱基的t/g突变、第101130个碱基的t/g突变、第101139个碱基的g/t突变、第101157个碱基的t/g突变、第101166个碱基的t/g突变、第104944个碱基的g/a突变和第105020个碱基的t/c突变。
优选地，按照genbank登录号nc_021440.1的序列上的第29448、29449、37238、37239、37240、39801、97585、100523、100730、100732、101076、101085、101103、101112、101121、101130、101139、101157、101166、104944和105020个碱基的顺序，
以“火炬松武32”作为父系祖先的后代或“火炬松武32”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、t、a、g、t、t、t、t、t、t、g、t、t、g和t；
以“火炬松p043”作为父系祖先的后代或“火炬松p043”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、a、g、t、t、t、t、t、g、g、t、t、g和t；
以“火炬松p040”作为父系祖先的后代或“火炬松p040”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、a、g、t、t、t、t、t、t、t、t、t、g和t；
以“火炬松p051”作为父系祖先的后代或“火炬松p051”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、t、a、g、t、t、g、t、t、t、g、t、t、g和t；
以“火炬松w26”作为父系祖先的后代或“火炬松w26”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、a、g、t、t、g、t、t、t、g、t、g、g和t；
以“火炬松w16”作为父系祖先的后代或“火炬松w16”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、a、g、t、t、g、t、g、t、t、t、t、g和t；
以“火炬松w11”作为父系祖先的后代或“火炬松w11”的21个分子标记的碱基依次为：t、t、t、t、t、a、g、g、a、g、t、t、t、t、t、t、t、t、t、a和t；
以“火炬松q6”作为父系祖先的后代或“火炬松q6”的21个分子标记的碱基依次为：t、t、a、a、a、a、t、g、a、g、t、t、t、t、t、t、g、t、t、g和c；
以“火炬松q13”作为父系祖先的后代或“火炬松q13”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、a、g、t、t、t、t、t、g、t、t、t、g和t；
以“火炬松n3”作为父系祖先的后代或“火炬松n3”的21个分子标记的碱基依次为：t、t、a、a、a、a、t、t、a、g、t、t、t、t、t、t、g、t、t、g和c；
以“火炬松g16”作为父系祖先的后代或“火炬松g16”的21个分子标记的碱基依次为：t、t、a、a、a、a、t、g、a、g、t、t、t、t、t、t、g、t、t、g和t；
以“火炬松g10”作为父系祖先的后代或“火炬松g10”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、a、g、t、t、g、t、t、t、g、t、t、g和t；
以“火炬松g01”作为父系祖先的后代或“火炬松g01”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、t、g、t、t、t、t、t、g、t、g、t、t、g和t；
以“火炬松a279”作为父系祖先的后代或“火炬松a279”的21个分子标记的碱基依次为：g、t、t、t、t、a、t、g、a、g、t、t、t、t、t、t、t、t、t、a和t；
以“火炬松a270”作为父系祖先的后代或“火炬松a270”的21个分子标记的碱基依次为：t、t、t、t、t、g、t、g、a、g、t、t、t、t、t、t、g、g、t、g和t；
以“火炬松a262”作为父系祖先的后代或“火炬松a262”的21个分子标记的碱基依次为：t、t、a、a、a、a、t、g、a、g、t、t、t、t、t、t、t、g、t、g和c；
以“火炬松a259”作为父系祖先的后代或“火炬松a259”的21个分子标记的碱基依次为：t、t、a、a、a、a、t、g、a、g、t、t、g、t、t、t、t、t、t、g和t；
以“火炬松a201”作为父系祖先的后代或“火炬松a201”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、g、t、t、g、t、t、t、t、g、t、t、g和t；
以“火炬松24”作为父系祖先的后代或“火炬松24”的21个分子标记的碱基依次为：t、t、t、t、t、a、t、g、a、g、g、g、g、g、t、g、g、t、t、g和t；
以“火炬松18”作为父系祖先的后代或“火炬松18”的21个分子标记的碱基依次为：t、g、t、t、t、a、t、t、a、g、t、t、t、t、t、t、t、t、t、g和t；
以“火炬松13”作为父系祖先的后代或“火炬松13”的21个分子标记的碱基依次为：t、t、a、a、a、a、t、g、a、g、t、t、t、g、g、t、g、t、t、g和t；
火炬松叶绿体基因组的分子标记在鉴定火炬松个体的父本来源中的应用也属于本发明的保护范围、
本发明还要求保护一个鉴定火炬松个体的种质或父系祖先的分子标记组合，包括叶绿体基因组上的21个分子标记中的一个或几个，所述21个分子标记为genbank登录号nc_021440.1的序列上的第29448、29449、37238、37239、37240、39801、97585、100523、100730、100732、101076、101085、101103、101112、101121、101130、101139、101157、101166、104944和105020个碱基，所述种质或父系祖先为“中国森林植物种质资源信息系统”中“种质中文名”为“火炬松武32”、“火炬松p043”、“火炬松p040”、“火炬松p051”、“火炬松w26”、“火炬松w16”、“火炬松w11”、“火炬松q6”、“火炬松q13”、“火炬松n3”、“火炬松g16”、“火炬松g10”、“火炬松g01”、“火炬松a279”、“火炬松a270”、“火炬松a262”、“火炬松a259”、“火炬松a201”、“火炬松24”、“火炬松18”或“火炬松13”中的一种或几种。
以上所述分子标记组合在鉴定火炬松个体的种质或父系祖先中的应用，所述种质或父系祖先为“中国森林植物种质资源信息系统”中“种质中文名”为“火炬松武32”、“火炬松p043”、“火炬松p040”、“火炬松p051”、“火炬松w26”、“火炬松w16”、“火炬松w11”、“火炬松q6”、“火炬松q13”、“火炬松n3”、“火炬松g16”、“火炬松g10”、“火炬松g01”、“火炬松a279”、“火炬松a270”、“火炬松a262”、“火炬松a259”、“火炬松a201”、“火炬松24”、“火炬松18”或“火炬松13”中的一种或几种。
一个鉴定火炬松个体的种质或父系祖先的引物组合，能够扩增以上所述分子标记组合，优选地，所述能够扩增以上所述分子标记组合包括核苷酸序列如seqnoid：1～20所述引物一条或几条。
以上所述引物组合在鉴定火炬松的种质或父系祖先、或检测以上所述分子标记组合，或鉴定火炬松的种质或父系祖先试剂盒中的应用，所述种质或父系祖先为“中国森林植物种质资源信息系统”中“种质中文名”为“火炬松武32”、“火炬松p043”、“火炬松p040”、“火炬松p051”、“火炬松w26”、“火炬松w16”、“火炬松w11”、“火炬松q6”、“火炬松q13”、“火炬松n3”、“火炬松g16”、“火炬松g10”、“火炬松g01”、“火炬松a279”、“火炬松a270”、“火炬松a262”、“火炬松a259”、“火炬松a201”、“火炬松24”、“火炬松18”或“火炬松13”中的一种或几种。
一个鉴定火炬松的种质或父系祖先的试剂盒，包括所述引物组合，所述种质或父系祖先为“中国森林植物种质资源信息系统”中“种质中文名”为“火炬松武32”、“火炬松p043”、“火炬松p040”、“火炬松p051”、“火炬松w26”、“火炬松w16”、“火炬松w11”、“火炬松q6”、“火炬松q13”、“火炬松n3”、“火炬松g16”、“火炬松g10”、“火炬松g01”、“火炬松a279”、“火炬松a270”、“火炬松a262”、“火炬松a259”、“火炬松a201”、“火炬松24”、“火炬松18”或“火炬松13”中的一种或几种。
一种叶绿体基因组的分子标记的筛选鉴定火炬松个体的种质或父系祖先的试剂盒：
一、组成
核苷酸序列如seqnoid：1～20所述引物，2×taqmastermix，ddh2o。
二、使用方法
(1)提取待测样本叶绿体dna；
(2)使用核苷酸序列如seqnoid：1～20所述引物对样本叶绿体dna进行pcr扩增，
其中，pcr反应体系如下：
扩增突变位点的pcr程序：94℃预变性4分钟，94℃变性20秒，56℃退火30秒，72℃延伸3分钟，变性、退火、延伸三个步骤循环35次，最后再72℃充分延伸10分钟。
(3)pcr产物检测
利用2％琼脂糖凝胶电泳检测pcr产物，根据判断是否为目的条带。
(4)sanger测序
采用abi3730xl测序仪对pcr产物进行测序。
与现有技术相比，本发明具有如下有益效果：
利用存在于21个火炬松种质间的21个叶绿体基因组snp特异性位点，实现了对不同个体的区分，并能够以此鉴定其父系祖先。本发明首先发现火炬松叶绿体基因组序列的多态性，发现了叶绿体基因组可以作为种群亲缘关系鉴定。通过对叶绿体基因组序列差异区段的比较，获得能够用于筛选鉴定火炬松个体的种质或父系祖先的鉴定的分子标记，并进一步获得检测分子标记的特异性pcr扩增引物，从而实现了利用pcr扩增技术对群体进行快速鉴定，在火炬松育种领域会发挥中大的作用。
具体实施方式
下面结合说明书和具体实施例对本发明作出进一步地详细阐述，所述实施例只用于解释本发明，并非用于限定本发明的范围。下述实施例中所使用的试验方法如无特殊说明，均为常规方法；所使用的材料、试剂等，如无特殊说明，为可从商业途径得到的试剂和材料。
实施例1叶绿体基因组的分子标记的筛选
一、实验方法
(1)叶绿体dna的提取
1、叶绿体的提取
①取松针约10g洗净擦干于4℃冰箱中黑暗处理24～48h。
②将黑暗处理过的针叶剪成大概1cm左右置入匀浆机，加入6倍体积预冷(4℃)buffera，使叶片完全被浸没。低速下匀浆4次，每次10秒左右，然后中速匀浆2次，每次10秒左右，若匀浆不充分则继续重复以上步骤，待匀浆充分后用纱布过滤浆液至一洁净烧杯中，收集滤液，弃滤渣；将滤液分装到50ml离心管中，4℃下2500rpm离心6min，弃上清。
③向沉淀中加30mlbufferb，充分悬浮沉淀，4℃下3000rpm离心6min，弃上清。
④向沉淀中加30mlbufferc，充分悬浮沉淀，4℃下3000rpm离心6min，弃上清。
⑤向沉淀中加15mlbufferd，充分悬浮沉淀，4℃下3000rpm离心6min，弃上清；得到的沉淀为叶绿体粗沉淀，用干净的牙签沾取微量沉淀，可在olympusbx43f荧光显微镜下观察叶绿体的形态与完整性。
2、叶绿体dna的提取
①向上述沉淀中加入20ml已预热至65℃的2×ctab提取缓冲液，于65℃恒温水浴锅中裂解细胞40min，期间上下颠倒离心管若干次。
②加入等体积的酚/氯仿/异戊醇(25:24:1)抽提，轻轻颠倒若干次充分混匀，4℃下10000rpm离心10min。
③把上清移入新的离心管中，然后加入等体积的氯仿继续抽提，轻轻颠倒若干次充分混匀,4℃下10000rpm离心10min。
④取上清，加入0.7倍体积的异丙醇和0.1倍体积的3mol/lnaac，颠倒若干次充分混匀，放入-20℃冰箱中静置过夜沉淀dna。
⑤4℃下10000rpm离心30min，弃上清，此时的沉淀为dna粗沉淀，用70％乙醇洗涤dna沉淀2次，无水乙醇洗涤一次，自然风干，加入400μlte溶解dna。
⑥加入20μlrnasei于dna溶解液中，轻轻混匀后在37℃金属浴lh酶解rna，再加入4μl蛋白酶k于37℃金属浴30min。
⑦利用天根dp320试剂盒纯化所得dna溶解液。
3、相关试剂的配置：
buffera：50mol/ltris,25mmol/ledta,1.25mol/lnacl,0.25mol/lbsa,ph＝3.6；高温高压灭菌，4℃保存。
bufferb:50mol/ltris,25mmol/ledta,1.25mol/lnacl,1mmol/ldtt,0.1％牛血清蛋白(使用前加)，ph＝8.0；高温高压灭菌，4℃保存。
bufferc:150mmol/lnacl,100mmol/ledta,ph＝8.0；高温高压灭菌，4℃保存。
bufferd:50mol/ltris,25mmol/ledta,ph＝8.0；高温高压灭菌，4℃保存。
1mol/l二硫苏糖醇(dtt)溶液：用20ml0.01mol/l乙酸钠溶液(ph5.2)溶解3.09gdtt，经过滤除菌后分装成小份于-20℃保存。
1％牛血清蛋白(100ml)：称取1gbsa，加水定容至100ml，过滤除菌，4℃保存。
2×ctab(200ml)：ctab4g,nacl16.364g,1mtris-hcl(ph＝8.0)20ml,0.5medta(ph＝8.0)8ml,加水定容至200ml；高温高压灭菌，室温保存。
5×tbe电泳缓冲液：0.45mol/ltris-硼酸，0.01mol/ledta；4℃保存。
6×loadingbuffer：0.25％(w/v)溴粉蓝，40％(w/v)蔗糖水溶液；4℃保存。
酚：氯仿：异戊醇(25:24:1)；避免光照氧化，现配现用。
3mnaac(500ml)：123.04g无水乙酸钠，醋酸调ph至5.2，定容至500ml；高温高压灭菌，4℃保存。
(二)高通量测序
1、文库构建
将在“中国森林植物种质资源信息系统”中“种质中文名”为“火炬松武32”、“火炬松p043”、“火炬松p040”、“火炬松p051”、“火炬松w26”、“火炬松w16”、“火炬松w11”、“火炬松q6”、“火炬松q13”、“火炬松n3”、“火炬松g16”、“火炬松g10”、“火炬松g01”、“火炬松a279”、“火炬松a270”、“火炬松a262”、“火炬松a259”、“火炬松a201”、“火炬松24”、“火炬松18”和“火炬松13”的21个火炬松种质的dna溶液根据dna量等量混合，混合后的dna溶液使用dna超声片段化技术，将检验合格的dna序列打断成约500bp大小的片段。经末端修复后,在修复产物3’末端加上a碱基，再连接上adaptor测序接头，回收纯化连接有adaptor接头的目的片段，pcr扩增后完成文库构建。经检验合格的文库采用hiseq2500/hiseq4000/miseq/novaseq/xten等高通量测序平台进行高通量测序。
2、高通量测序
dna测序采用illumina的hiseq2500/hiseq4000/miseq/novaseq/xten等高通量测序平台进行，为实现高质量的基因组序列测序和组装，对已构建的双末端文库进行pairend双末端测序，测序读长为双端150bp。要求测序q20(质量值在20以上的碱基所占的比例，即碱基测序错误率在1％以下)指标达到90％以上，q30(质量值在30以上的碱基所占的比例，即碱基测序错误率在0.1％以下)指标达到85％以上。
3、质量检测
为保证后续分析的可靠性，对测序数据进行必要的检测和数据筛选，检测的内容包括：测序数据量、测序数据质量、gc比例分布、测序正确率等；同时基于检测结果，对数据进行筛选，筛选的标准为：去除含有n或含有3个以上碱基质量值低于20的测序序列。
4、序列比对与变异检测
序列比对分析采用短片段测序比对软件bwa(版本：0.7.12-r1039)，将经过质量筛选后的双末端测序reads比对到火炬松叶绿体基因组(ncbi登录号：nc_021440.1)，比对采用align比对方法，比对采用默认参数，比对文件产生后，使用bwa软件的sampe方法定位测序数据在叶绿体基因组上的位置。根据测序数据比对结果，采用samtools的sort分析方法，对比对结果进行排序整理。并使用mpileup方法，分析检测序列比对结果中，叶绿体基因组中的单碱基突变(snp)，短片段插入缺失突变(indel)和简单串联重复突变(ssr)。
(三)高频突变区筛选
高频突变位点的筛选根据snp和indel位点在叶绿体基因组中的位置分布，以在一对pcr引物中扩增出多个snp或indel为筛选方向，以实现在尽量少的pcr反应中扩增出尽量多的分子标记为目标，从54个突变位点中选择21个突变位点用于分子标记。
二、实验结果
通过使用mpileup方法，初步得到54个突变位点(表1)。
表1:
(注：corrdiante：突变所在位置；ref：参考基因组的碱基类型；snv/indels：突变的类型，包括snv和indel；rate：突变类型所占的比例；gene：所在的基因区或基因间隔区。)
经过进一步筛选，得到21个在这21个火炬松种质中特异性好的的突变位点。使用pcr引物检测了这21个火炬松种植的样品，共检测出21个snp差异位点，检测出的位点与高通量测序检测位点一致，根据差异位点显示，可以区分不同的个体。
具体基因型如表2所示。
表2：
(pos表示位点所在叶绿体基因组的位置；snp为该位点参考)
实施例2对特定snp位点设计引物
一、实验方法
根据表2的突变位点的分布信息，对检测到snp位点的叶绿体基因组序列区段设为候选区段，对候选区段进行pcr引物设计。
二、实验结果
根据引物扩增的特异性和扩增效率，得到扩增表2的突变位点的引物，引物见表3。
表3：
实施例3一种叶绿体基因组的分子标记的筛选鉴定火炬松个体的种质或父系祖先的试剂盒
一、组成
核苷酸序列如seqnoid：1～20所述引物，2×taqmastermix，ddh2o。
二、使用方法
(1)提取待测样本叶绿体dna；
(2)使用核苷酸序列如seqnoid：1～20所述引物对样本叶绿体dna进行pcr扩增，扩增仪：苏州东胜兴业科学仪器有限公司etc811pcr扩增仪，
其中，pcr反应体系如下：
扩增突变位点的pcr程序：94℃预变性4分钟，94℃变性20秒，56℃退火30秒，72℃延伸3分钟，变性、退火、延伸三个步骤循环35次，最后再72℃充分延伸10分钟。
(3)pcr产物检测
利用2％琼脂糖凝胶电泳检测pcr产物，条带与marker对比，根据大小初步判断是否为目的条带。
(4)sanger测序
采用abi3730xl测序仪对pcr产物进行测序。测序引物为pcr扩增所用引物。
实施例4鉴定火炬松个体的种质
一、实验方法
按照实施例3中的试剂盒检测22个火炬松样品，以illumina高通量测序作为对照，检测表2中的突变位点。
其中，所有样本的父系祖先均属于在“中国森林植物种质资源信息系统”中“种质中文名”为“火炬松武32”、“火炬松p043”、“火炬松p040”、“火炬松p051”、“火炬松w26”、“火炬松w16”、“火炬松w11”、“火炬松q6”、“火炬松q13”、“火炬松n3”、“火炬松g16”、“火炬松g10”、“火炬松g01”、“火炬松a279”、“火炬松a270”、“火炬松a262”、“火炬松a259”、“火炬松a201”、“火炬松24”、“火炬松18”或“火炬松13”的22个火炬松种质。
二.实验结果
结果如表4所示。
结果显示检测到突变位点与高通量测序检测位点一致，因此本方法可靠、准确。