全基因组DNA甲基化分析表明胰腺上皮内瘤变病变源于腺泡细胞并在表观遗传上为癌变做好准备

Comprehensive DNA Methylation Analysis Indicates That Pancreatic Intraepithelial Neoplasia Lesions Are Acinar-Derived and Epigenetically Primed for Carcinogenesis

作者信息Emily K W Lo, Brian M Mears, H Carlo Maurer, Adrian Idrizi, Kasper D Hansen, Elizabeth D Thompson, Ralph H Hruban, Kenneth P Olive, Andrew P Feinberg
PMID36989344
发布时间2023-06
DOI10.1158/0008-5472.CAN-22-4052

实验完整度

研究包含多层级实验证据:激光显微切割获取四种纯化胰腺细胞进行全基因组亚硫酸氢盐测序,以bsseq和informME双重互补方法鉴定DMR,并通过GSEA、EpiDISH去卷积、体细胞突变验证等验证,体外和患者样本层面均进行了功能验证。

主要模型

人胰腺腺泡细胞(正常腺泡) 人胰腺导管细胞(非肿瘤导管) 人胰腺上皮内瘤变(PanIN)病变 人胰腺导管腺癌(PDAC)病变

重点核对

样本类型:四种细胞类型(腺泡、导管、PanIN、PDAC)来自18例患者的手术切除标本 DNA起始量:20–50 ng用于WGBS 测序平台:HiSeq 4000(2×151 bp)或NovaSeq 6000(2×150 bp) 亚硫酸氢盐转化率:≥98.8% DMR识别标准:小DMR(<5 kb)采用t统计量0.01/0.99分位数、CpG数≥10、平均甲基化差异≥25%;块DMR(≥5 kb)采用F统计置换检验,FWER≤0.01

摘要

胰腺导管腺癌(PDAC)被认为是由一系列体细胞突变的积累引起的,并且通常与胰腺上皮内瘤变(PanIN)病变相关。然而,关于人类PanIN和PDAC的起源细胞是腺泡细胞还是导管细胞仍存在争议。由于细胞类型身份在表观遗传上得以维持,胰腺肿瘤发生过程中的DNA甲基化变化可以为研究这个问题提供有力的视角。在这里,我们对18名患者的手术切除标本进行了激光捕获显微切割,以高纯度分离四种相关细胞类型的DNA用于全基因组亚硫酸氢盐测序:腺泡细胞、非肿瘤导管、PanIN病变和PDAC病变。使用两种互补的分析方法鉴定了差异甲基化区域(DMR):bsseq,它识别任何DMR,但对大型块状DMR特别有用;以及informME,它分析整个基因组的潜在能量景观,特别适用于识别差异甲基化熵。全局甲基化谱和块状DMR均明确表明PanIN起源于腺泡细胞。在基因水平上,PanIN病变表现出介于腺泡-导管之间的中间表型,类似于腺泡-导管化生。在97.6%的PanIN特异性DMR中,PanIN病变的甲基化水平介于正常和PDAC之间,从信息论的角度表明PanIN病变在表观遗传上已准备好进展为PDAC。因此,表观基因组分析补充了组织病理学,定义了向PDAC的分子进展。PanIN和PDAC病变之间的共同表观遗传谱系可能通过靶向异常甲基化的进展相关基因提供预防机会。

实验结论

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

研究问题
人类胰腺上皮内瘤变(PanIN)和胰腺导管腺癌(PDAC)的起源细胞是腺泡细胞还是导管细胞?
核心机制
PanIN病变在表观遗传上处于腺泡-导管化生(ADM)的中间状态,并表现出朝向PDAC的表观遗传预激活(epigenetic priming),其甲基化水平介于正常和PDAC之间。
主要证据
通过对18例患者样本的激光显微切割获取四种细胞类型进行WGBS,全局甲基化分析显示PanIN与腺泡聚为一组;在97.6%的PanIN特异性DMR中,PanIN甲基化水平为中间值;基因水平GSEA显示PanIN为腺泡-导管中间表型。
研究意义
表观基因组分析补充了组织病理学,有助于定义向PDAC的分子进展,并可能通过靶向异常甲基化的进展相关基因提供化学预防机会。

研究路径

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

1

样本获取与纯化

获取高纯度的四种胰腺细胞类型(腺泡、导管、PanIN、PDAC)用于全基因组甲基化分析。

对18例患者的手术切除标本进行激光捕获显微切割(LCM),提取DNA,并通过体细胞突变检测和EpiDISH去卷积验证纯度。

2

全基因组亚硫酸氢盐测序

获得全基因组DNA甲基化图谱,以比较四种细胞类型。

对提取的DNA进行WGBS文库制备和测序,使用NEBNext或Accel-NGS试剂盒,在HiSeq或NovaSeq平台测序。

3

全局甲基化分析

评估四种细胞类型的整体甲基化差异,并初步判断PanIN的起源。

进行主成分分析(PCA)、层次聚类、平均甲基化水平(MML)和归一化甲基化熵(NME)密度图分析,以及甲基化块的分析。

4

差异甲基化区域(DMR)鉴定

识别基因水平上细胞类型间差异甲基化的区域,分析PanIN与腺泡/导管的重叠。

使用bsseq进行小DMR(<5 kb)的成对比较,使用informME计算Jensen-Shannon距离识别差异甲基化区域。

5

基因集富集分析

评估PanIN病变的基因表型状态,验证其是否为腺泡-导管中间状态。

对DMR相关基因进行Enrichr过度表示分析和基于排名的GSEA(fgsea),使用MSigDB和自定腺泡/导管/ADM基因集。

6

PanIN特异性DMR和进展DMR分析

区分PanIN是前体病变还是终末状态,并识别与癌变相关的DMR。

通过置换检验识别PanIN特异性DMR,并筛选中间甲基化模式的DMR;同时识别neoplastic vs normal的进展DMR并进行富集分析。

研究方法

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

产品清单

实验环节名称品牌货号
PEN膜玻片Applied BiosystemsLCM0522
QIAamp DNA微量试剂盒Qiagen56304
Qubit双链DNA高灵敏度检测ThermoFisher ScientificQ32851
Lambda DNAPromegaD1521
NEBNext Ultra DNA文库制备试剂盒New England BioLabsE7370L
Accel-NGS甲基化测序DNA文库制备试剂盒Swift Biosciences30024
Covaris S2聚焦超声破碎仪Covaris--
NEBNext Ultra末端修复/dA加尾模块New England BioLabs7442L
NEBNext Illumina多重寡核苷酸New England BioLabsE7535L
NEBNext Ultra连接模块New England BioLabs7445L
EZ DNA甲基化快速试剂盒Zymo ResearchD5030
Kapa HiFi Uracil+ PCR系统Kapa BiosystemsKK2801
AMPure XP磁珠Beckman CoulterA63881
EZ DNA甲基化金试剂盒Zymo ResearchD5005
独特双索引引物Swift BioscienceX9096-PLATE
Agilent高灵敏度DNA试剂盒Agilent Technologies5067-4626
KAPA文库定量试剂盒KAPA BiosystemsKK4824
Illumina HiSeq 4000测序仪Illumina--
Illumina NovaSeq 6000测序仪Illumina--

关键环节

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

环节核对要点
样本获取与纯化
样本来源:18例患者;细胞类型:腺泡、导管、PanIN、PDAC;激光显微切割的切片厚度(5-9 μm)
阅读提示:Materials and Methods: Sample extraction and processing
DNA提取与定量
DNA起始量:20-50 ng;QIAamp DNA Micro Kit提取
阅读提示:Materials and Methods: Sample extraction and processing
全基因组亚硫酸氢盐测序
文库制备试剂盒(NEBNext或Accel-NGS);测序平台(HiSeq 4000或NovaSeq 6000);读长(2×151或2×150);亚硫酸氢盐转化率≥98.8%
阅读提示:Materials and Methods: Whole-genome bisulfite sequencing (WGBS), alignment, and quality control
DMR鉴定
小DMR阈值:t统计量0.01/0.99、CpG数≥10、甲基化差异≥25%;块DMR阈值:ns=500、h=2000、FWER≤0.01
阅读提示:Materials and Methods: Differential methylation analysis
基因集富集分析
基因集:MSigDB C2、C5及Schlesinger等(2020)来源的自定基因集;fgsea参数:minSize=10、maxSize=1000
阅读提示:Materials and Methods: Gene set enrichment analysis (GSEA)