与谷子穗长性状相关的分子标记及其检测引物和应用的制作方法

本发明涉及谷子育种
技术领域：
，尤其涉及与谷子穗长性状相关的分子标记及其检测引物和应用。
背景技术：
：我国是谷子的原产地和世界上栽培面积最大、产量最高的国家，产量约占全世界总量的80％。同时，我国也是谷子遗传资源数量最多、多样性最丰富的国家。穗长、株高等是影响谷子产量的重要因子，研究与控制谷子穗长等产量性状相关的基因位置及其功能等，对于指导谷子遗传育种具有重要的意义。单核苷酸多态性(singlenucleotidepolymorphism，snp)是基因组水平上由单个核苷酸碱基的变异引起的dna序列多态性，其在基因组中数量多，分布广泛，遗传稳定性好。随着基因测序技术的发展和成本的降低，对不同个体同一基因或基因片段进行直接测序和序列比较，可以确定碱基是否存在变异。因此，snp检测有利于基因分型，适用于快速和规模化筛查未知或已知snp与某遗传性状的关系。目前关于谷子数量性状相关的qtl(quantitativetraitlocus)定位和snp标记研究主要集中于ssr标记的开发利用等途经，而利用群体自交系基因组测序检测单核苷酸多态性(snp)分子标记的开发应用研究尚未有报道。因此，开展谷子数量性状snp分子标记的开发，并建立辅助选育体系，对于提高谷子产量，节约育种成本具有重要意义。技术实现要素：本发明提供一种与谷子穗长性状相关的分子标记及其检测引物和应用，可以预测谷子穗长，为实现谷子穗长性状的早期鉴定和筛选育种提供分子辅助技术支持。根据本发明的第一方面，本发明提供一种与谷子穗长性状相关的snp标记，该snp标记是pl-02-016，其位于谷子第2号染色体的序列第38395016bp位置，碱基为c/t。进一步地，上述pl-02-016位点所在的序列如seqidno：3所示，上述pl-02-016位点是seqidno：3所示序列自5’端起第148位碱基。进一步地，对上述snp标记，利用复合区间作图法，采用5％的总体显著水平，qtl检测的lod值为3.7556，表型变异解释率为3.35％。根据本发明的第二方面，本发明提供一种用于检测如第一方面的snp标记的引物对，包括：上游引物7_2f：5’-gtatgctccacgcccttta-3’(seqidno：1)和下游引物7_2r：5’-ttgcgattaccacttgatt-3’(seqidno：2)。根据本发明的第三方面，本发明提供一种用于检测如第一方面的snp标记的试剂盒，包括如第二方面的引物对，以及任选的用于pcr扩增的试剂成分，这些试剂成分可以包括pcr缓冲液、dntps、dna聚合酶等。根据本发明的第四方面，本发明提供一种检测如第一方面的snp标记的方法，采用如第二方面的引物对，对待检测的谷子基因组dna进行pcr扩增，并通过测序扩增产物来分析pl-02-016位点的碱基情况。根据本发明的第五方面，本发明提供如第二方面的引物对在检测如第一方面的snp标记中的应用。根据本发明的第六方面，本发明提供一种谷子穗长预测方法，采用如第二方面的引物对，对谷子基因组dna进行pcr扩增，并通过测序扩增产物来分析pl-02-016位点的碱基情况，进而预测谷子穗长。根据本发明的第七方面，本发明提供如第一方面的snp标记在谷子穗长预测、或谷子穗长性状的早期鉴定、或谷子分子辅助育种中的应用。根据本发明的第八方面，本发明提供如第二方面的引物对在谷子穗长预测、或谷子穗长性状的早期鉴定、或谷子分子辅助育种中的应用。本发明的有益效果是：本发明提供的与谷子穗长性状相关的snp标记可以用于谷子穗长性状的分子标记辅助育种，通过该snp标记可以预测谷子穗长，为实现谷子穗长性状的早期鉴定和筛选育种提供分子辅助技术支持，对于加快谷子品种的遗传选育和改良进程具有重要的理论和实践指导意义。本发明的特异性引物可以对谷子穗长进行良好分型，检测snp表达的差异。附图说明图1为本发明实施例中谷子父母本基因组dna电泳胶图，其中m泳道表示dl2000marker；1泳道表示张谷三号基因组dna；2泳道表示a2基因组dna。图2为本发明实施例中谷子父母本基因组pcr产物电泳胶图，其中a7-2表示以a2dna为模板、7-2为引物的pcr产物，37-2表示以张谷三号dna为模板、7-2为引物的pcr产物，m表示dnamarker的片段大小，从下至上依次为100bp、250bp、500bp、750bp、1000bp、2000bp。图3为本发明实施例中母本和父本snp标记pl-02-016所在的序列的比对结果图，snp标记pl-02-016位于第2号染色体，第38395016bp位置，母本碱基为c，父本碱基为t。具体实施方式下面通过具体实施方式结合附图对本发明作进一步详细说明。本发明实施例提供与谷子穗长性状相关的分子标记、引物及其应用。本发明实施例的谷子穗长紧密连锁snp标记pl-02-016，位于谷子第2号染色体的序列第38395016bp位置(参考基因组——yugu基因组上的碱基编号)，碱基为c/t。值得说明的是，本发明中bin是指bin作图(binmap)中，利用作图群体中全部的重组位点信息，获得构成重组事件的最小单位(recombinationbin)。利用得到的bin作为标记用于后续的连锁图谱构建和qtl定位。在本发明实施例中，遗传图谱(geneticmap或linkagemap)是指以遗传标记间重组率为基础构建的标记之间相对位置的线性排列图。基于高通量重测序技术通过对作图群体亲本和子代群体进行测序，将一定数量的连续snp作为判断子代重组位点的依据，得到每个子代的全基因组物理重组图谱。利用作图群体中全部的重组位点信息，获得构成重组事件的最小单位(recombinationbin)和bin图(binmap)，并利用得到的bin作为标记用于后续的连锁图谱构建和qtl定位。本发明实施例提供谷子父本“张谷三号”和母本“a2”材料(用谷子不育系1066a与谷子光温敏不育系821杂交后，经多代选育而成)、二者杂交的f1代材料及自交13代的f2群体441份材料。将各份材料在河北张家口种植，并记录和整理穗长等性状数据。将各份材料进行简并基因组重测序后，比对参考基因组拼接完成，获得snp标记。通过分析穗长性状数据，获得与之相关的snp标记，并通过克隆测序验证该snp标记。以下通过实施例详细说明本发明的技术方案和技术效果，应当理解，实施例仅是示例性的，用于说明本发明的可行性以及得到的结果，不能理解为对本发明保护范围的限制。实施例1：高通量测序和snp标记信息分析本实施例利用谷子父本“张杂谷3号”和母本“a2”材料、二者杂交获得的f1代材料及f1代自交13代后得到的f2群体441份材料，即重组自交系(recombinantinbredlines，rils)。将各份材料在河北张家口种植，并记录和整理穗长等性状数据。将各份材料进行简并基因组重测序后，比对参考基因组拼接完成，获得snp标记。通过分析穗长性状数据，获得与之相关的snp标记，并通过克隆测序验证该snp标记。本实施例利用radseq方法，提取每个样本个体(其中441个rils、2个亲本、1个f1个体)的基因组dna并进行dna质量检测，对合格的dna进行酶切，电泳回收dna片段，并加上接头进行簇(cluster)制备，最后上机测序。本实施例的具体步骤如下：(1)称取1.0g新鲜叶片，剪碎放入研钵，用液氮研磨后加入3ml1.5×ctab，研磨成匀浆转入15ml的离心管中，然后往研钵中加入1ml1.5×ctab冲洗，再转入离心管中。混匀后于65℃水浴30min，期间不时缓慢摇匀。其中1.5×ctab配方(1l)如下表1：表1成分用量ctab15g1mol/ltris.cl(ph为8.0)75ml0.5mol/ledta30mlnacl61.4g加去离子水定容至1l，使用前加入终浓度为0.2％(2ml)的巯基乙醇。(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，反应体系如下表2：表2(7)进行连接反应，反应体系如下表3：表3(8)每个样品各取1μl反应产物，加入一个新的离心管中，总体积12μl。每12个样品一组。(9)用3％回收胶电泳1h，eb染色后切取300-700bp大小片段。用qiaquickkit进行胶纯化回收，回收产物溶于30μleb溶液。(10)进行pcr反应，pcr反应体系如下表4：表4pcr反应程序如下：94℃预变性5分钟；94℃变性30秒，60℃退火30秒，72℃延伸40秒，运行10个循环；最后72℃延伸3分钟。(11)磁珠纯化，完成建库。具体纯化方法为：首先，pcr后加入1.2倍体积磁珠，静置10min。然后，置于磁力架上吸附，去除上清。然后，加入500μl70％乙醇洗涤两次。干热仪上蒸干后，加入15μleb溶液溶解5min。最后，磁力架上吸附，转移上清于1.5ml离心管中。(12)检测达到上机标准后上机测序。结果：对444个样本(其中441个rils、2个亲本、1个f1个体)进行酶切建库测序，得到75.99gb原始数据，平均每个个体171.54mb。将测序序列比对到参考基因组(即yugu基因组，可以从下列网址获取yugu基因组的序列：https://www.ncbi.nlm.nih.gov/genome/？term＝setaria+italica+(foxtail+millet))上，平均比对率86.326％，平均覆盖率8.226％，平均测序深度3.655x。通过测序和信息分析，获得亲本间有多态性的33771个snp标记。根据窗口滑动的方法，选取若干个snp为一个窗口，每次滑动一个snp确定每个窗口的基因型以及每个个体的交换位点。最后生成bin基因型。根据bin基因型数据，用mstmap软件构建遗传图谱，2022个bin定位到9条染色体上，再将遗传图谱数据导入到mapchart软件，整合一张基因组遗传连锁图谱。利用所构建的谷子高密度遗传图谱，采用复合区间作图分析(cim)，对谷子穗长性状表型进行qtl分析。qtl检测采用5％的总体显著水平，根据500次排列检验结果，确定穗长性状表型数据qtl分析的临界lod值，分析得到与穗长相关的预测bin区间和snp位点信息。结果显示，与谷子穗长性状相关的分子标记pl-02-016，位于遗传连锁图谱第2号染色体237.81厘摩(cm)，第201个bin至第202个bin(chr2_bin201-chr2_bin202)，利用复合区间作图法，采用5％的总体显著水平，qtl检测的lod值为3.7556，表型变异解释率为3.35％，加性效应值(additiveeffect，a)为-0.8594。实施例2：snp标记验证根据预测的bin区间和snp位点，比对参考基因组，得到相关区域基因序列，选取snp位点前后300bp左右，设计开发snp标记引物，以父本和母本材料dna为模板进行pcr扩增。选择引物扩增正常、pcr产物符合预测大小的扩增产物，回收产物并进行测序，选择父本和母本材料扩增基因序列有snp位点差异的标记引物。预测的bin标记中包含多个snp位点，根据pcr结果对这些位点进行筛选。首先，根据pcr产物回收后的测序结果，选择父本和母本材料扩增基因序列有snp位点差异的标记。具体步骤如下：(1)按照实施例1中步骤(1)至(4)，用ctab法分别提取父母本基因组dna。(2)用0.8％的琼脂糖凝胶检测基因组dna，谷子父母本基因组dna电泳胶图如图1所示。(3)将得到的父母本基因组dna存于-20℃备用。(4)分别以提取的父本和母本的基因组dna为模板，利用7_2f：5’-gtatgctccacgcccttta-3’(seqidno：1)和7_2r：5’-ttgcgattaccacttgatt-3’(seqidno：2)扩增引物进行pcr扩增。pcr反应体系如下表5：表5试剂用量无菌水20.2μl10×缓冲液(含mg2+)2.5μldntps(25mm)0.15μltaq酶(5u/μl)0.15μl正向引物(10μmol/l)0.5μl反向引物(10μmol/l)0.5μl模板1.0μl总体积25μlpcr反应程序如下：94℃预变性5分钟；94℃变性30秒，60℃退火30秒，72℃延伸40秒，运行35个循环；最后72℃延伸3分钟。pcr扩增产物经纯化后于4℃下保存。各pcr扩增产物取部分进行1％琼脂糖凝胶电泳检测，结果如图2所示。通过对pcr扩增产物测序可知，引物7_2f和7_2r扩增得到的母本扩增产物的测序序列如序列表中seqidno：3所示，pl-02-016位点是seqidno：3所示序列自5’端起第148位碱基。母本扩增产物和父本扩增产物的测序序列比对结果如图3所示，其中标灰的碱基表示pl-02-016位点的碱基。然后，根据样品穗长数据筛选：父本穗长为26cm，母本穗长为19cm，441份样品中，最小值17cm，最大值47cm，穗长长度最短的50个样品中，该snp与母本相同的大于70％，穗长长度最长的50个样品中，与父本snp相同的样品数多于80％。以上内容是结合具体的实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属
技术领域：
的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。sequencelisting<110>张家口市农业科学院<120>与谷子穗长性状相关的分子标记及其检测引物和应用<130>16i23502<160>3<170>patentinversion3.3<210>1<211>19<212>dna<213>用于检测pl-02-016的上游引物<400>1gtatgctccacgcccttta19<210>2<211>19<212>dna<213>用于检测pl-02-016的下游引物<400>2ttgcgattaccacttgatt19<210>3<211>281<212>dna<213>pl-02-016位点所在的序列<220><221>snp位点<222>(148)..(148)<223>n是c或t<400>3ccctttagcagagcgattgcacgaaacaaaagcctcgtatccttcagccattatgacatc60tgcagacaggtcctgcctcgacaatttggtctcctacatacaacgcagaaggcattgata120cagtttggcaaatgaatcacaaaccaanccccatttcacccagttttactatttgcattt180gactgaacaaacgatccaataacatcccaccacagaaccgacaatttgatcagcattaaa240cgagatacactaatttttttctgaaatcaagtggtaatcgc281当前第1页12