《Biology》：Integrated Transcriptomic and Metabolomic Analyses Identify Candidate Transcription Factors Associated with Flavonoid and Coumarin Accumulation in Psoralea corylifolia
编辑推荐：
背景：补骨脂(Psoralea corylifolia)为广泛使用的传统药用植物，其干燥成熟果实为主药用部位。黄酮与香豆素为该物种主要生物活性成分。然而，补骨脂不同组织间活性成分的组织特异性代谢轮廓、基因表达模式及潜在调控因子仍知之甚少。方法：研究人员以补骨脂
背景：补骨脂(Psoralea corylifolia)为广泛使用的传统药用植物，其干燥成熟果实为主药用部位。黄酮与香豆素为该物种主要生物活性成分。然而，补骨脂不同组织间活性成分的组织特异性代谢轮廓、基因表达模式及潜在调控因子仍知之甚少。方法：研究人员以补骨脂根、茎、叶、花、果五类组织为材料，开展广泛靶向代谢组学(widely targeted metabolomics)与转录组学(transcriptomics)整合分析，并结合加权基因共表达网络分析(WGCNA)，筛选与黄酮和香豆素积累相关的共表达模块及候选调控因子。结果：不同组织在代谢组与转录组水平均呈现显著组织特异性，果实与其他组织差异最明显。异补骨脂查耳酮(isobavachalcone)、补骨脂黄酮(bavachin)和补骨脂素样黄酮(corylin)等特征活性成分在果实中特异性高积累。差异表达基因主要富集于苯丙烷生物合成(phenylpropanoid biosynthesis)、黄酮生物合成(flavonoid biosynthesis)和异黄酮生物合成(isoflavonoid biosynthesis)通路。WGCNA显示洋红色模块(magenta module)与果实组织及上述活性代谢物含量显著正相关。从该模块鉴定出6个候选转录因子，并基于转录因子–通路基因共表达网络（344条边，|r| ≥ 0.8，p < 0.05）与连接度指标分为三档：优先枢纽为Cluster_22013.0（C3H型锌指转录因子）与Cluster_21217.8（NF-YA；拟南芥同源NFYA9/AT3G20910）；Cluster_9299.0（Rcd1-like）与Cluster_20910.0（NAC；拟南芥同源NAC002/AT5G04410）为高连接正相关的候选；Cluster_31326.1（bZIP）与Cluster_11436.0（C2H2）为仅与通路基因呈负相关的候选。结论：该研究刻画了补骨脂组织特异性代谢与转录模式，鉴定出与果实黄酮和香豆素积累相关的候选共表达模块及转录因子，所优先排序的转录因子梯队为后续功能研究提供基础。
研究背景方面，补骨脂(Psoralea corylifolia)是豆科一年生直立草本，干燥成熟果实为法定药用部位，常用于皮肤病、骨质疏松等。其药理活性主要来自香豆素与黄酮，前者如补骨脂素(psoralen)和异补骨脂素(isopsoralen)，后者如异补骨脂查耳酮、补骨脂黄酮和补骨脂素样黄酮。香豆素与黄酮均源于苯丙烷途径(phenylpropanoid pathway)，以对香豆酰辅酶A(p-coumaroyl-CoA)为共同前体；经苯丙氨酸解氨酶(PAL)、肉桂酸4-羟化酶(C4H)和4-香豆酰辅酶A连接酶(4CL)后分别进入香豆素支与黄酮支。目前补骨脂中催化呋喃香豆素特征骨架修饰的基因及其组织特异性表达调控因子报道很少，补骨脂属乃至豆科中调控呋喃香豆素组织特异性积累的转录因子尚未被分离，属内比较转录组研究几乎空白，基因–代谢物关联网络缺失限制了活性成分组织差异候选基因的识别。因此，研究人员开展五组织整合多组学研究，以为补骨脂活性成分生物合成与调控网络提供候选资源与数据基础。该论文发表在《Biology》。
主要关键技术方法方面，研究人员于广东云浮云安都杨镇工业转移园60号地块采集同日龄一年生补骨脂成熟植株，设根、茎、叶、花、果五组织，每组织3个生物学重复，每个重复混采6株，共18株，采用阻断配对采样。代谢端用70%甲醇提取后超高效液相色谱–串联质谱(UPLC-MS/MS)广泛靶向代谢组学检测，依数据库匹配注释代谢物并以相对峰面积定量。转录端用CTAB结合pBIOZOL提总RNA，建链特异性mRNA文库，Illumina NovaSeq X Plus进行PE150测序；因无染色体级参考基因组，用Trinity从头组装、Corset聚类得到Unigene，RSEM算FPKM，DESeq2做差异表达，iTAK鉴定转录因子。再联合OPLS-DA筛差异积累代谢物，WGCNA筛与果实和活性代谢物相关的模块，并在目标模块内用Pearson相关构建转录因子–通路基因共表达网络以优先排序候选转录因子。
3.1. Metabolite Profiling of P. corylifolia：研究人员经广泛靶向代谢组学在五组织中鉴定2830个化合物，黄酮类最多，其次萜类、木脂素类、香豆素类；PCA与样本相关分析显示五组织代谢轮廓分明，果实与其他组织明显分离，证明组织特异性代谢强且重复性好。
3.2. Analysis of Key Metabolites in P. corylifolia：以果实为对照的四组比较中，果实富集差异积累代谢物占优势；423个共有果实富集代谢物经富集分析指向苯丙氨酸、酪氨酸和色氨酸生物合成及香豆素、黄酮类。热图显示补骨脂苷(psoralenoside)与补骨脂辛(bakuchicin)偏根，苦参黄苷(swertiajaponin)与羽扇豆黄酮(luteone)偏叶，异补骨脂查耳酮、补骨脂黄酮、补骨脂素样黄酮偏果。
3.3. Transcriptomic Characterization of P. corylifolia：转录组得127.51 Gb clean data，80,402个Unigene，4681个转录因子分89科。与其他组织比果实的对比中差异基因众多，4311个共有差异基因富集于苯丙烷代谢、次生代谢生物合成、异黄酮与黄酮生物合成；差异转录因子以AP2/ERF-ERF、bHLH、MYB-related、bZIP居多。
3.4. Integrated Transcriptomic and Metabolomic Analysis：研究人员将Unigene映射到ko00940、ko00941、ko00943、ko00944等通路，显示PAL家系偏根，4CL在茎叶花较高，部分CHS与CHI在果中偏高且与黄酮积累趋势一致；黄酮亚支中F3'H成员表达不一，异黄酮支中大豆黄素(daidzein)与金雀异黄素(genistein)偏根、苜蓿酚(coumestrol)偏果；香豆素端TOGT1偏果，CYP82C4在花与果高，而补骨脂素与花椒毒素(xanthotoxin)偏根，说明基因表达与代谢物积累存在相关趋势。
3.5. WGCNA Co-Expression Network Analysis and Screening of Hub Transcription Factors：WGCNA将基因分为12模块，洋红色模块与果实及补骨脂素样黄酮、补骨脂定(psoralidin)、苜蓿酚等正相关。从该模块筛转录因子并与黄酮/香豆素结构基因建|r| ≥ 0.8且p < 0.05的网络，得344条边、6个转录因子：Cluster_22013.0(C3H)、Cluster_21217.8(NF-YA)为枢纽；Cluster_9299.0(Rcd1-like)、Cluster_20910.0(NAC)为正相关高连接候选；Cluster_31326.1(bZIP)、Cluster_11436.0(C2H2)仅负相关，符合负调控因子特征。
讨论部分总结：研究人员强调果实为药用部位却具最独特代谢与转录轮廓，与药用植物活性成分富集于药用器官的规律一致；苯丙烷途径在转录与代谢水平同步富集只说明相关趋势，不等同基因–代谢物因果。WGCNA仅给相关性而非因果，洋红色模块是果实次生代谢候选模块；本研究枢纽为C3H与NF-YA，不同于豆科常被关注的MYB，NF-YA可结合CHS等基因启动子CCAAT元件，拓宽了豆科果实次生代谢转录因子谱系；bZIP与C2H2全负相关问题提示负调控因子也是未来靶点。讨论也指明局限：无染色体级基因组、从头转录组可能漏注或并源错分；阻断配对采样未把池作块效应但主结论稳定；仅生信相关无RT-qPCR、酵母单杂交、双荧光素酶或遗传转化验证；代谢组无同时跑标准品且为相对定量；WGCNA样本仅15个、scale-free拟合R2=0.79略低于0.80，统计力有限。未来应做RT-qPCR、酵母单杂交、双荧光素酶、本氏烟草瞬时表达及补骨脂稳定转化。
结论部分翻译：本研究构建了补骨脂多组织代谢组与转录组数据集，系统刻画黄酮与香豆素的组织特异性积累模式。通过整合加权基因共表达网络分析与转录因子–通路基因共表达策略，研究人员将洋红色模块鉴定为与果实特征性活性成分富集相关的候选模块。该模块中6个候选转录因子按网络拓扑与模块内连接度分为三档：Cluster_22013.0(C3H)与Cluster_21217.8(NF-YA)为候选枢纽，Cluster_31326.1(bZIP)与Cluster_11436.0(C2H2)与通路基因仅呈负相关，符合其他物种已报负调控因子模式。本工作属发现阶段候选筛选，所有转录因子与模块均基于共表达相关，尚未实验验证，应视为未来功能验证的候选而非已确认调控因子。研究所提供的多组学数据、候选共表达模块与优先转录因子，为补骨脂活性成分生物合成与调控的后续功能基因组研究提供资源与框架。
