甲基化扩增子测序
甲基化扩增子测序
1. 背景介绍
DNA甲基化在基因表达调控、疾病发生发展中扮演重要角色，特定基因或区域的甲基化异常与多种疾病密切相关。甲基化扩增子测序(Bisulfite Amplicon Sequencing，BSAS)技术通过设计特异性引物，对目标基因或区域进行定向扩增和深度测序，实现对候选区域甲基化状态的精准定量分析。
该技术结合了亚硫酸氢盐转化的单碱基分辨率优势和PCR扩增的靶向性，能够在单个CpG位点水平上检测甲基化差异。相比全基因组甲基化检测方法，扩增子测序具有成本低、通量高、检测深度深等特点，特别适合于大样本量的候选基因验证和临床标志物研究。
随着精准医学的发展，DNA甲基化生物标志物在肿瘤早筛、疾病诊断、预后评估和药物反应预测中展现出巨大应用价值。甲基化扩增子测序为这些临床应用提供了可靠、经济、高通量的检测解决方案，已成为表观遗传学转化研究的重要工具。
2. 甲基化扩增子测序检测原理
甲基化扩增子测序技术的核心原理是将亚硫酸氢盐转化技术与多重PCR扩增相结合，实现对特定基因组区域的靶向甲基化检测。首先对基因组DNA进行亚硫酸氢盐处理，未甲基化的胞嘧啶(C)被转化为尿嘧啶(U)，而5-甲基胞嘧啶保持不变。
检测流程包括：根据目标区域序列设计转化后特异性引物，引物设计需避开CpG位点以确保对甲基化和非甲基化模板的扩增效率一致；对转化后的DNA进行多重PCR扩增，可同时扩增8-96个目标区域；扩增产物经纯化后构建测序文库，进行高通量测序；通过生物信息学分析，将测序reads比对到参考序列，计算每个CpG位点的甲基化比例。
该方法的优势在于可实现超高测序深度(通常>5000X)，能够准确检测低频甲基化事件和样本间微小差异。单个扩增子长度通常为150-500bp，每个扩增子可包含数个至数十个CpG位点，提供区域内甲基化模式的详细信息。
3. 技术创新与优势
3.1 智能化引物设计系统
开发专业的引物设计平台，确保检测准确性：
(1) 避开CpG位点设计：引物结合区域避开CpG二核苷酸，确保对甲基化和非甲基化模板的等效扩增；
(2) 优化扩增效率：通过计算机模拟预测引物特异性和扩增效率，避免引物二聚体和非特异性扩增；
(3) 多重兼容设计：对于多重PCR panel，优化引物间相容性，确保所有靶标均衡扩增；
(4) 灵活panel配置：支持单个基因精细分析(tile设计)或多基因panel(8-96个靶标)检测。
3.2 超深度测序策略
采用高深度测序策略，实现精准定量：
(1) 每个扩增子测序深度>5,000X，部分关键位点可达10,000-50,000X；
(2) 高深度覆盖使甲基化定量误差<1%，能检测低至1-2%的甲基化差异；
(3) 适用于异质性样本(如肿瘤组织、循环肿瘤DNA)的低频甲基化检测；
(4) 支持单细胞水平甲基化异质性分析。
3.3 标准化数据分析流程
建立专业的生物信息学分析流程：
(1) 自动化质控：reads质量过滤、接头去除、低质量碱基修剪；
(2) 精准比对：针对亚硫酸氢盐转化序列优化的比对算法；
(3) 甲基化定量：单CpG位点甲基化水平计算，支持C/T转化率校正；
(4) 统计分析：组间差异检验、聚类分析、相关性分析等多种统计方法；
(5) 可视化输出：热图、棒图、散点图等多种可视化方式展示结果。
图1 甲基化扩增子检测流程和技术创新
4. 应用场景与服务优势
4.1 应用场景
甲基化扩增子测序技术在多个研究和临床领域具有广泛应用：
(1) 候选基因验证：验证全基因组甲基化筛查发现的候选位点，进行大样本量确证研究；
(2) 生物标志物开发：开发用于疾病诊断、预后评估的甲基化标志物panel；
(3) 临床检测应用：肿瘤早筛(如结直肠癌SEPT9甲基化检测)、液体活检(ctDNA甲基化分析)；
(4) 药物反应预测：检测与药物代谢、药物敏感性相关基因的甲基化状态；
(5) 队列研究：大规模人群队列中特定基因甲基化与疾病风险的关联研究；
(6) 治疗监测：评估去甲基化药物(如阿扎胞苷)的治疗效果；
(7) 印记基因研究：印记缺失综合征(如Prader-Willi综合征)的分子诊断。
(8) 基因治疗安全性评估:CAR-T细胞、基因编辑治疗产品中关键基因甲基化状态监测,确保治疗产品的表观遗传稳定性；
(9) 细胞治疗质控:干细胞治疗、iPSC衍生细胞产品的甲基化质量控制,评估细胞分化状态和安全性。
4.2 服务优势
(1) 高性价比：相比全基因组方法，成本降低90%以上，适合大样本量研究；
(2) 超高通量：单次可检测数百个样本，支持大规模临床队列研究；
(3) 灵活定制：支持单基因tile设计或多基因panel，根据研究需求灵活配置；
(4) 快速交付：标准服务周期3-4周，加急服务可缩短至2周；
(5) 质量保证：ISO质量管理体系认证，每批次包含阳性对照和阴性对照；
(6) 专业咨询：提供panel设计建议、数据解读和后续实验方案支持。
5. 甲基化扩增子测序分析报告示例
我们提供详细专业的甲基化扩增子测序分析报告，包含以下核心内容：
（1）样本质控信息：DNA浓度、纯度（A260/A280、A260/A230）、完整性检测结果（RIN评分）；
Sample_ID | DNA_Concentration_ng_uL | A260_A280_Ratio | A260_A230_Ratio | DNA_Integrity | RIN_Score |
|---|---|---|---|---|---|
Control_1 | 49.96 | 1.806 | 2.166 | Intact | 9.1 |
Control_2 | 96.06 | 1.887 | 2.042 | Intact | 8.9 |
Control_3 | 78.56 | 1.86 | 2.036 | Intact | 9.4 |
Exp_1 | 67.89 | 1.871 | 2.037 | Intact | 8.7 |
Exp_2 | 32.48 | 1.802 | 2.061 | Intact | 8.9 |
Exp_3 | 32.48 | 1.897 | 2.105 | Intact | 9 |
（2）测序数据统计：原始/有效reads数、Q30比例、测序深度、覆盖率；
Sample_ID | Raw_Reads | Clean_Reads | Q30_Percentage | GC_Content | Mean_Depth | Coverage_10X_Percent | Coverage_100X_Percent |
|---|---|---|---|---|---|---|---|
Control_1 | 8065725 | 9916182 | 94.72 | 48.06 | 5034 | 99.04 | 96.79 |
Control_2 | 11275709 | 10974675 | 93.8 | 48.92 | 7253 | 98.35 | 97.77 |
Control_3 | 8084654 | 10152991 | 92.05 | 48.96 | 6955 | 99.84 | 95.27 |
Exp_1 | 8953277 | 9870690 | 95.77 | 50.73 | 6585 | 99.47 | 95.59 |
Exp_2 | 9954354 | 9928388 | 94.25 | 50.44 | 6021 | 99.79 | 95.14 |
Exp_3 | 10688875 | 8028178 | 93.54 | 51.33 | 7613 | 99.7 | 95.98 |
（3）扩增子覆盖统计：每个扩增子的平均覆盖深度、中位覆盖深度、覆盖均一性；
Sample_ID | Amplicon | Mean_Coverage | Median_Coverage | Coverage_Uniformity |
|---|---|---|---|---|
Control_1 | Amplicon_1 | 4337 | 4678 | 0.958 |
Control_2 | Amplicon_1 | 4502 | 8293 | 0.914 |
Control_3 | Amplicon_1 | 7302 | 6037 | 0.864 |
Exp_1 | Amplicon_1 | 8737 | 4654 | 0.886 |
Exp_2 | Amplicon_1 | 7219 | 6711 | 0.863 |
Exp_3 | Amplicon_1 | 8682 | 5448 | 0.871 |
（4）亚硫酸氢盐转化效率：通过非CpG context的C-T转化率评估转化效率，通常>99%；
Sample_ID | Non_CpG_Sites_Analyzed | C_to_T_Conversion_Rate | CHH_Conversion | CHG_Conversion |
Control_1 | 62757 | 99.79 | 99.5 | 99.03 |
Control_2 | 69830 | 99.49 | 99.62 | 99.6 |
Control_3 | 67429 | 99.74 | 99.24 | 99.28 |
Exp_1 | 56893 | 99.46 | 99.4 | 99.17 |
Exp_2 | 68077 | 99.41 | 99.44 | 99.53 |
Exp_3 | 64373 | 99.59 | 99.56 | 99.45 |
（5）单CpG位点甲基化水平：每个CpG位点的甲基化百分比、覆盖深度、标准误差；
CpG_Site | Sample_ID | Methylation_Percentage | Coverage_Depth | Standard_Error |
|---|---|---|---|---|
Amplicon_1_CpG_1 | Control_1 | 59.66 | 1199 | 1.78 |
Amplicon_1_CpG_1 | Control_2 | 56.79 | 690 | 0.96 |
Amplicon_1_CpG_1 | Control_3 | 46.01 | 1846 | 0.75 |
Amplicon_1_CpG_1 | Exp_1 | 32.55 | 1065 | 1.54 |
Amplicon_1_CpG_1 | Exp_2 | 32.2 | 1811 | 1.51 |
Amplicon_1_CpG_1 | Exp_3 | 37.38 | 1461 | 0.88 |
Amplicon_1_CpG_2 | Control_1 | 38.15 | 1258 | 1.13 |
Amplicon_1_CpG_2 | Control_2 | 44.8 | 1420 | 1.42 |
Amplicon_1_CpG_2 | Control_3 | 29.99 | 1011 | 1.1 |
（6）区域平均甲基化：每个扩增子的平均甲基化水平及标准差；
Amplicon | Sample_ID | Mean_Methylation | Std_Dev | CpG_Count |
Amplicon_1 | Control_1 | 36.42 | 3.95 | 5 |
Amplicon_1 | Control_2 | 36.86 | 2.98 | 5 |
Amplicon_1 | Control_3 | 32.9 | 4.95 | 5 |
Amplicon_1 | Exp_1 | 51.24 | 2.45 | 5 |
Amplicon_1 | Exp_2 | 51.23 | 3.19 | 5 |
Amplicon_1 | Exp_3 | 41.84 | 4.92 | 5 |
（7）差异甲基化分析：对照组vs实验组的统计检验（t检验、Mann-Whitney U检验）、效应量（Cohen's D）、显著性判断；
Amplicon | Control_ Mean | Control _Std | Experimental_Mean | Experimental _Std | Delta _Methylation | T_Test_ Statistic | T_Test_ P_Value | Mann_ Whitney_U | Mann_ Whitney_P_Value | Cohens_D | Significance |
|---|---|---|---|---|---|---|---|---|---|---|---|
Amplicon_1 | 35.39 | 2.17 | 48.1 | 5.42 | 12.71 | -3.7681 | 0.0196 | 0 | 0.1 | 3.0766 | Yes |
Amplicon_2 | 32.01 | 0.62 | 61.3 | 1.08 | 29.29 | -40.7253 | 0 | 0 | 0.1 | 33.252 | Yes |
Amplicon_3 | 37.73 | 3.19 | 42.35 | 2.85 | 4.62 | -1.8712 | 0.1346 | 1 | 0.2 | 1.5279 | No |
（8）可视化结果：甲基化热图、棒图、散点图、箱线图等多维度展示；
图2 甲基化水平箱线图
图3 甲基化水平热图
（9）ROC曲线分析(如适用)：评估甲基化标志物的诊断效能，包括AUC、敏感性、特异性；
图4 甲基化扩增子ROC曲线
(10) 原始数据：提供每个样本的甲基化原始数据表格，支持客户自主深度挖掘。
6. 甲基化扩增子测序服务内容
服务流程 | 服务内容 |
|---|---|
项目咨询 | Panel设计咨询、引物设计策略讨论、样本量估算 |
引物设计与合成 | 专业引物设计、计算机模拟验证、引物合成与质检 |
样本质检 | DNA浓度、纯度、完整性检测 |
文库构建 | 亚硫酸氢盐转化、多重PCR扩增、文库纯化与质检 |
高通量测序 | DNBSEQ-T7平台PE150测序，每扩增子>5000X覆盖 |
数据分析 | 甲基化定量、差异分析、统计检验、数据可视化 |
报告交付 | 提供完整分析报告、原始数据、结果解读 |
增值服务 | Sanger测序验证、临床级检测方案开发、IVD试剂盒开发支持 |
7. 样本要求
类别 | 具体要求 |
|---|---|
DNA样本标准 | ·总量：≥500ng (标准)，最低可至50ng (少量样本模式)； ·浓度：≥20ng/μL； ·纯度：OD260/280=1.8~2.0， OD260/230≥1.8； ·完整性：基因组DNA主带清晰，无明显降解。 |
Panel信息 | ·提供目标基因/区域信息(基因名称、染色体坐标)； · 如有参考文献或已知CpG位点，请一并提供； ·支持8-96个靶标区域同时检测。 |
实验分组 | ·建议每组≥3个生物学重复；对于临床验证建议样本量≥30。 |
样本信息 | ·样本编号、类型(组织/血液/细胞)； ·分组信息(病例/对照，处理/未处理等)； ·物种信息及参考基因组版本。 |
特殊样本 | ·FFPE样本、cfDNA、古DNA等特殊样本请提前咨询，需定制优化方案。 |
*注：①对于临床样本或珍贵样本，建议先进行小规模预实验验证panel有效性后再进行大规模检测。②所有样本须符合上述质量标准，以确保检测结果的准确性和可靠性；③对于特殊样本类型，请提前与舒桐技术团队沟通（电话15336557985；产品订购/技术支持 service@generulor.com）。
8.参考文献
[1] Grunau, C., et al. (2001). Bisulfite genomic sequencing: systematic investigation of critical experimental parameters. Nucleic Acids Research, 29(13), E65.
[2] Tusnády, G. E., et al. (2005). BiSearch: primer-design and search tool for PCR on bisulfite-treated genomes. Nucleic Acids Research, 33(1), e9.
[3] Umer, M., & Herceg, Z. (2013). Deciphering the epigenetic code: an overview of DNA methylation analysis methods. Antioxidants & Redox Signaling, 18(15), 1972-1986.
[4] Korshunova, Y., et al. (2008). Massively parallel bisulphite pyrosequencing reveals the molecular complexity of breast cancer-associated cytosine-methylation patterns. Nature Methods, 5(7), 597-599.