本文描述了一种非破坏性、高通量的方法,用于自主进行田间光系统II量子效率、光谱反射率和植物构型的测量,从而实现农艺学和育种田间试验中大规模冠层光合作用表型分析。
本文描述了一种非破坏性、高通量的方法,用于自主进行田间光系统II量子效率、光谱反射率和植物构型的测量,从而实现农艺学和育种田间试验中大规模冠层光合作用表型分析。
光合作用不仅为植物生物量的生成提供能量,也为固氮等共生过程提供能量。尽管通过改善光能截获和收获指数来进一步提高主要作物产量的遗传潜力已基本耗尽,但自然存在或人工诱导的光合性状遗传变异仍为产量的进一步提升提供了可观的潜力。然而,由于在不断变化的田间条件下光合作用具有高度动态性,除非具备高空间和时间分辨率的数据,否则难以对光合性能进行有针对性的选择。为弥补这一空白,我们将一种光诱导荧光瞬变(LIFT)装置安装在一台自主田间机器人上,用于测量光系统II的量子效率(Fq'/Fm'),该参数已被证明与整体光合性能高度相关。LIFT方法利用高频重复的亚饱和闪光诱导最大荧光,可在最远1米的距离上在不到1毫秒内完成测量。该机器人以0.5 m s-1的速度行进,基于全球导航卫星系统(GNSS)坐标自主导航,覆盖整个试验田。光谱测量仪和立体红、绿、蓝(RGB)相机提供了有关与三维(3D)植物结构相关的性状的额外信息,例如叶片角度和目标叶片上的光照强度。由此生成的高时空分辨率的光合效率图谱,为农艺田间试验或植物育种苗圃中植物生长表现提供了详尽的信息。
光合作用不仅为植物生物量的产生和自身代谢的维持提供能量基础,也为豆科植物的固氮作用提供能量1 以及其他共生过程2任何作物的产量潜力均取决于冠层截获的光合光子通量密度(PPFR)的比例εi),这部分辐射中实际用于光合作用并转化为生物量的比例(εc),以及分配到收获产品中的生物量能量比例(εp)3。该 εp 可被设定为收获指数(HI),如果后者被定义为收获产物的生物量除以地上和地下部分的总生物量3以往的育种工作主要集中在提高 εi 和 εc 并且在这方面取得了巨大成功3例如,现代谷类作物和豆科作物基因型可达到 εi 约为90%和收获指数(HI)约为60%3,4,进一步提升的空间十分有限 εi 和 εp3,5相比之下,……的数值 εc 在C中观察到3 和 C4 作物产量仍很少超过三分之一 分别为理论最大值9.4%和12.3%。这表明若能进一步提升作物光合效率,作物生产力仍有巨大的提升潜力 εc 可通过基于自然遗传变异的选择和/或通过基因工程对光合作用装置的组分(如酶1,5-二磷酸核酮糖羧化酶/加氧酶(RuBisCO))进行靶向优化来提高3,5,6,7然而,已有研究表明 i光合效率在田间条件下具有高度动态性,因为植物会不断适应变化的光照及其他环境协变量8,9 ii在稳态条件下测得的光合性状(例如,在人工光照的室内生长室中)10在田间观察到的非稳态光合作用相比,表现出不同的遗传力模式11。这限制了 室内光合作用表型分析在筛选主要粮食作物高产基因型方面的实用性有限,因为这些作物在实际生产中通常是在室外田间种植的。因此,除非能够获得大量作物基因型的高空间和时间分辨率数据,否则难以对光合性状进行有针对性的选择。6,12.
以往在田间进行光合作用无损测量的尝试,主要依赖于手持式气体交换13,14,15 或脉冲振幅调制(PAM)型叶绿素荧光(ChlF)16,17测量设备。这些方法存在一个共同的缺点,即需要手动将待测叶片夹入或放置到测量室中18,19。此外,田间气体交换测量需要数分钟时间以确保测量室内的平衡状态。这使得测量过程费力且缓慢,难以实现足够的通量来广泛研究波动条件下的光合作用调控,尤其难以在育种项目中对数百甚至数千个基因型进行有意义的筛选19。与前述PAM叶绿素荧光测量设备17类似,光诱导荧光瞬变(LIFT)方法20 利用了以下两个事实:i)到达光系统II(PS II)中叶绿素的光能,可用于光合作用,也可通过非光化学淬灭(NPQ)以热的形式耗散,或以叶绿素荧光(ChlF)形式释放;ii)通过强光脉冲饱和使光合作用途径受阻,导致PS II下游电子受体被还原,从而引起叶绿素荧光相应增加。基于这种诱导产生的可变叶绿素荧光,我们可以在环境光照条件下计算PS II的光合量子效率(Fq'/Fm'),或在叶片测量前保持暗适应状态下测定PS II的最大量子效率(Fv/Fm)。
该参数 Fv/Fm 由Kitajima和Butler于20世纪80年代初建立21, Butler22,以及 Björkman 和 Demmig23,报道了最优值为 ≈ 0.83,在不同物种的非胁迫叶片中均如此。该参数的开发与应用为…… Fq'/Fm' (ΦPS II)作为光合作用性能诊断方法始于Genty等人24,表明其与CO2 非胁迫条件下的同化作用。Maxwell 和 Johnson25 和 Baker26 随后提供了一个应用于实践的框架 Fq'/Fm', Fv/F米 以及非光化学淬灭(NPQ),强调它们对环境胁迫的敏感性及其在识别光抑制等损伤症状方面的价值。Murchie 和 Lawson27 强调了潜力(及局限性) Fq'/Fm用于田间表型分析和作物改良。最近,Long 等人8 强调作物的光合作用发生在波动光照条件下,光系统II效率和非光化学猝灭(NPQ)调控的动态变化显著影响碳同化,凸显了结合环境因子的时间变异性来解读荧光参数的必要性。
如先前研究28所述,使用固定式 LIFT 传感器测得的光系统 II 量子效率(Fq'/Fm')与脉冲振幅调制法(PAM)测量结果(R2 = 0.89)以及 CO2 同化速率(R2 = 0.89)具有良好的相关性29,30,31。与此一致,建模分析生长季中Fq'/Fm' 的响应已被证明可用于较好地预测作物光能利用效率 εc,、作物生产力12,32 以及胁迫耐受性33。然而,PAM 类型的叶绿素荧光(ChlF)测量设备采用单次强光脉冲(饱和脉冲)抑制光合作用途径,并辅以弱光脉冲测量最小和最大叶绿素荧光(测量脉冲),整个过程约需 1 秒,因此不适合超高通量测量。相比之下,LIFT 传感器在 750 µs(= 0.00075 s)内发射一系列快速的 300 个高强度闪光(每个闪光强度为 40,000 µmol 光子 m-2 s-1),逐步实现光合作用途径的饱和,同时离散式测量叶绿素荧光产量20。该方法允许测量装置与目标叶片之间保持较大距离19 ,并可通过将传感器牢固安装在自主移动载具上,实现快速、高通量的田间测量32。以往的研究(例如32)在田间条件下的通量往往受限,或仅局限于温室实验。尽管这些方法成功鉴定了参与光合作用调控的候选基因12,但很可能遗漏了在田间条件下发生的重要调控过程。因此,我们假设将 LIFT 传感器安装在自主机器人上,可成为一种合适的工具,用于在田间对大量基因型材料进行Fq'/Fm' 测量(每天可在数千株植物或试验小区上实现每秒多次测量)。
1. 设置PPFR传感器和数据记录仪
2. 设置机器人和LIFT传感器组件
3. 进行LIFT测量
4. 数据整合与预处理
注意:R 代码可在 GitHub 上获取(https://github.com/beat2keller/lift_data_processing)

2025年6月12日至27日,在7个测量日中,共采集了36个大豆(Glycine max (L.) Merr.)育种品系的91,205次叶绿素荧光(ChlF)瞬态测量数据。其中,74,927个数据点可进行地块地理配准,58,916个数据点经过质量筛选后成功与光谱和气象数据匹配。测量设置及代表性数据概览见图1,包括仪器运行状态的照片(图1A)、12个大豆基因型的ChlF诱导动力学曲线(图1B)以及相应的光谱反射率曲线(图1C)。
地理参照数据表明,在两个试验田中,Fq'/Fm' 和 NDVI 均表现出显著的空间变异(图2)。这些空间格局部分与机器人的行进方向相关(图3),因此在后续建模中对这一因素进行了校正,以降低机器人遮荫对目标叶片可能造成的影响。
对光系统II实际光化学效率(Fq'/Fm')对入射光合光子通量密度(PPFR)的响应进行线性混合效应模型分析(图4A),能够解释大部分的变异(R
= 0.45,R
= 0.60)。提取的基因型特异性斜率(响应G:PPFR)在不同育种品系间表现出明显差异(图4B),其中多个品系的响应曲线相较于群体平均值更陡峭或更平缓(图4C)。
最后,基于RGB成像和匹配与立体三维重建算法(Matching And Stereo 3D Reconstruction, MASt3R)36(图 5)获得的冠层水平三维重建结果,展示了将基于LIFT的生理学测量与结构表型分析相结合以实现冠层三维结构表征的潜力。

图 1:大豆基因型光诱导荧光瞬变(LIFT)与光谱反射率测量的概述。(A)LIFT 设备运行时的示例图像。插图显示用于诱导叶绿素荧光(ChlF)的蓝色激发光束。(B)2025 年 6 月 27 日对 12 种大豆基因型测得的瞬态 ChlF 诱导曲线(均值 ± 标准差,每种基因型的 n 为 88 至 553,ntotal = 2,035)。(C)2025 年 6 月 25 日对相同基因型测得的相应叶片反射光谱(均值 ± 标准误,ntotal = 2,035)。颜色代表不同基因型;误差线表示重复测量之间的变异。请点击此处查看该图的放大版本。

图 2:光系统 II 量子效率与归一化植被指数(NDVI)值的空间分布。 在包含 120 个试验小区和 36 个育种品系的两个大豆试验田中,测量了(上图)光系统 II 量子效率(Fq'/Fm')和(下图)归一化植被指数(NDVI)值的空间分布(n = 58,916)。点表示测量位置,颜色表示观测值,形状表示机器人行进方向:西北(NW)、东北(NE)、西南(SW)和东南(SE)。黑色线条标示田间小区边界。黑色圆点表示田间机器人的全球导航卫星系统(GNSS)坐标。请点击此处查看该图的放大版本。

图 3:光系统 II 的量子效率。 箱线图显示了光系统 II 的量子效率(Fq'/Fm')以及在测量时段(9 - 15 时)内 680 nm 波长下入射辐照度的对数值,分别针对东北(NE)方向和西南(SW)方向的探测车行进方向(nNE-ward = 14,339,nSE-ward = 25,669,ntotal = 40,008)。机器人在上午朝 SW 方向行驶时遮挡了测量点,在下午朝 NE 方向行驶时亦如此。每个箱体表示四分位距(IQR),其中横线和须线分别表示中位数和 1.5 × IQR 的界限。请点击此处查看该图的放大版本。

图 4:光合作用响应斜率。(A)基于线性混合效应模型分析得出的两个大豆试验田中光合作用响应斜率(基因型与光合光子通量速率(G: PPFR)之间的相互作用)的空间变异情况。(B)各育种品系经校正的平均值及模型估计的光合作用响应斜率,其中以颜色突出显示感兴趣的基因型。(C)入射 PPFR 与突出显示基因型的光系统 II 的量子效率(Fq'/Fm')之间的关系,附拟合的平方根回归曲线。请点击此处查看该图的高清版本。

图 5:通过 LIFT 测量冠层结构。 显示了大豆品种 Gallec 的冠层结构。(A)光诱导荧光瞬变(LIFT)设备的红、绿、蓝(RGB)相机在筛选光合作用时采集了图像。(B, C)采用匹配与立体三维重建(MASt3R)算法36重建三维冠层结构。金字塔形符号表示估计的相机位置。请点击此处查看该图的放大版本。
方法的通量
LIFT 方法允许传感器与目标叶片之间的测量距离远大于以往的气体交换和PAM型叶绿素荧光(ChlF)测量设备18,19,从而能够在田间和室内环境中实现自动化的高通量测量12,32,37。这消除了 对冠层进行物理接触的需要。该方法的通量随机器人行驶速度的提高而增加。为确保叶绿素荧光(ChlF)产量在时间上的准确测量,在一次 Fq'/Fm' 测量期间所移动的距离必须远小于光照区域的直径(LIFT传感器为20 mm)。当机器人在激发过程中移动时,光照点也随之偏移,导致每次闪光脉冲产生的叶绿素荧光产量降低,最终在过高速度下造成 Fm' 和 Fq'/Fm' 的错误估算。在机器人速度为0.5 m s-1时,最后一次激发闪光所照射的区域中,仅有约2.4%的面积未被首次激发闪光覆盖,这意味着由于机器人运动引起的 Fm' 系统性测量误差在数值上可忽略不计(完整计算见GitHub:https://github.com/beat2keller/lift_data_processing)。以0.5 m s-1的速度,我们在一个40 × 36 m的试验田中每小时可测量300个样地(每个1.5 × 2 m),相当于每天可测量数千个样地。此外,我们的测量装置总重量相对较轻(200 kg),结构紧凑,可通过轻型车辆运输,便于开展多点试验,以评估不同环境下光合效率的表现38。
叶片光合作用与冠层光合作用
该方法中的LIFT测量在空间上仅限于冠层最上层的受光部分,这一部分贡献了总光合作用的约50%至70%39此外,每片叶子能够吸收的光量在很大程度上取决于其相对于太阳的角度,通常更接近垂直方向的叶片会增加叶面积指数,促进光的透射以及冠层光合作用40,41使用本方法中的LIFT传感器组件,可通过光谱仪的反射率测量数据,对3D冠层结构对单叶或整冠光合作用的影响进行建模。28 和/或三维冠层几何结构的显式重建42 基于立体RGB相机系统的图像。新型算法,例如MASt3R36,利用深度变换器学习密集且具有几何感知能力的特征对应关系,从而在宽基线或重复区域实现更鲁棒和精确的三维重建。然而,即使没有此类额外校正,其有效性 已证明,采用自动化高通量方法(如LIFT)可鉴定出更具生产力和抗逆性的作物品种11,43.
库限制
库限制,即由于植物无法利用其能够产生的光合产物量而导致光合作用被主动下调44,45,可能会掩盖不同基因型之间εc 的差异。随着大气中CO2 浓度升高,由库限制引起的光合作用下调可能更加显著45。我们的 方法仅在库限制导致Fq'/Fm'同时下降时才能检测到该现象;已有迹象表明这种情况确实 可能 发生46。无论如何,库限制的相关性似乎强烈依赖于作物种类和发育阶段,小麦在灌浆期受其影响远大于豆科作物45 ,因为后者的共生固氮作用构成了一个强大的额外库1,44,46。如果预计在进行LIFT测量期间某一特定作物的库限制具有高度相关性,则需要进一步研究,将测得的Fq'/Fm'值与同步获取的气体交换数据进行比较,并验证Fq'/Fm'可作为εc良好替代指标这一假设的有效性。
冠层高度的测量距离与变异性
本方法中使用的LIFT传感器组件设计用于对作物冠层在约60 cm距离处进行测量32,47。该方法所用的LIFT传感器对测量距离的微小差异具有较强的抗干扰能力,但通常较大的测量距离会导致Fq'/Fm'的数值略低47。为避免此类偏差,我们建议使用Fq'/Fm'对光合光子通量密度(PPFR)的响应(该指标对测量距离的变化更具鲁棒性),而非直接使用Fq'/Fm'的绝对值,如文献12 所述,和/或在数据分析过程中明确校正冠层高度的差异。
叶片的意外人为遮荫
光照强度增加通常会降低Fq'/Fm'3,8,12。因此,应避免操作人员或测量装置造成不必要的遮光3,8,12。为减少机器人造成的遮光,可设置其行进方向以完全避免遮光32 (但会因无生产性的返程而降低通量),或在交替方向上进行测量,再对遮光效应进行统计学校正。本研究中,Fq'/Fm' 的统计模型中的 Heading × Hour 项(见步骤 4.3.10)确保了机器人对目标叶片的遮光效应已被纳入考虑。
结论
总之,安装在自主机器人上的LIFT传感器能够在田间条件下实现光合效率的快速自动化测量,克服了以往方法在通量方面的限制。尽管冠层结构和库容限制等因素需要仔细考虑,但我们的结果表明,在农艺田间试验和植物育种苗圃中进行可靠且可扩展的光合作用筛选是可行的。高分辨率立体RGB相机系统与作为点传感器的LIFT相结合,将进一步提高测量精度:由人工智能驱动的叶片三维定位将能够选择具有相似太阳朝向的叶片进行采样,并考虑不同基因型特有的叶片倾角。从研究角度来看,比较具有不同冠层结构的基因型的光合作用,以及在开放式空气CO2 加富(FACE)条件下开展高通量LIFT测量,可能尤为有意义,这有助于更好地理解在田间条件下PS II下游过程及潜在的库限制如何影响Fq'/Fm'。在实际应用方面,从育种苗圃到农户田间开展光合性能和抗逆性的筛选,对于精准农业应用具有巨大潜力。
本研究中使用的机器人由Caterra AG公司开发,该公司源自苏黎世联邦理工学院作物科学研究团队,本文大多数作者均隶属于该研究团队。然而,作者们并未直接参与该机器人的开发,且声明不存在其他利益冲突。
感谢Caterra AG提供了机器人底盘、相关技术支持以及协助将其改装用于我们的测量工作。感谢苏黎世联邦理工学院作物科学研究小组的Nicola Storni,他在机器人上安装GNSS定位记录系统方面提供了支持。特别感谢DSP Delley Seeds的Christoph Barendregt和Agroscope的Claude-Alain Bétrix,为PhenoSoy项目试验提供了Agroscope/DSP合作项目的植物材料,该项目由联邦农业局(FOAG)资助。本方法是在“通过改良豌豆提升轮作系统的可持续性与氮素利用效率”(ECOPRot)项目中开发的,该项目由ETH世界粮食系统中心资助,资金来源于Bayer AG的捐赠以及苏黎世联邦理工学院职业种子基金奖。 部分工作使用了大语言模型以辅助编程和语言表达,所有输出内容均由作者核查、验证并最终确定。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 相机三脚架 | - | - | 任何可将PPFR传感器水平安装的类型均可使用。 |
| 汽车蓄电池(12 V / 80 Ah) | Mercedes-Benz Group AG, 斯图加特, 德国 | A 001 982 81 08 | - |
| 数据记录仪 | HOBO Data Loggers, 美国马萨诸塞州伯恩 | HOBO U30 站点, https://www.hobodataloggers.com.au/product/hobo-u30-usb-stand-alone-data-logger/ | - |
| 数据记录仪桌面客户端 | Onset Computer Corporation, 美国马萨诸塞州 | HOBOware 版本 3.7.28, https://www.onsetcomp.com/support/help-center/software/hoboware | - |
| GNSS天线套件 | ArduSimple, 安道尔城, 安道尔 | AS-STARTKIT-LR-L1L2-EUNH-00 | - |
| GNSS记录仪桌面客户端 | u-blox AG, 瑞士塔尔维尔 | u-center GNSS, [原文如此] 评估软件版本 25.03, https://content.u-blox.com/sites/default/files/documents/u-center-25.03_ReleaseNote_UBXDOC-304424225-19688.pdf | - |
| 笔记本电脑 | - | - | 任何可运行与GNSS及数据记录仪桌面客户端兼容的Windows(Microsoft Corporation, 美国雷德蒙德)操作系统的笔记本电脑均可使用。 |
| LIFT传感器组件桌面客户端 | Soliense Inc., 纽约, 美国 | 版本.2016.04 | - |
| LIFT传感器组件 (含240 V市电适配器) | Soliense Inc., 纽约, 美国 | LIFT-REM | 除LIFT传感器本身外,LIFT传感器组件还集成了一个RGB立体相机系统,该系统由FLIR Integrated Imaging Solutions Inc.(加拿大不列颠哥伦比亚省)[现为Teledyne FLIR LLC, 美国威尔逊维尔]提供的两台Blackfly S BFS-U3-50S5C相机组成,以及一台由Ocean Insight [现为Ocean Optics](美国奥兰多)提供的波长范围为400至800 nm、分辨率为0.46 nm的光谱仪。 机器人尺寸约为≈ 170 × 210 × 90 cm(长×宽×高),重量约为≈ 170 kg。LIFT传感器组件的尺寸约为≈ 23 × 34 × 59 cm(长×宽×高,以安装在机器人上的状态为准)。 |
| 逆变电源 (12 V直流转240 V 50 Hz交流) | Green Cell CSG S.A., 克拉科夫, 波兰 | INV08 | - |
| 移动电源 | - | - | 任何USB移动电源均可使用。 |
| PPFR传感器 | HOBO / LI-COR, 林肯, 美国 | S-LIA-M003 | - |
| 机器人控制网站 | Caterra AG, 瑞士奥普菲孔 | - | 链接由机器人制造商为每位客户单独提供。 |
| 机器人原型底盘 (含遥控器及 240 V市电电池充电器) | Caterra AG, 瑞士奥普菲孔 | Firefly(原型,未上市), https://caterra.org/en/technologie/ | 任何具有适当有效载荷能力、履带宽度和离地间隙的其他野外机器人亦可使用。 |
| 白色参考板 | - | - | 任何白色或灰色参考板均可使用。 |