材料和试剂

移液器吸头（TipOne，STARLAB，产品目录号：S1161-1800，S1182-1830和S1181-3810）
0.65毫升生物奥管微管（Diagenode，目录号：C30010011）
0.5ml和2ml DNA LoBind管（Eppendorf，目录号：0030108035和0030108078）
Micro Bio-Spin 6柱（Bio-Rad Laboratories，目录号：7326221）
1.5 ml Lobind管（Eppendorf，目录号：0030108051）
2 ml Lobind管（Eppendorf，目录号：0030108078）
DNeasy Blood＆amp;组织试剂盒（QIAGEN，目录号：69504）
100-bp DNA标记（Thermo Fisher Scientific，Invitrogen TM，目录号：15628019）
E-gel EX琼脂糖凝胶2％（Thermo Fisher Scientific，Invitrogen TM，目录号：G401002）
β-葡糖基转移酶（β-GT）和相关的反应缓冲液（New England Biolabs，目录号：M0357S）
DBCO-PEG4-生物素（Sigma-Aldrich，目录号：760749）
UDP-6-N3-Glc（Active Motif，目录号：55020）
DMSO（Sigma-Aldrich，目录号：D8418）
QIAquick核苷酸去除试剂盒（QIAGEN，目录号：28304）
Dynabeads M-280链霉亲和素（Thermo Fisher Scientific，Invitrogen TM，产品目录号：11205D）
NEBuffer 2（New England Biolabs，目录号：B7002S）
10x NEBuffer 4（New England Biolabs，目录号：M0357S）
ATP（10mM）（New England Biolabs，目录号：P0756S）
dNTP溶液混合物（New England Biolabs，目录号：N0447S）
T4 DNA聚合酶（New England Biolabs，目录号：M0203S）
DNA聚合酶I，Large（Klenow）片段（New England Biolabs，目录号：M0210S）
T4多核苷酸激酶（New England Biolabs，目录号：M0201S）
T4 DNA连接酶高浓度（New England Biolabs，目录号：M0202T）
无核酸酶的水（Thermo Fisher Scientific，Invitrogen TM，目录号：AM9937）
SCL-exo P7衔接子2个寡核苷酸的退火（5'Phos =磷酸化5'末端）：
P7 exo适配器反向：5'Phos-GTGACTGGAGTTCAGACGTGTGCTCTTCCGATC-OH 3'
P7 exo-接头正向：5'-OH-GATCGGAAGAGCACACGTCT-OH 3'
Phi29聚合酶（New England Biolabs，目录号：M0269S）
λ外切核酸酶（New England Biolabs，目录号：M0262S）
RecJ外切核酸酶（New England Biolabs，目录号：M0264S）
糖原（5mg / ml）（Thermo Fisher Scientific，Invitrogen TM，目录号：AM9510）
氯化钠（NaCl）（Acros Organics，目录号：AC207790050）
EtOH（100％）（VWR，目录号：20821.310）
SCL-exo P7引物：
5'OH-GACTGGAGTTCAGACGTGTGCT-OH 3'
Agencourt AMPure XP（Beckman Coulter，目录号：A63880）
Qiagen MinElute PCR纯化试剂盒（QIAGEN，目录号：28004）
SCL-exo P5接头：2个寡核苷酸的退火：
P5 exo适配器反向：5'OH-AGATCGGAAGAGCG-OH 3'
P5 exo-接头正向：5'-OH-TACACTCTTTCCCTACACGACGCTCTTCCGATCT-OH 3'
NEBNext高保真2x PCR主混合物（New England Biolabs，目录号：M0541S）
SCL-exo通用P5 PCR引物（* =硫代磷酸酯S-键）：
5'-OH-AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACG * A-OH 3'
SCL-exo指数P7 PCR引物（* =硫代磷酸酯S键）（指数序列来自TruSeq LT）：
索引2：
5'-OH-CAAGCAGAAGACGGCATACGAGAT ACATCG GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引4：
5'-OH-CAAGCAGAAGACGGCATACGAGAT TGGTCA GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引5：
5'-OH-CAAGCAGAAGACGGCATACGAGATCACTGT GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引6：
5'-OH-CAAGCAGAAGACGGCATACGAGAT ATTGGC GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引7：
5'-OH-CAAGCAGAAGACGGCATACGAGATGATCTG GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引12：
5'-OH-CAAGCAGAAGACGGCATACGAGAT TACAAG GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引13：
5'-OH-CAAGCAGAAGACGGCATACGAGAT TTGACT GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引14：
5'-OH-CAAGCAGAAGACGGCATACGAGATGGAACT GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引15：
5'-OH-CAAGCAGAAGACGGCATACGAGAT TGACAT GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引16：
5'-OH-CAAGCAGAAGACGGCATACGAGAT GGACGG GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引18：
5'-OH-CAAGCAGAAGACGGCATACGAGAT GCGGAC GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
索引19：
5'-OH-CAAGCAGAAGACGGCATACGAGAT TTTCAC GTGACTGGAGTTCAGACGTGTGC * T-OH 3'
注意（关于寡核苷酸）：
1. 所有寡核苷酸均由Sigma-Aldrich生产，通过HLPC纯化并以100μM终浓度重新悬浮于水中。
2. SCL-exo P7适配器和 SCL-exo P5适配器通过将4对体积的退火缓冲液（参见配方）中的互补寡核苷酸对混合并通过在95℃下加热5分钟退火，然后让其缓慢冷却至室温来获得。
3. 设计用于SCL-exo的寡核苷酸来自Illumina©2007-2012 Illumina，Inc.的P5和P7寡核苷酸序列。保留所有权利。只有Illumina客户创建的衍生作品才被授权与Illumina仪器和产品一起使用。所有其他用途都是严格禁止的。
安捷伦高灵敏度DNA试剂盒（安捷伦科技，产品目录号：5067-4626）
Qubit dsDNA HS分析试剂盒（Thermo Fisher Scientific，Invitrogen TM，目录号：Q32854）
EDTA（500mM，pH8.0）（AppliChem，目录号：A4892,0500）
HEPES（1M）（Gibco TM，目录号：15630056）
脱氧胆酸钠（Sigma-Aldrich，目录号：D6750）
NP-40，IGEPAL CA-630（Sigma-Aldrich，目录号：I8896）
氯化锂（LiCl）（Sigma-Aldrich，目录号：62476）
氯化镁六水合物（MgCl 2•6H 2 O）（Merck，目录号：442611）
硫酸铵（（NH 4）2 SO 4）（Merck，目录号：101217）
DTT（Sigma-Aldrich，目录号：D9779）
Tris（MP Biomedicals，目录号：04819638）
盐酸（HCl）（Sigma-Aldrich，目录号：H9892）
分子生物学甲酰胺（Sigma-Aldrich，目录号：F9037）
1x PBS（Fisher Scientific，目录号：BP399）
退火缓冲区（请参阅食谱）
RIPA缓冲区（请参阅食谱）
尼克修复缓冲液低DTT 10倍（见食谱）
TE缓冲液（见食谱）
洗脱缓冲液（见食谱）
装订＆amp;洗涤（B＆amp; W）缓冲液（见食谱）

设备

PIPETMAN经典TM吸管（Gilson，目录号：F123600，F144801，F123602和F123615）
Bioruptor Pico带水冷却器（Diagenode，产品目录号：B01060001和B02010003）
E-gel Power Snap电泳装置（Thermo Fisher Scientific，InvitrogenTM，产品目录号：G8100）
Qubit 3荧光计（Thermo Fisher Scientific，Invitrogen TM，目录号：Q33216）
冷冻离心机（Eppendorf，型号：5424 R）
热循环仪ProFlex PCR系统（Thermo Fisher Scientific，Applied Biosystems TM，目录号：4484073）
ThermoMixer C和Eppendorf ThermoTop（Eppendorf，产品目录号：5382000015和5308000003）
DynaMag-2 Magnet（Thermo Fisher Scientific，目录号：12321D）
Speed-Vac Savant（Thermo Fisher Scientific，目录号：DNA120-115）
2100生物分析仪（Agilent Technologies，型号：2100，目录号：G2939BA）
小型离心机（Bio-Rad Laboratories，目录号：1660603）

程序

使用QIAGEN DNeasy试剂盒提取基因组DNA，并通过超声破碎成300bp片段。酶β-葡糖基转移酶催化叠加葡萄糖到存在于gDNA片段中的5hmC。然后叠氮化物与生物素缀合物反应，允许将修饰的DNA固定在链霉亲和素包被的磁珠上（图1A）。在修复完毕后，将Illumina P7衔接子连接和切口修复，将捕获的DNA与5'→3'外切核酸酶λ和RecJf孵育。 λ外切核酸酶消化双链DNA的一条链，并在遇到珠子结合的生物素化5hmC时停止，而RecJ外切酶消化可能由未修饰污染物DNA消化产生的单链DNA λ外切核酸酶。从珠上洗脱后，DNA变性为单链DNA分子。接着是第二链合成，Illumina P5衔接子的连接，PCR扩增和Illumina测序。单末端测序从P5衔接头开始，鉴定λ外切核酸酶停止消化的位置及其最近的羟甲基化CpG（图1B）。

图1. SCL-exo过程概述。：一种。作为gDNA化学修饰的第一步，β-葡糖基转移酶催化叠氮化物葡萄糖从UDP-6-N3-Glc转移至5hmC。然后使用点击化学将生物素缀合物（DBCO-PEG4-生物素）添加至N3-Glc修饰的5hmC。 B. SCL-exo协议流程图。

为Illumina测序准备样品
我们强烈建议使用无RNA基因组DNA（gDNA）用于SCL-exo方案。我们通过使用QIAGEN DNeasy试剂盒并按照制造商的方案中所述添加RNaseA消化步骤来纯化目标RNA的无gDNA。来自任何类型的组织或培养细胞的无RNA的gDNA可以用于SCL-exo方案。但是，人们应该记住，组织间5hmC的总量差异很大，因此SCL-exo所需的gDNA的起始量可能因样本来源而异。当处理来自不同测试条件的样品时，我们强烈建议在超声处理后在每个样品中加入相同量的羟甲基化DNA标准物。覆盖此标准的读数数量可用于标准化样品之间的SCL-exo信号。
1. 使用Bioruptor Pico在0.65ml超声管中在10μl10mM Tris，pH 8中超声处理1μg感兴趣的gDNA以获得约300bp的DNA片段。超声波周期应设置为30秒/ 30秒。为确保正确和可重复的超声处理，我们建议进行3次超声处理循环，然后进行短时离心，然后再进行3次超声波处理，然后短时间离心，最后进行4次超声处理。
2. 通过在E-gel EX琼脂糖凝胶（2％）中运行100-bp DNA标记和0.5μl超声处理的gDNA（在19.5μl水中稀释）10分钟，可以快速检查超声处理效率。你应该获得约250-300 bp的DNA片段。
注意：程序步骤A3到A12来自我们的同事，稍作修改（ Szulwach et al。，2012 ，Bio-protocol）。
1. 将剩余的9.5μl超声处理的DNA与：2μl10x NEBβ-GT反应缓冲液（随β-GT酶提供）+0.68μlUDP-6-N3-Glc（3mM）+1μlNEB Beta- GT酶+ 6.8μl水。
2. 通过移液混合并在37℃的热循环仪中孵育1小时（无加热盖）。
3. 用微型离心机快速离心（5分钟，2000 em x em）。
4. 通过10倍稀释在DMSO中的30mM储备溶液制备DBCO-PEG4-生物素缀合物在DMSO中的3mM工作溶液。在-20°C储存。
5. 向来自步骤A5的DNA样品中加入1μlDBCO-PEG4-生物素缀合物工作溶液以达到150μM的终浓度。
6. 通过移液混合并在37℃的热循环仪中孵育1小时（无加热盖）。
7. 用微型离心机快速离心（5分钟，2,000 em x g），并用QIAquick Nucleotide Removal Kit清除反应。
  每个柱子至少要加30μl水注意：生物素化的DNA样品可以在-20°C下保存几天。
8. 每次用25μlDynabeads M-280链霉抗生物素洗涤3次，每次用100μl1×Binding＆amp;在0.5ml Lobind管中洗涤（B＆amp; W）缓冲液（见配方）。用磁力支架将珠子从缓冲液中分离出来，并将珠子重悬于30μl2×B＆amp; W缓冲液和140μl1×B＆amp; W缓冲液中。
9. 将30μlDNA洗脱物（来自步骤A9）添加到来自之前步骤的重悬的珠子中。 B＆amp; W缓冲液的最终浓度应为1倍。
10. 在室温下孵育30分钟注意：在此步骤中准备步骤A15的混合物。
11. 转移到2毫升Lobind试管中，用1毫升1x B＆amp; W缓冲液使用磁力架清洗珠子5次。
12. 用1 ml 10 mM Tris-HCl pH 8清洗2次。不要让珠子干燥。
13. 然后珠经历5次连续反应（在2ml Lobind管中，以900rpm在热混合器中搅拌）如下：
  末端修复：准备含10μlNEB2缓冲液（10x），10μlATP（10mM），1μldNTP（10mM），5μlT4 DNA聚合酶（3U /μl），1μl DNA聚合酶I大Klenow片段（5U /μl），5μlT4多核苷酸激酶（T4PNK）（10U /μl）和68μl无核酸酶的水。
  将混合物加入2毫升Lobind中的珠子中。
  在30°C孵育30分钟，在900转/分钟的搅拌器中搅拌
14. 用1ml RIPA缓冲液（参见配方）洗涤2次，用10mM Tris-HCl（pH 8）洗涤2次。除去最后的Tris洗涤液后，用微型离心机快速离心（5分钟，2,000 emg ）并将管放回磁性支架。去除Tris的痕迹。确保您对步骤A18，A20，A22和A24执行相同的操作。不要让珠子干燥。
15. P7适配器的结扎准备含10μlNEB2缓冲液（10x），10μlATP（10mM），15μlSCL-exo P7接头（10μM），1μlT4 DNA连接酶（2,000μg） U /μl）和65μl无核酸酶的水。将混合物加入2毫升Lobind管中的珠子中。
  。在25°C孵育1小时，在搅拌机中以900转/分钟的速度搅拌。
16. 用1ml RIPA缓冲液洗涤两次，用1ml 10mM Tris-HCl（pH 8）洗两次。
17. 尼克修复：
  准备一个混合物，包含1.5微升Phi29聚合酶（10 U /μl），10微升自制尼克修复低DTT缓冲液（10倍）（见食谱），1.5微升dNTP（10毫米）和87微升无核酸酶水。将混合物加入2毫升Lobind管中的珠子中。
  。在30°C孵育20分钟，在900转/分的搅拌器中搅拌
18. 用1ml RIPA缓冲液洗涤两次，用1ml 10mM Tris-HCl（pH 8）洗两次。
19. λ外切核酸酶消化：
  准备一个混合物，包含2μLλ外切核酸酶（5U /μl），10μLNEBλ外切核酸酶缓冲液（10x）和88μL无核酸酶水。将混合物加入2毫升Lobind管中的珠子中。
  。在37°C孵育30分钟，在900转/分钟的搅拌下在热混合器中。
20. 用1ml RIPA缓冲液洗涤两次，用1ml 10mM Tris-HCl（pH 8）洗两次。
21. RecJ外切核酸酶消化：
  准备包含1微升RecJ外切核酸酶（30 U /μl），10微升NEB2缓冲液（10x）和89微升无核酸酶水的混合物。将混合物加入2毫升Lobind管中的珠子中。
  。在37°C孵育30分钟，在900转/分钟的搅拌下在热混合器中。
22. 用1ml RIPA缓冲液洗涤两次，用1ml 10mM Tris-HCl（pH 8）洗两次。
23. 洗脱：
  在100μl洗脱缓冲液（参见配方）中将珠子在90°C下孵育5分钟，然后直接放在冰上冷却样品。
24. 将100μl洗脱液转移至新的1.5 ml Lobind试管中，加入300μl10 mM Tris-HCl，pH 8。
25. DNA沉淀：
  1. 加2μl糖原，16μlNaCl（5M）并充分混合。加入800μl100％乙醇并充分混合。
  2. 在-80°C孵育管至少30分钟（如果可能，过夜）。
  3. 在4℃以20000×gg离心30分钟。
  4. 小心取出上清液，不要扰动沉淀。
  5. 加入500μl70％EtOH。
  6. 在4℃下以20000×g克离心5分钟。
  7. 。清除上清液。
  8. 加入500μl的100％乙醇。
  9. 在4℃下以20000×g 离心5分钟
  10. 。清除上清液。
  11. 在45℃的Speed-Vac中干燥颗粒10-20分钟并重悬于20μl10mM Tris-HCl（pH 8）中。
  12. 纯化的DNA样品可以在-20℃保存一晚。转到步骤A28。
26. DNA变性：
  1. 将20μlDNA溶液转移至PCR管中，并在热循环仪中于95℃孵育DNA样品5分钟。
  2. 然后将管直接放在冰上以冷却样品。
27. 第二链合成：
  1. 向含有20ul DNA溶液的试管中加入下列试剂：20μl不含核酸酶的水，5μlSCL-exo P7引物（1μM）和5μlNEB Phi29反应物缓冲区（10x）。轻轻混合。
  2. 在热循环仪中，将样品在65℃孵育5分钟，然后在30℃孵育2分钟。暂停PCR程序。
  3. 立即加1μlPhi29聚合酶（10U /μl）和1μldNTP（10mM），轻轻混合。
  4. 重新启动PCR程序，并将样品在30℃的热循环仪中孵育20分钟，然后在65℃孵育10分钟。
28. DNA纯化：
  1. 向52μl样品中加入52μl室温Ampure珠（1体积）。
  2. 在室温下孵育15分钟。
  3. 将试管放在磁力架上，小心取出上清液。将试管放置在磁力架上，用新鲜制成的80％乙醇洗涤珠子两次（加入第一次乙醇洗涤后至少等待30秒）。
  4. 用微型离心机离心（2000g×g 5秒），将管放回磁力架上，除去剩余的乙醇。
  5. 将管子在磁力架上打开并让其干燥10-15分钟。
  6. 加入22μl室温10 mM Tris-HCl，pH 8，从磁力架上取下试管并充分混合。确保所有的珠子重新悬浮并潮湿。
  7. 从磁力架上取下试管并在室温下孵育3分钟。
  8. 将珠子放回磁力架，一旦它们充分包装，小心吸取20μlDNA洗脱液，并将其放入新的PCR管中。
29. SCL-exo P5接头的连接：
  1. 在PCR管中，将以下试剂加入到20μlDNA样品中：22.5μl无核酸酶水，1.5μlSCL-exo P5衔接子（10μM），5μlNEB T4 DNA连接酶缓冲液（10x）和1μlT4 DNA连接酶（2,000U /μl）。轻轻混合。
  2. 在热循环仪中，在25°C孵育60分钟，然后在65°C孵育10分钟。
30. DNA纯化：
  向50μl样品中加入50μl室温Ampure珠（1体积），并按照步骤A30进行。所得到的20μl洗脱DNA溶液用于最终PCR。
31. PCR扩增：
  1. 在PCR管中，制备混合物，其中含有4μl不含核酸酶的水，25μlNEBNext高保真PCR主混合物（2x），0.5μlSCL-exo通用P5 PCR引物（25 μM）和0.5μlSCL-exo索引P7 PCR引物（25μM）（选择您感兴趣的指标）。加入20μlDNA样品并轻轻混合。
  2. 将管放入热循环仪并运行以下程序：
    98℃30秒
    然后18个循环：98℃10秒，65℃30秒，72℃30秒
    72°C 5分钟
    4°C永远
32. DNA纯化：
  向50μlPCR样品中加入50μl室温Ampure珠（1体积），并按照步骤A30进行。你应该得到一个20μL的SCL-exo库。
33. 使用Qubit和dsDNA高灵敏度试剂盒测量DNA浓度。使用Agilent BioAnalyzer检查库质量（请参见图2）。如果有接头或引物污染，建议重新进行Ampure纯化（1体积的珠子用于1体积的DNA）。将库合并为复用。获得足够的索引复杂性，以便索引排序成功。如果您有任何疑问，请联系您的测序设备。
34. 提交Illumina单端测序MiSeq / GAII / HiSeq到测序设备。
  
  图2. SCL-exo文库的质量控制A. SCL和SCL-exo文库的琼脂糖凝胶电泳。通过省略外切核酸酶消化步骤获得SCL文库。请注意，来自SCL文库的DNA片段平均比SCL-exo文库长100 bp。 B.一组SCL-exo文库的BioAnalyzer电泳图谱。按照制造商的方案，在Agilent高灵敏度DNA芯片上运行1μl的SCL-exo文库。 DNA文库长度应在200到400 bp之间。注意不存在接头二聚体峰（大约120bp）和缺少PCR引物（大约50bp）是很重要的。如果这些污染物存在，我们建议按照步骤A30重新进行Ampure纯化（1体积的Ampure珠粒用于1体积的DNA文库）。
生物信息学鉴定SCL-exo fastq文件中的羟甲基化CpGs
我们构思并实施了一种生物信息学方案，以通过测序平台一式三份地从SCL-exo fastq文件中鉴定羟甲基化的CpG。该协议涉及以下步骤：
1）使用程序 SolexaQA （Cox et。，2010）对序列读数进行修剪和过滤。
2）使用程序Bowtie （Langmead et al。，2009）分别将高质量读数映射到基因组的每条链上，以便为正向和反向股线。山姆文件是文本文件，包含序列读取及其相关的基因组定位（如果有的话），并且可以解析以识别映射基因组上独特位置的读段。
3）使用我们的python程序 generate-SCL-exo signal-from-sams 直接读取sam文件中的序列，为每个重复项创建一个羟甲基化的CpG信号（假发）文件。该程序计算唯一重叠任何给定CpG的读取次数，并将这些值以CpG或基准对分辨率存储到信号（假发）文件中。使用基因组浏览器可以将假发文件可视化，例如 IGB （Nicol et al。，2009）。
注意：我们所有的python程序都可以在以下网址找到： https：// mycore。 core-cloud.net/index.php/s/4gyZ9dLTqgo86dt 。
4）使用python程序 generate-SCL-exo consensus-signal ，通过检索存在于三个重复中的至少两个中的共有CpG二核苷酸来鉴别推定的羟甲基化胞嘧啶。
5）使用峰值调用算法（ generate-SCL-exo peaks ）确定在5hmC中显着富集的CpG二核苷酸集合，并且具有良好选择的阈值。
每个步骤的详情如下：
1. 修剪和过滤顺序读取
  对于羟甲基化的CpGs，应该只保留高质量的读数。因此，我们使用SolexaQA程序（Cox et。，，2010）修剪和过滤SCL-exo fastq文件中的读取内容。该程序需要两个参数：质量阈值和最小长度。首先，从读数中去除质量低于质量阈值的所有测序核苷酸。其次，短于最小长度的读取被删除。我们使用值20作为最小核苷酸测序质量，对应于10 -2（或任何给定核苷酸上发生测序错误的概率为1％）的p值和17作为最小读数长度。修剪是通过进入SolexaQA目录并在Linux下键入来实现的：
  
  perl DynamicTrim.pl fastq -h quality -d。
  
  其中是fastq文件的路径和文件名， quality 是质量值（，例如，20）。这将生成一个修剪过的fastq文件 fastq.trimmed 。然后通过输入：
  来实现过滤
  perl LengthSort.pl fastq.trimmed -l < minlength
  
  其中 minlength 是保留读取的最小长度（，例如，17）。
2. 将过滤读数映射到基因组的两条链上
  可以使用Bowtie（Langmead et al。，2009）分别将保留的高质量读数映射到基因组的正向和反向链上，具有以下参数：
  
  -m 1 --sam --strata - 最好的--norc [或--nofw]
  
  ：指定可用于映射过程的计算机处理器的数量;
  长度：将读取长度考虑在内以将读取映射到基因组上;
  nb _ 不匹配：允许的不匹配数量;
  -m 1表示我们只保留将基因组映射到独特位置的读段;
  --norc（分别是--nofw）基因组反向链（或正向链）不用于映射。
  
  请注意，最初需要对基因组fasta文件进行索引（请参阅 Bowtie用户指南 / a>）。
  阅读必须分别映射到正向和反向链上，为每个链产生一个sam文件。使用Linux命令与 Bowtie 启动映射：
  
  ./bowtie -p 处理器 --best -l 28 -n 2 -m 1 --sam --strata --norc基因组fastq ＆gt; fw - sam
  
  将读取文件 fastq 映射到索引基因组文件的正向链（通常为.hs文件），以生成正向链fw-sam文件，以及： ./bowtie -p processors --best -l 28 -n 2 -m 1 --sam --strata --nofw genome fastq＆gt;用于将读取文件fastq映射到索引基因组文件的反向链上，从而生成反向链rv-sam文件。
3. 解析单链sam文件以CpG或基准对分辨率生成假发文件
  1. 首先根据5-羟甲基胞嘧啶的存在捕获DNA片段，然后用外切核酸酶修剪。由于胞嘧啶在CpG背景中大部分是羟甲基化的，因此5hmC阳性读数应该在每个序列开始的几个核苷酸内富集于CpG中（Sérandour et al。，2016）。如果有的话，SAM文件包含序列读取及其相关的基因组定位。我们的python程序（ generate-SCL-exo signal-from-sams ）解析了正向和反向链接的sam文件，并依次考虑了唯一映射到基因组上的所有读取。它检查每个本地化读取是否在其序列的前几个核苷酸内呈现CpG。通常，位于读数开始处的10个碱基对的长窗口用于证明CpG的存在。读取不显示窗口内的任何CpG将被丢弃。在窗口内显示两个或多个CpG的读数被保留在一边（它们的CpG将被存储在不同的文件中），因为不可能确定地确定哪个CpG被羟甲基化。
  2. 当在窗口内找到一个CpG时，就可以确定其精确的基因组坐标（从sam文件提供的读取位置和读取内的CpG位置），并存储在内存中的哈希表类型结构中。第一次遇到CpG位置时，值为1与基因组位置相关联。如果CpG位置已经包含一个值，那么该值将增加1，从而有效地存储覆盖该特定位置的读取次数。请注意，程序计算与sam文件关联的链：在反向链sam文件中，必须从右向左读取读取序列，而CpG坐标必须进行调整（程序会自动解释这个坐标）。一旦所有读数都被解析，散列表默认以CpG分辨率存储在假发文件中：每个CpG的基因组位置与其重叠读数的数目相关联。默认情况下，我们的程序将叠加在两条链上发现的给定CpG的重叠数相加。可以使用该程序，以便它为CpG的每个链胞嘧啶产生不同的信号值，将每个胞嘧啶位置与每个链重叠的读数数目相关联，从而以碱基对分辨率产生假发文件。 >
  3. 我们的程序使用以下Linux命令启动：
    
    python generate-SCL-exo.py signal-from-sams fw-sam rv-sam window SCL-exo-wig [resolution]
    
    fw-sam 和 rv-sam 分别指定正向和反向链接的sam文件的文件名;
    window 代表用于识别羟甲基化CpG二核苷酸的窗口的长度（通常为10），以碱基对表示;
    分辨率是一个可选参数，取值为2或1，具体取决于是否将以CpG或基准对分辨率生成 SCL-exo-wig 文件（默认如果参数未指定，则为2）。
4. 确定在三次重复中至少有两次发现的共识CpGs
  我们提供了一个比较三个SCL-exo信号文件的Python程序 generate-SCL-exo consensus-signal ，并返回一个假发文件，其中包含至少两个三次重复中鉴定的羟甲基化CpG，连同它们平均信号。在三个文件中的至少两个文件中，识别出的CpG位置的值必须大于最小阈值 min-threshold 。这种共识信号的例子可以在图3中找到。
  该功能可以这样调用：
  
  python generate-SCL-exo.py consensus-signal SCL-exo-wig1 SCL-exo-wig2 SCL-exo-wig3 min-threshold
  
  图3.来自P19胚胎癌细胞的mm8 chr11区域中的SCL-exo信号的综合基因组浏览器视图。为了比较，显示了Input-seq（P19细胞的基因组DNA）和SCL-seq（无外切核酸酶步骤）。
5. 确定SCL-exo假发文件中显着富集的羟甲基化CpGs
  使用峰值调用算法来确定在SCL-exo信号文件内显着富集5hmC的CpG。该算法在SCL-exo假发文件中查找相邻基因组位置，该位置显示出两个CpG坐标高于预定义阈值的信号值（Sérandour et。，2016）。
  该程序既可以用于共识SCL-exo信号文件，也可以用于SCL-exo假发单独复制。
  Python程序 generate-SCL-exo 峰会将一个 SCL-exo-wig 文件和预定义的阈值结合起来，所有符合上述限制的CpG头寸。
  使用Linux命令启动CpG或基准对分辨率下的 SCL-exo-wig 文件上的峰值调用：
  
  python generate-SCL-exo.py peaks SCL-exo-wig threshold 。

数据分析

有关数据处理和分析的信息可以在原始研究文章中找到：
https://doi.org/10.1186/s13059-016-0919-y 。

笔记

有关步骤A15至A24的注意事项：
1. 10mM Tris-HCl（pH 8）中磁珠的磁性更小。在去除第二次Tris洗涤过程中将管保持在磁力架上，以避免珠损失。然后短暂旋转试管，将其放回磁力架上并除去残留的Tris缓冲液。应该小心地进行Tris洗涤以消除任何可能对随后的酶促反应有害的洗涤剂痕迹。
2. 不要让链霉亲和素珠干透。
  在洗涤之前几分钟准备酶混合物。
关于步骤A30，A32和A34中的安瓿珠纯化：
请注意，任何剩余的EtOH痕迹都会抑制下一个酶促反应。

食谱

退火缓冲区
10mM Tris，pH 8
50 mM NaCl
1 mM EDTA
RIPA缓冲区
50 mM HEPES，pH 7.6
1 mM EDTA
0.7％脱氧胆酸钠
1％NP-40
0.5 M LiCl
尼克修复缓冲区低DTT（10x）
100mM MgCl 2 2/2 500mM Tris-HCl，pH7.5
100 mM（NH4）2 SO4 4 /
10 mM DTT
TE缓冲液（pH 7.4）
10 mM Tris
1 mM EDTA
HCl
洗脱缓冲液
95％甲酰胺
10mM EDTA，pH 8
装订＆amp;洗涤（B＆amp; W）缓冲液（2x）
10 mM Tris-HCl（pH 7.5）
1 mM EDTA
2 M NaCl

致谢

我们感谢M. Bizot和G. Palierne提供技术援助。这项工作由La Ligue Contre le Cancer，CancéropoleGrand Ouest，CNRS和雷恩大学1资助。作者声明没有竞争利益。

参考

Cox，M.P.，Peterson，D.A。和Biggs，P.J。（2010）。 SolexaQA：Illumina第二代测序数据的快速质量评估。 BMC Bioinformatics 11：485.
Langmead，B.，Trapnell，C.，Pop，M.和Salzberg，S.L。（2009）。短DNA序列与人类基因组的超快速记忆效率比对 基因组生物学 10（3）：R25。
Nicol，J.W.，Helt，G.A.，Blanchard，S.G.，Jr.，Raja，A。和Loraine，A.E。（2009）。综合基因组浏览器：用于分发和探索基因组数据集的免费软件。 Bioinformatics 25（20）：2730-2731。
Sérandour，A.A.，Avner，S.，Mahe，E.A.，Madigou，T.，Guibert，S.，Weber，M.and Salbert，G。（2016）。通过化学标记和核酸外切酶消化进行5-羟甲基胞嘧啶的单CpG解析图谱将进化上不受保护的CpG鉴定为TET目标。 Genome Biol 17:56。
Szulwach，K.E.，Song，C.X.，He，C.and Jin，P。（2012）。 5-hydroxymethylcytosine（5-hmC）specific enrichment。 Bio-protocol 2（15）。