‘84K’杨组氨酸激酶基因PaHK3b 的克隆及功能分析

鲁俊倩，武 舒，钟姗辰，张伟溪，苏晓华，张冰玉

(林木遗传育种国家重点实验室，国家林业和草原局林木培育重点实验室，中国林业科学研究院林业研究所，北京 100091)

细胞分裂素是植物生长过程中的一类重要调控激素，其主要通过双组分系统进行信号感知与传递，参与芽的分化、根的生长、种子发育等过程，且在低温、干旱及高盐等非生物胁迫的逆境响应中也发挥重要作用[1-4]。

植物双组分系统由组氨酸激酶（Histidine Kinase，HK）蛋白、组氨酸磷酸转移（Histidine Phosphotransfer，HPt）蛋白和反应调节（Response Regulator，RR）蛋白组成[5-6]。拟南芥（Arabidopsis thaliana）组氨酸激酶由AHK2、AHK3和AHK4（又名CRE1）编码，与细胞分裂素结合对下游的AHP（Arabidopsis HPts）蛋白和ARR（Arabidopsis Response Regulators）调节蛋白进行调节，进而传递植物激素（细胞分裂素、乙烯）和非生物胁迫信号[2,7-11]。研究表明，组氨酸激酶基因在植物的生长发育中发挥重要的调控作用，如AHK4对植物根的发育具有明显的调控作用[12]，AHK3则对植物叶片的衰老和细胞分化起主导作用[13-15]，AHK2和AHK3共同调节拟南芥的种子萌发和芽的生长[16]。另外，组氨酸激酶基因在植物应对低温、干旱及高盐等非生物胁迫过程中也发挥重要作用[3-4]。Kang 等[17]发现，拟南芥ahk2和ahk3单突变体比野生型个体抗旱性明显增强，并且ahk2/ahk3和ahk3/ahk4双突变体比相应的单突变体抗旱。低温胁迫时，ahk2/ahk3和ahk3/ahk4双突变体比野生型个体具有较高的抗冻性[18]，ahk2和ahk3单突变体对干旱、盐胁迫、低温和强光等非生物胁迫的抗性比野生型植株显著增强[3,18]。另外，通过萌发试验发现，ahk2、ahk3及ahk4单突变体对外源ABA 高度敏感，外源ABA 对拟南芥组氨酸激酶基因具有负调控作用[3]。通过ahk2/ahk3双突变体和野生型植株全基因组表达分析比较发现，组氨酸激酶基因在植物ABA 调控的抗逆反应及非ABA 调控的抗逆反应中均起重要作用[19]。在毛果杨（P.trichocarpaTorr.&Gray）基因组中鉴定出1 个AHK2同源基因PtHK2、2 个AHK3同源基因PtHK3a和PtHK3b，研究表明，毛果杨PtHK2、PtHK3a和PtHK3b在杨树形成层发育过程中具有重要调控作用[20]。然而，组氨酸激酶基因是否参与林木抗逆反应尚未见报道。

银腺杨‘84K’（P.alba×P.glandulosa‘84K’）是我国从韩国引进的白杨派优良品种，生长快、材质好、抗性强、适应性广，是优良的绿化树种、生态树种和用材树种，且因其易于组培及遗传转化，成为林木基因工程的理想材料。本研究以银腺杨‘84K’为材料，克隆了拟南芥AHK3同源基因（PaHK3b）的启动子及全长CDS（coding sequence），对其启动子元件及蛋白结构域进行了分析，并对PaHK3b基因在高温、干旱、盐胁迫等非生物胁迫及不同植物激素处理下的表达进行了qPCR 检测；同时，通过原核表达初步确定其生物学功能。该研究为利用组氨酸激酶基因进行杨树抗逆分子改良奠定了基础。

1 材料与方法

1.1 试验材料

银腺杨‘84K’。

1.2 试验方法

1.2.1PaHK3b基因的启动子克隆及序列分析 根据Phytozome 网站上公布的毛果杨基因组信息，查找基因上游启动子序列，设计PaHK3b基因启动子区域特异性引物PaHK3b-pro-F 和PaHK3b-pro-R（表1）。采用植物基因组DNA 提取试剂盒（TIANGEN，北京）提取‘84K’杨叶片总DNA，用常规PCR 方法进行目的片段扩增，用1.5%琼脂糖凝胶进行电泳检测，用胶回收试剂盒（Axygen，美国）回收并纯化目的片段，并连接到克隆载体pMDTM19-T Vector（Takara，日本）转化至大肠杆菌DH5α（TIANGEN，北京），经过蓝白斑培养基（Amp+）筛选，挑取阳性单克隆送至生物公司（中美泰和，北京）进行序列测定。利用在线软件plantCARE（http://bioinformatics.psb.ugent.be/webtools/plantcare/html/）分析该基因启动子中含有的顺式作用元件。

1.2.2PaHK3b基因CDS 的克隆及序列分析 采用EASY spin Plus 植物RNA 快速提取试剂盒（Aidlab，北京）提取‘84K’杨叶片总RNA，并用PrimeScript RT reagent Kit with gDNA Eraser(Perfect Real Time)（takara，日本）试剂盒进行反转录，合成单链cDNA。根据Phytozome 网站上公布的毛果杨基因组信息，设计PaHK3b基因全长cDNA 的PCR 扩增引物PaHK3b-F 和PaHK3b-R（表1），以单链cDNA 为模板进行PCR 扩增，扩增程序为95℃ 5 min；95℃ 30 s，57℃ 30 s，72℃ 3 min，38 个循环；72℃ 7 min；4℃保温。用1.5%琼脂糖凝胶进行电泳检测，用胶回收试剂盒（Axygen，美国）回收纯化目的片段，并连接到克隆载体PLB-Vector（TIANGEN，北京）转化至大肠杆菌DH5α，涂含抗性（Amp+）平板，37℃过夜培养。挑取阳性单克隆送至生物公司（中美泰和，北京）进行序列测定。利用在线软件ExpasyProtParam（http://web.expasy.org/protparam/）对PaHK3b基因序列特征进行分析，使用GOR4 在线工具(https://npsa-prabi.ibcp.fr/cgi-bin/secpred_gor4.pl)进行蛋白的二级结构进行预测，利用NCBI 在线工具（https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi）进行蛋白结构域分析。

1.2.3 ‘84K’杨非生物胁迫处理和激素处理 将带有顶端的‘84K’杨组培苗嫩茎（2～3 cm）切下，转接到生根培养基 (1/2MS+0.02 mg·L−1NAA+0.05 mg·L−1IBA)中，培养于温度24℃、光周期16 h/8 h（光照/黑暗）、光照强度为50 μmol·m−2·s−1人工气候培养室，28 d 后，选取生长状态一致的组培苗，转移到装有5 mL 1/2MS 液体培养基的玻璃管（直径4 cm，高20 cm）中继续培养。培养5 d后，进行非生物胁迫处理和植物激素处理。非生物胁迫处理分别为42℃高温、0℃低温、200 mmol·L−1NaCl 和5% PEG6000 处理，植物激素为脱落酸（ABA）、6-苄氨基嘌呤（6-BA）、吲哚丁酸（IBA）、赤霉素（GA3）及水杨酸（SA），浓度均为10 μmol·L−1。NaCl、PEG6000 及植物激素处理方法为将药剂添加到1/2MS 液体培养基中，以培养于1/2MS 培养基中的组培苗为对照。各处理的时间均为3 h，每个处理3 个生物学重复，处理后取成熟叶片，液氮速冻保存于 −80℃超低温冰箱。

表1 相关引物序列及产物大小Table 1 The sequences and PCR product size of primers used in this study

1.2.4PaHK3b基因的表达分析 提取‘84K’杨根、茎、叶以及不同处理植株叶片总RNA，反转录为cDNA。设计PaHK3b基因实时定量PCR（qRT-PCR）引物PaHK3b-q-F 和PaHK3b-q-R（表1）。将合成的cDNA 稀释10 倍作为实时定量PCR 模板，参照TB Green TM Premix Ex TaqTMII（TliRNaseH Plus）（TaKaRa，日本）试剂盒说明配制反应体系：TB Green Premix TaqII（TliRNaseH Plus，2 × ）10 μL，Primer-F(10 μmol·L−1) 0.8 μL，Primer-R(10 μmol·L−1) 0.8 μL，cDNA 模板2 μL，ddH2O 补足至20 μL。用Roche Light Cycle 480Ⅱ型荧光定量PCR 仪（Roche，瑞士）进行qRTPCR 反应，反应程序为：预变性95℃30 s，变性95℃ 5 s，退火60℃ 30 s，40 个循环；溶解曲线为95℃ 5 s，65℃ 1 min。以Actin为内参基因，采用2−ΔΔCT算法计算PaHK3b基因的相对表达量[21]。利用Excel2010 和Spass23.0 软件对数据进行整理制图和差异分析。

1.2.5PaHK3b基因原核表达载体的构建 根据TIANGEN 公司在线无缝克隆引物设计工具（http://123.56.75.19）设计含酶切位点的原核表达载体引物pET-28a-PaHK3b-F 和pET-28a-PaHK3b-R（表1），PCR 扩增‘84K’杨PaHK3b基因的cDNA 序列，回收扩增产物。原核表达载体pET-28a（索莱宝，北京）经Hind III（Thermo Fisher Scientific，美国）单酶切，胶回收纯化，根据EasyGeno 快速重组克隆试剂盒（TIANGEN，北京）操作流程，将二者回收产物进行重组反应，将重组质粒转化至大肠杆菌DH5α 感受态细胞中，经抗性（Kan+）筛选，挑取单克隆进行PCR 验证，并送至生物公司测序验证，提取阳性菌株质粒，获得重组载体（pET-28a-PaHK3b）并转化至表达菌株BL21（DE3）（TIANGEN，北京）。

1.2.6 重组菌E.coliBL21（pET-28a-PaHK3b）的盐及干旱胁迫处理 将转化后的表达菌株BL21（DE3）接种于5 mL 含有50 mg·L−1的卡那霉素的LB 液体培养基，37℃摇床震荡培养至OD6000.6～0.8，加入终浓度为0.5 mmol·L−1的IPTG 诱导4 h，稀释菌液104倍，取5 μL 至含有0、50、100、150、200 mmol·L−1NaCl LB 固体培养基上涂抹直径约为1 cm的9 个点，超净台吹干菌液，37℃倒置过夜培养[22]。同时，取诱导后的产物1 mL 加入到1 mL 含有PEG6000 的LB 液体培养基中使PEG 终浓度为5%（W/V），放置37℃摇床震荡培养，0～7 h 内每隔1 h 测1 次OD600值，每个时间点测3 个重复[23-25]。

2 结果与分析

2.1 ‘84K’杨PaHK3b 基因CDS 和启动子的获得及生物信息学分析

通过PCR 扩增、克隆得到‘84K’杨PaHK3b基因全长CDS。CDS 全长为3 060 bp，编码1 019 个氨基酸，蛋白分子量为113 564.85 kDa，等电点为6.53，蛋白不稳定指数为37，为稳定蛋白；脂肪族指数为92.31，蛋白疏水性平均值为−0.088。蛋白主要以无规则卷曲为主，α-螺旋与延伸链则散布在蛋白中，其中，无规则卷曲416 个(40.82%)、α-螺旋377 个(37%)、延伸链则有226 个(22.18%)。蛋白结构域分析结果表明，PaHK3b 蛋白具有典型的细胞分裂素受体结构域，包含CHASE (Cyclaseand histidine kinase associated sensing extracellular)结构域、HisKA (His kinase A domain)结构域及REC(CheY-homologous receiver domain)结构域，其中，CHASE 结构域使细胞分裂素与其相应的受体结合，HisKA 结构域感知信号和His 残基自身磷酸化，REC 结构域与磷酸基团相结合，并向下游的信号蛋白进行传递。

克隆得到全长1 600 bp 的‘84K’杨PaHK3b基因启动子序列，并采用plantCARE 对其进行生物信息学分析。结果表明：PaHK3b启动子不仅具有TATA-box 和CAAT-box 核心基本元件，还包括G-box、GA-motif、TCT-motif 等光响应元件。另外，PaHK3b启动子区域还包含多个与逆境响应及激素调控相关的顺式作用元件，如低温响应元件LTR、防御与胁迫响应元件TC-rich repeats、赤霉素响应元件GARE-motif、水杨酸响应元件TCAelement 等（表2），表明该基因可能参与杨树逆境胁迫及激素信号响应。

表2 PlantCARE 预测PaHK3b 启动子的主要顺式作用元件Table 2 The main cis-elements in PaHK3b promoter predicted by PlantCARE

2.2 ‘84K’杨PaHK3b 基因器官表达差异性

采用qRT-PCR 技术对‘84K’杨PaHK3b基因在根、茎、叶的表达进行检测。结果表明：该基因在根、茎、叶中均有表达，以根表达量作为参照，茎表达量低于根组织，为根表达量的82%，叶片的表达量高于根，为根表达量的1.63 倍（图1）。以上结果表明，PaHK3b基因在‘84K’杨叶片表达量最高。

图1 ‘84K’杨PaHK3b 基因在根、茎、叶片中的表达差异Fig.1 Expression difference of PaHK3b gene in root，stem and leaf of ‘84K’ poplar

2.3 非生物胁迫下‘84K’杨PaHK3b 基因的表达情况

高温、低温、盐胁迫和干旱胁迫条件下，‘84K’杨PaHK3b基因表达量均明显高于对照，并且温度胁迫下PaHK3b基因表达量高于盐胁迫和干旱胁迫下的表达量（图2），其中，42℃高温胁迫、0℃低温胁迫时，PaHK3b基因表达量分别为对照的2.67、2.61 倍；200 mmol·L−1NaCl 胁迫、5%PEG6000 胁迫时，PaHK3b基因表达量分别为对照的2.28、1.87 倍。以上结果说明，PaHK3b基因参与了非生物胁迫响应，在杨树的逆境胁迫反应过程中起作用。

图2 非生物胁迫处理下‘84K’杨叶片PaHK3b 基因的表达差异Fig.2 Expression difference of PaHK3b in leaf of ‘84K’poplar under abiotic stresses

2.4 ‘84K’杨PaHK3b 基因对外源生长素及细胞分裂素处理的表达反应

图3 表明：PaHK3b基因对不同植物激素的响应有差异。生长素IBA 处理时，PaHK3b基因表达量与对照相似，而其他激素处理时，PaHK3b基因表达量均不同程度的低于对照，其中，6-BA、ABA、GA3、SA 处理时，PaHK3b基因表达量分别为对照的35%、39%、61%、66%。可见，外源赤霉素、脱落酸、细胞分裂素及水杨酸均能抑制杨树PaHK3b基因的表达，表明PaHK3b基因参与植物激素响应。

图3 不同植物激素处理下‘84K’杨叶片PaHK3b 基因的表达差异Fig.3 Expression difference of ‘84K’PaHK3b under several plant hormone treatments

2.5 重组菌E.coli BL21（pET-28a-PaHK3b）的抗逆检测

利用斑点法对PaHK3b基因的抗盐性进行了检测。结果表明：在未添加NaCl 的LB 固体培养基上，转入pET-28a-PaHK3b重组质粒的大肠杆菌BL21 菌株（重组菌株）与转入pET-28a 空载体的菌株（对照菌株）生长速度基本一致。在含有50、100、150 mmol·L−1的NaCl LB 固体培养基上，重组菌株与对照菌株生长均受到抑制，且随着NaCl 浓度的增加菌落生长减少；但重组菌生长速度较对照快，菌落多且菌斑大。当培养基中NaCl 浓度为200 mmol·L−1时，2 种菌株均停止生长（图4）。因此，PaHK3b基因通过原核表达能够提高大肠杆菌菌株的抗盐能力。

图4 大肠杆菌E.coli BL21（pET-28a）和E.coli BL21（pET-28a-PaHK3b）菌株耐盐性差异Fig.4 Difference in salt tolerance of E.coli BL21（pET-28a）and E.coli BL21（pET-28a-PaHK3b）

利用PEG 模拟干旱胁迫对PaHK3b基因的抗旱性进行评价，结果表明：在LB 液体培养基中，重组菌株的生长速度较对照菌株要快；在添加5%PEG6000 的LB 液体培养基中，重组菌株和对照菌株生长均受到抑制。培养1 h 时，重组菌株和对照菌株生长没有差异，之后转入pET-28a-PaHK3b重组质粒的菌株较转入pET-28a 空载的大肠杆菌生长速度快，培养7 h 后，含重组质粒菌株的OD600值是初始的1.8 倍，含空载质粒菌株的OD600值是初始的1.1 倍（图5）。因此，PaHK3b基因原核表达不仅能够提高大肠杆菌菌株生长速度，而且还能够提高其抗脱水能力。

图5 大肠杆菌E.coli BL21（pET-28a）和E.coli BL21（pET-28a-PaHK3b）菌株在LB 及LB-5%PEG6000 液体培养基生长曲线Fig.5 Growth curve of E.coli BL21（pET-28a）and E.coli BL21（pET-28a-PaHK3b）in LB and LB-5%PEG6000 liquid medium respectively

3 讨论

组氨酸激酶基因在植物激素信号响应及非生物逆境胁迫响应中发挥重要作用[3,26]。如拟南芥组氨酸激酶AHK3参与其对干旱、高盐及冷胁迫的抗逆反应[3-4]，转玉米（Zea maysL.）组氨酸激酶基因ZmHK9的拟南芥抗旱性明显增强[27]。本研究中，在高温、低温、高盐和PEG 模拟干旱等处理条件下，PaHK3b基因表达量均增加，其中，基因上调程度呈高温>低温>盐胁迫>干旱胁迫的模式；在6-BA、ABA、GA3 及SA 处理条件下，PaHK3b基因表达量均下降，基因下调程度呈6-BA>ABA>GA3>SA 的模式。因此表明，PaHK3b基因能够不同程度的参与植物激素信号响应及非生物胁迫信号响应。另外，SA 是植物防卫反应的重要内源植物激素，作为一种信号分子，在植物生物胁迫响应过程中发挥重要作用[28-29]。ABA 在植物抗逆胁迫反应及在生物胁迫和非生物胁迫中均具有重要调控作用。因此，PaHK3b基因可能通过ABA 及SA 途径的负调节提高了杨树对逆境胁迫的耐受性。

通过分析PaHK3b基因CDS 序列发现，‘84K’杨组氨酸激酶PaHK3b 蛋白具有典型的CHASE 细胞分裂素受体结构域。在植物外源激素6-BA 处理下，PaHK3b基因下调表达，表明该基因能够感知细胞分裂素信号并参与响应，这与组氨酸激酶基因具有CHASE 结构域，通过该结构域与细胞分裂素的结合，感知植物激素信号并进行传导是一致的[7,10]。对PaHK3b基因启动子序列分析发现，该区域包含低温响应、防御与胁迫响应、赤霉素响应、水杨酸响应等多个参与胁迫响应及激素调控的顺式作用元件。这与组氨酸激酶基因启动子区域有多个植物激素信号响应及非生物逆境胁迫相关的顺式作用元件是一致的[30-31]。另外，结合在植物激素6-BA、ABA、GA3、SA 和非生物胁迫高温、低温、高盐、PEG 处理下，PaHK3b基因均不同程度参与响应的表达，表明PaHK3b基因参与植物激素信号响应及非生物胁迫响应与其基因结构密切相关。

已有研究表明，通过构建基因的原核表达载体，并利用大肠杆菌表达系统可以对基因的功能进行初步鉴定。如将小麦（Triticum aestivumL.）脱水蛋白编码基因——DHN14基因的原核表达载体转入大肠杆菌后，大肠杆菌菌株对重金属胁迫及过氧化氢胁迫的耐受性显著增强，推测小麦脱水蛋白基因DHN14在其响应非生物胁迫的过程中起作用[23]。转入谷子（Setaria italicaL.）类受体蛋白激酶基因SiRLK35原核表达载体的大肠杆菌菌株生长状态较阴性对照好，同时获得的转SiRLK35基因水稻植株对盐胁迫的耐受性高于对照[22]。以上研究说明，原核表达结果能对克隆基因功能进行有效鉴定。在本研究中，转入PaHK3b基因原核表达载体的菌株在50～150 mmol·L−1NaCl 的LB 固体培养基及5% PEG6000 的LB 液体培养基中生长状态均好于转入空载体的对照菌株。因此可得，杨树PaHK3b基因能够促进大肠杆菌的生长，同时能够提高其耐盐和抗旱能力，推测该基因可能在杨树生长及抗逆境胁迫的过程中发挥作用。

4 结论

本研究从‘84K’杨中克隆了组氨酸激酶基因PaHK3b及其启动子，生物信息学分析表明，该基因编码蛋白具有典型的细胞分裂素受体结构域，其启动子区域有多个参与胁迫响应和激素信号转导的顺式作用元件。在高温、低温、高盐和PEG 模拟干旱等非生物胁迫下，PaHK3b基因表达量增加，而在6-BA、ABA、GA3 及水杨酸等植物激素处理条件下，其表达量下降。原核表达PaHK3b基因可以提高受体大肠杆菌的生长及抗旱耐盐能力。初步生物学功能研究表明，该基因参与杨树植物激素信号响应，并在其抗逆境胁迫过程中发挥重要调控作用。本研究结果为深入研究PaHK3b基因在激素信号转导及逆境胁迫下杨树生长发育的调控作用提供线索，为杨树抗逆分子育种及品种改良奠定基础。