-
细胞分裂素是植物生长过程中的一类重要调控激素，其主要通过双组分系统进行信号感知与传递，参与芽的分化、根的生长、种子发育等过程，且在低温、干旱及高盐等非生物胁迫的逆境响应中也发挥重要作用[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. trichocarpa Torr. & 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检测；同时，通过原核表达初步确定其生物学功能。该研究为利用组氨酸激酶基因进行杨树抗逆分子改良奠定了基础。
HTML
-
银腺杨‘84K’。
-
根据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/ ）分析该基因启动子中含有的顺式作用元件。基因 Gene 引物序列 Primer sequence (5'→3') 产物长度 Product length/bp PaHK3b-F ATGAGTTTGCTTCATGTATTTGGG 3 060 PaHK3b-R TAAATTCGCACTGGAAGTTGGC PaHK3b-pro-F GTGGCTTTGGGAGGAAGGAAT 1 600 PaHK3b-pro-R CCACATCTTGCCACCATCACC PaHK3b-q-F ATGAGTTTGCTTCATGTATTTGG 243 PaHK3b-q-R CTATAAATTCGCACTGGAAGTTG Actin-F AAACTGTAATGGTCCTCCCTCCG 193 Actin-R GCATCATCACAATCACTCTCCGA pET-28a-PaHK3b-F cgagctccgtcgacaagcttATGAGTTTGCTTCATGTATTTGGG 3 060 pET-28a-PaHK3b-R tcgagtgcggccgcaagcttTAAATTCGCACTGGAAGTTGGC Table 1. The sequences and PCR product size of primers used in this study
-
采用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 ）进行蛋白结构域分析。 -
将带有顶端的‘84K’杨组培苗嫩茎（2~3 cm）切下，转接到生根培养基(1/2MS+0.02 mg·L−1 NAA+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−1 NaCl和5% PEG6000处理，植物激素为脱落酸（ABA）、6-苄氨基嘌呤（6-BA）、吲哚丁酸（IBA）、赤霉素（GA3）及水杨酸（SA），浓度均为10 μmol·L−1 。NaCl、PEG6000及植物激素处理方法为将药剂添加到1/2MS液体培养基中，以培养于1/2MS培养基中的组培苗为对照。各处理的时间均为3 h，每个处理3个生物学重复，处理后取成熟叶片，液氮速冻保存于 −80℃超低温冰箱。
-
提取‘84K’杨根、茎、叶以及不同处理植株叶片总RNA，反转录为cDNA。设计PaHK3b基因实时定量PCR（qRT-PCR）引物PaHK3b-q-F和PaHK3b-q-R（表1）。将合成的cDNA稀释10倍作为实时定量PCR模板，参照TB Green TM Premix Ex TaqTM II（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，瑞士）进行qRT-PCR反应，反应程序为：预变性95℃30 s，变性95℃ 5 s，退火60℃ 30 s，40个循环；溶解曲线为95℃ 5 s，65℃ 1 min。以Actin为内参基因，采用2−ΔΔCT算法计算PaHK3b基因的相对表达量[21]。利用Excel2010和Spass23.0软件对数据进行整理制图和差异分析。
-
根据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，北京）。
-
将转化后的表达菌株BL21（DE3）接种于5 mL含有50 mg·L−1的卡那霉素的LB液体培养基，37℃摇床震荡培养至OD600 0.6～0.8，加入终浓度为0.5 mmol·L−1的IPTG诱导4 h，稀释菌液104倍，取5 μL至含有0、50、100、150、200 mmol·L−1 NaCl 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]。
1.1. 试验材料
1.2. 试验方法
1.2.1. PaHK3b基因的启动子克隆及序列分析
1.2.2. PaHK3b基因CDS的克隆及序列分析
1.2.3. ‘84K’杨非生物胁迫处理和激素处理
1.2.4. PaHK3b基因的表达分析
1.2.5. PaHK3b基因原核表达载体的构建
1.2.6. 重组菌E.coli BL21（pET-28a-PaHK3b）的盐及干旱胁迫处理
-
通过PCR扩增、克隆得到‘84K’杨PaHK3b基因全长CDS。CDS全长为3060 bp，编码1019个氨基酸，蛋白分子量为113564.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结构域与磷酸基团相结合，并向下游的信号蛋白进行传递。
克隆得到全长1600 bp的‘84K’杨PaHK3b基因启动子序列，并采用plantCARE对其进行生物信息学分析。结果表明：PaHK3b启动子不仅具有TATA-box和CAAT-box核心基本元件，还包括G-box、GA-motif、TCT-motif等光响应元件。另外，PaHK3b启动子区域还包含多个与逆境响应及激素调控相关的顺式作用元件，如低温响应元件LTR、防御与胁迫响应元件TC-rich repeats、赤霉素响应元件GARE-motif、水杨酸响应元件TCA-element等（表2），表明该基因可能参与杨树逆境胁迫及激素信号响应。
顺式作用元件
Cis-element起始位点（方向）
Start site/bp（direction）序列
Sequence特性
CharacteristicARE 1 150（+） AAACCA 厌氧诱导 Essential for the anaerobic induction GARE-motif 346（−） TCTGTTG 赤霉素响应 Gibberellin-responsive element LTR 318（−） CCGAAA 低温响应 Low-temperature responsiveness TC-rich repeats 88（+） ATTCTCTAAC 防卫及胁迫响应 Defense and stress responsiveness TCA-element 288（−）/818（+） CCATCTTTTT 水杨酸响应 Salicylic acid responsiveness GATA-motif 171（−） GATAGGA 光响应元件 Light responsive element GATA-motif 439（+） GATAGGG TCCC-motif 1511（+） TCTCCCT TCT-motif 1 462（+）/1 551（−） TCTTAC GT1-motif 677（+） GTGTGTGA GT1-motif 1 409（−）/1 353（−） GGTTAA GA-motif 706（+） ATAGATA G-box 9（−） CACGAC CAAT-box 619（−）/575（−）/993（−）/103（−） CCAAT 启动子和增强子区域常见的顺式作用元件
Common cis-acting element in promoter and enhancer regionsTATA-box 27（+）/195（+）/150（+）/1 067（−） TATA 核心启动子
Core promoter element around −30 of transcription start148（−）/1 065（+）/193（−） TATAAA 147（−） TATAAAT 136（+） ATATAA 1 109（−） ccTATAAAaa 192（−） TATAAAA Table 2. The main cis-elements in PaHK3b promoter predicted by PlantCARE
-
采用qRT-PCR技术对‘84K’杨PaHK3b基因在根、茎、叶的表达进行检测。结果表明：该基因在根、茎、叶中均有表达，以根表达量作为参照，茎表达量低于根组织，为根表达量的82%，叶片的表达量高于根，为根表达量的1.63倍（图1）。以上结果表明，PaHK3b基因在‘84K’杨叶片表达量最高。
Figure 1. Expression difference of PaHK3b gene in root，stem and leaf of ‘84K’ poplar
-
高温、低温、盐胁迫和干旱胁迫条件下，‘84K’杨PaHK3b基因表达量均明显高于对照，并且温度胁迫下PaHK3b基因表达量高于盐胁迫和干旱胁迫下的表达量（图2），其中，42℃高温胁迫、0℃低温胁迫时，PaHK3b基因表达量分别为对照的2.67、2.61倍；200 mmol·L−1 NaCl胁迫、5%PEG6000胁迫时，PaHK3b基因表达量分别为对照的2.28、1.87倍。以上结果说明，PaHK3b基因参与了非生物胁迫响应，在杨树的逆境胁迫反应过程中起作用。
Figure 2. Expression difference of PaHK3b in leaf of ‘84K’ poplar under abiotic stresses
-
图3表明：PaHK3b基因对不同植物激素的响应有差异。生长素IBA处理时，PaHK3b基因表达量与对照相似，而其他激素处理时，PaHK3b基因表达量均不同程度的低于对照，其中，6-BA、ABA、GA3、SA处理时，PaHK3b基因表达量分别为对照的35%、39%、61%、66%。可见，外源赤霉素、脱落酸、细胞分裂素及水杨酸均能抑制杨树PaHK3b基因的表达，表明PaHK3b基因参与植物激素响应。
Figure 3. Expression difference of ‘84K’PaHK3b under several plant hormone treatments
-
利用斑点法对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基因通过原核表达能够提高大肠杆菌菌株的抗盐能力。
Figure 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基因原核表达不仅能够提高大肠杆菌菌株生长速度，而且还能够提高其抗脱水能力。
Figure 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
2.1. ‘84K’杨PaHK3b基因CDS和启动子的获得及生物信息学分析
2.2. ‘84K’杨PaHK3b基因器官表达差异性
2.3. 非生物胁迫下‘84K’杨PaHK3b基因的表达情况
2.4. ‘84K’杨PaHK3b基因对外源生长素及细胞分裂素处理的表达反应
2.5. 重组菌E.coli BL21（pET-28a-PaHK3b）的抗逆检测
-
组氨酸激酶基因在植物激素信号响应及非生物逆境胁迫响应中发挥重要作用[3, 26]。如拟南芥组氨酸激酶AHK3参与其对干旱、高盐及冷胁迫的抗逆反应[3-4]，转玉米（Zea mays L.）组氨酸激酶基因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 aestivum L.）脱水蛋白编码基因——DHN14基因的原核表达载体转入大肠杆菌后，大肠杆菌菌株对重金属胁迫及过氧化氢胁迫的耐受性显著增强，推测小麦脱水蛋白基因DHN14在其响应非生物胁迫的过程中起作用[23]。转入谷子（Setaria italica L.）类受体蛋白激酶基因SiRLK35原核表达载体的大肠杆菌菌株生长状态较阴性对照好，同时获得的转SiRLK35基因水稻植株对盐胁迫的耐受性高于对照[22]。以上研究说明，原核表达结果能对克隆基因功能进行有效鉴定。在本研究中，转入PaHK3b基因原核表达载体的菌株在50~150 mmol·L−1 NaCl的LB固体培养基及5% PEG6000的LB液体培养基中生长状态均好于转入空载体的对照菌株。因此可得，杨树PaHK3b基因能够促进大肠杆菌的生长，同时能够提高其耐盐和抗旱能力，推测该基因可能在杨树生长及抗逆境胁迫的过程中发挥作用。
-
本研究从‘84K’杨中克隆了组氨酸激酶基因PaHK3b及其启动子，生物信息学分析表明，该基因编码蛋白具有典型的细胞分裂素受体结构域，其启动子区域有多个参与胁迫响应和激素信号转导的顺式作用元件。在高温、低温、高盐和PEG模拟干旱等非生物胁迫下，PaHK3b基因表达量增加，而在6-BA、ABA、GA3及水杨酸等植物激素处理条件下，其表达量下降。原核表达PaHK3b基因可以提高受体大肠杆菌的生长及抗旱耐盐能力。初步生物学功能研究表明，该基因参与杨树植物激素信号响应，并在其抗逆境胁迫过程中发挥重要调控作用。本研究结果为深入研究PaHK3b基因在激素信号转导及逆境胁迫下杨树生长发育的调控作用提供线索，为杨树抗逆分子育种及品种改良奠定基础。
DownLoad: