整合转录组与代谢组分析揭示多花黄精生长发育调控的关键基因与代谢途径

《Frontiers in Plant Science》：Integrated transcriptome and metabolome analysis reveals key genes and metabolic pathways regulating growth and development in Polygonatum cyrtonema Hua

　　

摘要
多花黄精（*Polygonatum cyrtonema* Hua）是一种药食同源的传统中药，其药用成分具有重要的临床价值。由于市场需求量巨大，目前主要依靠人工栽培进行生产。为获得优质的多花黄精，研究人员对一年生和三年生多花黄精进行了转录组和代谢组测序，以探究其生长调控机制及品质改良相关的关键基因。本研究共鉴定出1,957个差异表达基因（differentially expressed genes, DEGs）和163个差异表达代谢物（differentially expressed metabolites, DEMs）。整合转录组与代谢组分析表明，多花黄精的生长调控可能主要与鞘脂代谢、苯丙烷生物合成以及淀粉和蔗糖代谢有关。数据提示，双向糖转运蛋白SWEET14表达升高可能促进蔗糖向库器官转运，蔗糖可能被β-呋喃果糖苷酶水解，一方面为植物生长提供能量，另一方面促进果糖积累。此外，L-苯丙氨酸丰度的升高可能与次生代谢物增加有关，这可能为根茎中年龄依赖性生长和代谢物分配提供代谢基础。研究人员观察到鞘脂代谢相关基因的下调可能反映了多花黄精的多年生生长习性，即第一年较慢的生长可能促进鞘脂介导的根系发育。然而，需要强调的是，这些推断基于转录组与代谢组的相关性数据，仍需功能性验证以建立因果关系。

多花黄精生长发育调控的整合组学解读：关键基因与代谢通路
作为药食同源的传统中药，多花黄精（*Polygonatum cyrtonema* Hua）的药用部位为根茎，含有丰富的多糖、甾体皂苷和黄酮等生物活性天然产物，其药理活性已被广泛研究，包括抗氧化、抗炎、保肝、抗动脉粥样硬化、抗疲劳、调节肠道菌群、改善化疗诱导的恶病质、抗抑郁及免疫调节等。随着市场需求的持续增长，野生资源已无法满足供应，且自然栖息地因土地开发而退化，目前主要依靠人工栽培。然而，栽培药材的品质和疗效往往不如野生资源。已有研究多集中在加工方法、活性成分提取分离、药理学及临床应用等方面，基于生长调控机制开展分子生物学层面的种质资源改良研究相对较少。已有证据表明，多花黄精中多糖、皂苷等活性成分的含量随植株年龄而变化，为确保药用品质，通常在其生长3至5年时采收。然而，调控其年龄依赖性生长与品质形成的分子机制尚不清楚，系统性地从转录和代谢层面阐明这一问题，对指导优质栽培和遗传改良具有重要意义。
基于上述背景，研究人员以一年生和三年生多花黄精为材料，依据植物“源–库”理论设计了整合转录组与代谢组分析方案：对作为光合源器官的叶片进行转录组测序，同时对作为主要库器官的根茎进行代谢组检测，以揭示叶片中基因表达变化与根茎中代谢物积累之间的关联。研究共鉴定出1,957个差异表达基因（DEGs）和163个差异表达代谢物（DEMs），并通过GO和KEGG富集分析，确定了鞘脂代谢、苯丙烷生物合成以及淀粉和蔗糖代谢三条与生长调控密切相关的核心代谢途径，提出了一批潜在的关键候选基因，为多花黄精种质改良提供了数据支撑和理论基础。
在关键技术方法方面，研究人员使用了以下主要技术手段：1）基于Illumina平台的转录组测序，对叶片样本进行mRNA测序，经从头组装获得unigenes并进行六大数据库（NR、GO、KEGG、Swiss-Prot、eggNOG、Pfam）功能注释；2）基于UHPLC-MS/MS的非靶向代谢组学方法，对根茎样本进行代谢物提取、分离与鉴定，结合主成分分析（PCA）和正交偏最小二乘判别分析（OPLS-DA）筛选差异代谢物；3）采用qRT-PCR对随机选取的10个DEGs进行表达验证，以18S rRNA为内参基因，确认转录组数据的可靠性。样本来源于中国安徽省池州市植物研究基地，由同一生态区遗传稳定的野生种群种子繁育，并在标准化栽培条件下生长，每个处理组设置三个生物学重复。
研究结果部分如下：
1. 不同年份多花黄参的基因变化。通过对一年生（A组）和三年生（C组）叶片进行转录组测序，分别获得40,249,718和40,671,892条高质量clean reads。经从头组装获得164,565条转录本和69,373条unigenes，平均长度900 bp，N50为1,286 bp。功能注释结果显示，NR数据库注释率最高（47.79%），其次为eggNOG（45.94%）。以|log2(fold change)|>1且FDR<0.05为标准，共鉴定出1,957个DEGs，其中上调728个、下调1,229个。GO富集分析显示，胞外区域、催化活性和膜组分显著富集。KEGG富集分析揭示DEGs显著富集于苯丙烷生物合成、淀粉和蔗糖代谢以及鞘脂代谢等通路，其中SWEET14（双向糖转运蛋白）、β-呋喃果糖苷酶和苯丙氨酸解氨酶（PAL）被确定为关键候选基因。qRT-PCR验证结果与转录组数据趋势一致，确认数据可靠。
2. 不同生长年份多花黄精的代谢组变化。基于UHPLC-MS/MS的非靶向代谢组学共注释到1,842个代谢物。PCA显示组内重复性好、组间分离明显。以p<0.05且VIP≥1为标准，共筛选出163个DEMs，其中上调82个、下调81个，差异主要集中在有机酸及其衍生物、苯环型化合物和有机氧化合物等类别。KEGG富集分析显示DEMs显著富集于不同环境中的微生物代谢、氨酰-tRNA生物合成、氰氨基酸代谢和氨基酸生物合成等途径。值得注意的是，L-苯丙氨酸在三年生根茎中显著积累，这一发现与转录组中苯丙烷生物合成通路的共富集相呼应，提示三年生植株中该通路可能被协同增强。
3. 转录组与代谢组整合分析。将DEGs和DEMs共同参与的代谢途径进行分析，共富集到鞘脂代谢、苯丙烷生物合成和淀粉与蔗糖代谢三条与生长调控相关的途径。九象限图相关性分析显示，L-苯丙氨酸的积累与苯丙烷生物合成途径相关基因表达呈显著正相关（Pearson相关系数>0.8），SWEET14和β-呋喃果糖苷酶等淀粉和蔗糖代谢基因形成明显的相关簇，鞘脂代谢途径中丝氨酸棕榈酰转移酶（SPT）和PAL等关键基因与差异代谢物形成相关性网络。
讨论部分指出，L-苯丙氨酸的显著积累可能是根茎成熟和防御的代谢适应性表现，为多种酚类化合物的合成提供底物；鞘脂代谢相关基因在三年生植株中的下调可能反映了多年生生长策略——第一年通过较慢的生长促进鞘脂介导的根系建植；SWEET14的上调和β-呋喃果糖苷酶的表达可能共同促进了蔗糖向根茎的运输和果糖积累。研究结论翻译如下：本研究调查了一年生和三年生多花黄精的叶片转录组和根茎代谢组，整合分析表明若干代谢途径和候选基因可能与年龄依赖性生长调控相关。具体而言，三年生根茎中L-苯丙氨酸丰度升高可能与次生代谢物多样性增加有关，可能为根茎成熟和代谢物再分配提供代谢基础；双向糖转运蛋白SWEET14可能促进蔗糖向库器官转运，β-呋喃果糖苷酶可能将蔗糖水解为果糖以供根茎积累；而鞘脂代谢相关基因在老龄植株中的下调可能反映了多花黄精的多年生生长策略。研究人员强调，所有这些发现均基于转录组与代谢组的相关性数据，候选基因的因果作用有待功能研究验证。该研究为深入探究多花黄精生长与代谢物积累的分子机制提供了数据支撑和理论基础。