有针对性的DNA甲基化分析通过新一代测序

摘要

亚硫酸氢盐扩增子测序 （BSAS） 是一种定量目标基因组区域中胞嘧啶甲基化的方法。该方法在下一代测序之前使用亚硫酸氢盐转化与靶区的 PCR 扩增配对，以产生碱基特异性水平的 DNA 甲基化绝对定量。

摘要

表观遗传过程在基因表达控制中的作用已经为人所知多年。胞嘧啶残基的 DNA 甲基化对表观遗传学研究特别感兴趣，因为它已被证明是基因表达的持久和动态调节因子。由于缺乏高通量、定量准确的方法，检查健康和疾病表观遗传变化的努力受到了阻碍。随着新一代测序 （NGS） 技术的出现和普及，除了现有的基因组和转录组学方法外，这些工具现在还被应用于表观基因组学。对于胞嘧啶甲基化的表观遗传学研究，其中已经确定了感兴趣的区域，例如特异性基因启动子或 CpG 岛，并且需要以高定量准确性检查大量样品，我们开发了一种称为亚硫酸氢盐扩增子测序 （BSAS） 的方法。该方法将亚硫酸氢盐转化与目标区域的靶向扩增、转座子体介导的文库构建和台式 NGS 相结合。BSAS 提供了一种快速有效的方法，可一次分析多达 96 个样品中高达 10 kb 的目标区域，大多数具有基本分子生物学技能的研究小组都可以进行该方法。结果提供了具有碱基特异性的胞嘧啶甲基化的绝对定量。BSAS 可应用于任何 DNA 来源的任何基因组区域。该方法可用于感兴趣靶标区域的假设检验研究以及确认在全基因组甲基化分析中鉴定的区域，例如全基因组亚硫酸氢盐测序、减少代表性亚硫酸氢盐测序和甲基化 DNA 免疫沉淀测序。

引言

自的天然存在的DNA的修饰中胞嘧啶甲基¹的形式中，第一报告它已经过了半个世纪。胞嘧啶甲基化是其中基环上的5 ^个碳具有一个甲基（5-MC），最常在哺乳动物基因组中的CpG二核苷酸基序中产生的修饰的胞嘧啶核苷酸碱基。的5-MC在基因的启动子功能的存在通常与转录抑制有关，而没有与转录活性^相关联。

巨大进展已在我们的DNA甲基化在发展^3，后生轮廓^4，癌症发病^5,6，和其他一些研究领域transgeneration传播中的作用的理解。许多的这些进步已经走在过去的几年中作为新的方法已经被开发来分析DNA甲基化^7。同的出现和下一代测序（NGS）技术的普及，现在有许多方法来跨越整个基因组或基因组⁸的大区域轮廓的DNA甲基化。这些方法打开发现研究，量化DNA甲基化模式和DNA甲基化的差异的可能性。除了这些假设产生的方法，有针对性/假设检验的DNA甲基化的定量方法显著需求。

焦磷酸测序，甲基化特异性PCR，和亚硫酸氢盐转化的DNA直接Sanger测序一直是最常用的方法为目标区域的分析（ 即 ，单一的基因的启动子区域或CpG岛^）9-12。所有这些方法都依赖于亚硫酸氢盐转化率，其中，非甲基化胞嘧啶的转化成尿嘧啶的化学脱氨反应，而甲基化胞嘧啶的保持完好^13,14。的bisulfite- PCR扩增转换后的DNA结果更换尿嘧啶与胸腺嘧啶¹⁵允许差异甲基化读出的基差。而非常有用，所述限制于这些方法包括低定量准确，短的读取长度，以及低样品通量。为了解决这些限制我们开发我们称之为亚硫酸氢盐测序扩增子^（BSAS）16的方法。该方法的目标是能够分析所关注目标的基因组区域中的大量的高定量准确的样品。

BSAS使用的现有方法的元素（亚硫酸氢盐转化率和区域特异性PCR扩增），并将其与简单下一代文库构建（座子介导的^）17和台式NGS ^18（ 图1）。这种方法提供了一个更定量和更高的吞吐量的方法，研究在任何感兴趣的区域中的胞嘧啶甲基化。在这里，我们提出了方法进行详细描述，并进行故障排除和质量检查BSAS过程的方法。

方案

从组织1.核酸分离

注：联合隔离DNA和RNA从同一组织样品提供了机会，收集的DNA进行分析后生和配对的RNA为基因表达分析。这里给出的例子中是从实验性的组织，但对不同类型的样品或其它分离方法的各种途径柱纯化的形式都可以使用。的主要要求是对高度纯化的核酸而不污染蛋白质或有机溶剂。

1. 选择一块新鲜冷冻或新鲜组织进行机械均化并放置在冰上。存储或转移组织，以2.0毫升的圆底离心管中。
2. 添加1个5毫米的不锈钢珠含组织各管。添加适当量的裂解缓冲液，根据制造商的细胞数或组织重量的建议。
3. 使用珠磨机机械均质化，在均化组织30赫兹进行30秒。
   注：频率和机械均质化的持续时间依赖于组织类型。较软的组织将彻底均匀使用推荐的设置，而更严格的组织要求设置的优化。
4. 使用硅胶离心柱在室温下制造商的方案进行DNA和RNA的分离。
   1. 洗脱的RNA在50μl不含RNA酶的水并置于冰上。使用RNA通过定量PCR基因表达分析。
   2. 洗脱在100μlTE缓冲液中的DNA。使用洗脱液以增加DNA的浓度和产量重复洗脱的DNA断列。利用DNA进行有针对性的甲基化定量分析。
5. 使用标准的分光光度计吸光度在230nm，260nm处，280处和320处评估的DNA和RNA的质量。
   注：A260 / A280的比例应为〜1.8-2高质量的DNA。 A260 / A280的比例应为〜2的高品质的RNA。过量吸光度在230nm印度语的存在茨污染物的分离步骤。
6. 质量通过根据制造商的说明使用毛细管电泳芯片进一步检查的RNA。
   注：高品质的RNA会产生RNA完整性数（RIN）> 8。在除了rRNA的峰的电泳额外峰的存在表明样品污染。
7. 根据制造商的协议量化DNA和RNA的荧光测定法。
   注意：使用荧光检测是更为敏感和特异单独量化核酸比分光光度法。
8. 商店中分离的DNA和RNA在-80℃。

2.目标识别和引物设计

1. 确定感兴趣的基因组区域（S）以被BSAS基于适当的参照基因组序列分析。差异表达基因的基因启动子甲基化是定量的共同目标。
2. 在硅片亚硫酸氢盐转换准备对于测序后比对ED的参考基因组的读取通过改变.fasta序列文件在文本编辑器。在5'-3'方向，取代非CpG的胞嘧啶与胸腺嘧啶的（ 图2）。
3. 设计底漆集从亚硫酸氢盐转化的DNA扩增感兴趣的区域。
   1. 选择并复制感兴趣的未转化的区域，在5'-3'方向，成亚硫酸氢特异性PCR引物设计程序（ 图3）。
      注：最优亚硫酸氢盐PCR扩增子长度为每扩增子250-400碱基对作为亚硫酸氢盐处理的DNA片段，它是难以扩增大> 400bp的区域。如果，例如，1 kb的区域是感兴趣，多个引物对，可设计成覆盖该区域。设计扩增子为等于bp的尺寸，便于池。避免扩增长度<250基点，因为这些可能是长度不足以库生成。最佳的引物长度为BSAS是>每个引物20基点。
   2. 前锋选择_{55-65℃，Tm}值范围和最大_Tm值之间的差异和反向1-2℃的引物。
   3. 选择的引物对，最好覆盖感兴趣的区域。不要选择含有CpG位点或者直接相邻CpG位点的引物，因为这会导致偏见的PCR反应。使用标准的PCR引物，可以从任意数量的学术或商业供应商订购。
   4. 重建在RNA酶/ DNA酶的游离水冻干引物以100微米的工作原液。
4. 商店的PCR引物在-20℃下。稀释100μM的引物股票10μM工作股票PCR反应。

3.重亚硫酸盐特异PCR优化

注：对于基因组DNA的亚硫酸氢盐转化，可用一个号码为亚硫酸氢盐转化不同的商业试剂盒。选择最适合的计划实验套件或协议。

1. 基因组DNA之间的200毫微克〜2微克使用。为优化实验，运行多个转化反应，以便有足够的亚硫酸氢盐转化的DNA用于多个基站的PCR反应。
2. 使用小（ 如 10微升）洗脱的亚硫酸氢盐转化的DNA量保持高浓度的DNA。
   注：1-2微升亚硫酸氢盐转化的DNA足以用于优化PCR反应，如果1微克基因组DNA，用于转换输入。
3. 通过放大亚硫酸氢特异性PCR亚硫酸氢盐转化的DNA。 BS-PCR反应需要使用Taq酶能够放大亚硫酸氢盐转化的DNA。
   1. 组装为单个扩增子的靶扩增优化以下反应。对于多个样品，组装在一个96孔PCR板的反应。
      25微升2倍反应缓冲液
      0.5微升dNTP混合物
      5微升10μM正向引物（最终1μM）
      5微升10μM反向引（最终1μM）
      2微升模板亚硫酸氢盐转化DNA
      0.4微升（5 U /微升）DNA聚合酶
      12.1微升DNA酶自由水（最终反应体积50微升）
   2. 密封PCR板用适当的粘合剂或热封膜。
   3. 将反应在适当的热循环仪，使用下列循环条件与热盖：
      1. 执行初始变性，在95℃下进行10分钟。
      2. 变性在95℃下30秒。
      3. 被用于退火30秒，在引物的特异性T _米 。在开始的退火温度少数℃以下的底漆优化反应的_Tm值。
         注：最佳的退火温度也可使用梯度热循环仪，以测试一个温度范围内进行测定。
      4. 执行一个延伸72℃30秒。延长更长的扩增子的延长时间。
      5. 重复步骤3.3.3.2-3.3.3.4为35个循环，以TAL的初步优化。更高的循环数目可能是必要的，但通常应该避免，以防止克隆扩增，将创建反应工件。
      6. 执行一个最终延伸在72℃7分钟。
      7. 保持反应在4℃。
   4. 可视化扩增通过PAGE电泳。或者，根据制造商的协议使用的毛细管电泳的DNA芯片。
      1. 混合24微升含6微升5×加载染料在冰上和旋涡的PCR反应。
      2. 使1升的1×TBE缓冲液通过混合双蒸₂ 200毫升5×TBE的运行缓冲液和800毫升O.稀释DNA阶梯，每孔0.5微克DDH _{2 O}的负载和染料的终浓度。
      3. 负载井30微升PCR反应或阶梯。
      4. 运行凝胶在200 V的〜45分钟。停止凝胶运行，当第一染料前沿到达凝胶的末端。
      5. 染色的凝胶用5微克/毫升溴化乙锭通过混合5微升的10毫克/毫升在100ml双蒸₂ O的覆盖整个凝胶，并让孵育〜5分钟。
      6. 图像通过在482纳米的激发波长的UV光盒或多峰成像器的凝胶。
      7. 在参考了梯子，大小的PCR扩增子，并通过寻找多个频带比预期的扩增子大小等（ 图4A和C）确定该反应的特异性。
   5. 曾经为BS-PCR最佳条件已经确定，对实验样品进行BS-PCR。
      1. 净化从剩余的引物和其他反应组分的扩增子。执行二氧化硅柱纯化，SPRI珠粒，或基于凝胶的清理和定量。
      2. 洗脱扩增子在30微升TE缓冲液中。
      3. 量化使用荧光检测使用制造商的协议扩增。
      4. 如果存在，每个样品的多个扩增子，它们集中在每体积重量相等在-20℃下量和存储。

4. NGS文库制备和图书馆QC

注：BSAS NGS文库制备采用了简化的转座子介导的协议与双索引（见材料清单上的文库制备试剂盒的选择详细信息）。这种方法提供了一个非常迅速，高通量路径用于文库构建已验证并显示了小至在亚硫酸氢盐测序应用^16,19没有偏压。双索引允许对样品比单一索引的较高多路复用。其他图书馆编制的风格可能是可能的，但还没有经过测试。

1. 准备从BS PCR扩增NGS库在96孔PCR板。稀释扩增至0.2纳克/μL作最后输入的1纳克，或5微升0.2纳克/微升稀释质量。执行清理使用SPRI珠生成库。用50-90微升SPRI珠隔离库瓦特第i此协议。的SPRI珠粒的量取决于文库的靶的尺寸和所需的测序反应循环数。
   注：使用微孔板摇床或板状涡旋重悬珠。
   1. 洗脱库脱珠转移到新的96孔板，并密封与储存在-20℃的热封膜。
2. 确定库（ 图4B和D）的大小和浓度。
   1. 运行库上的毛细管电泳的DNA根据制造商的协议上浆芯片。
      1. 使用高灵敏度的毛细管电泳的DNA检测，以确定检测到的库的峰的平均粒径和摩尔浓度的根据制造商的协议的估计。
   2. 量化库通过qPCR，使用被设计针对在生成的文库，并使用标准的具有已知浓度的接头序列的引物。
      1. 运行的qPCR重根据制造商的协议使用绝对定量的设置上的实时仪器操作。
      2. 使用库的大小和摩尔定量从定量PCR来计算摩尔浓度的库。
3. 稀释库利用10毫摩尔Tris-CL与0.1％Tween-20的pH为8.5 4纳米。商店库在96孔PCR板在-20℃下密封带热封膜。

5.测序文库使用台式序

1. 解冻并根据制造商的说明制备试剂。
2. 冰上解冻4纳米库。
3. 使在1.5ml微量离心管中加入1ml 0.2N的NaOH中。
4. 在等体积量池4nM的库，如果多个库被使用。
5. 稀释，并根据制造商的说明，以所期望的最终摩尔浓度变性汇集库。
   1. 不断稀释变性冰上，直到准备汇集库加载的试剂盒。
6. 生成样品表含样品名称和相应的索引信息，并指定只生成.fastq文件。
   1. 在序相应的样品表文件夹中加载样品板。
7. 添加共有600微升稀释变性库试剂盒对应得很好。
8. 运行的测序反应来完成。

6.甲基化定量数据分析

注：可用于NGS数据分析，包括商业和开源的多个程序。上运行的程序的细节可以与每个特定的包中找到。一般说明在随该软件包的特定命令的每一步都给出。 （见材料清单在这个协议中使用的软件的详细信息）。

1. 进口.fastq.gz文件到合适的本身quencing分析管道用于读修剪保留质量分数（Qscores）。 （在这个例子程序中选择“导入”，选择用来产生读测序方法，选择.fastq.gz读取导入文件，并保留Qscores从.fastq文件读通过取消在“常规”放弃质量得分“微调应用选项​​“，选择一个位置为导入读取，并选择”完成“。）
   1. 修剪并保留读取使用的NGS数据分析管道以下设置。 （对于此程序中选择“修剪序列”下“NGS核心工具”，然后选择读取进行修剪。）
      1. 只保留读取包含≥Q30分数。 （对于此程序中选择“裁剪使用质量得分”，并设置上限0.001）。
      2. 只保留读取包含≤1暧昧的核苷酸。 （对于此程序中选择“修剪暧昧核苷酸”，并设置1。'模棱两可的最大数量“选择一个位置来保存修剪读取，并选择”完成“。）
   2. 地图修剪读取从2.2 在硅片亚硫酸氢盐转化基因参考。 （对于示例软件选择“地图读取参考”下“NGS核心工具”，选择修剪读取被映射，并选择“下一步”，选择要转换的参考序列，并选择在“参考屏蔽'无遮拦”，选择“下一个'。）
      1. 分配最高刑罚为不匹配，插入和删除。设置参数，使得所读取的长度的100％被要求以90％的序列同一性映射。例如软件，分配的不匹配，插入和删除一个比分3。分配用于读出长度映射到参考的最低分数的分数为1，并指定为0.9的分数之间映射的读取和参照同一性的最低分数。
      3. 用于对准输出生成每个样本作为一个独立的文件中。例如软件，选择“创建独立的读取映射”下“输出选项”和“保存”“结果处理”项下。选择一个保存位置读取映射，并选择“完成”。
   3. 运行变量调用工作流上的映射读取。对于例如软件，选择“概率变检测”下“重测序分析”中，选择读映射进行分析，并选择“下一步”。选择“忽略非特异性匹配”下“读过滤器”，然后选择“下一步”。
      1. 确保序列处于DEP次方≥千X通过在'意义'设置'最小覆盖“到1000。
      2. 确保该变种呼叫存在于正向和反向的读取。例如软件，选择“需要在这两个前沿存在和反向看台上'下'变过滤器”，然后选择“下一步”。
      3. 出口对齐，变种，叫数据作为注解表。例如软件，选择“创建注释表”下“输出选项”，然后选择一个位置来保存变量表文件。选择“完成”。
2. 使用在从变体表文件感兴趣区域注明CpG部位的变体频率计算5-MC百分比。胞嘧啶在CG的基准C中的频率是5-MC百分比。
   注：使用变种，呼吁高度甲基化的CpG的算法，更换的CpG胞嘧啶与胸腺嘧啶。这将identiFY映射胞嘧啶的作为变体，导致胞嘧啶频率的CpG的。
3. 通过首先计算C的百分比不是在基准（C）的CpG部位确定亚硫酸氢盐转化效率（BSC）。乘非CpG的C'S由总T的映射（T _MP）的数量由全碳的映射在参考T的（C _MPT）的比例，以及鸿沟。 100减去计算出的​​频率是库（等式1）的亚硫酸氢盐转化效率（BSC）。



公式1亚硫酸氢盐转化效率。
注意：虽然哺乳动物基因组包含非常低量的非CpG的胞嘧啶甲基化（用干细胞²⁰的异常可能），植物基因组含有一个显著量。为了确定硫酸氢转换效率在这些referencê基因组，一个合适的途径是测序线粒体基因组或叶绿体基因组或已知的非甲基化的基因存在于植物的基因组的部分^21，并如上所述计算的转换效率。亚硫酸氢盐转化效率应该在大多数情况下≥98％，与未转化的C'S由于非CpG甲基化可能。

结果

BSAS读取正确对准到转换后的参考序列将类似于图5。CpG二核苷酸可以清楚地被识别和甲基化状态可以通过观察在所述映射中读取CpG部位碱呼叫来估计。例如，对于甲基化控制，0％甲基化将导致所有的读取映射到CpG部位含有T的（ 图5A）。 100％甲基化的控制，将导致所有的读取映射读取到CpG部位含有C的（ 图5B）。

胞嘧啶频率（C / c + T）在CpG位点的定量产生原始样品中的甲基化频。从全基因组的甲基化对照（n = 3 /甲基化比率）产生的代表性标准曲线示出了在各个甲基化率（ 图6）的甲基化的定量的非线性，以及在定量精度。作为该方法中的总的一个例子，RNA和DNA共 - 异从小鼠小脑和视网膜（N = 4 /组）迟来。视紫红质的表达，选择性地表达在视网膜组织，使用qPCR测量。表达仅在视网膜（ 图7A）进行检测。使用BSAS，CpG甲基化水平进行定量的视紫红质启动子区。相比，在视网膜（P <0.001，参数t检验）（ 图7B）<15％整个启动子区域的甲基化的累积水平> 80％，在小脑。从站点特定的甲基化定量BSAS甲基化定量分析的好处，并且甲基化水平可以在任何给定的基因组区域的CpG一个特定的基础上进行比较。相比，视网膜样品（P <0.001，参数t检验在各CpG位点）（ 图7C）横跨视紫红质启动子CpG甲基化水平显著高于小脑样本。

1highres.jpg“/>
图1：BSAS方法的示意图 。在该方法中，基因组DNA是亚硫酸氢盐转化成修改非甲基化胞嘧啶向尿嘧啶。接着，在PCR期间这些尿嘧啶变为胸腺嘧啶。 PCR扩增是针对感兴趣的区域，高度丰富为只是这些序列。将所得的PCR扩增子通过一个简单的过程tagmentation制成双索引库。随后文库进行测序在工作台上的下一代序器和测序读数映射到胞嘧啶的是在一个碱基特异性的方式来确定的，在硅片转换参考序列和百分比的甲基化。

图2： 在亚硫酸氢硅片的关心区域的转换的兴趣从任何基因组参考的任何区域可SELEC特德在硅片亚硫酸氢盐转化。非CpG的胞嘧啶替换为胸腺嘧啶的在参考序列。的CpG的胞嘧啶保持在亚硫酸氢盐转化的胞嘧啶参考。

图3：引物放置亚硫酸氢盐PCR引物被设计针对在硅片亚硫酸氢盐转化的参考序列。引物被设计成使得它们不重叠的CG二核苷酸或设计邻近的CG二核苷酸。

图4：高，质量差的PCR扩增和测序文库的例子代表电泳痕迹和凝胶显示（A）的转座子介导的库代扩增的理想选择用一个高点。浓度产物。从该扩增子产生的（B）中示出了库的高浓度，甚至尺寸分布。（C）的质量差的PCR扩增子可以是低的大小，浓度和包含多个产品和/或引物二聚体。这些质量低劣的扩增子将导致（D）中失败，文库产生具有低浓度和低粒度分布。

图5：从甲基化对照测序输出实例。 （A） 在硅片转换的参考序列与CpG位点以红色突出显示。 0％甲基化控制映射到读取选定区域显示胸腺嘧啶（贪婪亮点）映射在CpG位点。（B）100％的甲基化控制映射读取显示胞嘧啶的（蓝色高亮显示）映射在CpG位点。

图6：在一系列的甲基化甲基化控制定量全基因组的甲基化控制从0％混合在甲基化率（N = 3 /比）至100％用BSAS定量。绘制定量预计与量化表明甲基化控制定量的线性度。在所有的控制，有高精度的甲基化定量分析。

图7：从组织样品配对DNA甲基化和基因表达分析的实施例。 （A）中从小鼠小脑和视网膜组织视紫红质的相对mRNA表达。（B）的视紫红质的平均启动子甲基小脑和视网膜的DNA（*** P <0.001，参数t-检验）之间定量BSAS。 C）CpG位点特异的甲基化水平的跨越视紫红质启动子（-244至+6，相对于转录起始位点[TSS]）在小鼠小脑和视网膜的DNA。

讨论

BSAS允许集中，准确的DNA甲基化与定量两个样本数量和指标数量的高吞吐量的能力。此外，该库一代使用转座子介导的tagmentation需要较少的输入DNA，更快，包含了比传统的剪切和随后的适配器连接技术更少的步骤。这些改进在现有的方法允许任何感兴趣的目标进行精确基层DNA甲基化分析。此外，BSAS为确认已使用的全基因组亚硫酸氢钠²²或²³捕捉方法确定的感兴趣区域（差异甲基化区域（的DMR），CpG ​​岛，监管区）的新方法。采用正交确认的方法和更精确的定量方法，提高了表观基因研究的严格分析。实现了与BSAS高阅读深度允许在一个窄的置信区间准确定量，再质保的义务在精确的定量分析^16。此外，为了在单次实验运行许多样品的能力可以增加对整个基因组的方法确认，这将增加功率统计样本大小;目前许多后生研究的一个弱点。重要的是，BSAS提供一个碱基特异性CpG甲基化的定量，其允许可检验的假设为后生研究的产生。假设如在特定的反应元件增加甲基化将导致一个特定的转录因子降低的结合。 BSAS，结合配对mRNA表达数据，提供了用于在一个单一的一块组织或细胞群的获得两种表观遗传调节和表达的方法。调控和表达的配对测量，也增加了科学严谨和结果的影响。

在BSAS协议中的关键步骤包括引物设计，扩增优化和生成库/QC。 PCR偏见介绍，如果引不正确的设计和放大，在不同的效率^8。使用甲基化的标准，如酶促产生100％和0％的胞嘧啶甲基化的标准，可以用于确定程度，如果有的话，甲基化定量的偏差，并且是一个适当的方法控制量化甲基化¹⁶时。类似地，应注意优化PCR反应，用于产生单个高质量扩增子时服用。如果没有扩增子是本使用在协议详述所建议的准则，优化步骤包括增加的PCR循环数或bsDNA输入量。如果有多个PCR产物，包括引物二聚体，优化步骤包括减小bsDNA输入量，降低了PCR引物的摩尔浓度的输入，提高退火温度，或降低PCR循环数。这些优化过程的一个或一个组合产生足够的结果，以产生一个单一的集锦。H-质量的扩增子。另外，重新设计引物对引物不产生一个扩增子的好办法。对于这些地区的地方重新设计不适合或没有可能，退化引物包括胞嘧啶的和胸腺嘧啶的在正向引物CpG位点和腺嘌呤的和鸟嘌呤的在反向引CpG位点可能工作。然而，这些引物组还需要进一步优化，以确保无PCR的偏差发生，这个协议的范围之外的处理。质量库是成功的非常重要的。确保QC措施进行了测序前确定库的规模，质量和数量。测序反应是容易发生故障与低质量的库的输入。偏压中的甲基化的定量可以通过确保甲基化频率存在于正向和反向测序读数来缓解。此外，基地调整至CpG位点的质量分数应≥Q30高协nfidence的定量。虽然其他人已经表明以前很少在亚硫酸氢盐转化的DNA ¹⁹ tagmentation反应没有偏差，保证上述指标可以确认低偏置甲基化定量。

BSAS目前在CpG位点的甲基化测量，但是，这种方法的未来的扩展将允许的CpG 5- HMC的成对测量用加入了如引入钾高钌酸铵重亚硫酸盐转化²⁴之前实验步骤5-MC和5- HMC区分。这将增加BSAS效用的影响通过允许配对甲基化，既5-MC和5- HMC，和表达数据，以确定DNA甲基化的基因表达的调节作用。转座子介导的文库制备的PCR扩增子并导致降低测序深度在扩增区域的端部。因此，高的兴趣的区域的设计应该驻留在扩增子的中间。怎么样以往，因为BSAS生成测序阅读深度> 1000的X，5-MC定量精度邻近扩增区域的端部测量仍然很高。

材料

本文使用的材料清单

姓名	公司	目录编号	评论
安捷伦 2100 生物分析仪	安捷伦	G29393AA
安捷伦 RNA 6000 纳米试剂盒	安捷伦	5067-1511
安捷伦高灵敏度 DNA 试剂盒	安捷伦	5067-4626
台风 9200	Amersham/GE Healthcare 生命科学
Agencourt AMPure XP - PCR 纯化	贝克曼库尔特	A63880	5 ml
PowerPac基本电源Bio	Rad	164-5050
twin.tech PCR板96	Eppendorf	951020362	半裙边
热封膜	Eppendorf	0030127838
热封箔	Eppendorf	0030127854
热封机	Eppendorf	5390000.024
0.1-10 μl 吸头	Eppendorf	0030073.002
2-200 μl 吸头	Eppendorf	0030073.045
50-1,000 μl 吸头	Eppendorf	0030073.100
MixMate	Eppendorf	5353000.014
1.5 ml 微量离心管	Eppendorf	0030121.023
离心机 5424	Eppendorf	5424000.010
2.0 ml 微量离心管	Eppendorf	022363352
Nextera XT DNA 样品制备试剂盒	Illumina	FC-131-1024	24 rxn
Nextera XT Index Kit	Illumina	FC-131-1001	24 indexes
MiSeq Desktop Sequencer	Illumina
MiSeq 试剂盒 v2	Illumina	MS-102-2002	300 次循环
KAPA SYBR FAST 通用 qPCR 试剂盒	KAPA Biosystems	KK4824
Quant-iT PicoGreen dsDNA 检测试剂盒	Life Technologies	P7589
Quant-iT RiboGreen RNA 检测试剂盒	Life Technologies	R11490
ProFlex 96 孔 PCR 系统	Life Technologies	4484075
Novex 6% TBE DNA 凝胶	Life Technologies	EC6261BOX	10 孔
Novex TBE 电泳缓冲液	Life Technologies	LC6675	5x （1 L）
Novex 高密度 TBE 样品缓冲液	Life Technologies	LC6678	5x （10 ml）
1 Kb Plus DNA 分子量标准	Life Technologies	10787-018
Xcell SureLock Mini-Cell	Life Technologies	EI0001
超纯 10 mg/ml 溴化乙锭	Life Technologies	15585-011
7900HT 快速实时 PCR 系统	Life Technologies	4329001	384 孔模块
DynaMag-96 侧裙边	Life Technologies	12027
SpectroMax M2 多模式微孔板检测仪	Molecular Devices/VWR	89429-532
AllPrep DNA/RNA Mini Kit	Qiagen	80204
不锈钢珠	Qiagen	69989	5 mm
CLC 基因组学工作台	Qiagen		甲基化管道软件
QIAQuick PCR 纯化试剂盒	Qiagen	28104	50 rxn
TissueLyser II	Retsch/Qiagen	85300
无核酸酶水	Sigma	W4502
TRIS 盐酸盐	Sigma	PHG0002-100G
Tween-20	Sigma	P9416-50ML
NanoDrop 2000	Thermo Scientific
384 孔全裙边板，标准	Thermo Scientific	AB-1384
绝对 qPCR 密封	Thermo Scientific	AB-1170
EZ DNA 甲基化-Lightning 试剂盒	Zymo Research	D5030	50 rxn
ZymoTaq DNA 聚合酶	Zymo Research	E2001	50 rxn