整合代谢组与转录组鉴定酸橙果皮发育期黄酮生物合成候选调控因子及MYB11–CHS2启动子结合

《Food Chemistry: Molecular Sciences》：Integrated Metabolomic and transcriptomic analysis identifies candidate regulators of developmental flavonoid biosynthesis in sour Orange Peel

　　

整合代谢组与转录组分析鉴定酸橙果皮发育期黄酮生物合成的候选调控因子。
•对八个发育阶段的整合代谢组与转录组分析揭示，代谢模式由早期糖基化黄酮（glycosylated flavonoids）向晚期多甲氧基黄酮（polymethoxylated flavonoids, PMFs）发生时间性转变。
•早期下调的UDP-糖基转移酶（UGTs）与持续活跃的O-甲基转移酶（OMTs）之间的底物竞争，驱动PMFs的选择性积累，尽管总黄酮含量整体下降。
•共表达网络分析优先筛选出一个黄酮相关turquoise模块，该模块包含35个核心结构基因和66个候选转录因子（transcription factors, TFs）。
•酵母单杂交（yeast one-hybrid, Y1H）实验证实，转录因子MYB11、MYBH和发育调控因子ANT直接结合CHS2启动子，从而建立了一条连接果实形态发生与特化黄酮代谢的新调控联系。
•核定位的MYB11被验证为核心转录激活因子，其直接结合并上调入口基因CHS2。

研究背景与问题：酸橙（Citrus aurantium L.）为多年生木本植物，具有药用和食用价值，其果实和果皮富含生物活性次生代谢物，尤其是黄酮类化合物（flavonoids），具有抗氧化、抗炎和抗肿瘤等活性。尽管其经济与药理意义重要，果皮黄酮组成随发育变化的转录过程仍不完全清楚。黄酮生物合成起源于苯丙烷途径（phenylpropanoid pathway），核心酶PAL、CHS、F3H等相对保守，但转录调控复杂且具物种特异性；苹果和葡萄等物种的调控架构不能直接外推至酸橙。代谢组与转录组整合已用于解析复杂代谢物生物合成架构，但酸橙尚无类似整合分析，尤其对医学价值较高的PMFs高丰度特征及其转录调控架构尚未探索。因此，研究人员对酸橙果皮八个发育阶段开展整合多组学研究，假设发育变化伴随黄酮代谢与转录调控的协同重编程；通过代谢组、RNA测序和加权基因共表达网络分析（weighted gene co-expression network analysis, WGCNA）构建候选调控网络，并以Y1H测试候选TFs与CHS2启动子的互作，进一步对MYB11进行亚细胞定位和双荧光素酶报告（dual-luciferase reporter）验证。该研究意义在于提供候选调控框架，并区分由酵母和植物报告实验支持的互作与由多组学推断的关联。

关键技术与样本来源（≤250字）：样本来自5年生酸橙树，种植于江西宜春湿地公园常规露天果园，2023年5—8月按果实大小和颜色设S1—S8八个阶段，每半月采样；初采5个独立生物重复，选取总黄酮含量最相似的3个作为最终重复。关键技术包括UPLC-QTOF-MS/MS靶向黄酮代谢组分析、Illumina HiSeq X Ten RNA-seq（150 bp paired-end，比对甜橙参考基因组）、WGCNA、FIMO启动子cis元件预测与Cytoscape网络可视化、定量实时PCR（quantitative real-time PCR, qRT-PCR）验证、Y1H、MYB11-GFP亚细胞定位及双荧光素酶报告实验。

研究结果：
3.1 酸橙果皮发育过程中的形态学特征与黄酮谱：通过监测S1—S8果实形态、果体积、果肉和果皮生物量及果皮/果重比，结合AlCl₃法测定总黄酮含量（total flavonoid content, TFC）和UPLC-QTOF-MS/MS靶向黄酮分析，发现TFC由S1的442.07±5.85 mg/g降至S8的91.20±1.98 mg/g；检出13个黄酮代谢物，11个差异积累黄酮（differentially accumulated flavonoids, DAFs）分为两簇，簇1糖基化黄酮早期高后降，簇2 PMFs逐渐增加并于S6达峰后缓降；鉴定出rhoifolin、naringin、neohesperidin、didymin及sinensetin等PMFs，提示发育过程中由糖基化向多甲氧基化转变，但剖面本身不能建立UGT/OMT差异活性或前体—产物关系。
3.2 不同发育阶段酸橙果皮转录组分析：RNA-seq质量较高，获得23,002个注释基因和1,903个新基因；主成分分析显示阶段分离，qRT-PCR验证的12个黄酮相关基因与转录组数据高度相关（r>0.8）；鉴定9,691个差异表达基因（differentially expressed genes, DEGs；log₂FC≥1，FDR<0.05），早期上调基因占优后转为下调；DEGs分为六簇，簇6于S1达峰后下降、簇3先升后降，与黄酮积累趋势时间相似；簇6富集于苯丙烷生物合成和类黄酮生物合成，PAL、CHS、F3H等富集，提供早期黄酮途径转录活性较高的线索。
3.3 黄酮生物合成结构基因表达动态及其与代谢物积累相关性：检查45个苯丙烷/黄酮途径相关候选基因，其中35个位于turquoise模块并保留于候选网络；LOC102617232和LOC102622674注释为咖啡酰辅酶A O-甲基转移酶样（caffeoyl-CoA O-methyltransferase-like, CCoAOMT-like）基因，早期S1—S2高表达后总体下降，不平行于PMFs晚期积累；该结果仅建立注释性通路归属和发育转录模式，未鉴定催化底物或直接贡献PMF生成。
3.4 酸橙果皮TFs表达模式：鉴定812个差异表达TFs，分属67个家族，约半数家族在S1—S2高表达，MYB、NF-YC和bHLH（RB）等FPKM较高；Mfuzz聚类中簇5含189个TFs，其表达趋势与簇1 DAFs及DEG簇6时间相似，故保留为候选调控因子，但功能需实验验证。
3.5 WGCNA鉴定黄酮相关基因模块：WGCNA将转录组分为12个模块，turquoise模块含2,831个基因，与黄酮动态最相关，与naringin、didymin、neohesperidin及TFC强正相关（r>0.9），与sinensetin等4种PMFs显著负相关；该模块基因多在S1高表达，并与DEG簇6重叠；KEGG富集于苯丙烷生物合成、次生代谢物生物合成和植物激素信号转导。UGT表达异质，CCoAOMT-like转录早期高后降，而PMFs于S6达峰；GA生物合成相关基因早期高后降，与部分苯丙烷基因相似，但仅提示关联，不建立通量变化或因果关系。
3.6 黄酮生物合成转录调控网络构建：将turquoise模块与TF簇5取交集，得116个候选TFs，其中66个MYB、bHLH、zf-HD、ERF等家族成员保留，并与35个黄酮结构基因显著共表达，作为候选枢纽基因；CHS2和UGT79A6与6个代表性TFs表达趋势平行，R²>0.7，支持TF—基因表达相关性。
3.7 TF结合CHS2启动子实验验证：以CHS2为优先对象，Y1H使用242-bp CHS2启动子片段、100 mM 3-AT抑制背景，发现MYB11、MYBH和ANT支持选择性培养基上生长，bHLH89未检测到结合；MYB11-GFP显示主要核定位并与核mCherry共定位；双荧光素酶实验中SK-MYB11与0800-LUC-CHS2共浸润较启动子对照显著提高LUC/REN比值，空报告无激活。结果支持MYB11为候选CHS2调控因子；MYBH和ANT仅为启动子结合候选，未测定位及激活/抑制；上述实验不建立内源CHS2表达或酸橙黄酮积累改变。

讨论与结论：讨论部分总结认为，发育期TFC下降与上游苯丙烷/黄酮基因转录降低一致，但转录水平不能证明生物合成通量；PMFs晚期积累未伴随两个CCoAOMT-like基因表达升高，家族层面先例不能直接转移至这两个基因，PMF丰度还可能反映合成、周转、运输和储存差异，均未直接量化；turquoise模块含35个候选结构基因和66个候选TFs，MYB11获启动子结合、核定位及报告激活支持，MYBH和ANT仅为启动子结合候选，ANT与发育调控的联系尚待定位、报告、遗传和蛋白互作实验；GA相关转录早高后降与部分上游基因相似，但转录不能反映内源GA浓度或信号活性；bHLH89在所用Y1H条件下未检测到结合。研究结论翻译：本研究结合八个酸橙果皮发育阶段的代谢组与转录组分析，记录了糖基化黄酮和多甲氧基化黄酮的对比积累模式。WGCNA鉴定出一个黄酮相关turquoise模块，包含35个候选通路相关结构基因和66个候选转录因子。在酵母单杂交条件下，MYB11、MYBH和ANT结合受试CHS2启动子片段，MYB11还在本氏烟草（N. benthamiana）中显示核定位并提高CHS2启动子报告活性。两个参考注释的CCoAOMT-like基因LOC102617232和LOC102622674的表达模式得到表征，但其酶功能及晚期PMF积累的直接转录解释均未建立。这些结果为进一步功能评估提供候选，而未赋予无支持的催化或调控作用。