李新，顾晓川,2，龙海波，彭焕，黄文坤，彭德良



禾谷孢囊线虫果胶酸裂解酶新基因的鉴定与表达特征分析

李新1，顾晓川1,2，龙海波3，彭焕1，黄文坤1，彭德良1

（1中国农业科学院植物保护研究所植物病虫害生物学国家重点实验室，北京100193；2中国热带农业科学院橡胶研究所，海南儋州571737；3中国热带农业科学院环境与植物保护研究所农业部热带作物有害生物综合治理重点实验室，海口571101）

【】禾谷孢囊线虫（）是一种严重危害麦类作物的重要植物病原线虫，对农业生产造成巨大的经济损失，然而其致病机制和有效防控方法还有待进一步研究。通过克隆禾谷孢囊线虫的果胶酸裂解酶新基因，并对其表达特性进行分析，为后续探究的基因功能及其与寄主的互作提供理论依据，并为探讨禾谷孢囊线虫防控途径提供新思路。【】采用同源克隆结合RACE技术从禾谷孢囊线虫中克隆出一个新的果胶酸裂解酶基因；采用DNAMAN、Clustal、SignalP 4.0 Server和GSDS等相关生物信息学软件和在线工具分析该基因的核苷酸和氨基酸序列，并使用MEGA 5.0构建系统进化树；采用原位杂交和半定量PCR方法确定该基因的在禾谷孢囊线虫中的表达部位及其在线虫不同龄期中的表达情况。【】从禾谷孢囊线虫中成功克隆出一个果胶酸裂解酶基因（GenBank登录号GQ998895)，该基因cDNA全长1 717 bp，包含一个长度为1 563 bp的开放阅读框，编码一个长度为521个氨基酸残基的蛋白，其理论分子量为57.5 kD，理论等电点为8.52。从线虫基因组DNA中扩增获得长度为7 199 bp的基因组全长，基因结构显示分析发现，基因组包含14个外显子和13个内含子，除第3个内含子剪接位点是GC-AG外，其余12个内含子都符合真核生物基因剪接位点GT-AG规则。同源比对结果表明，预测蛋白Ha-PEL-1的C端与大豆孢囊线虫果胶酸裂解酶HG-PEL-1、甜菜孢囊线虫果胶酸裂解酶HS-PEL-1均有67%的一致性和83%的相似性；此外，其N端信号肽后比报道的其他植物寄生线虫果胶酸裂解酶多出一段长度为254个氨基酸残基的序列，这段序列中，靠近N端的184个氨基酸残基与数据库中的蛋白均无相似性，而靠近C端有70个氨基酸残基（Lys205—Glu274）与韦塞尔斯布朗病毒NS5蛋白（注册号3ELD）的甲基转移酶区域有32%的一致性和47%的相似性。氨基酸序列分析发现，预测蛋白Ha-PEL-1包含一个长度为20个氨基酸残基的信号肽和4个果胶酸裂解酶第3家族（PL3）的高度保守区域以及多个保守的半胱氨酸残基；系统进化分析发现，及其他已报道的线虫果胶酸裂解酶基因与细菌和真菌来源的PEL聚在一个大的分支中；原位杂交结果显示主要在禾谷孢囊线虫亚腹食道腺中表达；半定量RT-PCR确定在寄生前和寄生后的2龄幼虫中大量表达。【】通过对禾谷孢囊线虫中一个新果胶酸裂解酶基因的克隆和表达特征分析，揭示该基因与禾谷孢囊线虫的侵染和寄生过程密切相关。

禾谷孢囊线虫；果胶酸裂解酶；基因克隆；原位杂交；发育表达分析

0 引言

【研究意义】小麦禾谷孢囊线虫（）是一种严重危害小麦（）、大麦（）和燕麦（）等禾谷类作物的重要病原线虫[1]。自1874年在德国首次发现后，现已在全球40多个国家和地区发生与危害[2]。中国1989年在湖北省天门县岳口镇首次发现该病原线虫[3]，随着农业机械化和小麦跨区联合收割，小麦孢囊线虫病发生范围迅速扩大，危害程度日趋严重，目前已扩散至全国16个省（市、自治区）的小麦种植区[4]，对小麦等麦类作物的生产造成了严重的经济损失[5]。研究表明，由植物线虫食道腺细胞产生、并通过口针穿刺分泌到寄主植物组织中的纤维素酶、果胶酸裂解酶等一系列细胞壁降解酶类，在线虫的侵染过程中，能够降解和软化细胞壁，从而有利于线虫的寄生和迁移[6]。对重要的细胞壁降解酶——果胶酸裂解酶（pectate lyase，PEL）基因进行克隆和表达特性分析，可为后续探究的基因功能及其与寄主的互作机制提供理论依据，为从基因和蛋白水平上对禾谷孢囊线虫进行防治打下基础。【前人研究进展】2000年，Popeijus等首次从马铃薯金线虫（）中发现果胶酸裂解酶基因[7]，随后陆续从大豆孢囊线虫（）、甜菜孢囊线虫（）、马铃薯金线虫、根结线虫（spp.）、松材线虫（）和燕麦真滑刃线虫（）等多种植物线虫中分离获得该基因[8-16]。基因组数据分析预测显示，松材线虫、北方根结线虫（）、南方根结线虫（.）分别含有15、22和30个果胶酸裂解酶基因[17-19]。在大豆孢囊线虫研究中，先后发现了两类结构和功能具有显著差异的果胶酸裂解酶基因，Bakhetia等[20]和Peng等[10]先后采用RNAi技术沉默了大豆孢囊线虫的两个果胶酸裂解酶基因（和），结果发现，当被沉默后，大豆孢囊线虫在寄主体内的雌雄比发生变化，但是对侵染寄主根系的线虫数量没有影响；而沉默后，寄主根系内的线虫数量和雌成虫数均明显下降，说明两类结构和功能不同的果胶酸裂解酶在线虫侵染过程发挥着不同的作用。【本研究切入点】在前期研究中，禾谷孢囊线虫中已分离克隆出多个与植物细胞壁降解相关的-1,4-内切葡聚糖酶基因（-1,4-endoglucanase gene）[21-22]、扩展蛋白基因（expansion gene）[23-24]和纤维素结合蛋白基因（CBP）[25]，但关于果胶酸裂解酶的研究未见报道。【拟解决的关键问题】从禾谷孢囊线虫中克隆获得一个果胶酸裂解酶基因，对其序列特征、表达部位和发育动态进行系统研究，为进一步解析该类基因的功能打下基础。

1 材料与方法

试验于2015年9月至2016年12月在中国农业科学院植物保护研究所完成。

1.1 线虫材料

含有禾谷孢囊线虫成熟孢囊的土样采集于河南省濮阳市南乐县。采用蔗糖离心法分离获得孢嚢，并将挑取的孢囊在4℃条件下保存8周以上，然后放于16℃条件下中孵化2龄幼虫[26]，参考亓晓莉等[27]方法，挑取单条线虫进行分子检测，确定为禾谷孢囊线虫后，将其他2龄幼虫接种到感病寄主小麦温麦19的根系中，进行扩大培养，从根系和土壤中再次分离出来的孢囊，在4℃低温刺激后，16℃孵化2龄幼虫用于下一步研究。

1.2 试剂

Trizol Reagent、SuperScriptTMFirst-strand Synthesis System for RT-PCR反转录试剂盒、5′ RACE System for Rapid Amplification of cDNA Ends Version 2.0 kit、Dynbeads®mRNA DIRECTTMKit购自Invitrogen公司；3′-Full RACE Core Set Ver.2.0购自Takara公司；DIG High Prime DNA Labeling and Detection Starter Kit I、PCR DIG Probe Synthesis Kit、蛋白酶K购自Roche公司；RQ1 RNase-Free DNase、pGEM-T easy vector购自Promega公司；DNeasy Blood & Tissue Kit购自QIAGEN公司；DNA Marker、DH5感受态细胞、凝胶回收试剂盒购自北京天根生化有限责任公司；其他常规试剂购自Sigma等公司。

1.3 方法

1.3.1 核酸提取 收集新鲜孵化的禾谷孢囊线虫2龄幼虫置于1 mL的Trizol溶液中，参考彭焕等[28]方法提取总RNA，总RNA用RQ1 RNase-Free DNase进行去DNA处理，纯化后产物置于-80℃保存。第一链cDNA的合成使用Invitrogen公司的SuperScriptTMFirst-Strand Synthesis System for RT-PCR试剂盒，反转录产物置于-20℃保存。参照DNeasy Blood & Tissue Kit试剂盒说明书提取禾谷孢囊线虫基因组DNA。

1.3.2 RACE 扩增3′端和5′端序列 根据已报道的植物寄生线虫果胶酸裂解酶蛋白序列，参照Takara公司3′-Full RACE Core Set Ver.2.0试剂盒引物设计说明，采用Primier 5.0软件，设计3′ RACE反应的简并引物（JPL-O和JPL-I，表1），以2龄幼虫总RNA为模板，进行3′末端扩增，具体步骤参照试剂盒说明书，然后对扩增片段进行回收、连接、转化、送生工生物工程（上海）股份有限公司进行测序。再根据已测序获得的基因3′末端片段序列，参考Invitrogen公司5′ RACE System for Rapid Amplification of cDNA Ends，Version 2.0 kit试剂盒引物设计说明，设计的5′ RACE特异引物（GSP-PL1、GSP-PL2和GSP-nPL4，表1），以2龄幼虫总RNA为模板，进行5′末端扩增，具体步骤参照试剂盒说明书，对扩增片段进行回收、连接、转化、测序。

1.3.3 cDNA全长和基因组DNA克隆 根据已测序获得的基因cDNA拼接全长序列，设计全长引物（PEL-L1和PL-LA2，表1），以第一链cDNA和基因组DNA为模板，分别进行PCR扩增。反应参数：94℃预变性5 min；94℃变性30 s，50℃退火30 s，72℃延伸2 min；35个循环；72℃延伸10 min；对扩增片段进行回收、连接、转化、测序。

1.3.4 序列分析 分别使用DNAMAN和Clustal软件进行核苷酸序列的翻译和氨基酸序列的比对分析，使用在线工具（http://www.expasy.ch/tools/ pi_tool.html/）进行蛋白质等电点和分子量的预测，使用SignalP 4.0 Server（http://www.cbs.dtu.dk/services/ SignalP/）预测蛋白质前体信号肽[29]，使用CBS Prediction Servers（http://www.cbs.dtu.dk/services/）对氨基酸序列中可能存在的糖基化位点和磷酸化位点进行预测，使用基因内含子分析工具（http://gsds.cbi.pku. edu.cn/）分析基因组结构[30]，使用MEGA 5.0构建系统进化树[31]。

1.3.5 原位杂交 原位杂交方法参考De Boer等[32]的方法进行，收集新孵化的禾谷孢囊线虫2龄幼虫用DEPC水清洗3次，4%多聚甲醛固定24 h，用手术刀将虫体切断后，采用甲醇和丙酮进行通透，蛋白酶K进行消化，再用1倍杂交液重悬并置于-20℃保存。根据全长序列设计特异性引物Yw-F和Yw-R（表1），以cDNA第一链为模板加入DIG标记的dNTP进行正、负链探针合成。将合成的DIG标记探针与预处理的禾谷孢囊线虫参考DIG High Prime DNA Labeling and Detection Starter Kit Ⅰ说明书进行杂交和洗脱，停止显色后，在Leica DM2500显微镜下进行观察拍照。

1.3.6 发育表达分析 参考Long等[22]方法，通过控制接种时间结合酶裂解处理分离获得5种不同发育阶段的禾谷孢囊线虫（侵染前2龄幼虫、侵染后2龄幼虫、3龄幼虫、4龄幼虫和雌虫）。采用Dynbeads®mRNA DIRECTTMKit提取不同龄期禾谷孢囊线虫mRNA，并取等量mRNA（200 ng）分别反转录合成第一链cDNA。分别以5个不同发育阶段线虫的第一链cDNA和小麦根系的第一链cDNA为模板，采用特异引物RT-F和RT-R（表1）进行PCR 扩增，30个循环，同时以引物ActinF和ActinR（表1）扩增禾谷孢囊线虫的作为阳性对照。扩增产物在1.5%琼脂糖凝胶中电泳分离检测。

表1 本研究中的引物序列

2 结果

2.1 核酸提取及Ha-pel-1克隆

采用Trizol法提取的禾谷孢囊线虫2龄幼虫总RNA，电泳分离检测可见总RNA的28S、18S和5S 3条清晰的带，其OD260/280值为1.9—2.0，说明蛋白质和DNA去除较彻底，RNA纯度较高，可以用于后续试验。电泳检测提取的禾谷孢囊线虫基因组DNA，获得大小为23 kb左右的完整条带，无拖尾和降解，完整性较好，可用于下游试验。

通过同源克隆策略，采用3′ RACE技术扩增获得长度为280 bp的禾谷孢囊线虫果胶酸裂解酶基因cDNA 3′末端片段（图1-A），BLAST比对分析表明，该片段与已报道的线虫果胶酸裂解酶基因具有较高的同源性。在此基础上，采用5′ RACE技术获得长度为1 499 bp的cDNA 5′末端片段（图1-B）。为验证拼接序列的准确性，采用全长引物从禾谷孢囊线虫cDNA模板中扩增获得到长度1 682 bp，包含完整ORF的cDNA全长序列（图1-C）；采用分步克隆和序列拼接，从禾谷孢囊线虫基因组DNA中扩增获长度为7 199 bp的基因组全长。

1：PCR产物Product of PCR；M1：DL 2000 plus marker；M2：DL 2000 marker

2.2 Ha-pel-1序列分析

序列分析结果表明，（GenBank登录号GQ998895）cDNA全长为1 717 bp（不包含多聚腺苷酸尾）（图2），其中5′端非编码区（UTR）为98 bp，3′端非编码区为56 bp，开放阅读框为1 563 bp，编码一个长度为521个氨基酸残基的蛋白，其理论分子量为57.5 kD，理论等电点为8.52。SignalP预测发现，Ha-PEL-1蛋白N端含有一个长度为20个氨基酸残基的信号肽，剪接位点在Gly20和Arg21之间。糖基化位点预测显示，该蛋白在Asn38存在一个N位糖基化位点。基因结构显示分析发现，基因组含有14个外显子和13个内含子，除第3个内含子的剪接位点为GC-AG外，其余12个内含子都符合真核生物基因剪接位点GT-AG规则，其中相位0的内含子有5个，相位1的内含子有6个，相位2的内含子有2个。

图2 Ha-pel-1基因组DNA结构分析

2.3 Ha-PEL-1氨基酸序列分析及系统进化分析

NCBI中BLASTp比对分析发现，Ha-PEL-1预测蛋白C端的247个氨基酸残基与已报道的大豆孢囊线虫果胶酸裂解酶HG-PEL-1、甜菜孢囊线虫HS-PEL-1均有67%的一致性和83%的相似性，与大豆孢囊线虫HG-PEL-2和马铃薯孢囊线虫GR-PEL-1的相似性分别为82%和81%；其N端信号肽后的254个氨基酸残基与GenBank蛋白数据库中的蛋白无相似性，而与RCSB Protein Data Bank比对后发现，这段氨基酸序列后70个氨基酸残基（Lys205—Glu274）与RCSB蛋白质数据库中韦塞尔斯布朗病毒NS5蛋白（注册号3ELD）甲基转移酶区域的71个氨基酸残基具有32%的一致性和47%的相似性。

多序列比对和保守结构域分析表明，预测蛋白Ha-PEL-1具有果胶酸裂解酶第3家族（PL3）的4个保守区域和多个保守的半胱氨酸残基，同时，在Ha-PEL-1上发现一个带负电荷的保守天冬氨酸和两个带正电荷的保守赖氨酸，可以确定Ha-PEL-1是果胶酸裂解酶第3家族的新成员。

采用ME法构建了Ha-PEL-1及42个已报道的植物线虫、真菌和细菌来源的果胶酸裂解酶的系统发育树。结果显示，所有植物线虫来源的果胶酸裂解酶聚为3个独立的分支，其中燕麦滑刃线虫AA-PEL-1和AA-PEL-2、松材线虫BX-PEL-1和BX-PEL-2、拟松材线虫（）BM-PEL-1和BM-PEL-2、马铃薯孢囊线虫HG-PEL-2、甜菜孢囊线虫HS-PEL-2和大豆孢囊线虫HG-PEL-5聚为一个分支，命名为Nematode PEL-1；南方根结线虫MI-PEL-1和爪哇根结线虫（）MJ-PEL-1聚在另一个分支，命名为Nematode PEL-2；而禾谷孢囊线虫Ha-PEL-1与大豆孢囊线虫HG-PEL-1和HG-PEL-2、甜菜孢囊线虫HS-PEL-1、马铃薯孢囊线虫GR-PEL-1、南方根结线虫MI-PEL-2和MI-PEL-3聚为一个分支，命名为Nematode PEL-3。其中第1分支Nematode PEL-1和第2分支Nematode PEL-2又与芽孢杆菌（sp.）、番茄细菌性溃疡病菌（）等的果胶酸裂解酶聚在一个大分支上，第3分支Nematode PEL-3与菊欧式杆菌（）、纤维堆囊菌（）等细菌聚为一支（图3）。

2.4 Ha-pel-1组织定位

采用原位杂交方法明确了在禾谷孢囊线虫中的表达部位。地高辛标记的负链探针杂交结果表明，在禾谷孢囊线虫的两个亚腹食道腺细胞中有明显杂交信号（图4-A），而作为对照的正链探针无杂交信号（图4-B），由此确定在禾谷孢囊线虫的亚腹食道腺细胞中特异表达。

2.5 Ha-pel-1发育表达分析

采用半定量RT-PCR方法对在禾谷孢囊线虫5个不同发育阶段中的表达差异进行了分析。结果表明，主要在禾谷孢囊线虫寄生前2龄幼虫和寄生后2龄幼虫中表达，而在寄生后期和雌成虫中表达量很低甚至检测不到，内标基因在各个龄期中均有表达，在小麦根系DNA的阴性对照中未检测到和的表达（图5）。

浅色部分：线虫PEL；深色部分：真菌PEL；无底色部分：细菌PEL；下划线：Ha-PEL-1

3 讨论

植物细胞壁是寄主植物抵御外界病原物侵染的物理屏障，其主要成分为纤维素、果胶质和半纤维素等，降解和软化植物细胞壁是植物线虫成功侵染和寄生的关键。目前的研究结果表明，植物线虫能从食道腺细胞中分泌一系列的细胞壁修饰酶类，共同作用于植物细胞壁，从而利于线虫侵染和寄生[6,33-34]，其中果胶酸裂解酶在植物细胞壁降解过程中发挥着不可替代的作用。本研究从禾谷孢囊线虫中分离出一个果胶酸裂解酶基因，并对其序列结构、系统进化关系及表达特性进行了系统研究。

结构域分析发现，Ha-PEL-1氨基酸序列C端有一个长度为247个氨基酸残基的果胶酸裂解酶区域，比报道的其他植物寄生线虫果胶酸裂解酶氨基酸序列多出一段长度为250个氨基酸左右的序列。这段序列中，靠近N端的184个氨基酸残基与数据库中的蛋白均无相似性，而靠近C端有70个氨基酸残基与韦塞尔斯布朗病毒NS5蛋白（注册号3ELD）的甲基转移酶区域有32%的一致性和47%的相似性。韦塞尔斯布朗病毒是一种危害人类的病毒，其NS5蛋白参与了病毒复制，该蛋白的甲基转移酶区域能够参与RNA的加帽反应，催化RNA的甲基化[35]。到目前为止，在其他病原线虫的果胶酸裂解酶中均未发现该结构域。在Ha-PEL-1中，果胶酸裂解酶结构域与类似甲基转移酶结构域结合在一起，该段与韦塞尔斯布朗病毒NS5蛋白具有相似性的氨基酸序列是否也具有甲基转移酶功能，有待于今后进一步研究。

A：Ha-pel-1负链探针杂交In situ hybridization of digoxigenin-labelled antisence probes of Ha-pel-1；B：Ha-pel-1正链探针杂交In situ hybridization of digoxigenin-labelled sence probes of Ha-pel-1；S：口针Stylet；M：中食道球Metacorpus；G：亚腹食道腺Subventral oesophageal gland region

1：侵染前2龄幼虫Pre-parasitic 2nd stage juveniles；2：侵染后2龄幼虫Parasitic 2nd stage juveniles；3：3龄幼虫Parasitic 3rd stage juveniles；4：4龄幼虫Parasitic 4th stage juveniles；5：雌虫Adult females；6：阴性对照Negative control

序列同源比对结果显示，禾谷孢囊线虫Ha-PEL-1与大豆孢囊线虫HG-PEL-1、HG-PEL-2，甜菜孢囊线虫HS-PEL-1和马铃薯孢囊线虫GR-PEL-1的相似度高达81%—83%；进化分析也发现上述果胶酸裂解酶序列被聚类到同一个分支中，连同南方根结线虫MI-PEL-2和MI-PEL-3，聚为线虫果胶酸裂解酶的第3分支Nematode PEL-3；两个根结线虫果胶酸裂解酶聚在一起单独形成线虫果胶酸裂解酶的第2分支Nematode PEL-2；其他的线虫果胶酸裂解酶则聚为第1分支Nematode PEL-1，推测这些植物线虫PEL可能具有3个不同的进化起源。此外，这3个线虫果胶酸裂解酶进化分支均与细菌或真菌来源的果胶酸裂解酶聚在一起，在植物寄生线虫的其他细胞壁降解酶中也曾发现类似现象，如纤维素酶和扩展蛋白等[36]，由上述结果推测植物线虫PEL的起源可能与其他微生物的基因水平转移有关，而植物寄生线虫分泌细胞壁降解酶也正是线虫与真菌和细菌之间存在基因水平转移的重要证明[36]。综上所述，本研究结果为植物寄生线虫果胶酸裂解酶的起源过程提供了新的依据。

原位杂交结果发现，在禾谷孢囊线虫的亚腹食道腺细胞中特异表达，这与其他的植物寄生线虫的PEL，如HG-PEL-1和HG-PEL-2、HS-PEL-1和HS-PEL-2、GR-PEL-1等的表达部位一致[7-16]，此外，Ha-PEL-1蛋白N端具有信号肽，推测Ha-PEL-1是由禾谷孢囊线虫亚腹食道腺细胞合成，并通过口针分泌到寄主体内；同时，禾谷孢囊线虫亚腹食道腺细胞也是其他细胞壁降解酶类的合成场所[21-24]。

突破植物细胞壁的物理屏障是植物线虫成功寄生的关键，在禾谷孢囊线虫侵染过程中，食道腺细胞会分泌一系列细胞壁降解酶类如果胶酸裂解酶、-1,4-内切葡聚糖酶[21-22]和扩展蛋白[23-24]等分泌物，协同口针穿刺的机械压力来完成对寄主的侵染和早期寄生。禾谷孢囊线虫在5个不同龄期的表达分析显示，该基因在寄生前2龄幼虫和寄生后2龄幼虫中大量表达，在寄生后期和雌成虫中表达量非常低或者不表达；Long等研究发现禾谷孢囊线虫-1,4-内切葡聚糖酶Ha-eng-1a、Ha-eng-2和Ha-eng-3及扩展蛋白Ha-EXPB1、HaEXPB2也均在侵染和寄生早期大量表达，由此推测可能在禾谷孢囊线虫侵染植物寄主的过程中与上述两类蛋白共同发挥作用[21-24]。

本研究从禾谷孢囊线虫中克隆获得一个果胶酸裂解酶基因，并对其序列特征、组织定位和发育表达进行了系统分析，结果可为明确在禾谷孢囊线虫与寄主互作过程中的功能研究打下基础，并为探究植物寄生线虫果胶酸裂解酶的起源提供了新的证据，后续将利用RNAi等技术，进一步对进行功能验证；此外，Ha-PEL-1中与韦塞尔斯布朗病毒NS5蛋白甲基转移酶区域相似的氨基酸序列是否也具有甲基转移酶功能，仍有待进一步研究。

4 结论

从禾谷孢囊线虫中克隆获得一个新的果胶酸裂解酶基因，该基因cDNA全长1 717 bp，编码一个长度为521个氨基酸残基的蛋白Ha-PEL-1；基因组全长为7 199 bp，包含14个外显子和13个内含子。对其序列特征和表达特性进行系统研究表明，在Ha-PEL-1中，类似甲基转移酶结构域与果胶酸裂解酶结构域结合在一起；Ha-PEL-1属于果胶酸裂解酶第3家族，根据进化分析推测其起源可能与其他微生物的基因水平转移有关；主要在禾谷孢囊线虫寄生前和寄生后的2龄幼虫的亚腹食道腺中大量表达，揭示该基因与禾谷孢囊线虫的侵染和寄生过程密切相关。

[1] Meagher J W. World dissemination of the cereal-cyst nematode () and its potential as a pathogen of wheat., 1977, 9(1): 9-15.

[2] Rivoal R, Nicol J M. Past research on the cereal cyst nematode complex and future needs. Cereal cyst nematodes: status, research and outlook//Antalya, Turkey, 2009: 3-10.

[3] 陈品三, 王明祖, 彭德良. 我国小麦禾谷孢囊线虫(Wollenweber)的发现与鉴定初报. 中国农业科学, 1991, 24(5): 89.

Chen P S, Wang M Z, Peng D L. Preliminary report of identification on cereal cyst nematode of wheat in china., 1991, 24(5): 89. (in chinese)

[4] Cui J K, Huang W K, Peng H, Liu S M, Wang G F, Kong L A, Peng D L. A new pathotype characterization of Daxing and Huangyuan populations of cereal cyst nematode () in China., 2015, 14(4): 724-731.

[5] Peng D L, Nicol J M, Li H M, Hou S Y, Li H X, Chen S L, Ma P, Li H L, Riley I T.Current knowledge of cereal cyst nematode () on wheat in China. Cereal cyst nematodes: status, research and outlook//. Antalya, Turkey, 2009: 29-34.

[6] Bohlmann H, Sobczak M. The plant cell wall in the feeding sites of cyst nematodes., 2014, 5(1): 89.

[7] Popeijus H, Overmars H, Jones J, Blok V, Goverse A, Helder J, Schots A, Bakker J, Smant G. Degradation of plant cell walls by a nematode., 2000, 406(6791): 36-37.

[8] De Boer J M, McDermott J P, Davis E L, Hussey R S, Popeijus H, Smant G, Baum T J. Cloning of a putative pectate lyase gene expressed in the subventral esophageal glands of., 2002, 34(1): 9-11.

[9] 彭焕, 彭德良, 黄文坤, 贺文婷, 胡先奇. 大豆孢囊线虫果胶酸裂解酶基因的克隆与分析. 中国农业科学, 2012, 45(5): 854-866.

Peng H, Peng D L, Huang W K, He W T, Hu X Q. Molecular cloning and analysis of a novel pectate lyase genefrom soybean cyst nematode., 2012, 45(5): 854-866. (in chinese)

[10] Peng H, Cui J K, Long H B, Huang W K, Kong L A, Liu S M, He W T, Hu X Q, Peng D L. Novel pectate lyase genes ofplay key roles in the early stage of parasitism., 2016, 11(3): e0149959.

[11] Vanholme B, Van Thuyne W, Vanhouteghem K, De Meutter J, Cannoot B, Gheysen G. Molecular characterization and functional importance of pectate lyase secreted by the cyst nematode., 2007, 8(3): 267-278.

[12] Kudla U, Milac A L, Qin L, Overmars H, Roze E, Holterman M, Petrescu A J, Goverse A, Bakker J, Helder J, Smant G. Structural and functional characterization of a novel, host penetration-related pectate lyase from the potato cyst nematode., 2007, 8(3): 293-305.

[13] Doyle E A, Lambert K N. Cloning and characterization of an esophageal-gland-specific pectate lyase from the root-knot nematode., 2002, 15(6): 549-556.

[14] Huang G, Dong R, Allen R, Davis E L, Baum T J, Hussey R S. Developmental expression and molecular analysis of twopectate lyase genes., 2005, 35: 685-692.

[15] Kikuchi T, Shibuya H, Aikawa T, Jones J T. Cloning and characterization of pectate lyases expressed in the esophageal gland of the pine wood nematode., 2006, 19(3): 280-287.

[16] Karim N, Jones J T, Okada H, Kikuchi T. Analysis of expressed sequence tags and identification of genes encoding cell- wall-degrading enzymes from the fungivorous nematode., 2009, 10: 525.

[17] Kikuchi T, Cotton J A, Dalzell J J, Hasegawa K, Kanzaki N, McVeigh P, Takanashi T, Tsai I J, Assefa S A, Cock P J A, Otto T D, Hunt M, Reid A J, Sanchez-Flores A, Tsuchihara K, Yokoi T, Larsson M C, Miwa J, Maule A G, Sahashi N, Jones J T, Berriman M. Genomic insights into the origin of parasitism in the emerging plant pathogen., 2011, 7(9): e1002219.

[18] Opperman C H, Bird D M, Williamsond V M, Rokhsar D S, Burke M, Cohn J, Cromera J, Diener S, Gajan J, Graham S, Houfek T D, Liu Q L, Mitros T , Schaff J, Schaffer R, Scholl E, Sosinski B R, Thomas V P, Windham E. Sequence and genetic map of: A compact nematode genome for plant parasitism., 2008, 105(39): 14802-14807.

[19] Abad P, Gouzy J, Aury J M, Castagnone-Sereno P, Danchin E G J, Deleury E, Perfus-Barbeoch L, Anthouard V, Artiguenave F, Blok V C, Caillaud M C, Coutinho P M, Dasilva C, De Luca F, Deau F, Esquibet M, Flutre T, Goldstone J V, Hamamouch N, Hewezi T, Jaillon O, Jubin C, Leonetti P, Magliano M, Maier T R, Markov G V, McVeigh P, Pesole G, Poulain J, Robinson-Rechavi M, Sallet E, Ségurens B, Steinbach D, Tytgat T, Ugarte E, Van Ghelder C, Veronico P, Baum T J, Blaxter M, Bleve-Zacheo T, Davis E L, Ewbank J J, Favery B, Grenier E, Henrissat B, Jones J T, Laudet V, Maule A G, Quesneville H, Rosso M N, Schiex T, Smant G, Weissenbach J, Wincker P. Genome sequence of the metazoan plant-parasitic nematode., 2008, 26(8): 909-915.

[20] Bakhetia M, Urwin P E, Atkinson H J. qPCR analysis and RNAi define pharyngeal gland cell-expressed genes ofrequired for initial interactions with the host., 2007, 20(3): 306-312.

[21] Long H B, Peng H, Huang W K, Wang G F, Gao B L, Moens M, Peng D L. Identification and molecular characterization of a new-1,4-endoglucanase gene () in the cereal cyst nematode., 2012, 134(2): 391-400.

[22] Long H B, Peng D L, Huang W K, Peng H, Wang G F. Molecular characterization and functional analysis of two new-1,4-endoglucanase genes () from the cereal cyst nematode., 2013, 62(4): 953-960.

[23] Long H B, Peng D L, Huang W K, Liu Y K, Peng H. Identification of a putative expansin gene expressed in the subventral glands of the cereal cyst nematode, 2012, 14(5): 571-577.

[24] Liu J, Peng H, Cui J K, Huang W K, Kong L A, Clarke L J H, Jian H, Wang G L, Peng D L. Molecular characterization of a novel effector expansin-like protein fromthat induces cell death in., 2016, 6: 35677.

[25] 顾晓川, 彭德良, 彭焕, 龙海波, 王高峰, 黄文坤, 何月秋. 禾谷孢囊线虫()纤维素结合蛋白基因()的克隆和序列分析. 植物病理学报, 2011, 41(3): 240-246.

Gu X C, Peng D L, Peng H, Long H B, Wang G F, Huang W K, He Y Q. Molecular cloning and sequencing of cellulose binding protein gene () from the cereal cyst nematode ()., 2011, 41(3): 240-246. (in chinese)

[26] Liu Y K, Huang W K, Long H B, Peng H, He W T, Peng D L. Molecular characterization and functional analysis of a new acid phosphatase gene () from., 2014, 13(6): 1303-1310.

[27] 亓晓莉, 彭德良, 彭焕, 龙海波, 黄文坤, 贺文婷. 基于SCAR标记的小麦禾谷孢囊线虫快速分子检测技术. 中国农业科学, 2012, 45(21): 4388-4395.

Qi X L, Peng D L, Peng H, Long H B, Huang W K, He W T. Rapid molecular diagnosis based on SCAR marker system for cereal cyst nematode., 2012, 45(21): 4388-4395. (in chinese)

[28] 彭焕, 彭德良, 黄文坤. 甘薯茎线虫-1,4内切葡聚糖酶基因() cDNA全长的克隆与序列分析. 农业生物技术学报, 2009, 17(6): 1035-1041.

Peng H, Peng D L, Huang W K. Molecular cloning and sequence analysis of a new-1, 4-endoglucanase gene () from plant-parasitic nematodeon sweetpotato in China., 2009, 17(6): 1035-1041. (in Chinese)

[29] Petersen T N, Brunak S, Von H G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions., 2011, 8(10): 785-786.

[30] 郭安源, 朱其慧, 陈新, 罗静初. GSDS：基因结构显示系统. 遗传, 2007, 29(8): 1023-1026.

Guo A Y, Zhu Q H, Chen X, Luo J C. GSDS: a gene structure display server., 2007, 29(8): 1023-1026. (in chinese)

[31] Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods., 2011, 28(10): 2731-2739.

[32] De Boer J M, Yan Y, Smant G, Davis E L, Baum T J.hybridization to messenger RNA of., 1998, 30(3): 309-312.

[33] Vanholme B, De Meutter J, Tytgat T, Van Montagu M, Coomans A, Gheysen G. Secretions of plant-parasitic nematodes: a molecular update., 2004, 332: 13-27.

[34] 彭德良, 郑经武, 廖金铃, 万方浩. 重要植物线虫致病相关基因研究进展//彭友良. 中国植物病理学会2006年学术年会论文集. 北京: 中国农业科学技术出版社, 2006: 221-229.

Peng D L, Zheng J W, Liao J L, Wan F H. Advanced on putative patasitim genes of plant-parasitic nematodes and prospective//Peng Y L.. Beijing: China Agricultural Science and Technology Press, 2006: 221-229. (in Chinese)

[35] Bollati M, Milani M, Mastrangelo E, Ricagno S, Tedeschi G, Nonnis S, Decroly E, Selisko B, Lamballerie X D, Coutard B, Canard B, Bolognesi M. Recognition of RNA cap in the Wesselsbron virus NS5 methyltransferase domain: implications for RNA-capping mechanisms in Flavivirus., 2009, 385(1): 140-152.

[36] Haegeman A, Jones J T, Danchin E G J. Horizontal gene transfer in nematodes: a catalyst for plant parasitism?, 2011, 24(8): 879-887.

（责任编辑 岳梅）

Identification and expression analysis of a new pectate lyase genefrom

LI Xin1, GU XiaoChuan1,2, LONG HaiBo3, PENG Huan1, HUANG WenKun1, PENG DeLiang1

(1Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193;2Rubber Research Institute, Chinese Academy of Tropical Agricultural Sciences, Danzhou 571737, Hainan;3Key Laboratory of Pests Comprehensive Governance for Tropical Crops, Ministry of Agriculture, Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101)

【】The cereal cyst nematode () is one of the important plant parasitic nematodes which seriously threatened cereal crops and caused huge economic losses to agricultural production in China. However, its pathogenic mechanism and effective prevention and control methods still need to be further studied. The objective of this study is to provide a theoretical basis for further study on the gene function ofand its interaction with host plants, and to give new ideas for the control strategies of cereal cyst nematode based on the cloning and expression analysis of a new pectate lyase genefrom.【】A novel pectate lyase genewas cloned fromusing homology cloning combined with RACE technology, and its nucleotide sequence and amino acid sequence were analyzed by related bioinformatics softwares and online tools, such as DNAMAN, Clustal, SignalP 4.0 Server and GSDS. a phylogenetic tree was also constructed using MEGA 5.0. the tissue localization and developmental expression characteristics ofwere analyzed byhybridization and semi-quantitative PCR method. 【】A novel pectate lyase gene (, GenBank accession number GQ998895) was cloned successfully from.was 1 717 bp in length which contained a 1 563 bp open reading frame (ORF) encoding a protein of 521 amino acid residues. The molecular weight ofencoding protein was 57.5 kD and isoelectric point was 8.52. The full length of genomic sequence ofwas amplified from the nematode genome DNA which contains 7 199 bp. Gene structure analysis showed that thegenome contains 14 exons and 13 introns, except for the 3rd intron splice sites are GC-AG, the other 12 introns are in line with the rules of the eukaryotic gene splicing site GT-AG. The results of homologous comparison showed that the C-terminal sequence of the putative Ha-PEL-1 had a 67% identity and a similarity of 83% with that of soybean cyst nematode HG-PEL-1 and beet cyst nematode HS-PEL-1. In addition, after the end of N-terminal signal peptide, the putative Ha-PEL-1 had a sequence of 254 amino acid residues more than other reported plant parasitic nematodes pectate lyases. In this sequence, 184 amino acid residues closing to the N-terminal had no similarity with protein database, while 70 amino acid residues (Lys205-Glu274) closing to the C-terminal had an identity of 32% and a similarity of 47% with the methyltransferase domain of Wesselsbron virus NS5 (Registration No. 3ELD). The amino acid sequence analysis revealed that the predicted protein contained a signal peptide of 20 amino acid residues, as well as 4 highly conserved regions and several conserved cysteine residues characteristic of class Ⅲ pectate lyases (PL3). A phylogenetic analysis revealed thatand other nematodes pectate lyase genes are gathered in a large branch with bacterial and fungal sources PEL.hybridization analyses showed that the transcripts ofwere mainly expressed in the two subventral gland cells of. a semi-quantitative RT-PCR analysis confirmed that its transcriptions were highly expressed at the pre-parasitic and parasitic 2nd stage juveniles.【】A new pectate lyase genefrom, closely related to the infection and parasitic process of cereal cyst nematode, was found and analyzed.

; pectate lyase; gene cloning;hybridization; developmental expression analysis

10.3864/j.issn.0578-1752.2017.19.009

2017-03-30；接受日期：2017-05-08

国家“973”计划（2013CB127502）、国家公益性行业（农业）科研专项（201503114）

李新，E-mail：lixin9745@163.com。通信作者彭德良，Tel：010-62815611；E-mail：dlpeng@ippcaas.cn