睡莲花器官发生的时空转录组图谱为基部被子植物花瓣的演化起源提供见解

Spatiotemporal transcriptome atlas of water lily floral organogenesis provides insight to the evolutionary origin of petals in basal angiosperms

作者信息Chang Liu, Ge Wang, Binbin Cao, Yonglong Li, Yamao Chen, Chunce Guo, Ji Qi
PMID42972450
发布时间2026-08-05
DOI10.1186/s13059-026-04229-1

实验完整度

高

空间转录组测序五张切片、RNA原位杂交验证及跨12种植物系统发育基因组学分析,含功能与机制验证

主要模型

睡莲(Nymphaea colorata)花芽发育四个阶段的空间转录组切片 蝴蝶兰(Phalaenopsis aphrodite)花芽对照比较样本 睡莲花器官的bulk RNA-seq数据

重点核对

花芽取自幼莲育种基地的N. colorata,覆盖四个发育阶段,共五张纵向切片(S1、S2、S3、S4-1、S4-2) 空间转录组测序在BMK S1000平台完成,cDNA于Illumina Novaseq平台以2×150 bp双端策略测序 组织与空间域识别:S1-S3用STEEL(S1参数--gini_read=5、--gini=0.1;S2/S3 --gini_read=10、--gini=0.1;20个主成分),S4用Seurat两轮聚类 MADS-box四聚体模型基于AGL6、AP3、PI、AG、SEP表达量以整数线性规划建模,并对表达值±20%扰动重复1000次评估稳健性 RNA原位杂交针对NcAGL6、NcAP3、NcPI、NcAG1、NcAG2在花器官决定关键阶段进行验证

摘要

背景:花瓣是花朵的关键演化创新,重塑了植物-传粉者相互作用,并促成了被子植物在陆地生态系统中的优势地位。然而,由于以往研究主要依赖形态学推断,且花器官发生过程中空间分辨的分子证据有限,花瓣的演化起源仍存争议。结果:本文对基部被子植物睡莲(Nymphaea colorata)发育中的花芽进行空间转录组测序,以解析花瓣形成所依赖的潜在基因调控格局。分析揭示出四个分生细胞群呈现出意想不到的层级结构,这些细胞群在花基部持续存在直至心皮形成,并沿近端-远端轴表现出分生潜力递减的梯度。其中一群细胞发育为雄蕊和花瓣(内轮花被片)。发育轨迹重建进一步揭示了从雄蕊原基到发育中花瓣的连续细胞命运转变。这些转变由MADS-box四聚体组成的定量变化驱动,逐步降低外轮花器官的生殖身份。涉及数千个组织偏好性表达基因的发育调控模块的空间定位表明,花瓣形态发生整合了额外的类叶遗传程序,最终将花瓣与雄蕊区分开来。结论:综上所述,这些结果揭示了花瓣通过招募种子植物之前已建立的祖先调控网络而从生殖器官演化衍生而来,为开花植物中花瓣状器官的发育基础和演化创新提供了新见解。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
基部被子植物睡莲中花瓣(内轮花被片)的发育基础及其演化起源是什么。
核心机制
内轮花被片与雄蕊起源于共同的分生组织祖细胞群,MADS-box四聚体组成的定量(剂量)变化逐步降低生殖身份并促进花瓣身份,随后招募极性、光合等类叶调控模块。
主要证据
五张花芽纵向切片的空间转录组测序结合轨迹分析、MADS-box四聚体整数线性规划建模、RNA原位杂交验证,以及跨12种绿色植物系统发育基因组学比较。
研究意义
揭示花瓣通过招募种子植物出现之前已建立的祖先调控网络而从生殖器官演化衍生,为被子植物花瓣状器官的发育基础与演化创新提供新见解。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

花芽取样与空间转录组测序

构建睡莲早期花器官发生的高分辨率时空基因表达图谱。

采集四个发育阶段的N. colorata花芽,固定、蔗糖处理、OCT包埋并液氮冷却异戊烷速冻,切片后经HE染色、透化、荧光cDNA合成与成像,在BMK S1000平台进行空间转录组测序,cDNA在Illumina Novaseq平台测序。

2

组织与细胞类型识别

为组织注释提供基于转录组的框架,界定花器官各组织域。

整合空间聚类、标记基因表达与组织学信息,S1-S3采用STEEL,S4采用Seurat两轮聚类,结合已知标记基因与HE染色组织学特征注释细胞类型。

3

基部花组织分生细胞群异质性分析

解析花基部多个分生细胞群的组织方式与发育轨迹。

提取S1-S3中位于花基部、高表达NcKNAT1、NcSTM、NcKNAT2且低表达分化基因的簇,用Seurat整合后以Monocle进行轨迹分析(DDRTree降维),重建M1–M4四个分生细胞群的发育轨迹。

4

MADS-box表达与四聚体组成分析

探究内轮花被片与雄蕊身份转变的分子基础。

比较雄蕊1–4与内轮/外轮花被片的基因表达相似性,检测SEP、AGL6、AP3、PI、AG1、AG2空间表达,并以整数线性规划模型在最大化四聚体总量约束下估计各器官四类四聚体(SeIT、PeIT、StIT、CaIT)的比例。

5

RNA原位杂交验证

独立验证花器官身份决定关键MADS-box基因的空间表达模式。

对NcAGL6、NcAP3、NcPI、NcAG1、NcAG2在花器官决定关键阶段进行RNA原位杂交,使用地高辛标记约30 bp反义与正义寡核苷酸探针。

6

内轮花被片与雄蕊分化轨迹及调控模块分析

阐明内轮花被片与雄蕊分化所招募的调控模块。

收集雄蕊1–4与内轮花被片细胞用Monocle 2进行轨迹分析,鉴定假时间相关基因并分为三个模块,进行GO富集与激素信号共表达网络分析。

7

跨物种系统发育基因组学比较

探究花器官发生相关基因的演化起源与保守性。

对12种绿色植物进行BLASTP同源鉴定与PhyloMCL同源群预测,用Dollo简约法推断基因家族在各分类节点的获得/丢失,重构最大似然树。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
BMK S1000平台BMK--
Illumina Novaseq平台Illumina--
石蜡Sigma-AldrichP3683
Leica RM2235切片机LeicaRM2235
蛋白酶KPINUOFEIPN0012
NBT/BCIP显色底物BOSTERAR1023
Nikon倒置荧光显微镜NIKONECLIPSE Ts2

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
花芽取样与空间转录组测序
取样地点(海南睡莲育种基地)、四个发育阶段划分、五张切片编号与固定包埋流程;RIN合格标准的具体数值文中未明确说明
阅读提示:Methods之Plant tissue dissection and spatial transcriptome sequencing;Fig. 1A;Additional file 1: Fig. S1
组织与细胞类型识别
S1-S3使用STEEL及S4使用Seurat的参数与聚类分辨率;各阶段最终分组数
阅读提示:Methods之Tissue identification from spatial transcriptomic datasets;Additional file 2: Table S2
分生细胞群轨迹分析
用于轨迹分析的cluster选取标准与Monocle参数
阅读提示:Methods之Trajectory analysis of floral tissues;Fig. 2B-2C
MADS-box四聚体建模
四类四聚体公式、最大化目标函数及±20%扰动1000次的稳健性检验设置
阅读提示:Methods之Construction of MADS box tetramer dynamic model;Fig. 3G
RNA原位杂交
固定时间、蛋白酶K浓度与消化时间、探针长度及杂交温度(杂交温度具体数值文中未明确说明)
阅读提示:Methods之RNA in situ hybridization;Additional file 1: Fig. S10;Additional file 2: Table S14
跨物种比较与系统发育分析
12个物种组成、同源群预测与Dollo简约法参数、RAxML模型与bootstrap次数
阅读提示:Methods之Ortholog identification and reconstruction of phylogenetic trees;Fig. 5A