RRBS甲基化测序
RRBS甲基化测序
1. 背景介绍
DNA甲基化是最重要的表观遗传修饰之一，在基因表达调控、胚胎发育、基因组印记和疾病发生等生物学过程中发挥关键作用。RRBS(Reduced Representation Bisulfite Sequencing，简化代表性亚硫酸氢盐测序)技术通过酶切富集基因组中CpG高密度区域，结合亚硫酸氢盐转化和高通量测序，以经济高效的方式获得全基因组范围内的甲基化信息。
RRBS技术于2005年首次报道，相比全基因组亚硫酸氢盐测序(WGBS)，RRBS只需较低的测序深度即可覆盖基因组中大部分CpG岛，特别适合于启动子区域和基因调控区域的甲基化分析。该技术已被广泛应用于发育生物学、肿瘤学、神经科学等领域的表观遗传学研究。
随着表观遗传学在精准医学中的重要性日益凸显，DNA甲基化标志物在疾病诊断、预后评估和治疗监测中展现出巨大应用价值。RRBS技术为大规模样本的甲基化筛查提供了理想的解决方案，已成为表观基因组学研究的重要工具。
2. RRBS检测原理
RRBS技术的核心原理是利用MspI限制性内切酶识别CCGG位点对基因组DNA进行片段化。由于CpG二核苷酸在基因组中分布不均匀，MspI酶切位点主要富集于CpG岛等CpG高密度区域，从而实现对这些功能重要区域的靶向富集。
检测流程包括：首先对高质量基因组DNA进行MspI酶切消化，产生富含CpG的DNA片段；经末端修复和接头连接后，进行片段大小筛选（通常为40-220bp和220-320bp两个区段）；随后对筛选的DNA片段进行亚硫酸氢盐转化处理，未甲基化的胞嘧啶被转化为尿嘧啶，而甲基化的胞嘧啶保持不变；最后通过PCR扩增和高通量测序获得甲基化信息。
通过生物信息学分析，将测序reads比对到参考基因组，根据C-T转化率计算每个CpG位点的甲基化水平。RRBS技术能够以单碱基分辨率检测甲基化状态，覆盖约10-12%的基因组区域，包含60-70%的CpG岛，主要集中在启动子及基因调控元件等功能重要区域。
图1 RRBS检测原理图
3. RRBS技术创新与优势
3.1 RRBS技术优势
RRBS甲基化技术和其他甲基化测序技术相比存在一定优势，特别是在成本、覆盖范围、测序深度、样本通量和周期等5个维度表现优异，是性价比最优的选择，特别适合：
（1）需要平衡成本和覆盖范围的项目；
（2）大规模样本(>50个)的队列研究；
（3）启动子和CpG岛区域的精细分析；
（4）肿瘤标志物的发现和初步验证。
图2 五种甲基化检测技术对比
3.2 优化的酶切策略
采用MspI酶切富集策略，具有以下优势：
(1) 精准靶向CpG岛：MspI识别CCGG序列，该序列在CpG岛中高度富集，实现对功能重要区域的有效捕获；
(2) 覆盖启动子区域：约85%的基因启动子区域包含CpG岛，RRBS能有效覆盖这些调控元件；
(3) 片段大小可控：通过双size selection策略，富集40-220bp和220-320bp片段，确保测序质量和覆盖度。
3.3 高效的文库构建体系
通过优化的文库构建流程,建立高效实验体系：
(1) 优化的亚硫酸氢盐转化条件，确保>99%的转化效率，同时最大限度保护DNA不被降解；
(2) 低循环数PCR扩增策略，减少PCR偏好性对甲基化定量的影响；
(3) 严格的质控体系，每个步骤进行质检，确保文库质量符合测序要求。
3.4 专业的生物信息学分析
整合完善的生物信息学分析流程，实现：
(1) 精确的甲基化水平计算：单碱基分辨率的甲基化定量分析；
(2) 差异甲基化区域(DMR)鉴定：识别样本间显著差异的甲基化区域；
(3) 功能注释和通路分析：整合基因功能和信号通路数据库，解读甲基化变化的生物学意义。
4. 应用场景与服务优势
4.1 应用场景
RRBS技术在多个研究领域具有广泛应用：
(1) 肿瘤表观遗传学研究：发现肿瘤特异性甲基化标志物，鉴定抑癌基因启动子异常甲基化；
(2) 发育生物学：追踪胚胎发育过程中的甲基化重编程，研究细胞分化的表观遗传调控；
(3) 疾病机制研究：探索复杂疾病（如糖尿病、心血管疾病、神经退行性疾病）的甲基化异常；
(4) 环境表观遗传学：评估环境因素（如饮食、污染物、应激）对DNA甲基化模式的影响；
(5) 大规模人群研究：进行多样本甲基化筛查,发现疾病相关的表观遗传变异。
4.2 服务优势
(1) 技术成熟：10年以上RRBS技术服务经验，已完成数千例样本分析；
(2) 质量保证：ISO9001质量管理体系认证，确保实验全流程标准化；
(3) 性价比高：相比WGBS、测序成本降低60-70%，适合大规模样本筛查；
(4) 分析专业：提供从基础分析到高级分析的全套生物信息学服务；
(5) 快速交付：标准服务周期4-6周，加急服务可缩短至3周。
5. RRBS分析报告示例
我们提供全面专业的RRBS分析报告，包含以下核心内容：
（1）数据质量评估：测序数据量统计、质量评分分布、GC含量分析、接头污染检测等；
Sample | Total_Reads | Clean_Reads | Clean_Rate_% | Q20_% | Q30_% | GC_Content_% |
Control_1 | 20219110 | 19691316 | 97.39 | 96.55 | 95.12 | 50.39 |
Control_2 | 19570006 | 18625606 | 95.17 | 98.6 | 94.4 | 50.83 |
Control_3 | 18787201 | 18254771 | 97.17 | 98.82 | 92 | 51.97 |
Treatment_1 | 19396025 | 18731569 | 96.57 | 97.3 | 93.16 | 50.45 |
Treatment_2 | 18654811 | 18267027 | 97.92 | 96.7 | 92.36 | 50.47 |
Treatment_3 | 21474675 | 20782598 | 96.78 | 96.14 | 94.43 | 48.68 |
（2）比对统计：唯一比对reads比例、覆盖深度分布、转化率统计等关键指标；
Sample | Clean_Reads | Mapped_Reads | Mapping_Rate_% | Bisulfite_Conversion_Rate_% | Average_Coverage_Depth |
Control_1 | 19691316 | 14667835 | 74.49 | 99.38 | 26 |
Control_2 | 18625606 | 13381069 | 71.84 | 99.5 | 25.3 |
Control_3 | 18254771 | 13525544 | 74.09 | 99.39 | 30.2 |
Treatment_1 | 18731569 | 13199593 | 70.47 | 99.67 | 34.4 |
Treatment_2 | 18267027 | 13508151 | 73.95 | 99.25 | 27 |
Treatment_3 | 20782598 | 13602682 | 65.45 | 99.36 | 33.3 |
（3）CpG覆盖分析：总体CpG覆盖数量、CpG岛覆盖统计、启动子区域覆盖情况；
Sample | Total_CpG_Covered | CpG_Islands_Covered | CpG_Islands_Coverage% | Promoter_CpG_Covered | Promoter_Coverage% |
Control_1 | 3774329 | 20264 | 66.97 | 188571 | 69.93 |
Control_2 | 3955808 | 20027 | 67.05 | 185258 | 69.39 |
Control_3 | 4128776 | 21099 | 71.04 | 210535 | 71.15 |
Treatment_1 | 3956551 | 20773 | 71.21 | 185056 | 70.68 |
Treatment_2 | 4115270 | 20842 | 68.66 | 201518 | 68.78 |
Treatment_3 | 3610687 | 19923 | 71.35 | 186374 | 70.87 |
（4）全基因组甲基化水平：染色体甲基化分布、CpG/CHG/CHH context分析、基因组元件甲基化模式；
① 染色体甲基化分布：
Sample | Chromosome | Methylation_% | CpG_Count |
Control_1 | chr1 | 73.02 | 65708 |
Control_2 | chr1 | 71.75 | 275732 |
Control_3 | chr1 | 76 | 239407 |
Treatment_1 | chr1 | 68.78 | 201836 |
Treatment_2 | chr1 | 70.08 | 71677 |
Treatment_3 | chr1 | 66.84 | 171790 |
② CpG/CHG/CHH context与基因组元件甲基化模式分析：
Sample | Overall_Methylation_% | CpG_Context_% | CHG_Context_% | CHH_Context_% | CpG_Island_% | CpG_Shore_% |
Control_1 | 72.48 | 76.79 | 1.31 | 0.62 | 33.71 | 56.04 |
Control_2 | 74.42 | 76.69 | 0.82 | 0.36 | 27.28 | 52.27 |
Control_3 | 73.53 | 75.25 | 0.72 | 0.36 | 28.38 | 57.43 |
Treatment_1 | 69.09 | 72.92 | 1.46 | 0.43 | 29.97 | 51.01 |
Treatment_2 | 69.51 | 70.15 | 0.78 | 0.75 | 27.4 | 49.45 |
Treatment_3 | 70.02 | 72.28 | 0.74 | 0.66 | 28.68 | 54.32 |
（5）差异甲基化DMR分析：差异甲基化DMR统计、火山图和热图可视化；
DMR_ID | Chromosome | Start | End | Length_bp | CpG_Count | Log2_Fold_Change | P_Value | Regulation |
DMR_0001 | chr15 | 93525563 | 93526304 | 741 | 19 | -0.937 | 2.77E-13 | Hypomethylated |
DMR_0002 | chr10 | 82251088 | 82252819 | 1731 | 13 | 0.331 | 1.23E-09 | Hypermethylated |
DMR_0003 | chr11 | 92893487 | 92893797 | 310 | 6 | -0.481 | 1.57E-08 | Hypomethylated |
DMR_0004 | chr8 | 39660995 | 39662914 | 1919 | 7 | -1.05 | 9.76E-12 | Hypomethylated |
DMR_0005 | chr18 | 17618153 | 17618886 | 733 | 19 | -0.484 | 5.35E-08 | Hypomethylated |
图3 差异甲基化DMR火山图
图4 差异甲基化DMR热图
（6）差异甲基化CPG分析：差异甲基化CpG位点统计、火山图；
CpG_ID | Chromosome | Position | Log2_Fold_Change | P_Value | Q_Value | Regulation |
CpG_00001 | chr20 | 70198496 | 0.647 | 2.25E-04 | 7.84E-04 | Hypermethylated |
CpG_00002 | chr22 | 91693120 | -1.275 | 1.55E-13 | 3.70E-13 | Hypomethylated |
CpG_00003 | chr19 | 62065507 | -0.532 | 8.37E-14 | 1.33E-13 | Hypomethylated |
CpG_00004 | chr16 | 9288515 | 0.345 | 1.45E-16 | 4.87E-16 | Hypermethylated |
CpG_00005 | chr4 | 57998429 | 0.707 | 2.70E-19 | 6.59E-19 | Hypermethylated |
图5 差异甲基化CpG火山图
（7）候选基因分析：重点基因甲基化模式展示、基因结构注释。
图6 重点基因甲基化水平展示
6. RRBS检测服务内容
服务流程 | 服务内容 |
项目咨询与评估 | 根据研究目的制定个性化检测方案,进行项目报价 |
样本接收与质检 | 严格按标准对DNA样本进行质检,确保符合建库要求 |
RRBS文库构建 | MspI酶切、片段筛选、亚硫酸氢盐转化、文库扩增 |
高通量测序 | 文库质检合格后进行PE150测序,人20M reads,小鼠10M reads |
生物信息学分析 | 甲基化水平计算、DMR鉴定、功能注释与通路分析 |
专业报告交付 | 提供标准化分析报告,含技术解读与咨询服务 |
增值服务 | 个性化分析、候选位点验证、多组学联合分析 |
7. 样本要求
类别 | 具体要求 |
DNA样本标准 | ·总量:≥1μg (Qubit定量); ·浓度:≥50ng/μL; ·纯度:OD260/280=1.8~2.0; ·完整性:主带清晰，无明显降解(需提供琼脂糖凝胶电泳图)。 |
实验分组要求 | ·建议每组≥3个生物学重复，以确保统计学显著性。 |
样本信息 | ·样本类型及名称; ·物种信息(人、小鼠、大鼠或其他); ·实验设计和分组信息。 |
其他样本类型 | ·可提供组织(>50mg)或细胞(>2×10⁷)样本，我们提供DNA提取服务。 |
*注:①所有样本须符合上述质量标准，以确保检测结果的准确性和可靠性；②对于特殊样本类型，请提前与舒桐技术团队沟通（电话 15336557985；产品订购/技术支持 service@generulor.com）。
8. 参考文献
[1] Meissner, A., et al. (2005). Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Research, 33(18), 5868-5877.
[2] Gu, H., et al. (2011). Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nature Protocols, 6(4), 468-481.
[3] Smith, Z. D., & Meissner, A. (2013). DNA methylation: roles in mammalian development. Nature Reviews Genetics, 14(3), 204-220.
[4] Boyle, P., et al. (2012). Gel-free multiplexed reduced representation bisulfite sequencing for large-scale DNA methylation profiling. Genome Biology, 13(10), R92.