小麦是全球最重要的粮食作物之一，为了加速优良品种的培育，科学家们发展出了一种名为“双单倍体（Doubled Haploid, DH）”的技术。这项技术能够快速获得基因完全纯合的植株，从而节省大量育种时间和成本。实现DH技术的一条重要途径是“小孢子胚胎发生”，通俗地说，就是让原本注定要发育成花粉（雄配子体）的未成熟细胞，在特定条件下“改道”，转变为一颗可以发芽生长的胚胎。然而，就像并非所有人都会弹钢琴一样，不同的小麦品种对这种“细胞改道”的反应能力差别巨大，大部分普通小麦品种对此表现得十分“迟钝”，这严重限制了该技术在小麦育种中的广泛应用。为什么有些小麦品种能高效地产生胚胎，而另一些却几乎“无动于衷”？这背后隐藏的分子密码是什么？
为了揭开这个谜题，一个研究团队在《Plant Cell Reports》杂志上发表了一项重要研究。他们以两个表现截然不同的春小麦品种——“Nanda”（高响应）和“Sadash”（难响应）为研究对象，系统地比较了它们在小孢子胚胎发生三个关键阶段（第0天单核期小孢子、第7天原胚、第21天胚状结构）的基因表达全景图。研究人员通过整合差异表达基因（DEG）分析、功能富集、同源基因特异性分析和机器学习等方法，不仅描绘了成功胚胎发生的转录组景观，更重要的是，揭示了造成品种间巨大差异的核心分子机制，并找到了潜在的调控“开关”。
为了开展这项研究，研究人员主要应用了以下几项关键技术：首先，他们建立了Nanda和Sadash两个基因型的小孢子培养体系，并在三个关键发育时间点（D0， D7， D21）收集了生物样本。其次，利用RNA-seq（RNA测序）技术对总共18个样本（2品种×3阶段×3生物学重复）进行了转录组测序。然后，使用生物信息学流程（包括使用iDEP， DESeq2等软件）进行差异表达基因（DEG）鉴定、功能富集（GO分析）和主成分分析（PCA）。此外，研究还特别针对小麦的六倍体特性，进行了同源基因（Homoeolog）表达动力学分析，以考察A、B、D三个亚基因组基因的表达贡献变化。最后，应用机器学习算法（MLSeq包， VoomNSC分类器）对转录组数据进行分类和特征选择，以挖掘潜在的生物标志物基因。
研究结果
高响应基因型Nanda的胚胎发生效率
研究首先在表型上确认了Nanda和Sadash的巨大差异。在诱导培养后，Nanda平均每穗能产生1002个胚状结构（Embryo-Like Structures, ELS），而Sadash仅产生25个。在再生植株中，Nanda有51%是绿苗，而Sadash仅有2%。这确立了Nanda作为高响应基因型的研究价值。
胚胎发生过程中的转录组动态与基因型变异
主成分分析（PCA）显示，样本能清晰地按发育阶段和基因型聚类，表明转录组响应具有阶段特异性和基因型特异性。差异表达分析发现，在从D0到D7的胚胎发生启动阶段，Nanda发生了比Sadash更剧烈、更广泛的转录重编程，显示出更强的转录可塑性。
同源基因表达动力学
由于普通小麦是由A、B、D三个祖先基因组融合而成的六倍体，其多数基因以三个一组的“同源基因”形式存在。研究发现，在胚胎发生过程中，大部分同源基因表达是平衡的。然而，一个关键发现是，在难响应的Sadash中，B亚基因组基因被显著抑制，并且这种抑制状态比在Nanda中更稳定。这些被持续抑制的B亚基因组基因中，有超过300个同时出现在品种比较的差异表达基因列表中，它们在Sadash中表达上调。这表明B亚基因组同源基因的表达抑制可能是导致Sadash胚胎发生能力低下的一个重要原因。
驱动胚胎发生进程的分子机制
通过对Nanda胚胎发生过程中上调基因的功能分析，研究人员揭示了成功“细胞改道”需要协同激活的多个生物学过程。在启动阶段（D0到D7），与表观遗传调控相关的核小体组装、染色质重塑和组织过程被显著富集。同时，涉及细胞骨架（如微管组织）、细胞壁动力学、代谢途径（如谷胱甘肽代谢、多糖代谢）以及激素信号（如生长素信号）的基因也大量上调。这表明胚胎发生起始于深刻的表观遗传和细胞结构重组，并伴随着代谢转换和信号激活。在后续发育阶段（D7到D21），代谢和信号通路持续活跃，而染色质相关活动则转向下调。
导致Nanda高胚胎发生效率的机制
比较Nanda和Sadash在同一发育阶段的基因表达差异，研究人员发现了阶段特异性的调控模式。在启动的关键期（D7），Nanda相较于Sadash显著上调了包括核小体组装、染色质重塑在内的表观遗传过程，以及多种代谢和应激反应通路。到了胚胎发育后期（D21），与细胞壁合成（如纤维素代谢）和特定氨基酸代谢相关的基因在Nanda中更为活跃。这揭示了Nanda的高效率源于其在正确的时间点，更强烈地启动了驱动胚胎发生的核心调控和代谢网络。
与胚胎发生相关的转录因子
转录因子是基因表达的“指挥官”。研究鉴定出超过1500个在胚胎发生过程中差异表达的转录因子。其中，bHLH、WRKY、ERF、bZIP、GRAS、ARF等家族成员在早期被激活。研究特别关注了一些已知的胚胎发生“主调控因子”，如WRI/PLT/AIL（属于AP2家族）、WOX、GRF等，它们在Nanda的胚胎发生过程中持续上调，可能对建立和维持胚胎发育程序至关重要。
富含差异表达基因的基因组区域
通过基因组扫描分析，研究人员发现了一些染色体区域在品种比较中持续富集着差异表达基因，这些区域可能包含了导致两个品种胚胎发生能力差异的结构变异或调控元件。值得注意的是，一个位于1B染色体上约230 Mbp的区域，在三个比较阶段都富集了上调的差异表达基因，这可能是一个关键的基因组“热点”。
利用机器学习发现胚胎发生候选基因
为了从海量数据中精准锁定最具鉴别力的基因，研究人员应用了机器学习模型。结果显示，仅有两个基因被模型判定为最有效的分类特征。其中之一是组蛋白去乙酰化酶（Histone Deacetylase， HDAC）基因TRAESCS1D02G454400，该基因在Nanda的三个阶段均稳定上调。组蛋白去乙酰化酶是重要的表观遗传调控因子，通过移除组蛋白上的乙酰基团来影响染色质结构和基因表达。该基因位于一个富含差异表达基因的基因组区域内，并被预测为潜在的生物标志物，这强烈暗示其在调控小麦小孢子胚胎发生效率中扮演着关键角色。
研究结论与重要意义
这项研究通过对高、低响应小麦品种的系统性转录组比较，构建了一幅关于小麦小孢子胚胎发生的详细分子图谱。研究结论强调，胚胎发生的效率是由多层面因素协同决定的：
- 1.
表观遗传调控的核心作用：成功的胚胎发生起始于大规模、协调的表观遗传重编程，特别是染色质结构的动态变化。高响应基因型Nanda在关键时期更有效地激活了这些过程。
- 2.
亚基因组调控的复杂性：在小麦这种多倍体作物中，同源基因的表达平衡至关重要。研究首次揭示了B亚基因组基因的持续性抑制与低胚胎发生效率相关，为理解多倍体作物中基因剂量效应与性状调控提供了新视角。
- 3.
基因型特异性基因组变异的基础性影响：研究发现了多个与胚胎发生能力相关的基因组富集区域，表明品种间的遗传结构差异（如结构变异）是导致其不同响应能力的根本原因之一。
- 4.
多通路协同的网络化调控：胚胎发生不是由单个基因或通路决定的，而是应激感知、激素信号、细胞骨架重组、代谢重编程等多条通路协同作用的结果。高响应基因型展现出了更优的协同与时机控制。
该研究的重要意义在于，它超越了以往仅在一个响应品种内描述变化模式的研究，通过直接比较，揭示了决定“响应”与“难响应”的关键分水岭。研究所鉴定出的关键生物学过程、候选转录因子、基因组热点以及机器学习筛选出的生物标志物基因（特别是组蛋白去乙酰化酶基因TRAESCS1D02G454400），为未来通过分子标记辅助选择、基因编辑或表观遗传调控手段来改良难响应小麦品种，从而释放小孢子胚胎发生技术在小麦双单倍体育种中的全部潜力，提供了宝贵的路线图和具体靶点。这不仅有助于加速小麦育种进程，也对理解植物细胞全能性和发育可塑性这一基础生物学问题具有重要价值。
