本发明属于作物遗传育种领域，涉及小麦植株千粒重判断标记及其应用。
背景技术：
小麦是全球约35-40％人口的主食，中国是全球最大的小麦生产国和消费国，常年产量约占全球总产17％。提高小麦产量对我国及全球粮食安全具有非常重要的战略意义。单位面积穗数、穗粒数和千粒重是构成小麦产量的三要素，其中，千粒重具有相对较高的遗传力，可以分解为粒长、粒宽、粒厚等基本要素。因此，鉴定千粒重相关位点，开发与其紧密连锁的分子标记将有助于小麦高产聚合育种。
千粒重属于数量性状，受主效-微效多基因共同控制，其相关qtl分布于小麦21条染色体，其中，主效qtl位于1a、3a、4b、5a、6a、7a、7b、7d，部分候选基因已被同源克隆，如位于6as染色体的tagw2；位于2a染色体的tacwi等。但是，千粒重的遗传机理较复杂，不同材料，遗传背景不同，控制千粒重的主效基因也存在差异。鉴定并开发出更多千粒重相关基因及其分子标记，不仅有利于加快多基因聚合育种的进程，同时有助于进一步阐明小麦产量形成的分子机制。
技术实现要素：
本发明提供一种小麦植株千粒重判断标记，该标记为jhmfei，序列如seqidno.1所示。
本发明还提供一种小麦植株千粒重判断标记引物，该引物分别为
jhmfei-f：ggaggaaacaaagtggtc，
jhmfei-r：gagaagaggaggtcgtag。
本发明还提供标记在小麦千粒重判断中的应用，该基因具有jhmfei-a类型和jhmfei-b类型，携带有jhmfei-a类型的小麦植株比携带jhmfei-b类型的小麦植株千粒重更高。
本发明还提供标记引物在小麦千粒重判断中的应用，具体方法为以所述引物对小麦基因组进行扩增，当扩增产物片段大小为489bp，则该小麦属于携带有jhmfei-a类型；当扩增片段大小为399bp和90bp，则该小麦属于携带有jhmfei-b类型；携带有jhmfei-a类型的小麦植株比携带jhmfei-b类型的小麦植株千粒重更高。
上述标记引物在小麦千粒重判断中的应用，其中，(1)所述pcr扩增体系为0.8ul2.5mmdntp,0.2ul10um标记引物jhmfei-f，0.2ul10um标记引物jhmfei-r1ul10×buffer,0.5utaqdnapolymerase,100ngdna模板，用双蒸馏水补齐到10ul；本发明中使用的是easytap品牌。
(2)扩增过程为94℃变性5min；94℃变性30s；60℃退火30s；72℃延伸30s，38个循环。
本发明还提供上述判断标记的获得方法，该方法包括如下步骤：
(1)提取小麦基因组dna；
(2)随机对小麦千粒重、粒长和粒宽性状进行测定；
(3)酶切方案设计和简化基因组测序；
(4)遗传图谱构建：利用slaf-seq方法，共得到在两个亲本之间存在多态性的slaf标签，构建亲本遗传图谱；采用highmap软件分析slaf标签在染色体上的线性排列，并估算相邻标签之间的遗传距离，最终得到遗传图谱；
(5)1b染色体千粒重主效qtl检测：利用软件在小麦1b染色体检测到一个粒重主效qtl，确定物理区间；本发明中采用软件qtlicimapping4.1检测小麦qtl。
(6)caps标记获得：以两个亲本为对照，根据163个家系的千粒重表型值，选择30个高千粒重家系的籽粒和30个低千粒重家系的籽粒，分别提取这60个家系的dna，然后等量混合构建高粒重池和低粒重池，用于660ksnp芯片检测；根据上述(5)中鉴定的粒重主效qtl物理位置，筛选该区间内在亲本以及两池间都存在一致突变的snp，利用软件查找酶切位点以及能够识别该位点的专一性限制性内切酶，采用软件在酶切位点上下游各1000bp设计引物，获得caps标记jhmfei；本发明所述的660ksnp芯片是委托北京博奥生物有限公司进行。本发明是利用primerpremier5软件设计引物。
(7)标记与植株千粒重相关性分析：利用连锁群亲本品种和自然群体进行主效位点与千粒重的相关性分析，验证了jhmfei对粒重具有显著的效应，为调控粒重的主效qtl，可以作为小麦植株千粒重标记。
上述判断标记获得方法，其中，连锁群体为京411/红芒春21，自然群体由369份小麦品种所组成。因此，步骤(4)共得到6718个slaf标签，获得的遗传图谱总图距为1310.38cm。步骤(5)所述的qtl位于marker24473338-marker236689383之间，物理区间为647895499bp-651868580bp。
上述判断标记获得方法，其中，步骤(3)所述酶切方案设计包括如下过程：利用酶切预测软件对参考基因组进行电子酶切预测，最终确定使用rsai酶，酶切片段长度选择在414-544bp的序列范围；其中，选择最适酶切原则：①位于重复序列的酶切片段比例尽可能低；②酶切片段在基因组上尽量均匀分布。
上述判断标记获得方法，其中，步骤(3)所述简化基因组测序包括如下步骤：将b过程获得的酶切大片段和datp在37℃条件下进行3’端加a处理；连接测序接头；pcr扩增；纯化；混样；切胶选取目的片段，测序；同时选用拟南芥(arabidopsisthalianaecotypecolumbia)作为对照进行测序。
上述判断标记获得方法，其中，步骤(3)还包括确定参考基因组方法，具体步骤如下：根据小麦的基因组大小以及gc含量等信息，最终选取小麦(triticumaestivumlinn.)a组基因组作为参考基因组进行酶切预测；其中，a组基因组为3.92gb，gc含量45.38％。
有益效果：
发明人利用千粒重差异较大的两个品种京411(千粒重42.54g)和红芒春21(千粒重18.73g)为亲本构建了ril群体(f10代)，两个亲本在千粒重上差异较大，有利于寻找控制千粒重的主效qtl，采用slaf-seq技术对两个亲本以及163个家系成员进行简化基因组测序，在全基因组水平上共开发出6718个多态性snp标记，且结合多年千粒重表型数据，在1b染色体上检测到控制千粒重的主效qtl，并获得与其紧密连锁的分子标记(marker24473338-marker24708085)，位于647-651mb之间。为进一步缩短该区间，继续以两个亲本为对照，根据163个家系成员的千粒重表型值，选择30个高千粒重家系和30个低千粒重的家系，分别提取各自基因组dna后等量混合，构建高低混合池，进行660k芯片扫描分型(委托北京博奥晶典有限公司完成)，获得在亲本以及两池间差异一致的snp标记用于加密目标区段。
本发明通过的简化基因组测序(specific-locusamplifiedfragmentsequencing，slaf-seq)，即利用生物信息学方法寻找合适的酶对基因组进行酶切，构建一定大小的插入片段文库，对其进行高通量测序，在全基因组范围内鉴定多态性标记(snps)。鉴定出来的大量多态性snps标记用于构建高密度遗传图谱以及定位控制千粒重的主效qtl。对于已经鉴定到的目标主效qtl，进一步利用小麦660k芯片分型上述两个亲本(京411和红芒春21)，以及所构建的高低混合池，寻找目标区间内(647-651mb)在两个亲本以及高低混合池间差异一致的snp标记，并将其转化成caps标记，用于针对性加密该目标区间，达到主效qtl精细定位目的。最终，确定了小麦植株千粒重判断标记，获得该标记的等位基因位点，不同基因类型与小麦千粒重存在显著相关性。根据标记基因获得的标记引物，使得规模化获得高千粒重小麦生产获得可能。本发明所述标记和标记引物尤其在小麦千粒重判断中为首次公开。
附图说明
图1caps标记jhmfei在京411和红芒春21中的酶切扩增产物电泳检测图谱。其中，a：京411和红芒春21扩增产物酶切前电泳检测图谱；b：京411和红芒春21扩增产物酶切后电泳检测图谱；1：京411，2：红芒春21
具体实施方式
下面结合实施例对本发明做进一步说明，下列实施例中未注明具体条件的实验方法，通常按照本领域的公知手段。
实施例1
(一)小麦基因组dna提取
用于小麦基因组dna提取的试验材料为亲本京411和红芒春21、京411/红芒春21群体(f10)163份家系材料、分别由来自京411/红芒春21群体的30份千粒重高和30份千粒重低的家系组成的两个高低混池、以及369份小麦品种。其中，京411/红芒春21群体163份家系材料用于构建高密度连锁图谱以及主效qtl定位；亲本京411和红芒春21，以及分别由来自京411/红芒春21群体的30份千粒重高和30份千粒重低的家系组成的高低混池用于660k芯片的扫描，主效qtl区间内差异snp的鉴定；369份小麦品种用于验证caps标记与千粒重的关系。369份小麦品种为市场所购，任意选择即可。
具体方法为：1、将小麦单籽粒磨碎，取0.1g加入0.7mlsds提取液(0.1mtris-hcl(ph＝8.5)，0.1mnacl，0.05medta(ph＝8.0)，2％sds)在60℃恒温下裂解45min，间歇震荡。
2、在4℃12000rpm条件下，离心10min。
3、取上清液加入等体积的酚：氯仿：异戊醇(25:24:1)上下颠倒，避免分层。
4、在4℃12000rpm条件下，离心10min。
5、取上清液加入等体积的酚：氯仿：异戊醇(25:24:1)上下颠倒旋转数次。
6、在4℃12000rpm条件下，离心10min。
7、取上清液加入等体积的异丙醇于-20℃冰箱静置30min。
8、在4℃12000rpm条件下，离心10min。
9、用70％的乙醇洗涤两次。
10、在4℃12000rpm条件下，离心10min。
11、沉淀自然风干，4℃存，备用。
(二)千粒重测定
利用万深种子分析仪对随机数取的小麦种子进行千粒重、粒长和粒宽性状测定，重复3次，计算平均值。
(三)简化基因组测序和酶切方案设计
1、参考基因组确定：选取中国春小麦基因组作为参考基因组进行酶切预测。
2、酶切方案确定：利用酶切预测软件对参考基因组进行酶切预测，选择最适酶切方案。
3、实验流程：根据选定的最适酶切方案，对检测合格的各样品基因组dna用限制性内切酶haeiii(newenglandbiolabs,neb)，进行酶切。对得到的酶切片段(slaf标签)用klenowfragment(3′→5′exo–)(neb)和datp在37℃下进行3′端加a处理、连接dual-index测序接头、pcr扩增(pcr扩增引物：faatgatacggcgaccaccgarcaagcagaagacggcatacg)、纯化(agencourtampurexpbeads(beckmancoulter,highwycombe,uk))、混样、切胶选取目的片段，文库质检合格后用illuminahiseqtm2500进行测序。为评估酶切实验的准确性，选用拟南芥(arabidopsisthalianaecotypecolumbia)作为对照(control)进行测序。
4、信息分析流程：利用dual-index对测序得到的原始数据进行识别，得到各个样品的reads。过滤测序reads的接头后，进行测序质量和数据量的评估。通过control数据评估酶切效率，以此判断实验过程的准确性和有效性。
(四)caps标记开发
1、引物设计：采用primerpremier5软件设计caps引物，由上海生物工程有限公司合成。
caps引物命名为jhmfei，分别为jhmfei-f：ggaggaaacaaagtggtc，
jhmfei-r：gagaagaggaggtcgtag。
2、pcr反应体系：0.8ul2.5mmdntp,0.4ul10um引物，1ul10*easytapbuffer,0.5ueasytap,100ngdna模板，用双蒸馏水补齐到10ul。
3、pcr反应程序：94℃变性5min；94℃变性30s；60℃退火30s(退火温度视引物而定)；72℃延伸30s，38个循环。
4、pcr产物检测：酶切产物利用2％琼脂糖凝胶电泳进行分析检测，gelstain荧光染料染色，bio-rad凝胶成像系统扫描拍照(图1)。
获得caps标记，具体序列如seqidno.1所示，在90位置具有c/t突变，90位置具有c突变为jhmfei-a类型；90位置具有t突变为jhmfei-b类型；
(五)caps标记与粒重的关系
同时利用京411/红芒春21群体和369份小麦品种组成的自然群体验证caps标记jhmfei与千粒重tgw(thousand-grainweight，tgw))的关系(表1和2)，结果表明，在6年环境下(2011、2012、2014、2015、2016、2017)，携带jhmfei-a类型的品种千粒重与携带jhmfei-b类型的品种的千粒重差异达到极显著水平(p<0.01)，且在京411/红芒春21群体中，该标记能够解释5.7-16.1％的表型变异(表2)。在自然群体中，携带jhmfei-a类型的品种的千粒重与携带jhmfei-b类型的品种的千粒重差异也达到极显著水平(p<0.01)，jhmfei-a的小麦植株比携带jhmfei-b类型的小麦植株千粒重更高。说明jhmfei对粒重影响较显著，为调控粒重的主效qtl。
表1caps标记jhmfei在ril群体和自然群体中与粒重的关系
表2caps标记jhmfei对千粒重的效应
可以知道，上述实施例仅为了说明发明原理而采用的示例性实施方式，然而本发明不仅限于此，本领域技术人员在不脱离本发明实质情况下，可以做出各种改进和变更，这些改进和变更也属于本发明的保护范围。
序列表
<110>安徽农业大学
<120>一种小麦植株千粒重判断标记及其应用
<160>1
<170>siposequencelisting1.0
<210>1
<211>489
<212>dna
<213>小麦(triticumaestivuml)
<220>
<221>misc_feature
<222>(90)
<223>niscort
<400>1
ggaggaaacaaagtggtctctctagtctactgccactccacttggtttggctgtctgttc60
gtttgctatacttaaagcgtgtaaaactgnaattggttattggttgttgtttagtcaggt120
accctcagtagtagtatatcccgtcatctgttttactgtaattcccacagcaataagttg180
tttttagaacgaaggctcaagtcgacccctgctttcaataaaaaagctatcaaccgacca240
ggattacagcatcaaccacaaaaacaaaaagacaaacctgctaagagatacatagtgctc300
tgcagctagcaagcttcatatggtctacaagcaaagatgaagaaaaactaaatcacagga360
acaacttgaaggcgatgaatccctccaccaagcgcactccaccacttagtgacaccaact420
agcaaattaagagtaacacattgtcatctccatgggctccggcaaattcagctacgacct480
cctcttctc489