Brassica rapa田间站点特异性蛋白质组重塑揭示跨纬度气候适应与作物可塑性

《BMC Plant Biology》：Field site-specific proteome remodeling under contrasting environmental conditions in Brassica rapa

背景：植物可塑性使植物能够适应异质环境，但作物长期气候驯化的分子机制仍认识不足。在此，研究人员调查了Brassica rapa在三个对比田间站点中的蛋白质组可塑性，这些站点代表亚北极（68°N）、温带（59°N）和亚热带（44°N）气候相关条件。研究人员栽培了涵盖三个作物类型的12份种质，并进行蛋白质组分析（proteome analysis）和蛋白质-代谢物相关分析（protein-metabolite correlation analysis）。结果：主成分分析（PCA）和聚类分析显示田间站点对蛋白质组模式具有强烈印记。功能富集分析（functional enrichment analysis）鉴定出站点特异性蛋白网络：亚北极条件有利于参与光合作用（photosynthesis）、碳水化合物代谢（carbohydrate metabolism）和糖酵解（glycolysis）的蛋白富集，支持能量捕获和冷驯化（cold acclimation）；温带条件诱导与呼吸作用（respiration）和氧化还原稳态（redox homeostasis）相关的蛋白；亚热带条件诱导参与分支特异性代谢（clade-specific metabolism）、尤其是硫苷生物合成（glucosinolate biosynthesis, GLS）的蛋白，并伴随氨基酸和氧化还原通路。基于蛋白质组数据，对三个田间站点的关键适应响应包括细胞呼吸（cellular respiration）蛋白如NADH脱氢酶亚基（NADH dehydrogenase subunits）的调节，以及植物激素（phytohormone）调控的变化：细胞分裂素（cytokinins, CKs）在亚北极富集，赤霉素调控蛋白（gibberellin-regulated proteins）在温带富集，生长素相关蛋白（auxin-associated proteins）在亚热带富集。大白菜（Chinese cabbage）种质A247表现出最高蛋白质组可塑性，而根用芜菁（root turnip）种质表现出较低可塑性。沿纬度气候梯度的反应规范分析（reaction-norm analyses）揭示平均蛋白Z分数（protein Z-scores）的基因型和作物类型依赖轨迹，显示不同栽培品种间存在对比性适应策略。结论：B. rapa在不同田间站点表现出强蛋白质组灵活性（proteomic flexibility），表明不同的代谢调整。这些发现提示蛋白质组谱分析（proteomic profiling）有助于鉴定适应性状，并可指导作物育种策略以提高未来气候变化下的韧性（resilience）。

研究背景：作物适应未来环境是作物研究的重要挑战，高温及伴随干旱是限制产量的主要环境因素。应用组学技术虽加深了对热胁迫响应分子基础的理解，但耐热品种培育进展有限，原因包括性状的多基因复杂性，以及大量数据来自温室或生长箱的短期胁迫处理，这种简化方法阻碍了知识向田间的转化。田间植物持续受到胁迫，并常同时暴露于多种胁迫约束；尽管田间试验劳动强度、成本和环境波动更大，但只有纳入真实田间适应作物，才能连接分子机制与田间表现。纬度梯度可在自然条件下提供大时间与地理尺度信息，是研究植物适应变化环境的方法学设置。长期以来可塑性多在表型水平研究，少数扩展至转录组；移植实验可检验可塑性范围和限度。Brassicaceae约占世界蔬菜产量10%，Brassica rapa包括大白菜、小白菜和芜菁等。其消费与癌症、心血管和神经退行性疾病风险负相关，保护作用常归因于硫苷（GLS，Brassicales几乎特有的次生代谢物）降解产物及类胡萝卜素、多酚等。此前100份B. rapa种质在2021年三个站点（亚北极、温带、亚热带）种植并鉴定抗虫种质；本研究用覆盖三个亚种的12份核心种质在相同条件下研究纬度梯度上蛋白质组可塑性范围。研究目的是鉴定田间条件下全局生理适应并检查哪些蛋白贡献站点特异性驯化，为预测气候变化对作物影响及改良生产系统韧性和产量提供分子背景。

研究人员开展的研究及意义：研究人员于2021年在三个田间站点种植12份B. rapa种质，进行叶片蛋白质组分析并与平行代谢物数据整合，鉴定站点相关蛋白质组重塑、关键通路和种质可塑性差异；结果表明B. rapa在亚北极、温带和亚热带环境中表现显著蛋白质组可塑性，细胞呼吸、植物激素代谢和GLS生物合成发生气候依赖重编程，为识别适应性状和指导气候韧性育种提供依据。

关键技术方法（不超过250字）：研究人员于2021年在三个站点（亚北极68°58′N、温带59°43′N、亚热带44°36′N）开展田间多地点实验，材料为12份B. rapa种质，覆盖大白菜7份、小白菜2份、根用芜菁3份，来自VIR种质库，三站点各设3个生物学重复、每样本5株叶片混合；采用nanoflow LC-MS/MS进行蛋白质组分析、数据库搜索与定量，并结合平行代谢物数据进行蛋白-代谢物相关和多组学整合；GLS水解产物用GC-MS分析；另在控制条件下对小白菜设置灌溉梯度并用UHPLC-DAD-ToF-MS分析GLS。

研究结果：
田间站点环境变异：通过记录2021年生长季三站点温度、降水、光周期和成熟期，发现温度和降水随纬度升高而降低，光周期亚北极最长（22–24 h），温带18 h，亚热带15 h，成熟期分别约8、12、14周；虫害评分显示P. xylostella叶损温带最高、亚热带次之、亚北极最低，M. brassicae温带最高、亚热带次之、亚北极无损伤；作物类型偏好为趋势。由此得出站点环境与植食者压力存在差异。

站点相关条件诱导B. rapa叶片蛋白质组变化：通过叶片蛋白质组分析鉴定并定量5640个蛋白；主成分分析（PCA，降维可视化）显示站点强影响，前两主成分解释30%变异；偏最小二乘判别分析（PLS-DA，监督分类）验证站点分离，Q2=0.873，置换p=0.024；联合PCA与代谢物显示协调分离；K-means聚类（按表达模式分组）得到10个表达簇；功能富集显示亚北极富集光合、细胞呼吸、碳水化合物代谢和蛋白生物合成，温带富集细胞呼吸、蛋白转运和氧化还原稳态，亚热带富集分支特异代谢、氧化还原和氨基酸代谢；关键蛋白与下游代谢物相关。由此得出站点条件驱动蛋白质组重编程。

细胞呼吸相关蛋白对环境的适应：基于MapMan注释鉴定198个细胞呼吸相关蛋白；41个糖酵解蛋白中胞质糖酵解在亚北极高丰度，甲基乙二醛/乙二醛降解在亚热带高；11个丙酮酸氧化和41个三羧酸循环（TCA cycle，线粒体碳代谢）蛋白多呈种质特异；氧化磷酸化（oxidative phosphorylation，线粒体ATP生成）复合体I的多个NADH脱氢酶（NADH dehydrogenase）亚基在温带高，亚热带B8高，亚北极铁硫蛋白4高；复合体III2、IV和V无站点依赖。由此得出呼吸链调节具站点特异性。

植物激素代谢蛋白反映站点适应：激素相关蛋白在各站点富集；脱落酸（ABA，胁迫响应激素）蛋白在亚北极和亚热带较高，PYR1亚北极特异或高；乙烯蛋白亚北极较低；茉莉酸（JA，防御激素）亚热带略高；细胞分裂素（CK，延缓衰老等）亚北极高；赤霉素（GA，促生长）温带高；油菜素内酯（BR，胁迫耐受）仅部分种质亚热带积累；生长素（IAA，调控生长）亚热带最高。由此得出乙烯、细胞分裂素、赤霉素和生长素等响应长期站点条件。

硫苷代谢相关蛋白受站点条件差异影响：GLS代谢蛋白亚热带富集；脂肪族GLS相关15个蛋白温带最高，吲哚/苯环族GLS相关14个蛋白亚热带最高；GLS代谢物分布与蛋白一致，脂肪族GLS在亚北极和温带为主，吲哚/苯环族GLS亚热带为主；磺化相关8个蛋白亚热带最高、亚北极次之、温带最低；31个GLS分解蛋白无站点特异。由此得出GLS组成受站点影响而分解较稳定。

灌溉不影响小白菜GLS积累：在控制条件下对小白菜设置水涝及100%、80%、60%、40%灌溉处理，检测9种GLS（5种脂肪族、3种吲哚、1种苯环族），总GLS为7.9±2.3至10.4±0.6 μmol g^?1干物质，总GLS和三类GLS均无显著差异。由此排除水分对GLS积累的显著影响。

根用芜菁种质蛋白质组可塑性低于叶菜种质：通过计算种质跨站点蛋白质组可塑性，范围0.660–0.741；大白菜A247（0.741）、A270（0.732）和小白菜A195（0.722）较高；根用芜菁A1300（0.660）、A1424（0.669）、A1409（0.676）最低。由此得出作物类型间存在可塑性趋势。

讨论与结论：讨论指出田间多站点不能完全分离温度、土壤、微生物和光等单因素，需每类气候设置多个站点；但与短期冷/热胁迫研究无普遍重叠，可能反映长期发育中的代谢重编程。亚北极糖酵解和线粒体功能增强可补偿低温酶活和ATP产生，NADH脱氢酶铁硫蛋白4与冷适应相关；激素串扰中乙烯与细胞分裂素呈负相关，ABA/PYR1、赤霉素与ABA负互作，生长素促生长；GLS总量和组成受温度/光周期影响，脂肪族与吲哚GLS和昆虫偏好相关，灌溉无显著影响；根用芜菁低可塑性可能反映储藏器官缓冲。研究结论译为：本研究表明B. rapa在亚北极、温带和亚热带环境中表现显著蛋白质组可塑性，反映细胞呼吸、植物激素代谢和GLS生物合成的气候依赖重编程；冷适应伴随糖酵解和线粒体功能相关蛋白丰度增强及细胞分裂素和ABA信号转变；温带条件促进赤霉素调控蛋白丰度及生长过程；亚热带条件诱导生长素相关响应和吲哚GLS积累；蛋白质组谱分析为解析地理梯度植物适应分子机制提供假设并有助于气候韧性育种。