棉花半乳糖基转移酶基因GhGalT1启动子的克隆及表达分析

秦丽霞 李 静 张换样 李 盛 竹梦婕 焦改丽 吴慎杰

山西省农业科学院棉花研究所, 山西运城 044000

糖基转移酶(glycosyltransferase, GT, EC 2.4.x.y)广泛存在于植物体内, 负责催化小分子化合物(激素、次级代谢物等)的糖基化, 将活性糖基从供体转移到受体(小分子化合物), 从而改变这些化合物的生物化学性质[1-2]。按照序列同源性及催化底物特性及糖苷键立体化学性质将GTs分为94个不同的家族(http://www.cazy.org/fam/acc_GT.html)[3]。越来越多的证据显示GTs在调节植物激素和植物抗逆反应中发挥重要作用。

在拟南芥中过量表达 UGT84B1导致转基因植株表现出生长素缺陷表型, 植株矮化, 多分枝, 叶片皱缩, 根系失去向地性生长[4-5]。UGT74E2在体外能糖基化生长素IBA, UGT74E2过量表达转基因拟南芥植株表现出矮化、多分枝, 抗盐、抗旱性增强[6]。在烟草中过量表达玉米素合成相关糖基转移酶ZOG1, 导致转基因比非转基因材料需多添加 10倍的玉米素才能诱导愈伤组织形成芽[7]。拟南芥UGT76C1、UGT76C2、UGT85A1、UGT73C1 和UGT73C5参与调控植物对细胞分裂素的响应[8-9]。UGT73C5和 UGT73C6催化油菜素内酯的糖基化,在拟南芥中过量表达 UGT73C5植株表现出油菜素内酯缺陷表型, 同时伴随着高水平的油菜素内酯糖苷物的积累, 而在突变体株系中则未检测到油菜素内酯糖苷物[10-12]。从拟南芥分离出的叶片伤口处检测到茉莉酸糖苷的富集, 表明 GTs在激素调节和病原菌防御之间也起重要作用[13]。UGT73B1、UGT73B2、UGT73B3和UGT71C1功能缺失突变体增强了植物对氧化胁迫的抗性[14-16]。在拟南芥中过量表达水杨酸合成糖基转移酶 AtSGTl导致植株体内水杨酸及其糖苷水平降低, 进而增加对半活体营养侵染型病原菌假单胞菌的敏感性[17]。虽然以上研究已表明糖基转移酶参与植物耐逆过程, 但对其具体作用机理尚不清楚, 糖基转移酶参与植物耐逆过程的分子机制有待深入探索。

非生物胁迫(干旱、土壤盐渍化、低温冷害等)是限制农业生产、威胁作物产量和品质的重要因素,而传统育种方法并不能精准有效地提高植物耐受逆境能力。因此, 通过植物基因工程技术结合传统育种手段, 使抗逆基因在目标作物中特异表达, 以提高作物对逆境胁迫的抗性, 逐渐成为作物育种中的一个重要手段。启动子作为启动目标基因表达的关键调控元件, 主要有组成型启动子、组织特异型启动子和诱导型启动子[18]。虽然组成型启动子在启动外源基因表达中具有高效、广谱和稳定性强等优点,但是由于高效、持续表达目标基因容易造成一些负面影响(如植株生长异常等), 而诱导型启动子是在特定条件下启动目标基因的表达, 很大程度上缓解了由于大量异源蛋白的高表达对植物生长造成的伤害[19-20]。

棉花的生长和发育过程经常受到高盐、干旱等非生物胁迫因素的影响, 导致产量降低。到目前为止, 在棉花中参与调节激素平衡及抗逆胁迫方面的糖基转移酶尚未见报道。我们之前从棉花cDNA文库中分离了一个棉花糖基转移酶 GhGalT1基因[21],为研究其在植物逆境应答中的功能, 本研究分离克隆了GhGalT1基因起始密码子ATG上游539 bp的启动子序列。通过PlantCARE软件分析发现其包含许多与逆境、激素等相关的顺式作用元件, 推测该启动子受干旱、热、脱水、激素等胁迫诱导表达, 为此我们构建了含有该启动子的重组表达载体pBI121-pGhGalT1-GUS, 并通过浸花法转化拟南芥对其功能进行了验证, 以期丰富植物诱导型启动子种类, 为植物抗逆遗传育种改良提供理论依据。

1 材料与方法

1.1 植物材料

利用本实验室保存的陆地棉(Gossypium hirsutum L.)品种Coker 312克隆GhGalT1基因的启动子,拟南芥(Arabidopsis thaliana, 哥伦比亚型)用于启动子重组质粒的转化与功能验证。

1.2 棉花基因组DNA的提取

选取饱满的Coker 312棉花种子, 以70%乙醇处理1 min, 10%过氧化氢(H2O2)处理1.5 h, 无菌水冲洗 3～4次, 随后加入适量无菌水浸泡至种子根尖露白。去掉种皮, 置1/2 MS培养基上萌发, 28ºC培养5～7 d后, 选取生长健壮的幼苗移栽至土壤中, 待幼苗长出2片真叶时, 取2 g左右幼嫩叶片利用CTAB法提取基因组DNA[22]。

1.3 GhGalT1基因启动子的克隆及顺式作用元件分析

以GhGalT1编码区序列为探针, Blast搜索中国农业科学院棉花研究所网站的陆地棉(Gossypium hirsutum)基因组序列数据库(http://cgp.genomics.org.cn/page/species/blast.jsp), 将该基因上游2 kb左右序列作为目标区域。

利用在线顺式作用元件分析工具 PLACE(http://www.dna.affrc.go.jp/PLACE/signalup.htm)和PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)分析启动子序列, 预测顺式作用元件。

以 Coker 312的基因组 DNA为模板, 使用PrimeSTAR HS DNA Polymerase高保真酶, PCR扩增 GhGalT1基因的启动子目的片段, 命名为pGhGalT1, 琼脂糖凝胶电泳检测扩增条带大小。扩增体系为 5×Prime STAR 缓冲液 (Mg2+plus) 4 μL,dNTPs (各 2.5 mmol L–1) 4 μL, 模板 DNA (500 ng μL–1) 1 μL, 引物 pGhGalT1-F (10 μmol L–1) 1 μL,pGhGalT1-R (10 μmol L–1) 1 μL (表 1), PrimeSTAR HS DNA 聚合酶 (2.5 U μL–1) 0.2 μL, 用 ddH2O 补至20 μL。反应程序为 98ºC 10 s, 55ºC 15 s, 72ºC 1 min,30个循环; 72ºC 10 min。将扩增得到的序列连接至pBluescript(pSK)载体上测序验证, 与棉花基因组数据库序列比对。

表1 本试验中所用引物序列Table 1 Primers used in this study

1.4 GhGalT1启动子载体构建、拟南芥转化及阳性植株检测

用 Sal I和 BamH I分别酶切植物表达载体pBI101和扩增的 GhGalT1基因启动子目的片段pGhGalT1后, 进行连接, 获得重组质粒 pBI101-pGhGalT1-GUS (图1), 转化农杆菌GV3101菌株。

采用农杆菌介导的浸花法转化拟南芥[23]。经卡那霉素抗性筛选后获得 T1代转基因拟南芥。取 T1代转基因植株幼苗叶片, 利用 SDS法提取基因组DNA, 并以此为模板通过 PCR扩增检测阳性植株,扩增引物为 pGhGalT1-F和 pGhGalT1-R、GUS-F和GUS-R和抗性标记基因NPTII-F和NPTII-R (表1), 产物大小分别为539、1812和795 bp。选取PCR检测呈阳性的单株收获种子并加代, 再经卡那霉素抗性筛选及PCR检测得到T3代纯合转基因阳性株系。

图1 pBI101-pGhGalT1-GUS表达载体构建示意图Fig. 1 Map of construction of the pBI101-pGhGalT1-GUS expression vector

1.5 GUS组织化学染色

取野生型和 T3代转基因拟南芥植株种子, 经10% NaClO溶液表面消毒后, 分别点播于MS培养基上, 4ºC春化处理2 d, 于恒温光照培养箱(16 h光照/8 h 黑暗, 22～24ºC)培养, 分别取生长 5、10 和 15 d幼苗、茎、花和角果、成熟叶片等组织进行GUS组织化学染色[24]。

选取长势一致的幼苗浸入 GUS染液(含 10 mmol L–1EDTA, 100 mmol L–1NaH2PO4·H2O, 0.5%Triton X-100, 1 mmol L–1K3Fe(CN)6, 1 mmol L–1K4Fe(CN)6·3H2O, 0.1% X-Gluc)中, 于 37ºC 染色 6～8 h, 除去染液。再用70%乙醇室温脱色3～5 h, 其间更换 70%乙醇 2～3次, 直至各组织的色素完全退去。取样品在Leica MZ16f体视显微镜下观察, 照相。

取适量生长10 d的转基因拟南芥幼苗分别移至含 150 mmol L–1NaCl、300 mmol L–1甘露醇、100 μmol L–16-BA、50 μmol L–1茉莉酸甲酯(MeJA)和500 nmol L–1油菜素内酯(BL)的MS液体培养基中处理6 h后, 进行GUS染色分析, 本试验重复3次。

1.6 GUS酶活测定

GUS活性的荧光检测方法参照Jefferson[25]。将生长10 d幼苗分别用MS液体和添加有150 mmol L–1NaCl、300 mmol L–1甘露醇、100 μmol L–16-BA、50 μmol L–1MeJA 或 500 nmol L–1BL 的 MS 液体培养基诱导6 h。取0.1 g样品于预冷的研钵中, 加0.5 mL预冷的研磨缓冲液(100 mmol L–1K3PO4, pH 7.8,1 mmol L–1EDTA, 1% Triton X-100), 迅速研磨成匀浆后转入预冷的1.5 mL离心管中, 2130×g, 4ºC离心10 min, 上清液即为GUS蛋白粗提取液。以4-MUG为作用底物进行酶促反应后, 反应产物在365 nm的激发光下测定455 nm的荧光强度(HITACHI F-2500),考马斯亮蓝法测定蛋白质含量。

1.7 GUS基因的表达分析

取适量生长10 d携有启动子序列的拟南芥幼苗分别移至含 150 mmol L–1NaCl、300 mmol L–1甘露醇、100 μmol L–16-BA、50 μmol L–1MeJA 或 500 nmol L–1BL的MS液体培养基中处理6 h后收取幼苗材料。用RNA-Solv Reagent (Omega)提取幼苗总RNA, 经DNase (RQ1 RNase-Free DNase, Promega)消化后, 用逆转录酶(M-MLV Reverse Transcriptase,Promega)合成cDNA。应用实时荧光定量RT-PCR技术(MJ Research, Opticon 2), 以 cDNA 为模板, 用GUS基因特异性引物GUS-RTF和GUS-RTR (表1)和 Real-time PCR Master Mix (TOYOBO)进行 PCR反应。以拟南芥的 Actin2作为内参基因, 目标基因每一个循环的扩增都被 SYBR-Green荧光检测, 具体步骤及分析参照Li等[26]的方法, 重复3次, 统计分析试验结果。

2 结果与分析

2.1 GhGalT1启动子的克隆及顺式作用元件分析

靠近GhGalT1编码区开放阅读框5′上游539 bp(图2)区域序列具有TATA-box、CAAT-box等真核生物启动子基本调控元件, 此外还包含多种逆境响应功能元件, 如干旱诱导的 MYB结合元件 MBS、脱水响应元件 MYCCONSE、热应答元件 HSE、防御与胁迫应答元件 TC-rich repeats、MeJA响应元件CGTCA-motif、调控生物节律元件以及多种光响应相关元件等(图3和表2)。

图2 棉花GhGalT1基因启动子的克隆Fig. 2 Amplified fragment of GhGalT1 gene promoter in cotton

图3 GhGalT1启动子序列及预测的顺式作用元件Fig. 3 Sequence and predicted cis-acting elements of GhGalT1 promoter

2.2 GhGalT1启动子载体构建及拟南芥转化与阳性植株的鉴定

为了研究GhGalT1启动子的活性及其表达模式,将该启动子扩增产物pGhGalT1 (图2)及植物表达载体pBI101分别用Sal I和Bam HI双酶切后连接, 经PCR扩增和酶切检测验证获得正确的重组质粒(图4-A和 B)。通过农杆菌介导法转化拟南芥, 获得转基因植株。经卡那霉素抗性筛选和PCR检测呈阳性的单株收获种子并加代(图5, 图6-A～C), 获得T3代纯合转基因阳性株系(图6-D)。

表2 GhGalT1启动子主要顺式作用元件及预测的功能Table 2 The cis-acting elements and predicted functions in the GhGalT1 promoter sequence

图4 GhGalT1启动子载体的构建Fig. 4 Vector construction of GhGalT1 promoter

图5 T1代转基因拟南芥阳性植株的卡那霉素抗性筛选Fig. 5 Kanamycin resistance screening of T1 generation transgenic positive plants

图6 转基因拟南芥阳性植株的PCR检测Fig. 6 Positive detection in transgenic Arabidopsis lines by PCR

2.3 GUS组织化学染色

图7表明GUS信号主要在生长5、10和 15 d转基因植株幼苗的主根及侧根的根尖表达, 在 10 d以后的子叶及莲座叶边缘也有GUS表达, 而在茎、花、角果等部位未检测到 GUS信号, 表明靠近GhGalT1编码区开放阅读框上游的 539 bp的GhGalT1启动子区域可以驱动GUS基因表达, 是启动子的活性区域。

图7 GhGalT1启动子转基因拟南芥的GUS活性检测Fig. 7 Histochemical assay of GUS activity in pGhGalT1:GUS transgenic Arabidopsis

2.4 不同胁迫条件下的GUS染色及定量分析

Genevestigator网站(https://genevestigator.com/gv/)分析显示, 与 GhGalT1基因同源的拟南芥AtGalTs基因受盐、干旱、激素(6-BA、MeJA、BL等)的诱导。将生长 10 d的幼苗移至分别添加 150 mmol L–1NaCl、300 mmol L–1甘露醇、100 μmol L–16-BA、50 μmol L–1MeJA 和 500 nmol L–1BL 的 MS液体培养基中处理6 h, 然后进行GUS染色分析。图 8显示, 与无胁迫对照相比, 经 NaCl、甘露醇和MeJA处理后, 转基因株系根中GUS信号明显减弱;用BL处理后, 根中GUS信号显著增强; 而在6-BA处理的小苗根中GUS信号完全消失, 子叶和莲座叶中的GUS信号却明显增强。

图8 盐、干旱和激素处理GhGalT1启动子转基因拟南芥的GUS活性检测Fig. 8 Histochemical assay of GUS activity of GhGalT1 promoter in transgenic Arabidopsis under stresses

进一步定量确定启动子的表达活性, 对盐、干旱、激素(6-BA、MeJA、BL)处理过的 10 d幼苗进行GUS酶活测定, 同时提取总RNA后, 对GUS基因进行实时荧光定量 PCR分析其表达量。结果经NaCl、甘露醇和MeJA处理后, GUS活性及GUS基因表达量显著降低, 6-BA、BL处理后GUS活性及GUS基因表达量显著升高(图9), 表明GhGalT1基因的启动子受盐、渗透胁迫和激素(6-BA、MeJA、BL等)的诱导。

图9 胁迫条件下转基因拟南芥GUS定量分析Fig. 9 Quantification of GUS activity and GUS gene in transgenic Arabidopsis under different stresses

3 讨论

真核生物基因表达受染色体水平上的基因活化调节、转录水平、转录后水平和翻译水平的调控, 外源基因在受体内的表达受到诸多因素的影响, 如插入位置、拷贝数、启动子类型等, 但很大程度上依赖于外源基因的启动子[27]。目前, 研究启动子的方法主要有生物信息学分析及试验分析两种。利用生物信息学分析启动子序列, 预测其中的顺式作用元件种类、数目及位置, 可以为进一步确定启动子的功能奠定基础。试验分析主要有瞬时表达、5′端片段缺失分析[28]、凝胶阻滞[29]、酵母单杂交[30]等方法。本研究采用生物信息学分析和启动子融合GUS基因瞬时表达相结合的方法, 研究GhGalT1基因启动子的重要顺式作用元件和表达特性。

通过生物信息学分析发现, 在 GhGalT1启动子区域存在多种重要的非生物胁迫应答元件, 例如响应干旱、防御与胁迫等逆境的元件。Genevestigator网站(https://genevestigator.com/gv/)分析也显示与GhGalT1基因同源的拟南芥 AtGalTs基因受盐、干旱、激素(6-BA、MeJA、BL等)的诱导。GUS组织化学染色表明邻近GhGalT1基因开放阅读框5′端上游的539 bp启动子区段是其活性核心区域; 同时生物信息学分析结果也显示该活性区域中包含了响应干旱和渗透胁迫的顺式作用元件。我们对转539 bp片段启动子拟南芥进行胁迫处理, 结果发现在NaCl、甘露醇和 MeJA处理后, 转基因株系中 GUS信号、GUS酶活性及GUS表达量明显降低。在BL处理后, 其活性及表达量显著增加, 而在6-BA处理后, GUS信号由根尖转移至子叶及莲座叶边缘部位且GUS活性及表达量显著升高。推测 GhGalT1基因的启动子响应盐、干旱和激素胁迫应答, 可能是通过结合干旱诱导型MYB结合元件MBS、热应答元件 HSE、防御与胁迫应答元件 TC-rich repeats、MeJA响应元件CGTCA-motif来调节的。

利用诱导性启动子在特定条件下驱动相应基因的高效表达, 增强植物对不良环境条件的抵抗性,近年来逐渐成为研究热点。目前已有一些关于植物诱导型启动子的研究报道, 如马铃薯损伤诱导型启动子Wun1[31]、柑橘冷诱导型启动子pPtcor8[32]、马铃薯光诱导型启动子Sgt1p[33]、拟南芥冷诱导型启动子CBF3[34]、拟南芥逆境诱导型启动子rd29A[35]、棉花受多逆境诱导表达的GhWRKY64基因启动子[36]、小麦干旱和渗透胁迫诱导型启动子TaPP2AbB″-α[37]等。这些能够同时受多种因素诱导表达的启动子被称为多元调控启动子, 目前能够广泛应用于转基因研究的多元调控启动子仍然比较少, 而本研究所克隆的棉花GhGalT1启动子受干旱和高盐、激素诱导表达, 属于多逆境诱导型启动子, 为了将该启动子更精确地用于棉花遗传改良, 后续的工作将对该启动子进行缺失片段分析, 确定对盐、干旱及不同激素响应的顺式作用元件及功能区, 有望进一步应用于植物遗传育种改良研究。

4 结论

从陆地棉品种Coker 312基因组DNA中克隆了GhGalT1基因的启动子 pGhGalT1, 长度为 539 bp,该区域包含响应干旱、热、脱水、防御与胁迫应答的顺式作用元件MBS、HSE、MYCCONSE、TC-rich repeats和CGTCA-motif。NaCl、甘露醇或MeJA处理后, 转基因株系中 GUS活性显著降低, BL或6-BA处理后, GUS活性显著升高, 验证了启动子活性区的功能。本研究结果为选用启动子改良作物提供了依据。

[1]Campbell J A, Davies G J, Bulone V V, Henrissat B. A classification of nucleotide-diphospho-sugar glycosyltransferases based on amino acid sequence similarities. Biochem J, 1997, 326:929–939

[2]Lao J, Oikawa A, Bromley J R, McInerney P, Suttangkakul A,Smith-Moritz A M, Plahar H, Chiu T Y, Gonzalez Fernandez-Nino S M, Ebert B, Yang F, Christiansen K M, Hansen S F,Stonebloom S, Adams P D, Ronald P C, Hillson N J, Hadi M Z,Vega-Sanchez M E, Loque D, Scheller H V, Heazlewood J L.The plant glycosyltransferase clone collection for functional genomics. Plant J, 2014, 79: 517–529

[3]Coutinho P M, Deleury E, Davies G J, Henrissat B. An evolving hierarchical family classification for glycosyltransferases. J Mol Biol, 2003, 328: 307–317

[4]Jackson R G, Kowalczyk M, Li Y, Higgins G, Ross J, Saiidberg G, Bowles D J. Overexpression of an Arabidopsis gene encoding a glucosyltransferase of indole-3-acetic acid: phenotypic characterisation of transgenic lines. Plant J, 2002, 32: 573–583

[5]Jackson R Q, Lim E K, Li Y, Kowalczyk M, Saiidberg G, Hoggett J, Ashford D A, Bowles D J. Identification and biochemical characterization of an Arabidopsis indole-3-acetic acid glucosyltransferase. J Biol Chem, 2001, 276: 4350–4356

[6]Tognetti V B, Van Aken O, Morreel K, Vandenbroucke K, van de Cotte B, De Clercq I, Chiwocha S, Fenske R, Prinsen E,Boerjan W, Genty B, Stubbs K A, Inze D, Van-Breusegem F.Perturbation of indole-3-butyric acid homeostasis by the UDP-glucosyltransferase UGT74E2 modulates Arabidopsis architecture and water stress tolerance. Plant Cell, 2010, 22:2660–2679

[7]Martin R C, Mok D W, Sniels R, Van Oiickelen H A, Mok M C.Development of transgenic tobacco harboring a zeatin O-glucosyltransferase gene from Phaseolus. In Vitro Cell Dev Biol Plant, 2001, 37: 354–360

[8]Wang J, Ma X M, Kojima M, Sakakibara H, Hou B K.N-glucosyltransferase UGT76C2 is involved in cytokinin homeostasis and cytokinin response in Arabidopsis thaliana. Plant Cell Physiol, 2011, 52: 2200–2213

[9]Wang J, Ma X M, Kojima M, Sakakibara H, Hon B K. Glucosyltransferase UGT76C1 finely modulates cytokinin responses via cytokinin N-glucosylation in Arabidopsis thahana. Plant Physiol Biochem, 2013, 65: 9–16

[10]Poppenberger B, Fujioka S, Soeno K, George G L, Vaistij F E,Hiranuma S, Seto H, Takatsuto S, Adam G, Yoshida S, Bowles D J. The UGT73C5 of Arabidopsis thaliana glucosylates brassinosteroids. Proc Natl Acad Sci USA, 2005, 102: 15253–15258

[11]Suzuki H, Fujioka S, Takatsuto S, Yokota T, Murofushi N, Sakurai A. Biosynthesis of brassinolide from teasterone via typhasterol and castaserone in cultured cells of Catharanthus roseus. J Plant Growth Regul, 1993, 13: 21–26

[12]Husar S, Berthiller F, Fujioka S, Rozhon W, Khan M, Kalaivanan F, Elias L, Higgins G S, Li Y, Schuhmacher R, Krska R, Seto H,Vaistij F E, Bowles D, Poppenberger B. Overexpression of the UGT73C6 alters brassinosteroid glucoside formation in Arabidopsis thaliana. BMC Plant Biol, 2011, 11: 51

[13]Glauser G, Boccard J, Rudaz S, Wolfender J L. Mass spectrometry-based metabolomics oriented by correlation analysis for wound-induced molecule discovery: identification of a novel jasmonate glucoside. Phytochem Anal, 2010, 21: 95–101

[14]Lim C E, Ahn J H, Lim J. Molecular genetic analysis of tandemly located glycosyltransferase genes, UGT73BI, UGT73B2, and UG17383, in Arabidopsis thaliana. J Plant Biol, 2006, 49:309–314

[15]Lim C E, Choi N J, Kim A, Lee S A, Huang Y S, Lee C H, Lim J.Improved resistance to oxidative stress by a loss-of-function mutation in the Arabidopsis UGT71C1 gene. Mol Cells, 2008, 25:368–375

[16]Kim A, Heo J O, Chang K S, Lee S A, Lee M H, Lim C E, Lim J.Overexpression and inactivation of UGT73B2 modulate tolerance to oxidative stress in Arabidopsis. J Plant Biol, 2010, 53:233–239

[17]Song J T, Koo Y J, Seo H S, Kim M C, Choi Y D, Kim J H.Overexpression of AtSGTl, an Arabidopsis salicylic acid glucosyltransferase, leads to increased susceptibility to Pseudomonas syringae. Phytochemistry, 2008, 69: 1128–1134

[18]李田, 孙景宽, 刘京涛. 植物启动子研究进展. 生物技术通报,2015, 31(2): 18–25 Li T, Sun J K, Liu J T. Research advances on plant promoter.Biotechnol Bull, 2015, 31(2): 18–25 (in Chinese with English abstract)

[19]Li F, Han Y Y, Feng Y N, Xing S C, Zhao M R, Chen Y H,Wang W. Expression of wheat expansin driven by the RD29 promoter in tobacco confers water-stress tolerance without impacting growth and development. J Biotechnol, 2013, 163:281–291

[20]Pino M T, Skinner J S, Park E J, Jeknic Z, Hayes P M,Thomashow M F, Chen T H. Use of a stress inducible promoter to drive ectopic AtCBF expression improves potato freezing tolerance while minimizing negative effects on tuber yield. Plant Biotechnol J, 2007, 5: 591–604

[21]Qin L X, Rao Y, Li L, Huang J F, Xu W L, Li X B. Cotton GalT1 encoding a putative glycosyltransferase is involved in regulation of cell wall pectin biosynthesis during plant development. PLoS One, 2013, 8: e59115

[22]Li X B, Lin C, Cheng N H, Liu J W. Molecular characterization of the cotton GhTUB1 gene that is preferentially expressed in fibers. Plant Physiol, 2002, 130: 666–674

[23]Clough S J, Bent A F. Floral dip: a simplified method for Agrobacterium mediated transformation of Arabidopsis thaliana.Plant J, 1998, 16: 735–743

[24]Qin L X, Li Y, Li D D, Xu W L, Zheng Y, Li X B. Arabidopsis drought-induced protein Di19-3 participates in plant response to drought and high salinity stresses. Plant Mol Biol, 2014, 86:609–625

[25]Jefferson R A. Assaying chimeric genes in plants: the Gus gene fusion system. Plant Mol Biol Rep, 1987, 5: 387–405

[26]Li X B, Fan X P, Wang X L, Cai L, Yang W C. The cotton Actin1 gene is functionally expressed in fibers and participates in fiber elongation. Plant Cell, 2005, 17: 859–875

[27]Wang X, Tan Y P, Zhou J, Wang C T, Liu X Q. Expression of a tobacco glycosyltransferase gene driving promoter in transgenic tobacco. Agric Sci Technol, 2010, 11: 83–85

[28]Luo K, Zhang G, Deng W, Luo F, Qiu K, Pei Y. Functional characterization of a cotton late embryogenesis-abundant D113 gene promoter in transgenic tobacco. Plant Cell Rep, 2008, 27:707–717

[29]Singh H, Sen R, Baltimore D, Sharp P A. A nuclear factor that binds to a conserved sequence motif in transcriptional control elements of immunoglobulin genes. Nature, 1986, 319: 154–158

[30]Li J J, Herskowitz I. Isolation of ORC6, a component of the yeast origin recognition complex by a one-hybrid system. Science,1993, 262: 1870–1874

[31]贾笑英, 向云, 张金文, 王蒂. 马铃薯损伤诱导型启动子Wun1基因的克隆及其GFP表达活性. 分子植物育种, 2006, 4:333–338 Jia X Y, Xiang Y, Zhang J W, Wang D. Cloning of potato wound-inducible promoter Wun1 and its GFP express activity.Mol Plant Breed, 2006, 4: 333–338 (in Chinese with English abstract)

[32]周潇, 姜航, 屈汉金, 邓子牛, Gentile A, 龙桂友. 柑橘冷诱导基因及其启动子表达载体构建与瞬时表达分析. 果树学报,2015, 32: 353–358 Zhou X, Jiang H, Qu H J, Deng Z N, Gentile A, Long G Y. Construction of plant vectors with promoter and cold-induced genes in citrus and transient expression verification. J Fruit Sci, 2015,32: 353–358 (in Chinese with English abstract)

[33]魏桂民, 张金文, 王蒂, 张俊莲, 陆艳梅, 高宜峰. 马铃薯Sgt1基因启动子的结构及功能分析. 中国生物化学与分子生物学报, 2013, 29: 969–977 Wei G M, Zhang J W, Wang D, Zhang J L, Lu Y M, Gao Y F.Promoter analysis of potato Sgt1 gene. J Biochem Mol Biol, 2013,29: 969–977 (in Chinese with English abstract)

[34]郭新勇, 程晨, 张选, 祝建波. 拟南芥冷诱导型启动子 CBF3驱动 IPT基因在烟草中的表达. 西北农业学报, 2012, 21:123–131 Guo X Y, Cheng C, Zhang X, Zhu J B. Expression of IPT gene linked with cold-induced promoter CBF3 from Arabidopsis thaliana in tobacco. Acta Agric Boreali-Occident Sin, 2012, 21:123–131 (in Chinese with English abstract)

[35]杨春霞, 陈英, 黄敏仁, 李火根. 拟南芥逆境诱导型启动子rd29A的克隆及活性检测. 南京林业大学学报(自然科学版),2008, 32: 6–10 Yang C X, Chen Y, Huang M R, Li H G. Cloning of stress-inducible promoter rd29A from Arabidopsis thaliana and its activity detection in transgenic tobacco. J Nanjing For Univ(Nat Sci Edn), 2008, 32: 6–10 (in Chinese with English abstract)

[36]杜皓, 丁林云, 何曼林, 蔡彩平, 郭旺珍. 受多逆境诱导表达的 GhWRKY64基因启动子克隆与功能分析. 作物学报, 2015,41: 593–600 Du H, Ding L Y, He M L, Cai C P, Guo W Z. Cloning and functional identification of promoter region of GhWRKY64 induced by multi-stresses in cotton (Gossypium hirsutum). Acta Agron Sin,2015, 41: 593–600 (in Chinese with English abstract)

[37]扆珩, 李昂, 刘惠民, 景蕊莲. 小麦蛋白磷酸酶 2A 基因TaPP2AbB″-α启动子的克隆及表达分析. 作物学报, 2016, 42:1282–1290 Yi H, Li A, Liu H M, Jing R L. Cloning and expression analysis of protein phosphatase 2A gene TaPP2AbB″-α promoter in wheat.Acta Agron Sin, 2016, 42: 1282–1290 (in Chinese with English abstract)