《Horticulturae》：Integrated Metabolome and Transcriptome Analyses Reveal the Developmental Program of Anthocyanin Pigmentation in Lavender (Lavandula angustifolia Mill.) Flowers
编辑推荐：
研究人员通过整合非靶向代谢组学（untargeted metabolomics）、RNA-seq转录组学与系统发育分析，结合花发育过程的形态学与显微观察，解析了狭叶薰衣草（Lavandula angustifolia Mill.）‘Xinxun 2’花部花青素
研究人员通过整合非靶向代谢组学（untargeted metabolomics）、RNA-seq转录组学与系统发育分析，结合花发育过程的形态学与显微观察，解析了狭叶薰衣草（Lavandula angustifolia Mill.）‘Xinxun 2’花部花青素（anthocyanin）呈色的发育程序。花青素含量呈现“起始—峰值—平台”轨迹，在盛花期（full flowering）达到高峰；然而在末花期（late flowering）色素总量稳定时，花冠仍持续向紫蓝方向偏移，表明最终色相不仅由色素数量决定。分子层面，飞燕草素分支关键酶F3′5′H（黄酮类3′,5′-羟化酶）是紫罗兰表型的主效决定因子；结构基因在初花期（early flowering）集中表达，而晚期修饰基因AOMT（花青素O-甲基转移酶）在盛花期被诱导，竞争性的黄烷-3-醇分支酶ANR同步上调；两个R2R3-MYB第六亚组转录因子MYB90与MYB114被识别为候选激活因子。该研究首次建立薰衣草花部色素沉积的阶段分辨分子框架，并确定盛花期为色素提取最优窗口，为花色分子育种提供具体靶标。
论文解读：薰衣草花青素呈色发育程序的多组学解析
研究背景与立项逻辑
狭叶薰衣草（Lavandula angustifolia Mill.）的紫罗兰至紫色花序是其观赏与经济价值的核心，但长期以来其花色形成的分子基础并不清晰。唇形科（Lamiaceae）植物虽普遍积累飞燕草素（delphinidin）衍生、并被芳香族与脂肪族酰基修饰的花青素，但薰衣草特异的结构基因组合、修饰逻辑与转录调控网络仍属空白。已有染色体水平基因组与多组学研究多聚焦于萜类与腺毛发育，花色仅停留于基础色谱调查与单个LaCHS克隆层面，缺乏把“代谢物积累—结构基因表达—上游转录因子”放在同一发育时间轴上读取的阶段分辨框架。花青素最终色相由B环羟基化分流、糖基化/甲基化/酰基化修饰、液泡转运与液泡微环境（pH、共色素、金属络合）共同决定；其中MYB–bHLH–WD40（MBW）复合体是保守主控模块。研究人员据此假设：薰衣草紫调源于飞燕草素类花青素发育期积累，由阶段特异结构基因与修饰基因协同，并由R2R3-MYB激活因子驱动。
研究总体结论与意义
研究人员以‘Xinxun 2’为材料，跨芽期、初花期、盛花期、末花期采样，联合非靶向代谢组（UPLC-MS/MS）与RNA-seq，重构了“起始—峰值—平台”的呈色发育程序：色素总量盛花期达峰，末花期不降解但色相向紫蓝微调；飞燕草素分支（F3′5′H）主导紫调，早期骨架基因与晚期修饰基因（AOMT及BAHD酰基转移酶、糖基转移酶）时间分区表达；ANR代表竞争性的黄烷-3-醇分流在盛花期被激活；MYB90与MYB114作为第六亚组R2R3-MYB枢纽候选调控因子。该工作首次给出薰衣草花色阶段分辨分子框架，并把盛花期定义为色素提取与育种干预的最优窗口，且因内源已具功能性F3′5′H，蓝化育种无需整套异源飞燕草素通路转入，只需靶向修饰与调控基因编辑。
主要关键技术方法
研究以新疆伊犁‘Xinxun 2’五年生田间植株为样本队列，设芽期、初花期、盛花期、末花期，代谢组每阶段四生物学重复、转录组三生物学重复。技术主线为：非靶向代谢组学（LC-MS/MS）定性定量花青素与黄酮衍生物；RNA-seq转录组组装与DESeq2差异表达分析；基于KIPEs3多证据框架（序列相似度+功能残基保留≥75%+最大似然系统发育）鉴定结构基因与修饰基因；MYB_annotator结合系统树界定R2R3-MYB亚组；Pearson相关与Cytoscape构建基因—代谢物共表达网络；RT-qPCR用La_Actin内参验证趋势。
研究结果
3.1 花发育形态与花青素动态
通过立体显微与细胞显微观察，结合pH差示法测花青素含量，研究人员发现：芽期无色，初花期花冠顶端淡紫，盛花期深紫、含量约0.32 mg/g FW达峰，末花期维持平台并略偏蓝。色素定位于表皮细胞液泡；末花期色移不依赖新合成，而可能与液泡微环境（pH/共色素）有关。
3.2 花青素积累的代谢组特征
PCA显示三花期明显分离；EF→FF差异代谢物显著上调，黄酮/花青素类在FF与LF高积。飞燕草素、矢车菊素、天竺葵素衍生物随发育上升，证明紫蓝表型由多类修饰花青素协同而非单一色素决定。
3.3 转录组与结构基因鉴定
从基因组预测蛋白中鉴定出苯丙烷—黄酮核心结构基因（PAL、C4H、4CL、CHS、CHI、F3H、F3′H、F3′5′H、DFR、ANS）及分支酶FLS、FNS II、ANR与arGST。核心骨架基因在EF达峰；AOMT、3AT、5MAT、F3GT等修饰基因与产物积累同向。R2R3-MYB系统树锁定第六亚组MYB90（La16G01215）与MYB114（La02G01755）为高表达候选激活子。
3.4 共表达网络与关键调控子
对20个基因与24个核心代谢物做Pearson网络：上游骨架基因关联苯丙烷前体，晚期修饰基因特异关联花青素衍生物。MYB90与MYB114处网络中枢，与多种花青素强正相关，符合MBW模块中激活型MYB特征。
3.5 花青素途径重构
整合序列、功能残基与系统发育，重构花冠途径：F3′5′H分流至飞燕草素，AOMT依次甲基化生成矮牵牛素/锦葵素型；3AT、5MAT、F3GT负责酰化、丙二酰化与糖基化；arGST参与液泡隔离；ANR竞争消耗底物流向黄烷-3-醇。表达架构呈“早合成—晚修饰—并行分流”。
讨论与结论翻译
讨论指出，薰衣草紫调不是色素量的简单反映，而是发育编程、空间分隔与层级调控的涌现属性：液泡隔离+表皮乳突光学界面构成呈色基础；末花期蓝移可由液泡碱化、分子内/间共色素、金属络合解释，但需直接测pH与共色素验证。F3′5′H强诱导使薰衣草天然具备蓝化底盘，AOMT甲基化与酰化共同把吸收峰向蓝移。MYB90/MYB114为候选主控，但直接转录激活仍待双荧光素酶、酵母单杂交与基因编辑确认。
结论部分译文：
通过跨三个花期的非靶向代谢组、RNA-seq、系统发育与表型整合分析，本研究确立：薰衣草紫—紫色调是发育编程且空间分区的进程，而非色素丰度的简单映射。花青素含量呈“起始—峰值—平台”轨迹，盛花期达峰；末花期色素量稳定却出现细微紫蓝偏移，该变化可能反映液泡微环境（pH和/或共色素）的内部理化调整，而非显著视觉信号。飞燕草素分支（F3′5′H）是紫调主效决定因子；途径时间分区表现为结构基因初花期达峰、AOMT盛花期诱导、竞争性黄烷-3-醇酶ANR同期上调；两个第六亚组R2R3-MYB（MYB90、MYB114）为该程序的候选激活子。盛花期被确定为色素提取最优窗口；鉴于薰衣草已具功能性F3′5′H分支，通过靶向改造修饰与调控基因获得真蓝花冠，比转入整套异源飞燕草素通路更可行。所述调控关系仍为相关性，需功能验证。该研究为薰衣草花色提供了阶段分辨分子框架，并为从机制走向应用的花色育种给出具体靶点。
