本申请涉及谷子育种技术领域，特别是涉及一种与谷子码数性状相关的snp标记及其检测引物和应用。
背景技术：
我国是谷子的原产地和世界上栽培面积最大、产量最高的国家，产量约占全世界总量的80％。同时，我国也是谷子遗传资源数量最多、多样性最丰富的国家。
谷穗由穗轴和众多谷码组成，码数和码粒数等是影响谷子产量的重要因子，研究与控制谷子码数、码粒数等产量性状相关的基因及其位置功能等，对于指导谷子遗传育种具有重要的意义。
单核苷酸多态性(singlenucleotidepolymorphism，snp)是基因组水平上由单个核苷酸碱基的变异引起的dna序列多态性，其在基因组中数量多，分布广泛，遗传稳定性好。随着基因测序技术的发展和成本的降低，对不同个体同一基因或基因片段进行直接测序和序列比较，就能够确定碱基是否存在变异。因此，snp检测有利于基因分型，适用于快速和规模化地筛查未知或已知snp与某遗传性状的关系。
目前关于谷子数量性状相关的qtl定位和snp标记研究主要集中于ssr标记的开发利用等途经，而利用群体自交系简并基因组测序检测多态性的snp分子标记的开发应用研究尚未有报道。因此，开展谷子数量性状snp分子标记的开发，并建立辅助选育体系，对于提高谷子产量，节约育种成本具有重要意义。
技术实现要素：
本申请的目的是提供一种新的与谷子码数性状相关的snp标记及其检测引物和应用。
为了实现上述目的，本申请采用了以下技术方案：
本申请的一方面公开了一种与谷子码数性状相关的snp标记，该snp标记位于第二条染色体第27578869bp至第28178402bp的bin标记内，并且，snp标记与谷子码数紧密连锁。
需要说明的是，本申请经过研究发现了一段与谷子码数紧密连锁的序列，这段序列内的snp标记直接影响谷子码数性状，这段序列即第二条染色体第27578869bp至第28178402bp的bin标记。
进一步的研究发现，第二条染色体第27578869bp至第28178402bp的bin标记中，尤其是snp标记的碱基为t或c，与谷子码数紧密连锁。
本申请的一种实现方式中，具体披露了其中一个snp标记，即cn-02-349位点，cn-02-349位点位于第二条染色体的序列第28175349bp位置，cn-02-349位点所在的序列如seqidno.1所示，cn-02-349位点是seqidno.1所示序列中自5’端起第66位碱基。
需要说明的是，第二条染色体第27578869bp至第28178402bp的bin标记与谷子码数紧密连锁，其中可能包含众多snp标记，cn-02-349位点只是本申请的一种实现方式中证实的与谷子码数性状相关的snp标记，在本申请的指导下，不排除在该bin标记内还有其它的snp标记。
优选的，本申请与谷子码数性状相关的snp标记，对snp标记，利用复合区间作图法，采用5％的总体显著水平，qtl检测的lod值为5.5675，表型变异解释率为4.95％。
本申请的另一面公开了一种检测本申请的snp标记的引物对，该引物对的上游引物为seqidno.2所示序列，下游引物为seqidno.3所示序列；
seqidno.2：5’-cgcaccaaactagccaccg-3’
seqidno.3：5’-caatgagcgagcagggagg-3’。
本申请的另一面公开了一种检测本申请的snp标记的试剂盒，该试剂盒中包括本申请的引物对。
优选的，本申请的试剂盒中还包括用于pcr扩增的试剂。
本申请的另一面公开了一种检测本申请的snp标记的方法，包括采用本申请的引物对，或者本申请的试剂盒，对待测的谷子基因组dna进行pcr扩增，通过pcr扩增产物的测序结果，获得snp标记的碱基情况。
本申请的另一面公开了一种预测谷子码数的方法，包括采用本申请的引物对，或者本申请的试剂盒，对预测对象的谷子基因组dna进行pcr扩增，通过pcr扩增产物的测序结果，获得本申请的snp标记的碱基情况，根据snp标记的碱基情况预测谷子码数。
本申请的再一面公开了本申请的snp标记在谷子码数预测、谷子码数早期鉴定或谷子分子辅助筛选育种中的应用。
本申请的再一面公开了本申请的引物对或本申请的试剂盒在谷子码数预测、谷子码数早期鉴定或谷子分子辅助筛选育种中的应用。
由于采用以上技术方案，本申请的有益效果在于：
本申请与谷子码数性状相关的snp标记，可以用于谷子码数性状的分子标记辅助育种，通过该snp标记检测可以预测谷子码数，为实现谷子码数性状的早期鉴定或筛选育种提供了科学依据，对加快谷子品种的遗传选育或改良进程具有重要的理论和实践指导意义。本申请检测snp标记的引物对，可以对谷子码数进行分型，检测snp表达的差异。
附图说明
图1是本申请实施例中谷子父母本基因组dna电泳图，其中，第一泳道为m泳道，表示dl2000marker，第二泳道即标注为1的泳道表示父本“冀张谷一号”基因组dna，第三泳道即标注为2的泳道表示母本“a2”基因组dna；
图2是本申请实施例中谷子父母本基因组dna采用引物对21_2f和21_2r进行扩增的pcr产物电泳图，其中，第一泳道即标注a21-2的泳道，表示以母本“a2”基因组dna为模板的pcr产物，第二泳道即标注为321-2的泳道，表示以父本“冀张谷一号”基因组dna为模板的pcr产物，第三泳道为m泳道，表示dl2000marker；
图3是本申请实施例中母本和父本snp标记中cn-02-349位点所在的序列的比对结果图，cn-02-349位点在第二条染色体第27578869bp至第28178402bp的bin标记中，母本碱基为t，父本碱基为c。
具体实施方式
本申请提供了与谷子码数相关的bin标记区域的snp标记、检测引物及其应用。谷子码数紧密连锁snp标记cn-02-349，位于第二条染色体第27578869bp至第28178402bp的bin标记内。
需要说明的是，本申请中bin是指bin作图(即binmap)中，利用作图群体中全部的重组位点信息，获得构成重组事件的最小单位(recombinationbin)。利用得到的bin作为标记用于后续的连锁图谱构建和qtl定位。
本申请中，遗传图谱(geneticmap或linkagemap)是指以遗传标记间重组率为基础构建的标记之间相对位置的线性排列图。基于高通量重测序技术通过对作图群体亲本和子代群体进行测序，将一定数量的连续snp作为判断子代重组位点的依据，得到每个子代的全基因组物理重组图谱。利用作图群体中全部的重组位点信息，获得构成重组事件的最小单位(recombinationbin)和bin图(即binmap)，并利用得到的bin作为标记用于后续的连锁图谱构建和qtl定位。
本申请提供谷子父本“冀张谷一号”和母本“a2”材料，二者杂交的f1代材料，及自交13代的f2群体441份材料。将各份材料在河北张家口种植，并记录和整理码数等性状数据。将各份材料进行简并基因组重测序后，比对参考基因组拼接完成，获得snp多态性标记。通过码数性状数据，获得与之相关的snp标记，并通过克隆测序验证snp标记。
下面通过具体实施例和附图对本申请作进一步详细说明。以下实施例仅对本申请进行进一步说明，不应理解为对本申请的限制。
实施例一高通量测序和snp标记信息分析
本例利用谷子父本“冀张谷一号”和母本“a2”材料，二者杂交的f1代材料，及自交13代的f2群体441份材料，即重组自交系(recombinantinbredlines，缩写rils)。将各份材料在河北张家口种植，并记录和整理码数等性状数据。将各份材料进行简并基因组重测序后，比对参考基因组拼接完成，获得snp标记。通过分析码数性状数据，获得与之相关的snp标记，并通过克隆测序验证snp标记。
本例利用radseq法，提取每个样本个体的基因组dna，共计444个样本，其中包括其中441个rils、2个亲本、1个f1个体。对提取的基因组dna进行质量检测，对合格的dna进行酶切，电泳回收dna片段，并加上接头进行cluster制备，最后上机测序。
本例的具体步骤如下：
(1)称取1.0g新鲜叶片，剪碎放入研钵，用液氮研磨后加入3ml1.5×ctab，研磨成匀浆转入15ml的离心管中，然后往研钵中加入1ml1.5×ctab冲洗，再转入离心管中，混匀后于65℃水浴30min，期间不时缓慢摇匀。
其中1l的1.5×ctab配方包括：15g的ctab、1mol/l的ph8.0的tris.cl75ml、0.5mol/l的edta30ml、nacl61.4g，加去离子水定容至1l；使用前加入约2ml的巯基乙醇，使其终浓度为0.2％。
(2)水浴完成后，冷却至室温，加入等体积氯仿/异戊醇，轻轻混匀，至下层液变为深绿色；其中氯仿/异戊醇的体积比为氯仿:异戊醇为24:1。
(3)4200rpm离心10min，将上层水相移到新的15ml离心管，加2倍体积预冷的无水乙醇，混合静止5min。于-20℃放置30min沉淀dna。
(4)4200rpm离心10min，弃掉上清，加入1ml75％乙醇洗涤沉淀1次，倒置离心管干燥dna，加入50μlte溶解dna。
(5)检测dna的浓度，并用水调整为20ng/μl。
(6)采用psti酶进行酶切，打断基因组dna，psti酶的酶切反应体系为30μl，包括：基因组dna20μl、psti酶1μl、10×psti缓冲液3μl、nuclease-freewater6μl。酶切反应体系配制好后，在37℃恒温孵育至少1小时。
(7)加接头反应，在打断的基因组dna片段上添加测序接头，加接头反应的体系为40μl，包括：psti酶切产物30μl、测序接头1μl、10×连接buffer1μl、t4dnaligase2μl、ratp0.4μl、ddh2o5.6μl。
(8)每个样品各取1μl反应产物，加入在一个新的离心管中，每12个样品一组，每组的总体积为12μl。
(9)用3％回收胶电泳1h，eb染色后切取300-700bp大小片段。用qiaquickkit进行胶纯化回收，回收产物溶于30μl的eb溶液。
(10)进行pcr反应，采用测序接头的通用引物对切记回收的片段进行pcr扩增，pcr反应体系25μl，包括：10×buffer2.5μl、25mm的dntps0.25μl、5u/μl的高保真酶0.25μl、10μm的测序接头正向引物0.5μl、10μm的测序接头反向引物0.5μl、切胶回收产物20μl、无菌水1μl。
pcr反应条件为，94℃预变性5min，然后进入10个循环：94℃变性30秒、60℃退火30秒、72℃延伸40秒，循环结束后72℃延伸3min。
(11)磁珠纯化，完成建库。
具体纯化方法为：首先，向pcr产物中加入1.2倍体积磁珠，静置10min。然后，置于磁力架上吸附，去除上清。然后，加入500μl70％酒精洗涤两次。干热仪上蒸干后，加入15μleb溶液溶解5min。最后，磁力架上吸附，转移上清于1.5ml离心管中。
(12)检测达到上机标准后上机测序。
结果：对444个样本进行酶切建库测序，得到75.99gb原始数据，平均每个个体171.54mb。将测序序列比对到参考基因组上，平均比对率86.326％，平均覆盖率8.226％，平均测序深度为3.655×。
其中，444个样本包括441个rils、2个亲本、1个f1个体；本例的参考基因组即yugu基因组，可以从下列网址获取yugu基因组序列：
https://www.ncbi.nlm.nih.gov/genome/？term＝setaria+italica+(foxtail+millet)。
通过测序和信息分析，获得亲本间有多态的33771个snp标记。根据窗口滑动的方法，选取若干个snp为一个窗口，每次滑动一个snp确定每个窗口的基因型以及每个个体的交换位点，生成bin基因型和bin图。根据bin基因型数据，用mstmap软件构建遗传图谱，2022个bin定位到9条染色体上，再将遗传图谱数据导入到mapchart软件，整合一张基因组遗传连锁图谱。
利用所构建的谷子高密度遗传图谱，采用复合区间作图分析(cim)，对谷子码数性状表型进行qtl分析。qtl检测采用5％的总体显著水平，根据500次排列检验结果，确定码数性状表型数据qtl分析的临界lod值，分析得到与码数相关的预测bin区间和snp位点信息。
结果显示，与谷子码数性状相关的snp标记位于遗传连锁图谱第二条染色体的194.66厘摩(cm)，第153个bin至第154个bin(即chr2_bin153-chr2_bin154)，也就是第二条染色体第27578869bp至第28178402bp的bin标记内；其中，cn-02-349位点位于第二条染色体的序列第28175349bp位置；这些与谷子码数性状相关的snp标记的碱基都是t或c。利用复合区间作图法，采用5％的总体显著水平，qtl检测的lod值为5.5675，表型变异解释率为4.95％，加性效应值(additiveeffect,a)为-3.9644。
实施例二snp标记验证
根据预测的bin标记和snp位点，对比参考基因组，得到相关区域基因序列，选取snp位点前后300bp左右，设计开发snp标记引物，以父本和母本材料dna为模板进行pcr扩增。选择引物扩增正常、pcr扩增产物复合预测大小的扩增产物进行切胶回收，并进行测序，选择父本和母本材料扩增基因序列有snp位点差异的标记引物。以f1代和f2代材料为模版，进行snp标记引物pcr扩增和结果测序，选择符合基因分型和与表型性状相关的标记引物。
首先，根据pcr产物回收后的测序结果，选择父本和母本材料扩增基因序列snp位点差异的标记。具体步骤如下：
(1)按照实施例一中步骤(1)至(4)，用ctab法分别提取父本和母本基因组dna。
(2)用0.8％的琼脂糖凝胶电泳检测提取的父本和母本基因组dna，检测结果如图1所示，其中，第二泳道即标注为1的泳道表示父本“冀张谷一号”基因组dna，第三泳道即标注为2的泳道表示母本“a2”基因组dna。可见父本和母本基因组dna提取成功。
(3)将提取的父本和母本基因组dna存储于-20℃备用。
(4)分别以提取的父本和母本基因组dna为模板，采用引物对21_2f和21_2r对基因组dna进行pcr扩增；引物对的上游引物21_2f为seqidno.2所示序列，下游引物21_2r为seqidno.3所示序列；
seqidno.2：5’-cgcaccaaactagccaccg-3’
seqidno.3：5’-caatgagcgagcagggagg-3’。
pcr反应体系25μl，包括：10×buffer2.5μl、25mm的dntps0.15μl、5u/μl的taq酶0.15μl、10μm的上游引物0.5μl、10μm的下游引物0.5μl、基因组dna1.0μl、无菌水20.2μl。
pcr反应条件为，94℃预变性5min，然后进入35个循环：94℃变性30秒、60℃退火30秒、72℃延伸40秒，循环结束后72℃延伸3min。
pcr扩增产物经纯化后于4℃保存。各pcr扩增产物取部分进行1％琼脂糖凝胶电泳检测，结果如图2所示。图2是引物对21_2f和21_2r进行扩增的pcr产物电泳图。图2的结果显示，设计的引物对能够对父本和母本的基因组dna扩增出符合预期的靶标片段。
本例分别对引物对21_2f和21_2r扩增的父本和母本基因组dna的扩增产物进行测序，引物21_2f和21_2r扩增得到的母本扩增产物的测序序列如seqidno.1所示序列，cn-02-349位点是seqidno.1所示序列中自5’端起第66位碱基。父本扩增产物和母本扩增产物的测序序列比对结果如图3所示，其中背景标灰的碱基表示cn-02-349位点的碱基，母本碱基为t，父本碱基为c。
然后，根据样品码数数据筛选：父本码数为80，母本码数为67。441份样品中，最小值20，最大值126；码数最大的50个样品中，与父本snp相同的样品数多于60％，码数最小的50个样品中，与母本snp相同的样品数多于80％。
以上内容是结合具体的实施方式对本申请所作的进一步详细说明，不能认定本申请的具体实施只局限于这些说明。对于本申请所属技术领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干简单推演或替换。
sequencelisting
<110>张家口市农业科学院(河北省高寒作物研究所)
<120>与谷子码数性状相关的snp标记及其检测引物和应用
<130>17i25623
<160>3
<170>patentinversion3.3
<210>1
<211>240
<212>dna
<213>cn-02-349位点所在序列
<400>1
cgcaccaaactagccaccggtttgtcgatcgatctctgccgccgctactccatcaatctg60
atgtgttgattggaatgtcttggactgatcaactaatcattaagttgagcgcgcggcggc120
acttaactcgccgcacgcatgccttgtggtgacgcctgcggcatgcatatatgcatgggc180
cacacggcacggcggcggcggccaatggccatggtccacgaccacgacccccatcaacca240
<210>2
<211>19
<212>dna
<213>人工序列
<400>2
cgcaccaaactagccaccg19
<210>3
<211>19
<212>dna
<213>人工序列
<400>3
caatgagcgagcagggagg19