《Metabolites》：Integrated Transcriptomic and Metabolomic Analyses Reveal Organ-Specific Molecular and Metabolic Features in Nitraria sibirica
编辑推荐：
：背景/目的： 白刺（Nitraria sibirica Pall.）是一种具有生态和植物化学价值的盐生灌木，但其器官水平的分子组织尚不清楚。方法： 研究人员整合了来自同一田间种群采集的根、茎、叶和果实的配对RNA测序与非靶向液相色谱-质谱联用（LC-MS）代
：背景/目的： 白刺（Nitraria sibirica Pall.）是一种具有生态和植物化学价值的盐生灌木，但其器官水平的分子组织尚不清楚。方法： 研究人员整合了来自同一田间种群采集的根、茎、叶和果实的配对RNA测序与非靶向液相色谱-质谱联用（LC-MS）代谢组学数据。结果： 研究人员定量了25,263个基因和2143个注释代谢物条目，四个器官在两组学数据中均清晰分离。在六组成对比较中，鉴定出12,160个非冗余差异表达基因（DEGs）和1446个非冗余差异积累代谢物（DAMs）；交集筛选进一步解析出3955个器官偏好性基因和539个器官偏好性代谢物。叶片含有最大且最多样化的偏好性代谢物集，并表现出光合作用、叶绿体、色素和碳固定特征。茎以次生细胞壁、细胞骨架、碳水化合物转化、苯丙烷类和木脂素相关模式为特征。果实优先积累有机酸、氨基酸相关代谢物、糖类和糖苷，并表现出能量代谢和黄酮类修饰特征。根含有的偏好性代谢物较少，其特征是信号传导、保护性响应、RNA加工和表面屏障代谢。加权基因共表达网络分析（WGCNA）为每个器官鉴定出一个显著的器官相关模块，跨组学分析将这些模块与偏好性代谢物联系起来，同时优先筛选候选基因-代谢物关联。结论： 这些数据定义了白刺匹配的四器官分子和代谢分区，并为功能验证和靶向代谢物验证提供了可测试的候选对象。
论文解读文章
一、研究背景与科学问题
白刺（Nitraria sibirica Pall.）是蒺藜科（Nitrariaceae）白刺属的一种多年生盐生灌木，广泛分布于欧亚干旱区的草原、半荒漠、荒漠及盐碱生境中。该物种通过其多年生灌木结构和较强的环境适应性，在减少风蚀、稳定沙面、截留土壤颗粒及维持荒漠植被方面发挥重要生态功能，是研究盐生植物适应机制和干旱区生态恢复的潜在宝贵本土资源。生理学研究表明，适度盐处理可促进幼苗生长，Na+和K+转运调控有助于盐胁迫下的离子稳态。此外，白刺果实含有黄酮类、生物碱、酚酸、花青素和多糖等生物活性成分，具有抗氧化、抗补体、益生元等生物活性；叶片亦含有酚苷、苯丙烷类、木脂素和黄酮类化合物。然而，以往研究主要聚焦于胁迫处理、生境或发育对比、基因组资源及特定基因家族，尚未在同一田间条件下对根、茎、叶、果实四大主要器官进行匹配的分子分区解析。根和茎作为生态和结构上重要的器官，其化学和分子特征研究相对不足。此外，白刺根系还是荒漠药用寄生植物锁阳（Cynomorium songaricum）的重要寄主，解析根系特异性分子和代谢特征兼具胁迫适应生物学和干旱区药用植物生态背景的双重意义。因此，开展匹配的四器官转录组与代谢组整合分析，有助于明确器官间的组成性分子分区，为功能验证和靶向代谢物分析提供资源。
二、研究设计与技术方法
研究人员于2026年6月白刺盛果期，在新疆玛纳斯县六户地镇土炮营国家沙漠公园（44.707438° N, 86.152036° E）采集样本。选取9株生长状态相似、无可见病害的健康成熟植株，分为3个不重叠组，每组3株的等量鲜重混合为1个复合生物学重复，每器官3个重复，共12个样本。主要技术方法包括：非靶向超高效液相色谱-高分辨质谱（UHPLC-HRMS）代谢组学分析，采用C18色谱柱和正负离子模式全扫描，基于内标校正和自建谱库进行代谢物注释；RNA测序采用DNBSEQ-T7RS平台进行双端150 bp测序，经fastp质控、HISAT2比对至白刺参考基因组、StringTie组装和DESeq2差异表达分析；功能注释与富集分析采用GO、KEGG等数据库；加权基因共表达网络分析（WGCNA）构建共表达模块并关联器官性状；整合分析采用配对比较的基因-代谢物相关、联合KEGG通路分析及模块-代谢物关联。
三、主要研究结果
3.1 多组学数据集概览与质量评估。 RNA测序获得2.586亿条干净读段，Q30值96.13%–97.18%，基因组比对率90.00%–90.86%。代谢组学提取60,606个特征，经四级鉴定体系注释2143个条目，QC样本间Spearman相关系数0.9977–0.9985，内标RSD中位数0.64%，表明数据质量可靠。
3.2 不同器官的整体代谢谱。 主成分分析显示四器官代谢谱清晰分离，PC1和PC2累计解释51.2%方差。脂质和类脂分子为最大类别（37.9%），其次为有机酸及其衍生物（14.4%）、苯丙烷类和聚酮类（13.1%）。器官内相关性高于器官间相关性。
3.3 不同器官的整体转录组谱。 25,263个基因的PCA显示器官依赖性结构，PC1和PC2累计解释48.73%方差。表达趋势分析将19,053个基因分为9个簇，呈现茎、叶、根、果实特异性表达最大值。iTAK预测1303个转录因子，MYB家族最大（160个成员）。
3.4 器官间差异积累代谢物。 六组比较整合获得1446个非冗余DAMs，脂质类占33.8%。器官偏好性DAMs共539个：根35个、茎87个、叶303个、果实114个。叶偏好性代谢物集最大且化学多样性最高，果实偏好性代谢物以有机酸、氨基酸相关代谢物、糖类和糖苷为特征。
3.5 差异基因表达与功能富集。 六组比较鉴定12,160个非冗余DEGs，其中3955个器官偏好性基因（根708、茎385、叶1935、果实927）。光合作用是唯一在所有六组比较中均显著富集的KEGG通路。叶偏好性基因富集于光合作用、色素代谢和碳固定；茎偏好性基因富集于碳水化合物转化和结构代谢；根偏好性基因涉及信号传导和表面屏障代谢；果实偏好性基因涉及蛋白质加工、能量代谢和黄酮类生物合成。
3.6 器官相关基因共表达模块。 WGCNA将19,053个基因分配至38个模块，鉴定出4个显著器官相关模块：darkred（根）、darkolivegreen（茎）、blue（叶）、antiquewhite1（果实），共含11,192个基因。各模块功能特征分别为RNA加工与防御信号、细胞骨架与次生细胞壁、叶绿体与光合功能、能量代谢与黄酮生物合成。共筛选出32个根、26个茎、128个叶和53个果实候选枢纽基因。
3.7 转录组与代谢组整合分析。 配对整合筛选出高幅度基因-代谢物关联，联合KEGG分析鉴定出黄酮和黄酮醇生物合成、类胡萝卜素生物合成等协调富集通路。模块-代谢物关联显示darkred、darkolivegreen和blue模块与对应器官偏好性代谢物显著富集。枢纽基因-代谢物整合获得363个根、458个茎、7595个叶和546个果实阳性关联对，代表性关联如CCR1-落叶松脂素、PSBY-角黄素、CYP75B2-山柰酚-3-O-半乳糖苷等。
四、研究结论与意义
该研究首次在同一田间种群中建立了白刺根、茎、叶、果实四器官匹配的转录组和代谢组分子分区图谱。研究定义了各器官独特的分子和代谢特征：叶片以光合和色素代谢为主导，茎以结构代谢和苯丙烷类为特征，果实富集有机酸和氨基酸相关代谢物，根则以信号传导和表面屏障代谢为特征。整合分析优先筛选出大量候选基因-代谢物关联，为后续功能验证和靶向代谢物分析提供了可测试的候选对象。该研究不仅深化了对盐生灌木器官水平分子组织的认识，也为白刺作为干旱区生态恢复物种和药用植物资源的开发利用提供了重要的组学资源。
