CN108950039B

CN108950039B - Dna条形码、引物、试剂盒、方法和应用

Info

Publication number: CN108950039B
Application number: CN201710351860.3A
Authority: CN
Inventors: 徐平; 唐蜀昆; 田飞; 施佳辉; 周文婧; 高林瑞; 高慧英; 职晓阳; 丁章贵; 贾鳗
Original assignee: Yunnan Dayi Microbial Technology Co ltd; Menghai Tea Industry Co ltd
Current assignee: MENGHAI TEA INDUSTRY Co.,Ltd.; Yunnan Dayi Microbial Technology Co., Ltd
Priority date: 2017-05-18
Filing date: 2017-05-18
Publication date: 2021-08-13
Anticipated expiration: 2037-05-18
Also published as: CN108950039A

Abstract

本发明属于物种和菌株鉴定领域，具体涉及一种DNA序列、一种用于鉴别食腺嘌呤节孢酵母菌株的DNA条形码、引物、试剂盒、方法和应用。该DNA条形码来源于食腺嘌呤节孢酵母TMCC 70007菌株的基因组，并且包含根据本发明所述的DNA序列，并且该DNA条形码的长度为500bp‑2300bp，优选为500bp‑1500bp，更优选为500bp‑1300bp。该DNA条形码能够准确鉴别普洱茶发酵菌株食腺嘌呤节孢酵母TMCC 70007，可以快速准确地将该菌株从易混淆的菌株或同物种内其他菌株中鉴定出来。

Description

DNA条形码、引物、试剂盒、方法和应用

技术领域

本发明属于物种和菌株鉴定领域，具体涉及DNA条形码、引物、试剂盒、方法和应用。

背景技术

普洱茶是一种具有云南地理标识的后发酵茶，其采用大叶种晒青毛茶作为原料，经过一系列工艺制作而成。传统的普洱茶制作过程为：采摘的新鲜茶叶经揉捻、晒青、除杂、潮水、渥堆、晾干、筛分、压制成型、包装出厂。在普洱茶的生产中，渥堆发酵过程是普洱茶品质形成的主要因素，在此过程中，茶叶中诸如茶多酚、咖啡因以及一些多糖物质等内含成分发生了巨大变化，成就了普洱茶的特殊风味、口感、品质以及各种保健功效。

传统的普洱生产中，湿热的环境激活茶叶本身所含的酶，把茶叶中所含的一部分内含成分转化成微生物能利用的物质，微生物大量滋生，产生丰富的胞内酶和胞外酶，催化茶叶中包含的成分发生一系列转变，再加上湿热等其他因素，逐渐形成不同品质的普洱茶。不同的产地，因为其微生物群落结构的不同，风味及品质也各具特色。

普洱茶除了其独特的风味和文化外，其诸如减肥、降血糖血脂、预防和改善心血管疾病、抗衰老、抗癌、消炎、助消化以及养胃等保健功效也备受人们关注。近年来，普洱茶越来越受欢迎，逐渐增加的市场需求拉动了普洱茶产业的发展，促进了云南地区经济增长。

然而，普洱茶也面临着一些困境，例如产品质量不稳定、存在一定安全隐患、生产周期长、劳动力投入过高、微生物数量超标和螨虫滋生等。随着人们生活水平的提高，消费者对食品的卫生、安全等问题日益重视。近几年食品质量安全事件频发，茶叶质量的安全问题也日益突出，除了农药的残留问题之外，渥堆发酵过程中微生物也有可能成为影响茶叶质量的重要因素。

目前多数厂家的普洱茶生产仍然是半自然人工渥堆的经验式发酵，虽然渥堆过程中以食腺嘌呤节孢酵母(Blastobotrys adeninivorans)等优势共性微生物为主的群落相对稳定，但产品的稳定性还存在较大的提升空间。生产过程中也不可避免的存在着一定的安全隐患，为了进一步获得消费者的青睐和市场的认可、突破外贸壁垒、提升普洱茶企业的市场竞争力，必须实现普洱茶生产工艺的人工可控化、清洁化以及高效化，并不断开发新产品、延伸产业链。要做到这些，就必须革新技术，发明一系列安全、清洁、高效、人工可控自动化的普洱茶新工艺，才能确保普洱茶产业的健康发展，给国家和人民带来长远利益。

随着普洱茶科研的进一步深入，越来越多的人们开始了人工接种发酵普洱茶的探索，目前已经有许多菌株应用于普洱茶人工可控发酵。近年由于市场需求加大，许多规模较小的厂家对普洱茶缺乏系统深入的研究，大肆滥用许多他人专利进行牟利，比如，按照一些接种发酵普洱茶的专利方法，采用一些菌种进行普洱茶发酵等一系列方法，进行普洱茶加工生产等经营活动，大大损害了一些具有相关专利的普洱茶生产企业的经济利益。

为确保普洱茶发酵微生物种质资源得到有效保护，防止滥用普洱茶发酵微生物的侵权行为发生，解决在普洱茶人工可控发酵过程中菌种侵权时难以举证的难题，开发对普洱茶发酵用菌株进行快速精准鉴别的方法势在必行，亟需建立其分子鉴定方法，开发通过DNA条形码技术准确鉴别与区分与其相似菌株的方法，以达到通过形态特征与分子数据分析相结合的方法来解决工业生产菌株鉴别和鉴定问题，为普洱茶有控工业发酵工艺开发和资源保护提供理论依据，促进普洱茶产业健康繁荣发展。

发明内容

为了克服形态学鉴定用于普洱茶发酵生产的食腺嘌呤节孢酵母菌株所存在的不足，本发明提供了一种用于鉴别食腺嘌呤节孢酵母菌株的DNA条形码、引物、试剂盒、方法和应用，从而对食腺嘌呤节孢酵母TMCC 70007菌株进行快速鉴别和区分，可为普洱茶新发酵工艺提供快速评估，防止发酵过程中其他杂菌的干扰，以及为普洱茶人工可控发酵工艺及菌种非法滥用提供举证方法和依据。

为了实现上述目的，本发明采用如下技术方案：

(1)一种分离的或合成的DNA序列，其核苷酸序列如SEQ ID No.2所示。

(2)一种用于鉴别食腺嘌呤节孢酵母菌株的DNA条形码，该DNA条形码来源于食腺嘌呤节孢酵母TMCC 70007菌株的基因组，并且包含根据(1)所述的DNA序列，并且该DNA条形码的长度为500bp-2300bp，优选为500bp-1500bp，更优选为500bp-1300bp。

(3)根据(2)所述的DNA条形码，其核苷酸序列如SEQ ID No.4所示，优选如SEQ IDNo.1所示。

(4)一种用于扩增(2)所述的DNA条形码的引物对。

(5)根据(4)所述的引物对，其正向引物的核苷酸序列与食腺嘌呤节孢酵母TMCC70007菌株基因组中的这样的序列相同：该序列为从所述TMCC 70007菌株基因组中如SEQID No.2所示的核苷酸序列的第1位的上游1000bp到如SEQ ID No.2所示的核苷酸序列的第1位的区域中的序列，并且该正向引物的长度为20-30bp；其反向引物与所述TMCC 70007菌株基因组中的这样的序列反向互补：该序列为从所述TMCC 70007菌株基因组中如SEQ IDNo.2所示的核苷酸序列的最后一位到如SEQ ID No.2所示的核苷酸序列的最后一位的下游1000bp的区域中的序列，并且该反向引物的长度为20-30bp。

(6)根据(4)所述的引物对，所述正向引物和反向引物的核苷酸序列分别如下所示：

正向引物：5’-CACTGCAACGGTCGCATAAG-3’；

反向引物：5’-GGTAACGGGCCACAGATTCA-3’。

(7)一种用于鉴别食腺嘌呤节孢酵母菌株的试剂盒，包含根据4(4)至(6)所述的引物对。

(8)一种用于鉴别食腺嘌呤节孢酵母菌株的方法，包括以下步骤：

a)提供待测菌株的基因组DNA；

b)以步骤a)所述的基因组DNA为模板，利用根据(4)所述的引物对进行PCR扩增，得到PCR产物；

c)电泳检测PCR产物，如果没有目标条带，则判定待测菌株不是食腺嘌呤节孢酵母TMCC 70007菌株，如果有目标条带，则进行步骤d)；

d)对所得PCR产物进行测序，得到待测核苷酸序列；将所述待测核苷酸序列与(2)所述的DNA条形码的核苷酸序列进行同源性比对，若同源性在97％以上，则判定待测菌株是食腺嘌呤节孢酵母TMCC70007菌株。

(9)根据(2)所述的DNA条形码在鉴别食腺嘌呤节孢酵母菌株中的应用。

(10)根据(4)所述的引物对在鉴别食腺嘌呤节孢酵母菌株中的应用。

本发明与现有技术相比具有以下优点和积极效果：

1、本发明采用蛋白质基因组学技术从食腺嘌呤节孢酵母TMCC70007菌株的基因组中发现了漏注释的蛋白编码基因AraC，并且发现该基因的蛋白编码框(SEQ ID No.2)为食腺嘌呤节孢酵母TMCC70007菌株的基因组所特有。本发明进一步通过仔细的研究和比较分析，将该特有的DNA序列开发为了DNA条形码，作为鉴别普洱茶工业发酵生产菌株食腺嘌呤节孢酵母TMCC 70007的有效工具。该条形码为食腺嘌呤节孢酵母TMCC 70007菌株基因组中包含上述特有的蛋白编码框的一段500-2300bp的序列，该条形码序列能实现食腺嘌呤节孢酵母种内菌株的快速准确鉴别与区分。

2、本发明进一步发现，相比于其他基因，AraC基因中的序列(例如SEQ ID No.1)具有通用、易扩增、易比对的特点，其在食腺嘌呤节孢酵母种内不同菌株之间差异明显。

3、本发明建立了普洱茶工业发酵生产菌株食腺嘌呤节孢酵母TMCC 70007的标准基因序列和样品鉴别方法，相比于传统的形态学鉴定方法，显著提高了鉴定效率。该方法对于样品的完整程度要求低，鉴别指标能够量化，这对于及时判定普洱茶及其种质资源提供了有效的依据。此外，进一步加入了形态学混淆种并运用基于聚类分析的系统发育树法设立鉴定规则进行了鉴定，其鉴定的可靠性和准确性也大大优于常规的分子鉴定方法，弥补了该类普洱发酵茶生产菌株基于DNA条形码技术的菌株鉴定的空白。

附图说明

图1示出了支持AraC编码基因编码LEGDYMEEWLR和TYAMSEIGAHFSSQR的肽谱匹配(Peptide Spectrum Matches，PSMs)图谱；其中图1A为LEGDYMEEWLR的质谱图；图1B为TYAMSEIGAHFSSQR的质谱图；

图2示出了肽段所在区域的蛋白质编码框的mRNA序列及其编码的蛋白质序列对应图，图中灰色部分标出了肽段LEGDYMEEWLR和TYAMSEIGAHFSSQR；

图3示出了AraC的部分蛋白编码区的扩增产物核苷酸序列，其中前后两端有下划线的区域即为引物所在区域，灰色背景序列为本发明用蛋白质基因组学技术发现的新肽段所对应的蛋白编码框，前后两个终止子分别为TAA和TAG；

图4示出了SEQ ID NO.3通过KEGG-BLASTP分析后得到的与之具有同源性的序列比对图，图4A中，Sbjct＝dti:Desti_1638 3-羟基-3-甲基戊二酰辅酶A合成酶；在图4B中，Sbjct＝ach:Achl_3576转录调控因子，AraC家族；

图5示出了条形码序列SEQ ID NO.1通过NCBI-BLASTN得到的与之具有同源性的序列比对图，图5A上方的值表示比对得分；5B中的Sbjct表示沙柳支孢瓶霉(Cladophialophora psammophila)CBS110553假定蛋白部分mRNA(登录号：XM_007742452.1)；

图6示出了食腺嘌呤节孢酵母TMCC 70007的条形码序列SEQ ID NO.1与食腺嘌呤节孢酵母LS3菌株的同源序列的比对图；

图7示出了对利用本发明的引物得到的待测菌株PCR扩增产物进行琼脂糖凝胶电泳的结果；

图8示出了SEQ ID NO.1序列与待测菌株PCR扩增序列的比对图。

具体实施方式

以下通过具体实施方式的描述并参照附图对本发明作进一步说明，但这并非是对本发明的限制，本领域技术人员根据本发明的基本思想，可以做出各种修改或改进，但是只要不脱离本发明的基本思想，均在本发明的范围之内。

本发明采用蛋白质基因组学技术从食腺嘌呤节孢酵母TMCC70007菌株的基因组中发现了漏注释的蛋白编码基因AraC，并且发现该基因的蛋白编码框(SEQ ID No.2)为食腺嘌呤节孢酵母TMCC70007菌株的基因组所特有。因此，该特有的DNA序列可以用于开发鉴别普洱茶工业发酵生产菌株食腺嘌呤节孢酵母TMCC 70007的DNA条形码。

因此，本发明的一个方面提供了一种分离的或合成的DNA序列，其核苷酸序列如SEQ ID No.2所示。

考虑到DNA条形码应当具有合适的长度和足够的菌株间特异性，本发明进一步通过仔细的研究和比较分析，基于上述特有的DNA序列获得了可准确、有效地鉴别食腺嘌呤节孢酵母TMCC 70007菌株的DNA条形码。该条形码为食腺嘌呤节孢酵母TMCC 70007菌株基因组中包含上述特有的蛋白编码框的一段500-2300bp的序列，该条形码序列能实现食腺嘌呤节孢酵母种内菌株的快速准确鉴别与区分。

具体而言，本发明另一方面提供了一种用于鉴别食腺嘌呤节孢酵母菌株的DNA条形码，该DNA条形码来源于食腺嘌呤节孢酵母TMCC70007菌株的基因组，并且包含本发明所述的核苷酸序列如SEQ ID No.2所示的DNA序列，并且该DNA条形码的长度为500bp-2300bp，优选为500bp-1500bp，更优选为500bp-1300bp。

由于SEQ ID No.2为食腺嘌呤节孢酵母TMCC 70007菌株基因组所特有，因此包含其的、长度如上所述的DNA条码的特异性更加得到保证。另一方面，通常，DNA条形码的长度大于500bp时，也能够满足易扩增、易比对的操作要求。当DNA条形码的长度过长时(例如大于2500bp)，对于扩增操作而言较不理想。

在本发明的一个优选的实施方案中，所述DNA条形码序列如SEQ ID No.4所示，更优选如SEQ ID No.1所示。

在本文中，术语“漏注释基因”是指物种完成基因组测序以后，采用基因预测软件(例如GeneMark、Augustus、Glimmer等)未能将该基因预测出来，这些基因一般表达量不高、在特定的条件下才表达，所以在研究中难以被发现。

术语“DNA条形码技术(DNA barcoding)”是指用基因组内一段标准的、短的DNA片段来鉴定物种的一项分子鉴定新技术，其可以快速、准确的进行物种鉴定。

术语“六框翻译”是蛋白质组学和基因组学中的已知术语，简述其原理为，DNA编码蛋白质时，采用三联体密码子编码蛋白质，给定一条DNA序列，有3种编码可能性，加上其互补链上的3种编码可能，共有6种编码可能性(+1，+2，+3，-3，-2，-1)。

本发明的另一个实施方案提供了一种用于扩增本发明所述的DNA条形码的引物对。

优选地，其正向引物的核苷酸序列与食腺嘌呤节孢酵母TMCC70007菌株基因组中的这样的序列相同：该序列为从所述TMCC70007菌株基因组中如SEQ ID No.2所示的核苷酸序列的第1位的上游1000bp到如SEQ ID No.2所示的核苷酸序列的第1位的区域中的序列，并且该正向引物的长度为20-30bp；其反向引物与所述TMCC70007菌株基因组中的这样的序列反向互补：该序列为从所述TMCC70007菌株基因组中如SEQ ID No.2所示的核苷酸序列的最后一位到如SEQ ID No.2所示的核苷酸序列的最后一位的下游1000bp的区域中的序列，并且该反向引物的长度为20-30bp。

本发明的引物对是基于本发明所发现的DNA条形码所设计引物对，只要能够从食腺嘌呤节孢酵母TMCC 70007菌株基因组中扩增出所述DNA条形码即可。例如，引物对可以为用于扩增SEQ ID No.4或SEQ ID No.1所示的DNA条形码的正向引物和反向引物。基于给定的模板DNA设计引物是本领域的常规技术。

在更优选的实施方案中，所述正向引物和反向引物的核苷酸序列分别如下所示：

AraC-F：5’-CACTGCAACGGTCGCATAAG-3’(SEQ ID No.5)；

AraC-R：5’-GGTAACGGGCCACAGATTCA-3’(SEQ ID No.6)。

本发明的引物能实现对所述DNA条形码序列的特异性扩增。

本发明还提供了一种用于鉴别食腺嘌呤节孢酵母菌株的试剂盒，包含根据本发明所述的引物对。

在另一实施方案中，所述试剂盒还包含根据本发明所述的DNA条形码。所述DNA条形码可存在于记录介质上。该记录介质例如为光盘。

本发明的又一方面提供了一种用于鉴别食腺嘌呤节孢酵母菌株的方法，包括以下步骤：

a)提供待测菌株的基因组DNA；

b)以步骤a)所述的基因组DNA为模板，利用根据本发明所述的引物对进行PCR扩增，得到PCR产物；

c)电泳(例如琼脂糖凝胶电泳)检测PCR产物，如果没有目标条带，则判定待测菌株不是食腺嘌呤节孢酵母菌株TMCC 70007，如果有目标条带，则进行步骤d)；

d)对所得PCR产物进行测序，得到待测核苷酸序列；将所述待测核苷酸序列与本发明所述的DNA条形码的核苷酸序列进行同源性比对，若同源性在97％以上，则判定待测菌株是食腺嘌呤节孢酵母TMCC 70007菌株。

在本发明的实施方案中，术语“同源性”和“同一性”可以相互交换使用。

在本发明的具体实施方案中，所述PCR扩增的程序为：1)94℃-96℃预变性8分钟；2)94℃-96℃变性45秒，55℃-57℃退火45秒，72℃延伸1分15秒，其中该程序2)进行30-35个循环；3)72℃延伸10分钟。

在本发明的又一实施方案中，所述方法还包括将测序结果得到的待测序列与本发明的DNA条形码进行聚类分析(例如系统发育树)，如果待测序列与DNA条形码聚类在一起，则确定待测菌株为食腺嘌呤节孢酵母菌株TMCC 70007。

在本发明一个具体的实施方案中，利用本发明的引物对对提取自待鉴别的菌株的基因组DNA进行PCR扩增，然后进行琼脂糖凝胶电泳检测。基于检测是否存在PCR产物来鉴别菌株：如果待鉴定的菌株没有扩增出相应的目标条带，则说明该菌株不是TMCC 70007；如果扩增出相应的目标条带，则证明该菌株可能是TMCC 70007。为了进一步进行鉴别，对PCR产物进行测序，将DNA测序结果与DNA条形码序列进行同源性比对，获得序列间的相似性(即同源性)，如果序列同源性小于97％，则判定待测菌株不是食腺嘌呤节孢酵母菌株TMCC70007。如果序列同源性大于或等于97％，则判定待测菌株是食腺嘌呤节孢酵母菌株TMCC70007。

如果进行聚类分析，例如系统发育树，则将所述DNA条形码和各待鉴定菌株的DNA测序结果(即待测序列)一起应用MEGA 5.1或PAUP软件构建NJ系统发育树。如果待鉴定菌株的待测序列与食腺嘌呤节孢酵母TMCC 70007的DNA条形码聚类在一起，则鉴定为食腺嘌呤节孢酵母TMCC 70007。

在此所使用的术语“聚类”是指经系统发育树分析后，处于同一个分支，且进化距离相同。

本发明还提供了本发明所述的DNA条形码在鉴别食腺嘌呤节孢酵母菌株中的应用。

本发明还提供了本发明所述的引物对在鉴别食腺嘌呤节孢酵母菌株中的应用。

实施例

下面结合具体的实施例对本发明做进一步说明。实施例中所用方法在未具体说明的情况下，均采用常规方法和已知工具。

实施例1：AraC基因和DNA条形码的获得

1、利用高覆盖蛋白质组技术，采用pFind和pAnno软件对普洱茶工业发酵菌食腺嘌呤节孢酵母TMCC 70007进行了蛋白质组的深度覆盖研究，并对其基因组进行了注释编码基因的验证，发现了漏注释的基因AraC。为了发现新的蛋白编码区，利用系统蛋白质基因组学中的六框翻译(Six Frame Translation)策略，获得了食腺嘌呤节孢酵母TMCC 70007基因组数据的六框翻译数据库，穷举了基因组的6种编码可能(+1，+2，+3，-1，-2，-3)，并且将该核酸序列称为“六框翻译核酸序列”，蛋白序列称为“六框翻译蛋白序列”。通常而言，六框翻译核酸序列是从一个终止子到下一个终止子的序列，在本文中，也将其称为“蛋白编码框”。利用这个数据库，采用pFind和pAnno软件对质谱数据进行了新肽段和新蛋白质的鉴定。

经鉴定，发现了一条现有技术中食腺嘌呤节孢酵母TMCC 70007注释基因组中未发现的2条肽段LEGDYMEEWLR和TYAMSEIGAHFSSQR，质谱谱图如图1所示。

如图1所示，LEGDYMEEWLR和TYAMSEIGAHFSSQR二级质谱谱图(MS2)的匹配好，信号强，结果可信。

2、根据新肽段所在的位置，以前一个终止密码子和后一个终止密码子包括的区域为界，得到六框翻译核酸序列(蛋白编码框)SEQ ID NO.2。SEQ ID NO.2序列对应的mRNA及其所编码的氨基酸序列如图2所示，并且该氨基酸序列为SEQ ID NO.3。

3、为了进一步确定该编码基因(AraC基因)的编码起始位点和终止位点，将上述蛋白编码框向上游和下游分别扩展1000bp，采用GeneMark.hmm和Augustus进行基因预测，参考物种选用粟酒裂殖酵母(Schizosaccharomyces pombe)，在LEGDYMEEWLR和TYAMSEIGAHFSSQR两个肽段所处的蛋白编码框区域均没有预测到该蛋白编码基因的存在，这说明这个蛋白编码区不容易被当前的基因预测软件注释到。而采用蛋白质基因组学手段，直接鉴定基因表达产物蛋白质的肽段来系统揭示漏注释基因是一个切实可行的技术手段。

4、将AraC基因的蛋白编码框表达的蛋白序列(即SEQ ID NO.3)进行NCBI-BLASTP分析，发现该蛋白序列在NCBI数据库中没有明显相似的序列。经KEGG-BLAST检索分析后，发现其与其他物种(脱硫念珠菌(Desulfomonile tiedjei)DSM 6799)的AraC家族仅具有35％的相似度，但是该相似序列的序列太短，其为转录调控因子。

表1经KEGG-BLASTP分析后与AraC蛋白相似度较高的序列

AraC蛋白序列经过KEGG-BLASTP分析的同源性最高的2条序列的同源区段对比结果如图4所示。

5、将该AraC基因的蛋白编码框序列(即SEQ ID NO.2)进行NCBI-BLASTN和KEGG-BLASTN分析，未发现该基因的同源序列。这进一步佐证了该基因的蛋白编码框序列可以用于开发普洱茶发酵菌株食腺嘌呤节孢酵母TMCC 70007的潜在的DNA条形码序列。

6、为了得到长度更加合适并且具有足够特异性的DNA条形码，选择了SEQ ID NO.2上游1000bp到下游1000bp之间的区域(即序列SEQ ID NO.4，这个序列大于2000bp，如SEQID NO.4所示。

对SEQ ID NO.4进行NCBI-BLASTN分析，检索Nucleotide collection(nr/nt)，结果发现该序列仅有1％(约30多bp)与沙柳支孢瓶霉(Cladophialophora psammophila)CBS110553假定蛋白部分mRNA(登录号：XM_007742452.1)具有同源性(见表2)，说明该序列具有作为DNA条形码的潜力。对SEQ ID NO.4进行NCBI-BLASTN分析，检索全基因组鸟枪法序列(Whole-genome shotgun contigs，WGS)数据库，食腺嘌呤节孢酵母(taxid:409370)，结果发现其与食腺嘌呤节孢酵母LS3菌株全基因组序列(登录号：CBZY010000003.1)中部分序列具有同源性，同一性为84％(见表3)。所以尽管TMCC 70007与LS3菌株属于同一物种，但在这个区段上，序列的差异已经足够将TMCC 70007鉴别开来，这进一步说明该序列可以作为特异性鉴别该菌株的DNA条形码序列。

表2

表3

综上所述，SEQ ID NO.2上游1000bp到下游1000bp之间的序列(SEQ ID NO.4)可以作为鉴别食腺嘌呤节孢酵母菌株的DNA条形码。

7、进一步地，基于SEQ ID NO.4，采用NCBI引物设计工具设计PCR引物，扩增产物须包括蛋白编码框(SEQ ID NO.2)序列。由此得到一系列适合鉴别食腺嘌呤节孢酵母菌株的DNA条形码。其中一对正反引物序列分别为：

AraC-F：5’-CACTGCAACGGTCGCATAAG-3’；

AraC-R：5’-GGTAACGGGCCACAGATTCA-3’。

该对引物理论扩增长度为1262bp，引物位置及扩增序列见图3所示，所扩增的DNA序列如SEQ ID NO.1所示。从长度和特异性方面考虑，并且经过下面的验证，其为用于鉴别食腺嘌呤节孢酵母菌株的最佳DNA条形码序列。

8、对SEQ ID NO.1序列进行NCBI-BLASTN分析发现，该DNA序列仅有一小部分序列仅有一部分序列与沙柳支孢瓶霉CBS 110553菌株的假定蛋白部分mRNA的30bp左右的序列同源性较高，但同源序列较短，不足条形码序列SEQ ID NO.1的3％，如表4及图5所示，这表明我们在食腺嘌呤节孢酵母TMCC 70007检测到的这个DNA条形码序列特异性。

表4条形码序列SEQ ID NO.1序列NCBI-BLASTN结果

9、目前，食腺嘌呤节孢酵母物种中，仅有食腺嘌呤节孢酵母LS3菌株的基因组已经被报道(Kunze et al.Biotechnology for Biofuels2014,7:66)。经Local-BLASTN分析后，发现该基因序列SEQ ID NO.1在LS3菌株基因组中存在同源序列。经DNAMAN分析后，两条序列的同源性为79.03％，尽管TMCC 70007与LS3菌株为同一个物种(即食腺嘌呤节孢酵母)，但在该DNA条形码序列的1262个碱基中，约有250个位点是不一样的，结果如图6所示。这进一步证明了利用该DNA条形码序列可以有效地区别TMCC 70007菌株与种内的其他菌株。

实施例2.利用DNA条形码鉴别菌株

依据待测样品的扩增结果和与食腺嘌呤节孢酵母TMCC 70007菌株SEQ ID NO.1的序列同源性，来判定待测样品是否为普洱茶工业应用菌株食腺嘌呤节孢酵母TMCC 70007。

(1)菌株来源

表5选用的相关菌株信息

(2)分别提取菌株DNA：采用OMEGA E.Z.N.A.^TM的Yeast DNA试剂盒提取酵母基因组DNA，并用灭菌的去离子水将样品的DNA浓度稀释到0.5μg/μL。

(3)扩增DNA片段，进行聚合酶链式(PCR)反应，所用引物序列分别为：

AraC-F：5‘-CACTGCAACGGTCGCATAAG-3’；

AraC-R：5’-GGTAACGGGCCACAGATTCA-3’。

PCR反应体系为50μL，PCR试剂为Thermo Scientific^TM的Taq DNA Polymerase(重组)：ddH₂O 37.7μL、MgCl₂ 5μL，dNTPs 4μL，正向引物1μL，反向引物1μL，Taq DNA聚合酶0.3μL、DNA模板1μL，不含染料。扩增程序为：94℃预变性8分钟；接下来94℃变性45秒，56℃退火45秒，72℃延伸1分15秒，一共进行32-35个循环；最后72℃延伸10分钟。

(4)扩增产物检测：用1.0％的琼脂糖凝胶、1×TBE电泳液电泳检测，使用DNAmarker检测PCR片段大小。若待测菌株没有扩增条带，说明该菌株不是食腺嘌呤节孢酵母TMCC 70007；若出现明显清晰的条带，且无杂带，送生物测序公司进行DNA片段测序。

(5)该引物仅能在食腺嘌呤节孢酵母TMCC 70007中实现扩增，在同物种的食腺嘌呤节孢酵母CBS 8244^T、CBS 7350、CBS 7370、CBS 8335和Blastobotrysraffinosifermentans CBS 6800^T中则不能实现扩增。结果如图7所示。该PCR引物的理论扩增序列为1262bp，和预期相符。为了进一步验证扩增的DNA的序列，进行了测序和同源序列比较。

(6)对于有条带的序列的测序结果，首先用软件Chromas检查测序后得到的序列峰图质量，确定峰图质量达到数据分析的要求之后，用DNASTAR软件包中的SeqMan进行正反向序列的拼接。将测序结果进行手工校对、序列拼接。如果待测菌株DNA片段与食腺嘌呤节孢酵母TMCC 70007标准DNA条形码的同源性在97％以上，即可判断所述的待测菌株可能为食腺嘌呤节孢酵母TMCC 70007菌株。例如，经过DNAMAN，将菌株TMCC 70007测序得到序列与DNA条形码序列SEQ ID NO.1比对后完全一致，进一步证明了这种方法的可靠性。结果如图8所示。

SEQUENCE LISTING

<110> 勐海茶业有限责任公司

<120> DNA条形码、引物、试剂盒、方法和应用

<130> FI-164053-59:52/C

<160> 6

<170> PatentIn version 3.5

<210> 1

<211> 1262

<212> DNA

<213> 食腺嘌呤节孢酵母（Blastobotrys adeninivorans）

<400> 1

cactgcaacg gtcgcataag cgatgctaaa aaaaaagatt taggcatcag ttaacgccct 60

agattaaatt tagccacatt agatgcgagt aatggtggct cgaggaatgg gttcaccccc 120

tttgcaaaca agacttagac cctgaccacc gtctccacgg gggctgccct ccaattagct 180

cattagctca ggtgagctcg ctccgccatg gtgccgctta gccactgagc caagccactt 240

ttaggtgcct ggttcctgaa caagcctgct aagaaaatca gtggagaagg atgtgtcgag 300

ggcctgaaag acgtcatgcg ctgtcgatca gcgggttccc ctacggcaaa gggcttatac 360

gggtgcgacg catggttaca actataaata attgccaatt ccgaggacca ccaccactta 420

ttgacaatgg tcgaggccca gacgacatgt gcggtttgtg aatggtgcga agggaacctg 480

tcgtatccta cgggaaagac atatgcaatg agcgagattg gtgctcactt tagctctcag 540

cggcatgtta gcaattggtt tacgtgggta aagtatctac gggagcaaag actagagggg 600

gattacatgg aggagtggct tcggcagcat ggtctatctc gagttcatca tcaggcgatc 660

aaaggaagga ggggtcatta gcactcttac cattaggtgc tagtgatttg ttttgccttt 720

tcgtttttgt tttcatttgc atttgcattg acgccctcat ttttattttc attttcgcct 780

ttgtctttat gatcctcatg atcctcacga tcctcattag cctcattggc ctcattagcc 840

tcattggagc tattggaacc attggtgcca ttggagccac tggagtcacg atcctcattg 900

gccttgttgg agccatttga gtctgcagtg gttctggtga ggtacacatt ctctattgca 960

tttttgttgg cagttttgtc ggtacgttct gccaattcta acaacctggc ttcttcctga 1020

agacgcttgg tcgacatcca tcctcttccg gggatgaaca agcggcggta ttgcttgcgc 1080

tccccttctt caagcttggc tctgatttca ctactggtta attgtacggg aggtctttcg 1140

tcccctggtt gggcactagt gctggtacta ctagatatgg cctccagttt agcacccaat 1200

tgtgtaccat gagtgtagga aaaagatgtc ttttcttgaa cttgaatctg tggcccgtta 1260

cc 1262

<210> 2

<211> 297

<212> DNA

<213> 食腺嘌呤节孢酵母（Blastobotrys adeninivorans）

<400> 2

taaataattg ccaattccga ggaccaccac cacttattga caatggtcga ggcccagacg 60

acatgtgcgg tttgtgaatg gtgcgaaggg aacctgtcgt atcctacggg aaagacatat 120

gcaatgagcg agattggtgc tcactttagc tctcagcggc atgttagcaa ttggtttacg 180

tgggtaaagt atctacggga gcaaagacta gagggggatt acatggagga gtggcttcgg 240

cagcatggtc tatctcgagt tcatcatcag gcgatcaaag gaaggagggg tcattag 297

<210> 3

<211> 97

<212> PRT

<213> 食腺嘌呤节孢酵母（Blastobotrys adeninivorans）

<400> 3

Ile Ile Ala Asn Ser Glu Asp His His His Leu Leu Thr Met Val Glu

1 5 10 15

Ala Gln Thr Thr Cys Ala Val Cys Glu Trp Cys Glu Gly Asn Leu Ser

20 25 30

Tyr Pro Thr Gly Lys Thr Tyr Ala Met Ser Glu Ile Gly Ala His Phe

35 40 45

Ser Ser Gln Arg His Val Ser Asn Trp Phe Thr Trp Val Lys Tyr Leu

50 55 60

Arg Glu Gln Arg Leu Glu Gly Asp Tyr Met Glu Glu Trp Leu Arg Gln

65 70 75 80

His Gly Leu Ser Arg Val His His Gln Ala Ile Lys Gly Arg Arg Gly

85 90 95

His

<210> 4

<211> 2292

<212> DNA

<213> 食腺嘌呤节孢酵母（Blastobotrys adeninivorans）

<400> 4

actaccgtac catcttccga atcttcgtcc gactggtact ctggctgatc ctctggctgg 60

tcctccgaat ggtcctctga ctgatcctcc gcatggtcgt caccttcggt atcgcggtaa 120

agctttatca ctccatggcc aagatacctt gatttgagct ggcccatagg accgctggga 180

tcttctgtca tgtcaagagg aaatacgtgg ataggccctg cacgccaatc cttggacgta 240

tggaggttca cgtactctaa cggtagcgat ttggctggtt caaactggcc tttacggaag 300

cgatggtggt ctggtctagg aagagatgcc cctttagcgc gctcaggagt aaaaatatcg 360

tcgtggttca cggccgcaac aatccagggc tcggggtcgt ccagtatgag gtcgataaca 420

atcgtatagc tgccagaagc agcttgcatc accgatcaat atccagtcaa tgtccaggag 480

aatgattaga ttgacgatag atgataaggt gtaaacgtga tctttatcag ctgatgataa 540

tcatttgata gtatatcagg tgattgcacg gctcctcaac ggttttgttg gtgcagcagc 600

cgtgtgacag acgtcactgc aacggtcgca taagcgatgc taaaaaaaaa gatttaggca 660

tcagttaacg ccctagatta aatttagcca cattagatgc gagtaatggt ggctcgagga 720

atgggttcac cccctttgca aacaagactt agaccctgac caccgtctcc acgggggctg 780

ccctccaatt agctcattag ctcaggtgag ctcgctccgc catggtgccg cttagccact 840

gagccaagcc acttttaggt gcctggttcc tgaacaagcc tgctaagaaa atcagtggag 900

aaggatgtgt cgagggcctg aaagacgtca tgcgctgtcg atcagcgggt tcccctacgg 960

caaagggctt atacgggtgc gacgcatggt tacaactata aataattgcc aattccgagg 1020

accaccacca cttattgaca atggtcgagg cccagacgac atgtgcggtt tgtgaatggt 1080

gcgaagggaa cctgtcgtat cctacgggaa agacatatgc aatgagcgag attggtgctc 1140

actttagctc tcagcggcat gttagcaatt ggtttacgtg ggtaaagtat ctacgggagc 1200

aaagactaga gggggattac atggaggagt ggcttcggca gcatggtcta tctcgagttc 1260

atcatcaggc gatcaaagga aggaggggtc attagcactc ttaccattag gtgctagtga 1320

tttgttttgc cttttcgttt ttgttttcat ttgcatttgc attgacgccc tcatttttat 1380

tttcattttc gcctttgtct ttatgatcct catgatcctc acgatcctca ttagcctcat 1440

tggcctcatt agcctcattg gagctattgg aaccattggt gccattggag ccactggagt 1500

cacgatcctc attggccttg ttggagccat ttgagtctgc agtggttctg gtgaggtaca 1560

cattctctat tgcatttttg ttggcagttt tgtcggtacg ttctgccaat tctaacaacc 1620

tggcttcttc ctgaagacgc ttggtcgaca tccatcctct tccggggatg aacaagcggc 1680

ggtattgctt gcgctcccct tcttcaagct tggctctgat ttcactactg gttaattgta 1740

cgggaggtct ttcgtcccct ggttgggcac tagtgctggt actactagat atggcctcca 1800

gtttagcacc caattgtgta ccatgagtgt aggaaaaaga tgtcttttct tgaacttgaa 1860

tctgtggccc gttacccata gtatttcgta tggtgcttcc cgcagaggcc cctacggtat 1920

ttcgtatggt ggtccctggg gtgttccctg gggtggttcc tacggtattt cgtgtggtgg 1980

tggtattggt atcaattgca aatggatttt tataattcac atctcgaggg ctgtgagggc 2040

tgtgcgaaga ggaggaatcc ttgtttaaca gtttccttat agattccatc catctctgat 2100

ggccctctct tgtgccaatg tcgtacactc cgtgctcagt tctcagaatg ttgctgcgct 2160

tatgtctggt aaagttgggc ttgtatacat ctgaattgtt ggactggaac tggctcattg 2220

gcgggggata aaagtcatac atttgctggc ccaattgagg gattcgcgac agtcgatcag 2280

atacctgaaa ta 2292

<210> 5

<211> 20

<212> DNA

<213> 人工序列

<400> 5

cactgcaacg gtcgcataag 20

<210> 6

<211> 20

<212> DNA

<213> 人工序列

<400> 6

ggtaacgggc cacagattca 20

Claims

1.一种用于扩增用于鉴别食腺嘌呤节孢酵母TMCC 70007菌株的DNA条形码的引物对，所述正向引物和反向引物的核苷酸序列分别如下所示：

正向引物：5’-CACTGCAACGGTCGCATAAG-3’；

反向引物：5’-GGTAACGGGCCACAGATTCA-3’；

其中，所述DNA条形码来源于食腺嘌呤节孢酵母TMCC 70007菌株的基因组，并且核苷酸序列选自如SEQ ID No.4所示的序列并且包含如SEQ ID No.1所示的序列。

2.一种用于鉴别食腺嘌呤节孢酵母TMCC 70007菌株的试剂盒，包含根据权利要求1中任一项所述的引物对。

3.一种用于鉴别食腺嘌呤节孢酵母TMCC 70007菌株的方法，包括以下步骤：

a)提供待测菌株的基因组DNA；

b)以步骤a)所述的基因组DNA为模板，利用根据权利要求1所述的引物对进行PCR扩增，得到PCR产物；

c)电泳检测PCR产物，如果没有目标条带，则判定待测菌株不是食腺嘌呤节孢酵母TMCC70007菌株，如果有目标条带，则进行步骤d)；

d)对所得PCR产物进行测序，得到待测核苷酸序列；将所述待测核苷酸序列与所述的DNA条形码的核苷酸序列进行同源性比对，若同源性在97％以上，则判定待测菌株是食腺嘌呤节孢酵母TMCC 70007菌株。

4.根据权利要求1所述的引物对在鉴别食腺嘌呤节孢酵母TMCC 70007菌株中的应用。