摘要本文对3种不同来源样品的4株椰酵假单胞菌的16S rDNA的序列进行了测定，以确定其系统发育位置及与其它菌种的遗传进化关系。实验中设计了7条引物，利用聚合酶链式反应(Polymase chain Reaction,PCR)扩增4株细菌的16S rDNA，通过扩增产物的克隆测序或直接测序，获得了它们的16S rDNA序列，将测得的序列与伯克霍尔德氏菌属中其它菌种的16S rDNA序列进行了聚类分析，得到了系统发育进化树状图。树状图及不同菌种之间的同源性比值：椰酵假单胞菌(P.cocovenenans subsp.farinofermentans)与唐菖蒲伯克霍尔德氏菌(Burkholderia gladioli)、椰毒伯克霍尔德氏菌(Burkholderia cocovenenans)的亲缘关系非常近，与其它菌种的亲缘关系相对较远，椰酵假单胞菌、唐菖蒲伯克霍尔德氏菌和椰毒伯克霍尔德氏菌可归属于一个独立的系统发育系。

中图分类号R155.31TS 207.4Q939.110.1

Squencing and analysis of 16S rDNA sequences for

P.cocovenenans subsp farinofermentans

Jiao Zhenquan,Liu Xiumei,Yang Ruifu,Meng Zhaohe

institute of Nutrition and Food Hygiene,Chinese Academy of Preventive Medicine beijing100050,China

in order to fix the phylogenetic taxonomic position of p.cocovenenans subsp.farinofermentans and describe its relationships with other closed species,seven primers were designed to amplify and sequence the 16S rDNA of 4 strains of it by using a clonal sequencing or direct sequencing method.The16S rDNA sequence of the 4 strains were aligned with the 16S rDNA sequences of other species in genus Burkholderia,and the polygenetic tree was produced by using a Clustal V and treeview software.The results showed that P.cocovenenans subsp.farinofermentans was in a high homology with Burkholderia gladioli and burkholderia cocovenenans， and they formed an independent phylogenetic clustal of genus Burkholderia.

Key words: P.cocovenenans subsp.farinofermentans,16S rDNA sequcence,phylogeny

椰毒假单胞菌酵米面亚种(Pseudomonas cocovenenans subsp.farinofermentans,简称椰酵假单胞菌)是我国学者于1977年在东北酵米面中毒食物中发现的一种新的食物中毒菌(曾暂命名为酵米面黄杆菌^［1］)。米酵菌酸和毒黄素是该菌产生的两种主要毒素，其中前者是引起食物中毒和死亡的主要毒性代谢产物，肝、脑、肾等主要实质性脏器是该毒素作用的靶器官。椰毒假单胞菌酵米面亚种可污染鲜银耳、酵米面、醋凉粉、汤圆粉和马铃薯粉等多种食品而变质发生食物中毒^［2］。

近10多年来，我们研究室在有关椰酵假单胞菌的菌体抗原结构，生化鉴定和检测方法及米酵菌酸结构的鉴定、产毒机制和消毒，去毒等方面进行了广泛和深入的研究。刘秀梅在美国ATCC(America type Culture Collection)所做的研究结果显示：椰假单胞菌在生物学性状上与唐菖蒲伯克霍尔德氏菌更为接近。邱茂锋等在利用rDNA指纹图对椰酵假单胞菌、椰毒伯克霍尔德氏菌和唐菖蒲伯克霍尔德氏菌进行对比研究时发现：唐菖蒲伯克霍尔德氏菌的图谱谱形与椰酵假单胞菌的RT9型相同，而椰毒伯克霍尔德氏菌与椰酵假单胞菌在谱形上则存在差别，缺乏椰酵假单胞菌所共有的唯一杂交带，说明椰酵假单胞菌与唐菖蒲伯克霍尔德氏菌的亲缘关系更为接近^［3］。

近几年来，细菌的分类学发展迅速，已经从单一的表型分类发展为系统分类、多相分类；即利用表型、遗传型和系统发育信息来区分细菌，进而阐述细菌的遗传进化关系。自从Woese根据16S rRNA序列将所有生命分成3个域后^［4］，作为细菌分类的系统发育框架基础16S rRNA，以其作为“细菌化石”的高度保守特点而被越来越多的分类学家所重视，它使人们可以进一步了解细菌的生命进化关系。

本研究的目的就是分析来源于3种不同样品的椰酵假单胞菌的16S rDNA序列，以确定它与椰毒伯克霍尔德氏菌和唐菖蒲伯克霍尔德氏及与伯克霍尔德氏菌属中其它菌种的遗传进化关系。

1材料和方法

1.1实验菌株

本实验选用椰毒假单胞菌酵米面亚种(Pseudomonas cocovenenans subsp.farinofermentans)的4个样品中的分离株进行16S rDNA序列测定：HN2y(CCTCCAB9725)1989年于银耳中分离，Co14 于1977年酵米面中分离，Sx8801于1987年由醋凉粉中分离，90-3于1990年由酵米面中分离。

1.2实验方法

1.2.1DNA的制备

菌株总DNA的提取：菌株用LB培养基培养至对数生长期，以8000r/min离心收集菌体，用SDS溶液进行裂解，再用酚/氯仿/异戊醇抽提和RNA酶消化，最后用乙醇沉淀，溶于TE(10mmol/L tris-HCl,pH7.6及1mmol/L EDTA,pH8.0)中。

1.2.216S rRNA引物的设计

根据16S rRNA高度保守的特点，从其保守区域寻找引物位点，利用Oligo4.1软件设计引物，共合成7条引物。

正向引物：

primer 1:5'-AGA GTT TGA TCA TGG CTC AG-3'(对应于大肠杆菌的第8—27bp)，

primer 2:5'-GTA GTC CAC GCT GTA AAC GAT G-3'(对应于大肠杆菌的第774—795bp)反向引物：

primer 6:5'-CTA CGT TTA CAG CGT GGA CTA C-3'(对应于大肠杆菌的第795—774bp)

primer 7:5'-TAC GGT TAC CTT GTT ACG AC-3'(对应于大肠杆菌的第1504—1485bp)

其它的3条引物序列分别为：

primer 3:5’-CTA CGG GAG GCA GCA GTG GG-3'

primer 4:5'-CAG CAG CCG CGG TAA T-3'

primer 5:5'-AAA CTT AAA GGA ATT GAC GG-3'。

1.2.316S rRNA基因扩增及序列测定

以总的DNA为模板，用引物P1，P7经PCR反应扩增出总的16S rDNA，PCR反应条件为：预变性94℃ 3min；变性94℃ 1min，退火55℃1min，复性72℃ 1min；最后延伸72℃ 10min，反应35个循环。对于HN2y和Co14,以总的16S rDNA为模板，用引物P1，P6及P2，P7进行PCR扩增，扩增产物用PGEM-T载体进行钝端连接后转化大肠杆菌JM109，经蓝白斑筛选及质粒大小比较后，以T7和SP6为测序引物利用A₃₇₃ dNA自动测序仪进行测序。对于Sx8801和90-3则用总的16S rDNA为测序模板，用P1，P3，P4，P2及P5为测序引物在A₃₇₃ dNA自动测序仪直接进行测序。

1.2.416S rDNA序列分析

将伯克霍尔得氏菌属(Burkholderia)中其它菌种的16S rDNA序列以GenBank/EMBL/DDBJ核酸序列库中检索出来^［5］，它们包括唐菖蒲伯克霍尔德氏菌(B.gladioli)、椰毒伯克霍尔德尔菌(B.cocovenenan)、洋葱伯克霍尔德氏菌(B.cepacia)、须芒草伯克霍尔德氏菌(B.andropogonis)、石竹伯克霍尔德氏菌(B.cariophylli)、植物伯克霍尔德氏菌(B.plantarii)、范氏伯克霍尔德氏菌(B.vandii)、越南伯克霍尔德氏菌(B.vienamiensis)、范氏伯克霍尔德菌(B.vandii)、越南伯克霍尔德尔菌(B.vienamiensis)、吡咯菌素伯克霍尔德氏菌(B.pyrrocinia)、鼻疽伯克霍尔德氏菌(B.mallei)、类鼻疽伯克霍尔德氏菌(B.pseudomallei)、泰国伯克霍尔德氏菌(B.thailandensis)、格氏伯克霍尔德氏菌(B.glathei)、吩嗪伯克霍尔德氏菌(B.phenazinium)荚壳伯克霍尔德菌(B.glumae)及B.graminis,将这些菌株的16S rDNA序列与测得的四株椰酵假单胞菌的16S rDNA序列利用Clustal Ⅴ软件进行排列^［6］，选用1307bp(相当于大肠杆的145-1452bp位置)进行比较，选取大肠杆菌(Escherchia coli,E.coli)作为系统发育进化树的根部，用Neighbor-Joining方法构建系统发良进行树状图^［7］，用Treeview3.2软件绘制系统发育进化树。各菌株的16S rDNA基因序列的核酸序列库中的索取号分别是：

b.gladioli X67038,B.cocovenenan U96934,B.cepacia U96927,B.andropogoni X67037,B.cariophylli X 67039,B.plantarii U96933,B.vandii U96932,B.vienamiensis u96928,B.pyrrocinia U96930,B.mallei,B.pseudomallei BPU91839,B.thailandensis bSU91838,B.glathei Y17052,B.phenaxinium U9636,B.glumae U96931,B.graminis u96940,E.coli J01859.

2结果与讨论

2.1序列测定

我们利用8条引物，应用PCR扩增产物克隆后测定及直接测序两种方法获得了椰酵假单胞菌HN2y,Co14,90-3和Sx8801四个菌株的16S rDNA序列，大小分别为1484bp,1499bp,1352bp和1369bp,序列见GenBank。

2.2序列分析将椰酵假单胞菌的4株菌的16S rDNA序列与伯克霍尔得氏菌属中其它菌种的16S rDNA序列利用Clustal V软件进行比较、聚类，用Treeview软件获得了系统发育进化树状图，如附图所示。不同菌种之间的同源性比值及核苷酸替代率Knuc见表2^［8］