《Horticulturae》：ABA-Responsive Peach PpMYB6 Enhances Freezing Tolerance and Restricts Plant Growth: PpCBF2 as a Direct Transcriptional Target
编辑推荐：
桃（Prunus persica）的生产及其栽培边界的北扩受到频繁极端气候事件的严重制约。基于两个独立平行处理（冷胁迫和ABA施用）的时间序列转录组学分析及qRT-PCR验证，研究人员鉴定出多个对冷胁迫和ABA共同响应的候选转录因子基因（PpERF48、PpE
桃（Prunus persica）的生产及其栽培边界的北扩受到频繁极端气候事件的严重制约。基于两个独立平行处理（冷胁迫和ABA施用）的时间序列转录组学分析及qRT-PCR验证，研究人员鉴定出多个对冷胁迫和ABA共同响应的候选转录因子基因（PpERF48、PpERF017、PpMYB6、PpWRKY46、PpWRKY40和PpbHLH35）。其中，通过生物信息学分析以及对过表达PpMYB6的转基因桃愈伤组织和拟南芥（Arabidopsis thaliana）株系的进一步表征，揭示了PpMYB6在调节植物生长和赋予低温胁迫耐受性中的双重作用。酵母单杂交和双荧光素酶报告基因实验证实，PpMYB6直接结合并激活PpCBF2启动子。酵母双杂交文库筛选鉴定出DWARF8作为候选互作蛋白，提出了一个工作假说，即PpMYB6可能通过赤霉素途径负调控营养生长。总之，本研究表征了一个与桃耐寒性和生长平衡相关的新分子模块，为培育耐寒品种和砧木提供了有前景的候选基因。
一、研究背景与研究目的
低温胁迫（包括<20℃的冷害胁迫和<0℃的冰冻胁迫）是限制桃产业进一步扩张及其北向栽培边界的主要非生物限制因素。内源植物激素脱落酸（ABA）作为经典的胁迫激素，在介导植物对冷、旱、盐碱胁迫的耐受性中发挥关键作用。外源ABA处理已被证实能有效缓解胁迫诱导的过氧化损伤并维持光合系统稳定性。尽管已有研究表明ABA通过多层机制（包括激素互作、层级转录因子级联和蛋白磷酸化）调控冷胁迫反应，但冷信号与ABA信号通路如何交互调控桃树低温耐受性的分子机制仍不清楚。为此，研究人员通过对桃幼苗叶片进行时间序列RNA测序（RNA-Seq），旨在解析冷胁迫和ABA处理各自独立诱导的全基因组转录程序，鉴定共同响应的核心转录因子，并重点表征了PpMYB6的功能及其下游靶标，为培育耐寒桃品种提供分子候选基因。该论文发表于《Horticulturae》。
二、主要关键技术方法
- 1.
时间序列转录组测序（RNA-Seq）：以耐寒桃品种‘世纪之星’一年生盆栽苗为材料，设置0℃冷胁迫处理和50 mg·L-1 ABA叶面喷施处理两个独立平行实验组，分别在0、1、2、4小时采集叶片进行RNA-Seq，共获得157.01 Gb高质量数据。
- 2.
差异表达基因（DEG）鉴定与功能富集：采用DESeq2鉴定DEGs，并通过Gene Ontology（GO）和KEGG通路富集分析解析其生物学功能。
- 3.
实时定量PCR（qRT-PCR）验证：对6个共响应转录因子进行表达模式验证。
- 4.
生物信息学与启动子分析：利用ExPASy、SWISS-MODEL、PlantCARE等数据库分析PpMYB6蛋白理化性质、保守结构域及启动子顺式作用元件。
- 5.
转基因功能验证：构建PpMYB6过表达载体，分别转化桃愈伤组织和拟南芥，并在-5℃（愈伤组织）和-4℃（拟南芥）条件下评估抗冻性。
- 6.
酵母单杂交（Y1H）与双荧光素酶报告基因实验：验证PpMYB6对PpCBF2启动子的直接结合与转录激活。
- 7.
酵母双杂交（Y2H）文库筛选：以截短的PpMYB6为诱饵筛选桃cDNA文库，鉴定互作蛋白。
三、研究结果
3.1 转录组测序与质量控制：获得高质量测序数据，GC含量超过45%，Q30比例持续高于93.12%，为后续分析奠定了坚实基础。
3.2 样本相关性分析与主成分分析（PCA）：各处理组内生物学重复间决定系数（R2）超过0.95，表明重复性极高。PCA显示ABA处理样品主要分布在PC1负轴，冷胁迫样品则随时间沿PC1正轴偏移，表明两种处理对桃全转录组产生显著且不同的影响。
3.3 DEG鉴定与表达动态：ABA处理4小时后DEG总数最多（2962个），冷胁迫4小时后次之（1898个）。维恩图分析发现1105个DEGs同时对冷胁迫和ABA处理做出响应。层次聚类将DEGs分为四个簇，其中簇b基因在两种处理下均随时间持续上调，代表了保守的胁迫适应模块。
3.4 GO功能富集分析：ABA处理后，早期DEGs富集于防御反应、胁迫应答和转录因子结合等过程，后期转向跨膜运输和细胞稳态维持。冷胁迫下，早期以激素信号和内源刺激应答为主，后期扩展至多糖代谢、转录调控和细胞外基质重塑，体现了从早期感知到后期代谢重塑的转变。
3.5 KEGG通路富集分析：ABA处理下，富集通路从早期的基础代谢合成转向后期的信号转导与胁迫应答通路。冷胁迫则呈现相反的时序变化，即早期集中于植物激素信号转导和MAPK级联，后期转向次级代谢生物合成。两者共同显著富集于苯丙烷生物合成和类黄酮生物合成等次级代谢通路。
3.6 qRT-PCR验证：对六个共响应转录因子（PpERF48、PpERF017、PpMYB6、PpWRKY46、PpWRKY40和PpbHLH35）进行qRT-PCR验证，其表达模式与RNA-Seq数据高度一致，其中PpMYB6在冷胁迫2小时表达量达到峰值（约为对照的4倍）。
3.7 PpMYB6生物信息学分析：PpMYB6蛋白由368个氨基酸组成，分子量42.22 kDa，等电点5.90，位于桃4号染色体。其N端15-118位氨基酸含有典型的R2R3-MYB和SANT保守结构域。启动子区域密集分布ABA响应元件（ABRE）和干旱/冷响应元件（DRE/CRT），证实其受ABA和冷信号的双重调控。
3.8 转基因愈伤组织和拟南芥表型分析：过表达PpMYB6的桃愈伤组织在-5℃处理后鲜重仅下降5.8%（野生型下降51.2%），抗冻性显著增强，但常温下生长速率略慢。在拟南芥中，过表达株系在-4℃处理后存活率（77.6%-81.9%）显著高于野生型（56.1%），且常温下植株高度（约15-16 cm）显著矮于野生型（23.5 cm），证实PpMYB6在增强抗冻性的同时抑制营养生长。
3.9 PpMYB6直接结合并激活PpCBF2启动子：Y1H实验证明PpMYB6可直接结合PpCBF2启动子。双荧光素酶报告基因实验显示，共表达PpMYB6与PpCBF2启动子驱动的LUC报告基因后，荧光强度及LUC/REN比值显著升高，确认PpCBF2是PpMYB6的直接转录靶标。
3.10 Y2H筛选PpMYB6互作蛋白：由于全长PpMYB6具有自激活活性，研究人员采用N端截短片段（1-130 aa）作为诱饵，从桃cDNA文库中筛选出24个非冗余候选互作蛋白。其中，DWARF8（一种DELLA蛋白）是赤霉素（GA）信号通路的关键转录抑制因子，暗示PpMYB6可能通过GA通路负调控生长。
四、讨论与结论总结
讨论部分：本研究发现冷胁迫和ABA处理共享大量重叠的DEGs（1105个），并在苯丙烷生物合成、类黄酮生物合成、MAPK级联和植物激素信号转导等通路中共富集，表明这些通路是果树保守的抗冷机制。值得注意的是，两种刺激引发了不同的早期动力学响应：冷胁迫优先触发信号转导级联，而后进行代谢重塑；而ABA作为化学信使则直接进入细胞启动转录级联，加速下游代谢合成。PpMYB6在桃愈伤组织和拟南芥中均表现出增强抗冻性但抑制生长的双重表型，这与ABA协调胁迫耐受性与生长权衡的已知功能一致。通过Y1H和双荧光素酶实验确认PpMYB6直接激活核心冷信号组分PpCBF2，为抗冻表型提供了分子联系。Y2H筛选到的候选互作蛋白NCED1（ABA合成限速酶）、DWARF8（DELLA蛋白，GA通路抑制因子）和SRC2（转录共激活因子）为进一步解析PpMYB6平衡冷耐受与生长的分子网络奠定了基础。
结论部分：总之，整合的时间序列转录组分析和实验验证确定了PpMYB6是一个被ABA信号和冷胁迫共同激活的核心转录因子基因。通过转基因桃愈伤组织和拟南芥的功能表征表明，PpMYB6增强了抗冻性，同时抑制了营养生长。在机制上，PpMYB6直接结合并反式激活PpCBF2启动子，确立了PpCBF2作为冷信号通路中直接下游靶标的地位。此外，蛋白质互作组筛选鉴定了候选伙伴，特别是DELLA蛋白DWARF8，提出了一个工作模型，即PpMYB6通过赤霉素信号串扰来抑制生长。总体而言，本研究破译了一个平衡环境适应力与生长的新型分子节点，为培育耐寒桃品种和砧木提供了一个卓越的遗传候选基因。
