《Plants》：Sucrose-Related Transcriptomic and Metabolic Profiles Are Associated with Latex Yield Variation in Three Hevea brasiliensis Clones
编辑推荐：
摘要专业翻译
来自巴西橡胶树（Hevea brasiliensis）的天然橡胶是一种不可替代的工业原料。橡胶产量潜力是否受蔗糖库—蔗糖分配限制仍未完全阐明。本研究整合了来自三个产量对比鲜明（分别命名为H、M和L）的橡胶树无性系的叶片、树皮和胶乳的转录组学与靶向
摘要专业翻译
来自巴西橡胶树（Hevea brasiliensis）的天然橡胶是一种不可替代的工业原料。橡胶产量潜力是否受蔗糖库—蔗糖分配限制仍未完全阐明。本研究整合了来自三个产量对比鲜明（分别命名为H、M和L）的橡胶树无性系的叶片、树皮和胶乳的转录组学与靶向糖代谢组学数据。转录组分析揭示了组织特异性功能分化：叶片富集光合碳固定相关通路，树皮富集物质运输相关通路，胶乳富集防御与再生相关通路。与L无性系相比，H无性系并未表现出光合通路转录富集程度的增强；相反，H无性系表现出更高的转化酶（invertase）和蔗糖合酶（sucrose synthase）基因转录丰度，同时伴随己糖水平升高，提示蔗糖水解与己糖周转可能存在基因型差异。代谢组学分析显示，L无性系叶片中蔗糖积累量较高，而胶乳中蔗糖浓度相对较低。加权基因共表达网络分析（WGCNA）鉴定出MEpink模块与叶片蔗糖含量高度相关。该模块富集于硫胺素代谢和氰氨基酸代谢通路，并包含候选转录因子WRKY40、MYC2和NAC002。研究结果表明，这些无性系间的产量差异与不同的蔗糖分配和利用模式相关。这些多组学特征为高产橡胶树育种提供了分子资源。
论文解读
一、研究背景与科学问题
天然橡胶（顺式-1,4-聚异戊二烯）凭借其独特的弹性、耐磨性和抗撕裂强度，在航空航天、交通运输、医疗防护和国防等关键领域具有不可替代的战略地位。与合成橡胶相比，天然橡胶具有更优异的动态力学性能，全球市场需求持续攀升。目前，巴西橡胶树（Hevea brasiliensis）几乎是天然橡胶唯一的商业化来源，种植区域主要集中在东南亚、非洲和南美洲的热带地区。然而，橡胶生产面临多重制约因素，包括病害、遗传变异有限和气候适应性受限等。作为一种多年生热带乔木，橡胶树育种周期长达30至40年，加之近几十年气候变率和病害频发，全球天然橡胶供需失衡问题日益加剧。因此，深入阐明橡胶树产量形成的分子调控机制，培育高产优质新品种，对保障全球天然橡胶供应链的安全稳定具有重要的理论价值和现实意义。
蔗糖在植物生理中具有核心地位，是叶片通过韧皮部向茎、果实和根等库组织转运的主要糖类运输形式。橡胶树在割胶采胶过程中承受剧烈的生理胁迫——胶乳含有40%的干胶颗粒，采胶形成了人工碳水化合物库，与维持、生长和繁殖所需的天然碳水化合物储备形成竞争。研究表明，割胶可增加橡胶树树干内的淀粉储备，但以牺牲径向生长为代价。可溶性糖尤其是蔗糖，在树皮中充当胶乳生产的中间介质，在木材中充当淀粉合成的底物。蔗糖浓度与胶乳产量直接相关，介导蔗糖进入乳管细胞的蔗糖转运蛋白基因HbSUT3的表达在乙烯处理和割胶机械刺激下均上调，而碱性/中性转化酶将蔗糖裂解为己糖，供后续胶乳生产利用。
尽管前人研究已探索橡胶产量形成的机制，但大多数研究集中于单一组织（主要是胶乳）或单一组学层面，对叶片、树皮和胶乳之间协调的糖代谢系统表征仍然有限。针对这一空白，研究人员选取了三个遗传背景不同、田间实测胶乳产量存在显著差异的橡胶树无性系（命名为H、M和L），通过整合转录组学和靶向糖代谢组学，结合加权基因共表达网络分析（WGCNA），系统考察了组织特异性蔗糖分配及其潜在的分子调控网络。
二、主要关键技术方法
研究人员从中国农业部景洪橡胶树种质资源圃的80份橡胶树种质中筛选出三个产量经单因素方差分析（ANOVA）和Tukey HSD检验确认存在统计学显著差异的无性系，标记为H、M和L。所有树木于2008年定植，2022年8月初采用s/4d3割胶制度开始割胶并持续三个月，于2022年8月、9月和10月末用量筒测定胶乳产量。每个无性系每种组织设置三个生物学重复，共采集27份样本（叶片、树皮、胶乳）。研究采用了以下核心技术手段：基于Illumina HiSeq? 2000平台的RNA测序（RNA-Seq）转录组分析，以FPKM定量基因表达水平，使用DESeq包进行差异表达分析；基于气相色谱-质谱联用（GC-MS）的靶向糖代谢物检测，采用内标法和外标校准曲线进行绝对定量；加权基因共表达网络分析（WGCNA）构建共表达模块并关联模块特征基因与糖代谢物；以及实时荧光定量PCR（qRT-PCR）验证关键基因表达。
三、研究结果
3.1 不同产量橡胶树无性系的比较转录组分析
研究人员对H、M和L三个无性系的树皮、叶片和胶乳进行转录组测序，构建了27个cDNA文库。主成分分析显示不同组织样本间清晰分离，聚类热图进一步表明叶片整体转录活性最高。采用BH校正p<0.05且|log?(倍数变化)|>1的严格阈值，鉴定出1674个统计学确认的差异表达基因（DEGs）；另以未校正名义p<0.05和|log?(倍数变化)|>1获得12500个探索性候选基因用于下游假设生成分析。京都基因与基因组百科全书（KEGG）富集分析显示组织特异性分布模式：叶片候选基因显著富集于光合作用、碳代谢和糖酵解/糖异生通路；树皮主要富集蛋白质降解与运输、次生代谢物生物合成和细胞壁代谢；胶乳主要富集类黄酮生物合成、氨基酸代谢和硫代谢。这些结果表明叶片的光合碳同化特征塑造源特性，树皮呈现与物质运输和代谢相关的特征，胶乳则表现与防御和再生一致的表达谱。
3.2 靶向糖组学揭示组织和无性系间糖谱差异
研究人员采用基于GC-MS的广泛靶向代谢组学分析，在叶片、树皮和胶乳中共检测到27类碳水化合物代谢物，鉴定出24个差异丰度代谢物（DAMs）。质控样本中超过80%的峰变异系数小于0.3，确认分析平台重复性良好。检测到的糖代谢物中单糖占比最高（78%），二糖和三糖分别占18%和4%。层次聚类显示叶片富含多种糖类，而树皮和胶乳糖含量普遍较低。蔗糖浓度方面，各组织内三个无性系间经多重检验校正后均未达到统计学显著差异，但存在一致的数值趋势：L无性系叶片和树皮中蔗糖数值最高，M无性系胶乳中蔗糖浓度数值最高。
3.3 WGCNA鉴定与糖代谢调控相关的共表达模块
研究人员通过WGCNA分析共鉴定出15个共表达模块。模块-样本相关性分析显示显著的組織特异性，其中MEpink、MEtan、MEblue、MEgreenyellow、MEgreen和MEpurple与叶片样本正相关而与树皮和胶乳负相关。模块-代谢物关联分析显示多个模块与关键糖代谢物显著正相关。在全27个样本中，MEpink和MEtan与蔗糖含量正相关（r=0.76和r=0.63）。考虑到组织身份可能造成的混杂，研究人员进行了组织分层和偏相关分析：叶片内MEpink与蔗糖仍显著相关（r=0.82，BH校正p=0.038），而MEtan在叶片内相关性中等但不显著（r=0.62，p=0.072）；控制组织和无性系后，MEpink（偏r=0.76，p=2.7×10??）和MEtan（偏r=0.59，p=0.003）均与蔗糖保持正相关。这表明MEpink是稳健的叶片相关蔗糖关联模块。
3.4 核心模块的功能与表达分析
对MEpink和MEtan模块基因进行KEGG通路富集分析，结果显示这些基因显著富集于硫胺素代谢和氰氨基酸代谢通路。表达热图显示这些基因在叶片中整体高表达，且多数基因在L无性系叶片中的表达强度略高于H和M无性系。通过将各组织无性系间差异表达候选基因与MEpink和MEtan模块基因取交集，叶片中鉴定出71个重叠基因，包括多个候选转录因子（如WRKY40、WRKY46、WRKY53、MYB44、MYC2、NAC002）和胁迫响应基因（如HSP70、HSP83A、LOX3.1），这些基因普遍在L叶片中高表达；树皮中鉴定出62个重叠基因，多数在M树皮中高表达；胶乳中鉴定出21个重叠基因，包括WRKY53、MYB44和NAC002等候选转录因子，在H和M胶乳中表达较高。
3.5 蔗糖代谢通路的整合多组学分析
研究人员整合转录组和代谢组数据系统分析了蔗糖代谢通路。基因表达层面，醛缩酶编码基因（ALFP3、ALFC2、ALFC6、ALFP2）在H和M叶片中高表达，蔗糖磷酸合酶和蔗糖磷酸磷酸酶基因（SPP1、SPSA1、SPSA3、SPSA4）在H叶片和H胶乳中高表达，而转化酶基因（INVE、INVC、INVA、CINV2）和蔗糖合酶基因（SUS、SUS2、SUS6）在H和M无性系的叶片、树皮和胶乳中均高表达。代谢物层面，L无性系叶片和树皮中蔗糖浓度数值较高，M无性系胶乳中蔗糖数值较高；葡萄糖和果糖在H叶片中更丰富，在L树皮和L胶乳中也升高。L无性系树皮和胶乳中积累的葡萄糖和果糖未被消耗，提示L无性系可能并非缺乏己糖供应，而是下游己糖加工发生改变。
3.6 qRT-PCR验证
研究人员选取三个糖代谢关键基因进行qRT-PCR验证。所有基因均表现组织特异性表达，叶片中转录丰度最高，树皮居中，胶乳最低。基因23505和基因22397在叶片中表现显著表达差异，基因16100在树皮中检测到显著差异，与RNA-Seq数据一致。qRT-PCR结果与转录组数据总体吻合，验证了RNA-Seq结果的可靠性。
四、讨论与结论
本研究表明，橡胶树各组织存在明确的功能分化：叶片富集光合碳固定通路，树皮富集物质运输通路，胶乳富集防御和再生通路。H无性系并未比L无性系表现出更高的光合通路富集程度，说明这些基因型间的产量差异与碳分配模式和基因表达谱差异相关。L无性系叶片中蔗糖数值较高但蔗糖裂解酶表达相对较低，而H和M无性系转化酶和蔗糖合酶表达较高且伴随叶片蔗糖较低和己糖积累较高，这一协调趋势与H无性系增强蔗糖水解为己糖以为橡胶生物合成提供碳前体的假说一致。M无性系胶乳中蔗糖浓度最高但产量表现不佳，说明高蔗糖底物丰度本身不足以驱动高橡胶生物合成。MEpink模块与叶片蔗糖含量稳健相关，其富集的硫胺素代谢通路可能参与碳代谢平衡的补偿性调控。候选转录因子WRKY40、MYC2和NAC002在L叶片中上调，与碳胁迫响应的推测一致。研究结论指出，这些多组学数据揭示了基因型产量变异与不同蔗糖分配和利用特征之间的相关性，为高产橡胶树育种提供了可检验的假说和候选分子资源，但驱动产量变异的因果机制和基因功能仍需后续功能实验验证。
