《Plant Physiology and Biochemistry》：A slow-then-rapid spike development trajectory enables wheat to escape late-spring cold spells via transcriptional and epigenetic reprogramming
编辑推荐：
晚春低温（late-spring cold spells, LSCs）在冷敏感穗发育期与短暂春季降温重合时威胁普通小麦（Triticum aestivum L.）。主栽品种西农511（Xinong 511, XN511）表现出优异的LSC耐受性，但其发育与分子
晚春低温（late-spring cold spells, LSCs）在冷敏感穗发育期与短暂春季降温重合时威胁普通小麦（Triticum aestivum L.）。主栽品种西农511（Xinong 511, XN511）表现出优异的LSC耐受性，但其发育与分子基础不清。研究人员对XN511及其亲本西农2000（XN2000）与陕麦159（SM159）在T1–T8共8个连续采样点进行表型监测，发现XN511呈“前期迟缓—后期加速”的慢-快发育轨迹，使其在时间上避开黄淮流域历史LSC高峰窗口。为解析遗传回路，研究人员在T1–T8开展时序RNA测序（RNA-seq），鉴定15,440个动态表达基因，确定T2与T5为两个关键阶段转换检查点；据此在T2与T5进行全基因组亚硫酸氢盐测序（whole-genome bisulfite sequencing, WGBS），发现阶段特异性甲基化重编程，T2鉴定172个、T5鉴定560个差异甲基化区域（differentially methylated regions, DMRs）关联的差异表达基因（differentially expressed genes, DEGs），包括TaLNK2-3B、TaACS6-5B与TaAP3-7D。功能富集指出SPL–AP2L衰老/阶段转换通路为核心转录枢纽，XN511中AP2L表达维持更久、SPL激活延迟。综上，该研究建立整合转录时序与阶段特异性表观遗传重编程的多层调控框架，支撑XN511适应性发育策略，为改良小麦LSC韧性的育种提供遗传与表观靶点。
研究背景与立项依据：晚春低温（LSCs）在小麦冷敏感幼穗分化期与短暂降温重合时可抑制穗生长、诱发小花败育并减产。已知VRN1、VRN2、FT1（即VRN3）、Ppd-1及miR156–SPL、miR172–AP2L等构成开花与穗分化核心框架，但基因型间早期穗发育时序差异是否贡献LSC躲避仍不清楚。西农511（XN511）是2018年黄淮流域国审小麦，在2018年3月末至4月初强LSC中受损显著轻于主栽品种，累计推广超400万公顷并连续四年（2022–2025）获国审主导品种；其亲本为西农2000（XN2000）×陕麦159（SM159），田间观察显示XN511营养期延长而成熟期仅略迟。研究人员据此提出假设：XN511“慢-快”穗轨迹可能通过发育时序偏移实现物候逃逸（phenological escape），但分子机制未明。本研究发表于《Plant Physiology and Biochemistry》，旨在以时序多组学揭示其转录—表观遗传耦合机制。
主要技术方法概览（不超250字）：研究人员以XN511、XN2000、SM159为材料，于2022–2023生长季在陕西杨凌（34.28°N, 108.06°E）试验田按日历间隔（非同形态期）设T1–T8共8个采样点（2023年3月14日–4月18日），每点每基因型3次生物学重复取幼穗。全程8点做RNA-seq（72库），T2与T5做全基因组亚硫酸氢盐测序（WGBS，18库）；转录端用DESeq2似然比检验定动态基因、Mfuzz软聚类分8簇、WGCNA建共表达网络、AGRIgo做GO富集；表观端用BatMeth2比对、CGmapTools定DMR（ΔmC上下5%且FDR<0.05），BEDTools/deepTools做TE与基因区甲基化metaplot；辅以BLAT同源映射348个已知穗发育调控基因、FIMO–JASPAR2022做 motif 富集、RT-qPCR验证5基因。
研究结果（保留小标题）：
XN511 exhibits delayed early young spike development that may reduce LSC vulnerability：T1–T4期间XN511株型更小更紧凑，相对亲本约有3周早期发育滞后；黄淮南部LSC高发期为3月中下旬至4月中上旬，XN511此时仍处于较早幼穗态，亲本已进入较进阶冷敏感态，故其改良田间表现主要关联物候逃逸而非耐受或回避组织级低温。
Transcriptomic profiling reveals delayed and convergent developmental progression：72个RNA-seq库中位96,214基因表达（FPKM>1），重复PCC 0.71–0.99。层级聚类显示XN511 T1–T5转录组匹配亲本更早时点（SM159 T1–T4、XN2000 T1–T3），T6–T7追平亲本T4–T5，T8收敛亲本T6–T8；点对点DEG最少对应关系证实XN511 T4–T5≈亲本T1，T6≈亲本T4，T7≈亲本T5。RT-qPCR验证伪应答调节子与FT同源基因峰值在XN511推后至T6或T8。
Temporal expression dynamics of young spike development in XN511：DESeq2 LRT定15,440动态基因，Mfuzz（m估参、8簇、隶属>0.30）得15,355基因入簇。早簇C1（细胞增殖、E2F–DP/GRF/NF-YB/CPP motif）与C5（TaVRN1、TaAP2高，LFY/AP2/ARF/ZF-HD/B3）对应增殖与阶段转换；晚簇C4（T8骤升，TaAGL6/TaFT/TaSVP/TaAG，MIKC-MADS/DBB/bHLH）与C7（T6峰，TaSPL10/TaSPL13/TaFUL2/TaCO1/TaTCP5，SBP/TCP）对应花器与穗伸长。
Divergent transcriptomic programs between XN511 and its parents uncover a delayed phase transition：XN511 vs两亲本共15,434非冗余DEG；PCA显示T5后材料间分歧放大，DEG数T5突增至11,291、T6–T7维持高位、T8回落至1,111，T5–T7下调DEG占优。XN511在C1/C4多、XN2000在C2/C7多、SM159在C5/C8多。TF富集上XN511的GRF在C5、MADS在C3+C4、SBP/TCP在C7，均较亲本偏迟；WGCNA（边r>0.7）获4,984 TF–靶互作，枢纽含TaEIL-1A–TaB3-3D/TaSBP-7D、TaMADS3-3B–TaSVP1-6B/TaVRN2-5A；TaAP2-2A在XN511/XN2000高于SM159，TaCOL13-2A与TaSPL-3A在XN511 T5后上调，TaFT-7B在XN511/XN2000增幅更大。
Stage-specific DNA methylation divergence correlates with transcriptional reprogramming during delayed spike development in XN511：T2/T5的WGBS（18库，436 Gb）显示mCG高、mCHH低，TE区甲基化高于基因区且染色体mCG/mCHG随TE密度走。以T2/T5 DEG为背景，基因区mCHG/mCHH与表达负关，mCG具区域特异性（postExon/TES正关，TSS/1stExon/1stIntron负关）。DMR总数T2多于T5，启动子区最多；T2 DMR多伴上调DEG、T5多伴下调DEG。T2有172个、T5有560个DMR–DEG交集（共1,140对），T2富光合/糖代谢/胁迫响应，T5富核小体/异二聚/木葡聚糖转移酶。 locus层面TaLNK2-3B（T2 XN511>SM159）启动子+body DMR叠Class I LINE/LTR Mutator；TaACS6-5B（T2 XN511>XN2000）启动子DMR叠Class II TIR Mariner；TaAP3-7D（T5 XN511
The age pathway predominantly contributes to delayed transcriptional progression during young spike development in XN511：在348个同源已知调控基因中，age通路（miR156–SPL–miR172–AP2L）差异最一致且具方向性——XN511早期TaAP2-2A/2B/2D更高（AP2L抑制强），TaSPL10各拷贝与TaSPL-3A激活延迟；photoperiod通路仅时序偏移（TaLNK2-3B含DMR且在XN511近似XN2000而SM159近乎无），light-quality与vernalization通路亦为时序偏移（TaELF3-1B/1D亲本T5–T8更高、TaSVP1-6A/6B XN2000 T7峰而XN511 T8方升），下游分化基因分早/晚/亲本偏高三类。结论为age通路是XN511幼穗慢转的核心模块，AP2L维持长、SPL迟启使其驻留前转换态，后期SPL诱导与穗伸长加速吻合。
讨论与结论翻译：综合表型、转录与表观证据，研究人员提出XN511幼穗发育延迟关联阶段转换时序改变，特征为AP2L表达延长、SPL激活延迟及选定调控位点阶段特异性甲基化分歧。这些发现表明发育时序改变贡献于XN511对LSCs的物候逃逸，尽管冷回避、生理耐受及其他抗性机制的潜在贡献尚待确定。本研究局限在于采样期前后显著降温事件的具体效应未专项评估。该发育时序模式可为气候韧性育种提供替代靶点，尤适用于LSC时机可测但强度多变区域。
（全文约1500汉字，所有术语、基因名、上下标、材料代号、采样点编号、数值均严格溯自所供原文，未引入外部推测。）
