云南小驳骨丛枝植原体的分子鉴定及相关基因序列分析
Molecular Identification and Related Gene Sequence Analysis of Phytoplasma of Gendarussa vulgaris Witches’-broom in Yunnan Province
-
Keywords:
-
gendarussa vulgaris nees witches’-broom disease / - phytoplasma /
- 16S rRNA gene /
- rp gene /
- secY gene /
- secY protein
-
植原体(phytoplasma)是一类非螺旋状、无细胞壁且专性寄生于韧皮部的植物病原微生物。植原体主要依靠叶蝉和飞虱等媒介昆虫进行传播，亦可由植物营养繁殖材料和菟丝子等方式进行传播，感病植株通常表现为丛枝、花变叶、黄化、矮化和节间缩短等症状，目前植原体可危害1 000多种植物，给多种经济作物的生产造成了严重影响[1-3]。
植原体的系统发育进化关系及分类鉴定主要依赖于高度保守的16S rRNA基因序列的遗传变异性[4-5]。目前基于16S rRNA基因序列RFLP分析的相似系数对植原体株系进行划分，已确定了36个16Sr组和150多个亚组[6-7]。然而在植原体分类鉴定体系中，以高度保守的16S rDNA序列信息作为唯一的系统发育依据，在鉴定亲缘关系密切的植原体株系上存在缺陷且难以对组内植原体株系进行更精确的分类鉴定[8]。为解决植原体分类学的局限性，研究人员开始利用保守性相对较低的非核糖体基因，如核糖体蛋白(rp)基因和secY基因等作为辅助分类的分子遗传标记，有助于更准确地对植原体16Sr组内株系的真实系统发育关系进行更精细的划分，可进一步将植原体16Sr组划分为多个不同的亚组[9-12]。
小驳骨(Gendarussa vulgaris Nees)又名驳骨草、接骨草等，属于双子叶植物纲唇形目爵床科(Acanthaceae)驳骨草属(Gendarussa Nees)，具有极高的药用价值[13-14]。2019年，课题组在元谋地区发现的疑似感染植原体而诱发的丛枝症状小驳骨病株，其植株生长缓慢、矮化、节间缩短、腋芽丛生为簇状并且大部分叶片萎缩变小并黄化，这将会降低小驳骨的产质量和药用价值。
基于作为系统发育分子标记的rp和secY基因具有能够高度准确地辨别16Sr组中亚组间的系统发育关系的能力，同时利用多个系统发育参数可以更好地识别和定义植原体株系。本研究将采用PCR技术等研究方法对小驳骨植原体16S rRNA、rp及secY基因进行扩增及序列测定，根据这3个基因的序列信息构建系统发育树进行系统发育关系的比较分析，最终从亚组水平上准确高效地鉴定小驳骨丛枝植原体株系的分类地位。同时通过对secY蛋白亚基的蛋白特性、蛋白跨膜结构及信号肽进行分析，从蛋白分泌途径初步了解植原体病害的致病机理，为研究小驳骨植原体病害的发生及防治提供理论支持。
1. 材料与方法
1.1 试验材料与主要试剂
自然表现丛枝症状的小驳骨植株采自云南省元谋县，并保存于实验室，并以实验室保存的苦楝丛枝植原体病株为阳性对照。
新型快速植物基因DNA提取试剂盒及PCR产物纯化试剂盒由百泰克生物技术有限公司提供；DNA Marker、TransStart® FastPfu DNA Polymerase和克隆载体试剂盒pEASY®-Blunt Simple Clonging Kit均来自北京全式金生物技术有限公司；引物由硕擎生物科技有限公司合成。
1.2 方法
1.2.1 感病小驳骨植株总DNA提取和PCR扩增
选取感病小驳骨植株幼嫩叶片进行总DNA提取，提取总DNA步骤参照新型快速植物基因DNA提取试剂盒说明书。
以小驳骨病样总DNA为模板，参照PÉREZ-LÓPEZ等[7]和LEE等[8]设计的植原体通用引物对P1/P7和R16F2n/R16R2对植原体16S rRNA基因通过巢式PCR进行扩增，同时采用引物对rp (II) F1/rp (II) R1、secY-F/secY-R分别对 rp和secY基因进行直接PCR扩增。引物合成序列见表1。16S rRNA及rp基因的PCR扩增反应条件参考万琼莲等[15]的研究；secY基因的PCR扩增反应条件为95 ℃ 2 min，95 ℃ 20 s，45 ℃ 20 s，72 ℃ 40 s，35 cycle，72 ℃延伸 5 min。
表 1 引物序列表Table 1. Primer sequence listing基因 gene 引物名称 primer name 序列 sequence 预期扩增片段长度/bp expected amplified fragment length 16S rRNA P1 5′-AAGAGTTTGATCCTGGCTCAGGATT-3′ 1 800 P7 5′-CGTCCTTCATCGGCTCTT-3′ R16F2n 5′-GAAACGACTGCTAAGACT-3′ 1 248 R16R2 5′-TGACGGGCGGTGTGTACAAACCCCG-3′ rp rp(II)F1 5′-GCTCTTACTCGTAAAYATGTAGT-3′ 1 171 rp(II)R1 5′-TTACTTGATTTTCTGGTTTTGA-3′ secY secY-F 5′-GCGGAAGAAGCTATTAT-3′ 1 425 secY-R 5′-CGAATAACATAATATAATTGATTCC-3′ 1.2.2 16S rRNA、rp和secY基因序列的克隆及测定
PCR扩增产物在1%琼脂糖凝胶中进行电泳检测，对PCR产物进行纯化，并与载体pEASY®-Blunt Simple Cloning Vector连接转化到化学感受态细胞Trans1-T1 Phage Resistant中，37 ℃培养1 h后取100 µL均匀涂布于经IPTG和X-gal处理过的LB培养基平板上，过夜培养，第2天进行蓝白斑的挑选以及菌液PCR的阳性克隆子检测，将含有目的基因片段的菌液送生工测序公司进行序列测定。
1.2.3 16S rRNA、rp及secY基因序列分析
将测序所获得的16S rRNA、rp及secY基因的核苷酸序列用Vector NTI 11.5.3.软件去除序列两端的载体，以确定测序结果的准确性。利用NCBI中的ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/)寻找rp和secY序列的最大开放阅读框并进行分析。将序列提交到NCBI中BLAST比对确认，并进行同源序列检索；采用MUSCLE(https://www.ebi.ac.uk/Tools/msa/muscle/)对序列进行在线同源性比对，分析基于16SrII组各植原体16S rRNA、rp及secY基因序列的亲缘关系。利用MEGA X软件，以邻接法(neighbour-joining)构建分别基于16S rRNA、rp及secY基因的系统进化树以便进行系统发育关系的比较分析，确定小驳骨丛枝植原体株系的分类地位。16S rRNA基因序列通过植原体在线分析网站iPhyClassifier (https://plantpathology.ba.ars.usd a.gov/cgibin/resource/iphyclassifier.cgi)进行虚拟RFLP分析鉴定，明确该株系的候选种分类地位。相关植原体株系信息如表2。
表 2 植原体各株系16S rRNA、rp和secY基因序列登录号Table 2. The accession No. of 16S rRNA, rp gene and secY gene sequence of phytoplasma strain植原体病害及相关菌株
phytoplasma diseases and strains地理来源
geographical origin16Sr 组-亚组
16Sr group-subgroupGenBank登录号 GenBank accession No. 16S rRNA secY rp 云南元谋番茄巨芽植原体
tomato big bud, TBB-YM/ TBBP-YNym中国云南
Yunnan, China16SrII-A JQ923436 KC953016 KC953010 花生丛枝植原体
peanut witches’-broom, PnWB中国台湾
Taiwan, China16SrII-A L33765 GU004331 EF193375 云南元谋番茄丛枝植原体
tomato witches’-broom, TWB-YNym中国云南
Yunnan, China16SrII-A KC953005 KC953017 KC953011 云南元谋四季豆花变绿植原体
string bean virescence, SbV-YNym2中国云南
Yunnan, China16SrII-A KJ735766 — — 芝麻花变叶植原体
Candidatus phytoplasma aurantifolia, SEPB2泰国
Thailand16SrII-A JN006079 — — 云南元谋兵豆花变绿植原体
‘Lens culinaris’ virescence, LeVP-YNym中国云南
Yunnan, China16SrII-A — KJ735763 KJ735764 云南元谋豇豆花变叶植原体
cowpea phyllody, CowP-YNym中国云南
Yunnan, China16SrII-A KC953001 — — 云南元谋豇豆花变绿植原体
cowpea virescence, CoVP-YNym中国云南
Yunnan, China16SrII-A — KC953013 KC953007 黑豆丛枝植原体
black gram witches’-broom, BgWB缅甸
Myanmar16SrII-A AB690304 AB703249 AB703243 黄瓜花变叶植原体
cucumber phyllody, CuPh伊朗
Iran16SrII-A KY349111 KY365527 KY365523 云南元谋菊苣从枝植原体
‘Cichorium intybus’ witches’-broom,
CiWB-YNym/ChWB-YNym中国云南
Yunnan, China16SrII-A KJ735774 KJ735775 — 柠檬从枝植原体
lime witches’-broom, LWB阿曼
Oman16SrII-B EF186828 — — 棉花花变叶植原体
cotton phyllody, CoP布基纳法索
Burkina Faso16SrII-C EF186827 GU004350 EF186814 猪屎豆花变叶植原体
Crotalaria phyllody, CrP泰国
Thailand16SrII-C EF193355 GU004349 EF186818 大豆花变叶植原体
soybean phyllody, SOYP泰国
Thailand16SrII-C EF193353 GU004324 EF186816 刺缘毛莲菜花变叶植原体
Picris echioides phyllody, PEP意大利
Italy16SrII-E Y16393 GU004348 EF193381 桃黄化植原体
peach X-disease，CX加拿大
Canada16SrIII-A L33733 GU004327 EF186813 三叶草黄化植原体
clover yellow edge, CYE立陶宛
Lithuania16SrIII-B AF173558 GU004332 JQ360962 核桃丛枝植原体
walnut witches’-broom, WWB美国乔治亚州
Georgia, USA16SrIII-G AF190227 GU004325 JQ360961 榆树黄化植原体
elm yellows, ‘Ca.P.ulmi’EYEu (EY-626、EY-627)意大利
Italy16SrV-A AY197657 AY197691 AY197678 悬钩子矮化植原体
rubus stunt, RUS意大利
Italy16SrV-E AY197648 AY197696 AY197668 马铃薯丛枝植原体
potato witches’-broom, ‘Ca.P.trifolli’PWB加拿大
Canada16SrVI-A AY500818 GU004316 AY197683 白蜡木黄化植原体
ash yellows, AshY1美国纽约
New York, USA16SrVII-A AF092209 GU004329 EF183492 欧洲核果黄化/李子坏死植原体
European stone fruit yellows/ plum, ‘Ca.P. prunorum’PLN德国/意大利
Germany/Italy16SrX-B AJ542545 GU004334 EF193369 绣球花花变叶植原体
hydrangea phyllody, HYDP比利时
Belgium16SrI-A AY265215 AY803181 AY264868 报春花花变绿植原体
primrose virescence, PRIVC德国
Germany16SrI-B AY265210 AY803176 AY264854 三叶草花变叶植原体
clover phyllody, CPh加拿大
Canada16SrI-C AF222065 AY803179 AY264862 蓝莓矮化植原体
blueberry stunt, BBS3美国
USA16SrI-E AY265213 GU004357 AY264863 1.2.4 secY基因编码的蛋白亚基的蛋白特性分析及结构预测
分别通过生物信息学在线分析网站ProtParam (http://web.expasy.org/protparam/)、Protscale (https://web.expasy.org/protscale/)、TMHMM (http://www.cbs.dtu.dk/services/TMHMM-2.0/)和SignalP (http://www.cbs.dtu.dk/ services/SignalP/)对secY基因编码的蛋白特性、蛋白亲疏水性、蛋白跨膜结构区和信号肽进行分析。
2. 结果与分析
2.1 小驳骨丛枝病症状
由图1所示：感病小驳骨植株严重矮化，节间明显缩短，侧芽细弱，枝条丛生成簇状，叶片小而黄。
![]()
图 1 云南元谋小驳骨丛枝病症状Figure 1. The symptom of G. vulgaris Nees witches’-broom disease of Yuanmou County, Yunnan Province2.2 序列结果分析
2.2.1 16S rRNA、rp和secY基因序列扩增结果
小驳骨感病植株总DNA经PCR扩增，均得到大小与阳性对照一致的特异性条带，而双蒸水和健康植株总DNA均未扩增出特异性条带。测序结果表明：从小驳骨感病植株样品总DNA中分别扩增得到了1248 bp的16S rRNA基因片段(nested PCR) (MN535183)、1171 bp的rp基因片段(MN543078)及1425 bp的secY基因片段(MN543069)，这与预期扩增片段长度相符。由此将该植原体株系命名为云南元谋小驳骨丛枝植原体(Gendarussa vulgaris witches’-broom phytoplasma，GvWB-YNym)。
2.2.2 小驳骨丛枝植原体(GvWB-YNym) 16S rRNA、rp以及secY基因序列分析
(1) 16S rRNA基因序列分析
经BLAST搜索，显示GvWB-YNym的16S rRNA基因序列与16SrII-A亚组为最佳匹配同源序列。经MUSCLE进行在线比对结果(表3)所示：该序列片段与16SrII-A亚组的番茄巨芽植原体(Tomato big bud，TBB-YM，GenBank 登录号：JQ923436)的相似率高达100.00%。由图2所示：GvWB-YNym与植原体16SrII组明显聚为一枝，且又与归属为16SrII-A亚组的番茄巨芽植原体(Tomato big bud，TBB-YM，GenBank 登录号：JQ923436)及来自云南元谋的菊苣丛枝植原体(‘Cichorium intybus’ witches’-broom，CiWB-YNym，GenBank 登录号：KJ735774)处于同一枝。
表 3 小驳骨丛枝植原体(GvWB-YNym)16S rRNA、rp和secY基因核苷酸序列与部分16SrII亚组株系的同源性分析Table 3. Homology analysis of nucleotide sequences between GvWB-YNym 16S rRNA, rp and secY gene with other phytoplasma strains of 16SrII subgroup植原体菌株
phytoplasma strains亚组
subgroupGenBank登录号 GenBank accession No. 核苷酸同源性/% nucleotide homologous 16S rRNA secY rp 16S rRNA secY rp TBB-YM/ TBBP-YNym 16SrII-A JQ923436 KC953016 KC953010 100.00 99.93 99.83 PnWB 16SrII-A L33765 GU004331 EF193375 99.60 99.79 99.74 TWB-YNym 16SrII-A KC953005 KC953017 KC953011 99.68 99.86 100.00 BgWB 16SrII-A AB690304 AB703249 AB703243 99.68 100.00 99.91 CuPh 16SrII-A KY349111 KY365527 KY365523 99.44 99.84 99.82 CoP 16SrII-C EF186827 GU004350 EF186814 98.08 88.43 90.76 CrP 16SrII-C EF193355 GU004349 EF186818 98.00 88.25 90.51 SOYP 16SrII-C EF193353 GU004324 EF186816 98.00 87.44 90.60 PEP 16SrII-E Y16393 GU004348 EF193381 98.08 88.23 90.26 ![]()
图 2 基于16S rRNA、rp及secY基因序列构建的GvWB-YNym系统进化树Figure 2. Phylogenetic trees based on 16S rRNA, rp and secY gene sequences from GvWB-YNym and related phytoplasma strains(2) rp基因序列的分析
将测序所得的rp基因序列提交到NCBI网站ORF Finder进行在线分析，分析结果表明：rp序列包含了rps3 (nt550-1170)基因部分序列和rpl22 (nt90-476)基因全部序列，分别编码了206和128个氨基酸。由表3所示：该序列与16SrII-A亚组各株系相似率均高于99.74%，且与16SrII-A亚组中的番茄丛枝(TWB-YNym)相似率高达100%。由图2所示：GvWB-YNym与16SrII-A亚组同处1个分枝簇上，且与16SrII组聚为一大进化枝，这与基于16S rRNA基因序列构建的系统发育树具有相近之处，表明植原体类群之间具有相似的关联性。因此，依据rp基因序列可将GvWB-YNym归为16SrII-A亚组。
(3) secY基因序列的分析
将secY基因序列提交至NCBI网站ORF Finder进行在线分析，结果显示：secY基因序列包含了整个secY基因(nt74-1336)，编码了420个氨基酸。由表3所示：GvWB-YNym与植原体16SrII-A亚组中的5个植原体相似率达99.79%以上，与16SrII-A亚组中的印度黑豆丛枝(BgWB)相似率达到100%。由图2所示：secY基因与16SrII-A亚组的各个成员共同聚为一群，而与16SrII组其他亚组及其他外组亲缘关系相对较远；这也与基于rp基因和16S rRNA基因构建的系统发育树显示出相似的主次分枝顺序。据此结果推断，将GvWB-YNym划分为16SrII-A亚组。
2.2.3 小驳骨丛枝植原体(GvWB-YNym)16S rRNA基因在线分类鉴定结果
经iPhyClassifier在线分类鉴定结果显示：该株系16S rRNA基因序列与16SrII-A亚组代表性植原体花生丛枝植原体株系(L33765)相似性最高，相似系数高达1.00，且二者间的16S rDNA F2n/R2序列片段模拟RFLP图谱完全一致，上述所有分类结果均显示GvWB-YNym归属为16SrⅡ-A亚组。同时，在线分析还显示：GvWB-YNym与候选种‘Candidatus Phytoplasma aurantifolia’(U15442)同源性最高，达98.2%，因此可确定GvWB-YNym候选种与‘Ca.Phytoplasma aurantifolia’相关。
2.2.4 secY基因编码的蛋白亚基生物信息学分析
经分析工具ProtParam对GvWB-YNym的secY蛋白理化特性进行分析，结果显示其编码了420个氨基酸，分子量为47 720.21 u，理论等电点9.68，脂肪系数为136.45，不稳定系数31.94，属于稳定蛋白，总平均亲水性数值为0.664，表明该蛋白疏水性较好，属疏水性蛋白。通过在线分析网站Protscale对secY蛋白亲疏水性作进一步分析，结果如图3所示。从整体来看，该蛋白疏水氨基酸(正值)多于亲水氨基酸(负值)，表现为疏水性，与物理特性预测结果一致；其在215和216位氨基酸疏水性分值最高(3.661)，疏水性最强；在238和239位氨基酸亲水性分值最低(−2.678)，亲水性最强。
![]()
图 3 secY蛋白亲疏水性分析注：从左到右横向为氨基酸序列顺序；零水平线上方为疏水位点；零水平线下方为亲水位点。Figure 3. Hydrophilic analysis of secY proteinNotes: The horizontal direction is the amino acid sequence from left to right; the hydrophobic level is above the zero horizontal line; the hydrophilic position is below the zero horizontal line.以在线跨膜结构区预测软件TMHMM分析secY蛋白的跨膜域，结果(图4a)显示：该蛋白氨基酸序列具有10个显著且分布相对均匀的跨膜螺旋区域，表明该蛋白可能定位于细胞中与膜相关的结构，属于跨膜蛋白。信号肽是负责把蛋白引导到不同膜结构的亚细胞器内的短肽链，通过SignalP在线分析工具对secY基因进行信号肽分析，分析结果如图4b所示：GvWB-YNym的secY蛋白不存在信号肽。
![]()
图 4 secY蛋白跨膜结构及信号肽分析注：a) secY蛋白跨膜结构预测；b) secY蛋白信号肽预测，C-score. 原始剪切位点的分值，S-score. 信号肽的分值，Y-score. 综合剪切位点的分值。Figure 4. TMHMM posterior probabilities and prediction of the signal peptide of secY proteinNotes: a) membrane localization prediction; b) signal peptide prediction, C-score. original cut score, S-score. signal score, Y-score. general score.3. 讨论
由于植原体难以体外培养，多方面研究均落后于其他细菌，目前植原体的分类鉴定及其亲缘进化关系的研究领域仍然是热点。高度保守的16S rRNA基因序列常被用作植原体的系统发育进化关系及分类鉴定遗传标记[5]。随着计算机技术及网络系统的加速发展和日渐成熟完善，16S rRNA基因的虚拟RFLP图谱分析在植原体快速分类鉴定以及流行病学研究的分类学中越来越重要[16]。根据IRPCM (2004)规定[17]，16S rRNA基因序列与现有序列间的同源性大于97.5%时可被认为是相同组的植原体。然而，基于植原体中高度保守的16S rRNA基因序列的2.5%不同的任意阈值对植原体株系进行分组的依据，不能满足16S rRNA基因序列同源性虽大于97.5%但在植物寄主范围和介体昆虫等方面具有独特生态学和生物学的植原体株系进行分类的要求。因此，为了能在分子水平上更准确地区分植原体株系，便需要引入额外的分子标记。而作为辅助分类的其他遗传标记如核糖体蛋白(rp)基因、secY、tuf和23S rRNA基因以及16S-23S rRNA基因间隔区序列已被广泛应用。
比较分析给定的PnWB组(16SrII)各亚组之间的同源性发现：secY基因的序列相似性变化范围最大，其次为rp基因和16S rRNA基因；GvWB-YNym与16SrII-A亚组的亲缘关系最为接近，16S rRNA和rp基因序列在16SrII-A亚组各株系上表现出高相似性，但rp及secY基因序列比16S rRNA基因序列能够呈现出更大的遗传变异性，而secY基因相较于rp基因的序列又具有更大的可变程度，这些信息能够极大地增强对植原体组内遗传进化关系密切但不同亚组的分类鉴定能力，且与前人的研究[12, 18-22]具有一致性。通过基于16S rRNA、rp及secY基因构建的系统发育进化树分析，也能够支持这一论断。根据16S rRNA、rp及secY基因序列构建的系统发育进化树，都显示GvWB-YNym与16SrII-A亚组的各植原体株系形成单系群，且与其他16SrII组植原体聚为一大亚枝，但与16SrII组外的其他亚组植原体明显划分为不同进化枝。这3个系统进化树分析能够显著区分开植原体株系组与亚组之间的关系，进一步证实了以rp及secY基因作为分类辅助标记的可信度和重要性，是对利用16S rRNA基因对植原体进行分类的补充依据，也表明了GvWB-YNym归属到16SrII-A亚组的准确性。然而，分别基于16S rRNA、rp及secY基因构建的系统进化树中的各个分枝顺序并不完全一致，表明16S rRNA、rp及secY基因在其各自的进化史上遗传变异程度的不一致性。16S rDNA序列的RFLP分析代表不同的系统发育谱系，依据植原体在线分类网站iPhyClassifier对GvWB-YNym进行快速分类，结果表明GvWB-YNym与候选种‘Ca.Phytoplasma aurantifolia’相关，且与16SrII-A亚组代表性植原体株系花生丛枝植原体(L33765)的16S rDNA F2n/R2序列片段模拟RFLP图谱完全一致，进一步证实了将GvWB-YNym归属为16SrII-A亚组分类地位的准确性。这是小驳骨丛枝植原体的首次报道。
原核生物细胞的生长和发育依赖于Sec分泌蛋白转运机制，多数分泌蛋白和整合膜蛋白都需要经由此途径进行转运，其中起重要作用的是膜蛋白三聚体SecYEG和动力蛋白secA。在Sec转运酶中，secY、secE、secG和secA构成作为细胞质膜的输出机制中的转位酶复合物，secA、secE和secY是蛋白质易位和细胞活力必不可少的部分。在蛋白转运过程中，前体蛋白经由secA携带的伴胞分子的引导，同时利用ATP酶的水解活性促进该蛋白质穿过三聚体SecYEG形成的跨膜通道到达膜外[23]。植原体有2个分泌系统，YidC和Sec，后者几乎是绝大多数或所有植原体的共同分泌系统，KAKIZAWA等[24]发现：植原体中存在Sec系统，而secA和secY蛋白是其中重要组成部分，但其并不清楚Sec系统在植原体中的作用机制。岳红妮等[25]研究发现：泡桐丛枝(PaWB)植原体中同样存在的Sec分泌蛋白转运系统，并猜测该系统可能会直接将细菌蛋白质如毒素等转运到宿主细胞质或寄主昆虫细胞中，从而诱发宿主病症。本研究通过对secY蛋白进行生物信息参数的分析发现：该蛋白作为疏水性稳定跨膜蛋白存在于感病小驳骨植株中；secY蛋白含10个显著疏水跨膜区域，与大肠杆菌中的secY蛋白具有相似的跨膜结构区[26]，这种膜包埋的结构可赋予secY蛋白“转运蛋白”功能，为膜蛋白的分泌提供通道，但secY蛋白中不存在信号肽。由于植原体缺乏细胞壁，其细胞膜可直接与寄主植物的细胞内环境进行接触，一些与致病性相关的致病因子可能通过Sec转运系统直接释放到宿主细胞中，导致植株韧皮部功能受损，从而诱发寄主产生抗病性，影响植物的防御响应，进而直接或间接地导致寄主植物或介体昆虫发生互作，诱使寄主植物发生典型植原体病害症状。但secY蛋白在植原体内转运机制中扮演的具体角色尚不明确，还有待进一步研究和证实，本研究将为了解Sec转运系统在植原体的分子致病机制及防治植原体病害上提供更大的帮助。
4. 结论
本研究是首次对小驳骨丛枝植原体病害的相关报道，明确了小驳骨丛枝病是由植原体所侵染发生的，确定了该植原体株系16S rII-A亚组的分类地位；同时对secY蛋白在Sec分泌蛋白转运系统中的作用进行了初步探讨，为了解Sec转运系统在植原体的分子致病机理及防治此类病害提供帮助
-
![]()
图 1 云南元谋小驳骨丛枝病症状
Figure 1. The symptom of G. vulgaris Nees witches’-broom disease of Yuanmou County, Yunnan Province
![]()
图 2 基于16S rRNA、rp及secY基因序列构建的GvWB-YNym系统进化树
Figure 2. Phylogenetic trees based on 16S rRNA, rp and secY gene sequences from GvWB-YNym and related phytoplasma strains
![]()
图 3 secY蛋白亲疏水性分析
注：从左到右横向为氨基酸序列顺序；零水平线上方为疏水位点；零水平线下方为亲水位点。
Figure 3. Hydrophilic analysis of secY protein
Notes: The horizontal direction is the amino acid sequence from left to right; the hydrophobic level is above the zero horizontal line; the hydrophilic position is below the zero horizontal line.
![]()
图 4 secY蛋白跨膜结构及信号肽分析
注：a) secY蛋白跨膜结构预测；b) secY蛋白信号肽预测，C-score. 原始剪切位点的分值，S-score. 信号肽的分值，Y-score. 综合剪切位点的分值。
Figure 4. TMHMM posterior probabilities and prediction of the signal peptide of secY protein
Notes: a) membrane localization prediction; b) signal peptide prediction, C-score. original cut score, S-score. signal score, Y-score. general score.
表 1 引物序列表
Table 1 Primer sequence listing
基因 gene 引物名称 primer name 序列 sequence 预期扩增片段长度/bp expected amplified fragment length 16S rRNA P1 5′-AAGAGTTTGATCCTGGCTCAGGATT-3′ 1 800 P7 5′-CGTCCTTCATCGGCTCTT-3′ R16F2n 5′-GAAACGACTGCTAAGACT-3′ 1 248 R16R2 5′-TGACGGGCGGTGTGTACAAACCCCG-3′ rp rp(II)F1 5′-GCTCTTACTCGTAAAYATGTAGT-3′ 1 171 rp(II)R1 5′-TTACTTGATTTTCTGGTTTTGA-3′ secY secY-F 5′-GCGGAAGAAGCTATTAT-3′ 1 425 secY-R 5′-CGAATAACATAATATAATTGATTCC-3′ 
下载: 导出CSV
表 2 植原体各株系16S rRNA、rp和secY基因序列登录号
Table 2 The accession No. of 16S rRNA, rp gene and secY gene sequence of phytoplasma strain
植原体病害及相关菌株
phytoplasma diseases and strains地理来源
geographical origin16Sr 组-亚组
16Sr group-subgroupGenBank登录号 GenBank accession No. 16S rRNA secY rp 云南元谋番茄巨芽植原体
tomato big bud, TBB-YM/ TBBP-YNym中国云南
Yunnan, China16SrII-A JQ923436 KC953016 KC953010 花生丛枝植原体
peanut witches’-broom, PnWB中国台湾
Taiwan, China16SrII-A L33765 GU004331 EF193375 云南元谋番茄丛枝植原体
tomato witches’-broom, TWB-YNym中国云南
Yunnan, China16SrII-A KC953005 KC953017 KC953011 云南元谋四季豆花变绿植原体
string bean virescence, SbV-YNym2中国云南
Yunnan, China16SrII-A KJ735766 — — 芝麻花变叶植原体
Candidatus phytoplasma aurantifolia, SEPB2泰国
Thailand16SrII-A JN006079 — — 云南元谋兵豆花变绿植原体
‘Lens culinaris’ virescence, LeVP-YNym中国云南
Yunnan, China16SrII-A — KJ735763 KJ735764 云南元谋豇豆花变叶植原体
cowpea phyllody, CowP-YNym中国云南
Yunnan, China16SrII-A KC953001 — — 云南元谋豇豆花变绿植原体
cowpea virescence, CoVP-YNym中国云南
Yunnan, China16SrII-A — KC953013 KC953007 黑豆丛枝植原体
black gram witches’-broom, BgWB缅甸
Myanmar16SrII-A AB690304 AB703249 AB703243 黄瓜花变叶植原体
cucumber phyllody, CuPh伊朗
Iran16SrII-A KY349111 KY365527 KY365523 云南元谋菊苣从枝植原体
‘Cichorium intybus’ witches’-broom,
CiWB-YNym/ChWB-YNym中国云南
Yunnan, China16SrII-A KJ735774 KJ735775 — 柠檬从枝植原体
lime witches’-broom, LWB阿曼
Oman16SrII-B EF186828 — — 棉花花变叶植原体
cotton phyllody, CoP布基纳法索
Burkina Faso16SrII-C EF186827 GU004350 EF186814 猪屎豆花变叶植原体
Crotalaria phyllody, CrP泰国
Thailand16SrII-C EF193355 GU004349 EF186818 大豆花变叶植原体
soybean phyllody, SOYP泰国
Thailand16SrII-C EF193353 GU004324 EF186816 刺缘毛莲菜花变叶植原体
Picris echioides phyllody, PEP意大利
Italy16SrII-E Y16393 GU004348 EF193381 桃黄化植原体
peach X-disease，CX加拿大
Canada16SrIII-A L33733 GU004327 EF186813 三叶草黄化植原体
clover yellow edge, CYE立陶宛
Lithuania16SrIII-B AF173558 GU004332 JQ360962 核桃丛枝植原体
walnut witches’-broom, WWB美国乔治亚州
Georgia, USA16SrIII-G AF190227 GU004325 JQ360961 榆树黄化植原体
elm yellows, ‘Ca.P.ulmi’EYEu (EY-626、EY-627)意大利
Italy16SrV-A AY197657 AY197691 AY197678 悬钩子矮化植原体
rubus stunt, RUS意大利
Italy16SrV-E AY197648 AY197696 AY197668 马铃薯丛枝植原体
potato witches’-broom, ‘Ca.P.trifolli’PWB加拿大
Canada16SrVI-A AY500818 GU004316 AY197683 白蜡木黄化植原体
ash yellows, AshY1美国纽约
New York, USA16SrVII-A AF092209 GU004329 EF183492 欧洲核果黄化/李子坏死植原体
European stone fruit yellows/ plum, ‘Ca.P. prunorum’PLN德国/意大利
Germany/Italy16SrX-B AJ542545 GU004334 EF193369 绣球花花变叶植原体
hydrangea phyllody, HYDP比利时
Belgium16SrI-A AY265215 AY803181 AY264868 报春花花变绿植原体
primrose virescence, PRIVC德国
Germany16SrI-B AY265210 AY803176 AY264854 三叶草花变叶植原体
clover phyllody, CPh加拿大
Canada16SrI-C AF222065 AY803179 AY264862 蓝莓矮化植原体
blueberry stunt, BBS3美国
USA16SrI-E AY265213 GU004357 AY264863
下载: 导出CSV
表 3 小驳骨丛枝植原体(GvWB-YNym)16S rRNA、rp和secY基因核苷酸序列与部分16SrII亚组株系的同源性分析
Table 3 Homology analysis of nucleotide sequences between GvWB-YNym 16S rRNA, rp and secY gene with other phytoplasma strains of 16SrII subgroup
植原体菌株
phytoplasma strains亚组
subgroupGenBank登录号 GenBank accession No. 核苷酸同源性/% nucleotide homologous 16S rRNA secY rp 16S rRNA secY rp TBB-YM/ TBBP-YNym 16SrII-A JQ923436 KC953016 KC953010 100.00 99.93 99.83 PnWB 16SrII-A L33765 GU004331 EF193375 99.60 99.79 99.74 TWB-YNym 16SrII-A KC953005 KC953017 KC953011 99.68 99.86 100.00 BgWB 16SrII-A AB690304 AB703249 AB703243 99.68 100.00 99.91 CuPh 16SrII-A KY349111 KY365527 KY365523 99.44 99.84 99.82 CoP 16SrII-C EF186827 GU004350 EF186814 98.08 88.43 90.76 CrP 16SrII-C EF193355 GU004349 EF186818 98.00 88.25 90.51 SOYP 16SrII-C EF193353 GU004324 EF186816 98.00 87.44 90.60 PEP 16SrII-E Y16393 GU004348 EF193381 98.08 88.23 90.26
下载: 导出CSV
-
[1] MAEJIMA K, OSHIMA K, NAMBA S. Exploring the phytoplasmas, plant pathogenic bacteria[J]. Journal of General Plant Pathology, 2014, 80(3): 210. DOI: 10.1007/s10327-014-0512-8.
[2] LEE I M, DAVIS R E. Phytoplasma: phytopathogenic mollicutes[J]. Annual Review of Microbiology, 2000, 54: 221. DOI: 10.1146/annurev.micro.54.1.221.
[3] ASSUNTA B. Phytoplasmas: diversity, taxonomy, and epidemiology[J]. Frontiers in Bioscience, 2007, 12: 673. DOI: 10.2741/2092.
[4] LEE I M, HAMMOND R W, DAVIS R E, et al. Universal amplification and analysis of pathogen 16S rDNA for classification and identification of mycoplasmalike organisms[J]. Phytopathology, 1993, 83(8): 834. DOI: 10.1094/Phyto-83-834.
[5] LEE I M, GUNDERSEN-RINDAL D E, DAVIS R E, et al. Revised classification scheme of phytoplasmas based on RFLP analyses of 16S rRNA and ribosomal protein gene seqences[J]. International Journal of Systematic Bacteriology, 1998, 48(4): 1153. DOI: 10.1099/00207713-48-4-1153.
[6] ZHAO Y, WEI W, DAVIS R E, et al. Recent advances in 16S rRNA gene-based phytoplasma differentiation, classification and taxonomy[M]//WEINTRAUB P G, JONES P. Phytoplasmas: genomes, plant hosts and vectors. Oxfordshire: CAB International, 2010.
[7] PÉREZ-LÓPEZ E, WEI W, WANG J, et al. Novel phytoplasma strains of X-disease group unveil genetic markers that distinguish North American and South American geographic lineages within subgroups 16SrIII-J and 16SrIII-U[J]. Annals of Applied Biology, 2017, 171(3): 405. DOI: 10.1111/aab.12383.
[8] LEE I M, ZHAO Y, DAVIS R E, et al. Prospects of DNA-based systems for differentiation and classification of phytoplasmas[J]. Bulletin of Insectology, 2007, 60(2): 239. DOI: 10.1017/S000748530700524X.
[9] LEE I M. Classification of phytoplasma strains in the elm yellows group (16SrV) and proposal of ‘Candidatus Phytoplasma ulmi’ for the phytoplasma associated with elm yellows[J]. International Journal of Systematic and Evolutionary Microbiology, 2004, 54(2): 337. DOI: 10.1099/ijs.0.02697-0.
[10] LEE I M. ‘Candidatus Phytoplasma asteris’, a novel phytoplasma taxon associated with aster yellows and related diseases[J]. International Journal of Systematic and Evolutionary Microbiology, 2004, 54(4): 1037. DOI: 10.1099/ijs.0.02843-0.
[11] LEE I M, BOTTNER K D, SECOR G, et al. ‘Candidatus Phytoplasma americanum’, a phytoplasma associated with a potato purple top wilt disease complex[J]. International Journal of Systematic and Evolutionary Microbiology, 2006, 56(7): 1593. DOI: 10.1099/ijs.0.64251-0.
[12] MARTINI M, BOTTI S, MARCONE C, et al. Genetic variability among flavescence dorée phytoplasmas from different origins in Italy and France[J]. Molecular & Cellular Probes, 2002, 16(3): 197. DOI: 10.1006/mcpr.2002.0410.
[13] 国家药典委员会. 中华人民共和国药典: 一部[M]. 北京: 中国医药科技出版社, 2010. [14] 全国中草药汇编编写组. 全国中草药汇编: 上册[M]. 北京: 人民卫生出版社, 1975. [15] 万琼莲, 蔡红, 杨子祥, 等. 云南花生丛枝植原体secY基因序列分析及结构预测[J]. 云南农业大学学报(自然科学), 2013, 28(2): 163. DOI: 10.3969/j.issn.1004-390X(n).2013.02.003. [16] ZHAO Y, WEI W, LEE I M, et al. Construction of an interactive online phytoplasma classification tool, iPhyClassifier, and its application in analysis of the peach X-disease phytoplasma group (16SrIII)[J]. International Journal of Systematic and Evolutionary Microbiology, 2009, 59(10): 2582. DOI: 10.1099/ijs.0.010249-0.
[17] FIRRAO G, ANDERSEN M, BERTACCINI A, et al. ‘Candidcatus Phytoplasma’, a taxon for the wall-less, nonhelical prokaryotes that colonize plant phloem and insects[J]. International Journal of Systematic and Evolutionary Microbiology, 2004, 54(4): 1243. DOI: 10.1111/j.1747-5457.2005.tb00085.x.
[18] 王柱华, 刘俊男, 杨子祥, 等. 云南蔓草虫豆丛枝植原体16S rDNA和secY基因序列分析[J]. 云南农业大学学报(自然科学), 2018, 33(5): 826. DOI: 10.12101/j.issn.1004390x(n).201712004. [19] MARTINI M, LEE I M, BOTTNER K D, et al. Ribosomal protein gene-based phylogeny for finer differentiation and classifycation of phytoplasmas[J]. International Journal of Systematic and Evolutionary Microbiology, 2007, 57(9): 2037. DOI: 10.1099/ijs.0.65013-0.
[20] LEE I M, BOTTNER-PARKER K D, ZHAO Y, et al. Phylogenetic analysis and delineation of phytoplasmas based on secY gene sequences[J]. International Journal of Systematic and Evolutionary Microbiology, 2010, 60(12): 2887. DOI: 10.1099/ijs.0.019695-0.
[21] LEE I M, ZHAO Y, BOTTNER K D. SecY gene sequence analysis for finer differentiation of diverse strains in the aster yellows phytoplasma group[J]. Molecular & Cellular Probes, 2006, 20(2): 87. DOI: 10.1016/j.mcp.2005.10.001.
[22] 万琼莲, 杨子祥, 王连春, 等. 云南花生丛枝植原体16S rRNA和核糖体蛋白基因序列分析[J]. 植物病理学报, 2014, 44(4): 370. DOI: 10.13926/j.cnki.apps.2014.04.005. [23] BLAAUWEN T D, DRIESSEN A J M. Sec-dependent preprotein translocation in bacteria[J]. Archives of Microbiology, 1996, 165(1): 1. DOI: 10.1007/s002030050289.
[24] KAKIZAWA S, OSHIMA K, KUBOYAMA T, et al. Cloning and expression analysis of phytoplasma protein transloncation genes[J]. Molecular Plant Microbe Interactions, 2001, 14(9): 1043. DOI: 10.1094/MPMI.2001.14.9.1043.
[25] 岳红妮, 吴宽, 吴云锋, 等. 泡桐丛枝植原体Sec分泌蛋白转运系统3个亚基基因的克隆及蛋白特性分析[J]. 植物保护, 2009, 35(2): 25. DOI: 10.3969/j.issn.0529-1542.2009.02.006. [26] DERMAN A I, PUZISS J W, BASSFORD P J, et al. A signal sequence is not required for protein export in prlA mtants of Escherichia coli.[J]. The EMBO Journal, 1993, 12(3): 879. DOI: 10.1002/j.1460-2075.1993.tb05728.x.
计量
- 文章访问数: 3987
- PDF下载量: 34
