1、类黄酮生物合成调控网络中转录因子的筛选鉴定
本课题组前期已鉴定获得‘红阳’的全红果肉突变体 “H16”,并发现二者的类黄酮与花青素含量存在显著差异(Ye 等,2025)。本研究对果实发育全过程的类黄酮变化进行了重新定量分析,结果显示,在果实发育的 3 个时期中,“H16” 的总黄酮与总花青素含量均显著高于‘红阳’(图1,A-C)。为揭示其潜在分子机制,本研究基于课题组前期获得的‘红阳’与 “H16” 转录组数据(NCBI SRA 数据库,登录号 PRJNA1032890)(Zhang 等,2023),开展了加权基因共表达网络分析(WGCNA)(图S1,A)。层次聚类树状图结果显示,差异表达基因(DEGs)可聚为19个模块,其中棕色、绿色、青色与粉色模块与总花青素及两种花青素组分的含量呈高度相关(图1,D)。GO与KEGG富集分析表明,这些模块中的差异基因主要富集于果实发育、次生代谢物生物合成相关通路,其中类黄酮生物合成与苯丙烷生物合成通路的富集程度最为显著(图S1,B 和 C)。基于模块隶属度与模块内连接度较高的基因,本研究构建了蛋白质 - 蛋白质相互作用(PPI)网络(图1,E),并从中筛选得到一个参与类黄酮生物合成调控的潜在关键因子AcbHLH137,后续将对其开展功能验证。
图1 基于WGCNA的猕猴桃果肉花青素积累相关关键转录因子的鉴定
2、AcbHLH137的表达与花青素含量呈正相关
为进一步探究AcbHLH137在类黄酮生物合成中的功能,本研究检测了其在果实发育过程中的表达模式。qRT-PCR结果显示,AcbHLH137在 “H16” 中持续上调表达,且在果实成熟期的表达量显著高于‘红阳’,该结果与转录组测序数据一致(图2,A 和 B),提示 AcbHLH137可能正向调控类黄酮的生物合成。相关性分析结果表明,AcbHLH137的表达量与类黄酮、花青素含量均呈显著正相关,但与原花青素含量无显著相关性(图 2,C)。
图2 AcbHLH137的分子特征
本研究分别从‘红阳’与 “H16” 中克隆获得AcbHLH137的全长编码序列(CDS),二者序列完全一致,开放阅读框(ORF)长度为969 bp。蛋白序列比对结果显示,AcbHLH137与其他物种同源基因的N端序列相似性较低,但其bHLH结构域与C端区域具有高度保守性(图S2)。系统发育分析表明,AcbHLH137在进化上与bHLH XII亚族成员(如茶树 CsbHLH137、拟南芥 AtBPE)亲缘关系较近(图2,D)。亚细胞定位实验结果显示,AcbHLH137定位于细胞核与内质网中(图2,E)。
3、过表达AcbHLH137促进转基因猕猴桃植株及愈伤组织中花青素与原花青素的积累
为进一步明确AcbHLH137在类黄酮代谢中的调控作用,本研究通过农杆菌介导的遗传转化获得了 2 个AcbHLH137过表达转基因猕猴桃株系 OE-8 与 OE-19(图S3)。与野生型(WT)相比,OE-8 与 OE-19 的叶片呈现出明显的黄化与红化表型(图3,A)。转基因植株中AcbHLH137的表达量分别为野生型的 2.4 倍与 4.1 倍(图3,B)。与野生型相比,OE-8与OE-19的花青素与原花青素含量均显著升高(图3,C 和 D)。RT-qPCR 结果显示,转基因植株中类黄酮途径的多个结构基因表达显著上调,包括UFGT、ANR、F3H、F3'H、F3'5'H、CHI、FLS(图3,E)。
此外,本研究还在猕猴桃愈伤组织中过表达了AcbHLH137(图3,F;图 S4,A)。转基因愈伤组织(OE-5)的原花青素含量显著高于野生型,DMACA 染色颜色更深(图3,G 和 H)。qRT-PCR 分析表明,除AcLAR外,AcbHLH137过表达显著上调了转基因愈伤组织中类黄酮途径绝大多数结构基因的表达(图S4,B)。值得注意的是,AcANR在过表达植株(OE-8、OE-19)与愈伤组织(OE-5)中均显著上调,而AcLAR在转基因愈伤组织(OE-5)中则显著下调(图 S4,B)。
图3 AcbHLH137过表达促进花青素和原花青素的积累
4、AcbHLH137直接结合类黄酮生物合成基因的启动子
为进一步探究AcbHLH137是否直接调控类黄酮生物合成途径结构基因的表达,本研究开展了一系列分子实验。首先克隆获得AcCHI、AcDFR、AcLAR与AcANR的启动子区域(约 2000 bp),并根据顺式作用元件将其截为不同片段(图S5,A)。经筛选确认,这些启动子片段在适宜的金担子素(AbA)浓度下无自激活活性(图S5,B)。酵母单杂交(Y1H)、双荧光素酶(LUC)与凝胶迁移阻滞实验(EMSA)结果共同证实,AcbHLH137可特异性结合AcCHI启动子上的经典 E-box(CACGTG)基序,并正向激活其转录(图4,A-C和H),表明 AcbHLH137 可促进酚类物质向类黄酮的转化。
图4 AcbHLH137激活AcCHI、AcDFR,抑制AcLAR的表达
此外,AcbHLH137 还可结合AcDFR启动子上的非经典 E-box(CATTTG)基序,并增强其转录活性(图4,A、D、E、I)。由于 DFR 催化生成的产物是花青素与原花青素的共同前体,该结果提示 AcbHLH137 可能同时正向调控花青素与原花青素的生物合成。本研究还通过 Y1H 与 LUC 实验证实,AcbHLH137 可促进AcF3GT1的转录(图 S6)。
Y1H 与 EMSA 实验结果显示,AcbHLH137 可结合AcLAR启动子上的 E-box(CAATTG)基序,但无法结合AcANR的启动子(图4,A和 J)。双荧光素酶实验表明,AcbHLH137可抑制AcLAR的转录(图4,F 和 G),这与转基因愈伤组织中AcLAR表达下调的结果一致。
5、AcbHLH137与AcMYB12存在蛋白互作
为进一步解析AcbHLH137调控类黄酮代谢的分子网络,本研究利用酵母双杂交(Y2H)实验筛选了与其互作的潜在 MYB 转录因子。AcbHLH137 在 Y2H 系统中无自激活活性(图S7)。Y2H 实验证实,AcbHLH137 与 AcMYB12 存在直接的蛋白互作(图5,A)。双分子荧光互补(BiFC)实验进一步在细胞核中验证了二者的互作关系(图5,B)。此外,Pull-down实验结果显示,AcMYB12-HIS 融合蛋白可捕获 AcbHLH137-MBP融合蛋白,而单独的 HIS 标签无法捕获,证实二者在体外也存在物理互作(图5,C)。上述结果表明,AcbHLH137与AcMYB12 可形成蛋白复合体,共同调控类黄酮的生物合成。
图5 AcbHLH137与AcMYB12蛋白间的互作验证
本研究检测了AcMYB12在‘红阳’与 “H16” 中的表达模式,结果显示,随着果实发育,两个品种中AcMYB12的表达量均呈下降趋势(图S8,A)。蛋白序列分析表明,AcMYB12属于R2R3型MYB转录因子,在进化上与SG5亚族MYB成员MdMYB12亲缘关系较近(图S8,B和C);已有研究表明,MdMYB12 参与苹果原花青素的合成调控(Wang 等,2017)。亚细胞定位分析显示,AcMYB12 同时定位于细胞核与内质网(图S8,D)。
6、过表达AcMYB12促进原花青素积累
为探究 AcMYB12 在类黄酮生物合成中的功能,本研究在‘毛花’猕猴桃(A. eriantha)中获得了AcMYB12稳定过表达株系OE-6与OE-8。经 DMACA 染色后,转基因植株叶片的染色更深,表明其原花青素积累水平更高(图6,A);叶片未呈现红色表型,可能与该品种本身不积累花青素的特性有关。过表达株系OE-6与OE-8中AcMYB12的表达量是野生型的 40 倍以上(图6,B)。定量结果显示,过表达株系的原花青素含量显著升高(图6,C),且类黄酮途径结构基因的表达均显著上调(图6,D)。
图6 过表达AcMYB12增强原花青素的积累
同样,在另一组转基因愈伤组织株系(OE-1、OE-2)中(经 qRT-PCR 与凝胶电泳验证,图S9,A和B),DMACA 染色呈现红棕色(图 S9,C和D),原花青素含量同样升高(图S9,E),类黄酮途径基因也呈现不同程度的上调(图S9,F)。上述结果表明,AcMYB12 可促进猕猴桃中原花青素的生物合成。
7、AcMYB12 激活AcDFR、AcF3GT1与AcLAR的表达
进一步研究AcMYB12是否直接调控类黄酮途径结构基因的转录。Y1H 实验结果显示,AcMYB12 可直接结合AcDFR与AcF3GT1的启动子(图7,A);EMSA 实验证实,该互作是通过结合二者启动子上的 MER 基序(AACGGCC/AACCCCC)实现的(图7,H 和 J)。LUC 实验进一步表明,AcMYB12 可正向调控AcDFR与AcF3GT1的表达(图7,B、C、F、G),该调控模式与AcbHLH137 类似,提示 AcMYB12 也可能促进花青素的生物合成。但由于‘毛花’本身不合成花青素,因此在过表达株系中未检测到花青素含量的变化。
图7 AcMYB12激活合AcDFR、AcLAR和AcF3GT1的表达
鉴于AcMYB12过表达株系中原花青素含量升高,且AcLAR与AcANR表达上调,本研究进一步探究了 AcMYB12 与这些基因的调控关系。Y1H 与 LUC 实验结果显示,AcMYB12 可结合AcLAR启动子并激活其表达(图7,D 和 E),但无法结合AcANR的启动子。EMSA 结果证实,AcMYB12 可特异性结合AcLAR启动子上的 MYB CORE 基序(TAGTTA)(图7,I)。上述结果表明,AcMYB12 通过直接激活AcLAR的表达,正向调控原花青素的积累(图7,F 和 G)。
8、AcbHLH137与AcMYB12共同调控花青素与原花青素生物合成的平衡
由于AcbHLH137与AcMYB12均能调控类黄酮途径中AcDFR、AcF3GT1与AcLAR的表达,为明确二者共同作用的调控效应,本研究将 AcbHLH137、AcMYB12分别与携带GUS报告基因的AcDFR、AcF3GT1、AcLAR启动子载体共转化烟草叶片(图8,A)。
图8 AcbHLH137和AcMYB12协同调控AcDFR、AcF3GT1、AcLAR的表达
GUS 活性与染色实验结果显示,AcbHLH137与AcMYB12可协同激活AcDFR的表达,共表达时的GUS染色强度显著强于二者单独表达(图8,B 和 C);二者对AcF3GT1的调控也存在类似的协同效应(图8,F 和 G)。
对于AcLAR,单独表达AcbHLH137时 GUS 活性显著降低,表明其对AcLAR具有负调控作用,这与前文的实验结果一致(图8,D);而单独表达AcMYB12时则可显著激活AcLAR的表达(图8,E)。当两个转录因子共浸润时,GUS活性处于二者单独作用的中间水平,表明 AcbHLH137 与 AcMYB12对AcLAR的转录调控存在拮抗效应。