苹果红肉着色的气候驱动因素

《Journal of Agriculture and Food Research》：Climatic drivers of red-flesh coloration in apple

　　

近年来，对新苹果品种尤其是红肉品种的需求激增。果实红肉着色源于花青素（anthocyanin）积累，受遗传、发育和环境因素调控。尽管光照和温度等气候单因子效应已有研究，但其在果园条件下对采收期着色的组合与动态影响仍不清楚。为探究基因型（Genotype）×环境（Environment）互作，研究人员基于6个红肉基因型在5个气候差异站点上的非靶向多站点方法，结合图像分析与气候数据建立统计模型，以确定气候因子对红肉色强度在不同果实发育阶段的隐含作用、排序及组合效应。结果凸显了最高气温（Tmax）、昼夜温差（Tmax?Tmin）和相对空气湿度（Hum）在红肉着色中的关键作用，刻画了它们的组合效应并按发育阶段识别了最优气候范围。最后，研究人员通过等位基因依赖的气候响应揭示了显著的基因型×环境互作。

研究背景方面，苹果（Malus domestica Borkh.）全球产量巨大，红皮与红肉表型因关联成熟度、风味及抗氧化健康效益而备受育种关注。红肉性状由野生种Malus sieversii f. niedzwetzkyana引入，核心由MdMYB10基因R6等位（R6-Mdmyb10）启动子六串联重复驱动组成型花青素合成，但表型强度在基因型间差异显著，还受LG16位点（F/a/A等位）修饰。已有研究多聚焦光（辐射）或平均温等单因子恒定处理，对多因子组合、动态变化及按果实发育阶段（细胞分裂、细胞膨大、成熟）解析不足；在气候变化导致热浪与干旱频发的背景下，田间生产条件评估尤为迫切。为此，研究人员在法国IFO牵头IFORED联盟下，于法、意、西、美、澳5个气候迥异站点种植6个R1/R6背景红肉基因型（RF01–RF06，LG16分别为Aa/Aa/Aa/aF/AA/Aa），2023–2024两季采收赤道横切面扫描提取CIELab*中a_mean作花青素proxy，匹配开花至采收日气象站数据，旨在量化田间红肉强度变异、判定各发育阶段有利/不利气候区间并构建预测模型。

主要关键技术方法上，研究人员构建基因型×站点×年份混合效应模型（lme4的lmer）估计方差组分与遗传力H²；以全期气候聚合变量（均值、总和、方差）为预测元，a_mean为响应，分别用Random Forest（MDA重要性）与OPLS（VIP>0.8筛选）锁定核心气候驱动；针对RF01与RF06采用SpiceFP 2D_1/2扩展法，将两连续气候预测元与三分类发育阶段（Phase 1: 0–30 DAFB；Phase 2: 31–130 DAFB；Phase 3: 131 DAFB–采收）离散网格化，联合Lasso与Fused Lasso惩罚识别组合效应区，并导出LASSO约束农艺模型。样本队列为5站×6基因型×2年共60个基因型-环境组合，每组合2株树各8果。

研究结果部分，3.1节五站点气候对比显示西班牙最暖干（T_mean≈21°C，Hum≈60%，Rf<1mm），法国最凉湿（T_mean13.8–17.9°C，Hum 72–82%），PCA（PC1+PC4）清晰分离站点与年份，证实温辐射与湿风呈拮抗轴。3.2节六基因型等位构成表明RF05（AA）a_mean最高（41.01），RF04（aF）最低（18.08），验证LG16 A等位增强、F等位削弱着色。3.3节混合模型方差组分中站点方差64.23、基因型54.03、G×E 15.61，H²=0.27，即73%变异归非遗传（环境+G×E+残差），美国站点偏离正、西班牙负，证明环境限制显著。3.4节RF与OPLS一致选出Dif_period_T_max_T_min、Mean_T_max、Mean_Hum、Mean_Etp、Mean_Rad为核心；OPLS R²Y(cum)=0.498、Q²(cum)=0.49，气候解释约50%肉质红度。SpiceFP 2D_1/2在RF01/RF06上揭示T_max与Hum组合随阶段移动：细胞分裂期偏好中温中湿，膨大期需较大昼夜温差伴适中湿度，成熟期低温高湿区间正效应最强，高温（>30°C）显著负向。

讨论与结论翻译浓缩如下：研究人员指出，红肉着色并非单因子线性响应，而是T_max、昼夜温差、相对湿度在多发育阶段门限组合的结果；Mdmyb10 R6提供组成型基础，LG16等位设定上限，气候在田间将潜在遗传力压缩或释放。研究首次在跨半球多站点尺度给出按Phase划分的气候最优窗口，可用于红肉苹果选址（如规避西班牙暖干站点）、灌溉遮阴管理（膨大期控温增湿差）及气候适应育种（选A等位配中高昼夜温差区）。结论确认：最高气温、昼夜温差与相对湿度为关键驱动；其组合效应具阶段特异性；基因型通过对气候响应的等位差异实现显著G×E互作，该框架为果园管理与育种提供量化依据。