食管癌部分甲基化结构域和差异甲基化区域的综合分析揭示了细胞类型特异性和癌症特异性的表观遗传调控

Comprehensive analyses of partially methylated domains and differentially methylated regions in esophageal cancer reveal both cell-type- and cancer-specific epigenetic regulation

作者信息Yueyuan Zheng, Benjamin Ziman, Allen S Ho, Uttam K Sinha, Li-Yan Xu, En-Min Li, H Phillip Koeffler, Benjamin P Berman, De-Chen Lin
PMID37620896
发布时间2023-08-24
DOI10.1186/s13059-023-03035-3

实验完整度

包含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验证

摘要

背景:作为最常见的恶性肿瘤之一,食管癌有两种亚型,鳞状细胞癌和腺癌,起源于不同的细胞。区分细胞类型特异性的分子特征与癌症特异性的特征是具有挑战性的。结果:我们分析了来自两种亚型的45个食管肿瘤和非恶性样本的全基因组亚硫酸氢盐测序数据。我们开发了一种新的序列感知方法来识别大的部分甲基化结构域(PMDs),揭示了肿瘤样本中PMDs在甲基化水平和基因组分布上的深刻异质性。我们识别出与抑制性转录、染色质B区室和高体细胞突变率相关的亚型特异性PMDs。虽然这些PMDs的基因组位置在正常细胞中是预先建立的,但肿瘤中的丢失程度显著更高。我们发现细胞类型特异性的H3K36me2沉积可能是PMDs基因组分布的基础。在较小的基因组尺度上,我们为每种亚型识别出细胞类型特异性和癌症特异性的差异甲基化区域(DMRs)。通过在这些DMRs内进行结合基序分析,我们发现细胞类型特异性的转录因子HNF4A在正常细胞中维持其生成的结合位点,同时在食管腺癌中与新的伙伴(如FOSL1)合作建立新的结合位点。最后,利用泛组织单细胞和泛癌表观基因组数据集,我们证明了在食管癌中鉴定出的相当一部分细胞类型特异性PMDs和DMRs实际上是起源于相关细胞类型的其他癌症中共同出现的标记。结论:这些发现加深了我们对正常和恶性状态下各种基因组尺度DNA甲基化动态的理解,为细胞类型特异性和癌症特异性的表观遗传调控提供了新的机制见解。

实验结论

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

研究问题
食管癌两种亚型中,细胞类型特异性和癌症特异性的表观遗传调控(如PMDs和DMRs)如何区分,以及其潜在机制是什么?
核心机制
细胞类型特异性的H3K36me2沉积可能通过招募DNMT3A维持DNA甲基化,从而决定PMDs的基因组分布;在肿瘤特异性DMRs中,HNF4A与FOSL1等AP-1因子合作建立新的结合位点,调控肿瘤相关基因表达。
主要证据
基于WGBS数据鉴定PMDs和DMRs,结合ChIP-seq、ATAC-seq、单细胞数据及功能实验(HNF4A敲低导致FOSL1结合减少),验证了细胞类型特异性和癌症特异性的表观遗传特征。
研究意义
为理解食管癌中DNA甲基化动态提供了新见解,并提示PMDs和DMRs的甲基化水平可作为潜在生物标志物用于癌症诊断。

研究路径

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

1

WGBS数据分析和PMD鉴定

分析食管癌和正常组织的全基因组甲基化模式,并开发新的PMD识别算法。

对45个食管样本进行WGBS测序,使用开发的MMSeekR方法识别PMDs,并与MethylSeekR和MethPipe进行比较。

2

亚型特异性PMDs的表征

鉴定食管癌两种亚型特有的PMDs,并分析其功能相关性。

根据频率将PMDs分为共享、EAC特异性、ESCC特异性和共享HMDs,分析其甲基化水平、染色质区室、体细胞突变率及基因表达关联。

3

H3K36me2 ChIP-seq分析

验证H3K36me2修饰与PMDs分布的关系。

在EAC和ESCC细胞系中进行H3K36me2 ChIP-seq,并分析其在PMD类别中的信号强度。

4

DMR鉴定和功能分析

在去除PMDs效应后,鉴定亚型特异性和肿瘤特异性的差异甲基化区域,并分析其调控功能。

使用Dmrseq方法鉴定DMRs,进行功能富集分析、基序分析、ATAC-seq和H3K27ac ChIP-seq数据整合,并验证转录因子结合。

5

肿瘤特异性DMRs的实验验证

验证HNF4A和FOSL1在肿瘤特异性DMRs中的协同结合。

通过shRNA敲低HNF4A,使用ChIP-qPCR检测FOSL1在选定DMR区域的结合变化。

6

泛癌和单细胞分析

验证PMDs和DMRs的细胞类型特异性在不同组织和癌症类型中的普遍性。

分析TCGA泛癌甲基化和ATAC-seq数据,以及正常组织的单细胞ATAC-seq数据,评估PMDs和DMRs的细胞类型特异性。

研究方法

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

产品清单

实验环节名称品牌货号
RPMI-1640培养基Gibco--
胎牛血清Omega Scientific--
青霉素-链霉素Thermo Scientific--
EZ DNA甲基化金试剂盒Zymo Research--
KAPA HiFi HotStart Uracil + Ready MixKapa Biosystems--
H3K36me2抗体Cell Signaling Technology2901S
FRA1抗体Cell Signaling Technology5281
Protein G DynabeadsLife Technologies--
Covaris超声仪----
Illumina HiSeq 2500平台Illumina--
Illumina HiSeq平台Illumina--
pLKO-puro载体Addgene8453

关键环节

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

环节核对要点
细胞培养
细胞系来源(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