WGBS数据分析和PMD鉴定
分析食管癌和正常组织的全基因组甲基化模式,并开发新的PMD识别算法。
对45个食管样本进行WGBS测序,使用开发的MMSeekR方法识别PMDs,并与MethylSeekR和MethPipe进行比较。
Comprehensive analyses of partially methylated domains and differentially methylated regions in esophageal cancer reveal both cell-type- and cancer-specific epigenetic regulation
包含WGBS数据分析、PMD/DMR鉴定、ChIP-seq验证、功能实验(shRNA敲低、ChIP-qPCR)以及多组学数据整合分析等多个独立证据层级。
食管癌肿瘤及非恶性组织样本(ESCC、EAC/NGEJ、NESQ) 食管癌细胞系(ESO26、OE19、TE5、KYSE70) 头颈部鳞状细胞癌细胞系(Cal27、Det562、FaDu)
WGBS数据的样本量和覆盖度(45个食管样本,≥3或≥5覆盖度) PMD calling方法(MMSeekR vs. MethylSeekR、MethPipe)及参数(窗口大小、模型选择) DMR鉴定参数(q值<0.05,delta methylation>0.2) H3K36me2 ChIP-seq在EAC和ESCC细胞系中的信号 HNF4A敲低后FOSL1结合变化的ChIP-qPCR验证
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
分析食管癌和正常组织的全基因组甲基化模式,并开发新的PMD识别算法。
对45个食管样本进行WGBS测序,使用开发的MMSeekR方法识别PMDs,并与MethylSeekR和MethPipe进行比较。
鉴定食管癌两种亚型特有的PMDs,并分析其功能相关性。
根据频率将PMDs分为共享、EAC特异性、ESCC特异性和共享HMDs,分析其甲基化水平、染色质区室、体细胞突变率及基因表达关联。
验证H3K36me2修饰与PMDs分布的关系。
在EAC和ESCC细胞系中进行H3K36me2 ChIP-seq,并分析其在PMD类别中的信号强度。
在去除PMDs效应后,鉴定亚型特异性和肿瘤特异性的差异甲基化区域,并分析其调控功能。
使用Dmrseq方法鉴定DMRs,进行功能富集分析、基序分析、ATAC-seq和H3K27ac ChIP-seq数据整合,并验证转录因子结合。
验证HNF4A和FOSL1在肿瘤特异性DMRs中的协同结合。
通过shRNA敲低HNF4A,使用ChIP-qPCR检测FOSL1在选定DMR区域的结合变化。
验证PMDs和DMRs的细胞类型特异性在不同组织和癌症类型中的普遍性。
分析TCGA泛癌甲基化和ATAC-seq数据,以及正常组织的单细胞ATAC-seq数据,评估PMDs和DMRs的细胞类型特异性。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | RPMI-1640培养基 | Gibco | -- |
| 胎牛血清 | Omega Scientific | -- | |
| 青霉素-链霉素 | Thermo Scientific | -- | |
| WGBS | EZ DNA甲基化金试剂盒 | Zymo Research | -- |
| KAPA HiFi HotStart Uracil + Ready Mix | Kapa Biosystems | -- | |
| ChIP-seq | H3K36me2抗体 | Cell Signaling Technology | 2901S |
| FRA1抗体 | Cell Signaling Technology | 5281 | |
| Protein G Dynabeads | Life Technologies | -- | |
| Covaris超声仪 | -- | -- | |
| WGBS测序 | Illumina HiSeq 2500平台 | Illumina | -- |
| ChIP-seq测序 | Illumina HiSeq平台 | Illumina | -- |
| 基因敲低 | pLKO-puro载体 | Addgene | 8453 |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 细胞培养 | 细胞系来源(ESO26、OE19、TE5、KYSE70)及培养条件(培养基、血清、抗生素、CO2浓度) 阅读提示:Methods: Cell culture |
| WGBS | DNA起始量(3 μg)、lambda DNA用量(26 ng)、测序平台(HiSeq 2500) 阅读提示:Methods: Whole-genome bisulfite sequencing (WGBS) |
| ChIP-seq | 细胞数量(10 million)、交联条件(1% PFA 10 min)、抗体信息(H3K36me2、FRA1) 阅读提示:Methods: Chromatin immunoprecipitation sequencing (ChIP-Seq) and ChIP-qPCR |
| HNF4A敲低 | shRNA序列、载体(pLKO-puro)、敲低效率验证 阅读提示:Methods: Cell culture (shRNA knockdown part) |
| WGBS数据分析 | 比对工具(Biscuit)、甲基化提取条件(覆盖度≥5)、PMD calling参数(MMSeekR) 阅读提示:Methods: DNA methylation data analysis, Application of MMSeekR to WGBS data |
| DMR分析 | DMR鉴定参数(q<0.05, delta>0.2)、PMD掩蔽条件(2/3样本) 阅读提示:Methods: Identification and characterization of DMRs |