时空图谱分析界定胚胎肺形态发生中的上皮–间充质转化窗口

Spatiotemporal Profiling Defines the Epithelial and Mesenchymal Transition Window in Embryonic Lung Morphogenesis.

作者信息Huiwen Zheng, Jinpei Lin, Hanyi Li, Shijie Hao, Mengnan Cheng
PMID42346508
发布时间2026-06-01
DOI10.3390/jdb14020025

实验完整度

中

整合空转与单细胞数据,行解卷积、PCA、CellChat、NicheNet、pySCENIC、gsMap及cBioPortal基因组分析,以计算分析为主,无湿实验功能验证。

主要模型

小鼠胚胎肺组织空间转录组(E10.5–E18.5) 小鼠肺来源单细胞RNA测序数据集(26万细胞) TP53/HNRNP突变型肺腺癌患者基因组数据

重点核对

空间转录组bin50解卷积与CytoSpace单细胞注释对比,验证细胞类型空间分布合理性 五个间充质亚群(C0–C4)标记基因表达及比例随发育阶段变化 E13.5前后早期与晚期上皮/间充质差异表达基因谱 CellChat推断BMP、FGF、PDGF、ADGRG通路配体–受体互作(差异阈值log2FC>0.5,互作p<0.01) gsMap整合GWAS统计量与空间转录组评估肺癌性状空间关联

摘要

肺的器官发生由胚胎发生过程中动态的上皮–间充质相互作用所协调,然而控制假腺体期这些过程的基因调控程序和信号动力学仍未被完全理解。在本研究中,我们整合了跨胚胎发育阶段的空转和单细胞转录组数据,以系统表征肺发育过程中的上皮和间充质动力学。为了在空转数据中获得单细胞分辨率下更精细的细胞类型,我们开发了一种基于bin的解卷积策略,实现了高精度的细胞类型分配。随后,我们构建了肺发育的3D时空图谱,并阐明了肺形态发生过程中上皮–间充质成熟所依赖的分子调控机制。此外,我们分析了转录因子模块活性、细胞间通讯信号,并预测了下游靶基因,同时整合公开GWAS元数据将发育程序与肺癌相关特征联系起来。我们观察到假腺体期与晚期胚胎阶段之间存在明显的阶段特异性功能异质性。值得注意的是,E13.5成为一个关键的转化窗口,在此期间祖细胞状态向更成熟的细胞表型转变。我们重建了上皮–间充质相互作用,并揭示了跨发育阶段的配体–受体信号和转录网络的协调重塑。调控网络分析进一步鉴定出以Tbx3、Tbx5、Gli1、Gata4/5、Foxa1/2和Cebpa为中心的时间协调转录因子模块,它们共同协调了分支形态发生、上皮模式化和组织稳定。与肺癌全基因组关联数据的整合表明,胚胎肺祖细胞状态与肺癌相关转录程序存在强关联,特别是涉及上皮–间充质可塑性和RNA剪接通路。此外,TP53/HNRNP突变型肺腺癌显示出与细胞骨架重塑和祖细胞状态再激活相关的胚胎样分子特征。总之,我们的研究提供了一个时空分辨的胚胎肺发育框架,并确定了一个连接肺形态发生、调控网络重塑和癌症相关上皮可塑性的关键转化窗口。

实验结论

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

研究问题
胚胎肺假腺体期上皮–间充质相互作用及调控程序如何时空动态变化,其与肺癌相关特征是否存在关联。
核心机制
E13.5为上皮与间充质同步成熟的转化窗口,Tbx3驱动的早期程序逐步被Tbx5和Gli1依赖网络取代,协调细胞骨架重塑与ECM组织;胚胎肺祖细胞状态与肺癌相关转录程序(尤其上皮–间充质可塑性和RNA剪接)相关。
主要证据
基于小鼠E10.5–E18.5空转与单细胞数据构建时空图谱,识别五个间充质亚群及上皮/间充质阶段差异表达;CellChat、NicheNet和pySCENIC揭示信号与转录因子模块动态;gsMap关联肺癌GWAS性状,cBioPortal分析TP53/HNRNP突变LUAD基因组改变。
研究意义
提供时空分辨的胚胎肺发育框架,为理解肺疾病分子基础建立发育框架,并对再生医学和治疗干预策略具有潜在意义。

研究路径

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

1

空间转录组细胞分割与细胞类型注释

在单细胞分辨率下对胚胎肺空转数据进行准确细胞类型注释

使用StereoCell进行细胞分割,将bin50数据与CytoSpace结合,并开发基于最优传输的bin解卷积策略,将空间bin分配给参考单细胞类型

2

鉴定上皮与间充质成熟的时间转化点

确定上皮和间充质群体从祖细胞向晚期胚胎阶段转变的关键时间点

对间充质细胞进行PCA和UMAP聚类,注释五个间充质亚群;比较E13.5前后亚群比例和基因表达变化;对早期(E10.5–E13.5)与晚期(E18.5)细胞进行差异表达和功能富集分析

3

上皮–间充质信号通讯与转录因子调控网络分析

解析胚胎期至假腺体期细胞间信号重塑及配体–转录因子调控关系

使用CellChat推断细胞间通讯,比较胚胎期与假腺体期信号输入输出变化;使用NicheNet预测配体–靶基因调控关系;使用pySCENIC推断转录因子调控网络和regulon活性

4

肺癌相关性状的遗传信息空间映射

评估胚胎肺发育细胞状态与肺癌遗传风险的关联

使用gsMap整合GWAS汇总统计量与空间转录组数据,计算基因特异性得分和性状遗传力富集;提取E10.5–E18.5单细胞数据计算肺癌性状GSS并量化上皮/间充质状态评分

5

TP53/HNRNP突变型肺腺癌基因组改变分析

比较TP53突变与TP53/HNRNP共突变LUAD样本的基因组共改变景观

通过cBioPortal对265例TP53突变LUAD按HNRNP家族基因突变状态分层,进行比较基因组改变分析

研究方法

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

关键环节

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

环节核对要点
空间转录组细胞类型注释
空转数据来源与样本纳入(E12.5可识别肺组织,E13.5因缺ssDNA图像被排除);bin50分辨率;CytoSpace默认参数;参考单细胞每种类型最多500个子采样细胞
阅读提示:Methods 2.1 Data Source、2.3 Cell-Type Annotation;Figure S2
上皮–间充质转化窗口鉴定
间充质亚群C0–C4标记基因定义;模块评分阈值>0.2;早期(E10.5–E13.5)与晚期(E18.5)差异分析分组
阅读提示:Results 3.2;Figure 2;Figure S3
细胞通讯与配体–靶基因预测
CellChat v1.6.1默认参数及差异阈值log2FC>0.5;NicheNet配体筛选log2FC>0.5且调整p<0.05,配体–靶基因权重阈值0.5;pySCENIC默认参数
阅读提示:Methods 2.5–2.7;Figure 3;Figure 4
肺癌GWAS空间关联
GWAS数据集ukb-b-14521;S-LDSC与Cauchy组合检验参数;上皮/间充质状态评分基因集
阅读提示:Methods 2.9;Results 3.6;Figure 5;Figure S5
TP53/HNRNP突变LUAD分层
总样本265例,TP53组203例,HNRNP组62例;数据集来源TCGA PanCanAtlas
阅读提示:Methods 2.10;Figure S5E