第一作者联系方式: E-mail: ruiyang_liu2007@126.com **同等贡献(Contributed equally to this work)
富含油酸的菜籽油具有重要的经济价值, 使得高油酸育种及油酸形成机制成为热点。油酸脱氢酶基因( FAD2基因)是控制油酸含量的关键酶基因。本文针对 BnFAD2-C5基因展开研究, 根据油菜和甘蓝的同源性, 克隆了1257 bp启动子序列, 利用 GUS和 GFP作为报告基因分别构建含有不同片段长度的启动子和内含子的缺失载体并转化拟南芥, 经GUS染色检测发现-319 ~ -1 bp为该研究中最小启动子; 采用Western blot技术分析启动子和内含子不同区域的功能, 发现 BnFAD2-C5启动子区域-1257 ~ -1020 bp和-319 ~ -1 bp能够诱导报告基因在转基因拟南芥种子发育中期高效表达, BnFAD2-C5内含子具有增强启动子转录水平的功能, 该功能主要由+631~ +1033 bp区域调控。
High oleic rapeseed breeding and the formation mechanism of oleic acid have become a central issue after finding the important economic value of rapeseed oil with high oleic acid. The fatty acid dehydrogenase gene ( FAD2) is a key enzyme gene to control oleic acid content, but the regulation of FAD2 gene is not well understood. According to the homology between rapeseed and oleracea, the BnFAD2-C5 promoter sequence of 1257 bp was cloned. Promoter and intron of BnFAD2-C5gene were analyzed using β-glucuronidase (GUS) reporter and green fluorescent protein (GFP) reporter system to construct deleted vectors and transform Arabidopsis thaliana.Deletion analysis of BnFAD2-C5 promoter through GUS stainning revealed that -319 to -1 bp was the minimum promoter region. And deletion analysis of BnFAD2-C5 promoter and intron through GFP reporter system using western technique showed that -1257 to -1020 bp and -319 to -1 bp regions of BnFAD2-C5 promoter could induce expression of reporter genes effectively in transgenic Arabidopsisseed in the mid stage of seed development, while BnFAD2-C5intron could confer the enhancement of promoter’s function and the intron-mediated enhancement region was mainly located in +631 to +1033 bp.
菜籽油在我国食用成品油市场占有重要地位, 高油酸菜籽油不容易被氧化, 保鲜更持久[1, 2, 3], 具有重要的商业价值[4]。因此, 油酸含量的高低也成为评判油菜品质的重要指标, 同时, 高油酸育种及形成相关机理也成为目前各油料作物研究的热点[5, 6]。
目前在高油酸育种方面, 通过传统育种的手段, 官春云等[7]用60Co辐射诱变获得油酸含量高达93.5%的材料, 和江明等[8]利用EMS诱变油菜小孢子获得油酸含量为80.34%的材料。油酸脱氢酶基因(fatty acid desaturase, FAD2)作为控制油酸、亚油酸和亚麻酸含量的主要基因[9, 10], 其功能主要是将油酸磷脂胆碱去饱和成亚油酸磷脂胆碱。Stoutjesdijk等[10]通过共抑制法沉默甘蓝型油菜FAD2基因, 将油酸含量提高至89%; Peng等[11]通过RNAi技术将甘蓝型油菜FAD2基因和FAE1基因共同沉默, 使油酸含量提高至85.34%。然而, 与高油酸油菜育种相比, 油酸形成机制的研究则比较落后。详细探究脂肪酸代谢路径中的关键酶基因(FAD2基因)的调控机制, 将有助于充分利用FAD2基因为高油酸育种提供参考。2006年, Kim等[12]通过缺失分析的手段, 对芝麻FAD2基因(SeFAD2)上游调控序列分析研究, 发现在上游5′ 端非翻译区(untranslated region, UTR)内部存在一个1131 bp长度的内含子, 并且具有增强启动子功能的效应, 而启动子区域则主要含有一些起到负调控作用的顺式作用元件。同时证实脱落酸(abscisic acid, ABA)可以促进芝麻FAD2启动子的启动转录功能, 增强SeFAD2基因的表达, 并确定了ABA应答元件在启动子上的大致区域。Xiao等[13]通过对甘蓝型油菜BnFAD2-A5基因研究, 预测并分析了该启动子及内含子上的顺式作用元件, 通过检测转基因拟南芥中GUS活力, 证明BnFAD2-A5启动子属于组成型启动子, 其内含子同样具有增强基因表达的功能。BnFAD2-A5启动子的活性受ABA的诱导出现高表达现象, 并且定位到ABA在启动子上的顺式作用结合位点ABRE所处的位置。
本实验室前期已克隆到甘蓝型油菜BnFAD2基因的4个拷贝(BnFAD2-A5、BnFAD2-C5、BnFAD2-A1和BnFAD2-C1), 初步分析了BnFAD2-A5[13]和BnFAD2- C1[14]在转录水平的功能, 而BnFAD2-A1则是假基因(待发表), 因此本研究针对BnFAD2-C5基因展开研究, 克隆BnFAD2-C5启动子, 采用生物信息学分析启动子及内含子的顺式作用元件, 通过PCR技术构建启动子及内含子的缺失载体并转化拟南芥, 分析BnFAD2-C5启动子和内含子的功能, 旨在为高油酸油菜分子育种在表达水平提供理论基础。
甘蓝型油菜湘油15和拟南芥(Arabidopsis thaliana Col)由国家油料改良中心湖南分中心提供。大肠杆菌trans-T1购自北京全式金公司, pCAMBIA1303载体和农杆菌GV3101均由湖南省作物基因工程重点实验室陈信波教授赠予。
本实验室已通过RACE技术克隆了BnFAD2-C5基因的全长序列, 确定转录起始位点及转录终止位点, 且发现在5′ UTR内部含有一个内含子序列(图1)。在此基础上, 本研究根据油菜与甘蓝的同源性, 依据甘蓝数据库信息(http://ocri-genomics.org/bolbase), 在转录起始位点上游约1700 bp处设计正向引物PBn-F1 (5′ -AAATGAAATGAAAATCATGGTAGGTG -3′ ), 在该基因编码区内部设计反向引物CDS-R (5′ - GCTTGATGTTGTGTCGGTTTCAGACTT-3′ ), 通过PCR技术从甘蓝型油菜基因组DNA中扩增BnFAD2- C5启动子区域, 然后转化到大肠杆菌并测序。
利用PLACE (http://www.dna.affrc.go.jp/htdocs/ PLACE)和PlantCARE (http://www.intra.psb.ugent.be: 8080/PlantCARE)网站对BnFAD2-C5启动子及内含子序列预测分析。
缺失载体构建过程中所用引物见表1。采用PBn-F1和IBn-r引物从基因组中克隆启动子+内含子片段, 经Hind III/Nco I双酶切处理, 连接到pCAMBIA1303载体, 构建表达载体Bn-1, 导入大肠杆菌保存。以含有Bn-1载体的大肠杆菌菌液为模板, 分别以引物PBn-F1/PBn-R、PBn-F2/PBn-R、PBn-F3/PBn-R、PBn-F4/PBn-R、PBn-F5/PBn-R克隆启动子缺失片段, 经Hind III/Nco I双酶切处理后, 按照上述方法, 将其连接到pCAMBIA1303载体, 构建启动子缺失表达载体BnP-1、BnP-2、BnP-3、BnP-4和BnP-5 (图2)。
同样以含有Bn-1载体的大肠杆菌菌液为模板, 以IBn-F/IBn-r引物进行扩增; 以IBn-R1、IBn-R2、IBn-R3和IBn-R4为反向引物, IBn-F1为正向引物, 进行PCR扩增, 各取3 μ L PCR产物与3 μ L IBn-F/IBn-r产物进行融合PCR, 反应体系含缓冲液4 μ L、dNTPs 2 μ L、HS Taq聚合酶0.5 μ L、ddH2O 7.5 μ L; 反应条件为95℃ 5 min; 95℃ 50 s, 58℃ 50 s, 72℃ 4 min, 15个循环。取2 μ L融合产物作为模板, 以IBn-F1/IBn-r为引物进行PCR扩增, 扩增程序为95℃ 5 min; 95℃ 50 s, 58℃ 50 s, 72℃ 4 min, 25个循环; 72℃ 10 min。将融合产物克隆到大肠杆菌并测序。选择测序正确的菌株摇菌, 提取质粒并进行Hind III/Nco I双酶切, 然后将酶切后的产物纯化并连接到经Hind III/Nco I双酶切处理的pCAMBIA1303载体, 构建内含子缺失载体BnI-1、BnI-2、BnI-3和BnI-4 (图2)。以含有CaMV35S启动子的pCAMBIA1303载体为阳性对照, 去除CaMV35S启动子的作为阴性对照。
将构建好的表达载体通过冻融法导入农杆菌GV3101, 将转化后的农杆菌菌液涂到含有庆大霉素(GM)和卡那霉素(Kan)的YEB固体平板筛选, 挑选单菌落进行PCR检测, 并将阳性菌液保存于-80℃。在转化拟南芥前, 将含有目的表达载体的农杆菌画线接种于含有50 mg L-1GM+50 mg L-1 Kan的YEB固体平板, 28℃倒置培养2 d, 挑取单菌落于10 mL含有50 mg L-1 GM+50 mg L-1 Kan的YEB液体培养基摇过夜, 取2 mL接种于50 mL含有50 mg L-1 GM+50 mg L-1 Kan的YEB液体培养基, 待其OD600为0.4~0.6时, 收集菌体(5000转 min-1, 5 min), 用5 mL 10 mmol L-1 MgSO4溶液重悬, 加入5 μ L乙酰丁香酮(200 μ mol L-1)于暗光下静置1.5 h。随后加入100 mL渗透培养基, 通过花序浸染法转化拟南芥。
参照Jefferson等[15]GUS组织化学检测的方法检测转基因拟南芥种子不同发育时期(4、7、10和13 d) GUS表达情况。以GFP为报告基因, 以Actin为内参, 通过Western blot技术对GFP和Actin显影, 利用Quality One软件将Western杂交斑点数字化, 以GFP含量比Actin含量所得的数值表示GFP蛋白的相对含量, 以此确定不同缺失载体启动报告基因GFP的功能大小。
根据油菜和甘蓝的同源性, 同源克隆了1257 bp启动子序列, 与甘蓝基因组序列仅有2个碱基的差异。利用PLACE和PlantCARE网站对BnFAD2-C5基因启动子和内含子区域预测分析(图3), 发现在启动子区域存在核心启动子必备的元件TATA box (-139 bp)和CAAT box (-61 bp), 同时还存在一些控制基因高效转录的元件(CCAATT box, GATA, 5′ UTR Py-rich), 激素应答元件(TCA motif、CGTCA motif、ABRE、MYB、MYC和GARE), 根毛特异性元件(root hair specific element, RHE), 光调控元件(GAG, G-boxGT-1 consensus合Dof), 花粉特异性表达元件(GTGA), 胚乳特异性表达元件(ACGT), 以及具有正、负调控功能的元件CACA。
在内含子区域, 同样存在TATA box和CAAT box, 分别位于+1241 bp和+1040 bp处。另外, 内含子区域也包含一些顺式作用元件, 如胚乳特异性表达原件(ACGT, Skn-1motif), 激素应答原件(MYB、MYC、GARE), 分生组织特异性原件(CAT box), 抗性原件(W-box), 光控原件(G-box、GT-1consensus), 热激蛋白原件(heatshockprotein element, HSE)以及Dof和E-box。
由GUS染色结果(图4)可知, 当BnFAD2-C5启动子缺失掉大部分区域, 仅保留-136 ~ +137 bp区域时(BnP-5), GUS染色结果与阴性对照相同, 该部分序列不能发挥启动子功能, 而BnP-1、BnP-2、BnP-3和BnP-4均有不同程度的着色。将含有BnP-1、BnP-2、BnP-3和BnP-4的转基因拟南芥发育中的种子进行蛋白相对定量检测, 将Western blot结果经Quality软件数字化后, 以GFP含量比Actin含量所得的数值表示GFP蛋白的相对含量, 结果如图5。BnP-1和BnP-4的GFP蛋白相对含量高于BnP-2和BnP-3, 而BnP-2和BnP-3的相对蛋白含量较为接近, 说明在缺失的-1257 ~ -1020 bp区间存在调控基因高表达的正调控元件, 在-581 ~ -319 bp区间存在抑制基因表达的负调控元件, 且主要是抑制-319 ~ -1 bp区间的顺式作用元件的功能发挥。BnP-1和BnP-4在种子发育中期出现表达高峰, 而BnP-2和BnP-3在整个种子发育过程中表达量变化较为平缓, 说明在-1257~ -1020 bp和-319~ -1 bp区间存在某些调控元件诱导报告基因GFP在种子发育中期高效表达。
由图5-A中BnP-1转基因拟南芥种子GFP蛋白相对含量与图6-A中Bn-1比较可知, BnFAD2-C5内含子具有增强启动子功能的效应, 为了研究发挥这一增强效应的内含子区域, 本研究将内含子区间缺失并转化拟南芥验证。分别采用GUS染色和Western blot方法检测转基因拟南芥不同发育时期(4、7、10和13 d)种子中报告基因GUS/GFP的含量。GUS染色结果显示, 随着内含子片段的逐步缺失, 染色有不同程度的变化(图7), 说明内含子不同区域发挥着一定的功能。利用Quality One软件将Western杂交斑点数字化后绘制曲线(图6), 发现当缺失+631~ +1033 bp区域后, GFP蛋白相对含量大幅下降, 当继续缺失+526 ~ +631 bp时, GFP相对蛋白含量几乎降至零, 但下降幅度减小, 说明在+631 ~ +1033 bp区域存在正调控元件促进基因的高效表达, 在+526 ~ +631 bp区域也存在正调控元件, 但功能较弱。随着进一步的缺失, BnI-3和BnI-4 GFP蛋白相对含量逐渐上升, 说明在缺失掉的+266 ~ +526 bp区间, 含有负调控元件抑制启动子的功能。因此, 推测内含子的增强效应主要取决于+631 ~ +1033 bp区域顺式作用元件的功能发挥。
BnFAD2-C5基因的全长序列已克隆, 发现长度
为1123 bp的内含子位于5′ UTR内部。类似的序列特征在玉米[16]、大豆[17]、拟南芥[12]、芝麻[12]和油菜[13]中均有报道, 且发现这些内含子通常具有增强效应。本文也证实该内含子具有增强BnFAD2-C5基因表达水平的功能。一般而言, 内含子增强效应由内含子剪切模式[18]、内含子的长度[19]以及内含子的特殊元件决定, 但是Parra等[20]认为内含子剪切模式不足以对启动子功能产生增强效应, Carola等[21]证实内含子的长度越长其增强效应越强。在本研究中, 缺失掉大部分内含子序列的BnI-4其GFP蛋白相对含量却比BnI-3和BnI-4有所增加, 那么, 认为内含子的长度并不是导致该内含子发挥增强效应的原因, 该内含子的增强效应极有可能是特殊的调控元件与启动子区域互作引起的。当内含子+631~ +1033 bp区域缺失后, GFP蛋白相对含量明显下降, 而随着进一步缺失并没有再次引起GFP蛋白相对含量的大幅下降(图5-A), 说明内含子增强效应区域主要存在于+631 ~ +1033 bp区域, 而不是整个内含子区间。Parra等[20]认为, 内含子的增强效应元件主要存在于内含子的5′ 端, 而不是内含子的3′ 端, 在5′ UTR和编码区也会大量出现。如果内含子含有大量CGATT拷贝, 则会大幅增加基因的表达量。在油菜BnFAD2- C5内含子的5′ 端+262 ~ +266 bp处存在一个CGATT, 并且BnI-4 GFP蛋白相对含量高于BnI-1、BnI-2和BnI-3, 所以当+266 ~ +631 bp区间的负调控元件被解除后, CGATT元件极有可能增强了基因的表达水平, 导致BnI-4中GFP蛋白相对含量增加。Kim等[12]发现芝麻FAD2基因的内含子增强效应区域为+791 ~ +1517 bp, 处于远离内含子5′ 端的位置, 并推测内含子中具有增强功能的元件与增强子的作用方式相近, 并不受距离的影响, 可以在远距离发挥增强作用。本研究中BnFAD2-C5内含子增强效应的主效元件位于+631 ~ +1033 bp区域, 也远离内含子5′ 端, 推测该内含子功能的发挥不受距离的影响, 发挥功能的元件可以远距离增强启动子的转录水平。
Rose[22]研究发现, 内含子增强效应主要是富含U的序列导致, 内含子的5′ 和3′ 剪切位点和分支点突变时, 会影响内含子的剪切, 但不一定会引起内含子增强效应的减弱。我们分析了+631 ~ +1033 bp区域, 其中存在38.7%的U序列, 同时还存在ABRE/MYB/MYC、E-box、GARE元件, 研究发现这些元件在油菜种子napA启动子和拟南芥中, 当植物受到干旱刺激和病原菌危害时能够增强相应启动子的功能, 诱导基因的高表达[23, 24]。
从图4-C可知, 含有BnP-1和BnP-4的转基因拟南芥种子发育至7 d时, 出现高表达现象, 具有典型诱导表达的特点, BnP-2和BnP-3则呈现组成型表达的特点, 在种子发育的各个时期蛋白含量变化微弱。因此我们认为-1020 ~ -319区域存在控制基因组成型表达的相关元件。在BnP-1的-1257 ~ -1020 bp和BnP-4的-319 ~ -1 bp中存在某些元件在种子发育至7 d时可以被激活, 诱导该基因的表达, 通过PLACE和PlantCARE网站预测可知, 在-1257~ -1020 bp和-319~ -1 bp分别含有2 Dof, 3 GT-1consensus、BoxII、TCA motif、GAG和GATA、G-box、3GT-1consensus、4 Dof、GAG、GC-motif、CGTCA motif、TATA-box和CAAT-box。GT-1consensus在病原菌和NaCl的诱导下会促进SCaM-4基因启动子转录[25], TCA motif是水杨酸应答元件, 抑制iMyAP基因启动子的转录[26]。Dof、BoxII、GAG是光诱导调控元件, 感受光刺激诱导基因的表达, Dof蛋白结合Dof调控位点增强cyPPDK基因和PEPC基因的转录[27]。CGTCA motif是茉莉酸应答元件, 当植物体内茉莉酸含量上升时, 诱导相关基因的表达[26, 28]。GATA元件, 具有增强cab基因启动子转录水平的功能[29]。CAAT-box能提高转录起始频率, 并且参与基因表达量的调控[30]。因此, 下一步工作可以针对-1257 ~ -1020 bp和-319 ~ -1 bp区域展开研究, 探索对BnFAD2-C5基因发挥诱导功能的主效区域及主效元件。
克隆了BnFAD2-C5启动子序列, 其-1257~ -1020 bp和-319~ -1 bp区域具有诱导基因表达的功能; BnFAD2-C5内含子具有增强转录水平的效应, 这一功能主要由+631~ +1033 bp区域调控。
The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|