线粒体 DNA 甲基化的精确检测方法

摘要

在这里, 我们提出了一个协议, 允许准确定量的线粒体 DNA (线粒体) 甲基化。在本协议中, 我们描述了一种酶消化 DNA 与BamHI结合的生物信息学分析管道, 可用于避免高估线粒体甲基化水平引起的线粒体二级结构。

摘要

利用亚硫酸氢钠的特性, 将非甲基化胞嘧啶转化为尿, 在单链 dna 环境下, 可以实现 dna 甲基化的量化。亚硫酸氢的测序可以是靶向 (使用 PCR) 或在整个基因组进行, 并提供绝对定量的胞嘧啶甲基化的单基分辨率。考虑到细胞核和线粒体 DNA 的独特性质, 特别是在二级结构中, 应该对亚硫酸氢嘧啶测序方法在线粒体中的 adaptions 进行研究。线粒体的二级和三级结构确实可以导致亚硫酸氢的测序工件导致误报, 因为不完全变性, 亚硫酸氢对单链 DNA 的接触较差。在这里, 我们描述了一个协议, 使用酶消化 DNA 与BamHI结合生物信息学分析管道, 以允许准确定量的胞嘧啶甲基化水平在线粒体。此外, 我们还提供了设计特定于线粒体的亚硫酸盐测序底漆的准则, 以避免将不理想的核线粒体片段 (NUMTs) 插入核基因组。

引言

线粒体基因组是一个圆形, 双链结构约16.5 公斤的基础 (kb) 长, 构成一个沉重的和轻的链。线粒体基因组存在于每个细胞内的多个拷贝中, 母系遗传, 并编码呼吸链复合体的基本成分¹。与细菌基因组相似, 与核基因组不同, 线粒体基因组组织在许多次和三级结构中, 例如在连续和超螺旋结构中 2, 这可以在排序过程中使访问变得困难.实验³。

在细胞核中, DNA 甲基化是一种广泛研究的表观遗传标记, 在许多过程中起着重要作用, 特别是在基因表达调控中。在哺乳动物基因组中, DNA 甲基化主要发生在 deoxycytidines 嘧啶环的5位位置上, 主要是在 CG dinucleotides (或 CpG) 上。胞嘧啶甲基化在体细胞基因组中的70% 都被发现, 占总 DNA 基数的 1%, 为⁴。DNA methylations 也被描述在非 CpG 的情况下, 如 CpA, CpT 和 CpC, 并存在于各种数量的核 DNA, 其价值高达25% 的甲基化 cytosines 在胚胎干细胞^5,6,7。

虽然核基因组的胞嘧啶甲基化被广泛接受, 但线粒体 DNA (线粒体) 甲基的存在仍然存在争议。首次研究线粒体甲基化是在培养细胞中进行的, 在那里线粒体甲基化是容易检测到的, 虽然在较低的水平与核 DNA⁸。在人类和小鼠细胞中, 线粒体甲基化也被检测到低水平 (2-5%)。使用依赖于5甲基胞嘧啶捕获 (如甲基化 DNA 免疫沉淀 (MeDIP) 和定量 PCR 的检测方法, 在各种鼠标和人和细胞^线(^9,10 ,) 中也发现线粒体甲基化.^11,12。在 ELISA 检测或质谱中使用 5-甲基胞嘧啶抗体, 从纯化的线粒体分数中检测出大量的 DNA 甲基化^水平13, 14, 15,¹⁶. 然而, 上述研究中的大多数化验方法都使用的技术不是为在单一碱基分辨率下的 DNA 甲基化提供绝对量化的。

定量和 resolutive dna 甲基化分析可以通过一种名为 "亚硫酸氢" 的技术来实现, 利用亚硫酸氢钠的特性将非甲基化胞嘧啶转化为尿的单链 DNA^上下文17.利用亚硫酸氢的测序, 一系列研究发现了胞嘧啶甲基化在不同层次的存在。甲基化线粒体在 D 环区域, 12S 或16S 区域在人的^{18,19,20,21,22,23}和老鼠²⁴中容易地被发现。然而, 组织和细胞具有令人着迷的变异, 在研究中总 cytosines 的1-20%。

与这些众多的研究相比, 只有少数研究, 包括从我们的小组, 有争议的存在线粒体甲基化^3,25,26,27或质疑的生物学相关性非常低的线粒体水平 (低于 2%)²⁸。最近, 我们报告了一个潜在的亚硫酸氢序列伪影在整个线粒体亚硫酸氢化序列³的观察。我们提供的证据表明, 线粒体 DNA 的二级结构可能导致亚硫酸氢亚硫酸盐测序中的误报, 从而高估甲基化水平。我们这里提供了一个协议, 以防止亚硫酸氢的伪影转换的线粒体。该协议使用简单的酶消化 DNA, 以扰乱线粒体二级结构, 并允许完全进入亚硫酸氢化后, 亚硫酸盐排序协议。此外, 我们还提供了一个伴随的生物信息学管道, 用于分析亚硫酸氢的测序。

研究方案

1. 限制酶治疗

1. 进行线性化线粒体 dna 通过治疗总人类 DNA 与限制酶 BamHI, 那削减在位置14258在人类线粒体脱氧核糖核酸。
   注: 对于小鼠 DNA, 使用限制酶 BglII 在相同条件下, 如下所示。
2. 对于要评估线粒体甲基化的每个样本, 准备一个0.2 毫升反应管, 并添加以下 mix:3 微克基因组 DNA (由荧光量化), 15 ul 缓冲 3, 3 ul BamHI 和水高达 150 ul
3. 在热循环仪中放置管4小时, 在37摄氏度。

2. 亚硫酸氢转化

1. 使用亚硫酸氢钠转换 BamHI 处理的 DNA, 如其他地方所述²⁹。
2. 使用 200-500 ng BamHI 处理的 DNA 为每个转换反应在总容量 20 ul。DNA 回收率为50-70%。
3. 洗脱 10 ul 洗脱缓冲液中的亚硫酸氢转化 DNA。
4. 用荧光法定量测定单链 DNA。
5. 可选: 在继续执行该协议之前, 将亚硫酸氢化转化后的 DNA 储存在-20 摄氏度。长期贮存, 将亚硫酸氢转化后的 DNA 储存在-80 摄氏度

3. 亚硫酸氢化测序底漆的设计

1. 用于亚硫酸氢的设计底漆使用在线工具 Methprimer³⁰和折半查找 31, 32, 以此进行排序.
   注: Methprimer 和折半查找允许搜索 cytosines 转化为 thymines 的基因组序列上的底漆结合序列, 类似于亚硫酸氢后处理。
2. 为优化和无偏放大, 避免在引物中的 CpG dinucleotides 和设计的扩增子大小之间的 100-300 bp。
3. 确保设计的底漆是特定的线粒体 DNA 使用折半查找的功能底漆搜索。插入已设计的亚硫酸氢转化底漆, 勾选亚硫酸氢框, 并包括参考基因组和 PCR 参数。
   注: 将显示可能的 PCR 产品列表。
4. 复制前1-5 个底漆序列, 并粘贴到一个排序形式的寡核苷酸合成。
   注意: 在实验室中已验证的人和鼠线粒体 dna 底漆序列的列表显示在表 1中。
5. 根据步骤4中的说明, 采用亚硫酸氢化凝胶电泳法进行的亚硫酸盐转化 PCR 试验底漆。对所有底漆对分别进行测试和优化, 确保琼脂糖凝胶上出现正确的扩增子尺寸。
   注: 如果需要复用底漆, 添加多个引物到一个单一的 pcr 反应和可视化 pcr 产品在琼脂糖凝胶电泳或, 更 resolutive 的方法, 如聚丙烯酰胺凝胶电泳, 以区别产品的相似大小。

4. 亚硫酸氢转化 PCR 和凝胶萃取

1. 为了放大感兴趣的区域, 使用热启动 Taq 聚合酶进行 PCR。简单地, 混合100的 DNA 转换, 500 µM 感和反感底漆, 1 µL dNTP 混合 (每 dNTP 10 µM) 和聚合酶在 7.5 U/反应在总容量50µL. 运行 PCR 以以下情况: 5 分钟在95°c (六十年代94°c; 六十年代55°c *; 六十年代72°c) 35 个周期; 10分钟72°c。单独或复用底漆。
2. 可选: 当复用引物, 使用 200 ng 转换 DNA 和以下循环条件: 5 分钟95°c (六十年代94°c; 九十年代55°c *; 九十年代72°c) 35 循环; 10分钟72°c。
   注: 温度可能因底漆熔化温度而异。
3. 主题 PCR 产品的凝胶电泳2% 琼脂糖在100伏特。
4. 在柔和的紫外光下, 用手术刀提取 PCR 产物并加入管内。不要将 PCR 产品暴露在过量的紫外线下, 以避免 DNA 损伤。避免曝光时间超过三十年代。
5. 使用凝胶纯化方法纯化 PCR 产品, 如前所述³。
6. 用荧光定量测定步骤4.5 中纯化的 DNA。继续步骤5或在20°c 存储样品。

5. 亚硫酸氢亚硫酸盐测序库准备

1. 对亚硫酸氢转化 PCR 产物进行库制备。
   可选: 在这个阶段多路 PCR 产品。
2. 使用多达 100 ng 的 PCR 产品进行端修复。添加3µL 端准备酶和6.5 µL 端修复反应缓冲器 (10x) 到55.5 µL PCR 产品在0.2 毫升管。放置和孵化管在热循环仪为30分钟在20°c 之后30分钟在65°c。
3. 结扎测序适配器通过混合最终修复 DNA 与15µL 钝/TA 连接酶混合, 2.5 µL 适配器和1µL 结扎增强剂。如果使用 < 100 ng PCR 产品作为输入, 稀释适配器10次。
4. 将混合 (步骤 5.3) 放在热循环仪中, 在20摄氏度孵化15分钟。暂停热循环仪, 并添加3µl 用户酶的管。混合并放置在热循环仪的管, 孵化15分钟在37摄氏度。
5. 尺寸选择适配器结扎的 dna, 使用固体相可逆固定 (SPRI) 珠的不同比例的珠子, DNA 取决于 PCR 产品大小。
6. 添加13.5 µL H₂O, 从步骤5.4 和55µL 悬浮 SPRI 珠中的适配器结扎 DNA。在室温下孵育5分钟, 在磁性支架上放置管。当溶液清晰时, 将上清液转移到新管, 丢弃含有珠子的管子。
   注: 上清包含适配器结扎的 DNA 和大不必要的碎片绑定到丢弃的珠子。
7. 添加25µL 悬浮 SPRI 珠到上清从步骤 5.6, 混合和孵化5分钟室温。当溶液清晰时, 将管子放在磁性支架上, 并丢弃上清液。
   注: 被丢弃的上清液含有不需要的 dna, 而适配器结扎的 dna 则绑定在珠子上。
8. 当管在磁性立场, 添加200µL 80% 乙醇 (新鲜准备), 以洗涤珠。孵化三十年代, 除去和丢弃上清。重复这一步骤, 共2洗涤。
9. 空气干珠5分钟, 而管是在磁性立场与一个开放的盖子。
10. 从磁铁中取出管, 洗脱23µL 洗脱缓冲液中的 DNA。室温孵育2分钟。
11. 将管子放在磁性支架上, 当溶液清晰时, 将2µL 转移到96井 pcr 板上进行预 pcr 放大 (步骤 5.14)。
12. 将其余的大约21µL 转移到一个新的0.2 毫升的 PCR 放大管 (步骤 5.17)。
    注: 确保不转移任何珠子, 因为珠子可以抑制酶反应的下一步。
13. 继续步骤5.14 或在20°c 存储样品。
14. 用 RT qPCR 进行预 PCR 放大, 以估计放大适配器结扎 DNA 所需的周期数。每反应, 添加以下在96井 PCR 板包含2µL DNA (步骤 5.11): 0.4 µL 指数底漆, 0.4 µL 通用 PCR 底漆, 7.2 µL H₂O, 10 µL RT-qPCR 主混合。
15. 将板材置于实时 PCR 机中, 并将板材置于以下情况下: 三十年代在98°c (十年代98°c; 七十五年代65°c) 20 个周期; 5分钟65°c。
16. 估计每个样本所需的放大周期数, 将一个 ct 值减去与 PCR 反应线性相对应的预放大 RT qPCR 的 ct 值。
17. 从步骤5.12 中放大的适配器结扎的 dna 混合:21 µL dna, 25 µL pcr 总混合, 1 µL 指数底漆, 1 µL 通用 pcr 底漆和2µL H₂O。在热循环仪, 主题每个样品到以下情况: 三十年代在98°c (十年代98°c; 七十五年代65°c) [从步骤 5.16] 周期计算的 Ct; 5分钟65°c。
    注意: 请记住, 每个样本只使用一个索引底漆, 以允许多路复用。索引在 PCR 期间入。
18. 在22µL 洗脱缓冲器中纯化扩增的 DNA SPRI 珠和洗脱。
19. 采用凝胶电泳或替代方法控制库质量。寻找正确的库峰值大小 (PCR 产品的大小加上适配器)。
20. 通过涂抹分析估计库产品的平均基对大小: 在高级全局设置中, 在 "涂抹分析" 下双击 "表"。从 100-1000 bp 定义区域, 然后单击确定。在区域表下, 将显示定义的区域, 并计算库的平均大小 (bp)。
21. 通过荧光量化库。继续步骤6或存储库在-20 °c。

6. 下一代测序

1. 使用公式计算库的 nM 浓度:
   
2. 稀释图书馆到同一毫微米浓度例如2毫微米 (选择4毫微米, 2 毫微米, 1 毫微米, 或0.5 毫微米取决于图书馆浓度)。汇集所有库, 以达到 2 nM 库的池。
3. 变性和稀释图书馆跟随测序仪器指南。简而言之: 变性 2 nM 池添加 0.2 N 氢氧化钠, 在室温下孵育5分钟. 稀释变性池, 最后稀释下午10点。变性和稀释一个商业可利用的控制图书馆到最后稀释 12.5 pM。到新的管, 混合下午10点池与 20% 12.5 pM 控制库。
   注: 亚硫酸氢转换引入非常低的复杂性。20% 控制库的峰值将导致更好的排序质量。
4. 使用深测序仪运行 150 bp 配对端测序。

7. 计算分析-估计甲基化水平

1. 生成. fastq 文件 (此处称为 sample1_R1.fastq.gz 和 sample1_R2.fastq.gz), 生成所有感兴趣的基因组的俾斯麦索引 (此处位于名为bismarkIndex的文件夹中), 并确保 chrM 的基因组序列可在 fasta 中使用。格式 (此处位于名为chrM的文件夹中)。
2. 预处理读取以删除引物引入的适配器和偏差: 使用工具 Trim_galore³³。从正向和反向读取的 5 ' 端移除两个核苷酸。
   trim_galore--clip_R1 2-clip_R2 2-trim1-配对 sample1_R1.fastq.gz sample1_R2.fastq.gz
3. 使用俾斯麦³⁴将预处理的读取与整个基因组对齐
   俾斯麦--燕尾/bismarkIndex/-1/sample1_R1_val_1.fq.gz-2/sample1_R2_val_2.fq.gz
   注意: 如果使用了另一个光刻, 请确保丢弃不能唯一对齐的读取。
4. 提取读取映射到染色体 M, 这里使用Samtools³⁵
   samtools 排序-l 0-o sample1_sorted. bam
   sample1_R1_val_1_bismark_bt2_pe. bam
   samtools 索引 sample1_sorted. bam
   samtools 视图-u sample1_sorted. bam chrM |samtools 排序-n o sample1_chrM. bam
5. 提取甲基化信息
   bismark_methylation_extractor-cytosine_report--全面的 genome_folder/chrM/sample1_chrM. bam
6. 使用. bedGraph、. CX_report 或 .txt 文件进行进一步分析。使用 sample1_chrM.CX_report.txt 文件进行可视化和测试。
7. 在预处理步骤中执行进一步的裁剪, 如果询问多个区域, 并且在映射过程中生成的 M 偏置图中可以看到底漆偏倚。有关详细信息, 请参阅俾斯麦文档。

8. 计算分析. 差异检验

1. 确保 CX_report 文件包含7列, 染色体 (这里 chrM), 位置, 链, 甲基化读数的数量, 非甲基化读取的数量, c 上下文和周围序列, 每 C 染色体 m。
2. 当 CX_report 覆盖 chrM 的整体, 检查子集的区域 (s) 被放大。
3. 对样本之间的差异使用非参数测试, 例如渔民对单个 C 的确切测试或整个区域的标志测试。
   注: 量化通常会接近0% 甲基化, 这就是为什么假设正常不合适的原因。
4. 同样地, 对于可视化, 可以想象分位数或计算二项比例置信区间。

结果

在研究线粒体甲基化时, 该协议的两个步骤至关重要。1) 二次结构的开放性和 2) 线粒体 DNA 特异引物的设计。

通过使用限制酶 BamHI (图 1) 消化人类基因组 dna, 线粒体 dna 结构将在核苷酸位置14258上被切割, 第二结构将被打开。

亚硫酸氢转化的可能损害与完整的?...

讨论

在这里, 我们提供了一种亚硫酸氢测序协议, 专门用于审问线粒体甲基化。与亚硫酸氢酶测序协议用于基因组 DNA 的区别在于利用了先前的限制酵素消化步骤和生物信息学分析, 排除了 NUMT 序列产生的假阳性。

在研究线粒体甲基化时, 我们提供了一种避免亚硫酸氢测序工件的协议。以前报告的亚硫酸氢亚硫酸盐序列工件导致了误报³⁷。蛋白酶 K 去除 DNA 副蛋白的...

材料

Name	Company	Catalog Number	Comments
BamHI	New England BioLabs	# R0136
EZ DNA methylation-lightning kit	Zymo Research	# D5030
Qubit ssDNA assay	Thermo Fisher Scientific	# Q10212
Qubit assay tubes	Thermo Fisher Scientific	# Q32856
HotStarTaq plus DNA polymerase kit	Qiagen	# 203603
QIAquick Gel Extraction Kit	Qiagen	# 28704
NEBNext Ultra DNA Library Kit for Illumina	New England BioLabs	# E7370S
NEBNext Multiplex Oligoes for Illumina	New England BioLabs	# E7335S
AMPure XP Beads	Beckman Coulter	# A63881
High Sensitivity DNA chip	Agilent	# 5067-4626
Qubit high sensitivity dsDNA assay	Thermo Fisher Scientific	# Q33230
MiSeq reagent kit v2 300 cycles	Illumina	# MS-102-2003
PhiX control v3	Illumina	# FC-110-3001
Sodium hydroxide	Sigma	# S5881
Thermal cycler C1000	Biorad	# 1851148
CFX96 Real-Time PCR detection system	Biorad	#1855195
Qubit Fluorometer	Thermo Fisher Scientific	# Q33226
Bioanalyzer 2100	Agilent	# G2939BA
MiSeq instrument	Illumina	# SY-410-1003