Summary
在这里, 我们提出了一个协议, 允许准确定量的线粒体 DNA (线粒体) 甲基化。在本协议中, 我们描述了一种酶消化 DNA 与BamHI结合的生物信息学分析管道, 可用于避免高估线粒体甲基化水平引起的线粒体二级结构。
Abstract
利用亚硫酸氢钠的特性, 将非甲基化胞嘧啶转化为尿, 在单链 dna 环境下, 可以实现 dna 甲基化的量化。亚硫酸氢的测序可以是靶向 (使用 PCR) 或在整个基因组进行, 并提供绝对定量的胞嘧啶甲基化的单基分辨率。考虑到细胞核和线粒体 DNA 的独特性质, 特别是在二级结构中, 应该对亚硫酸氢嘧啶测序方法在线粒体中的 adaptions 进行研究。线粒体的二级和三级结构确实可以导致亚硫酸氢的测序工件导致误报, 因为不完全变性, 亚硫酸氢对单链 DNA 的接触较差。在这里, 我们描述了一个协议, 使用酶消化 DNA 与BamHI结合生物信息学分析管道, 以允许准确定量的胞嘧啶甲基化水平在线粒体。此外, 我们还提供了设计特定于线粒体的亚硫酸盐测序底漆的准则, 以避免将不理想的核线粒体片段 (NUMTs) 插入核基因组。
Introduction
线粒体基因组是一个圆形, 双链结构约16.5 公斤的基础 (kb) 长, 构成一个沉重的和轻的链。线粒体基因组存在于每个细胞内的多个拷贝中, 母系遗传, 并编码呼吸链复合体的基本成分1。与细菌基因组相似, 与核基因组不同, 线粒体基因组组织在许多次和三级结构中, 例如在连续和超螺旋结构中 2, 这可以在排序过程中使访问变得困难.实验3。
在细胞核中, DNA 甲基化是一种广泛研究的表观遗传标记, 在许多过程中起着重要作用, 特别是在基因表达调控中。在哺乳动物基因组中, DNA 甲基化主要发生在 deoxycytidines 嘧啶环的5位位置上, 主要是在 CG dinucleotides (或 CpG) 上。胞嘧啶甲基化在体细胞基因组中的70% 都被发现, 占总 DNA 基数的 1%, 为4。DNA methylations 也被描述在非 CpG 的情况下, 如 CpA, CpT 和 CpC, 并存在于各种数量的核 DNA, 其价值高达25% 的甲基化 cytosines 在胚胎干细胞5,6,7。
虽然核基因组的胞嘧啶甲基化被广泛接受, 但线粒体 DNA (线粒体) 甲基的存在仍然存在争议。首次研究线粒体甲基化是在培养细胞中进行的, 在那里线粒体甲基化是容易检测到的, 虽然在较低的水平与核 DNA8。在人类和小鼠细胞中, 线粒体甲基化也被检测到低水平 (2-5%)。使用依赖于5甲基胞嘧啶捕获 (如甲基化 DNA 免疫沉淀 (MeDIP) 和定量 PCR 的检测方法, 在各种鼠标和人和细胞线(9,10 ,) 中也发现线粒体甲基化.11,12。在 ELISA 检测或质谱中使用 5-甲基胞嘧啶抗体, 从纯化的线粒体分数中检测出大量的 DNA 甲基化水平13, 14, 15,16. 然而, 上述研究中的大多数化验方法都使用的技术不是为在单一碱基分辨率下的 DNA 甲基化提供绝对量化的。
定量和 resolutive dna 甲基化分析可以通过一种名为 "亚硫酸氢" 的技术来实现, 利用亚硫酸氢钠的特性将非甲基化胞嘧啶转化为尿的单链 DNA上下文17.利用亚硫酸氢的测序, 一系列研究发现了胞嘧啶甲基化在不同层次的存在。甲基化线粒体在 D 环区域, 12S 或16S 区域在人的18,19,20,21,22,23和老鼠24中容易地被发现。然而, 组织和细胞具有令人着迷的变异, 在研究中总 cytosines 的1-20%。
与这些众多的研究相比, 只有少数研究, 包括从我们的小组, 有争议的存在线粒体甲基化3,25,26,27或质疑的生物学相关性非常低的线粒体水平 (低于 2%)28。最近, 我们报告了一个潜在的亚硫酸氢序列伪影在整个线粒体亚硫酸氢化序列3的观察。我们提供的证据表明, 线粒体 DNA 的二级结构可能导致亚硫酸氢亚硫酸盐测序中的误报, 从而高估甲基化水平。我们这里提供了一个协议, 以防止亚硫酸氢的伪影转换的线粒体。该协议使用简单的酶消化 DNA, 以扰乱线粒体二级结构, 并允许完全进入亚硫酸氢化后, 亚硫酸盐排序协议。此外, 我们还提供了一个伴随的生物信息学管道, 用于分析亚硫酸氢的测序。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. 限制酶治疗
- 进行线性化线粒体 dna 通过治疗总人类 DNA 与限制酶 BamHI, 那削减在位置14258在人类线粒体脱氧核糖核酸。
注: 对于小鼠 DNA, 使用限制酶 BglII 在相同条件下, 如下所示。 - 对于要评估线粒体甲基化的每个样本, 准备一个0.2 毫升反应管, 并添加以下 mix:3 微克基因组 DNA (由荧光量化), 15 ul 缓冲 3, 3 ul BamHI 和水高达 150 ul
- 在热循环仪中放置管4小时, 在37摄氏度。
2. 亚硫酸氢转化
- 使用亚硫酸氢钠转换 BamHI 处理的 DNA, 如其他地方所述29。
- 使用 200-500 ng BamHI 处理的 DNA 为每个转换反应在总容量 20 ul。DNA 回收率为50-70%。
- 洗脱 10 ul 洗脱缓冲液中的亚硫酸氢转化 DNA。
- 用荧光法定量测定单链 DNA。
- 可选: 在继续执行该协议之前, 将亚硫酸氢化转化后的 DNA 储存在-20 摄氏度。长期贮存, 将亚硫酸氢转化后的 DNA 储存在-80 摄氏度
3. 亚硫酸氢化测序底漆的设计
- 用于亚硫酸氢的设计底漆使用在线工具 Methprimer30和折半查找 31, 32, 以此进行排序.
注: Methprimer 和折半查找允许搜索 cytosines 转化为 thymines 的基因组序列上的底漆结合序列, 类似于亚硫酸氢后处理。 - 为优化和无偏放大, 避免在引物中的 CpG dinucleotides 和设计的扩增子大小之间的 100-300 bp。
- 确保设计的底漆是特定的线粒体 DNA 使用折半查找的功能底漆搜索。插入已设计的亚硫酸氢转化底漆, 勾选亚硫酸氢框, 并包括参考基因组和 PCR 参数。
注: 将显示可能的 PCR 产品列表。 - 复制前1-5 个底漆序列, 并粘贴到一个排序形式的寡核苷酸合成。
注意: 在实验室中已验证的人和鼠线粒体 dna 底漆序列的列表显示在表 1中。 - 根据步骤4中的说明, 采用亚硫酸氢化凝胶电泳法进行的亚硫酸盐转化 PCR 试验底漆。对所有底漆对分别进行测试和优化, 确保琼脂糖凝胶上出现正确的扩增子尺寸。
注: 如果需要复用底漆, 添加多个引物到一个单一的 pcr 反应和可视化 pcr 产品在琼脂糖凝胶电泳或, 更 resolutive 的方法, 如聚丙烯酰胺凝胶电泳, 以区别产品的相似大小。
4. 亚硫酸氢转化 PCR 和凝胶萃取
- 为了放大感兴趣的区域, 使用热启动 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。单独或复用底漆。 - 可选: 当复用引物, 使用 200 ng 转换 DNA 和以下循环条件: 5 分钟95°c (六十年代94°c; 九十年代55°c *; 九十年代72°c) 35 循环; 10分钟72°c。
注: 温度可能因底漆熔化温度而异。 - 主题 PCR 产品的凝胶电泳2% 琼脂糖在100伏特。
- 在柔和的紫外光下, 用手术刀提取 PCR 产物并加入管内。不要将 PCR 产品暴露在过量的紫外线下, 以避免 DNA 损伤。避免曝光时间超过三十年代。
- 使用凝胶纯化方法纯化 PCR 产品, 如前所述3。
- 用荧光定量测定步骤4.5 中纯化的 DNA。继续步骤5或在20°c 存储样品。
10分钟72°c。单独或复用底漆。5. 亚硫酸氢亚硫酸盐测序库准备
- 对亚硫酸氢转化 PCR 产物进行库制备。
可选: 在这个阶段多路 PCR 产品。 - 使用多达 100 ng 的 PCR 产品进行端修复。添加3µL 端准备酶和6.5 µL 端修复反应缓冲器 (10x) 到55.5 µL PCR 产品在0.2 毫升管。放置和孵化管在热循环仪为30分钟在20°c 之后30分钟在65°c。
- 结扎测序适配器通过混合最终修复 DNA 与15µL 钝/TA 连接酶混合, 2.5 µL 适配器和1µL 结扎增强剂。如果使用 < 100 ng PCR 产品作为输入, 稀释适配器10次。
- 将混合 (步骤 5.3) 放在热循环仪中, 在20摄氏度孵化15分钟。暂停热循环仪, 并添加3µl 用户酶的管。混合并放置在热循环仪的管, 孵化15分钟在37摄氏度。
- 尺寸选择适配器结扎的 dna, 使用固体相可逆固定 (SPRI) 珠的不同比例的珠子, DNA 取决于 PCR 产品大小。
- 添加13.5 µL H2O, 从步骤5.4 和55µL 悬浮 SPRI 珠中的适配器结扎 DNA。在室温下孵育5分钟, 在磁性支架上放置管。当溶液清晰时, 将上清液转移到新管, 丢弃含有珠子的管子。
注: 上清包含适配器结扎的 DNA 和大不必要的碎片绑定到丢弃的珠子。 - 添加25µL 悬浮 SPRI 珠到上清从步骤 5.6, 混合和孵化5分钟室温。当溶液清晰时, 将管子放在磁性支架上, 并丢弃上清液。
注: 被丢弃的上清液含有不需要的 dna, 而适配器结扎的 dna 则绑定在珠子上。 - 当管在磁性立场, 添加200µL 80% 乙醇 (新鲜准备), 以洗涤珠。孵化三十年代, 除去和丢弃上清。重复这一步骤, 共2洗涤。
- 空气干珠5分钟, 而管是在磁性立场与一个开放的盖子。
- 从磁铁中取出管, 洗脱23µL 洗脱缓冲液中的 DNA。室温孵育2分钟。
- 将管子放在磁性支架上, 当溶液清晰时, 将2µL 转移到96井 pcr 板上进行预 pcr 放大 (步骤 5.14)。
- 将其余的大约21µL 转移到一个新的0.2 毫升的 PCR 放大管 (步骤 5.17)。
注: 确保不转移任何珠子, 因为珠子可以抑制酶反应的下一步。 - 继续步骤5.14 或在20°c 存储样品。
- 用 RT qPCR 进行预 PCR 放大, 以估计放大适配器结扎 DNA 所需的周期数。每反应, 添加以下在96井 PCR 板包含2µL DNA (步骤 5.11): 0.4 µL 指数底漆, 0.4 µL 通用 PCR 底漆, 7.2 µL H2O, 10 µL RT-qPCR 主混合。
- 将板材置于实时 PCR 机中, 并将板材置于以下情况下: 三十年代在98°c (十年代98°c; 七十五年代65°c) 20 个周期; 5分钟65°c。
- 估计每个样本所需的放大周期数, 将一个 ct 值减去与 PCR 反应线性相对应的预放大 RT qPCR 的 ct 值。
- 从步骤5.12 中放大的适配器结扎的 dna 混合:21 µL dna, 25 µL pcr 总混合, 1 µL 指数底漆, 1 µL 通用 pcr 底漆和2µL H2O。在热循环仪, 主题每个样品到以下情况: 三十年代在98°c (十年代98°c; 七十五年代65°c) [从步骤 5.16] 周期计算的 Ct; 5分钟65°c。
注意: 请记住, 每个样本只使用一个索引底漆, 以允许多路复用。索引在 PCR 期间入。 - 在22µL 洗脱缓冲器中纯化扩增的 DNA SPRI 珠和洗脱。
- 采用凝胶电泳或替代方法控制库质量。寻找正确的库峰值大小 (PCR 产品的大小加上适配器)。
- 通过涂抹分析估计库产品的平均基对大小: 在高级全局设置中, 在 "涂抹分析" 下双击 "表"。从 100-1000 bp 定义区域, 然后单击确定。在区域表下, 将显示定义的区域, 并计算库的平均大小 (bp)。
- 通过荧光量化库。继续步骤6或存储库在-20 °c。
6. 下一代测序
- 使用公式计算库的 nM 浓度:
- 稀释图书馆到同一毫微米浓度例如2毫微米 (选择4毫微米, 2 毫微米, 1 毫微米, 或0.5 毫微米取决于图书馆浓度)。汇集所有库, 以达到 2 nM 库的池。
- 变性和稀释图书馆跟随测序仪器指南。简而言之: 变性 2 nM 池添加 0.2 N 氢氧化钠, 在室温下孵育5分钟. 稀释变性池, 最后稀释下午10点。变性和稀释一个商业可利用的控制图书馆到最后稀释 12.5 pM。到新的管, 混合下午10点池与 20% 12.5 pM 控制库。
注: 亚硫酸氢转换引入非常低的复杂性。20% 控制库的峰值将导致更好的排序质量。 - 使用深测序仪运行 150 bp 配对端测序。
7. 计算分析-估计甲基化水平
- 生成. fastq 文件 (此处称为 sample1_R1.fastq.gz 和 sample1_R2.fastq.gz), 生成所有感兴趣的基因组的俾斯麦索引 (此处位于名为bismarkIndex的文件夹中), 并确保 chrM 的基因组序列可在 fasta 中使用。格式 (此处位于名为chrM的文件夹中)。
- 预处理读取以删除引物引入的适配器和偏差: 使用工具 Trim_galore33。从正向和反向读取的 5 ' 端移除两个核苷酸。
trim_galore--clip_R1 2-clip_R2 2-trim1-配对 sample1_R1.fastq.gz sample1_R2.fastq.gz - 使用俾斯麦34将预处理的读取与整个基因组对齐
俾斯麦--燕尾/bismarkIndex/-1/sample1_R1_val_1.fq.gz-2/sample1_R2_val_2.fq.gz
注意: 如果使用了另一个光刻, 请确保丢弃不能唯一对齐的读取。 - 提取读取映射到染色体 M, 这里使用Samtools35
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 - 提取甲基化信息
bismark_methylation_extractor-cytosine_report--全面的 genome_folder/chrM/sample1_chrM. bam - 使用. bedGraph、. CX_report 或 .txt 文件进行进一步分析。使用 sample1_chrM.CX_report.txt 文件进行可视化和测试。
- 在预处理步骤中执行进一步的裁剪, 如果询问多个区域, 并且在映射过程中生成的 M 偏置图中可以看到底漆偏倚。有关详细信息, 请参阅俾斯麦文档。
8. 计算分析. 差异检验
- 确保 CX_report 文件包含7列, 染色体 (这里 chrM), 位置, 链, 甲基化读数的数量, 非甲基化读取的数量, c 上下文和周围序列, 每 C 染色体 m。
- 当 CX_report 覆盖 chrM 的整体, 检查子集的区域 (s) 被放大。
- 对样本之间的差异使用非参数测试, 例如渔民对单个 C 的确切测试或整个区域的标志测试。
注: 量化通常会接近0% 甲基化, 这就是为什么假设正常不合适的原因。 - 同样地, 对于可视化, 可以想象分位数或计算二项比例置信区间。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
在研究线粒体甲基化时, 该协议的两个步骤至关重要。1) 二次结构的开放性和 2) 线粒体 DNA 特异引物的设计。
通过使用限制酶 BamHI (图 1) 消化人类基因组 dna, 线粒体 dna 结构将在核苷酸位置14258上被切割, 第二结构将被打开。
亚硫酸氢转化的可能损害与完整的线粒体二级结构很可能高估线粒体 unconversion 率。通过亚硫酸氢化酶测序, 研究了在5个不同区域的线粒体基因组: 位移环 (d-环)、tRNA 苯丙氨酸和12S 核糖体中, 在未消化和消化的人骨骼肌细胞的总 DNA 中, 线粒体 unconversion 率。(tRNA-F+12S), 16S 核糖体 (16S), NADH 脱氢酶 5 (ND5) 和细胞色素 b (CYTB) mRNA 编码基因 (图 2)。Unconversion 率从所有被调查的区域 0-15.1% 不等于未消化的线粒体。然而, 当 DNA 在亚硫酸氢转换之前被消化时, 所有区域的 unconversion 速率降到最大值 1% (图 3)。因此, 说明在亚硫酸氢转化之前 BamHI 消化的重要性, 以避免高估线粒体甲基化水平。
在分离线粒体基因组时, 核 DNA 的污染是不可避免的, 因为几乎整个线粒体基因组都复制到了核基因组中, 它导致了线粒体基因组的核插入, 称为核线粒体片段 (NUMTs)36。为了确保所分析的 dna 甲基化水平与线粒体 DNA 相对应, 而不是 NUMTs, 对设计特定于线粒体基因组的引物具有最外层的重要性。图 4显示了如何检查底漆特异性和表 1显示人和老鼠线粒体的底漆序列。
| 名字 | 序列 | 位置 | 扩增子尺寸 | 退火温度 | 链 | ||
| 人类 | |||||||
| D 循环 | F: AAATCTATCACCCTATTAAC | ||||||
| R: GTGGAAATTTTTTGTTATGATGT | 6-298 | 292 | 55°C | 反 义 | |||
| D 循环 | F: CATAACAAAAAATTTCCACCAAAC | ||||||
| R: GGGAAAATAATGTGTTAGTT | 279-458 | 179 | 55°C | 反 义 | |||
| 十二年代 + TF | F: TTTATATAACTTACCTCCTC | ||||||
| R: GTGTTTGATGTTTGTTTTTTTTG | 577-765 | 188 | 55°C | 反 义 | |||
| 466 | F: AATAAATTTATAGGTTTTTAAATTATTAAAT | ||||||
| R: TAACTAATAAAATCTTAACATATACTACTC | 2763-2873 | 110 | 53°C | 感 | |||
| CYTB | F: GGTATTATTTTTTTGTTTGTAATTATAGTA | ||||||
| R: CCTCAAATTCATTAAACTAAATCTATCC | 15091-15243 | 152 | 55°C | 感 | |||
| 鼠标 | |||||||
| 426 | F: ACACATACAAACCTCCATAAAC | ||||||
| R: GAGGTATAATTTAGTTAAAT | 111-287 | 176 | 53°C | 反 义 | |||
| 466 | F: AAATTCCAATTCTCCAAACATAC | ||||||
| R: TTTGGATTTTTTTTTTAGGT | 1749-1896 | 147 | 53°C | 反 义 | |||
| CYTB | F: CTACAAAAACACCTAATAACAAAC | ||||||
| R: AGTGTATGGTTAAGAAAAGA | 14129-14307 | 178 | 55°C | 反 义 | |||
| 466 | F: AGAGAAATAGAGTTATTTTATAAATAAG | ||||||
| R: CTAAAAACAAAATTTTAAATCTTAC | 2576-2727 | 151 | 55°C | 感 | |||
| MDV | F: TTAAGTTAATTAGGTTTGATAATAGTGA | ||||||
| R: TAAAACCCTATTAAAAATAATATTCCTATA | 12659-12910 | 251 | 55°C | 感 | |||
| CYTB | F: TGGGTTTTTTTTAGGAGTTTGTTTA | ||||||
| R: CAATAAAAATACTCCAATATTTCAAATTTC | 14242-14504 | 262 | 55°C | 感 | |||
表 1: 人类和小鼠线粒体的引物序列.请注意, PCR 引物的退火温度因引物熔化温度的不同而异。
图 1: BamHI 消化或未消化的 DNA 的凝胶电泳.人骨骼肌的基因组 DNA 未经治疗或用 BamHI 治疗, 并在1.5% 琼脂糖凝胶上进行凝胶电泳。请单击此处查看此图的较大版本.
图 2: 人类线粒体基因组的地图.用亚硫酸氢酶测序法和 BamHI 限制酵素研究了线粒体基因组的区域。线粒体 DNA 包含13基因、2核糖体 rna 和22个转移 rna。缩写: PH1: 重链启动子 1;PH2: 重链提升 2;PL: 轻型链启动子;OH: 复制的起源从重的线;OL: 从光链复制的起源。请单击此处查看此图的较大版本.
图 3: 在亚硫酸氢转化之前的 BamHI 消化降低了人骨骼肌细胞的线粒体 unconversion 率.通过亚硫酸氢化酶测序, 在人体骨骼肌细胞线粒体基因组的五个不同区域 (N= 3) 审问了未消化和消化的线粒体 DNA unconversion 百分比。整个正方形代表 cpg 站点, 并且开放正方形代表非 CpG 站点。给出了最小最大区间的结果, 用符号检验检验了 unconversion 的显著差异。D 环路 (6-298) P= 1.83E-42;D 环路 (279-458): P= 4.55E-13;tRNA-F+12S: P= 8.38E-16;16S: P= 5.91E-06;ND5: P= 1.70E-05;CYTB: P= 2.69E-10。D 环 (6-298) 包括复制的起源和 tRNA-F+12S 包括重链启动子2。这些数据已发布到其他地方3。请单击此处查看此图的较大版本.
图 4: 使用折半查找验证对线粒体基因组的底漆特异性。A)亚硫酸氢转化底漆序列添加到折半查找功能底漆搜索 (http://bisearch.enzim.hu/？m=genompsearch) 中。亚硫酸氢盒被打勾, 并选择一个参考基因组。B)示例在感官和反义链上生成的潜在 PCR 产品。请单击此处查看此图的较大版本.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
在这里, 我们提供了一种亚硫酸氢测序协议, 专门用于审问线粒体甲基化。与亚硫酸氢酶测序协议用于基因组 DNA 的区别在于利用了先前的限制酵素消化步骤和生物信息学分析, 排除了 NUMT 序列产生的假阳性。
在研究线粒体甲基化时, 我们提供了一种避免亚硫酸氢测序工件的协议。以前报告的亚硫酸氢亚硫酸盐序列工件导致了误报37。蛋白酶 K 去除 DNA 副蛋白的不足已被描述为37。不完全变性或部分 reannealing 可能导致不完全转化的延伸, 可能是通过降低亚硫酸氢的进入单链 DNA37。后来的伪影可以在线粒体中发挥作用, 在这种情况下, 次级结构会损害对 cytosines 的接触。据我们所知, 在估计线粒体甲基化的情况下, 先在2016年的3、28中首次报告了在亚硫酸氢转化之前添加限制酶消化步骤。在我们的手太, 事先消化与单刀限制酶显着降低检测未转换 cytosines3,28。
在这里, 我们提供了一个生物信息学的管道, 以分析亚硫酸氢的测序结果, 并在整个线粒体基因组亚硫酸氢化序列的背景下检测亚硫酸氢测定序列。这种方法来源于我们的观察, 即线粒体的非转化率与测序深度呈反比关系3,28。这一观察表明, 线粒体的二级甚至三级结构会在测序库的构建过程中损害特定区域线粒体的碎片化, 也会损害亚硫酸氢钠对线粒体的充分利用。由于消化与一刀酶释放二级和三级结构, 我们的结果表明, 线粒体结构是氢氧化物的来源, 并验证本文所描述的方法, 以定量的胞嘧啶甲基化的线粒体。
线粒体制剂中核 DNA 的不可避免的污染是一个挑战。这导致序列读数中 NUMT 序列的污染, 这可能会在线粒体中对胞嘧啶甲基化的估计, 即使是在执行靶向的, 而不是全线粒体基因组亚硫酸氢的序列。为避免这种偏差, 设计独特的亚硫酸氢化物测序底漆可以排除在 NUMTs 的甲基化污染。NUMTs 跨越人类和小鼠线粒体基因组的某些区域, 具有100% 的同一性, 因此, 在某些线粒体区域设计独特的引物是不可能的。在执行全线粒体基因组亚硫酸氢测定序列时, 测序长读可以缓解线粒体和 NUMTs 之间的差别。将读数对准整个基因组, 只使用唯一映射到线粒体基因组的读数, 可以进一步确保检测线粒体甲基化。
最后, 这个协议可以在线粒体基因组之外使用, 例如, 在其他亚硫酸氢的情况下, DNA 的二级结构代表一个潜在的工件来源。可以研究甲基化水平和测序深度之间的关系, 以评估在复杂 DNA 序列的亚硫酸氢酶测序之前, 是否需要用限制酵素进行消化。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
作者没有竞争的经济利益来申报。
Acknowledgments
德基础代谢研究基金会是在哥本哈根大学的一个独立研究中心, 由来自德基金会的无限制捐赠部分资助。
Materials
| 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 |
References
- Falkenberg, M., Larsson, N. -G., Gustafsson, C. M. DNA replication and transcription in mammalian mitochondria. Annu Rev Biochem. 76, 679-699 (2007).
- Kolesar, J. E., Wang, C. Y., Taguchi, Y. V., Chou, S. -H., Kaufman, B. A. Two-dimensional intact mitochondrial DNA agarose electrophoresis reveals the structural complexity of the mammalian mitochondrial genome. Nucleic Acids Res. 41 (4), e58 (2013).
- Mechta, M., Ingerslev, L. R., Fabre, O., Picard, M., Barres, R. Evidence Suggesting Absence of Mitochondrial DNA Methylation. Front Genet. 8, 166 (2017).
- Ehrlich, M., Gama-Sosa, M. A., et al. Amount and distribution of 5-methylcytosine in human DNA from different types of tissues of cells. Nucleic Acids Res. 10 (8), 2709-2721 (1982).
- Lister, R., Pelizzola, M., et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 462 (7271), 315-322 (2009).
- Yan, J., Zierath, J. R., Barres, R. Evidence for non-CpG methylation in mammals. Exp Cell Res. 317 (18), 2555-2561 (2011).
- Patil, V., Ward, R. L., Hesson, L. B. The evidence for functional non-CpG methylation in mammalian cells. Epigenetics. 9 (6), 823-828 (2014).
- Nass, M. K. Differential methylation of mitochondrial and nuclear DNA in cultured mouse, hamser and virus trasnforemd hamser cells in vivo and in vitro methylation. J Mol Biol. 80 (1), 155-175 (1973).
- Shock, L. S., Thakkar, P. V., Peterson, E. J., Moran, R. G., Taylor, S. M. DNA methyltransferase 1, cytosine methylation, and cytosine hydroxymethylation in mammalian mitochondria. PNAS. 108 (9), 3630-3635 (2011).
- Bellizzi, D., D'Aquila, P., et al. The Control Region of Mitochondrial DNA Shows an Unusual CpG and Non-CpG Methylation Pattern. DNA Res. 20 (6), 537-547 (2013).
- Jia, Y., Li, R., et al. Maternal Low-Protein Diet Affects Epigenetic Regulation of Hepatic Mitochondrial DNA Transcription in a Sex-Specific Manner in Newborn Piglets Associated with GR Binding to Its Promoter. PLoS ONE. 8 (5), e63855 (2013).
- Jia, Y., Song, H., Gao, G., Cai, D., Yang, X., Zhao, R. Maternal Betaine Supplementation during Gestation Enhances Expression of mtDNA-Encoded Genes through D-Loop DNA Hypomethylation in the Skeletal Muscle of Newborn Piglets. J Agric Food Chem. 63 (46), 10152-10160 (2015).
- Infantino, V., Castegna, A., Iacobazzi, F., Spera, I. Impairment of methyl cycle affects mitochondrial methyl availability and glutathione level in Down's syndrome. Mol Genet Metab. 102 (3), 378-382 (2011).
- Chen, H., Dzitoyeva, S., Manev, H. Effect of valproic acid on mitochondrial epigenetics. Eur J Pharmacol. 690 (1-3), 51-59 (2012).
- Dzitoyeva, S., Chen, H., Manev, H. Effect of aging on 5-hydroxymethylcytosine in brain mitochondria. Neurobiol Aging. 33 (12), 2881-2891 (2012).
- Menga, A., Palmieri, E. M., et al. SLC25A26 overexpression impairs cell function via mtDNA hypermethylation and rewiring of methyl metabolism. FEBS J. 284 (6), 967-984 (2017).
- Frommer, M., McDonald, L. E., et al. A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. PNAS. 89 (5), 1827-1831 (1992).
- Byun, H. -M., Panni, T., et al. Effects of airborne pollutants on mitochondrial DNA Methylation. Part Fibre Toxicol. 10 (1), 1 (2013).
- Janssen, B. G., Byun, H. -M., Gyselaers, W., Lefebvre, W., Baccarelli, A. A., Nawrot, T. S. Placental mitochondrial methylation and exposure to airborne particulate matter in the early life environment: An ENVIRONAGE birth cohort study. Epigenetics. 10 (6), 536-544 (2015).
- Zheng, L. D., Linarelli, L. E., et al. Insulin resistance is associated with epigenetic and genetic regulation of mitochondrial DNA in obese humans. Clin Epigenetics. 7 (1), 1739 (2015).
- Byun, H. -M., Colicino, E., Trevisi, L., Fan, T., Christiani, D. C., Baccarelli, A. A. Effects of Air Pollution and Blood Mitochondrial DNA Methylation on Markers of Heart Rate Variability. J Am Heart Assoc. 5 (4), (2016).
- Bianchessi, V., Vinci, M. C., et al. Mitochondrion. MITOCH. 27 (C), 40-47 (2016).
- Wijst, M. G. P., van Tilburg, A. Y., Ruiters, M. H. J., Rots, M. G. Experimental mitochondria-targeted DNA methylation identifies GpC methylation, not CpG methylation, as potential regulator of mitochondrial gene expression. Sci Rep. 7 (1), 177 (2017).
- Wong, M., Gertz, B., Chestnut, B. A., Martin, L. J. Mitochondrial DNMT3A and DNA methylation in skeletal muscle and CNS of transgenic mouse models of ALS. Front cellr Neurosci. 7, 279 (2013).
- Dawid, I. B. 5-methylcytidylic acid: absence from mitochondrial DNA of frogs and HeLa cells. Science. 184 (4132), 80-81 (1974).
- Gama-Sosa, M. A. Levels and distribution of 5-methylcytosine in Chordate DNA. , 1-201 (1985).
- Hong, E. E., Okitsu, C. Y., Smith, A. D., Hsieh, C. -L. Regionally specific and genome-wide analyses conclusively demonstrate the absence of CpG methylation in human mitochondrial DNA. Mol Cell Biol. 33 (14), 2683-2690 (2013).
- Liu, B., Du, Q., et al. CpG methylation patterns ofhuman mitochondrial DNA. Nature Publishing Group. , 1-10 (2016).
- Donkin, I., Versteyhe, S., et al. Obesity and Bariatric Surgery Drive Epigenetic Variation of Spermatozoa in Humans. Cell Metab. 23 (2), 369-378 (2016).
- Li, L. C., Dahiya, R. MethPrimer: designing primers for methylation PCRs. Bioinformatics. 18 (11), 1427-1431 (2002).
- Tusnády, G. E., Simon, I., Váradi, A., Arányi, T. BiSearch: primer-design and search tool for PCR on bisulfite-treated genomes. Nucleic Acids Res. 33 (1), e9 (2005).
- Arányi, T., Váradi, A., Simon, I., Tusnády, G. E.
- Krueger, F. Trim Galore!. , Available from: http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ (2018).
- Krueger, F., Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 27 (11), 1571-1572 (2011).
- Li, H., Handsaker, B., et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 25 (16), 2078-2079 (2009).
- Ramos, A., Barbena, E., et al. Nuclear insertions of mitochondrial origin: Database updating and usefulness in cancer studies. Mitochondrion. 11 (6), 946-953 (2011).
- Warnecke, P. M., Stirzaker, C., Song, J., Grunau, C., Melki, J. R., Clark, S. J. Identification and resolution of artifacts in bisulfite sequencing. Methods (San Diego, Calif). 27 (2), 101-107 (2002).
