黄芩(
Scutellaria baicalensis)属大宗中药材,主要分布于温带半湿润半干旱地区,有较强的耐寒性和耐旱性
[1]。黄芩作为我国重要的药用植物资源,具有清热解毒和止血安胎等功效
[2]。至今已经从黄芩中分离出上百种化合物,黄芩苷、黄芩素等为其主要活性成分和特征性成分
[3]。然而在黄芩的生长发育过程中,常受链格孢菌(
Alternaria alternata)侵染而引发叶斑病
[4],降低叶片光合面积与效率,减少干物质积累,导致根部发育受阻及产量下降
[5]。目前,针对黄芩叶斑病的研究多集中在病原菌的鉴定层面
[4],防治手段较为基础,防治效果有限。因此,对叶斑病分子生物学方面的防御机制研究显得尤为重要。
黄酮类化合物在植物生长发育及抵御逆境胁迫中发挥重要作用
[6]。其作为重要的次生代谢产物,能有效缓解由环境胁迫诱发的活性氧损伤,维持细胞代谢稳态。黄酮类化合物的合成起始于苯丙氨酸的苯丙烷途径及随后的黄酮类特异性合成途径
[7]。其中一些关键酶基因如苯丙氨酸解氨酶(phenylalanine ammonia-lyase,PAL)、肉桂酸4-羟化酶(cinnamate4-hydroxylase,C4H)
[8]、4-香豆酸-CoA连接酶(4-coumaroyl-CoA ligase,4CL)
[9]及随后的CHS、CHI
[10]在黄酮类化合物代谢调控中起着核心作用。当
SICHS2表达下调时,番茄(
Solanum lycopersicum)中黄酮类生物合成受到显著抑制
[11];
AeCHS过表达可以提高拟南芥(
Arabidopsis thaliana)中类黄酮含量
[12];
F3'H上调增加了马铃薯(
Solanum tuberosum)中花青素和黄酮醇的含量
[13]。
绿原酸(chlorogenic acid,CGA)属于羟基肉桂酸衍生物家族,由咖啡酸和奎尼酸通过酯化反应组成
[14],是植物中含量最高的酚酸类代谢物,在杜仲(
Eucommia ulmoides)和金银花(
Lonicera japonica)中含量极为丰富
[15]。它不仅是人类膳食中重要的抗氧化剂,同时也是植物自身抵御逆境、调节生长的关键因子
[16]。植物体内绿原酸含量高低与植物防御反应之间存在密切关系,外源绿原酸可调节活性氧水平、能量代谢及苯丙烷通路,提高抗逆水平
[17],且绿原酸含量与苹果(
Malus domestica)褐变速度呈显著负相关
[18];绿原酸处理还能激活水杨酸信号通路,显著提高桃(
Prunus persica)果实对扩展青霉(
Penicillium expansum)的抗性
[19]。尽管对绿原酸在果实抗病方面的作用已有初步研究,但其通过调控次生代谢网络来增强植物抗病性的机制尚不清晰。
已有研究
[20]发现,与单作相比,黄芩与元宝枫(
Acer truncatum)复合种植能显著促进黄芩幼苗的生长,且叶绿素、脯氨酸和可溶性糖含量均显著升高,代谢物含量发生显著变化。为了进一步揭示其中的分子机制,本研究以林下经济植物黄芩为对象,研究外源绿原酸对病原菌胁迫下黄芩抗病性的影响,并利用转录组测序技术,鉴定黄酮类生物合成途径中受外源绿原酸影响的关键基因,探究绿原酸增强黄芩抗病性的分子机制,以期为林下经济植物黄芩抗病性机理研究提供新思路。
1 材料与方法
1.1 植物材料种植与处理
黄芩种子于恒温箱中37 ℃催芽,发芽后在恒温人工气候室中培养(光/暗周期16 h/8 h),35 d时对幼苗进行绿原酸处理:每株幼苗根灌100 mL 2 mmol·L-1绿原酸,对照组施加等体积的70%甲醇溶液。处理48 h后选择长势一致的黄芩幼苗36株(3个生物学重复,每个重复6棵苗)收集样品,立即放入液氮中速冻,-80 ℃保存备用。在北京海淀区(40°0′N,116°34′E)进行黄芩与元宝枫的复合种植。黄芩单植区(HQ)栽植区面积3 m×3 m,种植间距5 cm×5 cm,黄芩苗龄3个月。黄芩-元宝枫复合区(HQF)栽植区面积3 m×3 m,中心定植12年生元宝枫,黄芩种植间距5 cm×5 cm,黄芩苗龄3个月。
1.2 转录组测序分析
将经绿原酸处理后的黄芩叶片送至百迈客公司进行转录组测序,使用HISAT2软件将Clean Reads与黄芩参考基因组(GWHDEDD00000000)进行比对,获取Reads在参考基因组上的定位信息。然后利用StringTie对比对上的reads进行组装,重构转录组用于后续分析。以≥2且FDR<0.01为标准筛选差异表达基因用于后续分析。
1.3 黄芩叶片侵染方法
取大小一致的叶片,用注射器针头对叶片背面造成伤口。在链格孢菌(Alternaria alternata)菌板上加1 mL蒸馏水,用刮棒刮取孢子制成孢子液,取15 μL添加到伤口处,对照组用蒸馏水代替。28 ℃避光侵染5 d,利用ImageJ 1.54g软件计算叶片侵染面积。
1.4 叶绿素、丙二醛含量测定
参考代阳等
[20]的方法测定叶绿素、丙二醛(MDA)含量
。1.5 黄酮类化合物定量分析
参考代阳等
[20]的方法,对流动相洗脱程序进行小幅调整。采用安捷伦1260系统HPLC测定黄酮含量,条件:柱温30 ℃,流动相A为体积分数0.1%甲酸溶液,B为乙腈溶液,检测波长为280 nm,体积流速为1 mL·min
-1,进样量为20 μL。流动相洗脱程序:0~5.00 min,30%~35% B;5.00~10.00 min,35%~40% B;10.00~22.00 min,40%~50% B;22.00~25.00 min,50% B;25.00~30.00 min,50%~80% B;30.00~30.10 min,80%~30% B;30.10~35.00 min,30% B。
1.6 顺式作用元件预测
1.7 RNA提取和RT-qPCR
使用诺唯赞(北京)生物技术有限公司的RNA提取试剂盒提取黄芩总RNA,反转录合成cDNA,按照扩增程序进行实时荧光定量PCR,以2-ΔΔCT作为表达量的计算标准。内参基因引物见附表1。
1.8 代谢组学分析
将复合种植与单独种植的黄芩样品送至百迈客公司进行处理和检测,每组3个生物学重复,每个重复5棵苗,样品提取物采用超高效液相色谱-电喷雾电离-串联质谱(UPLC-ESI-MS/MS)系统进行分析,以
≥2且
P<0.05为标准进行差异代谢物的筛选,采用在线网址(
https://www.bioinformatics.com.cn/)进行热图绘制。
1.9 数据处理
采用GraphPad Prism 9软件对数据进行统计学处理。每个样本至少重复3次,采用Student’s t-test对对照组和处理组之间显著性进行分析。
2 结果与分析
2.1 黄芩与元宝枫复合种植后黄芩代谢物变化分析
采用UPLC-ESI-MS/MS技术检测了黄芩与元宝枫复合种植后黄芩代谢物变化,每个处理包括3个独立的生物学重复。结果分析显示,黄芩与元宝枫复合种植能促进黄芩有机酸类代谢物显著积累。聚类热图(
图1)显示,复合种植显著影响了黄芩中205种有机酸类代谢物,包括137种羧酸、8种苯甲酸、5种胆汁酸、6种羟基酸、2种酮酸,以及47种其他有机酸。其中107种有机酸表达上调、98种有机酸表达下调(
图1A),上调幅度较高的代谢物分别为绿原酸、(Z)-3,7-二甲基-2,6-辛二烯-1-醇乙酸酯、古液酸、2-氨基-2-脱氧异分支酸酯、4-乙基醚-3,4,5-三羟基苯甲酸(
图1B)。其中绿原酸含量高低与植物防御反应之间存在密切关系
[17],可能通过复合种植增加黄芩中绿原酸的积累,从而影响黄芩的抗病性和生长发育。
2.2 绿原酸处理对黄芩抗病性影响分析
为了探究绿原酸的积累对黄芩生长与抗病性的影响,用2 mmol·L
-1的绿原酸根灌35 d龄黄芩。正常条件下的绿原酸处理组与对照组表型无显著差异。链格孢菌侵染胁迫后,随着侵染时间增加,绿原酸处理减轻了病害对植物的损伤,主要表现为病斑面积占比的减小,对照组在侵染第5天时病斑面积占比为23.8%,而绿原酸处理组叶片的病斑面积占比显著降低至5.5%(图
2A、
2B)。与对照组相比,病原菌胁迫下,绿原酸预处理的黄芩MDA含量显著下降(
图2D),可以保护细胞膜免受损伤。施加病害胁迫后叶片叶绿素含量降低,但与对照组相比,绿原酸预处理后的植物叶片中叶绿素含量显著上升(
图2C)。以上结果表明,外源绿原酸有效减小了黄芩叶片病害面积,增强了黄芩的抗病性。
2.3 绿原酸处理下黄芩叶片的转录组测序分析
为了进一步阐明绿原酸提高黄芩抗病性的机制,对黄芩根施2 mmol·L
-1绿原酸,48 h后取样送至检测公司进行转录组测序,每组3个生物学重复,每个重复6棵苗,各组测序reads数量分别为2 015万、2 074万、2 044万、1 969万、2 023万、2 013万,总体Q20和Q30分别高于99%和96%,表明基因丰度和转录本长度足够,可用于进一步分析。以
≥2且FDR<0.01为筛选标准,共得到563个差异表达基因,其中200个基因上调表达、363个基因下调表达(图
3A、
3B)。
KEGG富集分析显示,对照组和绿原酸处理组下的DEG主要富集在与代谢相关的途径(
图3C)。包括苯丙烷生物合成、淀粉和蔗糖代谢、亚麻酸代谢、半乳糖代谢、类黄酮生物合成等代谢途径,其中苯丙烷生物合成中富集基因最多。已知苯丙烷生物合成是黄酮类物质合成的源头,因此,初步判断外源绿原酸可调节黄芩叶片中黄酮类物质代谢,该过程与黄芩抗病性密切相关。
2.4 绿原酸诱导黄芩叶片中黄酮类化合物的积累
结果显示,绿原酸处理显著影响黄芩叶片黄酮代谢途径相关基因表达。为了进一步探究绿原酸处理对黄芩中黄酮类物质代谢的具体影响,对叶片中汉黄芩苷、芹菜素、野黄芩苷、柚皮素4种黄酮类代谢物含量进行测定。高效液相色谱分析显示,绿原酸处理显著提高了叶片中柚皮素、芹菜素2种黄酮类物质含量,芹菜素质量分数增加了2.75倍,柚皮素质量分数增加了4.07倍,黄芩叶片中主要的黄酮类化合物野黄芩苷质量分数由0.96 mg·g
-1增加至1.51 mg·g
-1(
图4),这些代谢物对植物抗氧化和抗逆具有积极作用。因此,绿原酸处理可能通过调节黄酮类物质在黄芩体内的积累,特别是芹菜素和柚皮素的生物合成,来增强黄芩对病原菌胁迫的抗性。
2.5 类黄酮生物合成途径酶基因鉴定结果
绿原酸处理显著诱导了芹菜素和柚皮素的积累,为了进一步揭示其中的分子机制,集中分析了苯丙烷途径酶基因的表达变化,尤其是芹菜素、柚皮素生物合成相关基因,以FDR<0.01,
≥2为标准,共筛选到
SbPAL、
SbC4H、
Sb4CL1、
Sb4CL2、
SbFNS、
SbHCT1、
SbHCT2 7个酶基因表达上调(图
5A、
5B),均为绿原酸响应过程中的关键基因,其中
SbFNS表达上调最明显。因此,本研究认为绿原酸可能通过调控酶基因
SbPAL、
SbC4H、
Sb4CL1、
Sb4CL2、
SbFNS、
SbHCT1、
SbHCT2的表达,介导黄芩叶片中柚皮素、芹菜素的积累,从而提高黄芩的抗病性。
2.6 黄酮生物合成酶基因表达分析
为了验证转录组测序结果,RT-qPCR检测了7个酶基因的相对表达量,发现大多数基因与转录组中的表达趋势一致。外施绿原酸后,
Sb4CL1的表达量增加5.63倍,
Sb4CL2的表达量增加3.02倍,
SbHCT1的表达量增加3.50倍,
SbHCT2的表达量增加2.29倍,
SbFNS的表达量增加3.06倍,
SbPAL的表达量增加2.27倍,
SbC4H的表达量增加2.55倍(
图6)。
2.7 酶基因启动子元件分析及转录因子预测结果
基因的顺式作用元件是植物抗病相关基因的调控开关,转录因子是激活这些开关的关键调控蛋白。为了进一步研究绿原酸对代谢物合成的调节网络,利用PLANTPAN 4.0对7个酶基因的转录因子结合位点预测并用TBtools进行可视化显示(
图7A),筛选得到MYB、bHLH、WRKY、NAC、GRF 5个与抗病相关的转录因子结合位点,且在转录组数据中发现
SbWRKY1、
SbMYB、
SbGRF等基因的表达也受到绿原酸的调控。已有研究
[21]表明,转录因子可以和WD40蛋白相互作用形成复合物调控黄酮合成酶基因的表达,由此推断,绿原酸可能是通过上调转录因子的表达来提高黄酮类化合物含量,从而增加植株抗病性。
对7个黄酮合成酶基因的顺式作用元件分析表明,酶基因启动子区域富含多种与植物生长发育和抵御胁迫相关的顺式作用元件(
图7B),包括低温胁迫元件、光响应元件、防御和胁迫反应元件;同时,启动子中还含有多个与植物激素信号转导有关的调控元件,茉莉酸甲酯响应元件数量最多,
Sb4CH启动子中存在玉米素响应元件,
SbHCT2和
Sb4CL1启动子中存在赤霉素响应元件。以上结果表明,黄酮生物合成酶基因的表达受光信号、植物激素及逆境胁迫的显著影响,在植物胁迫响应中发挥关键性作用。
3 讨 论
黄芩作为重要的林下经济物种和药用植物,富含多种黄酮类化合物,却常受到根腐病和叶斑病等多种病害的侵害,生存和产量受到严重影响
[5]。当植物感染叶斑病后,叶绿素分解,进而出现光损伤、光合抑制,严重时叶片会枯萎脱落,直接导致减产并显著降低产品质量
[22]。目前有关黄芩根腐病的研究已较为系统深入,叶斑病虽田间发生普遍、危害性逐渐增强,但其致病菌近年才被相继报道,抗病机制研究相对滞后
[4]。
在植物生长发育过程中,次生代谢产物与植物的抗逆性直接相关
[6]。外源施加激素或次生代谢物对植物体内代谢物含量具有显著的诱导作用,施加ABA可以提高木豆(
Cajanus cajan)中染料木黄酮、染料木苷、对香豆酸的含量,增强木豆的抗旱性
[23]。根施绿原酸缓解了番茄的低温胁迫损伤
[24]。在本研究中,根灌绿原酸降低了真菌侵染对叶片的损伤,减小了病斑面积(
图2),黄酮类化合物柚皮素、芹菜素含量显著增加(
图4),说明黄芩的抗病性增强可能与黄酮类物质含量的增加有关,黄酮类化合物能够破坏病原菌的细胞膜结构,并且可以与植物细胞壁成分结合,增加细胞壁的厚度和强度
[25],但黄酮类化合物提高植物抗病性的分子机制还需要进一步研究。本研究分析了绿原酸处理后的转录组数据,苯丙烷生物合成途径、类黄酮生物合成途径、异黄酮生物合成途径都得到显著富集,并且柚皮素和芹菜素相关的合成酶基因表达量都显著增加(
图5)。
在本研究中,类黄酮相关合成酶基因表达上调,且
SbWRKY1、
SbMYB、
SbGRF等表达量显著上调。对7个表达上调的合成酶基因前2 000 bp启动子进行分析,发现存在多个MYB、bHLH、GRF等转录因子结合位点,因此推测绿原酸不仅是一种代谢产物,更可能作为一种外源诱导子激活植物的防御反应。大量研究证明转录因子在调控次生代谢产物积累中起着重要作用,黄酮合成酶基因的表达受到复杂的调控网络控制,番茄中SIbHLH95正调控
SIFLS的表达,抵抗灰霉病
[26];ScWRKY41直接结合
ScF3'H启动子结合区域增强黄酮类化合物从而提高马铃薯的抗寒性
[27];GmMAPK6-GmMYB4-MBW模块正向调控大豆(
Glycine max)异黄酮的合成
[28]。由此推测,外源绿原酸可能通过激活上游调控元件,诱导转录因子的转录,上调关键酶基因表达,从而提高次生代谢物含量,该推论有待后续试验进一步验证。
4 结论
本研究发现外源施加绿原酸提高黄芩的抗病性,病斑面积占比减小。黄芩中黄酮类物质含量丰富且在植物抗逆过程中发挥关键作用,高效液相色谱检测结果显示,根灌绿原酸后黄芩中柚皮素、芹菜素含量显著增加,并且转录组测序数据分析发现,差异表达基因主要富集在类黄酮生物合成相关途径,关键的黄酮合成酶基因
SbPAL、
SbC4H、
Sb4CL1、
Sb4CL2、
SbFNS、
SbHCT1、
SbHCT2表达均显著上调,随后通过RT-qPCR验证了基因的表达量,与转录组数据趋势一致。因此,本研究证明了外源绿原酸通过上调黄酮类生物合成途径中关键基因的表达,增加黄芩中黄酮类化合物含量,从而增强黄芩的抗病性(
图8),为林下经济植物黄芩抗病性机制研究提供了新思路。