本发明涉及一种基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法及引物，属于生物检测技术领域。
背景技术：
双孢蘑菇（agaricusbisporus），中文别名为蘑菇、洋蘑菇、白蘑菇。其含有丰富的蛋白质、多糖、维生素、核苷酸和不饱和脂肪酸，不仅营养丰富，肉质肥美，味道鲜美，而且热能低，具有很高的医疗保健作用，日益受到各国人民的喜爱。据中国食用菌学会统计，2015年全国双孢蘑菇总产量达337多万吨，其中双孢蘑菇工厂化日产达到417吨，说明双孢蘑菇在国内越来越受消费者们的青睐。
随着双孢蘑菇的种植业发展迅速，不仅和他们的栽培工艺先进有关，还和他们有适合本国工厂化生产的菌种有关。as2796是我国首次通过杂交育种的方法培养的双孢蘑菇新品种，并成为我国唯一的当家品种。但是，国内双孢蘑菇的工厂化品种大部分是进口国外的。没有真正的完全适合国内生产工艺的工厂化品种。这主要是因为双孢蘑菇育种者们对分离双孢蘑菇同核不育单孢菌株没有完全破解，制约了双孢蘑菇杂交育种的进程。通常育种者都是将收集到的担孢子进行出菇，如果不能出菇者将其视为单孢不育菌株。这个过程繁琐耗时耗力，并且不能判定得到的单孢菌株的交配型。
技术实现要素：
本发明目的在于提供一种基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法及引物。
为实现上述目的，本发明提供的技术方案是：一种基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的引物，所述引物包括：
abhd2f1：accttgtcttggtgagtggg；
abhd2r1：tctctacctcctctacagtt。
为实现上述目的，本发明提供的技术方案是：一种基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法，提取双孢蘑菇的基因组dna，然后以引物abhd2f1和引物abhd2r1进行pcr扩增，所得pcr扩增产物用内切酶进行酶切，接下来检测酶切产物；
aba1f1：accttgtcttggtgagtggg；
aba1r1：tctctacctcctctacagtt。
优选的技术方案为：所述pcr扩增反应的体系为：50ng/μl双孢蘑菇的基因组dna0.5μl，10μmol/l的引物2μl，0.1mmol/ldntp2μl，高保真酶pfu0.5μl，5μl的pcr缓冲液，用ddh2o补齐至50μl。
优选的技术方案为：所述pcr扩增反应的条件为：94℃预变性5min，94℃变性30s，55℃退火30s，72℃延伸1min，30个循环。
优选的技术方案为：酶切反应条件为：（1）a+因子的内切酶为0.5μl，a-因子的内切酶为0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11μl，总体积是20μl，37℃温育3小时；（2）a+因子的内切酶0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11.5μl，总体积是20μl，37℃温育3小时；（3）a-因子的内切酶为0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11.5μl，总体积是20μl，37℃温育3小时。
优选的技术方案为：如果pcr扩增产物只能被a+因子的内切酶酶切，则此菌株是同核不育单孢菌株的a+因子的交配型；如果pcr扩增产物只能被a-因子的内切酶酶切，则此菌株是同核不育单孢菌株的a-因子的交配型，如果pcr扩增产物既能被a+因子的内切酶酶切，也能被a-因子的内切酶酶切，则此菌株是双核体菌株。
优选的技术方案为：a+因子的内切酶选自内切酶xmni、bpmi、eco57mi；a-因子的内切酶选自内切酶clai、bspd、afⅲ、ndei、psti、bglⅲ、bstyi、bsaai、pmli。
优选的技术方案为：a+因子的内切酶为内切酶xmni，a-因子的内切酶为内切酶ndei。
优选的技术方案为：提取双孢蘑菇的基因组dna的方法包括下列步骤：用孢子收集器将成熟双孢蘑菇子实体上的孢子收集后，用ddh2o稀释成105cfu/ml浓度的孢子悬浮液，再取100μl孢子悬浮液涂布在pda平板上，再将在pda平板上长出来的单菌落转移到新的pda平板上，培养结束后直接挑取菌丝放在dna裂解液中，在80℃进行温浴5min后即得。
优选的技术方案为：提取双孢蘑菇的基因组dna的方法包括下列步骤：将双孢蘑菇菌株的菌丝在液体pda中25℃培养5-7天，收集菌丝，然后用质量浓度为1.5%的溶壁酶进行酶切，酶解温度为30℃，酶解时间是3h；将获得的原生质体用0.6m的甘露醇进行洗涤收集，将获得的原生质体稀释成105cfu/ml浓度的孢子悬浮液，再取100μl孢子悬浮液涂布在pda平板上，25℃培养；将原生质体萌发的单菌落菌丝放在dna裂解液中，在80℃进行温浴5min后即得。
由于上述技术方案运用，本发明与现有技术相比具有的优点是：
本发明利用pcr扩增获得的特异性片段，用特异性内切酶进行特异性酶切，通过检测酶切产物判断双孢蘑菇的类型，具有耗时短，操作简便，精准度高，无需出菇，成本低的优点。
附图说明
图1内切酶xmni和ndei区分孢子单核体电泳图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效。
实施例一：基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法及引物
本发明基于分子标记技术，建立一个简单快速鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法。
本发明是基于双孢蘑菇的交配型基因a因子序列，对双孢蘑菇的交配型基因的特征片段进行pcr扩增，扩增引物是seq.no.1和seq.no.2，扩增获得的特征片段在不同交配型的单核体中的长度都是462bp，我们将其中一个基因序列默认为交配型a+的序列，把另一个序列默认为交配型为a-的序列。
两个单核体的a+和a-因子进行序列比对，发现它们只有在snp位点上有差异。根据这些位点不同，a因子的两个序列的酶切位点不同，根据不同的内切酶位点进行酶切，根据酶切的图谱差异，筛选单核体菌株及其鉴别单核体的交配型。
这些不同的位点从a+因子到a-因子位置分别是：t99-a99，c102-g102，c155-t155，c195-g195，c206-g206，t247,-c247，a342-g342，c355-t355。
a+因子序列特有的酶切位点有：xmni195bp,bpmi/eco57mi206bp;a-因子特有的酶切位点有：clai/bspdi102bp,afⅲ155bp,ndei247bp,psti/bglⅲ/bstyi342bp,bsaai/pmli355bp。
其中内切酶xmni对a+因子特异性片段abd2酶切的片段大小分别为195bp和267bp，bpmi/eco57mi对a+因子特异性片段abd2酶切的片段大小分别为206bp和256bp；用clai/bspd酶切a-因子特异性片段abd2的大小为102bp和360bp,afⅲ酶切a-因子特异性片段abd2的大小为155bp和307bp,ndei酶切a-因子特异性片段abd2的大小为247bp和215bp，psti/bglⅲ/bstyi酶切a-因子特异性片段abd2的大小为342bp和120bp，bsaai/pmli酶切a-因子特异性片段abd2的大小为355bp和107bp。
利用获得的a+和a-因子序列上的酶切位点的不同，分别进行内切酶的组合，然后对pcr扩增获得的pcr产物进行酶切，再将获得的酶切产物5μl，进行1.5%的琼脂糖凝胶电泳，观察酶切的结果。
所述的pcr扩增条件是：94℃预变性5min，94℃变性30s，55℃退火30s，72℃延伸1min，30个循环，4℃保存。反应体系是：50ng/μl基因组dna0.5μl,10μmol/l的引物2μl，0.1mmol/ldntp2μl,保真性酶0.5μl，5μl的pcr缓冲液，用ddh2o补齐至50μl。
酶切反应条件为：a+因子的内切酶为0.5μl，a-因子的内切酶为0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11μl，总体积是20μl，37℃温育3小时。
一种基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法，提取双孢蘑菇的基因组dna，然后以引物abhd2f1和引物abhd2r1进行pcr扩增，所得pcr扩增产物用内切酶进行酶切，接下来检测酶切产物；
aba1f1：accttgtcttggtgagtggg；
aba1r1：tctctacctcctctacagtt。
优选的实施方式为：所述pcr扩增反应的体系为：50ng/μl双孢蘑菇的基因组dna0.5μl，10μmol/l的引物2μl，0.1mmol/ldntp2μl，高保真酶pfu0.5μl，5μl的pcr缓冲液，用ddh2o补齐至50μl。
优选的实施方式为：所述pcr扩增反应的条件为：94℃预变性5min，94℃变性30s，55℃退火30s，72℃延伸1min，30个循环。
优选的实施方式为：如果pcr扩增产物只能被a+因子的内切酶酶切，则此菌株是同核不育单孢菌株的a+因子的交配型；如果pcr扩增产物只能被a-因子的内切酶酶切，则此菌株是同核不育单孢菌株的a-因子的交配型，如果pcr扩增产物既能被a+因子的内切酶酶切，也能被a-因子的内切酶酶切，则此菌株是双核体菌株。
优选的实施方式为：a+因子的内切酶选自内切酶xmni、bpmi、eco57mi；a-因子的内切酶选自内切酶clai、bspd、afⅲ、ndei、psti、bglⅲ、bstyi、bsaai、pmli。
优选的实施方式为：a+因子的内切酶为内切酶xmni，a-因子的内切酶为内切酶ndei。
优选的实施方式为：提取双孢蘑菇的基因组dna的方法包括下列步骤：用孢子收集器将成熟双孢蘑菇子实体上的孢子收集后，用ddh2o稀释成105cfu/ml浓度的孢子悬浮液，再取100μl孢子悬浮液涂布在pda平板上，再将在pda平板上长出来的单菌落转移到新的pda平板上，培养结束后直接挑取菌丝放在dna裂解液中，在80℃进行温浴5min后即得。
本实施方式基于pcr扩增获得的特异性片段的snp位点的不同，交配型基因a+和a-因子的序列上的酶切位点不同。从筛选到的a+因子的酶切位点中选了xmni，其可以把a+因子序列酶切成195bp和267bp大小；我们从a-因子的酶切位点中选了内切酶ndei进行酶切，分别可以获得247bp和215bp大小序列。
将获得的酶切产物5μl，进行1.5%的琼脂糖凝胶电泳，观察酶切的结果，如图1所示，图1中，数字1-9代表挑选的9个孢子萌发的菌株，m代表dnamarker，从上到下条带大小分别是2000bp，1000bp,750bp,500bp,200bp和100bp，x代表内切酶xmni，n代表内切酶ndei。
酶切结果分为三种：只能被内切酶xmni切成195bp和267bp大小的序列，是交配型为a+；只能被内切酶ndei切开，获得247bp和215bp的序列，是交配型为a-;同时能被xmni和ndei切开的菌株为双核体菌株。
所述的pcr扩增条件是：94℃预变性5min，94℃变性30s，55℃退火30s，72℃延伸1min，30个循环，4℃保存。反应体系是：50ng/μl基因组dna0.5μl,10μmol/l的引物2μl，0.1mmol/ldntp2μl,保真性酶0.5μl，5μl的pcr缓冲液，用ddh2o补齐至50μl。
所述的酶切条件是：内切酶xmni0.5μl，内切酶psti、bglⅲ或者ndei0.5μl，hbuffer2μl，pcr产物6μl，ddh2o11μl，总体积是20μl,37℃温育3小时。
实施例二：基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法及引物
高保真酶pfu，dntps，酶buffer购自碧云天生物科技有限公司，内切酶xmni和xhoⅰ购自大连宝生，特异性引物aba1f1：accttgtcttggtgagtggg；aba1r1：tctctacctcctctacagtt。购自上海生工生物科技有限公司。
双孢蘑菇孢子收集来自上海联中食用菌合作社种植的菌株a15子实体。用包子收集器将成熟子实体上的孢子收集后，用ddh2o稀释成105浓度的孢子悬浮液，再取100μl涂布在pda平板上。再将在平板上长出来的单菌落转移到pda平板。
将分离得到菌株进行菌丝进行菌丝pcr，即直接挑取微量菌丝放在dna裂解液中，在80℃进行温浴5min，即可作为模板进行pcr。dna裂解液购自上海生工生物科技有限公司。
pcr扩增条件是：94℃预变性5min，94℃变性30s，55℃退火30s，72℃延伸1min，30个循环，4℃保存。反应体系是：50ng/μl基因组dna0.5μl,10μmol/l的引物2μl，0.1mmol/ldntp2μl,保真性酶0.5μl，5μl的pcr缓冲液，用ddh2o补齐至50μl。
将获得的pcr产物直接进行两种类型的单酶切，即分别xmni和bsaai进行酶切。
酶切条件是：内切酶ecorv为0.5μl，mbuffer的2μl,pcr产物6μl，ddh2o为11.5μl，总体积是20μl,37℃温育3小时。酶切的产物用1.5%的琼脂糖凝胶进行电泳，genefinder荧光染料进行染色并拍照。
内切酶xmni为0.5μl，mbuffer的2μl,pcr产物6μl，ddh2o为11.5μl，总体积是20μl,37℃温育3小时。酶切的产物用1.5%的琼脂糖凝胶进行电泳。内切酶bsaai为0.5μl，mbuffer的2μl,pcr产物6μl，ddh2o为11.5μl，总体积是20μl,37℃温育3小时。酶切的产物用1.5%的琼脂糖凝胶进行电泳。
经检测，bsaai只能被内切酶bsaai酶切，得到355bp和107bp片段。证明此菌株是同核不育单孢菌株的a-因子的交配型。经现有技术育种试验，结果相同。
实施例三：一种基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法及引物
提取双孢蘑菇的基因组dna的方法包括下列步骤：将双孢蘑菇菌株的菌丝在液体pda中25℃培养5-7天，收集菌丝，然后用质量浓度为1.5%的溶壁酶进行酶切，酶解温度为30℃，酶解时间是3h；将获得的原生质体用0.6m的甘露醇进行洗涤收集，将获得的原生质体稀释成105cfu/ml浓度的孢子悬浮液，再取100μl孢子悬浮液涂布在pda平板上，25℃培养；将原生质体萌发的单菌落菌丝放在dna裂解液中，在80℃进行温浴5min后即得。
其它部分同实施例二。
酶切条件为：酶切反应条件为：（1）a+因子的内切酶为0.5μl，a-因子的内切酶为0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11μl，总体积是20μl，37℃温育3小时；（2）a+因子的内切酶0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11.5μl，总体积是20μl，37℃温育3小时；（3）a-因子的内切酶为0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11.5μl，总体积是20μl，37℃温育3小时。
a+因子的内切酶选自内切酶eco57mi；a-因子的内切酶bglⅲ。
pcr扩增产物只能被a+因子的内切酶eco57mi酶切，采用现有技术对酶切产物进行测序，酶切产物的片段的长度分别为206bp和256bp，证明此菌株是同核不育单孢菌株的a+因子的交配型。经现有技术育种试验，结果相同。
实施例四：基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法及引物
其它部分同实施例三。
酶切条件为：酶切反应条件为：第一次实验：a+因子的内切酶0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11.5μl，总体积是20μl，37℃温育3小时；第二次实验：a-因子的内切酶为0.5μl，mbuffer2μl，pcr扩增产物6μl，ddh2o为11.5μl，总体积是20μl，37℃温育3小时。
a+因子的内切酶选自内切酶xmni；a-因子的内切酶bspd。
pcr扩增产物既能被a+因子的内切酶xmni酶切，也能被a-因子的内切酶bspd酶切，经现有技术对两种酶切产物进行测序，第一次实验的酶切产物的片段的长度分别为195bp和267bp，第二次实验的酶切产物的片段的长度分别为102bp和360bp，此菌株是双核体菌株。经现有技术育种试验，结果相同。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
序列表
<110>上海市农业科学院
<120>基于snp位点差异的鉴定双孢蘑菇同核不育单孢菌株及其交配型的方法及引物
<160>4
<210>1
<211>22
<212>dna
<213>人工序列
<400>1
accttgtcttggtgagtggg
<210>2
<211>22
<212>dna
<213>人工序列
<400>2
tctctacctcctctacagtt
<210>3
<211>730
<212>dna
<213>人工序列
<400>3
accttgtcttggtgagtgggatatgcattttgggcgaagtattttttcagaaaaggtgtc60
cattccttgaaaaacatcaataaataatcgctcaaggattccatcacttacaggcttgaa120
aactggccgctcatctgtagaactagcttgttttaagcgttcttgaactcgtatggaagc180
agcctggagaatggcttcctgccaactggtgatgtttcgagaataacgtgcatcgaggac240
acggtttatatgtttggcgagattgcgcagctcagtattggaggaatatttactgagccg300
tcggaagagtttttcatagtgacgctcgtaagtcgactgcaaatctttcacacgcgatgc360
aaagatgttcgcggcttcctgagactcttcgaaagcgcatcctagttgcacgagttgggc420
cgccagggatggtattgggggaaactgtagaggaggtagaga462
<210>4
<211>737
<212>dna
<213>人工序列
<400>4
accttgtcttggtgagtggggtatgcattttgggcaaagtattttttcaaaaaaggtgtc60
cattcctaaagaaacatcagtaaagaatcgctcaaggaatcgataacatacaggcttgaa120
agctggccgctcatctgtagaactagcttgtttagagcgttcttgaactcgcagggaagc180
agcctggagaatggtttcctgccaagtggtgatgtttcgagagtaacgtgcgtcgaggac240
atggttcatatgtttggcgagattgcgcaactcggtactggaggaatatttactgagccg300
tcggaaaagtttctcatagtgacgctcataagtcgactgcagatctttcacacgtgatgc360
aaatacgttcgcagcttcctgagactcttcgaaagtgcatcctagttgcacgagttgtgc420
caccggggatggtattgggggaaactgtagaggaggtagaga462