首页 > 分享 > 曼地亚红豆杉NPR1基因家族鉴定及表达模式分析

曼地亚红豆杉NPR1基因家族鉴定及表达模式分析

花匠小妙招
2025-03-25 22:02

红豆杉Taxus wallichianavar. chinensis (Pilger) Florin又称为紫杉、赤柏松，由于其能够提取天然抗癌药物紫杉醇而闻名于世。紫杉醇是高效的抗癌药物之一，它对多种癌症，如卵巢癌、乳腺癌、宫颈癌等的治疗有着明显的积极作用[1]，同时也逐渐涉及到对胃癌、食道癌等癌症的治疗[2-3]。紫杉醇虽然广受欢迎，但是其在红豆杉中含量极低，且红豆杉生长极为缓慢，因此造成红豆杉野生资源遭受严重破坏。

水杨酸是植物抗病反应调控中的主要参与者，同时也响应植物参与的其它非生物胁迫。病程相关基因非表达子1（non-expressor of pathogenesis- related genes 1，NPR1）作为植物激素水杨酸信号通路调控的重要组成部分，在水杨酸参与的植物系统获得抗性（systemic acquired resistance，SAR）反应中起着重要的传递作用[4]。同时，它也可以通过作用于某些转录因子来诱导抗性基因的表达，从而激活植物的免疫系统，提高植物的抗性[5-6]。此外，有研究也证实NPR1参与介导茉莉酸、乙烯、脱落酸等其他植物激素信号通路的交互作用[7-8]。

目前，NPR1作为水杨酸信号通路重要调控因子已得到广泛研究，其作用机制已部分阐明。水杨酸通过改变细胞氧化还原势使得多聚体NPR1蛋白解聚为单聚体进入细胞核[9-10]，NPR1在核内与其它转录因子互作，激活或增强防御基因的表达，从而启动多种防御反应[6,11]。同时，NPR1蛋白可以通过磷酸化、泛素化等修饰展现出不同的稳定性与功能[12]。Liu等[13]在拟南芥的基因组中发现了5个NPR1的同源基因，并且将它们分别命名为NPR2、NPR3、NPR4、NPR5和NPR6，之后Fu等[14]发现AtNPR3/4能够直接结合水杨酸发挥负调控作用，而Ding等[15]则发现AtNPR1、AtNPR3/4都可以直接与水杨酸相互作用，平行行使激活抗病基因的功能，协同激活免疫反应。随着NPR1及其同源基因和水杨酸相互调控机制逐步被揭示，NPR1在调控水杨酸信号转导过程中的重要性逐渐被人们认识。

研究发现，除了茉莉酸甲酯（methyl jasmonate，MeJA）以外，水杨酸也会影响红豆杉中紫杉醇的生物合成，从而影响紫杉醇的含量[16]。然而，相关分子机制研究较为少见。本研究通过对曼地亚红豆杉Taxus× media Rehder的NPR1基因家族进行生物信息学以及表达模式分析，为进一步研究水杨酸调控红豆杉紫杉醇合成的分子机制提供理论依据，并为后续利用细胞工程和基因工程提高紫杉醇终产量奠定理论基础。

1材料

曼地亚红豆杉3～5年生小苗种植于中国林业科学研究院科研温室内，并经中国林业科学研究院林业研究所杨艳芳研究员鉴定确认为曼地亚红豆杉Taxus× mediaRehder。曼地亚红豆杉愈伤组织保存于中国林业科学院研究院林木遗传育种国家重点实验室，利用B5培养基进行培养，28 d继代1次，培养条件为25 ℃，暗培养。转录组测序在安诺优达基因科技（北京）有限公司进行。

2 方法

2.1TmNPR1家族成员鉴定及筛选

从TAIR（https://www.arabidopsis.org/）数据库中获取拟南芥AtNPR1（AT1G64280.1）蛋白序列。将拟南芥AtNPR1序列作为query，与课题组前期测序获得的曼地亚红豆杉全长转录组数据进行Tblastn比对，E值＜1×10−10，获得TmNPR1候选基因核苷酸序列。将所得到的核苷酸序列通过在线网站ORF Finder（https://www.ncbi.nlm.nih. gov/orffinder/）来搜寻序列对应的完整开放阅读框（open reading frame，ORF），再通过CD-HIT软件对NPR1候选基因的ORF序列进行去冗余处理，调整cut-off参数为0.90，最终获得TmNPR1s基因核苷酸序列。

2.2TmNPR1家族系统进化关系分析

利用DNAMAN将所获得的TmNPR1s氨基酸序列、拟南芥NPR1氨基酸序列与中国红豆杉NPR1氨基酸序列[17]进行保守结构域的同源性比对。使用MEGA X软件采用临接法（neighbor-joining，NJ）对所获得的TmNPR1s蛋白序列和其他植物的NPR1蛋白序列（下载自Genbank数据库与TAIR数据库）进行进化树的构建，bootstrap为1 000次重复，其他参数为默认值。利用iTOL（https:// itol.embl.de/）在线网站对进化树进行美化。

2.3TmNPR1家族蛋白理化性质分析

使用ExPaSy-ProtParam tool（https://web. expasy.org/protparam/）网站对所获得的NPR1蛋白序列的氨基酸数量（number of amino acids）、相对分子质量（molecular weight）、等电点（pI）、不稳定系数（instability index）、脂肪系数（aliphatic index）和总平均亲水性（grand average of hydropathicity，GRAVY）等理化性质进行分析。

2.4TmNPR1家族保守结构域与保守基序分析

利用NCBI网站上的在线工具Batch CD-Search（https://www.ncbi.nlm.nih.gov/Structure/ cdd/wrpsb.cgi）对TmNPR1s基因的ORFs进行保守结构域预测分析，随后将预测结果通过DNAMAN和TBtoolsv1.108软件进行可视化处理[18]。利用MEME在线网站（https://meme-suite.org/meme/tools/meme）对TmNPR1s家族蛋白进行保守基序（Motif）分析[19]，最大motif检索数为10，其余为默认参数，将结果通过TBtoolsv1.108软件进行可视化处理[17]。

2.5TmNPR1s基因表达模式分析

从NCBI网站SRA数据库中下载云南红豆杉不同组织（根、茎、叶）的转录组数据（SRP127697）、低温处理（0℃处理2～3年生中国红豆杉幼苗12 h）后红豆杉的转录组数据（SRP096539）和MeJA处理（红豆杉愈伤组织分别经过MeJA处理0、0.5、3、24 h）后的曼地亚红豆杉转录组数据（SRP133888），以及本课题组前期测得的冠菌素（coronatine，COR）处理（在曼地亚红豆杉的B5培养基中加入1 μmol/L COR，分别于0 h、4 h、12 h、48 h、4 d和14 d时收集愈伤组织）后的转录组数据[20]。

利用所获得的TmNPR1s序列分别与以上转录组数据进行Blastn同源比对，一致性最高同源序列的FPKM值即为不同处理下TmNPR1s的表达量，并利用GraphPad Prism 8软件作图，进行基因差异表达分析。

3 结果与分析

3.1TmNPR1家族成员鉴定及其编码蛋白理化性质分析

本研究从曼地亚红豆杉全长转录组中鉴定出3个NPR1基因，通过Batch CD-Search在线工具预测分析发现3条NPR1的蛋白序列均含有BTB_POZ结构域与AnK保守域（图1），是典型的NPR1家族蛋白。将鉴定到的3条基因序列分别命名为TmNPR1～TmNPR3，以便后续分析。

对获得的3个TmNPR1基因编码蛋白进行理化性质分析，如表1所示，发现3个NPR1蛋白分别含有600、519和658个氨基酸残基；相对分子质量在65 020～75 170；3个蛋白的等电点分别为5.86、5.56和5.34，全部小于7，说明该家族蛋白全部偏酸性；不稳定系数为42.43～51.79，表明3个蛋白均不稳定；脂肪系数80.56～93.72，都在80以上，说明都为脂溶性蛋白；总平均亲水性分别为−0.201、−0.235和−0.498，均为负值，表明3个蛋白都为亲水蛋白。

3.2 TmNPR1家族系统进化及多序列比对分析

本研究对包括3个TmNPR1蛋白在内的共32个NPR1蛋白进行系统进化树构建。结果如图2所示，这些NPR1家族蛋白共分为3个亚家族，分别为Clade Ⅰ（AtNPR1/2 subfamily）、Clade Ⅱ（AtNPR3/4 subfamily）和Clade Ⅲ（AtNPR5/6 subfamily）。在3个TmNPR1蛋白中，TmNPR1、TmNPR2聚在Clade II亚家族，TmNPR3聚在Clade III亚家族。在拟南芥中，AtNPR5/6被证明并没有参与免疫反应，但能够调控叶片与花的发育[21]，TmNPR3与AtNPR5/6聚类在一起，说明TmNPR3可能对红豆杉植株生长形态起到调控作用。而TmNPR1/2则与Clade II类AtNPR3/4蛋白亲缘关系较近，表明TmNPR1/2可能具有与AtNPR3/4蛋白相似的功能。

进一步对3个TmNPR1蛋白（TmNPR1-3）序列以及6个拟南芥NPR1蛋白序列（AtNPR1-6）、1个中国红豆杉NPR1蛋白[17]（TcNPR1）进行多序列对比分析，结果如图3所示，所有序列均包含BTB_POZ结构域和锚蛋白重复结构域AnK，这2个结构域在红豆杉和拟南芥中较为保守，但TmNPR3的BTB_POZ结构域和AnK保守域位置与其他序列相比有较大差异（图1）；除AtNPR5/6和TmNPR3外，其余NPR1蛋白序列均包含NPR_like_C区域。此外，在寡聚物-单体相互转变过程中发挥重要作用的3个半胱氨酸残基（Cys82、Cys155和Cys216）在除TmNPR3外的9个NPR1蛋白中完全保守。同时，TmNPR1、TcNPR1和AtNPR3/4还具有EAR基序（VDLNETP），该基序在转录抑制方面发挥着重要的作用。

3.3TmNPR1家族蛋白保守基序分析

利用MEME在线网站对3个TmNPR1蛋白以及拟南芥、中国红豆杉的NPR1蛋白进行保守基序分析，结果如图4所示。所有NPR1蛋白均含有motif1、3、5，但TmNPR3只含有这3种motif，而且motif1、3在TmNPR3蛋白上的位置也与其他NPR1蛋白存在较大差异。此外，除AtNPR5、6与TmNPR3外，剩余7个NPR1蛋白均含有motif1～与TmNPR2具有极高的相似性，除去motif3、7、9在位置上有所偏移，其余motif类型、数量以及位置几乎没有差别，TmNPR2与TmNPR1、TcNPR1相比，也只是motif位置有所不同，可见，在不同物种中NPR1蛋白具有较高的保守性，但也存在motif类型、数量以及位置等差异。

3.4TmNPR1s基因不同组织表达模式分析

为进一步探讨3个TmNPR1基因的表达特性，本研究同时利用转录组数据对红豆杉3个组织（根、茎、叶）中TmNPR1～TmNPR3的表达模式进行分析。结果发现3个基因具有组织表达特异性（图5）。在根中，TmNPR1的表达量最高；在茎中，TmNPR1表达量略高于其他2个基因，但是明显低于该基因在根和叶中表达量；在针叶中，TmNPR1和TmNPR2表达量最高且两者非常相近，TmNPR3在针叶中表达量虽不及前2个基因，但是比其在根和茎中的表达量相对较高。

3.5TmNPR1s基因在低温胁迫表达分析

本研究对低温处理下红豆杉植株中TmNPR1-3的表达水平也进行了检测。结果如图6所示，3个基因在正常温度下表达水平较为接近，但是在低温诱导（LT）后，TmNPR1基因的表达与对照（NT）相比显著提升，而TmNPR2和TmNPR3基因在低温处理后的表达发生了降低，且远低于TmNPR1的表达水平。

3.6TmNPR1s基因在COR和MeJA处理下表达分析

如图7所示，经过COR处理后，TmNPR1表达量明显高于TmNPR2和TmNPR3，并且呈现随着处理时间的增加而升高进而再降低的趋势，在处理4 h时达到最大值，在处理4 d时表达量达到最低值，处理14 d时表达量基本同对照持平。TmNPR2和TmNPR32个基因在COR处理后表达量变化幅度较小，这表明TmNPR1对COR的处理响应更为积极。

本研究还对MeJA处理下曼地亚红豆杉细胞系中TmNPR1s基因的表达水平进行了分析。结果如图8所示，整体来看与COR处理下结果相似，三者中，TmNPR1基因表达水平最高，TmNPR2和TmNPR3表达量较低，且2个基因在处理后各个时间点表达量并没有显著变化。与COR处理不同的是，在MeJA处理后，TmNPR1基因在处理0～24 h呈现逐渐升高趋势，并在24 h表达量最高。

4 讨论

水杨酸是植物抗病反应中激活植物防御保护机制的信号分子[22]，而NPR1是水杨酸信号通路的重要调控因子，在植物中常常以家族形式存在，该基因家族在植物的抗病以及逆境胁迫响应等方面起着重要作用。本研究从曼地亚红豆杉全长转录组中筛选出3个与拟南芥AtNPR1同源性较高的NPR1基因，并发现其有2个进化分枝，TmNPR1/2与AtNPR3/4聚为一类，TmNPR3与AtNPR5/6为一类，该实验结果和鳄梨Persea americana的研究结果类似[23]，可以看出该家族基因在植物的进化种具有高度保守性。

拟南芥中NPR1家族成员，结构十分相似，但是功能却相差甚远。尽管AtNPR1和AtNPR3/AtNPR4都是水杨酸受体，但是它们在水杨酸诱导的防御基因表达的转录调控中起相反的作用。SA通过抑制AtNPR4的转录抑制活性和促进AtNPR1的转录激活活性，从而激活关键免疫调节因子的表达[15]。本研究通过结构域分析，发现3个TmNPR1蛋白的N端均含有BTB_POZ结构域和AnK保守域，但是TmNPR3的C端缺少NPR1_Like_C结构域。Motif分析结果也看出TmNPR3与TmNPR1和TmNPR2结构上存在显著差异。此外，motif分析结果显示，TmNPR1和TmNPR2与拟南芥的AtNPR1-3不论从motif数目还是位置上都一致，但进化树分析结果显示TmNPR1/2与AtNPR3聚在一起。Ding等[15]发现AtNPR3和AtNPR4蛋白含有EAR基序（VDLNETP），该基序在转录抑制方面发挥着重要的作用。多序列比对结果可以看出，TmNPR1/2也具有EAR基序。由此推测TmNPR1/2更有可能与AtNPR3的功能相似。上述结果也表明TmNPR1家族既存在较高的保守性，同时，不同成员之间也很有可能存在功能特异性。

NPR1家族成员还参与器官发育进程，研究发现拟南芥中的AtNPR5和AtNPR6基因，具有调控叶片及花器官生长发育的功能[21, 24]。在本研究组织表达量分析中，TmNPR1家族基因成员TmNPR1/2在针叶和/或根中高表达，TmNPR3在针叶中相对表达量较高。该结果与棉花、小麦NPR1基因表达结果相似，具有典型的组织表达特异性[25-26]。此外，进化树分析结果显示TmNPR3与AtNPR5/6聚在一起，因此，TmNPR3极可能在曼地亚红豆杉中调控针叶的形成，同时TmNPR1/2也可能参与了根和针叶器官的发育。

NPR1家族除了具有抗病以及影响器官发育的功能，还被证实与植物非生物胁迫相关。如在低温胁迫下，AtNPR1与热休克转录因子1（recombinant heat shock tranion factor 1，HSFAI）相互作用，提高拟南芥的抗寒性[27]。在高粱中，SbNPR1基因在PEG6000和NaCl处理0.5 h时表达达到高峰[28]，可见其可能也与干旱和高盐胁迫有关。因此，本研究利用转录组数据对曼地亚红豆杉NPR1基因家族的表达模式进行分析，冷胁迫、COR以及MeJA处理的结果均显示TmNPR1表达量最高，TmNPR2在COR处理14 d表达量较高。可见，三者中，TmNPR1响应激素信号以及低温诱导，而TmNPR3极有可能不参与低温胁迫。此外，本研究中来自曼地亚红豆杉的TmNPR1与前人报道的中国红豆杉的TcNPR1序列一致性达到88%，进化树分析也聚类在一起，推测两者为不同种红豆杉中的直系同源基因。研究发现TcNPR1强烈响应干旱和高盐逆境胁迫[17]，本研究中TmNPR1也能响应低温胁迫，这表明红豆杉TmNPR1基因能响应多种非生物胁迫。

本研究从曼地亚红豆杉转录组中鉴定出3个NPR1基因成员，发现TmNPR1家族基因编码的3个蛋白在N端均具有典型的BTB_POZ结构域和AnK保守域，TmNPR1家族有很大可能与其他植物的NPR1蛋白功能相似。蛋白功能结构域分析结果预示着TmNPR1家族成员之间可能存在功能的特异性。表达模式分析显示，3个TmNPR1基因在红豆杉不同组织中和不同处理下的表达量各不相同，TmNPR1基因可能在红豆杉体内响应多种非生物胁迫，TmNPR3基因则极可能参与红豆杉针叶的形成。该研究为进一步探索红豆杉NPR1基因功能提供理论支撑，并为后续深入挖掘水杨酸调控红豆杉紫杉醇合成的分子机制提供理论依据。

利益冲突所有作者均声明不存在利益冲突

参考文献（略）

来源：杜雨晴，袁菊红，王旭，张恺恺，陈段芬，邱德有，杨艳芳．曼地亚红豆杉NPR1基因家族鉴定及表达模式分析 [J]. 中草药, 2024, 55(1): 240-248.