《Frontiers in Plant Science》：Widely targeted metabolomic and transcriptomic analyses of the effects of blue polarized light and ordinary light on Dendrobium officinale
编辑推荐：
本研究通过转录组学与广靶代谢组学联合分析，系统揭示了蓝光偏振光（BP）与普通光源（白光W、蓝光B）在调控铁皮石斛（Dendrobium officinale）表型、转录表达及代谢谱方面的显著差异。研究发现，偏振光能促进石斛茎部转红，并显著上调类黄酮（flavonoid）和苯丙烷（phenylpropanoid）合成通路的关键基因（如F5H， CHS， F3H）与代谢物（如melanoside A）表达，提示其可作为一种有效的光调控手段，用于提升铁皮石斛中具有药用价值的次生代谢物产量。
1 引言
铁皮石斛（Dendrobium officinale）是中国传统名贵中药材，具有生津润肺、增强免疫力等多种功效。其药用成分主要包括多糖、类黄酮和联苄类化合物。光照是调控植物功能代谢物合成的重要因素。普通光源是指电磁波方向无序的光源，而偏振光源则具有特定方向的规律性，其发光机制与普通光源完全不同，可能为探索光对植物代谢影响的新机制提供了窗口。目前关于偏振光对植物影响的研究有限，特别是在铁皮石斛等药用植物中的分子机制尚不清楚。
2 结果
2.1 不同光照处理下铁皮石斛的生长状态
偏振光能显著影响铁皮石斛的表型。研究显示，在蓝光偏振光（BP）照射下，植株茎秆颜色偏红。白光（W）处理下的叶片面积最大，植株最矮且茎最粗。这些表型差异表明偏振光源与普通光源对铁皮石斛的形态建成有明显不同影响。
2.2 转录组数据测序与组装
研究构建了9个mRNA文库进行测序，测序数据质量高（Q30>97.32%），与铁皮石斛参考基因组的比对率（91.24–92.61%）和单一位点比对率（87.40–89.80%）均很高，证实了数据的可靠性。
2.3 不同光照处理下的差异表达基因分析
转录组分析共鉴定出4729个差异表达基因（Differentially Expressed Genes， DEGs）。其中，BP与W比较（BP-W）中有2448个DEGs，BP与B比较（BP-B）中有2065个DEGs，B与W比较（B-W）中有2763个DEGs。韦恩图显示有187个基因在三个比较组中都发生变化，而各组也有其特有的差异基因。多数DEGs表现为上调表达。
2.4 差异表达基因的GO富集分析
基因本体论（Gene Ontology， GO）富集分析表明，BP-W和BP-B比较组中的DEGs在细胞组分方面显著富集于微管（microtubules），在分子功能方面显著富集于UDP-糖基转移酶活性（UDP-glycosyltransferase activity）和微管结合（microtubule binding）。这提示这些通路在铁皮石斛响应偏振光与普通光差异中扮演重要角色。
2.5 铁皮石斛中光诱导保守基因的功能分类
KEGG（Kyoto Encyclopedia of Genes and Genomes）通路富集分析显示，不同比较组的DEGs富集的通路存在差异。BP-W组的DEGs主要富集于苯并噁嗪类生物合成、类黄酮生物合成等通路；BP-B组主要富集于玉米素（zeatin）生物合成、异黄酮生物合成等通路；B-W组则主要富集于亚油酸代谢、玉米素生物合成等通路。值得注意的是，植物激素信号转导、玉米素生物合成、苯丙烷生物合成等通路在三个比较组中均显著富集，表明它们是光响应的核心通路。
2.6 不同光照处理下的代谢物变化
代谢组学分析（OPLS-DA）表明，BP、B、W三组处理的代谢物谱存在显著分离。在BP-W比较中鉴定出1190个差异丰度代谢物（Differentially Abundant Metabolites， DAMs），BP-B比较中鉴定出812个DAMs，且上调代谢物数量普遍多于下调代谢物。这些特定代谢物（Specific Changed Metabolites， SCMs）可分为多酚、生物碱、类黄酮等18类。在BP-W和BP-B比较中，下调最显著的SCMs均包括13，21-二氢宽木酮（13,21-dihydroeurycomanone）、异牡荆素2’’-O-阿拉伯糖苷（isovitexin 2’’-O-arabinoside）和丹叶大黄素3’-O-葡萄糖苷（rhapontigenin 3’-O-glucoside）。
2.7 DEGs与SCMs的联合分析
使用双向正交偏最小二乘法（O2PLS）对转录组和代谢组进行整合分析，结果显示BP-W和BP-B比较组中两组学数据存在显著相关性。KEGG富集气泡图分析表明，BP-W组中DEGs富集前5的通路包括托烷、哌啶和吡啶生物碱生物合成、类黄酮代谢等；其DAMs富集前5的通路包括丙酮酸代谢、苯丙烷代谢、类黄酮代谢等。BP-B组中DEGs富集前5的通路涉及玉米黄质代谢、异黄酮生物合成等；其DAMs富集前5的通路涉及泛醌生物合成、异黄酮生物合成等。
对共相关基因和代谢物数量最多的前10条KEGG通路进行分析发现，BP-W和BP-B组均包含植物激素信号转导、苯丙烷生物合成、碳代谢等通路，且植物激素信号转导通路的相关基因和代谢物数量均位居第一。
2.8 代谢相关DEGs和SCMs的分析
上述分析表明，偏振光与普通光源之间在植物激素信号转导、玉米素生物合成、苯丙烷生物合成和类黄酮生物合成方面存在显著差异。
对植物激素信号转导相关DEGs的分析发现，生长素（auxin）信号通路中的TIR1、IAA、ARF等13个基因在BP处理下表达水平最高，B处理次之，W处理最低。而细胞分裂素（cytokinin）信号中的AHP、赤霉素（gibberellin）信号中的DELLA基因则在W处理下表达最高。脱落酸（abscisic acid）信号中的PYL在BP组表达最高。乙烯（ethylene）和水杨酸（salicylic acid）信号通路中的多个基因也在不同处理间表现出差异表达模式。
对苯丙烷和类黄酮生物合成通路的分析发现，上调的DEGs包括F5H（阿魏酸-5-羟化酶）、CHS（查尔酮合酶）、F3H（柚皮素-3-加双氧酶）、CYP75B1（类黄酮3’-单加氧酶）等；下调的DEGs包括4CL（4-香豆酰辅酶A连接酶）和I2’H（异黄酮/4’-甲氧基异黄酮2’-羟化酶）。上调的SCMs包括芹菜素（apiengin）、槲皮素（quercetin）、芦丁（rutin）等；下调的SCMs包括苯丙氨酸（phenylalanine）、酪氨酸（tyrosine）、牡荆素（vitexin）等。
2.9 不同光照处理下铁皮石斛DEGs的qPCR验证
qPCR验证了9个与苯丙烷和类黄酮生物合成相关的基因表达。结果显示，F5H、CHS、F3H、CYP75B1、DFR（双功能二氢黄酮醇4-还原酶）、PGT1（根皮苷合酶）和ANS（花青素合酶）的表达量在BP处理组最高，B组次之，W组最低。而4CL和I2’H的表达模式则完全相反，在W组最高，BP组最低。这表明偏振光可能通过正调控F5H、CHS等基因表达，负调控4CL、I2’H表达来影响铁皮石斛功能代谢物的积累。
3 讨论
3.1 铁皮石斛对蓝光偏振光与普通光源的植物激素信号转导响应存在显著差异
KEGG分析和共相关通路分析均表明植物激素信号转导和玉米素生物合成是关键差异通路。DEGs分析揭示了生长素、赤霉素、脱落酸和乙烯信号通路基因的显著差异表达。这种差异可能源于光信号对激素通路的特异性调控，以及偏振光物理特性对光受体激活效率的提升。偏振光因其方向一致性，能更高效地激活光敏色素（如phyB）和蓝光受体，加速钙离子流和活性氧（Reactive Oxygen Species， ROS）信号（如H2O2）的积累，进而更高效地调控激素响应因子。此外，偏振光穿透力强，在提升光合效率的同时也可能导致叶绿体活性氧（cpROS）激增，通过逆向信号与激素信号互作，优先启动类黄酮合成通路以应对光氧化损伤。铁皮石斛作为附生兰科植物，其ERF5–WOX4调控模块和ROS介导的激素通路对偏振光更为敏感，从而优化了其光保护机制。
3.2 铁皮石斛对蓝光偏振光与普通光源的苯丙烷代谢和类黄酮代谢响应存在显著差异
GO和KEGG分析均显示苯丙烷和类黄酮代谢通路显著富集。联合分析也鉴定出大量参与这两条通路的DEGs和代谢物。造成差异的主要原因在于：第一，偏振光对光信号传递的增强效应。其方向一致的电磁振荡特性能更高效地激活光受体，加速ROS等信号分子积累，通过MAPK（Mitogen-Activated Protein Kinase）级联反应优先启动类黄酮合成。第二，对代谢通路的差异调控。偏振光富含蓝光成分，能显著增强类黄酮合成关键上游酶CHS的活性，促进槲皮素等黄酮醇的合成。类黄酮作为光保护剂，在偏振光下加速积累以抵抗定向光辐射引起的氧化应激。而苯丙烷途径的起始酶PAL（苯丙氨酸解氨酶）受光调控但对光质不敏感，其产物主要用于结构支持，在不同光源下积累差异不大。第三，偏振光的生物学优势。它能更有效地穿透叶片，提高光合效率，为次生代谢提供更多碳骨架。类黄酮在偏振光下定向积累于细胞壁和液泡，形成局部高浓度保护层。第四，环境适应性的进化因素。铁皮石斛长期适应林冠下的偏振光环境，其类黄酮代谢通路对偏振光具有高敏感性，以优化光保护策略。
3.3 铁皮石斛响应蓝光偏振光与普通光在微管活性、UDP-糖基转移酶活性和微管结合方面存在显著差异
GO分析显示，BP-W和BP-B比较组中，与微管、UDP-糖基转移酶活性和微管结合相关的术语富集显著性极强。关于微管通路，蓝光偏振光作为一种特殊光学信号，可能因其偏振特性被石斛细胞中的特定光受体识别，触发不同于普通光的信号转导通路。植物细胞内的微管系统对光信号敏感，蓝光能重排皮层微管，而蓝光偏振光可能增强了这种重排的特异性、方向和动态过程，为细胞提供了额外的空间取向信息。UDP-糖基转移酶活性差异可能与蓝光偏振光通过独特的光信号转导，调节该酶基因转录，进而影响底物糖基化程度和修饰位点选择有关，最终影响植物对光环境的适应。微管结合通路差异则可能源于蓝光偏振光诱导的细胞内信号变化，对微管结合蛋白的功能和活性状态产生了不同的调控作用，例如通过磷酸化修饰改变其与微管的亲和力，进而影响细胞器运输等生理过程。
4 材料与方法
4.1 植物材料与光照处理
使用生长一致的铁皮石斛组培苗，在光培养箱中分别进行白光（W）、蓝光（B， 450 nm）和蓝光偏振光（BP， 450 nm， 90°）处理60天。光强为200 μmol·m-2·s-1，光周期12 h/d，湿度55%-60%，温度25±2°C。培养基为1/2 MS固体培养基。以W为对照，B和BP为实验组。
4.2 RNA提取与测序
对三种处理的植株进行高通量测序，每组三个生物学重复。使用Illumina HiSeq平台进行建库和测序。原始数据经质控后与铁皮石斛参考基因组（v1.0）进行比对。
4.3 差异表达基因分析
使用DESeq2软件进行差异表达基因分析。筛选标准为P值<0.01且log2(差异倍数)>1。
4.4 差异表达基因功能分析
使用ClusterProfiler软件，以GO和KEGG数据库为参考，对所有显著差异基因进行功能富集和通路富集分析。
4.5 miRNA和mRNA相对表达水平检测
使用Roche LightCycler 480平台，以肌动蛋白（actin）为内参基因，通过实时荧光定量PCR（qRT–PCR）检测mRNA相对表达水平，采用2-ΔΔCt法计算。
4.6 广靶代谢组检测
4.6.1 代谢物提取： 使用UPLC–ESI–MS/MS系统（Waters Acquity I-Class PLUS UPLC， Applied Biosystems 6500+ Q TRAP MS）分析样本提取物。
4.6.2 LC–MS/MS分析： 详细设置了液相色谱条件和质谱离子源参数，采用多反应监测（MRM）模式。
4.6.3 数据分析： 原始峰面积经总峰面积归一化后，进行主成分分析和Spearman相关性分析评估重复性。使用ropls软件包进行OPLS-DA建模和置换检验，结合差异倍数、P值和VIP值筛选差异丰度代谢物，标准为FC>1， P值<0.05且VIP>1。最后利用超几何分布检验计算代谢物在KEGG通路中的差异富集。
4.7 数据分析
基因表达和功能代谢物数据分析均包含至少3个生物学重复。使用SPSS V19.0进行单因素方差分析（ANOVA），使用GraphPad Prism 6.0和OmicShare在线软件作图。
