早期拟南芥种子发育的转录图谱揭示组织间协调机制

成功的种子发育对于被子植物繁殖至关重要，并且需要三种遗传上彼此不同的组织之间的协调：受精产物胚（embryo）和胚乳（endosperm），以及母体来源的种皮（seed coat）。目前，对于支撑种子组织特异性功能及组织间协调的转录程序，认识仍不完整。为解决这一问题，研究人员对授粉后3、5和7天的拟南芥（Arabidopsis thaliana）种子进行了单核RNA测序（single-nucleus RNA sequencing, snRNA-seq）。研究表征了所有主要的种子细胞或细胞核类型，进一步细化了胚乳中的转录状态，并绘制了细胞类型特异性基因上的选择作用特征。除其他发现外，分析结果揭示了参与油菜素内酯（brassinosteroid, BR）响应性转录因子激活相关基因的区室化分布、编码短分泌肽（short, secreted peptides, SSPs）基因在胚乳中的丰富表达，以及快速进化基因在胚乳和种皮亚型中的富集表达，从而阐明了种子基因的细胞类型特异性与物种特异性。

该研究发表于《Nature Plants》，围绕早期拟南芥种子发育过程中不同组织之间如何实现精细协调这一核心问题，构建了一个具有时间分辨率和细胞核分辨率的转录图谱。种子由胚、胚乳和母体种皮三类遗传背景不同的组织共同构成，三者在发育过程中必须严格同步，才能保证后代正常形成。既往研究已经指出，种皮发育依赖胚乳来源的生长素，胚形态建成又受早期种皮来源生长素调控，说明种子内部存在广泛而复杂的跨组织信号交流。然而，由于这些组织之间被细胞壁和膜结构分隔，且分别构成共质体隔离区（symplastic fields），因此其协调机制一直缺乏系统性的细胞类型解析。尤其是在胚乳区室化、合点端种皮（chalazal seed coat, CZSC）功能分工、短分泌肽（SSPs）介导的跨组织通讯以及快速进化基因在种子中的空间分布等方面，原有认识仍较为有限。因此，开展该研究的必要性在于，以高分辨率转录组手段重建早期种子各组织和细胞核类型的分子状态，从而解析种子组织特异性功能及其协调的潜在机制。

研究人员以拟南芥Col-0为材料，选取授粉后3、5、7天这一起始胚乳最为动态的发育窗口，建立了早期种子单核转录图谱。研究共获得54,210个细胞核转录谱，覆盖胚、胚乳、种皮以及少量未受精胚珠和珠柄组织，并对不同时间点数据进行了聚类、分层注释和整合分析。基于已知标记基因、差异表达分析、基因本体（Gene Ontology, GO）富集及杂交链式反应RNA荧光原位杂交（hybridization chain reaction RNA-fluorescence in situ hybridization, HCR RNA-FISH）验证，研究系统界定了大多数已知种子细胞类型，并在胚乳和下部种皮中解析出更精细的亚群。最终结果表明：早期拟南芥种子在空间上具有高度区室化的转录程序；合点端种皮不同亚型在营养转运、胼胝质代谢和激素相关过程中具有互补功能；珠孔端区域存在BR生物合成和应答相关基因的集中表达；胚乳中尤其是珠孔胚乳（micropylar endosperm, MCE）和合点胚乳（chalazal endosperm, CZE）富集大量SSP编码基因；此外，外周胚乳（peripheral endosperm, PEN）和CZE表达了多数呈正选择信号的种子基因。这些发现为理解种子内部跨组织信号、营养分配和基因演化提供了高分辨率框架。

主要技术方法方面，研究人员采用10x Genomics平台对授粉后3、5、7天拟南芥Col-0种子分离细胞核后进行snRNA-seq；每个时间点设置2个生物学重复，样本来源为每次10–15个角果中分离的约500–800粒种子。随后结合Seurat、Harmony、monocle3等进行数据整合、聚类注释、拟时序与模块评分分析，并通过差异表达、GO富集、相关性分析和正选择位点分析解析功能特征。空间定位方面，采用HCR RNA-FISH对关键标记基因进行原位验证。进化层面，利用OrthoFinder、codeml/PAML和InterProScan对单拷贝直系同源基因（single-copy orthologues, SCOs）进行正选择与结构域定位分析。

在研究结果方面，论文依次呈现了以下几个核心部分。

A transcriptional atlas of early Arabidopsis seed development
研究首先构建了早期拟南芥种子发育的单核转录图谱。通过对3、5、7 DAP样本进行snRNA-seq，研究人员在最高分辨率下分别鉴定出34、33和25个簇，并将其归并为跨时间点可比的层级注释体系。图谱共包含54,210个核，其中种皮占比最高，胚乳次之，胚所占比例较小。研究进一步验证了图谱完整性，检测到悬垂体（suspensor）等稀有胚细胞类型，但部分已报道的茎端分生组织（shoot apical meristem, SAM）细胞类型未被充分捕获。总体而言，该图谱覆盖了绝大多数已知的解剖学种子细胞类型，为后续功能解析提供基础。

A high-resolution census of the developing seed coat
在种皮部分，研究建立了发育中种皮的高分辨率细胞类型普查。研究识别出一般种皮（general seed coat, GSC）的五层结构，并重点解析了下部种皮这一高度异质区域。通过ARR22等标记，研究鉴定出18个ARR22⁺簇，涉及GSC、CZSC及合点增殖组织（chalazal proliferating tissue, CPT）。其中，CPT被进一步区分为“持续型”和“瞬时型”亚群。差异表达与GO分析显示，持续型CPT富集核苷分解代谢相关基因，符合该区正在发生程序性细胞死亡（programmed cell death, PCD）的特征；而两类CPT在细胞死亡触发机制和金属响应等方面存在差异。这说明下部种皮并非单一功能结构，而是由多个具有分子专属性的亚型组成。

Cell types in the placentochalazal region of the CZSC have complementary functions
在CZSC的胎座-合点区（placentochalazal region），研究发现不同细胞亚型承担互补功能。依据UMAMIT14、SWEET10、TET5和AAP6等标记，研究鉴定出多个CZSC亚型。TET5⁺亚群差异基因富集于“(1→3)-β-ᴅ-glucan（callose）代谢过程”，并特异表达CALS8等胼胝质生物合成相关基因，提示其主要参与胼胝质沉积和胞间连丝通透性调控。相对地，SWEET10⁺亚群富集磷酸盐转运和磷酸酶活性相关基因，如PHO1;H1及多个海藻糖磷酸酶基因，并伴随cwINV4表达，提示其可能是营养输入和糖利用调控的主要部位。由此可见，胎座-合点区不同细胞类型在母体资源输入过程中形成了功能互补。

BR biosynthesis, homeostasis and response genes show concentrated expression in the micropylar region of the seed
研究进一步分析了油菜素内酯相关基因的空间分布。模块评分结果显示，BR生物合成在ii2层最显著，其次为外珠被1层（outer integument 1, oi1）；而更高分辨率分析揭示，这种富集主要由一个珠孔端oi亚型驱动。该亚型以SWEET12表达为特征，并高表达催化BR生物合成最后两步的BR6OX1和BR6OX2，提示其是早中期种子发育中BR产生的关键部位。与此同时，BZR1/BES1家族转录因子在5 DAP的MCE中高度富集，而BR受体BRI1在该处较低，说明BR应答可能与邻近珠孔端种皮信号输入有关。BES1靶基因从3到7 DAP逐步上调，支持珠孔端区域BR信号向胚乳传播并参与局部发育调控的观点。

The micropylar-to-embryo surrounding region endosperm shift is characterized by UMAMIT and nitrate transporter gene expression
对于胚乳发育动态，研究指出5到7 DAP之间发生了显著转录转变。7 DAP时，GLIP6⁺的MCE，即此阶段的胚周区域（embryo surrounding region, ESR），成为最具转录独特性的群体，其特征是UMAMIT家族和硝酸盐转运体基因的上调。该区域已知与胚角质层和胚鞘形成相关，并受ZOU/RGE1调控。研究在ZOU⁺群体中识别出两个状态：一个表达调控胚鞘形成的KRS，另一个特异表达NAC074并富集NAC087，后者与胚乳程序性细胞死亡相关。这些结果表明，7 DAP珠孔端/胚周胚乳内部存在与屏障建成和细胞死亡相关的连续状态变化。

The developmental basis for transcriptional polarity in the CZE
围绕CZE，研究揭示了其转录极性的发育基础。CZE是不发生细胞化的特殊胚乳区域，位于母体-子代界面。借助大量CZE细胞核数据，研究重新审视了囊（cyst）、结节（nodule）和类结节（nodule-like）状态。RALFL3被鉴定为3 DAP某一CZE亚群的高度特异标记，HCR RNA-FISH显示其信号位于合点囊底部或遍布囊体；相对地，AT3G49307在囊顶部和结节中更强。基于此，研究提出合点囊内部存在顶端-基部转录轴，RALFL3和AT3G49307分别标记基部与顶端区域。拟时序分析以3 DAP PEN为起点，显示通向基部囊与其余CZE亚型的是两条不同分支，说明基部囊具有独立发育轨迹。结合相关性分析，基部囊与邻近母体组织ii1和CPT转录上更相似。研究因此提出RALFL3⁺细胞核可能是最早迁移至合点端的“奠基者”细胞核，后续核融合形成成熟囊结构。

Discrete families of SSPs are enriched in endosperm subtypes
在SSP方面，研究表明若干离散家族的SSP在胚乳亚型中显著富集。通过对拟南芥SSP进行系统检测和家族注释，研究发现防御素样（defensin-like, DEFL）、低分子量富半胱氨酸（low-molecular-weight cysteine-rich, LCR）和脂质转运（lipid transfer, LTP）是图谱中最主要的SSP家族。3 DAP MCE、3 DAP CZE和5 DAP CZE是SSP表达枢纽，无论在表达量还是独特SSP数量上均最突出。富集的SSP不仅来自具有细胞间信号功能的TPD、CLE和LCR家族，也包括具有抑制功能的果胶甲酯酶抑制蛋白（pectin methylesterase inhibitor, PMEI）家族。许多此类SSP在受精后上调，提示其很可能参与种子发育期而非仅限于受精过程的跨组织通讯。

Rapidly evolving single-copy orthologues are compartmentalized in the endosperm
在进化分析部分，研究检测了拟南芥与Arabidopsis lyrata、Arabidopsis arenosa、Capsella grandiflora共有的7,187个SCO中的正选择信号，发现141个种子基因具有显著M2a/M1a似然比检验（likelihood ratio test, LRT）结果，其中103个在种子不同簇间差异表达。总体上，胚乳亚型差异表达了最多的快速进化基因，而ii1′和ii2等种皮层则对少数此类基因呈现高表达富集。进一步结构域分析显示，359个高置信正选择位点中，相当比例位于分泌蛋白胞外结构域或信号肽中，其次位于内在无序区（intrinsically disordered regions, IDRs）。这些结果支持种子中尤其是胚乳富含快速进化基因，并提示与细胞外信号传导和蛋白互作相关的结构区域可能是种子快速进化的重要靶点。

在讨论部分，论文强调，该研究提供了一个覆盖早期种子主要组织和细胞核类型的高分辨率转录图谱，深化了对细胞间信号、功能区室化和基因多样化的理解。研究不仅补充了营养转运、胼胝质沉积和BR生物合成定位等方面的新证据，也通过高覆盖率数据识别出CZE中的稀有细胞核群体，从而对合点囊起源提出新的发育模型。对于SSP，研究指出过去图谱由于探针和分辨率限制，可能低估了其在种子发育中的作用；当前数据则显示MCE和CZE是SSP高度富集的信号界面。对于快速进化基因，研究确认胚乳中存在显著富集，并指出胞外结构域和IDRs中的正选择位点可能与种子内外信号转导有关。同时，作者也说明其分析限于SCO，因此未覆盖扩张基因家族成员。

论文结论部分可概括翻译为：研究人员构建了一个具有高分辨率注释的早期种子发育转录图谱，为未来研究早期种子发育奠定了基础。为此，研究团队还建立了在线资源，用于数据浏览与分析。