大豆半胱氨酸蛋白酶抑制剂基因GmCYS2的功能鉴定

柯丹霞 彭昆鹏 贾 妍 曾 硕 王英枝 张静怡



大豆半胱氨酸蛋白酶抑制剂基因的功能鉴定

柯丹霞*彭昆鹏 贾 妍 曾 硕 王英枝 张静怡

信阳师范学院生命科学学院 / 大别山农业生物资源保护与利用研究院, 河南信阳 464000

半胱氨酸蛋白酶抑制剂(cystatin, CYS)在结瘤, 根瘤发育和衰老过程中起重要作用。本研究克隆了1个大豆CYS家族基因, 氨基酸序列比对及进化树分析表明, GmCYS2与木豆() CYS相似性最高。在体外对该基因编码蛋白进行了表达和纯化, 重组蛋白GmCYS2的酶活性抑制实验显示, 该蛋白对L类和B类组织蛋白酶的抑制活性明显高于对H类组织蛋白酶, 并且在30 d根瘤提取物中的抑制活性高于在60 d根瘤提取物中。此外, 构建的植物表达载体, 通过百脉根毛根转化法获得转基因的超表达复合体植株。的超量表达增加了百脉根的结瘤数目, 并且上调共生标记基因的表达。结果说明GmCYS2蛋白具有一定的酶活性抑制作用, 并且正调控百脉根共生结瘤过程。

大豆; 半胱氨酸蛋白酶抑制剂; 半胱氨酸蛋白酶; 共生结瘤; 百脉根

半胱氨酸蛋白酶抑制剂CYS是一类小分子蛋白质, 能够可逆地抑制木瓜蛋白酶C1A和豆荚蛋白C13肽酶家族中的半胱氨酸蛋白酶活性[1]。它们的保守结构域包含一个甘氨酸残基和两个发夹环, 一个保守的QxVxG结构域和一个色氨酸残基[2]。此外, 在植物CYS中还普遍存在着另一个保守的LARFAV结构域([LVI]-[AGT]-[RKE]-[FY]-[AS]-[VI])[3]。CYS和它们的半胱氨酸蛋白酶之间的动态平衡在调控植物发育、植物响应生物胁迫以及非生物胁迫等过程中发挥重要作用, 如内源性蛋白质转换[4], 贮存蛋白的积累[5], 种子的萌发和成熟[6], 细胞程序性死亡[7], 非生物环境胁迫[8]以及植物防御螨[9]、真菌[10]的侵害等。

Tan等[11]利用已公布的苹果()基因组数据库, 系统鉴定了苹果基因组中所有的基因, 并详细分析了其全基因组特征, 检测了该类基因在干旱、高温、盐、脱落酸(ABA)等不同非生物胁迫下的表达模式。对单个基因的功能研究发现基因可能通过影响植株体内ROS的积累和根毛的发育从而增强转基因拟南芥的耐旱性[12]。基因参与ABA胁迫反应, 通过促进气孔关闭上调ABA和干旱相关基因的转录水平, 从而增强转基因拟南芥和苹果的抗旱能力, 能够与类受体蛋白激酶MpFER相互作用共同参与苹果干旱胁迫应答过程[13]。基因还能够延缓自然和应激条件下苹果叶片的衰老[14]。在拟南芥中的异位表达可以同时增强转基因拟南芥对盐胁迫和对衣霉素诱导的内质网胁迫的抗性[15]。最新的研究发现, 在拟南芥种子萌发过程中,基因的表达受热应激和外源ABA胁迫诱导, 过量表达该基因增强了拟南芥的耐高温能力和对外源ABA的不敏感性[16]。水稻基因沉默后, 转基因植株表现出较高的豆荚蛋白酶和木瓜蛋白酶活性, 幼苗生长受到促进, 碱胁迫下RNAi植株根的形成和茎的生长受到抑制, 说明基因具备双重功能, 同时参与植物发育过程和胁迫反应[17]。野生大豆GsCPI14通过与钙/钙调素结合的受体类激酶GsCBRLK相互作用增强植物的耐碱性[18]。此外, 拟南芥、和基因受到甜菜胞囊线虫感染后表达量上升[19]。苎麻BnCPI的体外重组蛋白可以有效抑制木瓜蛋白酶和无花果蛋白酶的蛋白酶活性, 还能抑制一些重要植物病原真菌的菌丝生长[20]。

综上,基因已经在许多植物如拟南芥、水稻、苹果、苎麻、大豆中得到了分离和鉴定, 其在胁迫反应中的作用研究得最细致和深入, 但关于该类基因在豆科植物结瘤固氮过程中的作用却知之甚少。随着大豆全基因组测序的完成和公布[21], 目前已经有20个大豆CYS家族基因被鉴定出来, RNA- seq数据显示, 这些基因在大豆的14个不同组织中均有表达[22], 其中7个基因在根瘤中具有转录活性[23]。Yuan等[24]对20个大豆CYS家族基因的序列结构特征、启动子调控元件、亲缘进化关系及组织表达特征等进行了系统详细的分析, 并检测了该类基因在大豆接种根瘤菌前后不同时期的根及根瘤中的表达情况, 以为候选基因, 利用发根农杆菌介导的百脉根毛根转化法, 深入研究了该基因在结瘤过程中的生物学功能。同时该研究指出,基因在接种后5 d和16 d根中表达水平明显升高, 在30 d根瘤中的表达量迅速上升, 为12 d根瘤的6倍, 在42 d和84 d根瘤中基因也能维持较高转录水平, 因此推测基因可能在根瘤的发生、发育和衰老过程中发挥功能。鉴于此, 本研究以为对象, 拟通过蛋白体外抑制活性的测定揭示其生化功能, 通过转基因的超表达百脉根植株共生结瘤表型的鉴定揭示其生物学功能, 为深入研究该类基因在共生信号转导途径中发挥的作用提供新的分子生物学证据, 为进一步寻找根瘤特异的半胱氨酸蛋白酶提供有益的线索。

1 材料与方法

1.1 试验材料

大豆Willimas 82 (W82)种子由中国科学院东北地理与农业生态研究所孔凡江研究员提供,由复旦大学生命科学学院王应祥教授提供百脉根() MG-20种子、百脉根中生根瘤菌菌株、植物表达载体pU1301由华中农业大学农业微生物国家重点实验室张忠明教授提供。

1.2 植物材料处理

大豆W82种子经表面灭菌后, 置无菌水润湿的滤纸上萌发, 于16 h/8 h光照/黑暗, 相对湿度70%, 28℃~无菌Fahraeus同样条件~–℃

1.3 RNA提取、RT-PCR及qPCR检测

参照RNA提取试剂盒(Invitrogen, USA)说明书, 提取大豆的总RNA。使用NanoDrop仪器检测RNA样品的浓度和纯度。按照TaKaRa公司反转录试剂盒操作说明得到cDNA第1链, RT-PCR扩增目的基因。按照TaKaRa公司PrimeScript RT reagent Kit操作说明进行荧光定量PCR检测, 根据相对定量法(2–ΔΔCt)公式计算结果。

1.4 序列分析

利用NCBI网站的Blastp工具寻找GmCYS2蛋白的同源蛋白, 用DNAMAN软件进行序列同源性比对和进化树分析。

1.5 重组蛋白GmCYS2的表达和纯化

根据NCBI网站公布的基因GenBank登录号为XP_003524913.1, 使用表1中F-GmCYS2- pro.和R-GmCYS2-pro.引物, 通过PCR扩增除去基因N端信号肽序列(1~78 bp)之外的cDNA片段, 并插入大肠杆菌表达载体pGS-21a。将测序正确的重组质粒pGS-21a-GmCYS2转化到大肠杆菌BL21 (DE3)中。将单个菌落接种到含有氨苄霉素的培养基中, 以37℃摇床培养。当细胞密度在600 nm达到OD = 0.6时, 加入IPTG诱导表达。诱导条件分别为15℃诱导16 h, 37℃诱导4 h。用SDS-PAGE和Western blot检测不同诱导条件下蛋白的表达情况。最终从37℃诱导4 h的细胞裂解液上清液中, 通过一步Ni柱纯化法得到GmCYS2蛋白。以BSA为标准品, 用Bradford法测定蛋白浓度, 通过考马斯亮蓝染色的SDS-PAGE凝胶扫描分析估计蛋白纯度。蛋白经0.22 μm滤膜过滤除菌, 为避免重复解冻, 同等分分装保存在–80℃冰箱备用, 存储缓冲液为50 mmol L–1Tris-HCl, 150 mmol L–1NaCl, 10% Glycerol, pH 8.0。

1.6 GmCYS2蛋白的体外酶活性抑制实验

取出–80℃保存的30 d和60 d大豆根瘤样品, 在研钵中充分研磨, 加入0.15 mol L–1的NaCl重悬浮, 在50 mmol L–1pH 6.0的磷酸钠与2 mmol L–1EDTA缓冲液中, 4℃孵育1 h, 离心15 min, 收集上清液。采用Protein Assay试剂盒(Bio-Rad公司)定量分析蛋白质含量。通过根瘤提取物对特异性荧光底物的水解作用来测定重组蛋白GmCYS2的抑制活力。荧光标记的底物被酶解, AMC (Ex/Em = 365 nm/ 465 nm)荧光基团得以游离, 产生亮蓝色荧光, 因此很容易被DAPI滤光片所检测。反应过程如下: 在试管中加入1 μmol L–1GmCYS2纯化蛋白和1 μg根瘤蛋白提取液, 25℃反应10 min, 随后分别加入蛋白酶特异底物Z-FR-AMC (组织蛋白酶L类)、Z-RR-AMC (组织蛋白酶B类)和Bz-FVR-AMC (组织蛋白酶H类), 终浓度为25 mmol L–1, 反应缓冲液为pH 8.0的HEPES, 反应体系为500 μL[4]。混合液在30℃孵育1 h后, 以同体积冰乙醇终止反应, 加去离子水3 mL。最后进行荧光测定, 以A值反映酶活性。所有的分析设置3组重复, 以不加抑制剂蛋白的试管作为空白对照。抑制率(%) = (A空白对照–A加抑制剂)/ A空白对照×100%。

表1 本研究中所使用的引物

小写字母区域为酶切位点序列。

The lowercases in sequences indicate the restriction enzyme cutting sites.

1.7 转GmCYS2基因的超表达百脉根植株的获得

根据NCBI网站公布的基因序列, 使用引物F-OX和R-OX (表1), 扩增基因全长序列并构建植物表达载体pU1301-GmCYS2, 重组质粒经冻融法转入发根农杆菌LBA1334菌株中备用。百脉根()的毛根转化方法如前所述[25-26]。简言之, 收集剪去下胚轴的幼苗子叶部分, 放入准备好的携带有重组质粒的农杆菌LBA1334菌悬液中, 侵染30 min, 并转移到Murashige和Skoog (MS)平板上共培养5 d。随后移至含250 mg L–1羧卞的MS培养基上, 直至毛根长出。标记每根毛根并用GUS染液鉴定。每个幼苗只留下1~2个转基因阳性毛根。提取阳性毛根总RNA, 反转录后以表1中F-OX和R-OX为引物扩增基因, 百脉根多聚泛素(,)基因作为内参, 携带有空载体(pU1301)阳性毛根的复合体植株作为阴性对照。将RT-PCR检测为阳性的复合体植株移入无菌沙盆中炼苗1周, 接种根瘤菌, 每天浇灌无菌Fahraeus无氮营养液。

1.8 复合体植株的共生表型鉴定及基因表达检测

接种30 d后, 对结瘤表型进行照相统计, 记录单株结瘤数和样本量()并计算每个植株(CK和GmCYS2-OX)的平均结瘤数, 重复2次。利用F-GmCYS2-rt和R-GmCYS2-rt引物, 荧光定量PCR检测超表达植株阳性毛根中的转录水平, 利用引物F-NIN-rt和R-NIN-rt, F-ENOD40-1-rt和R-ENOD40-1-rt, F-ENOD40-2-rt和R-ENOD40-2-rt分别检测结瘤标记基因、和的转录水平, 百脉根基因作为内参。重复试验3次。采用Microsoft Excel 2007和SPSS 13.0软件分析数据并绘制图表。

2 结果与分析

2.1 GmCYS2基因序列同源比对及进化树分析

序列分析表明, 完整的基因全长为393 bp, 克隆的基因序列与参考基因组序列完全一致。该基因编码130个氨基酸残基, N端包含26个氨基酸残基编码的信号肽。在NCBI网站上对基因的编码序列进行Blastp比对发现, 该序列与木豆(XP_020204069.1)、狭叶羽扇豆(XP_019445765.1)、绿豆(XP_014508215.1)、赤豆(XP_017435098.1)、鹰嘴豆(XP_004504436.1)、蒺藜苜蓿(XP_013446602.1)、落花生(XP_016196387.1)、蔓花生(XP_015958299.1)、水稻(XP_015635666.1)和烟草(XP_016482039.1)的CYS相似性分别为73%、67%、64%、60%、55%、54%、54%、54%、50%和47%。GmCYS2与它的同源蛋白都包含3个保守的结构域, 即N端的“G” 残基, 序列中央的“QxVxG”, C端的“W”残基, 另外还有1个可变的“LARFAV”序列(图1-A)。通过DNAMAN软件获得大豆GmCYS2蛋白与上述10种物种CYS蛋白的进化树(图 1-B), 大豆GmCYS2蛋白与木豆CYS蛋白处在同一进化分支上, 亲缘关系最近。

2.2 重组蛋白GmCYS2的表达和纯化

由于蛋白酶抑制剂在中可溶性较差, 经包涵体复性纯化的蛋白酶活性很难保证, 需添加大标签GST用于蛋白表达。此外, GmCYS2序列的N端第1~26位氨基酸编码信号肽, 在构建体外表达载体时需将此信号肽去掉。因此我们构建的重组蛋白序列如图2所示, 包含双标签His-GST以及GmCYS2 C端的104个氨基酸, 蛋白大小为40 kDa左右(图2)。蛋白诱导条件分别为15℃诱导16 h和37℃诱导4 h。用SDS-PAGE和Western blot方法检测不同诱导条件下重组蛋白在中的表达情况。在2种诱导条件下的全细胞裂解液中都检测到了蛋白的表达(图3-A、B中1和2泳道)。全细胞裂解液经超声波破碎、离心, 分别取上清液和沉淀进行SDS-PAGE分析显示, 15℃诱导16 h条件下, 目的蛋白条带主要出现在沉淀中, 上清液中几乎检测不到(图3-A、B中3和4泳道)。37℃诱导4 h条件下, 目的蛋白条带在上清液和沉淀中都存在(图3-A, B中5和6泳道)。因此, 重组蛋白最终从37℃诱导4 h的细胞裂解液上清液中, 通过一步Ni 柱纯化法获得(图3-C)。以BSA为标准品, 用Bradford法测定, 纯化蛋白的浓度为1.80 mg mL–1, 通过考马斯亮蓝染色的SDS-PAGE凝胶扫描分析估计蛋白纯度约为85%。

2.3 GmCYS2蛋白的体外酶活检测

在30 d和60 d根瘤蛋白提取液中, GmCYS2重组蛋白主要抑制L类和B类组织蛋白酶的活性, 抑制率在26.7%~39.1%之间, 而对H类组织蛋白酶的抑制活性较低(< 15%)。此外, GmCYS2蛋白在30 d根瘤提取物中的抑制活性明显高于在60 d根瘤提取物中(< 0.05)(图4)。说明GmCYS2蛋白在体外具有蛋白酶抑制活性, 且具有种类特异性和组织特异性。

2.4 过量表达GmCYS2对百脉根共生结瘤的影响

将基因在玉米()泛素启动子(ZmUb)的控制下表达于百脉根的转基因毛状根中, 分别获得pU1301-GmCYS2 (GmCYS2-OX)和空载体pU1301 (CK)复合体植株, 根瘤菌接种4周后记录结瘤表型。结果GmCYS2-OX毛状根长出更多的根瘤(图5-A), 每个植株的结瘤数从对照的5.8个增加到10.1个(图5-B和表2)。对单株结瘤数量分布的分析显示, 45%的GmCYS2-OX植株结瘤数目>10。与之相反, 对照中没有观察到结瘤数目>10的植株(图5-C)。转基因阳性毛根经GUS染色后呈现明显的蓝色, 阴性毛根仍为白色(图5-D), 进一步通过RT-PCR检测阳性毛根中基因的表达情况, 发现空载体对照中没有出现基因的转录信号, 而随机挑选的3个转基因复合体植株中均出现明显的转录信号(图5-E), 说明植物表达载体的构建是有效的。

2.5 转基因毛根中GmCYS2以及共生基因的转录水平检测

通过实时荧光定量PCR检测显示, 在超表达植株中mRNA的表达水平为对照的8.8倍(图6), 表明基因的过量表达可以引起植株结瘤数目的增加。与此同时, 以qPCR分析、和这3个与结瘤相关的基因在超表达植株中的表达水平, 进一步验证基因的过量表达对结瘤相关Marker基因的影响。结果表明这3种结瘤基因的表达水平, 特别是和基因, 与对照相比, 在GmCYS2-OX毛状根中显著增加(图6)。

图1 大豆GmCYS2与其他植物中同源蛋白的序列比对(A)及进化树分析(B)

图A中用黄色、红色、绿色和蓝色方框标出半胱氨酸蛋白酶抑制剂的主要保守结构域。图B中标尺代表遗传相似性, 指不同植物间同源蛋白进化关系的远近。

The main cystatin conserved motifs are in yellow, red, green, and blue boxes in panel A. The scale in panel B represents genetic similarity, indicating the proximity relationships among species.

图2 His-GST-GmCYS2融合蛋白的序列分析

氨基酸第5~10位: His标签; 第14~233位: GST标签; 第244~250位: TEV(烟草蚀纹病毒)切割位点; 251~355位: GmCYS2蛋白序列(去掉N端26个氨基酸编码的信号肽)。

5–10: His tag; 14–233: GST tag; 244–250: TEV cleavage site; 251–355: GmCYS2 protein sequence (Removing signal peptide of 26 amino acids at N end).

图3 GmCYS2蛋白的体外表达、抗体检测和蛋白纯化

A: SDS-PAGE分析蛋白在BL21(DE3)菌株中的表达情况; B: Western blot分析, 抗-His抗体检测; C: 纯化蛋白的SDS-PAGE分析。M1: 蛋白marker; M2: Western marker; M3: 蛋白预染marker; PC1: BSA (1 μg); PC2: BSA (2 μg); NC: 未诱导的全细胞裂解液。1: 15℃诱导16 h的全细胞裂解液; 2: 37℃诱导4 h的全细胞裂解液; NC1: 未诱导的细胞裂解液上清液; NC2: 未诱导的细胞裂解液沉淀; 3: 15℃诱导16 h的细胞裂解液上清液; 4: 15℃诱导16 h的细胞裂解液沉淀; 5: 37℃诱导4 h的细胞裂解液上清液; 6: 37℃诱导4 h的细胞裂解液沉淀; 7: 从37℃诱导4 h的细胞裂解液上清液中纯化的蛋白。

A: SDS-PAGE analysis of protein expression in BL21 (DE3) strain; B: Western blot analysis of protein with anti-His antibody; C: SDS-PAGE analysis of purified protein. M1: protein marker; M2: Western marker; M3: protein pre staining marker; PC1: BSA (1 μg); PC2: BSA (2 μg); NC: Non induced whole cell lysate. 1: Whole cell lysate induced under 15°C for 16 h; 2: Whole cell lysate induced under 37°C for 4 h; NC1: Supernatant of non-induced cell lysate; NC2: The precipitation of non-induced cell lysate; 3: Supernatant of cell lysate induced under 15°C for 16 h; 4: The precipitation of cell lysate induced under 15°C for 16 h; 5: Supernatant of cell lysate induced under 37°C for 4 h; 6: The precipitation of cell lysate induced under 37°C for 4 h; 7: Purified protein from supernatant of cell lysate induced under 37°C for 4 h.

表2 过表达GmCYS2复合体百脉根植株的结瘤数目统计

统计接种30 d后单个植株的平均结瘤数目, 括号中的代表植株总数, 试验重复2次(1和2)。

Mean number of nodules per plant at 30 d after rhizobial inoculation with SD are shown. The number () of scored plants is indicated in parentheses. Two independent experiments (1 and 2) were performed.

图4 重组蛋白GmCYS2在不同时期根瘤提取物中的抑制活性分析

样品分别为30 d和60 d根瘤提取物。反应液中加入纯化的GmCYS2蛋白, 另外分别加入组织蛋白酶L类、B类和H类3种不同底物, 试验重复3次。*< 0.05。

Samples were extracts of 30-day- and 60-day-old nodules. Reaction buffer added with purified GmCYS2 protein and the specific substrates cathepsin L-like, B-like, and H-like, respectively. Values are triplicate measurements. *< 0.05.

3 讨论

固氮效率与结瘤, 根瘤发育和衰老密切相关。

已有证据表明, 半胱氨酸蛋白酶在几种豆科植物根瘤发育和衰老过程中扮演重要角色, 如基因沉默紫云英推迟根瘤发育和类菌体衰老, 延长了紫云英根瘤寿命[27]。抑制半胱氨酸蛋白酶CYP15A延缓蒺藜苜蓿根瘤的衰老[28]。豌豆在不定型根瘤开始衰老时被激活[29]。在大豆中, 1983年最早有报道半胱氨酸蛋白酶在根瘤衰老时表达, 18个大豆半胱氨酸蛋白酶被发现在根瘤发育和衰老时有很强的转录活性[30]。然而, 这些半胱氨酸蛋白酶在根瘤发育和衰老过程中到底发挥怎样的功能？为什么如此多的半胱氨酸蛋白酶参与调控根瘤共生？这些问题目前并没有得到很好的阐明。先前的研究表明抑制半胱氨酸蛋白酶能够延迟根瘤衰老, 所以对它们的天然抑制剂CYS的研究将有助于解答这些疑问。大豆CYS家族基因的序列结构以及表达特征已经得到了详尽和细致的揭示, 但是对于单个基因具体功能的研究还十分有限, 特别是在共生结瘤信号途径中的功能及机制研究的报道很少。本研究克隆得到的基因与植物中大多数基因一样, 氨基酸序列中包含3个典型基序即N端的“G”残基, 中部的“QVVSG”序列和C端的“W”残基, 还有一个保守的“LGRFSV”序列(图1),这段序列在家族基因的不同分支中是高度变化的。

图5 GmCYS2基因的过量表达对百脉根共生结瘤的影响

A: 复合体植株接种根瘤菌30 d后的结瘤表型。毛状根表达pMUb:GmCYS2 (GmCYS2-OX)或空载体pU1301 (CK), Bar = 10 mm; B:单株复合体百脉根植株的平均结瘤数目; C:根据单株根瘤数划分的不同结瘤种类所占的相对百分比; D:阳性毛根的GUS鉴定; E: RT-PCR检测阳性毛根中基因的表达。M: DL2000 DNA marker; 1: 阳性质粒对照; 2: 空载体对照; 3~5: 复合体植株。

A: Transgenic plants at 30 days after inoculation withshowing hairy roots expressing pMUb: GmCYS2 (GmCYS2-OX) or the empty vector pU1301 (CK). Bar = 10 mm; B: Mean numbers of nodules per plant of; C: The relative ratios of various groups divided by nodule number per plant; D: GUS identification of positive hairy roots; E: Detection ofgene expression in positive hairy roots by RT-PCR; M: DL2000 DNA marker; 1: plasmid as positive control; 2: empty vector control; 3–5: composite plants.

图6 转基因毛根中GmCYS2以及共生基因的转录水平检测

荧光定量PCR检测以及共生基因和在超表达植株(GmCYS2-OX)和对照植株(CK)阳性毛根中的转录水平。*<0.05, **<0.01。

qPCR analysis of the transcript levels ofand symbiotic marker genes,, andin the GmCYS2-OX and control hairy roots. *<0.05, **<0.01.

CYS通过抑制半胱氨酸蛋白酶的活性, 在植物种子萌发、宿主免疫、高敏感细胞死亡、病原体及昆虫的防御反应等过程中发挥关键作用。本研究在体外表达和纯化了GmCYS2重组蛋白(图3), 酶活性抑制实验显示, GmCYS2蛋白对L类组织蛋白酶的抑制活性最强, 其次是B类组织蛋白酶。此外, GmCYS2蛋白在成熟根瘤中的抑制活性明显高于衰老根瘤(图4)。该实验不仅证实GmCYS2蛋白在体外具有蛋白酶抑制活性, 还进一步揭示该蛋白对抑制的半胱氨酸蛋白酶种类具有选择性, 更有可能在成熟根瘤而不是衰老根瘤中发挥功能。

基因的表达也受非生物胁迫和激素信号的影响[31]。启动子预测发现,启动子具有脱落酸、赤霉素、水杨酸、生长素和茉莉酸5种激素反应元件。有趣的是,启动子有一个结瘤特异因子结合位点元件, 该元件也存在于典型结瘤素基因大豆血红蛋白LBC3的启动子片段中。在接种根瘤菌的大豆根中表达量明显上升[24], 在根瘤衰老过程中表达量也是上升的[23]。这些研究结果表明,很可能参与大豆共生结瘤过程。水稻的1个基因在大豆中异位表达, 正调控大豆的结瘤和缺氮反应[32]。本研究利用农杆菌介导的百脉根毛根转化法, 对转基因的超表达百脉根植株共生表型进行记录及统计分析表明, 大豆基因在百脉根中的异位表达能够显著增加百脉根的结瘤数目(图5), 并且上调共生Marker基因的表达(图6)。结瘤起始基因是侵入线和根瘤原基形成必需的关键基因[33]。结瘤素基因也是根瘤起始和发育至关重要的Marker基因, 百脉根中有2个基因,和[34]。由此推测GmCYS2可能通过抑制特异的半胱氨酸蛋白酶活性, 进而在结瘤过程中发挥正调控作用。

4 结论

克隆了大豆半胱氨酸蛋白酶抑制剂基因, 该基因编码蛋白具有典型的CYS家族蛋白保守结构域。GmCYS2蛋白在体外具有蛋白酶抑制活性, 过量表达该基因能显著增加转基因百脉根的结瘤数目, 并且上调共生关键Marker基因的表达。推测GmCYS2可能通过抑制特异的半胱氨酸蛋白酶活性, 在结瘤过程中发挥正调控作用。该结果为揭示半胱氨酸蛋白酶抑制剂在共生结瘤过程中的功能及其调控机制提供了分子证据。

[1] Martinez M, Diaz-Mendoza M, Carrillo L, Diaz I. Carboxy terminal extended phytocystatins are bifunctional inhibitors of papain and legumain cysteine proteinases., 2007, 581: 2914–2918

[2] Stubbs M T, Laber B, Bode W, Huber R, Jerala R, Lenarcic B, Turk V. The refined 2.4 A X-ray crystal structure of recombinant human stefin B in complex with the cysteine proteinase papain: a novel type of proteinase inhibitor interaction., 1990, 9: 1939–1947

[3] Margis R, Reis E M, Villeret V. Structural and phylogenetic relationships among plant and animal cystatins., 1998,359: 24–30

[4] Martinez, M, Cambra I, Carrillo L, Diaz-Mendoza M, Diaz I. Characterization of the entire cystatin gene family in barley and their target cathepsin L-like cysteine-proteases, partners in the hordein mobilization during seed germination., 2009,151: 1531–1545

[5] Benchabane M, Schlüter U, Vorster J, Goulet M C, Michaud D. Plant cystatins., 2010, 92: 1657–1666

[6] Arai S, Matsumoto I, Emori Y, Abe K. Plant seed cystatins and their target enzymes of endogenous and exogenous origin., 2002,50: 6612–6617

[7] Belenghi B, Acconcia F, Trovato M, Perazzolli M, Bocedi A, Polticelli F, Ascenzi P, Delledonne M. AtCYS1, a cystatin from, suppresses hypersensitive cell death., 2003,270: 2593–2604

[8] Hwang J E, Hong J K, Lim C J, Chen H, Je J, Yang K A, Kim D Y, Choi Y J, Lee S Y, Lim C O. Distinct expression patterns of two Arabidopsis phytocystatin genes,and, during development and abiotic stresses., 2010,29: 905–915

[9] Carrillo L, Martinez M, Ramessar K, Cambra I, Castañera P, Ortego F, Díaz I. Expression of a barley cystatin gene in maize enhances resistance against phytophagous mites by altering their cysteine-proteases., 2011,30: 101–112

[10] Popovic M, Andjelkovic U, Burazer L, Lindner B, Petersen A, Gavrovic-Jankulovic M. Biochemical and immunological chara­cterization of a recombinantly-produced antifungal cysteine proteinase inhibitor from green kiwifruit ()., 2013,94: 53–59

[11] Tan Y, Wang S, Liang D, Li M, Ma F. Genome-wide identification and expression profiling of the cystatin gene family in apple (×Borkh.)., 2014,79: 88–97

[12] Tan Y X, Li M J, Ma F W. Overexpression of MpCYS2, a phytocystatin gene from(Willd.) Borkh., confers drought tolerance and protects against oxidative stress in Arabidopsis., 2015, 123: 15–27

[13] Tan Y X, Li M J, Yang Y L, Sun X, Wang N, Liang B W, Ma F W. Overexpression of MpCYS4, a phytocystatin gene from(Willd.) Borkh., enhances stomatal closure to confer drought tolerance in transgenic Arabidopsis and apple., 2017, 8: 33

[14] Tan Y X, Yang Y L, Li C, Liang B W, Li M J, Ma F W. Overexpression of MpCYS4, a phytocystatin gene from(Willd.) Borkh., delays natural and stress-induced leaf senescence in apple., 2017, 115: 219–228

[15] Tan Y X, Wei XY, Wang P, Sun X, Li M J, Ma F W. A phytocystatin gene from(Willd.) Borkh.,, confers salt stress tolerance and functions in endoplasmic reticulum stress response in Arabidopsis., 2016, 34: 62–75

[16] Song C, Kim T, Chung W S, Lim C O. The Arabidopsis phytocystatinenhances seed germination and seedling growth under heat stress conditions., 2017, 40: 577–586

[17] Christoff A P, Passaia G, Salvati C, Alves-Ferreira M, Margis- Pinheiro M, Margis R. Rice bifunctional phytocystatin is a dual modulator of legumain and papain-like proteases., 2016, 92: 193–207

[18] Sun X, Yang S, Sun M, Wang S, Ding X, Zhu D, Ji W, Cai H, Zhao C, Wang X, Zhu Y. A novelcysteine proteinase inhibitor GsCPI14, interacting with the calcium/calmodulin- binding receptor-like kinase GsCBRLK, regulated plant tolerance to alkali stress., 2014, 85: 33–48

[19] Labudda M, Różańska E, Szewińska J, Sobczak M, Dzik J M. Protease activity and phytocystatin expression inupon heterodera schachtii infection., 2016, 109: 416–429

[20] Yu Y, Zhang G, Li Z, Cheng Y, Gao C, Zeng L, Chen J, Yan L, Sun X, Guo L, Yan Z. Molecular cloning, recombinant expression and antifungal activity of, a Cystatin in Ramie (L.).(Basel), 2017, 8: 265

[21] Schmutz J, Cannon S B, Schlueter J, Ma J, Mitros T, Nelson W, Hyten D L, Song Q, Thelen J J, Cheng J, Xu D, Hellsten U, May G D, Yu Y, Sakurai T, Umezawa T, Bhattacharyya M K, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, Gill N, Joshi T, Libault M, Sethuraman A, Zhang X C, Shinozaki K, Nguyen H T, Wing R A, Cregan P, Specht J, Grimwood J, Rokhsar D, Stacey G, Shoemaker R C, Jackson S A. Genome sequence of the palaeopolyploid soybean., 2010,463: 178–183

[22] Severin A J, Woody J L, Bolon Y T, Joseph B, Diers B W, Farmer A D, Muehlbauer G J, Nelson R T, Grant D, Specht J E, Graham M A, Cannon S B, May G D, Vance C P, Shoemaker R C. RNA-Seq atlas of: a guide to the soybean transcriptome., 2010,10: 160

[23] Van Wyk S G, Du Plessis M, Cullis C A, Kunert K J, Vorster B J. Cysteine protease and cystatin expression and activity during soybean nodule development and senescence., 2014,14: 294

[24] Yuan S, Li R, Wang L, Chen H, Zhang C, Chen L, Hao Q, Shan Z, Zhang X, Chen S, Yang Z, Qiu D, Zhou X. Search for nodulation and nodule development-related cystatin genes in the genome of soybean ()., 2016, 7: 1595

[25] Ke D X, Li X Y, Han Y P, Cheng L, Yuan H Y, Wang L. ROP6 is involved in root hair deformation induced by Nod factors in., 2016, 108: 488–498

[26] 柯丹霞, 李祥永, 彭昆鹏, 韩雅彭. 百脉根基因启动子的克隆与表达分析. 信阳师范学院学报(自然科学版), 2018, 31: 46–51 Ke D X, Li X Y, Peng K P, Han Y P. Cloning and expression analysis of the promoter region ofgene of.(Nat Sci Edn), 2018, 31: 46–51 (in Chinese with English abstract)

[27] Li Y, Zhou L, Li Y, Chen D, Tan X, Lei L, Zhou J. A nodule- specific plant cysteine proteinase, AsNODF32, is involved in nodule senescence and nitrogen fixation activity of the green manure legume., 2008,180: 185–192

[28] Sheokand S, Dahiya P, Vincent J L, Brewin N J. Modified expression of cysteine protease affects seed germination, vegetative growth and nodule development in transgenic lines of, 2005,169: 966–975

[29] Vincent J L, Brewin N J. Immunolocalization of a cysteine protease in vacuoles, vesicles, and symbiosomes of pea nodule cells., 2000,123: 521–530

[30] Pfeiffer N E, Torres C M, Wagner F W. Proteolytic activity in soybean root nodules: activity in host cell cytosol and bacteroids throughout physiological development and senescence.,1983, 71: 797–802

[31] Zhang X, Liu S, Takano T. Two cysteine proteinase inhibitors from, AtCYSa and AtCYSb, increasing the salt, drought, oxidation and cold tolerance., 2008,68: 131–143

[32] Quain M D, Makgopa M E, Cooper J W, Kunert K J, Foyer C H. Ectopic phytocystatin expression increases nodule numbers and influences the responses of soybean () to nitrogen deficiency., 2015,112: 179–187

[33] 柯丹霞, 李祥永. 结瘤信号途径中相关调控蛋白的研究进展. 信阳师范学院学报(自然科学版), 2015, 28: 621–626Ke D X, Li X Y. Research progress of key regulatory proteins in nodulation pathway.(Nat Sci Edn), 2015, 28: 621–626 (in Chinese with English abstract)

[34] Ke D X, Fang Q, Chen C F, Zhu H, Chen T, Chang X J, Yuan S L, Kang H, Ma L, Hong Z L, Zhang Z M. Small GTPase ROP6 interacts with NFR5 and is involved in nodule formation in., 2012, 159: 131–143

更正: 发表于本刊2018, 44(2): 185–196, 朱莹等“R2R3-MYB转录因子GmMYB184 调节大豆异黄酮合成”一文中的图9和图10更正为:

图9 毛状根中沉默和过表达GmMYB184 的相对表达量

A: 毛状根中沉默4 的相对表达量; B: 毛状根中过表达4 的相对表达量。**和***分别表示转基因株系基因相对表达量与对照在<0.01 和<0.001 水平差异显著。

A: relative expression levels of RNAi ofin hairy roots; B: relative expression levels of over expression ofin hairy roots.

** and *** represent the relative expression levels ofbetween the transgenic lines and control lines were significantly different at the 0.01 and 0.001 probability levels, respectively.

图10 沉默和过表达GmMYB184 毛状根相对异黄酮含量

A: 沉默4 毛状根中相对异黄酮含量; B: 过表达4 毛状根中相对异黄酮含量。**表示转基因株系基因相对表达量与对照在<0.01 水平差异显著。

A: relative isoflavone contents of RNAi ofin hairy roots; B: relative isoflavone contents of over expression ofin hairy roots. ** represents the relative expression levels ofbetween the transgenic lines and control lines were significantly different at the 0.01 probability level.

Functional Characterization of Soybean Cystatins Gene

KE Dan-Xia*, PENG Kun-Peng, JIA Yan, ZENG Shuo, WANG Ying-Zhi, and ZHANG Jing-Yi

College of Life Sciences / Institute for Conservation and Utilization of Agro-bioresources in Dabie Mountains, Xinyang Normal University, Xinyang, Henan 464000, China

CYS (cystatin) plays an important role in nodulation, nodule development and senescence. In this study, we cloned a soybean CYS family gene, which had the highest similarity with pigeonpea () CYS, shown by amino acid sequence alignment and phylogenetic tree analysis. The gene encoding protein was expressed and purified. Enzyme acti­vity inhibition experiment of recombinant protein GmCYS2 showed that the cathepsin inhibitory activity of GmCYS2 protein on L and B classes was higher than that on H class, and cathepsin inhibitory activity in nodule extract from nodule at 30 d was gene­rally higher than that at the 60 d. In addition, the plant overexpression vector was constructed and the-overexpression composite plants were obtained byhairy root transformation method. Overexpression ofincreased the number of nodules inplants and up regulated the expression levels of symbiotic marker genes. These results suggest that GmCYS2 protein can inhibit the enzyme activity and positively regulate the symbiotic nodulation of.

;cystatin; cysteine protease; symbiotic nodulation;

2018-03-26;

2018-04-20.

10.3724/SP.J.1006.2018.01159

柯丹霞, E-mail: kdx_029@163.com

2017-11-23;

本研究由国家自然科学基金项目(31400213), 河南省科技攻关计划项目(182102110448), 信阳师范学院青年骨干教师资助计划项目(2015)和信阳师范学院“南湖学者奖励计划”青年项目(2016)。

This study was supported by the National Natural Science Foundation of China (31400213), Henan Province Science and Technology Research Projects (182102110448), the Funding Scheme for Young Core Teachers of Xinyang Normal University (2015), and Nanhu Scholars Program for Young Scholars of XYNU (2016).

URL: http://kns.cnki.net/kcms/detail/11.1809.S.20180419.1721.008.html