模型构建与预训练
利用跨细胞类型的染色质可及性数据学习调控语法。
基于213种人类胎儿和成体细胞类型的scATAC-seq数据,构建GET模型,通过自监督掩码预训练预测掩蔽调控元件的基序结合分数和可及性分数。
A foundation model of transcription across human cell types
包含模型预训练与微调、跨细胞类型预测验证、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突变特异转录程序的富集分析
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
利用跨细胞类型的染色质可及性数据学习调控语法。
基于213种人类胎儿和成体细胞类型的scATAC-seq数据,构建GET模型,通过自监督掩码预训练预测掩蔽调控元件的基序结合分数和可及性分数。
将预训练模型适配到基因表达预测,验证其泛化能力。
使用配对的多组学或独立scRNA-seq数据(153种细胞类型)微调GET,预测基因表达,并评估在未见过细胞类型上的表现。
评估GET对未见调控元件的活性预测能力。
使用K562 lentiMPRA数据,GET模型在未见lentiMPRA数据的情况下预测迷你启动子活性,并与Enformer对比。
从模型解释中鉴定调控基因表达的远端增强子和调控因子。
通过Jacobian归因方法计算区域和基序贡献分数,聚焦胎儿血红蛋白相关位点(BCL11A、NFIX、KLF1、HBG2),并与Enformer等模型比较。
利用模型解释推断转录因子之间的功能性相互作用。
基于基因-基序矩阵进行相关性和因果发现(LiNGAM),构建TF-TF相互作用网络,与STRING数据库比较。
通过蛋白质结构预测和实验验证预测的TF-TF相互作用。
使用AlphaFold2预测TF-TF结构,并选择TFAP2A-ZFX和PAX5-NR2C2进行Co-IP和BioID验证。
验证PAX5 G183S突变对转录程序的影响。
分析B-ALL患者样本(141例,包括4例G183S种系突变),比较PAX5-WT、PAX5缺失和G183S突变组的差异表达基因。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | DMEM培养基 | Gibco | 11965 |
| 胎牛血清 | HyClone | SH30070 | |
| 细胞裂解 | NP-40裂解液 | -- | -- |
| 蛋白酶和磷酸酶抑制剂混合物 | Sigma-Aldrich | PPC1010 | |
| 免疫共沉淀 | TFAP2A抗体 | Santa Cruz Biotechnology | sc-12726 |
| Laemmli上样缓冲液 | BioRad | 1610737 | |
| Tris-甘氨酸凝胶 | ThermoFisher | XP00100 | |
| PVDF膜 | Immobilon-P | IPVH00010 | |
| TFAP2A抗体 | ABclonal | A2294 | |
| ZFX抗体 | ThermoFisher | PA5-34376 | |
| β-肌动蛋白抗体 | Santa Cruz Biotechnology | sc-47778 | |
| SRF抗体 | Abclonal | A16718 | |
| β-肌动蛋白抗体 | Cell Signaling Technology | 4967 | |
| 邻近标记 | pCDNA3.1-MCS-13Xlinker-BioID2-HA质粒 | Addgene | 80899 |
| pCDH-GFP-puro载体 | System Bioscience | CD513B-1 | |
| 生物素 | Sigma-Aldrich | B4501 | |
| 蛋白酶和磷酸酶抑制剂 | Life Technologies | 78443 | |
| 苯zonase核酸酶 | Sigma-Aldrich | 70746-3 | |
| Pierce BCA蛋白定量试剂盒 | ThermoFisher Scientific | 23225 | |
| Dynabeads MyOne链霉亲和素C1磁珠 | Life Technologies | 65002 | |
| Western blot | 链霉亲和素-辣根过氧化物酶抗体 | Life Technologies | S911 |
| 抗PAX5抗体 | Cell Signaling | 8970 | |
| 抗HA抗体 | Cell Signaling | 3724 | |
| 抗NR2C2抗体 | Cell Signaling | 31646 | |
| 抗NCOR1抗体 | Cell Signaling | 5948 | |
| NRIP1-HRP抗体 | Santa Cruz Biotechnology | sc-518071 | |
| NR3C1抗体 | Cell Signaling | 12041 | |
| 结构分析 | 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 |