数据收集与配对筛选
构建用于分析转录因子结合差异的ChIP-seq配对数据集。
从GTRD数据库下载ChIP-seq数据,使用UniBind p值筛选高质量实验,通过Jaccard距离和层次聚类选择具有显著结合差异的配对。
TFscope: systematic analysis of the sequence features involved in the binding preferences of transcription factors
方法基于公开ChIP-seq数据,进行了大规模配对分析,包括细胞类型、处理和旁系同源TF比较,并验证了模型准确性,但属计算方法研究,缺乏湿实验功能验证。
人类细胞系(如HepG2、K562、MCF7等)的ChIP-seq数据 小鼠细胞系SKH1和U937的ChIP-seq数据 旁系同源TF FOSL1/FOSL2、AR/NR3C1的ChIP-seq数据
272对细胞类型ChIP-seq数据的筛选标准(UniBind p值) 79对处理ChIP-seq数据的筛选标准(同一研究、排除时程) JASPAR2020数据库基序 FIMO工具和峰值重定位设置 模型训练与测试集划分比例(70%训练,30%测试)
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
构建用于分析转录因子结合差异的ChIP-seq配对数据集。
从GTRD数据库下载ChIP-seq数据,使用UniBind p值筛选高质量实验,通过Jaccard距离和层次聚类选择具有显著结合差异的配对。
将ChIP-seq峰对齐到最可能的转录因子结合位点,以便提取特征。
使用JASPAR PWM通过FIMO扫描峰附近序列,选择最佳结合位点作为锚点,提取以锚点为中心的1kb序列,去除无基序序列,并平衡两类序列数量。
从序列中提取三类可能的判别特征:核心基序、核苷酸环境、辅因子。
使用TFscope-DM学习判别性PWM,TFscope-NE提取差异性的k-mer区域频率,TFscope-CF识别差异性的辅因子PWM和区域。
整合三类特征,训练分类模型以区分两个ChIP-seq实验的结合位点,并评估特征重要性。
将DM、NE、CF变量合并,使用LASSO正则化的逻辑回归模型训练,通过AUROC评估预测能力,并通过移除不同特征计算AUROC损失来衡量特征重要性。
验证TFscope在不同情境下的适用性和准确性,并总结结合差异的规律。
对272对细胞类型、79对处理和若干旁系同源TF进行全流程分析,比较不同模型(如仅用核心基序,不含某些特征)的效果,并通过聚类总结特征使用模式。
按研究目的归类文中使用的方法,便于定位所需技术。
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 数据选择 | ChIP-seq数据的来源(GTRD)、UniBind p值阈值(<1%)、Jaccard距离和层次聚类参数 阅读提示:Material and methods: Selection of 272 ChIP-seq pairs |
| 序列提取与比对 | JASPAR PWM版本选择、FIMO参数、1kb序列提取、平衡策略 阅读提示:Material and methods: Sequence extraction and alignment |
| 特征提取 | TFscope-DM的序列长度、TFscope-NE的k-mer长度和bin数、TFscope-CF的bin数和PWM库 阅读提示:Material and methods: TFscope-DM, TFscope-NE, TFscope-CF |
| 模型训练与评估 | 训练/测试集比例(70/30)、LASSO惩罚参数、模型变量选择标准(10个参数) 阅读提示:Material and methods: Measure of variable importance, Model learning |