用新型连续定量表型分析方法揭示蓝莓重复开花的遗传控制

《Scientia Horticulturae》：Unveiling the genetic control of repeat flowering in blueberry with a novel continuous quantitative phenotyping approach

　　

多年生作物的开花是受遗传、环境信号和植物活力影响的性状。在此，研究人员研究了蓝莓（Vaccinium corymbosum）中重复开花（repeat flowering, RF）的遗传和环境控制。RF是在一个来自重复开花与非重复开花品种（分别为‘Hortblue Petite’和‘Nui’）杂交的全同胞群体中测量的。利用纵向表型对第一次和第二次开花高峰的曲线下面积（area under the curve, AUC）进行建模。研究人员发现RF受灌木大小和活力的强烈影响，随后将其纳入开花曲线下面积模型。在4号和10号染色体的三个热点区域检测到与首次开花和RF相关的数量性状位点（quantitative trait loci, QTLs），并讨论了已知在温度和光周期控制下调控开花的相关基因。本文报告的表型分析方案和统计建模方法是研究多个遗传位点与环境变量对发育性状（如开花）复杂相互作用的有效策略。现在需要重复的多环境和多年测量实验设计来证实研究人员的结果，并进一步阐明蓝莓中RF的决定因素。

**研究背景**
蓝莓（Vaccinium corymbosum）是多年生落叶果树，其开花受遗传、环境信号（温度、光周期）和植物活力共同调控。重复开花（repeat flowering, RF）现象，即同一生长季内多次开花，在某些品种（如‘Hortblue Petite’）中出现，但遗传机制尚不明确。目前存在的问题包括：RF与植物大小和活力的关系未量化，传统表型方法难以捕捉RF的连续动态变化，且遗传位点与环境变量的交互作用复杂。因此，本研究旨在开发一种连续定量表型分析方法和统计模型，以揭示RF的遗传控制，并评估灌木状态对RF的影响。论文发表在《Scientia Horticulturae》。

**研究内容与结论**
研究人员利用重复开花品种‘Hortblue Petite’与非重复开花品种‘Nui’杂交产生的F1群体（共278个后代，来自正反交），在Motueka果园进行两年观察。通过纵向连续定量表型分析（每周计数开放花朵）和统计建模（计算开花曲线下面积AUC，并校正灌木大小与活力），定位了RF相关QTLs，并鉴定了候选基因。结论：RF受遗传和环境共同调控，灌木大小和活力对RF有显著影响；在4号和10号染色体上发现三个QTL热点区域，包含多个与开花调控相关的候选基因（如DEAH12、TPP、HD3A等）。该研究为复杂发育性状的分析提供了有效策略，对蓝莓育种和应对气候变化具有潜在价值。

**关键技术方法**
本研究采用纵向连续定量表型分析方案，在第二生长季对F1群体（278个后代，来自‘Hortblue Petite’בNui’杂交）每周进行开放花朵计数，并视觉评估灌木大小与活力（1-5级）。利用贝叶斯广义线性混合模型（负二项分布）对开花曲线进行拟合，计算AUC。通过polyqtlR软件进行QTL分析，基于18,195个SNP标记的连锁图谱（总长1,756.9 cM，覆盖12个连锁群）。在参考基因组V. caesariense W85-20 v2.0中搜索候选基因，并使用Tripal MegaSearch和Mercator4进行功能注释。

**研究结果**
**2.1 重复开花及相关性状的表型分析**
通过两年表型数据发现，RF存在与否的比例约为1:1，且两年间高度一致（χ2检验显著）。第一年开花与RF花数存在正相关（r2=0.50）。灌木大小和活力与第一年开花及RF花数显著相关（如灌木大小与第一年开花花数r2=0.52，与RF花数r2=0.47），表明植物状态影响RF表达。

**2.2 开花曲线下面积**
利用AUC量化RF强度。模型揭示了不同植株开花曲线形状和AUC值的差异：RF的起始时间和持续时间变异远大于第一年开花。时间归一化使曲线更平滑，AUC计算基于拟合曲线在原始时间尺度上的积分。

**2.3 灌木大小和活力对RF曲线的影响**
灌木大小和活力强烈影响RF的发生和强度。大多数高活力（评分4-5）大灌木表现出RF，而小弱灌木（评分1-2）则无RF。预测曲线显示RF初始快速上升后逐渐下降，且大小和活力类别间曲线形状相似。

**2.4 QTL定位结果**
在连锁群（LG）4、10和11上检测到与第一年开花和RF相关的QTL。在LG4上存在两个热点区域：一个对应第一年开花（25.1-41.5 Mbp），另一个包含RF和灌木大小/活力QTL（45.2-52.6 Mbp）。在LG10上存在一个RF热点（8.2-17.6 Mbp），涵盖所有RF性状的QTL。LG11上仅检测到一个第一年开花结束日期的QTL。此外，在LG6上检测到灌木大小QTL。

**2.5 灌木与RF QTL的遗传相关性**
LG4上的灌木QTL与RF QTL存在共定位（峰位置相距数cM），但未发现多效性。LG10上的RF QTL与灌木QTL无相关性。然而，LG4灌木活力QTL的等位基因影响大灌木中所有RF性状的分布，表明植物状态与RF之间存在遗传互作。

**2.6 候选基因**
在染色体4热点区域鉴定到DEAH12（RNA解旋酶，与开花相关）、OsZHD1（锌指同源域转录因子，调控光周期开花）、TPP（海藻糖-6-磷酸磷酸酶，参与T6P信号调控开花）和GGPPS（牻牛儿基牻牛儿基焦磷酸合酶，影响侧枝生长）等候选基因。在染色体10热点区域鉴定到HD3A（FT同源，成花素）、AFP3（长日光周期开花关键调控因子）、FLD（自主开花途径基因）、NAC52（转录抑制因子）和CLE相关蛋白等。GO富集分析显示，转移酶、氧化还原酶、RNA生物合成（DNA结合转录调控）和泛素-蛋白酶体系统等类别在两条染色体热点区域均显著富集。

**讨论与结论**
讨论部分总结：本研究通过连续定量表型分析和AUC建模方法，揭示了RF的遗传控制，并证明灌木状态是重要协变量。QTL热点中的候选基因提供了调控线索，特别是独立于光周期和低温的RF调控机制，对应对气候变化和延长产季具有潜力。但需更多重复实验（如多环境试验、定量测量灌木体积）和更广泛遗传背景验证QTLs。结论部分翻译如下：通过本研究应用的连续定量表型分析方案和开发的AUC建模方法，研究人员推进了当前对蓝莓RF遗传控制的认识。由于该性状的高度复杂性，包括其遗传决定因素以及与环境变量和植物生长状态的相互依赖性，需要额外的实验来进一步阐明该性状。特别是，需要跨多种环境和灌木大小与活力类别的重复研究来证实本研究的发现。例如，在受控条件下进行多环境试验，独立调控昼夜节律和温度，或在光周期和温度不同的地点进行试验，将有助于更好地理解它们对RF的单独和交互效应。此外，使用定量测量如植物体积和春梢、夏梢数量将比本研究所用的视觉评分系统更准确地评估灌木大小和活力，从而允许更精确的G×E模型。此外，在将已鉴定的QTL应用于育种计划之前，需要在更广泛的遗传背景（即来自不同亲本和/或种质群体的家系）中进行验证。尽管如此，本研究鉴定的候选基因集为RF的可能调控机制及其独立于低温和光周期的发生提供了见解，这突显了其在应对气候变暖和延长反季节生产方面的巨大潜力。通过精细定位缩小QTL区间以及通过RNA测序分析转录表达数据，未来有助于优先考虑候选基因，然后可通过qPCR或转基因等方法进行功能验证实验。