烟草Hsc70-2 的克隆及对马铃薯Y 病毒侵染烟草的促进作用

奥鹏网院作业

烟草Hsc70-2 的克隆及对马铃薯Y 病毒侵染烟草的促进作用
龚明月1，段啸天1，余婷婷1，王杰2，申莉莉2，李莹2，刘明宏3，李永亮4，吕洪坤5，章松柏1，杨金广2

（1 长江大学农学院，湖北荆州434025；2 中国农业科学院烟草研究所，山东青岛266101；3 贵州省烟草公司遵义市公司，贵州遵义563000；

4 云南省烟草公司保山市公司，云南保山678000；5 中国烟草总公司海南省公司海南雪茄研究所，海口571100）

摘要：【目的】马铃薯Y 病毒（Potato virus Y，PVY）是危害烟草、马铃薯、番茄、辣椒等主要农作物的重要病毒，给农业生产带来巨大损失。论文旨在克隆Hsc70-2 蛋白在本氏烟（Nicotiana benthamiana）中的基因序列并分析其生物学信息，研究NbHsc70-2 蛋白对PVY 侵染烟草的影响，为进一步解析PVY 的侵染机制提供理论依据。【方法】以本氏烟为材料，克隆NbHsc70-2 蛋白的基因编码序列（coding sequence，CDS），利用MEGA 6.0 进行多序列比对并构建系统进化树；利用实时荧光定量PCR（qRT-PCR）分析NbHsc70-2 在本氏烟各组织中的表达水平；采用在线软件BaCeILo 和SignalP 4.0 对NbHsc70-2 进行生物信息学分析；构建NbHsc70-2-RFP 融合蛋白确定该蛋白的亚细胞定位；研究PVY 处理对本氏烟叶片中NbHsc70-2 的影响；利用激光共聚焦观察PVY-GFP 处理对NbHsc70-2-RFP 定位的影响；构建NbHsc70-2 VIGS 沉默体系及瞬时表达载体，采用qRT-PCR 比较NbHsc70-2在沉默/过表达后对PVY 表达的影响。【结果】NbHsc70-2 编码649 个氨基酸，系统进化树分析表明，NbHsc70-2与普通烟NaHsc70-2 序列相似性最高，亲缘关系最近，属于热激蛋白家族，C 端具有热激蛋白家族的高度保守基序结构；qRT-PCR 分析表明，NbHsc70-2 在叶中表达量最高，在根和茎中的表达量较低；BaCeILo 预测及激光共聚焦显微镜观察均显示NbHsc70-2 定位在细胞质中，并且PVY-GFP 侵染本氏烟后NbHsc70-2-RFP 部分转移至细胞核，与PVY-GFP 共定位在细胞质及细胞核中；将NbHsc70-2 基因沉默载体pTRV::NbHsc70-2 导入本氏烟中7 d后相比于对照组，沉默组出现心叶皱缩及矮化现象；接种PVY-GFP 7 d 后通过手持紫外灯观察到沉默组中无明显绿色荧光出现，而对照组中PVY-GFP 系统侵染使得非接种叶片呈现荧光现象，沉默组荧光点数约为对照组28%；接种PVY 后，利用qRT-PCR 检测沉默NbHsc70-2 后1、3、5 d 的PVY CP 基因积累量，结果表明沉默NbHsc70-2后PVY CP 基因积累量下降，其中3 d 和5 d 与对照相比差异显著，基因表达水平分别为对照的14%和0.004%；将NbHsc70-2 过表达载体pEarleyGate100::GWC::NbHsc70-2 与PVY 同时导入本氏烟后，结果表明相比于对照组，过表达组48 h 和72 h 其PVY CP 基因积累量显著上升，基因表达水平分别为对照组的2.31、2.56 倍。【结论】PVY 侵染会引起NbHsc70-2 表达量上升；NbHsc70-2 是PVY 侵染本氏烟重要的组成成分，沉默该基因显著抑制PVY的表达，过表达NbHsc70-2 则显著提高PVY 的表达，即NbHsc70-2 的表达水平与PVY 复制呈正相关，NbHsc70-2蛋白促进PVY 对烟草的侵染。

关键词：本氏烟；NbHsc70-2 蛋白；马铃薯Y 病毒；基因沉默；过表达；侵染复制

0 引言
【研究意义】烟草马铃薯Y 病毒病又称“烟草脉坏死病”或“烟草脉带病”，是由马铃薯Y 病毒（Potato virus Y，PVY）侵染烟草引起的一种系统性侵染病害[1-2]。PVY 于1931 年由SMITH 首次发现[3]，属于马铃薯Y 病毒属（Potyvirus）成员，是危害烟草、马铃薯、番茄、辣椒等茄科作物的一种病毒。马铃薯Y 病毒属病毒为正链RNA 病毒，分布广泛，几乎在全世界都有发生，尤其在热带和亚热带地区更为流行，给农业生产带来巨大损失[4]。因此，研究PVY 的致病机制对于烟草马铃薯Y 病毒病的防治具有重要意义。【前人研究进展】正链RNA 病毒通过诱导宿主蛋白来促进其在受感染宿主中的高效复制。在模型RNA病毒——番茄丛矮病毒（Tomato bushy stunt virus，TBSV）致病机理的研究中发现，该病毒的复制酶复合物中含有热激蛋白70（heat shock protein 70，Hsp70），是TBSV 复制所必需的一种丰富的细胞质伴侣，这种胞质Hsp70 在TBSV 复制过程中发挥着多种作用，如影响病毒复制蛋白的亚细胞定位和膜插入，以及病毒复制酶的组装[5]。并且通过蛋白质组学分析证明，Hsp70 是TBSV 病毒复制复合体（virus replication complex，VRC）的一种组成成分[6]。热激蛋白（Hsp）是植物细胞内重要的分子伴侣，该蛋白在进化上高度保守，按照分子量大小的不同分为6个家族，即Hsp110s、Hsp90s、Hsp70s、Hsp60s、小分子sHsps 和泛素[7]，其中研究最多的是Hsp70家族。根据热激蛋白70 的表达类型，可分为组成型Hsc70 和诱导型Hsp70，前者在正常条件下就有本底表达，以维持细胞正常生命活动的需要[8-9]；而后者在受环境因素胁迫时诱导表达，以应付外界恶劣环境条件的威胁[10]；Hsc70 广泛分布于原核及真核生物细胞中，是维持生命活动的重要功能蛋白，具有多项生物学功能：协助新生蛋白的折叠，参与细胞内蛋白质的转运，参与免疫复合物的形成和分解及降解冗余蛋白[11-12]。【本研究切入点】Hsp70 蛋白家族在多种病毒和寄主互作中的作用已有大量研究报道，而Hsc70 对植物病毒尤其是对PVY 侵染复制的影响并不明确。前期研究发现本氏烟（Nicotiana benthamiana）在PVY 侵染后NbHsc70-2 mRNA 水平显著上升，推测NbHsc70-2 可能参与PVY 的侵染过程。【拟解决的关键问题】以本氏烟为材料，通过开展生物信息学分析、亚细胞定位、组织特异性表达和基因功能验证等研究，探索NbHsc70-2 在PVY 侵染复制中的作用，为烟草抗马铃薯Y 病毒病研究提供理论依据。

1 材料与方法
试验于2018—2019 年在中国农业科学院烟草研究所植物保护研究中心进行。

1.1 植物材料与病毒
供试植物材料为本氏烟，本氏烟种子点播于混合土中（土壤﹕草炭= 1﹕1），于25℃人工气候室中培养，光周期为14 h 光照/10 h 黑暗，相对湿度为70%，试验采用5—6 叶期烟苗。

PVY 毒源为中国农业科学院烟草研究所保存的N株系。取1 g PVY 野生型毒源植株的叶片，放入提前灭菌消毒的研钵中，加入40 mL PBS 缓冲液，研磨成浆并用纱布过滤除去叶片残渣获得悬浮液以备用；浸润PVY 病毒汁液于本氏烟烟苗（4 周）的第3—4 片真叶，每片叶200 µL。

PVY-GFP 属于PVY 坏死株系（PVYN）的侵染性克隆，是通过无义突变在PVY 侵染性克隆的P1 与HC-pro 间插入一个SacⅡ限制性内切酶酶切位点，并在该酶切位点中插入表达了外源绿色荧光蛋白（GFP）的重组病毒，该病毒能系统侵染本氏烟，病毒扩增至植株非接种部位而呈现荧光[13]。

1.2 NbHsc70-2 的克隆与分析
本氏烟RNA 提取按照TaKaRa RNAiso Plus（Total RNA 提取试剂）试剂盒手册进行，并反转录为cDNA（Vazyme），根据拟南芥中已报道的Hsc70 核苷酸序列在美国本氏烟草数据库（https://solgenomics.net/ organism/nicotiana_benthamiana/genome）中比对得到一个描述为Heat shock 70 kD protein 基因（基因号：Niben101Scf00449g06008.1），根据该基因序列设计特异性引物NbHsc70-2 F/R（表1），以本氏烟cDNA 为模板PCR 扩增NbHsc70-2；PCR 产物连接pEasy-Blunt Simple 载体并转化E. coil 感受态菌株（TaKaRa），阳性克隆委托派森诺生物科技有限公司测序；利用MEGA 6.0[14]进行核苷酸多序列比对分析后，基于最低的贝叶斯标准（Bayesian information criterion，BIC）确定最适合核苷酸替代模型为木村双参数模型（Kimura two-parameter model，K2P），应用于最大似然法（maximum-likelihood，ML）构建系统发育树。各分支的置信度通过1 000次自举抽样（bootstrapping）进行评估。

1.3 NbHsc70-2 蛋白的亚细胞定位
利用BaCelLo[15]和SignalP 4.0[16]进行NbHsc70-2的亚细胞定位和信号肽预测；根据NbHsc70-2 序列设计不含有终止密码子的带XbaI/KpnI 酶切位点的特异性引物NbHsc70-2-XbaI F/KpnI R（表1），并以本氏烟 cDNA 为模板进行 PCR 扩增，利用In-Fusion 技术（TaKaRa）将产物与Fu46 载体相连，构建NbHsc70-2::RFP 融合基因，再利用LR 同源重组技术（Invitrogen）将融合基因连接到pEarleyGate100载体，最终得到pEarleyGate100::NbHsc70-2::RFP 载体；将含有pEarleyGate100::NbHsc70-2::RFP 载体和pEarleyGate100::RFP 载体的 LBA4404 农杆菌（OD600=0.8）浸润本氏烟下表皮，培养于25℃人工气候室中，光周期为14 h 光照/10 h 黑暗，相对湿度为70%，72 h 后制片置于激光共聚焦显微镜（SP8，Leica）下观察明、暗视野下的表达情况。

表1 试验中采用的引物及序列
Table 1 Sequence of the primers used in this study

pagenumber_ebook=116,pagenumber_book=774
1.4 NbHsc70-2 沉默载体的构建
根据NbHsc70-2 基因序列，用Premier 5.0 设计一对含有限制性核酸内切酶EcoRI、KpnI 的特异性引物RNAi NbHsc70-2 F/R（表1），PCR 扩增大小为350 bp目的沉默片段，利用In-Fusion 技术（TaKaRa）将产物与pTRV2 载体相连，构建pTRV::NbHsc70-2 重组载体；将含有pTRV::NbHsc70-2 载体、pTRV::PDS 及pTRV00 空载体的LBA4404 农杆菌（OD600=0.8）注射本氏烟下表皮，7 d 后检测沉默效率，进行后续试验。

1.5 NbHsc70-2 瞬时过表达载体的构建
根据NbHsc70-2基因序列设计含有AhdI酶切位点的特异性引物NbHsc70-2-AhdI F/R（表1），并以本氏烟cDNA 为模板进行PCR 扩增，利用In-Fusion 技术（TaKaRa）将产物与 GWC 载体相连，构建GWC::NbHsc70-2 入门载体，再利用LR 同源重组原理（Invitrogen）将入门载体连接到pEarleyGate100载体，最终得到pEarleyGate100::GWC::NbHsc70-2表达载体；提前24 h 给本氏烟浸润PVY 病汁液，将含有pEarleyGate100::GWC::NbHsc70-2 表达载体和pEarleyGate100::GWC 载体的LBA4404 农杆菌（OD600=0.8）浸润本氏烟下表皮，48 h 后开始采样检测PVY CP 基因积累量。

1.6 qRT-PCR 与统计分析
根据不同样品中特定基因的基因序列设计各基因 的荧光定量检测引物（表1），同时设置Actin 作为内参。分别提取各取样时间点的处理组及对照组材料的总RNA（TaKaRa），反转录成cDNA（Vazyme）。以cDNA 为模板，按照ChamQTM Universal SYBR® qPCR Master Mix（Vazyme）试剂盒说明书进行操作，在Applied Biosystems 7500 real-time PCR 机型上进行扩增。反应体系：10 µL 2×ChamQ Universal SYBR qPCR Master Mix，0.4 µL Primer1（10 µmol·L-1），0.4 µL Primer2（10 µmol·L-1），2 µL cDNA，补蒸馏水至总体积20 µL。反应程序：95℃ 30 s；95℃ 5 s；60℃ 30 s；40 循环；95℃ 15 s；60℃ 1 min，95℃ 15 s。相对表达量采用2-ΔΔCt 法计算[17-18]，使用7500 Software v2.3 进行数据统计分析并使用GraphPad Prism6.0 绘图。P＜0.05 表示差异显著，P＜0.01 表示差异极显著。

2 结果
2.1 NbHsc70-2 克隆与序列分析
根据拟南芥中已报道的Hsc70 核苷酸序列在美国本氏烟草数据库中比对得到一个描述为Heat shock 70 kD protein 基因（基因号：Niben101Scf00449g06008.1）。为鉴定Niben101Scf00449g06008.1 基因，在NCBI（美国国立生物技术信息中心数据库）中比对发现Niben101Scf00449g06008.1 与NaHsc70-2、NsHsc70-2、NtHsc70-2 序列相似性最高，分别有69.07%、69.07%、68.25%核苷酸相似性，但NaHsc70-2、NsHsc70-2、NtHsc70-2 的编码序列（coding sequence，CDS）均有1 953 bp，而Niben101Scf00449g06008.1 基因的CDS 序列中只有1 443 bp，因此推测Niben101Scf00449g06008.1 基因CDS 序列部分缺失，通过比对NCBI 上已发表的Hsc70 序列发现该基因非常保守，设计特异性引物，通过PCR扩增得到Niben101Scf00449g06008.1 基因的完整CDS 序列；将该CDS 序列与NaHsc70-2、NsHsc70-2、NtHsc70-2 序列进行比对，发现分别具有94.49%、94.49%、93.78%核苷酸相似性，对应的蛋白质间分别具有97.38%、97.38%、97.23%氨基酸相似性。通过获取与Niben101Scf00449g06008.1 相似度高于96%的基因CDS 序列，并用Clustal W 进行核酸序列比对，利用ML 算法构建系统发育树，结果均表明Niben101Scf00449g06008.1 与NaHsc70-2 关系最近（图1），暗示Niben101Scf00449g06008.1 可能是Hsc70-2 在本氏烟中的同源基因。从图2 中可以看出，Niben101Scf00449g06008.1 与NaHsc70-2 蛋白氨基酸序列相似性最高，表明Niben101Scf00449g06008.1 是Hsc70-2 在本氏烟中的同源基因，命名为NbHsc70-2。

pagenumber_ebook=117,pagenumber_book=775
图1 NbHsc70-2 及其同源基因系统发育分析（ML 算法）
Fig. 1 Phylogenetic analyses of NbHsc70-2 and its homologues (ML algorithm)

2.2 NbHsc70-2 组织表达特异性
以Actin 作为内参基因，利用qRT-PCR 分别检测NbHsc70-2 在本氏烟根、茎、叶中的表达情况。结果显示NbHsc70-2 在叶中的相对表达量较高，在根和茎中的表达量很低，其中叶片相对表达量是根中的1.91倍（图3）。因此，后期试验均采集叶片样品进行研究。

2.3 PVY 侵染对NbHsc70-2 的影响
采用qRT-PCR 检测PVY 侵染1、3 d 后NbHsc70-2 mRNA 水平的变化，结果显示PVY 侵染后NbHsc70-2 mRNA 水平明显上升，与对照相比，处理组1 d 无明显差异，PVY 处理3 d 后NbHsc70-2 mRNA 水平上调，处理组是对照组的9.24 倍（图4），表明PVY 侵染会引起NbHsc70-2 mRNA 水平上调，暗示NbHsc70-2 可能在PVY 侵染过程中发挥重要作用。

利用BaCelLo 预测NbHsc70-2 蛋白的亚细胞定位，发现NbHsc70-2 蛋白定位于细胞质上的可能性高于其他部位。SignalP4.0 Server 信号肽预测结果显示该蛋白无信号肽，即NbHsc70-2 蛋白不属于分泌蛋白。为验证亚细胞定位的预测，将NbHsc70-2 蛋白与红色 荧光蛋白RFP 构建融合表达载体，注射本氏烟表皮瞬时表达进行亚细胞定位，激光共聚焦显微镜显示融合蛋白与Cytoplasm-FDA（细胞质Marker）定位重合，表明NbHsc70-2 蛋白定位在细胞质中（图5）。

pagenumber_ebook=118,pagenumber_book=776
图2 NbHsc70-2、NaHsc70-2、NtHsc70-2 及NsHsc70-2 蛋白氨基酸序列比对分析
Fig. 2 Amino acid sequences multiple alignment among NbHsc70-2, NaHsc70-2, NtHsc70-2 and NsHsc70-2

深蓝色表示氨基酸完全相同，浅蓝色表示有两个氨基酸不同，粉色表示有一个氨基酸不同 Dark blue indicates that the amino acids are identical, light blue indicates that there are two amino acids differences, and pink indicates that there is one amino acid difference

pagenumber_ebook=118,pagenumber_book=776
图3 NbHsc70-2 的组织表达特异性分析
Fig. 3 Analysis of tissue expression specificity of NbHsc70-2

取同一植株的不同组织样品进行分析Different tissue samples from the same plant were taken for analysis；误差线表示3 组生物学重复的标准差，**表示差异极显著（P＜0.01）。下同 Error bars indicate the standard deviations of three sets of biological replicates, and ** indicates that the difference is extremely significant (P＜0.01). The same as below

pagenumber_ebook=118,pagenumber_book=776
图4 PVY 侵染本氏烟后NbHsc70-2 水平的变化
Fig. 4 Gene accumulation of NbHsc70-2 after PVY infection in N. benthamiana

pagenumber_ebook=119,pagenumber_book=777
图5 PVY 不侵染和PVY 侵染NbHsc70-2 蛋白的亚细胞定位
Fig. 5 Subcellular localization of NbHsc70-2 protein by PVY (B) and not by PVY (A)

A：NbHsc70-2 蛋白的亚细胞定位，Cytoplasm-FDA 为细胞质Marker Subcellular localization of NbHsc70-2 protein, Cytoplasm-FDA is cytoplasmic Marker；B：PVY-GFP 侵染本氏烟对NbHsc70-2 蛋白定位的影响，箭头所指为细胞核The effect of PVY-GFP infecting N. benthamiana on NbHsc70-2 protein localization, the arrow is referred to the nucleus

利用PVY-GFP 融合表达载体与NbHsc70-2-RFP融合表达载体共浸润本氏烟叶片 48 h 后，观察PVY-GFP 侵染对NbHsc70-2 蛋白定位的影响。结果显示，PVY-GFP 侵染后，NbHsc70-2-RFP 融合蛋白的细胞质定位发生变化，部分转移至细胞核中，NbHsc70-2-RFP 融合蛋白与PVY-GFP 共定位在细胞质及细胞核中，推测NbHsc70-2 可能是PVY 病毒复制复合体的一种组成成分，参与胞内移动过程。

2.4 沉默NbHsc70-2 对PVY 积累量的影响
为了进一步研究NbHsc70-2 编码的Hsc70-2 蛋白在PVY 侵染本氏烟中的功能，利用病毒诱导基因沉默技术构建VIGS 沉默体系下调本氏烟中NbHsc70-2 的含量，分析寄主中NbHsc70-2 的含量降低后对PVY侵染本氏烟的影响。由于NbHsc70-2 具有高度的CDS核酸相似性，因此选择NbHsc70-2 相对不保守区设计沉默靶标。同时，根据不保守区设计qRT-PCR 引物以扩增NbHsc70-2，用于检测基因沉默效率（表1）。含NbHsc70-2 沉默重组质粒农杆菌浸润本氏烟7 d 后，对NbHsc70-2 沉默效率为（67.50±0.29）%（图6-B）。从图6-A 可以看出沉默NbHsc70-2 导致心叶皱缩及植株矮化，暗示NbHsc70-2 编码的NbHsc70-2 蛋白可能参与调节植株生长发育。

含不同重组质粒农杆菌浸润寄主植株7 d 后，将PVY-GFP 接种于TRV 浸润叶的上面第3 片叶，接种7 d 后，在手提式紫外灯下观察本氏烟中的绿色荧光情况；在TRV 对照组中PVY-GFP 系统侵染使得非接种叶片呈现荧光现象，而NbHsc70-2 沉默组中无明显绿色荧光出现（图6-C），荧光点数约是对照组的28%（图6-D）。表明本氏烟中NbHsc70-2 的mRNA 含量降低，可导致PVY-GFP 系统性侵染受到抑制。

为了进一步明确下调NbHsc70-2 表达后抑制PVY侵染本氏烟，浸润接种pTRV::NbHsc70-2 后7 d 接种野生型病毒PVYN，从RNA 水平检测病毒的相对累积量。从图7-A 可以看出沉默组与对照组的PVY CP 表达量都呈现上升趋势，但沉默组的上升趋势较对照组慢；在沉默组上接种PVY 后 1 d，CP 表达量较对照组已有下降；3 d 沉默组中的表达量是对照组的14%；5 d 沉默组中的表达量是对照组的0.004%。表明沉默NbHsc70-2 后，显著延缓植株内PVY CP 的mRNA 积累。

pagenumber_ebook=120,pagenumber_book=778
图6 沉默NbHsc70-2 表型分析
Fig. 6 Phenotype analysis of NbHsc70-2 silencing

A：沉默NbHsc70-2 7 d 后的烟草表型Tobacco phenotype after 7 days of silencing NbHsc70-2；B：沉默NbHsc70-2 7 d 后的沉默效率Silencing efficiency after 7 days of silencing NbHsc70-2；C：接种PVY-GFP 7 d 后的烟草表型Tobacco phenotype after 7 days of PVY-GFP inoculation；D：接种PVY-GFP 7 d 后的烟草相对荧光强度Relative fluorescence intensity of tobacco after 7 days of PVY-GFP inoculation
pTRV1 与pTRV2::NbHsc70-2 或空的pTRV2（作为阴性对照）通过根癌农杆菌共浸润培养3 周的本氏烟。侵染7 d 后统计NbHsc70-2 沉默效率并拍照pTRV1 and pTRV2::NbHsc70-2 or empty pTRV2 (as a negative control) were co-infiltrated with Agrobacterium tumefaciens for 3 weeks of N. benthamiana. The silencing efficiency of NbHsc70-2 was counted 7 days after infection and photos were taken

pagenumber_ebook=120,pagenumber_book=778
图7 沉默及过表达NbHsc70-2 对PVY 积累的影响
Fig. 7 Effects of NbHsc70-2 silencing and overexpression on PVY accumulation

A：沉默NbHsc70-2 后对PVY 积累的影响Effect of silencing NbHsc70-2 on PVY accumulation；B：过表达NbHsc70-2 后对PVY 积累的影响Effect of overexpression of NbHsc70-2 on PVY accumulation
pTRV1 与pTRV2::NbHsc70-2 或空的pTRV2（作为阴性对照）通过根癌农杆菌共浸润培养3 周的本氏烟。侵染7 d 后接种PVY 开始取样进行检测。pEarleyGate100:: NbHsc70-2 或空的pEarleyGate100（作为阴性对照）通过根癌农杆菌共浸润培养3 周的本氏烟，提前24 h 接种PVY。侵染48 h 后开始取样进行检测pTRV1 and pTRV2::NbHsc70-2 or empty pTRV2 (as a negative control) were co-infiltrated with A. tumefaciens for 3 weeks of N. benthamiana. Seven days after infection, PVY was inoculated and samples were taken for test. pEarleyGate100:: NbHsc70-2 or empty pEarleyGate100 (as a negative control) was inoculated with A. tumefaciens for 3 weeks and was inoculated with PVY 24 hours in advance. Sampling was taken 48 h after infection

2.5 过表达NbHsc70-2 对PVY 积累量的影响
为了进一步确定NbHsc70-2 编码的Hsc70-2 蛋白在PVY 侵染寄主过程中的作用，通过LR 反应将构建好的pGWC-NbHsc70-2 与pEarleyGate100 重组，构建NbHsc70-2 瞬时表达载体，以表达载体pEarleyGate100在本氏烟中过量表达 NbHsc70-2，提高寄主中NbHsc70-2 含量，分析PVY 侵染植株后RNA 水平上的病毒累积量。

从图7-B 可以看出经过表达载体处理48 h 后，过表达组的PVY CP 表达量较对照组已有上升，过表达组中的表达量是对照组的2.31 倍；处理72 h 后过表达组中的PVY CP 表达量是对照组的2.56 倍。表明过表达寄主中NbHsc70-2 的含量能显著增加植株内PVY CP 的积累。

3 讨论
植物热激蛋白70 是热激蛋白家族中保守、普遍表达的，家族中通常具有两个主要功能区：N 端核酸结合区（nucleotide binding domain，NBD）和C 端底物结合区（substrate binding domain，SBD）。N 末端核酸结合区具有结合并水解ATP 的活性中心，它与ATP的结合-水解过程调控其对多肽的亲和力[19]。C 端存在高度保守基序结构（motif），通常用于指示亚细胞定位信号或与其他伴侣或共同伴侣相互作用[20-21]。例如定位于叶绿体中的Hsp70 其motif 结构的氨基酸序列为 PECDVLDADFTDSK，定位于线粒体中的为PEAEYEEAKK，定位于内质网中的为HDEL，而定位于细胞质中的为EEVD[22]。本试验中，在美国本氏烟草数据库中筛选到的Niben101Scf00449g06008.1 基因经鉴定为NbHsc70-2，编码649 个氨基酸，与已报道的NaHsc70-2 蛋白的相似性达到97.38%，系统进化树也显示NbHsc70-2 与NaHsc70-2 基因亲缘性最高，并且该蛋白C 端存在EEVD 基序，亚细胞定位结果显示该蛋白也定位于细胞质中，这说明NbHsc70-2 是Hsc70-2 在本氏烟中的同源基因。

Hsc70 蛋白是一种结构型Hsp70，在所有的细胞内均可表达，是一种ATP 结合蛋白[23]，其在非生物胁迫下无明显表达差异，如热刺激后表达量不增加或减少增加；在生物胁迫中，如在病毒的侵染过程中，Hsc70 被认为影响病毒蛋白的复制、积累、运动和折叠。许多病毒将Hsc70 招募到膜复制工厂，如红花苜蓿坏死花叶病毒（Red clover necrotic mosaic virus，RCNMV）[24]、芜菁花叶病毒（Turnip mosaic virus，TuMV）[25]和中国小麦花叶病毒（Chinese wheat mosaic virus，CWMV）[26]。WANG 等[27]研究表明，宿主Hsc70是乙型肝炎病毒（Hepatitis B virus，HBV）复制的辅助因子，Hsc70 的表达水平与HBV 的复制呈正相关性；MATHIOUDAKIS 等[28]研究表明，Hsc70-3 mRNA水平及蛋白水平在番茄花叶病毒（Pepino mosaic virus，PepMV）侵染过程中被诱导增加；有研究表明，Hsp70 家族与VRC 相关，能增强病毒RNA 的复制；例如Hsp70 在番茄丛矮病毒（Tomato bushy stunt virus，TBSV）侵染复制中是TBSV VRC 的一种组成成分[5-6]；在TuMV 中Hsc70-3 是膜结合复制复合体的一部分[25]。此外Hsc70 还与病毒基因组成分相互作用，以增强病毒在宿主植株中的RNA 积累；Hsp72在丙肝病毒（Hepatitis C virus，HCV）复制中与复制蛋白NS5A、NS5B（RNA-dependent RNA polymerase，RdRp）互作[29]。Hsc70-3 在PepMV 复制中与PepMV CP 相互作用；并且Hsc70 还参与病毒细胞间转运，在苘麻花叶病毒（Abutilon mosaic virus，AbMV）中运动蛋白（movement protein，MP）与Hsc70 相互作用[30]，表明Hsc70 可能在病毒转运方面具有一定的功能；在PepMV 中Hsc70-3 与PepMV CP 相互作用复合物出现在细胞质和细胞核中，表明Hsc70-3 参与PepMV 的细胞内移动[28]；ALAM 等[31]研究表明，Hsc70-2 与纯化的黄瓜坏死病毒（Cucumber necrosis virus，CNV）病毒粒子结合，并在病毒颗粒的解体中发挥作用。

本研究发现PVY 侵染的本氏烟中NbHsc70-2 mRNA 水平上升，且引起部分NbHsc70-2 蛋白定位转移至细胞核中，NbHsc70-2 蛋白与PVY 共定位在细胞质及细胞核中；在AbMV 中MP 与Hsc70 相互作用，在PepMV 中也发现Hsc70-3 参与PepMV 的胞内移动；推测NbHsc70-2 可能参与PVY 的胞内移动过程。本研究利用VIGS 技术下调宿主中NbHsc70-2 含量，进一步发现该基因缺失后，能够显著抑制PVY CP 的积累量，表明NbHsc70-2 的表达水平与PVY 的复制呈正相关，沉默组与对照组相比下降幅度较大可能是由于NbHsc70-2 是PVY 侵染复制所需要的重要组成部分，NbHsc70-2 是PVY VRC 的一部分，PVY 侵染需要招募宿主NbHsc70-2 蛋白到膜复制工厂形成PVY病毒复制复合体，缺失NbHsc70-2 导致PVY 无法形成VRC 进而导致沉默组与对照组的差异较大；并且病毒的复制是快速增长的，这也是导致沉默组与对照组差异较大的原因；同时沉默NbHsc70-2 后7 d 植株出现皱缩及矮化症状，这也暗示了NbHsc70-2 在植物生长发育中具有重要的作用。但NbHsc70-2 功能的丧失可能造成了植株的一些生理缺陷，因此本研究进一步上调宿主中NbHsc70-2 的含量，结果表明过表达NbHsc70-2 能够提高PVY CP 的积累量，这与WANG等[27]的研究结果一致。

4 结论
NbHsc70-2 为Hsp70 家族蛋白，PVY 侵染后导致其基因表达上升。NbHsc70-2 是响应PVY 侵染的重要基因，正向调控PVY 病毒的复制，抑制该基因的表达可以提高植株对PVY 的抗性，该基因在烟草抗马铃薯Y 病毒病方面具有较大的应用价值，未来可将该蛋白作为一个抗病毒药剂靶标。

References

[1] 张秀丽. 烟草马铃薯Y 病毒综合防治技术. 黑龙江农业科学, 2008(3): 151. ZHANG X L. Comprehensive control technology of Potato virus Y (PVY) in tobacco. Heilongjiang Agricultural Sciences, 2008(3): 151. (in Chinese)

[2] 余春英, 张西仲, 王定福, 李余湘, 罗红香, 蒙祥旭, 曾宪佛. 黔南烟区烟草马铃薯Y 病毒病发病原因及防治措施. 安徽农业科学, 2010, 38(10): 5110-5112. YU C Y, ZHANG X Z, WANG D F, LI Y X, LUO H X, MENG X X, ZENG X F. Causes and prevention measures of tobacco Potato virus Y disease in Qiannan tobacco growing areas. Journal of Anhui Agricultural Science, 2010, 38(10): 5110-5112. (in Chinese)

[3] SMITH K M. On the composite nature of certain potato virus diseases of the mosaic group as revealed by the use of plant indicators and selective methods of transmission. Proceedings of the Royal Society of London Series B, Containing Papers of a Biological Character, 1931, 109(762): 251-267.

[4] 崔晓艳, 魏太云, 陈新. 马铃薯Y 病毒属病毒的细胞生物学研究进展. 中国农业科学, 2012, 45(7): 1293-1302. CUI X Y, WEI T Y, CHEN X. Advances in research on cell biology of Potyvirus. Scientia Agricultura Sinica, 2012, 45(7): 1293-1302. (in Chinese)

[5] WANG R Y, STORK J, NAGY P D. A key role for heat shock protein 70 in the localization and insertion of tombusvirus replication proteins to intracellular membranes. Journal of Virology, 2009, 83(7): 3276-3287.

[6] SERVA S, NAGY P D. Proteomics analysis of the tombusvirus replicase: Hsp70 molecular chaperone is associated with the replicase and enhances viral RNA replication. Journal of Virology, 2006, 80(5): 2162-2169.

[7] 任月, 韩莹琰, 李婷, 郝敬虹, 范双喜. 叶用莴苣热激蛋白90 (LsHsp90)基因的克隆及其在热激下的表达. 中国农业科学, 2013, 46(16): 3514-3522. REN Y, HAN Y Y, LI T, HAO J H, FAN S X. Molecular cloning and expression analysis of heat-shock-protein 90 (LsHsp90) gene from leaf lettuce (Lactuca sativa L.) under heat shock. Scientia Agricultura Sinica, 2013, 46(16): 3514-3522. (in Chinese)

[8] SUNG D Y, VIERLING E, GUY C L. Comprehensive expression profile analysis of the Arabidopsis Hsp70 gene family. Plant Physiology, 2001, 126: 789-800.

[9] SWINDELL W R, HUEBNER M, WEBER A P. Transcriptional profiling of Arabidopsis heat shock proteins and transcription factors reveals extensive overlap between heat and non-heat stress response pathways. BMC Genomics, 2007, 8: 125.

[10] BRIZARD J P, CARAPITO C, DELALANDE F, VAN DORSSELAER A, BRUGIDOU C. Proteome analysis of plant-virus interactome: Comprehensive data for virus multiplication inside their hosts. Molecular and Cellular Proteomics, 2006, 5(12): 2279-2297.

[11] 王利彬, 董爽爽, 王灿国, 程敦公, 李豪圣, 刘爱峰, 宋健民, 刘建军, 刘成, 张玉梅, 穆平, 赵振东, 曹新有. 黑麦热激蛋白ScHsp90-1 基因的克隆及表达分析. 山东农业科学, 2018, 50(3): 1-5. WANG L B, DONG S S, WANG C G, CHENG D G, LI H S, LIU A F, SONG J M, LIU J J, LIU C, ZHANG Y M, MU P, ZHAO Z D, CAO X Y. Cloning and expression analysis of heat shock protein gene ScHsp90-1 in rye. Shandong Agricultural Sciences, 2018, 50(3): 1-5. (in Chinese)

[12] WANG N B, ZHAO J, HE X Y, SUN H Y, ZHANG G P, WU F B. Comparative proteomic analysis of drought tolerance in the two contrasting Tibetan wild genotypes and cultivated genotype. BMC Genomics, 2015, 16: 432.

[13] SUN H J, SHEN L L, QIN Y X, LIU X W, HAO K Q, LI Y, WANG J, YANG J G, WANG F L. CLC-Nt1 affects Potato virus Y infection via regulation of endoplasmic reticulum luminal Ph. New Phytologist, 2018, 220(2): 539-552.

[14] TAMURA K, STECHER G, PETERSON D, FILIPSKI A, KUMAR S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Molecular Biology and Evolution, 2013, 30(12): 2725-2729.

[15] PIERLEONI A, MARTELLI P L, FARISELLI P, CASADIO R. BaCelLo: A balanced subcellular localization predictor. Bioinformatics, 2006, 22(14): e408-e416.

[16] PETERSEN T N, BRUNAK S, VON HEIJNE G, NIELSEN H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nature Methods, 2011, 8(10): 785-786.

[17] YE C, DICKMAN M B, WHITHAM S A, PAYTON M, VERCHOT J. The unfolded protein response is triggered by a plant viral movement protein. Plant Physiology, 2011, 156(2): 741-755.

[18] YE C M, CHEN S, PAYTON M, DICKMAN M B, VERCHOT J. TGBp3 triggers the unfolded protein response and SKP1-dependent programmed cell death. Molecular Plant Pathology, 2013, 14(3): 241-255.

[19] FLAHERTY K M, WILBANKS S M, DELUCA-FLAHERTY C, MCKAY D B. Structural basis of the 70-kilodalton heat shock cognate protein ATP hydrolytic activity. II. Structure of the active site with ADP or ATP bound to wild type and mutant ATPase fragment. The

Journal of Biological Chemistry, 1994, 269(17): 12899-12907.

[20] DAUGAARD M, ROHDE M, JAATTELA M. The heat shock protein 70 family: Highly homologous proteins with overlapping and distinct functions. FEBS Letters, 2007, 581(19): 3702-3710.

[21] MAYER M P, BUKAU B. Hsp70 chaperones: cellular functions and molecular mechanism. Cellular and Molecular Life Sciences, 2005, 62(6): 670-684.

[22] GUY C L, LI Q B. The organization and evolution of the spinach stress 70 molecular chaperone gene family. The Plant Cell, 1998, 10(4): 539-556.

[23] PEREZ-VARGAS J, ROMERO P, LOPEZ S, ARIAS C F. The peptide-binding and ATPase domains of recombinant hsc70 are required to interact with rotavirus and reduce its infectivity. Journal of Virology, 2006, 80(7): 3322-3331.

[24] MINE A, HYODO K, TAJIMA Y, KUSUMANEGARA K, TANIGUCHI T, KAIDO M, MISE K, TANIGUCHI H, OKUNO T. Differential roles of Hsp70 and Hsp90 in the assembly of the replicase complex of a positive-strand RNA plant virus. Journal of Virology, 2012, 86(22): 12091-12104.

[25] DUFRESNE P J, THIVIERGE K, COTTON S, BEAUCHEMIN C, IDE C, UBALIJORO E, LALIBERTE J F, FORTIN M G. Heat shock 70 protein interaction with Turnip mosaic virus RNA-dependent RNA polymerase within virus-induced membrane vesicles. Virology, 2008, 374(1): 217-227.

[26] YANG J, ZHANG F, CAI N J, WU N, CHEN X, LI J, MENG X F, ZHU T Q, CHEN J P, ZHANG H M. A furoviral replicase recruits host HSP70 to membranes for viral RNA replication. Scientific Reports, 2017, 7: 45590.

[27] WANG Y P, LIU F, HE H W, HAN Y X, PENG Z G, LI B W, YOU X F, SONG D Q, LI Z R, YU L Y, CEN S, HONG B, SUN C H, ZHAO L X, KREISWIRTH B, PERLIN D, SHAO R G, JIANG J D. Heat stress cognate 70 host protein as a potential drug target against drug resistance in hepatitis B virus. Antimicrobial Agents and Chemotherapy, 2010, 54(5): 2070-2077.

[28] MATHIOUDAKIS M M, VEIGA R, GHITA M, TSIKOU D, MEDINA V, CANTO T, MAKRIS A M, LIVIERATOS I C. Pepino mosaic virus capsid protein interacts with a tomato heat shock protein cognate 70. Virus Research, 2012, 163(1): 28-39.

[29] CHEN Y J, CHEN Y H, CHOW L P, TSAI Y H, CHEN P H, HUANG C Y, CHEN W T, HWANG L H. Heat shock protein 72 is associated with the hepatitis C virus replicase complex and enhances viral RNA replication. The Journal of Biological Chemistry, 2010, 285(36): 28183-28190.

[30] KRENZ B, WINDEISEN V, WEGE C, JESKE H, KLEINOW T. A plastid-targeted heat shock cognate 70 kDa protein interacts with the Abutilon mosaic virus movement protein. Virology, 2010, 401(1): 6-17.

[31] ALAM S B, ROCHON D. Cucumber necrosis virus recruits cellular heat shock protein 70 homologs at several stages of infection. Journal of Virology, 2016, 90(7): 3302-3317.

Cloning of Hsc70-2 and Its Promoting Effect on Potato virus Y Infection in Nicotiana benthamiana

GONG MingYue1, DUAN XiaoTian1, YU TingTing1, WANG Jie2, SHEN LiLi2, LI Ying2, LIU MingHong3,LI YongLiang4, LÜ HongKun5, ZHANG SongBai1, YANG JinGuang2

(1College of Agronomy, Yangtze University, Jingzhou 434025, Hubei; 2Tobacco Research Institute, Chinese Academy of Agricultural Sciences, Qingdao 266101, Shandong; 3Zunyi City Company, Guizhou Tobacco Company, Zunyi 563000, Guizhou; 4Baoshan City Company, Yunnan Tobacco Company, Baoshan 678000, Yunnan; 5Hainan Cigar Research Institute, Hainan Provincial Corporation, China National Tobacco Corporation, Haikou 571100)

Abstract:【Objective】Potato virus Y (PVY) is an important virus that harms major crops such as tobacco, potato, tomato, pepper and so on, and brings great losses to agricultural production. The objective of this study is to clone NbHsc70-2, analyze the bioinformatics of NbHsc70-2 protein, research the effect of NbHsc70-2 protein on PVY infection in Nicotiana benthamiana, and to provide a theoretical basis for further analyzing the infection mechanism of PVY.【Method】The gene coding sequence (CDS) of Hsc70-2 protein from N. benthamiana was cloned in pEarleyGate100 vector. Multi-sequence alignment and phylogenetic analysis were performed using MEGA 6.0 software. Fluorescence quantitative RT-PCR (qRT-PCR) analysis of NbHsc70-2 expression level was performed in different tissues of N. benthamiana. Bioinformatics analysis using online software BaCeILo and SignalP 4.0.was performed for construction of NbHsc70-2-RFP fusion protein and determination of the subcellular localization of the identified NbHsc70-2 protein. PVY-GFP was used to analyze the effect of PVY infection on NbHsc70-2 in N. benthamiana leaves and localization of this protein. Furthermore, NbHsc70-2 VIGS silencing system and transient expression vector were constructed, and the effect of the protein on PVY expression after silencing/overexpression was analyzed by comparing the results with qRT-PCR assay.【Result】NbHsc70-2 encodes a total of 649 amino acids. The phylogenetic analysis indicated that NbHsc70-2 belongs to the heat shock protein (HSP) family and has the highest similarity with the NaHsc70-2. The C-terminal has the highly conserved motif structure of the HSP family. qRT-PCR analysis showed that NbHsc70-2 had the highest expression level in leaves and low expression in roots and stems. Moreover, BaCeILo prediction and laser confocal microscopy depicted that NbHsc70-2 was localized in the cytoplasm, and after PVY-GFP infection to N. benthamiana, NbHsc70-2-RFP was partially transferred to the nucleus and colocalized with virus in the cytoplasm and nucleus. The induction of NbHsc70-2 gene silencing vector pTRV::NbHsc70-2 into N. benthamiana resulted in stunted growth of infected plants compared to the control at 7 dpi. Whilst, inoculation of silent plants with PVY-GFP didn’t reflect green fluorescence, whereas PVY-GFP inoculation resulted in the green fluorescence production from the leaves of plants kept as control, the number of fluorescent spots in the silent group was about 28% of that in control group. qRT-PCR analysis showed that further inoculation of PVY resulted in accumulation of PVY CP at 1, 3, and 5 d after silencing NbHsc70-2. Moreover, the level of PVY CP decreased after silencing NbHsc70-2, and the difference between 3 and 5 dpi was significant compared to the control, the gene expression level was 14% and 0.004% of that in the control, respectively. Contrarily, high level of PVY CP was observed when N. benthamiana plants were inoculated with NbHsc70-2 overexpression vector along with PVY. The detailed analysis showed that the accumulation of PVY CP was significantly increased at 48 h and 72 h, and the gene expression level was 2.31 and 2.56 times of that in the control group, respectively.【Conclusion】PVY infection causes increased expression of NbHsc70-2. NbHsc70-2 is an important component of PVY infection to N. benthamiana, silencing NbHsc70-2 significantly inhibited PVY expression, and overexpression of NbHsc70-2 significantly increased PVY expression, i.e., the expression level of NbHsc70-2 was positively correlated with PVY replication. NbHsc70-2 protein promoted PVY infection in tobacco.

Key words: Nicotiana benthamiana; NbHsc70-2 protein; Potato virus Y (PVY); gene silencing; overexpression; infection and replication

收稿日期：2019-08-19；

接受日期：2019-09-18

基金项目：烟草绿色防控重大专项（110101601024（LS-04））、云南省烟草公司科技项目（2018530000241014）、海南省烟草公司科技项目（201846000024056）

联系方式：龚明月，E-mail：576092384@qq.com。

通信作者章松柏，E-mail：yangtze2008@126.com。

通信作者杨金广，E-mail：yangjinguang@caas.cn

开放科学（资源服务）标识码（OSID）：

pagenumber_ebook=113,pagenumber_book=771
（责任编辑 岳梅）