在慢性淋巴细胞性白血病患者中使用甲基CpG结合域捕获法进行综合DNA甲基化分析

摘要

该工作描述了优化的甲基CpG结合结构域（MBD）测序方案和用于鉴定慢性淋巴细胞性白血病（CLL）患者中差异甲基化CpG富集区的计算管道。

摘要

长期非编码RNA（lncRNA）在癌症中的作用已经走在前列，因为在癌症发展和进展期间越来越需要了解其机械功能。尽管如此，在癌症中lncRNA和重复序列的全球表观遗传调控尚未得到很好的研究，特别是在慢性淋巴细胞性白血病（CLL）中。本研究着重于一种独特的方法:使用甲基结合结构域（MBD）蛋白进行基于免疫沉淀的双链甲基化DNA片段捕获，随后进行下一代测序（MBD-seq）。本研究使用属于两个预后亚组（5个IGVH突变样本+ 5个IGVH未突变样本）的CLL患者样本。与正常健康对照相比，分析显示5,800个高甲基化和12,570个低甲基化的CLL特异性差异甲基化基因（cllDMG）。重要的是，这些结果确定了几种CLL特异性差异甲基化的lncRNA具有潜在预后价值的蛋白质编码基因。该工作概述了MBD-seq和生物信息学管道的详细协议，用于使用CLL患者样本全面分析高CpG丰富地区的全球甲基化谱。最后，使用焦磷酸测序验证蛋白质编码基因和lncRNA，这是一种高度定量的方法来分析CpG甲基化水平，以进一步证实MBD-seq方案的发现。

引言

近年来，使用下一代测序技术分析全球DNA甲基化谱已越来越受欢迎。基于全基因组的甲基化测定，包括基于微阵列和非微阵列的方法，基于以下方法开发:基因组DNA的亚硫酸氢盐转化，甲基化敏感的限制酶消化和甲基化DNA的免疫沉淀，使用甲基CpG特异性抗体。

异常DNA甲基化是白血病和淋巴瘤的标志之一，包括慢性淋巴细胞性白血病（CLL）。早些时候，包括我们在内的几个组织表征不同CLL预后亚组和正常健康B细胞对照的DNA甲基化谱，使用基因组DNA的亚硫酸氢盐转化，然后是基于微阵列的方法或全基因组测序^{1,2,3 ， 4。基因组DNA的亚硫酸氢盐转化导致未修饰的胞嘧啶脱氨酸尿嘧啶，留下基因组中修饰的甲基化胞嘧啶。一旦转化，DNA的甲基化状态可以通过PCR扩增和使用不同的定量或定性方法如基于微阵列或全基因组亚硫酸氢盐测序（WGBS）的测序来测定。尽管基于亚硫酸氢盐转化的方法具有许多优点，并且广泛用于不同的癌症类型以分析DNA甲基化水平，但与该技术相关的一些缺点是。 WGBS测序允许单碱基对分辨率较低的DNA，是分析大量样品的最佳选择。然而，这种方法不能区分基因组中5mC和5mmC水平之间的修饰。此外，基于微阵列的方法不提供完整的c基因组超量。}

在我们实验室⁷最近的一项研究中，基于免疫沉淀的方法，而不是亚硫酸氢盐转化，用于在CLL患者和正常健康对照组中在全球范围内鉴定高度富含CpG的差异甲基化区域。在甲基CpG结合域（MBD）下一代测序（MBD-seq）中，双链片段化DNA的富集取决于CpG甲基化程度。该方法可以克服亚硫酸氢盐转化方法的缺点，并且还可以以无偏差和不依赖PCR的方式提供CpG甲基化的全基因组覆盖。此外，与基于亚硫酸氢盐转化的微阵列方法不同，MBD-seq可用于分析重复元件的甲基化状态，例如长散布核元素（LINE），短散点核元素（SINE），长末端重复序列（LTRS） 等等然而，与亚硫酸氢盐转化方法相比，MBD-seq协议需要相当大量的输入DNA。此外，测序的质量读数和数据取决于所用抗体的特异性，亲和性和质量。

目前的研究解释了详细的MBD-seq方案来丰富用于下一代测序的甲基化DNA。它使用商业甲基化DNA结合浓缩试剂盒（在材料表中列出 ），以及计算流程，用于显示和解释甲基化测序数据，以鉴定与正常健康对照相比的CLL特异性高和低甲基化区域。基本上，这种方法利用人类MBD2蛋白与甲基化CpG相互作用的MBD的能力来提取富含甲基化CpG的DNA，其次是甲基化DNA的高通量测序。

Access restricted. Please log in or start a trial to view this content.

研究方案

收集CLL样本的道德准许是2007-05-21，具有以下注册号:EPN Gbg dnr 239/07。所有CLL患者根据最近修订的标准⁸诊断，并在诊断时收集样品。在获得书面同意后，研究中的患者被纳入瑞典西部不同血液学部门。本研究仅选择了白血病细胞百分率≥70％的CLL外周血单个核细胞（PBMC）样本。

准备工作

1. 根据制造商的说明书，使用商业分离试剂盒（参见材料表 ），从CLL和正常的健康外周血样品中分离PBMC进行DNA提取。
2. 高压灭菌1.5 mL管。解冻甲基化DNA结合试剂盒中的所有试剂（参见材料表 ）/ LI>
3. 使用无DNA酶的水从试剂盒提供的5x储备洗涤缓冲液中制备10 mL 1x珠洗涤缓冲液以洗涤珠粒并稀释试剂盒提供的MBD蛋白质。
4. 事先获得或准备以下物品:3M乙酸钠，无水乙醇和70％乙醇用于DNA沉淀（ 材料表 ）。

2.基因组DNA提取和超声处理

1. 使用市售的DNA提取试剂盒（ 材料表 ）根据制造商的方案从患者和正常PBMC样品中分离基因组DNA。使用分光光度计在260nm定量基因组DNA。
   注:DNA提取柱容量为5-6万个细胞，最大;因此，对于超过500万个细胞的样品使用多于一个色谱柱是很重要的。将DNA洗脱两次，共计100μL10 mM Tris EDTA（TE）缓冲液。用于甲基的MBD-生物素蛋白用于MBD测序的矿工试剂盒不会与单链DNA结合。因此，保持洗脱的DNA冷冻或在4℃保持其双链性质是非常重要的。
2. 使用TE缓冲液（pH8）将每个样品的5μg基因组DNA稀释至总共200μL，最后浓度为25ng /μL。
3. 使用超声波仪（ 材料表 ）在特殊设计的管中进行超声处理，总共30个循环（每个循环30秒和30秒;每5个循环后，短暂旋转管以将样品收集到底部）。
   注意:这将产生150 bp和300 bp之间的片段DNA，这对于该方案是最佳的。
4. 检查所有样品的超声处理范围，然后通过使用DNA电泳成像设备，在市售的预制2％琼脂糖凝胶上运行1μL片段DNA样品和DNA大小梯度，进行下一步的MBD测序。可视化t他使用标准的紫外线透射仪DNA。

3.与MBD-生物素蛋白结合之前的珠粒制备

1. 从试剂盒提供的储备管中重新悬浮磁性链亲和素珠，轻轻上下移液以获得均匀的悬浮液。不要涡旋珠子或让它们干燥。
2. 将50μL珠（每个5μg片段化的DNA样品）放在分开的，干净的和标记的1.5mL管中。加入50μL的1x珠洗涤缓冲液，达到100μL的最终体积。
   注意:当使用小的0.5 mL聚合酶链反应（PCR）条纹用于洗涤珠子时，如试剂盒方案中所述，使用约150-200μL的每个洗涤缓冲液。然而，在1.5mL管的情况下，每管加入至少250μL至300μL洗涤缓冲液用于珠洗涤。
3. 将管置于磁性架上1分钟，以使所有的磁珠集中在管的内壁上面对磁铁。取出液体，不用珠子接触，使用200μL移液管。
4. 从磁性支架上取出管子，加入250μL的1x珠洗涤缓冲液，并用移液管轻轻混合。
5. 对所有样品重复步骤3.3和3.4至少4-5次，最后重悬于250μL的1x珠洗涤缓冲液中。把它们放在冰上。

4.将MBD-生物素蛋白与洗涤的珠结合

1. 加入35μLMBD蛋白（7μL1μgDNA样品）分离管，并使用1倍珠洗涤缓冲液将总体积达250μL。
2. 将250μL稀释的MBD蛋白加入到250μL洗涤的珠中，使其在室温下端对端旋转1小时。
3. 将珠粒和蛋白质混合1小时后，用磁性链亲和素珠清洗MBD蛋白质生物素。
   1. 将管置于磁性架上1分钟，取出t他使用移液器不接触珠子。加入250μL的1x珠洗涤缓冲液，并将管置于旋转混合器中5分钟。
4. 重复步骤4.3.1两次，最后将洗涤的MBD-生物素珠重悬于200μL的1倍珠洗涤缓冲液中，使珠准备好进行甲基化DNA捕获。
   注意:以这种方式洗2-3次，除去所有背景未结合的MBD蛋白质珠，并改善MBD蛋白偶联珠与片段化基因组DNA的有效结合。

5.将MBD-生物素珠与片段基因组DNA结合

1. 在干净的1.5 mL无DNA酶的管中，加入100μL5倍珠洗涤缓冲液和180μL片段化基因组DNA（步骤2.3）。使用无DNA酶的水，使最终体积达到500μL。
2. 向剩余的20μL片段化的基因组DNA中加入380μL无DNA酶的水并冷冻。使用这些样品作为输入DNA随后与最终洗脱的甲基化DNA样品一起控制和沉淀，（见步骤7.1）。
3. 将含有洗涤的MBD-生物素珠（步骤4.4）的管置于磁性架上1分钟，并除去液体，而不干扰珠粒。加入500μL在珠洗涤缓冲液中稀释的片段化基因组DNA。
4. 用石蜡膜密封所有的管子，并在4℃下以8-10rpm在端对端的旋转架上放置过夜。
   注意:DNA和生物素 - 珠粒结合反应可以在室温下进行1小时，但在4℃下保持过夜可以提高最终甲基化DNA的回收率。

6.去除未结合的DNA并从珠中洗脱甲基化的DNA

1. DNA和MBD-珠粒结合反应后，将管置于磁性架上1分钟，将所有的珠粒集中在管内壁上。
2. 用移液器取出上清液而不接触并将此未结合的DNA样品级分保存在冰上。
3. 向珠中加入200μL的1x珠洗涤缓冲液，并将管在室温下放置在旋转架上3分钟。将管放在磁性支架上并取出液体。重复洗涤两次以除去残留的未结合的DNA。
4. 最终洗涤后，加入200μL高盐洗脱缓冲液（2,000 mM NaCl），在试剂盒中洗脱DNA。
5. 将管在室温下放置在旋转架上15分钟。将它们放在磁性架上1分钟，并使用移液管小心地将上清液转移到新的干净的1.5 mL管中。
6. 加入200μL高盐洗脱缓冲液，并在室温下旋转管重复洗脱15分钟。将第二个洗脱液加入含有200μL第一个洗脱液的同一根管中。
   注意:最终的400μL总洗脱DNA现在可以进行乙醇沉淀，以提取纯化的适合下一代测序的DNA。

乙醇沉淀和甲基化DNA的富集

1. 加入1μL糖原（20μg/μL;包含在试剂盒中）; 40μL3M乙酸钠，pH 5.2;和800μL冰冷的无水乙醇至400μL洗脱的DNA，以及步骤5.2中制备的400μL的输入DNA样品。
2. 通过涡旋混合管并在-80℃下孵育过夜。
3. 以12,000 xg（最大速度）和4°C离心管。仔细弃去上清液，不会干扰沉淀，加入500μL70％乙醇，并涡旋管。
4. 在4℃下以最大速度再次离心15分钟，小心地用移液管除去上清液。在室温下以最大速度离心管1分钟，并用移液管吸头彻底清除残留的乙醇。
5. 在室温下将颗粒空气干燥5分钟rature。向DNA沉淀中加入10μL无DNA酶的水。进行甲基化DNA定量（步骤7.6），MBD测序（步骤7.7）和分析（第8和9节）。
6. 根据制造商的说明书，使用市售的荧光定量试剂盒（参见材料表 ）量化最终回收的甲基化DNA样品。
7. 注意:使用该方案，从每个样品中回收30-50ng的最终DNA。 DNA样本准备在干冰上发送用于下游文库构建和高通量MBD测序
8. 使用商业平台（ 材料表 ）进行DNA文库构建和高通量MBD测序，如参考文献⁹所述。
   注意:对于最终甲基化DNA的初始质量控制，使用所有样品的50 bp双末端测序进行文库制备。处理原始数据进行后序测序质量控制，a在步骤8.1中阐述，并且使用生物信息学方法和统计方法进一步处理它，如下所述（步骤8.2-9.6）。

8.生物信息学分析方法1:鉴定CLL相关差异甲基化区域（cllDMRs）

1. 使用可用的质量控制工具（如Trimmomatic ¹⁰或Cutadapt）清洁获得的49-bp读取（FASTQ格式）的适配器。使用FastQC工具包交叉检查修剪读取的质量:
   java -jar trimmomatic.jar SE SAMPLE_uncleaned.fastq SAMPLE.fastq ILLUMINACLIP:adapters.fasta:2:30:10
2. 将清洁的FASTQ文件与短读基因组对齐器Bowtie对齐参考基因组¹¹ 。指定参数如下:
   1. 允许最多两个不匹配（Bowtie参数: -v 2）。
   2. 将报告限制为每次读取六个最佳对齐（Bowtie参数: -m 6）进行控制用于多映射读取。
      bowtie -v 2 -a -m 6 HG19_INDEX -S SAMPLE.fastq> SAMPLE.sam
      注意:使用SAMtools将对齐后的各个样本生成的SAM文件转换为BAM。
      samtools视图-bS -o SAMPLE.bam SAMPLE.sam
3. 使用对齐样本（BAM）的ChIP-Seq（MACS）峰值呼叫者的基于模型的分析来预测来自CLL亚组^12的富集/甲基化区域。
   注意:比较I:使用输入样品作为对照组，将正常和CLL患者样品作为治疗组。此步骤是收集所有负峰（丰富输入/背景的峰）。比较二:使用正常样本组作为对照组，CLL患者样本组作为治疗组。获得的正峰是CLL高甲基化区，负峰是正常或差异甲基化区（DMR）的CLL低甲基化区。
   macs14 -t SAMPLE_TREATMENT.bam -c SAMPLE_CONTROL.bam - 格式BAM -g hs
4. 使用BED工具去除差异甲基化区域（DMR）中的比较I背景峰¹³ 。
   bedtools subtract -a -b
5. 预测富集DMR的重复元素（SINE-Alu，LINE 等 ）的百分比。
   1. 使用fastacmd通过参考基因组HG19的染色体坐标提取DMR的FASTA序列。
      fastacmd -d HG19_genome.fa -s染色体-L开始，结束-1 50000> DMRs.fasta
   2. 使用RepeatMasker命令行工具来预测DMR中FASTA序列中存在的重复元素的百分比。
      RepeatMasker -gc -gccalc -s -species human -html DMRs.fasta
6. 使用荷马"annotatePeaks.pl"与可用的Ensembl [PMC4919035] o注释预测的DMRr代码[PMID 22955987]转录本注释（蛋白质编码和非编码转录本）。
   annotatePeaks.pl DMRs.bed -gtf
   注意:这提供了关于不同基因区域（ 即启动子，外显子，内含子，3'UTR和5'UTR）和基因间区域的峰分布的信息。与DMR相关的转录物或基因称为差异甲基化基因（DMG）。
7. 使用富集工具与更新的或当前的功能数据库来查找CLL DMGs（仅蛋白质编码基因）丰富的功能¹⁴ 。
   注意:基于功能注释的基因集聚簇（GeneSCF） ¹⁴是一种基于实时的功能丰富分析工具，使用更新的KEGG和基因本体作为参考数据库。

9.生物信息学分析方法2:识别CLL相关显着差异甲基化区域（cll sigDMR）

1. 按照方法I（第8节）中的步骤8.1-8.5，并使用基于读数计数的差异浓缩分析，利用步骤9.2 - 9.6中详细介绍的可用工具，将更多的统计信息添加到分析流程中。
2. 使用"Subread"软件包¹⁵中的"featureCounts"量化在普通和CLL患者样本中映射到单独峰值或DMR的读数。
   subread / bin / featureCounts -Q 30 -F SAF -a DMRs.SAF -o DMRs_counts.table SAMPLE_TREATMENT.bam SAMPLE_CONTROL.bam
   注意:可以引入映射质量过滤器来避免量化不良品质的读取（ 例如，-Q 30）。为此步骤，准备获得的DMR的SAF文件。有关SAF文件格式的更多信息，请使用此链接http://bioinf.wehi.edu.au/featureCounts/。
3. 使用RAW读数计数表，其中包含Normal a中各个峰的读数d CLL患者样本组在edgeR中作为输入¹⁶ 。
   注意:比较正常与CLL患者样本组以找到sigDMRs。对于差异富集分析，请按照edgeR的指南详细分步说明（https://bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf）。
4. 过滤sigDMRs使用false discovery rate（FDR）和log-fold更改由edgeR ¹⁶预测。
5. 使用可用的Ensembl或Gencode转录本注释（蛋白质编码和非编码转录本），使用Homer"annotatePeaks.pl"注释预测的sigDMR。
   annotatePeaks.pl sigDMRs.bed HG19_genome.fa -gtf -CpG
   注意:这提供了关于不同基因区域（ 即启动子，外显子，内含子，3'UTR和5'UTR）和基因间区域的峰分布的信息。与si相关的转录物或基因gDMR称为sigDMG。
6. 使用更新或当前功能数据库的浓缩工具来查找CLL sigDMG（仅蛋白编码基因）富集的功能¹⁴ 。
   注:GeneSCF是功能丰富分析的实时工具之一，使用KEGG和Gene本体作为参考数据库。

Access restricted. Please log in or start a trial to view this content.

结果

最近对CLL患者和匹配的正常健康对照进行了MBD-seq，以鉴定CLL特异性差异超低和低甲基化基因⁷ 。用于分析CLL和正常健康样品产生的数据的实验和生物信息管道如图1A和1B所示 。这些分析确定了几个CLL特异性差异甲基化区域（cllDMR），与对照样品相比，IGHV突变和IGHV未突变样品显着高/低甲基化，p值<0.00001。

Access restricted. Please log in or start a trial to view this content.

讨论

MBD-seq是一种经济有效的基于免疫沉淀的技术，可用于研究具有完整全基因组范围覆盖的甲基化模式。 MeDIP seq（甲基化DNA免疫沉淀随后测序）和MBD-seq均导致富含CpG的甲基化DNA的富集。然而，与MeDIP seq ¹⁹相比，MBD-seq对结合高度富含CpG的区域具有更强的亲和力。使用甲基结合富集试剂盒，可以分别通过高盐和低盐缓冲液洗脱DNA将DNA分馏成高CpG和低CpG富集区。在该研究中，仅进行单?...

Access restricted. Please log in or start a trial to view this content.

材料

Name	Company	Catalog Number	Comments
Dneasy Blood and tissue kit	Qaigen	69504
Lymphoprep solution	A X I S-S H I E L D	1114544
Nano drop 2000	Thermo Fischersceintific
TE buffer pH 8	Sigma aldrich	93283
Bioruptor standard sonication device	Diagenode	UCD-200
TPX bioruptor tubes 1.5 mL	Diagenode	C30010010-300
3 M Sodium acetate	Diagenode	C03030002
E-gel iBase safe imager combo kit	Thermo Fischersceintific	G6465EU
E-gel 2% Agarose gels	Thermo Fischersceintific	G441002
Methylminer Methylated DNA enrichment kit	Thermo Fischersceintific	ME10025
Labquake Tube Shaker/Rotators	Thermo Fischersceintific	415110
Dynal MPC-S	Thermo Fischersceintific	A13346
Vortex mixer	VWR	12620-848
Absolute Ethanol	Any company
70% Ethanool	Any company
DNAse free water	Milli Q
DNA precipitant (3M sodium acetate)	Diagenode	C03030002
Safe seal 1.5 mL eppendorf tubes	Eppendorf	4036-3204
Qubit dsDNA HS Assay Kit	Thermo Fischersceintific	Q32851
Qubit 0.5 mL tubes	Thermo Fischersceintific	Q32856
Qubit	Thermo Fischersceintific	Q32866
Illumina Hiseq2000 Platform	Illumina
Water  Bath	Grant
Heat block	grant
Tube rotater	Labquake