橡胶树白粉菌（HO-73）启动子WY172 不同长度片段的克隆及表达活性分析

殷金瑶 王义 徐良向 朱利 王晨 刘文波 缪卫国

（海南大学植物保护学院 热带农林生物灾害绿色防控教育部重点实验室，海口 570228）

启动子是一段位于结构基因5′端上游的DNA序列，提供RNA 聚合酶识别和结合位点，是调控基因表达强弱的重要元件，通过参与基因的转录环节来调控基因表达的水平、部位及方式［1-2］。研究表明，通过自身启动子在植物中对外源基因进行功能表达，为植物生物技术中有益转基因的应用提供了可能，如提高植物胁迫耐受性等［3］。

目前，最经典的启动子是花椰菜花叶病毒（CaMV）35S 启动子［4］、玉米泛素（Ubiquitin）启动子［5］和水稻肌动蛋白（ACT1）启动子［6］等。但随着应用不断深入，上述几个启动子的缺陷也不断暴露出来，如因高效表达需要产生的大量能耗会对植物生长产生影响，引起植物生长缓慢［7］。同时因其表达具有持续性，不随时间和空间的变化而变化，在植物体内表达不具有特异性［8］，因而人们开始寻求更加高效、优质的启动子。

有研究表明，使用丝状真菌中的强启动子可以高效启动目的基因的转录。在丝状真菌中利用自身启动子起始转录外源基因的转录效率通常很高，但是一些外源蛋白的产量远远低于内源蛋白的产量［9］。有文献报道，启动子活性、翻译效率、蛋白糖基化、目标基因的拷贝数以及菌株自身蛋白酶的降解能力等均会影响蛋白的产量，其中以启动子转录水平的高低影响最为关键，使用强启动子可以实现目的基因的高效表达［10-11］。目前已报道的丝状真菌启动子主要是曲霉葡糖淀粉酶基因（glaA）启动子，构巢曲霉的3-磷酸甘油醛脱氢酶基因（gpdA）启动子，米曲霉的a-淀酶基因（amyB）启动子，瑞氏木霉纤维二糖水解酶基因（cbh1）启动子［12］。

橡胶树白粉菌（Oidium heveae）引起的橡胶树白粉病是橡胶生产中重要病害之一［13］。在前期研究中，我们首次报道了橡胶树白粉菌基因组［14］，构建了适于快速检测橡胶树白粉菌菌株HO-73 启动子的探针载体［15］，并通过将启动子与CaMV35S 启动子及ACT1 启动子的活性比较，初步鉴定了WY7［16］、WY51［17］、WY193［18］等能够用于调控双子叶植物和单子叶植物中的外源目的基因表达的橡胶树白粉菌内源启动子。本研究中通过对其基因组进行分析，预测筛选了一些疑似内源启动子的核心区域，并通过克隆得到了其中1 个疑似新的内源启动子WY172（BankIt2158599，另文发表）。本研究以WY172 上游2K 序列为研究对象，通过利用生物信息学软件对其进行检索并预测转录因子结合位点，通过构建不同长度疑似具有启动子活性片段的表达载体，观察报告基因的表达情况，探讨WY172 及4 个不同长度缺失片段的表达活性，旨在寻求更好的可服务于基因工程、应用于生产实际的优秀启动子。

1 材料与方法

1.1 材料

橡胶树白粉菌菌株HO-73 筛选、保存于本实验室，培养于橡胶苗中感品种热研7-33-97［19］上，培养室内温度为25℃恒温，光周期为14 h/10 h（光/暗），湿度范围88%左右。模式生物三生烟（Nicotiana tabacumcv.Xanthi nc）作为研究对象，培养于上述培养室中。

1.2 方法

1.2.1 WY172 启动子生物信息学分析 Promoter Scan 预测橡胶树白粉菌全基因组DNA，得到若干橡胶树白粉菌疑似内源启动子的核心区域，将其中一个命名为WY172（BankIt2158599），并利用在线数据库EPD 确认所验证的启动子是否为未报道的新启动子。利用在线数据库PROMO 对WY172 上游2K序列进行转录因子预测分析，利用在线生物信息学软件JASPAR 检索相应的转录因子结合位点，根据转录因子结合位点在序列上的分布位置对WY172上游2K 序列进行渐变缺失突变，用引物设计软件Primer Premier 5 设计5′端不同缺失的上游引物和3′端下游引物。同时通过plantCARE 软件对WY172 启动子上游2K 序列的作用元件进行分析。

1.2.2 WY172 启动子及不同长度片段的克隆 提取橡胶树白粉菌基因组DNA（OMEGA，D3390），以此为模板，通过上游特异引物172F、172Q1F、172Q2F、172Q3F、172QF 和共同的下游引物172NR进行PCR 扩增，获得5 种不同长度的启动子片段。

将PCR 产物经1.0%的琼脂糖凝胶电泳检测并切胶回收。将回收产物与T-easy 载体连接，进行TA克隆，并转化Trans5α 感受态细胞，涂布于Amp 抗性的LB固体培养基（100 μg/mL），37℃过夜培养16 h。挑取单菌落于Amp 抗性的LB 液体培养基中，37℃摇床震荡培养12 h 后进行菌液PCR 验证，验证正确的阳性克隆菌株送华大生物公司测序。测序正确的菌株用80%的甘油保菌。

1.2.3 WY172 启动子及不同长度片段瞬时表达载体的构建 选择HindIII/BamH I 分别酶切并回收pBI121 载体和1.2.2 中的含不同长度片段的克隆载体，进行连接，并转化Trans5α 感受态细胞，涂布于Kan 抗性的LB 固体培养基（100 μg/mL），37℃培养16 h。挑取单菌落于Kan 抗性的LB 液体培养基中，37℃摇床震荡培过夜后进行菌液PCR 验证，正确的阳性克隆菌株送华大生物公司测序，测序正确的菌株用80%的甘油保菌。

1.2.4 三亲杂交法转化农杆菌菌株 挑选含有目的基因的大肠杆菌菌株，在含有50 μg/mL 卡那霉素的LB 液体培养基中37℃震荡培养（带有CaMV35S 启动子的pBI121 载体作为阳性对照），至OD600nm为0.5时，与LBA4404、HB101 感受态等体积混匀，涂布于不含任何抗生素的LB 固体培养基上，28℃过夜培养。用接种针将长出的菌落转移到含有100 μg/mL利福霉素和50 μg/mL 卡那霉素的固体LB 培养基上，28℃培养3-4 d，长出单菌落。将长出的单菌落再次转移到含有100 μg/mL 利福霉素和50 μg/mL 卡那霉素的LB 固体培养基上，挑选单斑在含有50 μg/mL卡那霉素和100 μg/mL 利福霉素的LB 液体培养基中37℃震荡培养，进行菌液PCR 检测，正确的阳性克隆菌株送华大生物公司测序，测序正确的菌株用80%的甘油保菌。

1.2.5 GUS 染色 利用农杆菌介导的转化法（ATMT）侵染三生烟，将侵染后的烟草叶盘28℃共培养2 d后取样，每个片段启动子取6 片叶盘浸泡于GUS 染色液［20］中，37℃过夜培养，后将样品转入70%乙醇中脱色至叶片组织无色，观察。

1.2.6 GUS 酶活性测定 取各样品0.1 g 液氮研磨，提取植物总蛋白（CWBIO，Plant Protein Extraction Reagent），-80℃保存。GUS 蛋白浓度的测定采用BCA 法（Solarbio，BCA Protein Assay），使用多功能酶标仪Infinie 200 PRO 测定样品在562 nm 处的吸光值，最后根据标准曲线计算样品蛋白浓度。采用Jefferson 等［21］报道的方法测定GUS 蛋白活性，结果以生成的4-甲基伞形酮（4-methylumbelliferon，4-MU）的量与总蛋白的量和时间的比值（pmol 4-MUG/mg protein/min）来表示。实验设置3 次重复。最后使用SPSS 软件统计结果。

2 结果

2.1 WY172启动子生物信息学分析

在线数据库EPD 确认未发现与WY172 同源的启动子。对转录因子的预测分析结果显示WY172启动子上游2K 序列上共有转录因子结合位点17 个（表1），运用IBS 软件绘制基因图谱（图1），直观地展示出各个转录因子结合位点在基因上的位置，为后续缺失突变位置的确定提供参考。根据转录因子结合位点在序列上的分布位置，以不切断转录因子结合位点为原则，设计了4 个缺失突变片段，分别 为WY172Q（2250 bp）、WY172Q1（1551 bp）、WY172Q2（1281 bp） 和WY172Q3（871 bp）（ 图2）。结合前期研究得到的WY172 启动子（251 bp），通过Primer Premier 5 软件设计了6 个引物（表2）。PlantCARE 软件分析结果（表3）显示，其上游2K序列含有丰富的作用元件，包括光响应元件、脱落酸响应元件、生长素响应元件、茉莉酸甲酯响应元件、创伤诱导元件、厌氧感应元件、参与栅栏叶肉细胞分化的元件及玉米蛋白代谢调节的顺式作用元件等。

2.2 WY172启动子及不同长度片段的克隆

通过PCR 及产物胶回收，成功获得5 个带有酶切位点的不同长度片段，并构建了Teasy-WY172、Teasy-WY172Q、Teasy-WY172Q1、Teasy-WY172Q2、Teasy-WY172Q3 共5 个载体，菌液PCR 结果显示出与目的片段大小基本一致的条带（WY172Q 为2300 bp 左右，WY172Q1 为1500 bp 左 右，WY172Q2 为1300 bp 左 右，WY172Q3 为800 bp 左 右，WY172为250 bp 左右）。并经过华大生物公司测序及NCBI比对结果显示序列相似性为100%，说明成功克隆了WY172 启动子及不同长度片段。

表1 WY172 启动子上游2K 序列含有的转录因子结合位点

图1 WY172 启动子上游2K 序列转录因子分布图

表2 引物序列

2.3 WY172启动子及不同长度片段瞬时表达载体的构建

阳性克隆菌株的菌液PCR 结果显示获得与目的条带大小一致的片段，华大公司测序结果经过比对显示DNA 相似性为100%，获得了重组的植物表达 载 体pBI121-WY172、pBI121-WY172Q、pBI121-WY172Q1、pBI121-WY172Q2、pBI121-WY172Q3（以pBI121 载体作为阳性对照），同时通过Vector NT1绘制了植物表达载体示意图（图2-A），并利用IBS软件绘制了WY172 启动子不同长度片段特征示意图（图2-C）。将5 个重组植物表达载体通过三亲杂交法转化农杆菌，菌液PCR 检测结果显示出与目的基因大小一致的特异片段（图2-B）（WY172Q 为2300 bp 左右，WY172Q1 为1500 bp 左右，WY172Q2 为1300 bp 左 右，WY172Q3 为800 bp 左 右，WY172为250 bp 左右），经过华大生物公司测序及NCBI 比对结果显示序列相似性为100%，说明获得了含有各缺失片段启动子的农杆菌菌株。

表3 WY172 启动子上游2K 序列含有的作用元件

图2 不同长度启动子片段及载体构建

2.4 GUS蛋白染色

GUS 染色结果（图3）显示，所有菌株瞬时表达的叶盘都成功染出了蓝色，其中阳性对照菌株（pBI121 载体）侵染的叶盘染色最浅，而WY172 不同片段菌株侵染的叶盘染色深度均强于阳性对照。pBI121-WY172Q3 叶盘染色最深，颜色由深至浅排列是：pBI121-WY172Q3＞ pBI121-WY172Q2＞ pBI121-WY172Q1＞ pBI121-WY172＞ pBI121-WY172Q，总体呈现随长度增长GUS 染色颜色逐渐变浅的趋势，但WY172Q3 片段出现蓝色显著加深现象。结果表明，WY172 启动子的5 个不同长度片段均能启动GUS基因的表达，其中WY172Q3 的表达活性最强。

图3 不同长度启动子片段烟草瞬时表达的GUS 染色

2.5 GUS酶活性检测

为测定不同长度片段的酶活性，绘制了BSA 标准曲线和4-MU 标准曲线，R2值分别为0.9995 和0.9998，达到了测定要求（图4）。SPSS 软件统计酶活性分析结果（图5）显示，野生型几乎检测不到酶活性，而CaMV35S 启动子和5 个WY172 不同长度片段启动子均具有驱动GUS基因表达的酶活性，GUS 酶活性趋势为：WY172Q3＞ WY172Q2＞WY172Q1＞ WY172＞ WY172Q＞ CaMV35S，总体呈现随长度的增长活性逐渐递减趋势，但在WY172Q3处酶活性出现显著回升。酶活性检测结果与GUS 染色结果一致，再次证明WY172 启动子的5 个不同长度片段均具有启动子活性，WY172Q3 表达活性最强。

图4 BSA 标准曲线和4-MU 标准曲线

图5 不同长度启动子片段烟草瞬时表达的GUS 酶活性

3 讨论

通过对不同长度启动子的表达活性进行分析可以得知启动子发挥作用的能力［22-23］，为充分解读其对基因的调控表达过程及后续应用奠定基础。橡胶树白粉菌是一种专性寄生菌，分子研究严重落后，启动子方面的研究更是少之又少［24］。目前，已见报道的研究多集中于启动子的筛选和鉴定。本实验室前期通过不同卡那霉素浓度下启动子的耐受性大小来筛选强启动子，并以橡胶树白粉菌基因组为研究对象，鉴定出多个在单子叶和双子叶植物中均能表达的强启动子。在此研究的基础上，本研究聚焦于WY172 启动子不同长度片段驱动外源基因表达的能力，以WY172 启动子上游2K 序列为研究对象，进行渐变缺失突变，分析不同片段驱动GUS基因的表达活性。

本研究经过初步分析发现WY172 启动子上游2K 序列含有丰富的顺式作用元件，其中有参与创伤反应的WUN-motif 作用元件，厌氧感应的ARE 顺式作用元件，参与叶肉细胞栅栏组织分化的HD-Zip1元件，参与玉米蛋白代谢调节的O2-site 顺式作用元件，以及参与昼夜控制的circadian 顺式作用元件等功能作用元件，还有茉莉酸甲酯响应、赤霉素响应、生长素响应等激素作用元件，此外还有低温响应、光响应元件等非生物因素的作用元件。这表明WY172 启动子很可能在植物生长的各个阶段都能发挥作用，多样性的元件也为下一步高效表达外源基因的启动子的开发提供了可能，为WY172 启动子的进一步研究和开发奠定了基础。

GUS 染色结果及酶活性检测表明，在调控GUS基因表达方面，WY172 启动子的不同长度片段都具有比CaMV35S 启动子强的启动子活性。CaMV35S启动子是目前国际上应用最广、活性最强的启动子之一，但随着研究的深入，该类启动子在应用上的不足已经显现地越来越明显［25］。WY172 启动子的不同长度片段均具有比CaMV35S 启动子更强的活性，这为基因工程工具的开发提供另外一个可能性。在后续实验中将会对WY172 启动子的具体类型及应用范围等进行分析，以便为基因工程开发一个更有效的启动子工具。

顺式作用元件的种类和数量的不同，以及转录因子结合位点的差异，都有可能影响基因的表达效果。活性分析结果显示，WY172 不同片段的活性强弱关系为WY172Q3＞ WY172Q2＞ WY172Q1＞WY172＞ WY172Q＞ CaMV35S。WY172Q3 段 启 动 子表达活性升高可能是由于在251 bp-821 bp 区段含有转录激活元件，使得该区段具有显著增强效应；WY172Q2 段启动子区段也有增强效应但表达活性低于251 bp-821 bp 段的原因，一个可能是821 bp-1281 bp 区段含有的激活元件种类和数量都少于WY172Q3 段，也可能是821 bp-1281 bp 区段同时含有抑制表达元件，与251 bp-821 bp 区段中的一部分激活元件的功能相抵了；WY172Q1 段启动子及WY172Q 段启动子表达活性降低，可能是由于1281 bp-2250 bp 区段含有抑制表达元件，导致了抑制效应。未来，我们将会对具体的抑制、激活元件进行定位及探究，从而研究该启动子不同情况下的表达量变化。

我们的启动子来源于丝状真菌，这类真菌表达外源蛋白的量比较低［26］，为实现外源基因的高效表达，通常需要强大的启动子来驱动基因的转录，但目前针对丝状真菌启动子的研究与应用比较少，主要是木霉属（Trichoderma）和黑曲霉属（Aspergillus niger）的一些启动子［27］。除此之外，丝状真菌启动子大部分都只能应用于自身菌体中外源基因的表达，而来源于橡胶树白粉菌的WY172 启动子不仅具有丝状真菌启动子的优点，还可以应用于植物，这为其推广和应用提供了更大的可能。为后续对于其下游基因和调控方式的确定提供方向、奠定基础，同时从一定程度上也促进了橡胶树白粉菌的分子系统发育和形态学分析研究工作。

4 结论

本研究克隆了WY172 启动子及其4 个不同长度的缺失片段，并构建了植物表达载体，通过GUS染色和酶活性检测，对启动子表达活性进行了定性和定量分析。证明以上WY172 启动子5 个不同长度片段均具有启动子活性，均能驱动GUS基因的表达。GUS 染色结果显示WY172 启动子不同长度片段的蓝色强度均优于阳性对照CaMV35S 启动子，其中WY172Q3 驱动GUS基因表达的酶活性最高。且酶活性检测也与GUS 染色结果相符。不同片段的缺失对表达活性的影响不同，pBI121-WY172Q3 启动子表达活性最强。