欧洲红豆杉或曼地亚红豆杉的SNP分子标记及应用的制作方法

本发明涉及木材树种鉴定的分子生物学检测领域，具体地说，涉及欧洲或曼地亚红豆杉的snp分子标记，并利用该snp位点将上述两类树种从其他的红豆杉品种中鉴定出来的应用。

背景技术：

红豆杉，又叫紫杉，属于国家一级保护植物，是第四纪冰川以后遗留下来的世界珍稀濒危物种。红豆杉，既有高大的乔木，又有低矮的灌木，全世界约有几十种，曼地亚红豆杉就是其中之一。曼地亚红豆杉是红豆杉的其中一种，红豆杉(学名taxusmairei)是第四纪冰川时期孑遗植物，世界珍稀濒危物种，在全世界共有十一种，目前我国共有四个种和一个变种，即云南红豆杉、西藏红豆杉、东北红豆杉、中国红豆杉和南方红豆杉(变种)。曼地亚红豆杉是我国20世纪90年代中期从加拿大引种而来。现在四川、广西资源、山东等地有栽培。3年的曼地亚基本都是带花蕾的，雄性曼地亚的花蕾比较大比较明显，从8月份花蕾开始明显，一直到来年的3月左右开花，给雌株授粉。雌株花蕾比较小，和芽很相似，开花的时候是一个小水珠，没有花粉。

红豆杉有"植物活化石"和"国宝"之称。从中提取的紫杉醇是世界公认的抗癌药物，每公斤市场价在一千万美元以上，是世界上濒临灭绝的天然珍稀抗癌植物，属国家一级保护植物，至今已有250万年的历史。由于它的根、茎、叶、皮及种子含有紫杉醇，所以价格相当昂贵，不少国家将其列为"国宝"。美国、英国、俄罗斯、韩国和中国对紫杉进行多年研究，将紫杉醇从紫杉中提炼出来，用于临床治疗癌症。经验证明，紫杉醇具有独特的抗肿瘤机制和显著的抑制肿瘤的作用，被专家和学者确认为本世纪初国际上最有开发前途的抗癌药物。

形态学鉴定是欧洲红豆杉或曼地亚红豆杉鉴定中常采用的方法，但形状特征易受生境、气候、生理状况等的影响而常常导致主观辨别上的偏差，仅靠形态差异难以准确鉴定。

近年来，利用分子生物学技术进行物种鉴定已经在一些珍惜物种及检验检疫方面得到了应用，如何利用分子生物学技术建立高效、快速、灵敏且准确度高的鉴定方法将是欧洲红豆杉或曼地亚红豆杉领域的关键技术之一。

技术实现要素：

本发明提供了一种用于鉴定欧洲红豆杉或曼地亚红豆杉的snp分子标记，所述snp分子标记位于seqidno.1所示序列的第450位碱基，若第450碱基为c，则鉴定为欧洲红豆杉或曼地亚红豆杉。

所述snp分子标记在欧洲红豆杉或曼地亚红豆杉鉴定中的应用。

进一步的，从南方红豆杉、曼地亚红豆杉、欧洲红豆杉、西藏红豆杉、东北红豆杉和云南红豆杉中鉴定出欧洲红豆杉或曼地亚红豆杉。

本发明提供了一种用于鉴定欧洲红豆杉或曼地亚红豆杉的核苷酸序列，所述核苷酸序列包括seqidno.1，或者seqidno.1所对应的cdna或mrna，所述seqidno.1序列的第450碱基为c。

本发明提供了一种鉴定欧洲红豆杉或曼地亚红豆杉的试剂盒，包括用于特异性检测seqidno.1序列第450碱基为c存在与否的一种或多种引物，或引物和探针的组合。

进一步的，所述引物和探针选自：用于测序检测的特异性扩增引物和测序引物，引物序列分别为：

上游引物f1：5'-gcggtaggatcattgtcg-3'

下游引物r1：5'-tcctccgcttattgatatgc-3'

测序引物：5'-gcggtaggatcattgtcg-3'；

或者，用于荧光定量pcr检测的扩增特异性引物和探针，引物和探针序列分别为：

上游引物f2：5'-ccggtctgcaaagcagct-3'

下游引物r2：5'-ggaaccagtggacggtcg-3'

探针p2：5'-gattatttcgcgcctttggacg-3'；

或者，用于arms-pcr扩增的特异性扩增引物，引物序列为：

上游引物f3：5'-taatgccgtacagattatttcgtgc-3'

下游引物r3：5'-gagagccgagatatccgttgc-3'

进一步的，所述试剂盒用于从南方红豆杉、曼地亚红豆杉、欧洲红豆杉、西藏红豆杉、东北红豆杉和云南红豆杉中鉴定出欧洲红豆杉或曼地亚红豆杉。

所述的试剂盒还包括阳性质控品和阴性质控品。

本发明提供了一种制备鉴定欧洲红豆杉或曼地亚红豆杉试剂盒的方法，包括制备检测反应试剂的步骤，所述反应试剂包含用于特异性检测seqidno.1序列第450碱基为c存在与否的一种或多种引物，或引物和探针的组合。

本发明提供了一种用于鉴定欧洲红豆杉或曼地亚红豆杉的寡核苷酸，所述寡核苷酸用于特异性检测seqidno.1序列第450碱基为c存在与否，所述寡核苷酸选自：用于测序的特异性扩增引物序列和测序引物序列，分别为

上游引物f1：5'-gcggtaggatcattgtcg-3'

下游引物r1：5'-tcctccgcttattgatatgc-3'

测序引物：5'-gcggtaggatcattgtcg-3'；

或者，用于荧光定量pcr扩增特异性上游引物f2序列、下游引物r2序列和探针p2序列，分别为：

上游引物f2：5'-ccggtctgcaaagcagct-3'

下游引物r2：5'-ggaaccagtggacggtcg-3'

探针p2：5'-gattatttcgcgcctttggacg-3'；

或者，用于arms-pcr扩增的特异性上游引物f3序列和下游引物r3序列分别为：

上游引物f3：5'-taatgccgtacagattatttcgtgc-3'

下游引物r3：5'-gagagccgagatatccgttgc-3'

其中，所述探针5'端的荧光基团选自fam、thex或tet，3'端的淬灭基团选自bhq1或tamra。

本发明还提供了一种鉴定欧洲红豆杉或曼地亚红豆杉的方法，包括以下步骤：

提取待测植物的基因组dna；

以步骤(1)提取的基因组dna为模板，进行pcr扩增，检测seqidno.1序列第450碱基为c存在与否；

若经检测，snp位点为c时，则该待测植物鉴定为欧洲红豆杉或曼地亚红豆杉。

进一步的，pcr扩增所需的引物或探针选自上游引物f1、下游引物r1和测序引；或者，上游引物f2序列、下游引物r2序列和探针p2序列；或者，上游引物f3和下游引物r3。

本发明与现有技术相比具有以下优点和效果：本发明在对红豆杉大量基因信息进行分析比较的基础上，并根据序列位点多态性，设计欧洲红豆杉或曼地亚红豆杉种属特异引物和探针，建立了欧洲红豆杉或曼地亚红豆杉特异性的荧光定量pcr鉴定检测方法。本发明并不需要木材有完整的形态结构，只需从木料上取极少量木材样本，即可满足检测要求，操作简单。并且通过欧洲红豆杉或曼地亚红豆杉snp分子标记做出判断，鉴定过程客观、准确。克服了传统鉴定方法主要依赖于形态鉴定方法，并必须建立在木材有完整易辨识的形态结构的基础上的技术局限性。本发明采用arms技术利用特异引物对突变靶序列进行高精准pcr扩增放大，与此同时，利用探针对扩增产物进行检测，在实时荧光定量pcr平台上实现对欧洲红豆杉或曼地亚红豆杉检测的高特异性和高灵敏度。

附图说明

图1-1是利用本发明荧光定量pcr法鉴定阳性质控品(欧洲红豆杉叶片基因组dna)的结果

图1-2是利用本发明荧光定量pcr法鉴定阴性质控品(南方红豆杉叶片基因组dna)的结果。

图2是利用本发明利用荧光定量pcr法鉴定曼地亚红豆杉的结果。

图3是利用本发明利用荧光定量pcr法鉴定欧洲红豆杉的结果。

图4是利用本发明利用荧光定量pcr法鉴定南方红豆杉的结果。

图5是利用本发明利用荧光定量pcr法鉴定西藏红豆杉的结果。

图6是利用本发明利用荧光定量pcr法鉴定云南红豆杉的结果。

图7是利用本发明利用荧光定量pcr法鉴定东北红豆杉的结果。

图8是测序法鉴定曼地亚红豆杉的结果。

图9是测序法鉴定欧洲豆杉的结果。

图10是测序法鉴定南方红豆杉的结果。

图11是测序法鉴定西藏红豆杉的结果。

图12是测序法鉴定云南红豆杉的结果。

图13是测序法鉴定东北红豆杉的结果。

图14是arms-pcr法电泳结果：1泳道为质控品，2泳道为曼地亚豆杉；3泳道为欧洲红豆杉；4泳道为南方红豆杉；5泳道为西藏红豆杉；6泳道为云南红豆杉；7泳道为东北红豆杉；8泳道为100bpladdermaker。

图15六种不同品种红豆杉的序列对比图

具体实施方式

下面结合具体实施例，进一步阐述本发明。但并不因此将本发明限制在所述的具体实施例中。实施例中未说明的条件和方法，通常按照所属领域实验人员常规采用方法：譬如，sambrook等分子克隆实验手册(newyork:coldspringharborlaboratorypress，1989)中所述的操作技术规程，或按照制造厂商所建议的实验条件。

实施例1欧洲或曼地亚红豆杉snp分子标记的获得

1、根据ncbi数据库公开的its基因，设计相应的引物进行pcr扩增，采用上游引物f1：5'-gcggtaggatcattgtcg-3'(seqidno.2)，下游引物r1：5'-tcctccgcttattgatatgc-3'(seqidno.3)，和如下的pcr反应体系和pcr反应程序，以经鉴定的南方红豆杉(taxusmaireisyhu)、曼地亚红豆杉(taxusmadia)、欧洲红豆杉(taxusbaccata)、西藏红豆杉(taxuswallichianazucc.)、东北红豆杉(taxuscuspidatas.etz.)和云南红豆杉(taxusyunnanensischengetl.k.f)6个红豆杉属物种基因组dna为模板，进行pcr扩增，扩增产物经1.5％琼脂糖凝胶电泳检测并在紫外凝胶成像系统观察合格后进行测序。

其中，vazymeacetaqdna酶(货号p401-d1250u)。

按传统方法提取植物基因组dna，也可采用商品化dna提取试剂盒进行dna提取。例如，所述红豆杉属物种基因组dna的提取采用杭州百迈生物股份有限公司的植物基因组dna提取试剂盒(koing^tm通用型基因组提取试剂盒，货号f1612)，参照试剂盒说明进行操作，将提取的dna置于-20℃下保存备用2、序列比对及物种特异的分子标记的开发

将步骤1获得的6个红豆杉属物种的its基因序列，通过dnaman软件进行序列比对寻找合适的snp位点。序列比对结果如图15所示，经分析发现位于欧洲红豆杉和曼地亚红豆杉its基因内的seqidno.1所示序列第450位上的碱基为c，曼地亚红豆杉its基因内相应于seqidno.1所示序列上的碱基为c。而其它红豆杉属物种相应位置为t。即seqidno.1所示序列的第450位上的碱基为c是欧洲红豆杉或曼地亚红豆杉的snp位点(序列中方框所示位置)，seqidno.1核苷酸序列为：

实施例2欧洲或曼地亚红豆杉snp分子标记在荧光pcr扩增法鉴定欧洲或曼地亚红豆杉中的应用

本实施例是在获得欧洲或曼地亚红豆杉snp位点的基础上，针对包含所述snp位点的基因序列，利用primerpremier5.0分别设计出特异性上游引物f2、下游引物r2和探针p2，引物和探针由苏州金唯智生物科技有限公司合成。利用所述特异性引物和探针对南方红豆杉、曼地亚红豆杉、欧洲红豆杉、西藏红豆杉、东北红豆杉和云南红豆杉进行荧光pcr扩增检测。具体的，所述上游引物f2、下游引物r2和探针p2分别是：

上游引物f2：5'-ccggtctgcaaagcagct-3'(seqidno.4)

下游引物r2：5'-ggaaccagtggacggtcg-3'(seqidno.5)

探针序列p2：5'-gattatttcgcgcctttggacg-3'(seqidno.6)

探针5'端的荧光基团选自fam、thex、tet等，3'端的淬灭基团：bhq1、tamra等。

例探针为5'-fam-gattatttcgcgcctttggacg-bhq1-3'。

以6个红豆杉属物种基因组dna为模板，20μlpcr反应体系中包括：

pcr扩增程序：

包括以下步骤：1)提取待测植物的基因组dna；2)以待测植物的基因组dna为模板，利用引物和探针，pcr扩增红豆杉样品its基因；同时设置阳性和阴性对照组，所述阳性和阴性对照组分别以含有提取的欧洲红豆杉叶片基因组dna和南方红豆杉叶片基因组dna为模板进行扩增；3)采用荧光pcr法检测pcr扩增产物，扩增曲线有标准的s型曲线，且ct≤36，则待测植物为欧洲或曼地亚红豆杉物种。

本实施例2采用的南方红豆杉、曼地亚红豆杉、欧洲红豆杉、西藏红豆杉、东北红豆杉和云南红豆杉样品也是经现有方法鉴定不同与实施例1的另一组样品，按本发明所述方法对样品进行荧光pcr法检测，检测结果图1至7所示，其中图1-1为阳性质控品(欧洲红豆杉叶片)的扩增结果，图1-2为阴性质控品(南方红豆杉叶片)的扩增结果；其余样本中，仅欧洲红豆杉和曼地亚红豆杉出现s型曲线(如图2或图3所示)，且ct≤36,而其他红豆杉类没有出现s型曲线(如图4至7所示)，这表明本发明针对欧洲红豆杉和曼地亚红豆杉its基因snp位点具有极高的特异性，因此可以用于快速鉴定欧洲红豆杉和曼地亚红豆杉。

利用实施例1所述引物和测序方法，对本实施例2中的南方红豆杉、曼地亚红豆杉、欧洲红豆杉、西藏红豆杉、东北红豆杉和云南红豆杉样品进行鉴定，鉴定结果见图8至13。测序结果为只有欧洲红豆杉和曼地亚红豆杉在第450位上的碱基为c，而其它红豆杉属物种为t。

实施例3欧洲或曼地亚红豆杉snp分子标记在arms-pcr法鉴定西藏红豆杉中的应用

本实施例是在获得欧洲或曼地亚红豆杉snp位点的基础上，针对包含所述snp位点的基因序列，利用primerpremier5.0分别设计出特异性上游引物f3、下游引物r3，引物由苏州金唯智生物科技有限公司合成。利用所述特异性引物和探针对南方红豆杉、曼地亚红豆杉、欧洲红豆杉、西藏红豆杉、东北红豆杉和云南红豆杉样品进行荧光pcr扩增检测。具体的，所述上游引物f3序列、下游引物r3序列分别是：

上游引物f3：5'-taatgccgtacagattatttcgtgc-3'(seqidno.7)

下游引物r3：5'-gagagccgagatatccgttgc-3'(seqidno.8)

以6个红豆杉属物种基因组dna为模板，20μlpcr反应体系中包括：

arms-pcr反应程序：

包括以下步骤：1)提取待测植物的基因组dna；2)以待测植物的基因组dna为模板，利用引物和探针，arms-pcr扩增红豆杉its基因；3)采用电泳法检测pcr扩增产物，欧洲红豆杉和曼地亚红豆杉的电泳有明显的条带，而其他红豆杉则无扩增条带。

本实施例2采用的南方红豆杉、曼地亚红豆杉、欧洲红豆杉、西藏红豆杉、东北红豆杉和云南红豆杉也是经现有方法鉴定不同与实施例1的另一组样品，检测结果如图14所示，欧洲红豆杉和曼地亚红豆杉的扩增结果为有明显条带(第2和3泳道)；其他红豆杉的扩增结果为无扩增条带。这表明本发明的针对欧洲红豆杉和曼地亚红豆杉its基因snp位点具有极高的特异性，因此可以用于快速鉴定欧洲红豆杉和曼地亚红豆杉。

实施例2和3的鉴定结果汇总表见下表：

本申请实施例中所列反应体系，反应程序和检测步骤等只是一种具体的实施方式，本领域技术人员还可以根据实际需要对反应体系，反应程序等作出适应性的调整。

序列表

<110>杭州百迈生物股份有限公司

<120>欧洲红豆杉或曼地亚红豆杉的snp分子标记及应用

<130>201704

<150>2017108141463

<151>2017-09-11

<160>8

<170>siposequencelisting1.0

<210>1

<211>1099

<212>dna

<213>欧洲红豆杉(taxusbaccata)

<400>1

gaacgtgaggagagtgggcgtgttcggcacgtccggcgcgtctcctccccgtcccccgct60

cttgcgggagtttggactggagttccggatgcttttccctttgcaatgcctcgggggtcc120

cctccgctggctttctctcgggtgggcgttccggtcttcgcaagttgttctcctgcggtt180

cgggttattacaagggcgagtgccagtgtcgtctacagtgtggctatgctaggggctcgt240

tttcggctgcagtgggtgcaccagcgactctgttcgtctggcagtcctgctcctttgcgg300

tcgaagggtcccagcaatttccacctgctaggcggttccttgcaaacgcccgtggctgct360

ttgcagaccgttgaacatttgccggtctgcaaagcagctgcaaccccccgaagtgtcgcc420

gggcaggtaatgccgtacagattatttcgcgcctttggacgggtgcacctgcgtacgtcg480

ccgggcgggacgaccgtccactggttccgttctcatctgttcggtctccggcctttcgac540

actccggggcggggacaccctttctccgattttcccccacggggcatcgcctctggttgc600

aggggggcgtcgagtggaggtgcgtgtcataacactcaacgcatcggtgcggattgcacc660

aaggatctgataaatagctgtgcgggctgcgcgccatgcgcgcgtcggactcacggcgcc720

aattaaacacgactctcggcaacggatatctcggctctcgccacgatgaagaatgtagcg780

aaatgcgatacttagtgtgaattgcagaatcccgtgaatcatcgagtctttgaacgcaag840

ttgcgcccgaggcctcggccgagggcacgtctgcttgggcgtcgcactctaaaatcgact900

tcccccgttcctcggcgggaggatcggagatggctgtccgtgcccccagcggcgcggtcg960

gctcaatgggcacgagatttgtggccgatgtcaccatgagcggtggccaatgttgttcgg1020

cgttggattgtggcatggggtcaccgatttgggcggaactttactccttgttaatccgaa1080

tcttccccgggagggccgt1099

<210>2

<211>18

<212>dna

<213>人工序列(artificialsequence)

<400>2

gcggtaggatcattgtcg18

<210>3

<211>20

<212>dna

<213>人工序列(artificialsequence)

<400>3

tcctccgcttattgatatgc20

<210>4

<211>18

<212>dna

<213>人工序列(artificialsequence)

<400>4

ccggtctgcaaagcagct18

<210>5

<211>18

<212>dna

<213>人工序列(artificialsequence)

<400>5

ggaaccagtggacggtcg18

<210>6

<211>22

<212>dna

<213>人工序列(artificialsequence)

<400>6

gattatttcgcgcctttggacg22

<210>7

<211>25

<212>dna

<213>人工序列(artificialsequence)

<400>7

taatgccgtacagattatttcgtgc25

<210>8

<211>21

<212>dna

<213>人工序列(artificialsequence)

<400>8

gagagccgagatatccgttgc21