《Horticulturae》：Integrated Transcriptomic and Metabolomic Analysis Reveals Metabolic Associations Underlying Oil Accumulation in Macadamia Cultivars
编辑推荐：
澳洲坚果（Macadamia integrifolia）种仁含油量可达80%，但品种间脂质积累差异的分子基础仍不清楚。研究人员选用4个油脂性状差异显著的澳洲坚果品种（O.C、QA1、QA2、QA3），整合转录组学、代谢组学与营养分析，探究与种仁油脂积累相关的代
澳洲坚果（Macadamia integrifolia）种仁含油量可达80%，但品种间脂质积累差异的分子基础仍不清楚。研究人员选用4个油脂性状差异显著的澳洲坚果品种（O.C、QA1、QA2、QA3），整合转录组学、代谢组学与营养分析，探究与种仁油脂积累相关的代谢通路。表型评价显示各品种果实性状与营养成分存在显著差异，O.C粗脂肪含量最高，QA1粗蛋白含量最高。转录组分析鉴定出大量差异表达基因（DEGs），主要富集于甘油酯代谢（glycerolipid metabolism）、α-亚麻酸代谢及氨基酸生物合成通路。代谢组学共检测778种代谢物，脂类相关化合物占较大比例。整合分析揭示3条与油脂积累相关的核心通路：氨基酸生物合成、甘油酯代谢与亚油酸代谢。氨基酸生物合成通路中，各品种间氨基酸差异积累伴随ASNS、GLUL、MAT等基因表达变化；甘油酯代谢中，TAG合成相关基因（DGAT、PDAT）差异表达，LysoPC、LysoPE等磷脂相关代谢物含量变异；亚油酸代谢中，脂肪酸去饱和与延长基因（SAD、FAD2、FAD3、FAE1/KCS）差异表达，多种脂肪酸相关代谢物丰度呈品种依赖性变化。该研究系统刻画了澳洲坚果种仁粗脂肪含量与脂肪酸相关代谢物的转录-代谢变异，并筛选出可供功能验证与分子育种的候选基因。
研究背景与立项依据
澳洲坚果为山龙眼科（Proteaceae）多年生乔木，种仁干重含油率可达80%，以单不饱和脂肪酸油酸（C18:1）为主，具高氧化稳定性与心血管健康效益。现有主栽品种间种仁含油量及脂肪酸组成差异显著，但品种依赖型油脂积累与脂肪酸组成的分子基础尚未阐明，尤其脂质生物合成、脂肪酸修饰及基因-代谢物关联层面的认识存在空白。常规育种依赖表型选择，受环境扰动大且童期长达7–10年。为突破表型选择瓶颈，研究人员开展本研究，旨在从多组学层面解析品种间油脂积累差异的转录与代谢基础，为分子育种提供候选资源。论文发表于《Horticulturae》。
主要技术方法概述
研究人员于贵州兴义万峰林澳洲坚果试验园采集4个品种（O.C、QA1、QA2、QA3）成熟种仁，每品种3株生物学重复、每株混样10果。营养品质采用索氏提取测粗脂肪、凯氏定氮测粗蛋白、蒽酮比色测可溶性糖、茚三酮比色测总氨基酸。转录组基于Illumina NovaSeq 6000测序，比对至M. integrifolia HAES 741染色体级参考基因组（GCA_013358625.1），DESeq2筛选DEGs（|log2FC|≥1，FDR<0.05）。靶向代谢组采用UPLC-MS/MS（MRM模式）检测，VIP≥1且|log2FC|≥1且p<0.05定义差异积累代谢物（DAMs）。转录-代谢整合将DEGs与DAMs映射KEGG，对氨基酸生物合成、甘油酯代谢、亚油酸代谢三条通路做Z-score热图与共表达（Pearson |r|≥0.80，FDR<0.05），并以RT-qPCR验证6个DEGs。
研究结果
3.1 品种间表型与营养特征
通过对四品种果实形态与种仁营养测定，研究人员发现QA1与QA3单果重与横径较大，QA2较小，O.C居中；O.C粗脂肪与可溶性糖最高，QA1粗脂肪最低但粗蛋白最高，QA3粗脂肪与总氨基酸均偏高。结果表明果径大小与种仁含油率非正相关，营养储备分配具品种特异性。
3.2 转录组质量与DEGs
12个文库Clean reads 42.01–50.72 M，Q30>98%，比对率90.40%–94.90%；PCA与相关系数>0.90显示重复可靠。以O.C为对照，QA1 vs O.C得1101 DEGs，QA2 vs O.C得3704 DEGs（最多），QA3 vs O.C得847 DEGs；356个DEGs三者共有。186个差异转录因子（TFs）中bHLH、ERF、MYB、bZIP、WRKY占比最高。
3.3 DEGs功能富集
GO富集指向甘油三酯代谢、中性脂质合成、酰基转移酶活性等；KEGG显著富集甘油酯代谢、α-亚麻酸代谢、氨基酸生物合成、糖酵解、植物激素信号转导等，提示脂质与中心碳代谢协同差异。
3.4 代谢组与DAMs
778种代谢物中脂类主导；QA1/QA2/QA3 vs O.C分别获233/272/202个DAMs，103个三者共有。KEGG富集共同通路含亚油酸代谢、氨基酰-tRNA生物合成等，亚油酸代谢富集印证脂肪酸修饰参与品质形成。
3.5 脂质相关通路整合
3.5.1 氨基酸生物合成：天冬氨酸族（ASNS、metE、metK/MAT）、谷氨酸族（GLT1、GLUL）、丝氨酸族（serA/PGDH、PSAT、cysE）及支链氨基酸模块（ilvB、ilvE）伴随天冬氨酸、天冬酰胺、谷氨酸、丝氨酸、缬氨酸等DAMs协同变异，提示氨基酸代谢通过碳骨架与乙酰-CoA前体供给关联油脂积累。
3.5.2 甘油酯代谢：LysoPC(16:0/18:1/18:2)、LysoPE(16:0/18:2)、胆碱、磷酸胆碱等DAMs品种间差异明显；GPAT、LPAT、PAP、DGAT、PDAT及PLA、PLC、PLD差异表达，DGAT与PDAT表达模式与粗脂肪表型对应，提示其为粗脂肪变异候选。
3.5.3 亚油酸代谢：硬脂酸（C18:0）、油酸（C18:1）、亚油酸（C18:2）、γ-亚麻酸（C18:3）丰度品种间不同；SAD、FAD2、FAD3、FAE1/KCS差异表达且与对应脂肪酸积累趋势一致，表明去饱和与延长步骤的转录调控塑造脂肪酸组成。
3.6 RT-qPCR验证
6个选定DEGs的RT-qPCR与RNA-seq表达趋势一致（R2=0.924–0.995），佐证转录组可靠性。
讨论部分总结
研究人员指出果径与含油率解耦、储藏物质分配非简单此消彼长，而受碳氮分配与代谢能力调控。氨基酸通路通过TCA（三羧酸循环）节点供给前体关联油脂积累；甘油酯通路中DGAT、PDAT及LysoPC/LysoPE重塑可能参与TAG（三酰甘油）装配但并非直接因果；SAD-FAD2轴决定C18:1/C18:2比例，FAD3与FAE1/KCS体现高不饱和与超长链分支的品种策略。作者明确本研究局限为单时点单点采样、关联非因果、缺酶活与通量验证，需发育时序与基因功能实验补证。
结论部分翻译
本研究整合营养、转录组与代谢组分析，刻画了四个澳洲坚果品种的种仁油脂相关分子与代谢变异。O.C粗脂肪含量最高，三对比较中基因与代谢物均呈显著品种间分化。整合分析确定氨基酸生物合成、甘油酯代谢与亚油酸代谢为基因-代谢协同变异的核心通路；氨基酸代谢、甘油酯代谢及脂肪酸去饱和延长相关基因的差异表达，伴随LysoPC、LysoPE、C18:0、C18:1、C18:2、C18:3等代谢物变化。研究识别出与品种间粗脂肪及脂肪酸相关代谢物变异关联的候选基因、代谢物与通路，但未直接证实代谢通量、脂质周转或候选基因因果功能，后续需功能、酶活、发育与多环境研究加以验证。
