《Genes》：Time-Series Transcriptome Reveals Differential Regulation of Flavonoid Biosynthesis in Dendrobium huoshanense Under Musical–Acoustic Stimuli: Insights for Green Horticultural Production
编辑推荐：
背景：声音振动是一种新兴的物理信号，能够影响植物生长与次生代谢；然而，珍贵药用兰科植物对复杂声学刺激的全基因组转录响应仍未得到充分表征。方法：本研究对暴露于两种具有对比声学特征的音乐结构化声学刺激（《森林幻想曲》与《二泉映月》）及无音乐对照条件下的驯化霍山石斛
背景：声音振动是一种新兴的物理信号，能够影响植物生长与次生代谢；然而，珍贵药用兰科植物对复杂声学刺激的全基因组转录响应仍未得到充分表征。方法：本研究对暴露于两种具有对比声学特征的音乐结构化声学刺激（《森林幻想曲》与《二泉映月》）及无音乐对照条件下的驯化霍山石斛（Dendrobium huoshanense）组培苗进行了时序RNA测序（RNA-seq）；分别在1、6、12和24小时采集叶片样本，并采用实时荧光定量聚合酶链式反应（qRT-PCR）验证转录组数据。结果：两种声学处理触发了类黄酮生物合成的差异性转录重编程。在《森林幻想曲》处理下，早期生物合成基因查尔酮合酶（CHS）和柚皮素7-O-甲基转移酶（NOMT）被迅速上调；而《二泉映月》处理则持续抑制后期修饰基因（柚皮素3-双加氧酶（F3H）、类黄酮3′-单加氧酶（F3′H）、羟基肉桂酰辅酶A莽草酸/奎宁酸羟基肉桂酰转移酶（HCT）、咖啡酰辅酶A O-甲基转移酶（CCoAOMT））。基于Pearson相关性的网络分析鉴定出差异表达基因（DEG）衍生的枢纽结构基因（肉桂酸4-羟化酶（C4H）、F3′H、CCoAOMT）及共表达的转录因子候选基因（MYB、WRKY、NAC）。结论：上述结果表明，霍山石斛能够在转录水平上区分复杂的音乐信号，提示短期自然风格声学刺激有望成为提升这一珍贵药用兰科植物品质的园艺策略。
一、研究背景与意义
霍山石斛（Dendrobium huoshanense）是我国珍稀且特有的药用兰科植物，其茎部富含多糖、类黄酮、生物碱、酚类、菲类和联苄类化合物等多种活性成分，具有抗氧化、免疫调节和保肝等药理作用。其中，类黄酮类化合物是评价霍山石斛药用价值的核心品质标志物。然而，由于过度采集和生境破坏，野生霍山石斛资源已严重枯竭，人工栽培群体的次生代谢物含量差异较大，亟需开发低成本、无污染的物理品质提升技术。
已有研究表明，声波作为植物可感知的机械振动信号，能够改变基因表达和次生代谢物积累。特定声频可促进绿豆幼苗生长、提高番茄总酚和维生素C含量；月见草花在接触模拟传粉者的声音振动后数分钟内即增加花蜜糖浓度；不同音乐风格可诱导芝麻菜产生不同的生理响应，证明植物能够区分不同的声学信号。在分子层面，植物质膜机械敏感离子通道可将声音引起的机械张力转化为胞内Ca2?信号，进而通过钙调素和钙依赖性蛋白激酶（CDPK）介导的信号级联调控转录因子活性，最终改变下游代谢基因的转录。
然而，现有植物声学研究的刺激多为简单的单频纯音，关于药用植物如何响应具有不同声学特征和风格属性的复杂音乐信号的研究仍然匮乏。类黄酮生物合成相关酶如查尔酮合酶（CHS）、羟基肉桂酰辅酶A莽草酸/奎宁酸羟基肉桂酰转移酶（HCT）和咖啡酰辅酶A O-甲基转移酶（CCoAOMT）在很大程度上决定了霍山石斛的药用品质，但其对音乐声学刺激的转录响应尚未得到充分表征。为此，研究人员对暴露于两种具有高度分化声学特征的代表性音乐处理的霍山石斛组培苗进行了时序RNA-seq分析，旨在表征不同音乐刺激引发的全局转录组变化、解析类黄酮生物合成核心基因的差异时间表达模式，并通过共表达网络分析鉴定枢纽基因和候选转录因子，为探索基于音乐的物理栽培技术、推动这一珍贵药用兰科植物的无化学投入绿色园艺改良提供基础转录组学证据。
二、主要技术方法
研究人员以安徽六安亿湖生物科技有限公司提供的一年生驯化霍山石斛组培苗为材料，在25°C、65%相对湿度、1950勒克斯光照强度、12小时光/12小时暗的光周期条件下培养。设置三个处理组：无音乐对照组（CK）、循环播放《森林幻想曲》的森林音乐组（Forest）和循环播放二胡曲《二泉映月》的二泉音乐组（Erquan），播放强度均为60–70 dB。音乐仅在光照期连续播放12小时，暗期停止。在音乐处理开始后1、6、12、24小时采集叶片样本，每个时间点每个处理组随机选取6株幼苗混合为一个生物学重复，每组设置3个生物学重复，共获得36个有效测序样本。主要技术方法包括：利用fastp进行测序数据质控，STAR将clean reads比对至霍山石斛染色体级别参考基因组（登录号CNA0014590），featureCounts进行基因表达定量，DIAMOND进行NR、Swiss-Prot、KOG、GO和KEGG数据库功能注释，PlantTFDB进行转录因子家族注释；采用DESeq2进行差异表达基因（DEG）筛选；通过clusterProfiler进行GO和KEGG富集分析；基于Pearson相关系数构建共表达网络并依据节点度鉴定枢纽基因；采用qRT-PCR验证转录组数据的可靠性。
三、主要研究结果
3.1 测序数据质量控制和基因组注释统计
三个处理组的clean reads数量均衡，CK为3.47×10?–5.80×10?，Forest为3.51×10?–5.04×10?，Erquan为3.57×10?–5.04×10?。碱基质量高，Q20范围为99.65%–99.74%，最低Q30值达97.74%。GC含量在45.31%–46.27%之间。总基因组比对率为96.37%–96.89%，唯一比对率为87.44%–88.28%。共分析20,247个注释基因，其中20,241个在NR数据库有注释（99.97%），14,564个在PFAM（71.93%），10,011个在Swiss-Prot（49.44%），8853个在GO（43.72%），8363个注释为转录因子（41.30%），5405个获得KEGG注释（26.70%）。
3.2 样本重复性和全局转录组分离
生物学重复间的Pearson相关系数均高于0.94，表明实验重复性良好。PCA分析显示PC1解释20.6%的总方差，PC2解释16.2%。CK、Forest和Erquan样本在二维空间中清晰分离，并按处理类型和采样时间聚类。对前2000个最可变基因的层次聚类结果与PCA一致，音乐处理类型和处理持续时间是转录组差异的主要驱动因素，确认了实验分组的有效性。
3.3 差异表达基因的全局特征
在8个处理对对照的比较中，两种音乐处理均在12小时达到DEG数量峰值。Forest-12小时产生1271个DEG（489个上调，782个下调）；Erquan-12小时产生1168个DEG（294个上调，874个下调）。Erquan处理下下调基因比例更高。KEGG富集呈现处理特异性模式：苯丙烷和类黄酮生物合成途径在Forest-1小时对CK中显著富集；类黄酮生物合成途径在Erquan-24小时对CK中达到本研究最高富集显著性（padj = 1.96×10??）。
3.4 两类类黄酮生物合成基因的表达模式
在ko00941和ko00944注释的38个基因中，23个基因在所有样本中稳定表达。根据时间热图表达行为，类黄酮生物合成途径的核心酶分为A类和B类。A类基因（CHS、NOMT）仅在Forest处理下1小时和6小时显著上调，在Erquan处理下不变或受抑制。B类基因（F3H、F3′H、HCT、CCoAOMT）在Erquan处理下多个时间点持续下调，在Forest处理下保持稳定。值得注意的是，HCT在Erquan组所有四个时间点均下调。
3.5 时间点特异性持续差异表达基因
UpSet图分析显示，在所有四个时间点均存在差异表达的基因数量有限：Forest持续上调22个、持续下调11个；Erquan持续上调8个、持续下调23个。将至少在3个时间点显著差异的基因定义为持续DEG后，60个Forest特异性持续上调基因在嘌呤代谢途径显著富集（FDR = 0.019）；74个Erquan特异性持续下调基因在类黄酮生物合成途径显著富集（FDR = 0.038）。
3.6 qRT-PCR验证
选取A类和B类类黄酮生物合成模块中6个代表性结构基因（CHS、NOMT、F3H、F3′H、HCT、CCoAOMT）进行qRT-PCR验证，以TBP为内参基因。qRT-PCR获得的上下调方向与基于FPKM的RNA-seq趋势总体一致，尽管个别样本存在定量幅度偏差，证实转录组数据集可靠，可用于下游共表达网络构建等分析。
3.7 类黄酮生物合成核心基因的共表达网络
筛选出13个候选类黄酮核心基因用于网络构建。在|r| ≥ 0.8阈值下，获得266个显著共表达对（229个正相关，37个负相关），包含250个网络基因。依据节点度拓扑指标，C4H拥有最多的共表达伙伴（117个），其次是F3′H（58个）和CCoAOMT（31个），被鉴定为枢纽结构基因。网络内嵌28个来自18个家族的转录因子基因，其中WRKY（6个）、MYB（3个）和NAC（3个）是最丰富的家族。对全部250个网络基因的KEGG富集显示，类黄酮生物合成途径显著性最高（15个基因，FDR = 6.40×10?1?），其次为二苯乙烯类、二芳基庚烷类和姜辣素生物合成、苯丙烷生物合成和α-亚麻酸代谢。
四、讨论与结论
讨论部分指出，自然风格音乐《森林幻想曲》处理可在1–6小时内快速激活类黄酮生物合成上游基因，为园艺生产提供了实践启示：在幼苗期施加短期自然风格音乐刺激，可作为无化学投入的物理措施，潜在提升霍山石斛的类黄酮相关药用价值。相反，慢节奏、低频的《二泉映月》音乐可能诱导类应激转录响应并抑制多个后期类黄酮修饰基因。两种音乐处理的作用时间模式存在差异：Forest主要在前1–6小时产生瞬时诱导效应，而Erquan产生持续至24小时的转录抑制效应。共表达网络分析鉴定出C4H、F3′H和CCoAOMT为与MYB、WRKY和NAC转录因子候选基因共表达的核心基因。研究同时指出多项局限性：缺乏代谢物定量和基因功能验证；仅采用两种音乐风格，无法排除曲目特异性效应；仅设置单一12小时光照期音乐暴露时长；植物声学感知的完整信号转导级联仍需进一步阐明。未来研究将整合广泛靶向代谢组学和基因功能验证，系统阐明声学调控类黄酮生物合成的分子机制，并优化声辅助绿色栽培策略。
结论部分翻译：时序转录组分析表明，霍山石斛在暴露于不同风格的音乐声学刺激时，类黄酮生物合成发生分化的转录重编程。自然风格《森林幻想曲》在1–6小时内快速上调早期类黄酮生物合成基因（CHS、NOMT），而忧郁风格《二泉映月》持续抑制多个后期修饰基因（F3H、F3′H、HCT、CCoAOMT）。共表达网络分析确定C4H、F3′H和CCoAOMT为与候选MYB、WRKY和NAC转录因子共表达的枢纽结构基因。研究结果表明，短期自然风格声学刺激是提升霍山石斛品质的一种有前景的无化学投入园艺方法。未来工作需要补充代谢组学测量和基因功能验证，以确认音乐声学处理在实际园艺生产中的实际功效。
