亚硫酸氢盐测序

技术原理编辑本段

亚硫酸氢盐测序（Bisulfite Sequencing）的核心是亚硫酸氢盐对DNA中胞嘧啶的化学修饰。在碱性条件下，亚硫酸氢根离子（HSO3-）与胞嘧啶的5,6双键发生亲核加成，形成胞嘧啶-6-磺酸中间体；随后在酸性条件下脱氨基，生成尿嘧啶-6-磺酸；最后在碱性条件下脱磺酸基，尿嘧啶（U）释放。整个过程中，5-甲基胞嘧啶（5mC）由于C5位上的甲基阻碍了亲核加成，不发生上述反应。因此，经过亚硫酸氢盐处理的DNA中，未甲基化的C转化为U，而甲基化的C保持为C。

未修饰胞嘧啶经亚硫酸氢盐处理转化为U，扩增后读作T；5mC在处理中保留，扩增后读作C。

在后续PCR 扩增中，U与腺嘌呤（A）配对，被DNA 聚合酶掺入为胸腺嘧啶（T）；而5mC仍与鸟嘌呤（G）配对，保持为C。最终测序得到的序列中，原始未甲基化的C位点表现为T，原始甲基化的C位点表现为C。通过将测序读段（reads）与参考基因组比对，统计每个胞嘧啶位点上C和T的比例，即可计算出该位点的甲基化水平（C/(C+T)）。这种单碱基分辨率的定量能力是亚硫酸氢盐测序区别于甲基化敏感限制性内切酶、甲基化DNA 免疫沉淀（MeDIP）等方法的根本优势。

技术流程编辑本段

标准的亚硫酸氢盐测序流程包括以下关键步骤：

1. DNA提取与片段化： 提取高质量基因组DNA，避免RNA和蛋白质污染。对于全基因组亚硫酸氢盐测序（WGBS），通常通过超声打断或酶切将DNA片段化至200-500 bp。对于靶向测序，则通过杂交捕获或多重PCR富集目标区域。

2. 亚硫酸氢盐转化： 这是技术的核心步骤。常用试剂包括亚硫酸氢钠（NaHSO3）或亚硫酸氢铵，在pH 5.0-6.0、50-60°C条件下反应数小时。转化效率是关键质量指标，通常要求>99%。商业化试剂盒（如EZ DNA Methylation-Gold Kit）通过添加变性剂和抗氧化剂提高了转化效率和DNA回收率。近年来，酶法转化（如APOBEC3A介导的TET辅助吡啶硼烷测序）作为非亚硫酸氢盐替代方法出现，但亚硫酸氢盐转化仍是金标准。

3. 文库构建： 转化后的DNA为单链且富含U，需进行第二链合成或直接进行单链文库构建。传统方法采用随机引物合成第二链，但会导致序列组成偏差。目前主流采用亚硫酸氢盐测序专用文库构建策略，如PBAT（Post-Bisulfite Adaptor Tagging）或TAIL-PCR，以减少偏差并提高覆盖均一性。

4. 测序与数据分析： 根据通量需求选择Sanger测序（针对克隆片段）、焦磷酸测序（针对短区域）、或高通量测序（Illumina、MGI等平台）。数据分析流程包括：质量控制、去除接头和低质量碱基、将读段比对到经计算机模拟转化的参考基因组（如使用Bismark、BSMAP、BS-Seeker2等专用比对工具）、甲基化水平提取和差异甲基化分析。由于亚硫酸氢盐转化后DNA序列复杂度降低（C→T），比对时需采用三字母或四字母编码策略处理。

主要衍生技术编辑本段

全基因组亚硫酸氢盐测序（WGBS）： 覆盖全基因组所有胞嘧啶，是甲基化研究的金标准，但成本高、数据量大，且对于重复序列区域覆盖较差。

简化代表性亚硫酸氢盐测序（RRBS）： 利用MspI等限制性内切酶切割CCGG位点，富集CpG密集区域，以较低成本覆盖约1-3%的基因组，尤其适用于启动子和CpG岛分析。

靶向亚硫酸氢盐测序： 通过杂交捕获或多重PCR富集特定区域，适用于大样本队列的候选基因或全基因组关联研究（EWAS）验证。

亚硫酸氢盐克隆测序（BSP）： 将转化后的PCR产物克隆至载体，对单克隆进行Sanger测序，获得单分子甲基化模式，适用于短区域（<500 bp）的精细分析。

焦磷酸测序（Bisulfite Pyrosequencing）： 对短序列（<100 bp）进行定量甲基化检测，通量低但快速、定量准确，常用于临床样本的特定位点验证。

应用领域编辑本段

亚硫酸氢盐测序广泛应用于表观遗传学研究的各个层面。在基础研究中，它被用于绘制发育过程中DNA甲基化动态图谱，揭示印记基因、X 染色体失活和细胞重编程中的表观遗传调控机制。在癌症研究中，WGBS和RRBS被用于鉴定肿瘤抑制基因启动子的异常高甲基化及癌基因的低甲基化，为肿瘤分子分型和生物标志物开发提供依据。在植物科学中，该技术用于研究春化作用、杂种优势和胁迫响应中的DNA甲基化变化。此外，在微生物学中，亚硫酸氢盐测序用于研究细菌的甲基化组（methylome），揭示限制-修饰系统和DNA甲基化在毒力调控中的作用。

局限性与挑战编辑本段

亚硫酸氢盐测序的主要局限性包括：①DNA降解严重，亚硫酸氢盐处理会导致DNA大量断裂，限制了对微量样本和降解样本（如FFPE组织）的应用；②序列复杂度降低，C→T转化导致读段比对困难，尤其在重复序列区域；③无法区分5mC与5-羟甲基胞嘧啶（5hmC），因为两者在亚硫酸氢盐处理中均不发生转化；④对非CpG甲基化（如植物中的CHG和CHH）的检测需要更复杂的分析策略；⑤成本较高，WGBS对大规模样本的应用受限。针对这些问题，研究人员开发了氧化亚硫酸氢盐测序（oxBS-seq）以区分5mC和5hmC，以及酶法转化方法以克服DNA降解问题。

发展前景编辑本段

随着单细胞测序和长读长测序技术的发展，亚硫酸氢盐测序正朝着单细胞甲基化组和全基因组甲基化单倍型分析方向演进。单细胞亚硫酸氢盐测序（scBS-seq）和单细胞RRBS（scRRBS）已实现对单个细胞甲基化组的检测，为研究细胞异质性和发育谱系提供了新工具。同时，纳米孔测序平台通过直接识别碱基修饰（如5mC和5hmC），有望在保留长读长优势的同时避免亚硫酸氢盐转化带来的DNA损伤。然而，亚硫酸氢盐测序作为甲基化检测的“金标准”，在可预见的未来仍将是验证新型甲基化检测方法准确性的参照标准。