多甲氧基黄酮在柑橘属C. depressa果实发育期果皮中的动态代谢谱及生物合成解析

《Metabolites》：Dynamic Profiling of Polymethoxyflavones and Biosynthetic Insights in the Peel of Citrus depressa During Fruit Development

　　

背景/目的：C. depressa Hayata果皮富含多甲氧基黄酮（PMFs），尤以川陈皮素含量最为丰富，其具有多种生物活性。然而，PMF积累的发育动态及其分子调控机制尚不清楚。本研究旨在表征果皮发育过程中的PMF谱，并鉴定调控川陈皮素富集的关键基因与关键时期。方法：于5个发育阶段（S1–S5）采集C. depressa果皮样品。采用高效液相色谱（HPLC）定量川陈皮素含量。采用广泛靶向代谢组学分析PMF谱，并结合转录组测序与加权基因共表达网络分析（WGCNA）鉴定候选结构基因。采用实时荧光定量PCR（qRT–PCR）验证表达模式。结果：川陈皮素呈双峰积累模式，在S1和S4阶段出现峰值（S4阶段为4.637 mg/g鲜重）。代谢组学鉴定出68种PMFs，其中川陈皮素为优势成分（占总PMF峰面积信号的44.5–51.5%）。转录组与WGCNA分析揭示，OMT、CYP450和UGT基因在S4阶段协同上调，qRT–PCR证实其表达峰值与川陈皮素积累同步。结论：川陈皮素在C. depressa果皮中呈双峰积累模式，于S4阶段达到峰值（4.637 mg/g FW），该阶段被确定为最佳采收窗口。O-甲基转移酶（OMT）、细胞色素P450（CYP450）和UDP-糖基转移酶（UGT）基因在S4阶段的协同上调为C. depressa的定向利用与分子育种提供了候选靶标。

研究背景方面，C. depressa Hayata（台湾香檬/山橘仔）原产于中国台湾省及琉球群岛，为东亚特色柑橘种质资源，其果皮富含多甲氧基黄酮（PMFs），以川陈皮素为优势活性单体，在不同栽培品系油胞层中占总PMFs的56.74–64.77%，干燥白皮层中约达24 mg/g。PMFs因高度甲氧基化而具有增强的脂溶性与膜通透性，具备抗炎、抗氧化、抗肿瘤、神经保护及心血管保护等活性，其活性与甲氧基数目及位置密切相关；川陈皮素具有5,6,7,8,3′,4′-六甲氧基取代模式，经Nrf2/HO-1、NF-κB及MAPK通路发挥抗氧化、抗炎与抗凋亡作用。目前存在的问题是：C. depressa果皮PMF积累的发育动态尤其是川陈皮素，及其分子调控机制尚不清楚；尽管已报道其果皮川陈皮素含量极高，但积累确立的具体发育阶段与关键基因不明。研究目的即整合代谢组与转录组系统阐明C. depressa果皮PMFs（尤其川陈皮素）的发育积累动态并鉴定关键调控基因，为分子育种或采后调控提升功能成分提供理论基础，并为确定最佳采收期提供科学依据。

研究人员以福建省亚热带植物研究所标准化栽培基地（原引自台湾屏东，2015年7月引入）的C. depressa为材料，于2025年按花后60、90、120、150、180天（DAFB）采集S1（绿果膨大期）、S2（稳定生长期）、S3（初黄期）、S4（半黄期）、S5（全黄期）五个发育阶段果皮，每阶段3株6年生健康树、每株4个方位各2果混合为1个复合样，设3个生物学重复。主要关键技术方法包括：高效液相色谱（HPLC）定量川陈皮素（外标法，335 nm检测）；广泛靶向代谢组学（SCIEX QTRAP，多反应监测MRM模式，Metware自建数据库MWDB鉴定）；转录组测序（Illumina NovaSeq 6000，HISAT2比对、StringTie组装定量、DESeq2差异分析）结合加权基因共表达网络分析（WGCNA）筛选与川陈皮素相关的模块与候选基因；实时荧光定量PCR（qRT–PCR，以UBQ4为内参，2?ΔΔCT法）验证候选基因表达。

结果部分：（1）川陈皮素在C. depressa果皮中的积累动态：HPLC定量显示川陈皮素含量随发育动态变化，S4最高（4.637 mg/g鲜重），为S1的1.06倍且差异不显著；S2与S5显著低于S4，分别为S1的0.77倍与0.78倍；S3为3.685 mg/g，仅为S1的0.84倍；单因素方差分析显示五阶段间差异显著（p = 0.012），表明积累具发育阶段特异性，S4为最佳采收窗口，且S4近成熟期果实成熟指标明显高于S1（可溶性固形物8.5% vs. 5.0%；总糖61.0 vs. 23.2 mg/g FW）。（2）C. depressa果皮PMFs的积累模式：对68种PMFs经Z-score标准化后以欧氏距离和Ward最小方差法层次聚类，分为5种积累模式——早期优势型（11种，16.2%）、中期峰值型（10种，14.7%）、双峰型（13种，19.1%）、中期特异型（5种，7.4%）与后期爆发型（29种，42.6%）；川陈皮素属双峰型，于S1和S4出现两个积累峰、S3为低谷；pmp000107与mwsTZ409均在S4达峰。（3）不同积累模式下PMFs的结构修饰特征：68种PMFs按结构修饰分为全甲氧基黄酮（9种）、羟基化PMFs（44种）和PMF糖苷（15种）；模式3甲氧基数目最多（4.46 ± 0.97）、羟基数目最少（1.08 ± 1.12）、糖基化比例最低（7.7%）；模式5甲氧基较多（4.31 ± 1.07）但糖基化比例达34.5%；模式4羟基最多（3.20 ± 0.45）且未检出糖基化，表明积累模式与甲氧基化、羟基化及糖基化程度显著关联。（4）关键PMF组分相对峰面积比例动态：川陈皮素在发育过程中始终占比最大（占总PMF峰面积信号的44.5–51.5%）且比例稳定；3′,4′,5′,5,7-五甲氧基黄酮占17.3–21.9%，于S5达峰（21.9%）；甜橙素（9.5–11.2%）于S1达峰后渐降；PDP030604于S3升至4.9%，其余阶段低于3.2%；Wagp008065自S1的3.6%渐降；pmp000107与MWStz409于S4分别升至4.3%与6.0%，揭示以川陈皮素为核心、多组分协同变化的动态积累格局。（5）修饰酶基因表达动态及其与川陈皮素积累的同步性：转录组共鉴定282个修饰酶基因（166个CYP450、75个UGT、41个OMT）；OMT与CYP450呈S1和S4双峰高表达，S4期OMT表达约为S1的1.3倍，UGT主要在S4持续高表达（平均约为S1的3.8倍），与川陈皮素双峰积累及后期糖基化比例升高高度同步；15个文库组装获86,225条unigene，以S1为参照分别鉴定出S2、S3、S4、S5的差异表达基因4679、10,576、17,983和20,067个；WGCNA显示蓝色模块在S4特异高表达且与川陈皮素含量正相关最强（r = 0.76，p = 0.001），从中筛得4个OMT、2个UGT和4个CYP450候选基因，其表达峰均集中于S4，与川陈皮素S4积累峰高度一致。（6）qRT–PCR验证候选修饰酶基因表达：以UBQ4为内参，6个代表性候选基因表达趋势与RNA-seq基本一致，均在S4显著高表达或维持高水平；UGT基因cluster 5991.0上调最强，S4相对表达分别为S1、S2、S3、S5的31.15、11.08、2.72和6.59倍；CYP450基因cluster 19,819.0在S4相对表达最高（为S1的4.98倍）并自S1至S4渐升；OMT基因cluster 16,382.0于S4达峰、S3次之；cluster 34,645.3在S4为S1的2.79倍而S5显著降至S1的0.97倍；cluster 30,032.2与21,636.4在S3–S4高表达而S5显著下调，表明参与PMF结构修饰的OMT、CYP450和UGT基因在S4协同上调或维持高表达活性，与川陈皮素积累动态高度吻合。

讨论部分总结：研究人员报道C. depressa果皮川陈皮素丰富且于S4达峰；与广陈皮等随成熟度增加而下降、且多以干重（DW）计并依赖炮制富集或易降解的情况不同，本研究以鲜重（FW）计（鲜柑橘皮含水量通常>50%，换算干重后更高），S4含量约为四株冲绳香檬油胞层（1.29–1.70 mg/g FW）的3倍，且无需采后陈化或加工富集，据知为首次报道引种C. depressa川陈皮素的发育动态；S2–S3下降可能归因于果实快速膨大导致的稀释效应，S3–S4回升可能反映生物合成或转运活性增强（机制尚不清楚），该动态与脐橙等随成熟逐渐下降的报道不同，提示C. depressa可能具有独特的类黄酮生物合成调控机制。代谢组鉴定68种PMFs（9种全甲氧基黄酮、44种羟基化PMFs、15种糖基化PMFs），种类数超过既往报道，组成以羟基化和糖基化PMFs为主，区别于九里香（以全甲氧基苷元为主）及橘C. reticulata cv. Chachiensis与甜橙（PMF谱相对有限），显示明显的种特异性。川陈皮素相对峰面积比例44.5–51.5%与商业柑橘（50.64%）及甜橙（29.63%）趋势一致而区别于柚（5.12%），提示其优势地位在柑橘属中保守；C. depressa PMFs含2–7个甲氧基，以4或5个甲氧基的中度甲基化成分（34种，50.0%）为主，高甲氧基（≥6个）仅5种但包含核心成分川陈皮素，模式3甲氧基化较高而羟基化与糖基化较低的结构组合可能利于果皮PMF快速积累；糖基化产物仅占22.1%，S4期UGT与OMT共存高表达而产物仍以全甲氧基川陈皮素为主，提示O-甲基化可能在羟基底物竞争中占优。在分子基础方面，CYP450先催化黄酮骨架区域选择性羟基化提供甲基化位点，OMT再将羟基转化为甲氧基并经多个OMT分步协同构建多甲氧基结构，UGT subsequently催化糖基化增强稳定性与溶解度；同种中已分离5个CdFOMT基因且CdFOMT5对3-、5-、6-、7-羟基具广泛甲基化活性，CreOMT3/4/5过表达使总PMF提高11.6–16.9%、RNAi沉默降低12.9–18.6%，证实OMT为关键驱动基因；本研究OMT cluster 16,382.0与21,636.4、CYP450 cluster 19,819.0与34,645.3及UGT cluster 5991.0在S4协同高表达并与川陈皮素峰同步，提示S4为三类修饰酶协同激活的代谢窗口、川陈皮素生物合成在此阶段重新激活；S1峰可能由其他家族成员驱动，S4未显著超过S1可能与UGT cluster 5991.0强烈上调导致糖基化分流抵消生物合成有关，此推断尚待功能验证与糖苷产物定量分析。

研究结论部分翻译：研究人员系统表征了C. depressa果皮在果实发育过程中川陈皮素与PMFs的积累模式及分子表达谱。C. depressa果皮川陈皮素含量在S4阶段达峰（4.637 mg/g FW），无需陈化或加工即可获得较高含量，表明C. depressa可能是川陈皮素的天然原料来源。所鉴定68种PMFs涵盖多种结构类型（苷元、羟基化形式与糖基化形式），川陈皮素为优势成分，占总PMF峰面积信号的44.5–51.5%，PMF谱总体以中度甲基化形式为特征且糖苷比例较低（22.1%）。三类修饰酶基因（OMTs、CYP450s与UGTs）在S4协同高表达，其表达峰与川陈皮素积累峰吻合，提示S4是C. depressa果皮PMF生物合成与积累的关键调控窗口。尽管这些候选基因的功能及其对川陈皮素积累的具体贡献尚需验证，上述发现为通过优化采收期和调控关键基因提高C. depressa功能成分含量提供了理论基础。