马铃薯生长-防御资源分配响应的多设施高时间分辨率表型组学研究

《Plant Phenomics》：Multi-facility, high-temporal-resolution phenomics of growth–defence allocation responses in potato

　　

理解植物如何在生长与防御之间分配资源是植物生物学的一个基本问题。此类分配响应通常较为细微，需要精密的测量工具和比较性实验设计方能检测。本研究在两家独立设施中部署了每日监测的高通量表型组学框架，在无病原体条件下评估马铃薯（Solanum tuberosum）品种Désirée及其六个转基因株系的生长、形态与生理特性。研究人员对经抗性诱导化合物β-氨基丁酸（BABA）或亚磷酸钾（KPhi）处理的植株，以及携带单个Rpi基因（R2A、blb1）、三基因叠加（3R1、3R2）或通过RNAi敲低茉莉酸共受体StCOI1（coi1H1、coi1X5）的株系进行了表型分析，以比较抗性诱导剂与基因编码抗性在分配响应上的差异。在四周内，图像衍生性状揭示了两设施间一致的效果：BABA增加了株高并提高了绿色指数，而KPhi未对生长产生一致性改变；StCOI1沉默株系在株高上有所增加。终点酶谱分析显示，BABA和KPhi在各自设施内激活了不同的酶类，且KPhi始终表现出较低的表观生长成本；StCOI1沉默株系尽管株高表型一致，却表现出设施依赖性的酶构型；然而，具体的酶特征并未在两设施间保守，表明跨设施酶表型分析需明确表征当地代谢基线。总体而言，高通量表型组学在跨平台间解析了与抗性相关的细微分配响应，表明BABA和KPhi在特定环境下与不同的终点酶特征相关，其中KPhi的较低成本在环境中保守，但特征身份由处理和环境共同塑造。

一、研究背景与问题

植物在生长与防御之间的资源分配权衡是植物生物学和作物改良领域的根本性问题。植物可采用不同的防御策略，从代谢成本高昂的组成型防御到更为节能的诱导型防御。细胞表面定位的模式识别受体（PRRs）和胞内核苷酸结合富亮氨酸重复受体（NLRs，亦称R基因）分别激活PAMP触发免疫（PTI）和效应子触发免疫（ETI），两者相互依赖并可触发系统获得性抗性（SAR）。过去五十年的研究表明，防御与生长之间的权衡常由基因多效性效应和植物激素间的拮抗交叉对话引起。大多数R基因在感染期间并未被强烈转录诱导，表明与抗性相关的任何适合度成本更可能源于组成型受体维持和/或病原识别后触发的下游防御响应。理解哪些下游过程驱动这些成本，以及这些成本是机制特异性的还是共享的，直接关系到R基因堆叠和易感（S）基因利用等作物改良策略。此外，外源化合物如β-氨基丁酸（BABA）和亚磷酸钾（KPhi）可作为抗性诱导剂（RIs）激活诱导抗性（IR），但其相关的适合度成本难以量化。一个尚未充分探索的问题是，不同防御机制（化学诱导剂与基因编码抗性）是否参与相同或不同的代谢再分配途径。由于碳水化合物代谢酶位于碳通量的关键分支点，其活性特征可作为资源如何在生长与防御之间分配的生化读数。跨环境研究中，当地条件可能独立调节碳代谢，即使宏观表型在各地点一致，也可能产生环境特异的酶生物特征。因此，本研究利用两家北欧植物表型设施，结合时间序列成像与化学抗性诱导剂及机制多样的转基因马铃薯株系，在无病原体条件下直接比较抗性诱导剂与遗传抗性机制是否施加相似的适合度成本并参与相同的代谢再分配途径。

二、研究内容与结论

研究人员在芬兰赫尔辛基大学国家植物表型基础设施（NaPPI）和丹麦哥本哈根大学PhenoLab两个高通量植物表型设施中开展了实验。以马铃薯品种Désirée为共同参照，评估了六个转基因株系：两个过表达单个Rpi基因的株系（R2A、blb1）、两个携带三个Rpi基因堆叠的株系（3R1、3R2）以及两个StCOI1 RNAi敲低株系（coi1H1、coi1X5）。每周对Désirée植株喷施BABA或KPhi共四次。通过图像传感器追踪冠层面积、株高、形态性状和冠层颜色达四周，终点测量包括地上部生物量、块茎产量以及碳水化合物和抗氧化代谢关键酶的活性特征。研究发现BABA处理在两个设施中均显著增加株高并提高绿色指数，同时伴随冠层面积减小和干重降低；KPhi处理未对生长产生一致性改变，但表现出较低的表观生长成本。StCOI1沉默株系在两个设施中均表现出株高增加，但酶构型呈设施依赖性。酶谱分析显示BABA和KPhi在每个设施内与不同的酶类相关，但具体酶特征未跨设施保守。该研究首次在农业相关作物中通过多设施高通量表型组学框架系统解析了生长-防御分配响应，表明宏观生长-防御权衡在不同环境中可能由不同的代谢再分配途径支撑，为抗性诱导剂和R基因堆叠策略的作物应用提供了重要的表型学依据。

三、关键技术方法

本研究采用的主要技术方法包括：广义加性混合模型（GAMMs）用于分析时间序列表型数据，捕捉非线性时间趋势和组间特异性轨迹；RGB成像与多光谱成像系统分别用于获取冠层面积、株高和归一化差异植被指数（NDVI）等性状；叶绿素荧光（ChlF）成像用于评估光系统II最大量子产率（F_v/F_m）和荧光衰减比（Rfd_Lss）等生理参数；终点酶活测定涵盖八种抗氧化酶和十二种碳水化合物代谢酶；主成分分析（PCA）用于酶谱数据的降维与可视化。样本队列包括NaPPI和PhenoLab各54株植物，每处理或株系六个生物学重复。

四、研究结果

诱导抗性（IR）处理结果：BABA处理植株在两个设施中均显著高于对照，株高差异在NaPPI达+86.3 mm（38 DAT），在PhenoLab达+131 mm（34 DAT）；KPhi处理无显著差异。冠层面积方面，BABA处理植株在20-35 DAT显著小于对照，而KPhi处理在22-25 DAT显著大于对照。绿色指数方面，BABA处理植株在NaPPI（23-38 DAT）和PhenoLab（21-34 DAT）均显著高于对照。叶绿素荧光测量显示KPhi处理显著提高有效光系统II量子产率（QY_Lss），BABA处理显著降低Rfd_Lss而KPhi处理显著提高Rfd_Lss。酶学分析显示两种诱导剂在各自设施内形成不同聚类，但驱动分离的关键酶在两设施间完全不同。生物量方面，BABA处理在NaPPI显著降低干重，而KPhi处理无显著差异。

转基因株系结果：株高方面，blb1在两个设施中均显著低于对照；3R2在NaPPI（34-38 DAT）显著高于对照；coi1H1和coi1X5在两个设施中均在生长后期显著高于对照。冠层面积方面，coi1H1在19 DAT至生长结束显著低于对照。酶学分析显示coi1H1和coi1X5在两个设施中均紧密共聚类，但驱动因子完全不同（NaPPI为PFK，PhenoLab为AGPase、UGPase和PGM）。生物量方面，NaPPI中coi1X5干重显著降低，3R1块茎数显著增加，coi1X5块茎重显著增加；PhenoLab中携带Rpi基因的株系（3R1、3R2、R2A、blb1）干重显著降低，而coi1H1鲜重和干重均显著增加。

五、讨论与结论

讨论部分指出，高通量表型组学能够解析抗性相关分配响应的细微变化，BABA和KPhi在特定环境下与不同的终点酶特征相关，KPhi的较低成本在环境中保守但特征身份由处理和环境共同塑造。跨设施联合PCA显示设施相关的聚类，NaPPI和PhenoLab对照植株占据不同的代谢位置，表明当地代谢基线的差异。研究结论为：高通量表型组学在跨平台间解析了细微的抗性相关分配响应，表明BABA和KPhi在给定环境中与不同的终点酶特征相关，KPhi的较低成本在环境中保守，但特征身份由处理和环境共同塑造；StCOI1沉默株系尽管株高表型一致，却表现出设施依赖性的酶构型，表明跨设施酶表型分析需要明确表征当地代谢基线。