一种跨人类细胞类型的转录基础模型

A foundation model of transcription across human cell types

作者信息Xi Fu, Shentong Mo, Alejandro Buendia, Anouchka P Laurent, Anqi Shao, Maria Del Mar Alvarez-Torres, Tianji Yu, Jimin Tan, Jiayu Su, Romella Sagatelian, Adolfo A Ferrando, Alberto Ciccia, Yanyan Lan, David M Owens, Teresa Palomero, Eric P Xing, Raul Rabadan
PMID39779852
期刊Nature
发布时间2025-01
DOI10.1038/s41586-024-08391-z

实验完整度

高

包含模型预训练与微调、跨细胞类型预测验证、B细胞PAX5突变的功能实验(Co-IP、BioID)及患者样本转录组分析,证据层级多样。

主要模型

胎儿及成体人类细胞类型(213种,包括成红细胞、B细胞等) 肿瘤细胞系(K562) GBM患者样本 REH B-ALL细胞系和HeLa细胞

重点核对

GET模型在未见过细胞类型(如胎儿星形胶质细胞)上的基因表达预测准确性(Pearson r = 0.94) 在胎儿成红细胞中预测胎儿血红蛋白相关顺式调控元件(BCL11A、NFIX、KLF1、HBG2) PAX5 G183S突变对PAX5-NR2C2相互作用的影响(BioID/Co-IP实验) K562 lentiMPRA数据的零样本预测性能 B-ALL患者样本中G183S突变特异转录程序的富集分析

摘要

转录调控涉及调控序列和蛋白质之间的复杂相互作用,指导着所有生物过程。现有的转录计算模型缺乏泛化能力,无法准确外推到未见过的细胞类型和条件。在此,我们介绍了GET(通用表达转换器),一种可解释的基础模型,旨在揭示跨213种人类胎儿和成体细胞类型的调控语法。仅依赖染色质可及性数据和序列信息,GET在预测基因表达方面达到了实验级别的准确性,即使在先前未见的细胞类型中也是如此。GET还展现出对新测序平台和检测方法的显著适应性,能够在广泛的细胞类型和条件下进行调控推断,并揭示了普遍性和细胞类型特异性的转录因子相互作用网络。我们评估了其在调控活性预测、调控元件和调控因子推断以及转录因子物理相互作用识别方面的性能,发现其在预测慢病毒大规模平行报告基因检测读数方面优于当前模型。在胎儿成红细胞中,我们鉴定了先前模型遗漏的远端(大于1兆碱基)调控区域;在B细胞中,我们鉴定了一种淋巴细胞特异性的转录因子-转录因子相互作用,该相互作用解释了白血病风险易感种系突变的功能意义。总之,我们提供了一个通用且准确的转录模型,以及具有细胞类型特异性的基因调控和转录因子相互作用目录。

实验结论

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

研究问题
如何构建一个可泛化到未见细胞类型和条件的转录调控模型,以揭示细胞类型特异性的调控语法和转录因子相互作用?
核心机制
GET基于染色质可及性和序列信息学习调控语法,预测基因表达;通过解释模型,揭示了转录因子之间的因果相互作用,并鉴定出细胞类型特异性的调节机制,如PAX5和NR2C2核受体域的相互作用介导白血病风险突变的影响。
主要证据
在左出星形胶质细胞上的预测与观测表达Pearson相关系数为0.94;在胎儿成红细胞中鉴定出BCL11A增强子,并优于Enformer等模型;Co-IP实验显示TFAP2A和ZFX相互作用;BioID实验显示PAX5 G183S突变增加与NR2C2的相互作用;患者样本转录组分析显示G183S突变特异性基因富集于NR/3和PAX5+NR/3靶基因。
研究意义
GET提供了一个可泛化的转录模型和转录因子相互作用的目录,有助于识别调控元件和上游调控因子,并为理解疾病相关突变的机制提供线索。

研究路径

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

1

模型构建与预训练

利用跨细胞类型的染色质可及性数据学习调控语法。

基于213种人类胎儿和成体细胞类型的scATAC-seq数据,构建GET模型,通过自监督掩码预训练预测掩蔽调控元件的基序结合分数和可及性分数。

2

微调预测基因表达

将预训练模型适配到基因表达预测,验证其泛化能力。

使用配对的多组学或独立scRNA-seq数据(153种细胞类型)微调GET,预测基因表达,并评估在未见过细胞类型上的表现。

3

零样本调控活性预测

评估GET对未见调控元件的活性预测能力。

使用K562 lentiMPRA数据,GET模型在未见lentiMPRA数据的情况下预测迷你启动子活性,并与Enformer对比。

4

顺式调控元件鉴定

从模型解释中鉴定调控基因表达的远端增强子和调控因子。

通过Jacobian归因方法计算区域和基序贡献分数,聚焦胎儿血红蛋白相关位点(BCL11A、NFIX、KLF1、HBG2),并与Enformer等模型比较。

5

转录因子相互作用发现

利用模型解释推断转录因子之间的功能性相互作用。

基于基因-基序矩阵进行相关性和因果发现(LiNGAM),构建TF-TF相互作用网络,与STRING数据库比较。

6

结构预测与实验验证

通过蛋白质结构预测和实验验证预测的TF-TF相互作用。

使用AlphaFold2预测TF-TF结构,并选择TFAP2A-ZFX和PAX5-NR2C2进行Co-IP和BioID验证。

7

患者样本转录组分析

验证PAX5 G183S突变对转录程序的影响。

分析B-ALL患者样本(141例,包括4例G183S种系突变),比较PAX5-WT、PAX5缺失和G183S突变组的差异表达基因。

研究方法

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

产品清单

实验环节名称品牌货号
DMEM培养基Gibco11965
胎牛血清HyCloneSH30070
NP-40裂解液----
蛋白酶和磷酸酶抑制剂混合物Sigma-AldrichPPC1010
TFAP2A抗体Santa Cruz Biotechnologysc-12726
Laemmli上样缓冲液BioRad1610737
Tris-甘氨酸凝胶ThermoFisherXP00100
PVDF膜Immobilon-PIPVH00010
TFAP2A抗体ABclonalA2294
ZFX抗体ThermoFisherPA5-34376
β-肌动蛋白抗体Santa Cruz Biotechnologysc-47778
SRF抗体AbclonalA16718
β-肌动蛋白抗体Cell Signaling Technology4967
pCDNA3.1-MCS-13Xlinker-BioID2-HA质粒Addgene80899
pCDH-GFP-puro载体System BioscienceCD513B-1
生物素Sigma-AldrichB4501
蛋白酶和磷酸酶抑制剂Life Technologies78443
苯zonase核酸酶Sigma-Aldrich70746-3
Pierce BCA蛋白定量试剂盒ThermoFisher Scientific23225
Dynabeads MyOne链霉亲和素C1磁珠Life Technologies65002
链霉亲和素-辣根过氧化物酶抗体Life TechnologiesS911
抗PAX5抗体Cell Signaling8970
抗HA抗体Cell Signaling3724
抗NR2C2抗体Cell Signaling31646
抗NCOR1抗体Cell Signaling5948
NRIP1-HRP抗体Santa Cruz Biotechnologysc-518071
NR3C1抗体Cell Signaling12041
ChimeraX----

关键环节

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

环节核对要点
模型预训练
使用的ATAC-seq数据集(Domcke等和Zhang等)、细胞类型数量、预训练轮数、掩码比例、优化器超参数
阅读提示:Methods: ATAC-seq data processing, Training scheme
基因表达预测
微调所使用的配对数据、训练和评估细胞类型划分、评估指标(Pearson r等)、输入窗口大小
阅读提示:Methods: RNA-seq data processing, Model evaluation
零样本MPRA预测
K562 NEAT-seq数据、lentiMPRA数据集、每次插入基因组的位置数、预测活性计算方法
阅读提示:Methods: LentiMPRA zero-shot prediction
顺式调控元件鉴定
使用的模型版本(binary ATAC)、Jacobian归因方法、HbFBase数据、K562 CRISPRi基准
阅读提示:Methods: Identifying important regions and regulators, Enhancer–gene pair prediction
TF相互作用发现
基因-基序矩阵的构建方法、LiNGAM因果发现参数、与STRING比较的截止值
阅读提示:Methods: Causal discovery of regulator interaction
Co-IP实验
细胞系(HeLa)、抗体品牌和编号、裂解液成分、孵育条件
阅读提示:Methods: TFAP2A coimmunoprecipitation
BioID实验
细胞系(REH)、质粒构建、生物素浓度和处理时间、链霉亲和素磁珠品牌和用量
阅读提示:Methods: Proximity labelling assay to detect PAX5–NR2C2 interactions
患者样本分析
样本数量(141例)、CDKN2A缺失选择、比较组定义、差异表达和富集分析方法
阅读提示:Methods: 未见明确详细Methods,但可参考Results和Extended Data Fig. 11