《Journal of Genetics and Genomics》：A single-nucleus transcriptome atlas of soybean anthers
编辑推荐：
为解决大豆花药发育过程中细胞类型特异性转录调控机制不明确、传统RNA-seq方法分辨率有限的问题，研究人员利用单核RNA测序（snRNA-seq）技术，构建了涵盖大豆花药多个发育阶段的单细胞转录组图谱。该研究鉴定并表征了9种细胞类型，揭示了从二倍体细胞到单细胞小孢子（UM）再到二细胞小孢子（BM）的关键转换节点中的动态转录重编程，发现了生殖细胞（GC）和营养细胞（VC）在RNA转运与稳定性调控上的功能分化。通过功能验证，证实了基因OSD1A和PKSA在花粉发育和育性中的关键作用。该高分辨率图谱为解析大豆花药发育的转录调控网络提供了宝贵资源，并为通过操纵雄性不育性以推动大豆杂交育种奠定了理论基础。
“民以食为天”，粮食安全关乎国计民生。大豆，作为最重要的油料作物和植物蛋白来源，其产量直接影响着全球的油脂和蛋白供应。然而，大豆生产正面临着日益严峻的挑战，特别是来自环境变化的压力。在气候变化的背景下，预计到2050年，全球大豆产量可能会在高排放情景下减少约2%。为了应对这一挑战，提高大豆产量和培育优良品种变得至关重要。这其中，植物的“生儿育女”——有性生殖过程扮演着核心角色，而花药作为雄性生殖器官，其正常发育直接决定了花粉的活力和最终作物的产量。
过去，科学家们通过批量RNA测序（bulk RNA-seq）等技术，已经对大豆花药的基因表达有了一定的了解，也发现了一些与雄性不育相关的候选基因。然而，花药是一个由多种不同类型细胞（如绒毡层细胞、表皮细胞、花粉母细胞、小孢子等）组成的复杂器官，批量RNA测序就像将一锅“细胞浓汤”混合后进行检测，无法分辨每种细胞独特的“语言”（基因表达谱），也无法精准捕捉不同细胞谱系在发育过程中的动态变化。为了“倾听”每一个细胞的声音，单细胞RNA测序（scRNA-seq）技术应运而生，并在拟南芥、玉米等植物中取得了成功。但这种方法需要制备原生质体，会带来细胞应激，可能改变细胞的真实状态。另一种技术，单核RNA测序（snRNA-seq），由于基于细胞核悬液，避免了细胞壁消化的步骤，在植物生物学领域显示出更广泛的适用性和稳定性。
针对大豆这种重要豆科作物，其花药发育过程的高分辨率、细胞类型特异性转录组图谱一直缺失。为了填补这一空白，并为提高大豆雄性育性、助力杂交育种提供分子层面的路线图，研究人员开展了一项系统性研究。他们利用snRNA-seq技术，绘制了覆盖大豆花药多个发育阶段的单细胞转录组图谱，深入解析了其发育过程中的转录调控机制，并鉴定了关键的功能基因。相关成果发表在《Journal of Genetics and Genomics》上。
研究人员为开展此研究，主要运用了以下关键技术方法：以大豆栽培品种Williams 82 (Wm82)为材料，根据花蕾长度（1-5 mm）分组取样，覆盖了从花粉母细胞到二细胞小孢子的不同发育阶段。利用单核RNA测序（snRNA-seq）平台（10x Genomics）对分离的花药细胞核进行高通量测序。对获得的超过13万个细胞核的转录组数据进行生物信息学分析，包括细胞聚类、细胞类型注释、差异表达基因分析、基因本体（GO）富集分析、拟时序（pseudotime）分析和共表达网络分析（WGCNA）。通过RNA荧光原位杂交（RNA-FISH）验证了细胞类型特异性标记基因的表达。利用CRISPR/Cas9基因编辑技术构建了OSD1A基因的敲除突变体，并通过EMS（乙基甲磺酸）诱变筛选并结合图位克隆（bulked-segregant analysis, BSA）鉴定了PKSA基因的突变体，对这两个基因在花粉发育和育性中的功能进行了表型验证。
研究结果
1. Identification and annotation of cell clusters from snRNA-seq of soybean anthers（基于大豆花药snRNA-seq数据的细胞簇鉴定与注释）
研究人员通过对不同长度花蕾（代表不同发育阶段）的花药进行snRNA-seq，获得了130,605个高质量细胞核的数据。聚类分析产生了28个初始簇，最终整合注释为9种主要的细胞类型，涵盖了体细胞谱系（绒毡层细胞、表皮细胞、药室内壁细胞）和生殖细胞谱系（花粉母细胞、减数分裂细胞、四分体细胞、单细胞小孢子、营养细胞和生殖细胞）。通过已知标记基因的表达和基因表达相关性热图，清晰地分离了这些细胞类型。对细胞类型比例的分析显示，测序结果与花药组织切片观察到的细胞类型动态变化趋势一致，验证了细胞注释的可靠性。
2. Functional profiling of cell-type-specific genes during soybean anther development（大豆花药发育过程中细胞类型特异性基因的功能分析）
通过分析在特定细胞类型中上调表达的基因，研究人员共鉴定出4921个细胞类型特异性基因。基因本体（GO）富集分析揭示了与各细胞类型已知功能相符的特征。例如，花粉母细胞富集了有丝分裂相关通路；减数分裂细胞富集了DNA重组相关通路；单细胞小孢子富集了细胞极性建立相关通路；生殖细胞富集了基因组稳定性（如组蛋白H2A泛素化）相关通路。此外，多个细胞类型共享了如花粉壁形成、细胞通讯、mRNA出核转运等通路。RNA原位杂交实验验证了在四分体、单细胞小孢子、营养细胞和生殖细胞中特异性上调表达的部分标记基因，证实了细胞注释和标记基因的可靠性。
3. Developmental trajectory of reproductive cells in soybean anthers（大豆花药中生殖细胞的发育轨迹）
对生殖细胞谱系进行拟时序分析，清晰地描绘了从花粉母细胞到减数分裂细胞、四分体细胞、单细胞小孢子，最终分化为营养细胞和生殖细胞的发育轨迹。分析发现，在发育过程中，转录本丰度（以UMI计数为指标）呈现动态变化。从减数分裂I期到II期，转录本丰度急剧下降至34.7%。在单细胞小孢子阶段，转录本丰度维持在较低水平（约为花粉母细胞阶段的67.3%）。进入二细胞小孢子阶段后，营养细胞的转录本丰度呈上升趋势并超过了花粉母细胞阶段，而生殖细胞的丰度则维持在类似单细胞小孢子的低水平。对相邻发育阶段基因表达的比较发现，从二倍体减数分裂细胞到单倍体单细胞小孢子的转换，伴随着大量新基因的激活，尽管总体转录本丰度降低。
4. Transcriptome changes during the transition from meiocyte to bicellular microspore（从减数分裂细胞到二细胞小孢子转换过程中的转录组变化）
深入分析从减数分裂细胞经四分体到单细胞小孢子这一关键转换过程，研究人员发现了独特的转录组变化模式。在单细胞小孢子中，与转录、DNA模板相关的基因（如bZIP和MADS转录因子）被激活，但同时“翻译的负调控”和“rRNA加工”等通路也被富集，表明此时细胞核转录增强，但mRNA出核转运和随后的翻译过程可能受到抑制。
进入二细胞小孢子阶段后，营养细胞和生殖细胞的转录程序出现显著分化。与单细胞小孢子相比，两者共同下调的基因远多于共同上调的基因。生殖细胞特异性上调的基因富集在染色质修饰和“多聚腺苷酸化mRNA出核转运”、“核转录mRNA多聚A尾缩短”等通路，后者涉及CCR4-Not复合体组分（如NOT4B, NOT3），表明生殖细胞依赖加速的RNA转运和降解来实现RNA周转。相反，营养细胞上调了与翻译、核糖体生物发生相关的通路，并高表达稳定mRNA的多聚A结合蛋白（PABs），显示出活跃的转录和翻译活动。
5. Soybean generative-cell genes experience weak selection pressure（大豆生殖细胞基因经历较弱的选择压力）
进化选择压力分析显示，在花粉发育过程中，四分体特异性基因受到最强的纯化选择。与营养细胞特异性基因相比，生殖细胞特异性基因受到的选择压力更弱。这可能是因为生殖细胞完全包裹在营养细胞的细胞质中，二者存在广泛的细胞通讯，缓冲了选择压力。此外，在驯化相关的选择性清除区域中，单细胞小孢子特异性基因的比例最高，表明在驯化过程中，单细胞小孢子（花药发育的早期单倍体阶段）承受了最强的选择压力。
6. Tapetum and microspore interactions vary across developmental stages（绒毡层与小孢子的相互作用随发育阶段而变化）
研究人员对绒毡层细胞进行了深入分析。根据拟时序和分化状态预测，将绒毡层细胞分为早期（C13）和晚期（C9）两个阶段。早期绒毡层细胞高表达有丝分裂细胞周期过程相关基因，而晚期绒毡层细胞则高表达脂质生物合成、花粉壁组装和孢粉素生物合成过程相关基因，包括MS2、CYP703A2、CYP704B1、PKSA和PKSB等关键酶基因。
通过加权基因共表达网络分析（WGCNA），在绒毡层细胞中构建了基因调控网络，不仅验证了已知的花粉育性调控核心通路（DYT1-TDF1-AMS-MS188-MS1），还系统性地鉴定了与孢粉素生物合成相关的保守基因集。此外，还发现了多个与雄性不育相关的新候选基因，如Glyma.06G285100、Glyma.12G120700、Glyma.15G271600、Sec8 (Glyma.20G182800)和Glyma.08G183500等。
7. Functional analysis of OSD1A and PKSA validates the snRNA-seq-inferred expression patterns in soybean anther cell types（OSD1A和PKSA的功能分析验证了snRNA-seq推断的大豆花药细胞类型表达模式）
为验证snRNA-seq鉴定的关键基因，研究人员对OSD1A和PKSA进行了功能分析。OSD1A是一个细胞周期调控因子，在减数分裂II期和四分体中高表达。通过CRISPR/Cas9技术构建的osd1a敲除突变体，花粉活力下降约12%，约20%的花粉在减数分裂II期末期染色体分离失败，产生染色体数目异常的二分体，并导致植株结实率严重下降（>90%），植株矮化和成熟延迟。
PKSA是孢粉素合成途径中的聚酮合酶，在晚期绒毡层细胞中特异性高表达。通过EMS诱变和BSA图位克隆鉴定的pksa突变体（E738D），表现出完全的雄性不育。其花粉粒粘连、外壁网状结构被破坏、无活力，植株不能结荚。RNA原位杂交证实了PKSA在绒毡层细胞中的特异性表达。这些结果验证了snRNA-seq数据推断的细胞类型特异性表达模式的准确性，并证实了这两个基因在大豆花粉发育和育性中的关键作用。
结论与讨论
本研究利用单核RNA测序（snRNA-seq）成功构建了一个全面、高分辨率的大豆花药发育转录组图谱，涵盖了超过13万个细胞和9种不同的细胞类型。这张图谱不仅提供了详细的细胞特异性基因表达谱，还系统揭示了从二倍体孢子体到单倍体配子体转换过程中动态的转录重编程事件。研究发现了在单细胞小孢子阶段，尽管总体转录本丰度降低，却有大量新基因被激活，同时mRNA转运和翻译过程可能受到抑制，暗示这可能是一个为后续配子体基因表达快速启动做准备的“蓄力”阶段。进入二细胞小孢子阶段后，营养细胞和生殖细胞展现出截然不同的转录程序分工：营养细胞转向活跃的翻译机制，以支持花粉管生长和营养供应；而生殖细胞则专注于mRNA的快速出核转运和降解，可能与降低突变风险、保证遗传物质精确传递有关。进化分析进一步表明，生殖细胞特异性基因经历了比营养细胞基因更弱的纯化选择，而单细胞小孢子特异性基因在驯化过程中经历了最强的正向选择。
研究还深入解析了绒毡层细胞在不同发育阶段的功能转换及其与生殖细胞的互作，并通过基因共表达网络鉴定了包括已知核心通路和一系列新候选基因在内的雄性不育相关基因资源库。最后，通过对OSD1A和PKSA基因的功能验证，不仅确认了snRNA-seq图谱的可靠性，也直接证明了这些关键基因在调控大豆减数分裂、孢粉素合成及雄性育性中的决定性作用。
这项研究的意义在于，它首次为重要豆科作物大豆的花药发育提供了单细胞水平的分子图谱，极大地深化了对大豆雄性生殖发育调控网络的理解。所鉴定的细胞类型特异性标记基因、关键调控通路和候选雄性不育基因为后续的功能研究和分子设计育种提供了宝贵的靶点资源。特别是在全球粮食安全挑战和气候变化压力日益增大的背景下，利用这些研究成果，结合先进的基因编辑技术，有望精准创制理想的大豆雄性不育系，从而突破大豆杂交育种的瓶颈，充分利用杂种优势，为培育高产、稳产、抗逆的大豆新品种，保障未来全球蛋白质供应提供强有力的科技支撑。
