硬粒小麦类胡萝卜素色素沉积的物候降水窗口：网络内预测与跨环境基因型稳定性

《European Journal of Agronomy》：Phenological precipitation windows are associated with carotenoid pigmentation in durum wheat: Within-network prediction and genotype stability across contrasting environments

　　

后期降雨与硬粒小麦黄色指数降低相关。
随着作物发育，降雨与黄色指数的关联方向发生反转。
网络内降雨模型：留一站验证R²为0.73（置信区间CI 0.41–0.84）。
生长度日（GDD）可减少由日历驱动产生的虚假温度趋势。
‘Nachit’（斜率为0.74）是色素稳定性育种的候选亲本。

研究背景与问题：硬粒小麦（Triticum turgidum L. subsp. durum）是地中海盆地及以外地区意大利面、蒸粗麦粉和粗粒粉制品的原料；粗粒粉黄色是市场价值的关键品质性状，主要由胚乳类胡萝卜素（carotenoids）尤其叶黄素及其酯决定。黄色指数（yellow index, YI）受遗传与环境共同控制，单点试验遗传力中至高，但多环境分析中明显下降，基因型×环境互作（genotype × environment interaction, G×E）是重要变异来源。多环境网络环境组分常主导色素变异，但具体气候相关因子及其敏感发育阶段尚未解析；既有证据多基于最终籽粒组成，未区分哪一发育阶段对哪一气候变量敏感。日历月份对齐物候会引入噪声，而生长度日（growing degree days, GDD）对齐在硬粒小麦气候–颜色分析中应用较少；蒸汽压亏缺（vapour pressure deficit, VPD）作为大气蒸发需求指标受到关注，但与温度高度共线。摩洛哥国家硬粒小麦多环境试验网络跨4个生态区，为阶段解析和网络内预测提供较广推断空间。气候变化下降雨变率和后期热事件增加，使品质性状的气候关联解析具有必要性。研究人员据此识别关键发育阶段气候变量与YI关联，检验GDD对齐是否揭示阶段依赖效应，建立简约预测模型并进行网络内交叉验证，同时刻画降水–品质风险梯度和基因型敏感性；设计为观察性，报告环境关联而非因果效应。

研究开展、总体结论与意义：研究人员基于2017–2023年7个收获年、7个试验站和35个品种的网络数据，发现成熟期累积降水是YI最强负相关因子，营养期降水呈正相关，降雨–YI关联随发育发生方向反转；简约三特征降水模型在网络内留一站交叉验证中R²为0.73，但真正新地点预测尚未建立；‘Nachit’在采样充分的品种中斜率最低，是色素稳定性育种候选亲本；GDD物候对齐优于日历对齐；网络平均YI在七年内下降，但序列过短。该研究为网络内品质风险预警、品种推荐和色素稳定性育种提供了模板，并强调外部验证的重要性。

关键技术方法（不超过250字）：研究样本来源于摩洛哥国家硬粒小麦多环境试验网络（7个试验站、4个生态区、2016/17–2022/23共7个收获年、35个品种；28个站–年环境进入气候分析）。YI以全谷粉CIE b^*测定，2019–2020年由总类胡萝卜素浓度经校准转换。气候采用ERA5-Land与CHIRPS日数据，按从播种起GDD（T_base=0 °C）划分7个物候期。统计包括49个气候特征与YI的Pearson相关及Benjamini–Hochberg错误发现率（false discovery rate, FDR）控制、站去均值、偏相关、18天滑动窗口与4000次置换族系阈值、前向逐步选择、留一站交叉验证（leave-one-station-out, LOSO）/留一年交叉验证（leave-one-year-out, LOYO）及嵌套交叉验证、站分层bootstrap、组内相关系数（intraclass correlation coefficient, ICC）方差分解、Finlay–Wilkinson回归、分段回归与风险分箱，并开展探索性VPD与土壤分析。

研究结果：
3.1 黄色指数变异与环境主导：通过35品种×29环境矩阵的ICC方差分解，得出约66%的YI变异归因于站–年环境；收获年独特贡献42%，站独特贡献11%，站×年贡献8%；2018年网络均值最高（18.4 b^*），2023年最低（14.1 b^*），表明环境是主要变异来源。
3.2 物候期特异性气候–YI关联：通过对49个气候特征×物候期相关进行FDR控制，16个在q<0.05显著、3个在q<0.01显著；成熟期累积降水（ΣP）负相关最强（r=-0.706），营养期降水呈正相关；偏相关控制成熟期ΣP后该信号仍独立（偏r=-0.682），灌浆期最高温度（T_max）和VPD的关联明显衰减，表明存在水分时序不对称。
3.3 关键发育窗口：通过18天滑动窗口逐日相关和置换阈值分析，降水在播后130天（灌浆期）出现r=+0.71峰值，且对12–21天窗宽稳定；成熟期转为负相关，零交叉日148接近灌浆–成熟边界；热和VPD变量无单一关键期，说明降水敏感性具有明确阶段定位。
3.4 成熟期降水作为主导预测因子：通过线性回归得出成熟期累积降水与YI强负相关（r=-0.706，R²=0.50），每增加1 mm对应Δb^*=-0.101；原生b^*子集更强（r=-0.820，Δb^*=-0.134）；分段回归无显著断点；风险分箱中0–10 mm有10%环境低于b^*<16，10–20 mm为67%，>20 mm为75%，但区间宽，属描述性梯度。
3.5 VPD作为候选机制通路：通过相关和偏相关分析，灌浆期VPD与YI负相关（r=-0.509），但其与灌浆期T_max高度共线（r=0.947），控制其一后另一变量关联弱，表明VPD独立贡献无法与热胁迫分离。
3.6 预测模型与网络内交叉验证：通过前向逐步选择和交叉验证，三特征降水模型（成熟期ΣP、灌浆期最大日降水（P_max）、拔节/孕穗期ΣP）获得LOSO R²=0.73（95% CI 0.41–0.84；嵌套交叉验证0.66），原始尺度LOYO R²=0.56；全49特征岭回归LOSO R²<0，提示高维过拟合，结论为网络内可预测但外部未验证。
3.7 基因型敏感性：通过Finlay–Wilkinson回归，32个注册品种斜率差异显著（F(31, 641)=1.80，p=0.0055），斜率范围0.51–1.40；‘Nachit’斜率0.74、为采样充分品种中最低，是候选亲本但未通过32次比较校正；国产与进口欧洲品种平均斜率无差异证据（p=0.69），说明色素稳定性应直接选择。
3.8 气候趋势与品质展望：通过七年趋势和GDD/日历对比，网络平均YI下降0.56 b^* yr^?1（p=0.009）；成熟期降水不显著上升（+2.57 mm yr^?1）；日历法灌浆期T_max呈假升温（+0.22 °C yr^?1），GDD对齐后为+0.10 °C yr^?1且不显著；年内混合模型仍显示成熟期降水负关联，提示趋势解释需谨慎。
3.9 探索性土壤分析：通过SoilGrids站点黏粒含量与降水敏感性关联分析，未发现显著相关（ρ=0.31，p=0.54），仅为假设生成。
3.10 敏感性分析：通过原生b^*、剔除灌溉站等子集重复分析，核心成熟期降水关联和模型性能稳定（原生b^* r=-0.820；LOSO R²=0.71–0.74），时间交叉验证更可变，说明核心结论不依赖单一站点或测量方法。

讨论部分总结：讨论认为成熟期累积降水是最强环境相关因子，稳健性检验和偏相关支持其独立性；水分时序不对称是核心发现：营养期降雨与较高色素沉积相关，成熟期降雨与较低色素沉积相关。候选机制包括成熟期水分促进酶促、氧化或微生物降解，灌浆期VPD与热胁迫共线无法分离，病害、脂氧合酶及再湿等通路未被直接测量。GDD对齐明显优于日历对齐，日历模型空间性能差且产生虚假温度趋势。三特征降水模型属于网络内信号而非已验证的外部预测器，全特征岭回归显示过拟合。风险梯度可构成潜在决策支持框架，但区间宽且需外部验证及操作条件。育种上品种敏感性差异显著，‘Nachit’为候选亲本，国产与进口品种平均敏感性无差异，色素稳定性应作为直接选择目标。七年气候趋势过短，2023年网络均值低于商业阈值，成熟期降水增加不显著，长期趋势不能确立；局限包括品种特异抽穗期未纳入、站点数有限、成熟期边界未田间验证、病害/氮/辐射/湿度未测、全谷b^*与粗粒粉色并非同一性状。

研究结论翻译（依据原文Conclusions）：成熟期累积降水是硬粒小麦类胡萝卜素色素沉积的主导环境相关因子：在28个站–年环境中，后期降雨每增加1 mm，YI降低0.101个b^*单位（r=-0.706，p<0.0001）；在23个以b^*原生测量的环境中降低0.134个b^*单位（r=-0.820），解释环境水平方差的一半，并在控制灌浆期温度和VPD后仍存在（偏r=-0.682）。全发育期存在水分时序不对称：营养期降雨与较高色素沉积相关，成熟期降雨与较低色素沉积相关。简约三特征降水模型在网络内交叉验证中重现该信号（LOSO R²=0.73，95% CI 0.41–0.84；每个折内重新选择特征时为0.66；敏感性子集为0.71–0.74），在原始尺度且站截距仅由训练年份估计时LOYO R²=0.56；真正新地点预测仍有待建立。GDD对齐至关重要：日历模型空间失败（LOSO R²=0.21），并产生虚假灌浆期温度趋势（+0.22 °C yr^?1），GDD对齐后消失（+0.10 °C yr^?1，p=0.87），提示多地点品质研究不应使用日历月份替代物。32个注册品种间基因型敏感性显著不同（F(31, 641)=1.80，p=0.0055），范围为0.51–1.40；‘Nachit’在采样充分品种中斜率最低（0.74，95% CI 0.50–0.98），是培育色素稳定性的候选亲本，但无单个品种在32次比较校正后仍显著。国产与进口欧洲释放品种平均敏感性无差异证据（p=0.69），色素稳定性应直接选择而非假定随本地适应而来；整合分析与经验G×E方法也可能有益于品质性状选择。网络平均YI在七个收获年下降（?0.56 b^* yr^?1，p=0.009），同时成熟期降水不显著增加（+2.57 mm yr^?1，p=0.21）；七年序列太短，不能确立长期趋势。VPD–类胡萝卜素降解通路机制上合理（r=-0.509，p=0.006），但统计上与热胁迫不可分离（VPD–T_max r=0.947站去均值，0.927原始），需实验确认。三种潜在输出为：基于降水的风险梯度，成熟期降水<10 mm时低于商业阈值环境比例为10%，10–20 mm为67%，>20 mm为75%，区间分别为0–45%、30–93%、19–99%，数据支持的是梯度而非某一截断值；基于CHIRPS卫星降水和GDD追踪的网络内预测模型，可能构成决策支持框架，但当前试验网络外可迁移性尚未独立验证；品种敏感性排序可指导品种推荐。