技术概述

RRBS（Reduced Representation Bisulfite Sequencing）即简化代表性亚硫酸氢盐测序技术，是一种专门针对基因组中CpG岛区域进行高精度DNA甲基化分析的高通量测序方法。CpG岛是指基因组中富含CpG二核苷酸的区域，通常长度在500-2000bp之间，GC含量大于55%，这些区域在基因表达调控中发挥着至关重要的作用。

传统的全基因组甲基化测序（WGBS）虽然能够覆盖整个基因组，但成本高昂且数据量巨大，对于重点关注CpG岛区域的研究而言存在一定的资源浪费。RRBS技术通过限制性内切酶酶切的方式，特异性富集基因组中的CpG岛区域，在保证数据质量的前提下，大幅降低了测序成本和数据分析负担。

RRBS技术的基本原理是利用甲基化敏感的限制性内切酶（如MspI）对基因组DNA进行酶切，酶切位点主要位于CCGG序列处。由于CpG岛区域富含这种序列，因此酶切后可以特异性富集CpG岛片段。经过接头连接、亚硫酸氢盐转化、PCR扩增和测序等步骤，最终实现对CpG岛区域单碱基分辨率的甲基化检测。

亚硫酸氢盐处理是RRBS技术的核心环节，它能够将未甲基化的胞嘧啶（C）转化为尿嘧啶（U），而甲基化的胞嘧啶则保持不变。通过测序比对，可以准确判断每个CpG位点的甲基化状态，从而绘制出详细的甲基化图谱。

RRBS技术的优势在于其高覆盖率和单碱基分辨率。与甲基化芯片相比，RRBS不受探针设计的限制，能够发现新的甲基化位点；与MeDIP-seq相比，RRBS能够提供精确的甲基化水平定量信息。目前，RRBS已成为研究DNA甲基化修饰的主流技术手段，广泛应用于表观遗传学研究的各个领域。

CpG岛主要位于基因的启动子区域，约60%的人类基因启动子含有CpG岛。这些区域的甲基化状态与基因表达密切相关：当CpG岛发生甲基化时，通常会抑制基因转录；而去甲基化则往往伴随基因的激活表达。因此，通过RRBS分析CpG岛的甲基化模式，可以深入了解基因表达调控机制及其在疾病发生发展中的作用。

检测样品

RRBS CpG岛分析对样品的要求相对严格，高质量的DNA样品是获得可靠检测结果的基础。以下是常见的检测样品类型及其具体要求：

  • 动物组织样品：新鲜或冷冻保存的动物组织，要求组织量不少于50mg，建议使用液氮速冻后-80°C保存。组织样品需避免反复冻融，以防止DNA降解。常用的动物组织包括肝脏、脾脏、肾脏、脑组织、肿瘤组织等。
  • 植物组织样品：新鲜的叶片、根、茎等植物组织，建议取样后立即冷冻保存。植物样品由于细胞壁的存在，DNA提取难度较大，建议适当增加样品量。植物样品需特别注意多糖多酚对后续实验的影响。
  • 血液样品：采集的外周血需使用抗凝管保存，推荐使用EDTA抗凝管。血液样品应在采集后24小时内处理，或分离白细胞后冷冻保存。血液样品量建议不少于2mL全血。
  • 培养细胞样品：体外培养的细胞系，需收集处于对数生长期的细胞。细胞数量建议不少于10^6个，离心收集后用PBS洗涤，去除培养基残留。细胞沉淀可于-80°C保存或直接提取DNA。
  • 基因组DNA样品：已提取的基因组DNA样品，要求浓度不低于50ng/μL，总量不少于1μg。DNA样品需具备良好的完整性，OD260/280比值在1.8-2.0之间，无明显降解。建议使用琼脂糖凝胶电泳检测DNA完整性。
  • 临床病理样品：石蜡包埋组织（FFPE）样品也可用于RRBS分析，但需注意FFPE样品中DNA可能存在降解和化学修饰。建议对FFPE样品进行质量评估后再进行检测。新鲜冷冻病理组织样品效果更佳。

样品运输过程中需保持低温条件，推荐使用干冰或冰袋进行冷链运输。对于特殊样品类型，如微量样品、古DNA样品等，需根据具体情况进行评估和优化实验方案。

检测项目

RRBS CpG岛分析涵盖多个层面的检测项目，从基础的甲基化位点检测到综合的生物信息学分析，为研究者提供全面的甲基化研究数据：

  • CpG位点甲基化水平检测：对每个CpG位点的甲基化比例进行定量分析，甲基化水平以0-1之间的数值表示，0代表完全未甲基化，1代表完全甲基化。RRBS可实现单碱基分辨率的甲基化检测。
  • CpG岛区域甲基化分析：统计每个CpG岛区域的平均甲基化水平、甲基化密度分布等指标。识别高甲基化和低甲基化的CpG岛区域，分析其与基因启动子的对应关系。
  • 启动子甲基化分析：针对基因启动子区域的CpG岛进行专项分析，建立启动子甲基化与基因表达的相关性。识别差异甲基化启动子区域，预测其对基因表达的影响。
  • 差异甲基化区域（DMR）鉴定：比较不同样品组之间的甲基化差异，筛选具有统计学意义的差异甲基化区域。DMR分析是表观遗传学研究的核心内容，可为疾病标志物的发现提供线索。
  • 基因本体（GO）功能注释：对差异甲基化基因进行GO功能分类，分析甲基化变化涉及的生物学过程、分子功能和细胞组分。
  • 信号通路富集分析：利用KEGG、Reactome等数据库对差异甲基化基因进行通路富集分析，揭示甲基化调控相关的生物学通路。
  • 甲基化图谱绘制：可视化展示基因组甲基化分布情况，包括染色体水平甲基化分布图、CpG岛甲基化热图、启动子甲基化图谱等。
  • 甲基化与转录组关联分析：结合RNA-seq数据，分析甲基化状态与基因表达的关系，构建甲基化调控网络。

检测方法

RRBS CpG岛分析的实验流程包括样品准备、文库构建、上机测序和数据分析四个主要阶段：

一、基因组DNA提取与质量控制

使用适当的DNA提取试剂盒从样品中提取基因组DNA，通过NanoDrop或Qubit测定DNA浓度，琼脂糖凝胶电泳检测DNA完整性。DNA样品需满足浓度和纯度要求后方可进入后续实验流程。

二、限制性内切酶酶切

使用MspI限制性内切酶对基因组DNA进行酶切。MspI识别并切割CCGG序列，由于CpG岛富含该序列，酶切后可获得大量富含CpG位点的DNA片段。酶切反应需在适宜的温度和缓冲液体系中进行，确保酶切效率。

三、末端修复与接头连接

对酶切后的DNA片段进行末端修复，使其形成平末端。随后添加特异性接头，接头设计需考虑亚硫酸氢盐处理后的序列变化。接头连接效率直接影响后续文库构建的质量。

四、片段大小筛选

通过凝胶电泳或磁珠筛选的方法，回收特定大小范围的DNA片段（通常为40-220bp），这一步骤富集CpG岛区域片段，去除非特异性片段。

五、亚硫酸氢盐转化

这是RRBS实验的核心步骤。使用亚硫酸氢盐处理DNA，使未甲基化的胞嘧啶脱氨转化为尿嘧啶，而甲基化的胞嘧啶保持不变。转化效率直接影响甲基化检测的准确性，需选择转化效率高的试剂盒并严格按操作流程执行。

六、PCR扩增

对亚硫酸氢盐处理后的DNA片段进行PCR扩增，引入测序所需的标签序列。PCR条件需优化，既要保证扩增效率，又要避免过度扩增导致的偏好性。

七、文库质量检测与定量

使用Qubit对文库进行定量，通过qPCR检测文库的有效浓度。使用Agilent Bioanalyzer或TapeStation检测文库片段大小分布，评估文库质量。

八、高通量测序

文库质检合格后，在Illumina测序平台进行高通量测序。根据项目需求选择适当的测序深度，一般建议每个样品的测序数据量不低于5G。

九、生物信息学分析

测序获得的原始数据经过质量过滤、去除接头序列等预处理后，比对到参考基因组。随后进行甲基化位点识别、甲基化水平计算、差异甲基化区域分析等一系列生物信息学分析。

检测仪器

RRBS CpG岛分析实验过程涉及多种精密仪器设备，确保实验的准确性和可重复性：

  • Illumina NovaSeq 6000测序系统：目前最先进的高通量测序平台之一，具有高数据产出量、高测序质量和低错误率的特点。支持多种读长模式，可根据项目需求灵活配置。单次运行可产生数十Tb的测序数据，适合大规模甲基化研究项目。
  • Illumina NextSeq 550测序系统：中通量测序平台，适合中小规模的RRBS项目。具有操作简便、运行周期短的特点，适合对时效性要求较高的研究。
  • Agilent 2100 Bioanalyzer：用于DNA文库质量评估的微流控芯片电泳系统，可快速检测文库片段大小分布、浓度等指标，是文库质检的标准设备。
  • Agilent TapeStation系统：自动化电泳分析系统，用于DNA样品和文库的质量控制，具有高通量、自动化的特点。
  • Qubit 4荧光计：高灵敏度的DNA定量仪器，使用荧光染料特异性结合DNA，准确测定DNA浓度，不受RNA、蛋白质等其他物质的干扰。
  • StepOnePlus实时PCR系统：用于文库有效浓度定量，通过qPCR方法检测可用于测序的文库分子数，是文库定量金标准。
  • Bio-Rad T100 PCR仪：用于文库扩增和亚硫酸氢盐转化后DNA的PCR扩增，具有精确的温度控制能力。
  • Eppendorf离心机系列：包括高速冷冻离心机和微量离心机，用于样品离心、核酸沉淀、文库纯化等步骤。
  • NanoDrop分光光度计：用于DNA浓度和纯度的快速测定，通过紫外吸光度判断DNA质量。

所有仪器设备均需定期校准和维护，建立完善的仪器使用记录和质控体系，确保实验数据的可靠性和可追溯性。

应用领域

RRBS CpG岛分析在生命科学研究和医学领域具有广泛的应用，为揭示表观遗传学机制提供了强有力的技术支撑：

肿瘤表观遗传学研究

肿瘤发生发展过程中普遍存在DNA甲基化异常。抑癌基因启动子CpG岛的高甲基化导致基因沉默，而癌基因启动子的去甲基化则促进其表达。RRBS技术可系统筛选肿瘤相关甲基化标志物，为肿瘤早期诊断、预后评估和靶向治疗提供依据。在肺癌、肝癌、胃癌、结直肠癌、乳腺癌等多种肿瘤中，RRBS已成功鉴定出大量具有临床价值的甲基化标志物。

发育与分化研究

胚胎发育和细胞分化过程伴随大规模的甲基化重编程。通过RRBS分析不同发育阶段或分化状态的细胞，可绘制发育过程中的甲基化动态变化图谱，揭示干细胞多能性维持和定向分化的表观遗传调控机制。在诱导多能干细胞（iPSC）研究中，RRBS用于评估重编程过程中甲基化状态的恢复情况。

神经系统疾病研究

多种神经系统疾病与DNA甲基化异常相关，包括阿尔茨海默病、帕金森病、精神分裂症、抑郁症等。RRBS技术可分析脑组织或外周血中的甲基化改变，探索疾病相关的表观遗传学机制，寻找潜在的诊断标志物和治疗靶点。

心血管疾病研究

动脉粥样硬化、高血压、心肌病等心血管疾病的发病过程中存在甲基化调控异常。RRBS可分析血管组织、心肌组织或血液样品中的甲基化变化，研究环境因素（如饮食、压力）通过表观遗传学机制影响心血管健康的途径。

代谢性疾病研究

肥胖、2型糖尿病、非酒精性脂肪肝等代谢性疾病的发病与表观遗传学改变密切相关。RRBS可用于分析代谢器官（如肝脏、脂肪组织、胰岛）的甲基化变化，揭示营养和环境因素导致代谢紊乱的表观遗传学基础。

免疫性疾病研究

自身免疫性疾病（如系统性红斑狼疮、类风湿关节炎）患者的免疫细胞存在甲基化异常。RRBS可分析免疫细胞亚群的甲基化特征，探索免疫耐受打破的表观遗传学机制。

植物表观遗传学研究

在植物中，DNA甲基化参与基因表达调控、转座子沉默、基因组印记等生物学过程。RRBS可用于研究植物发育、逆境响应、杂种优势等过程中的甲基化变化，为作物遗传改良提供新思路。

药物表观遗传学效应评估

某些药物（如DNA甲基转移酶抑制剂）通过影响甲基化发挥治疗作用。RRBS可用于评估药物处理后的甲基化变化，阐明药物作用机制，为药物研发和临床应用提供参考。

常见问题

问题一：RRBS与全基因组甲基化测序（WGBS）有何区别？应如何选择？

RRBS和WGBS都是基于亚硫酸氢盐转化的甲基化检测技术，但存在明显差异。WGBS覆盖整个基因组，能够检测所有CpG位点，但测序成本高，数据量大。RRBS通过酶切富集CpG岛区域，虽然覆盖范围有所限制，但针对CpG岛的覆盖度可达80%以上，且成本仅为WGBS的1/3到1/2。如果研究重点是基因启动子区域的甲基化，或样品数量较多、预算有限，RRBS是更经济高效的选择；如果需要研究基因组非CpG岛区域（如基因体、增强子区域）的甲基化，则应选择WGBS。

问题二：RRBS检测对DNA样品有何特殊要求？

RRBS对DNA质量要求较高，主要包括以下几点：DNA总量建议不低于1μg，浓度不低于50ng/μL；DNA需保持良好的完整性，避免降解，可通过琼脂糖凝胶电泳检测；DNA纯度要高，OD260/280比值在1.8-2.0之间，不应含有蛋白质、RNA、有机溶剂等杂质。对于FFPE样品，由于DNA可能存在降解和化学修饰，需要评估DNA质量后再决定是否适合进行RRBS分析。此外，样品不应经过甲醛等化学固定处理，因为这些处理可能影响后续的亚硫酸氢盐转化效率。

问题三：RRBS的测序深度应如何设定？

RRBS的测序深度需根据研究目的和预算综合考虑。对于常规的甲基化谱图分析，每个样品产生5-10G的有效数据通常可以满足需求；如果进行精细的差异甲基化分析或稀有细胞群体研究，建议适当增加测序深度。一般来说，每个CpG位点获得10-30×的覆盖深度可获得较为准确的甲基化定量结果。测序深度的增加有助于提高甲基化检测的准确性和统计效力，但也会相应增加成本。

问题四：亚硫酸氢盐转化效率如何评估？

亚硫酸氢盐转化效率是影响RRBS数据质量的关键因素。通常通过以下方式评估转化效率：在实验中加入已知甲基化状态的对照DNA（如lambda DNA或pUC19质粒），这些对照DNA不含甲基化胞嘧啶，如果转化效率正常，这些DNA中的胞嘧啶应全部转化为尿嘧啶；也可分析内源性的非CpG胞嘧啶转化情况作为评估指标。一般要求亚硫酸氢盐转化效率达到99%以上。如果转化效率偏低，需要排查转化试剂、反应条件等因素。

问题五：RRBS数据分析的主要难点有哪些？

RRBS数据分析面临以下主要挑战：一是序列比对难度大，亚硫酸氢盐处理后DNA序列发生改变（C变为T），需要使用专门的比对算法（如Bismark、BSMAP）；二是重复序列区域的甲基化分析存在不确定性，因为短读长难以唯一比对到基因组；三是覆盖不均匀，由于酶切位点分布不均，某些区域的覆盖深度可能很高，而某些区域可能没有覆盖；四是批次效应校正，多个测序批次之间的数据需要校正批次差异。选择经验丰富的数据分析团队和成熟的分析流程对于获得可靠结果至关重要。

问题六：如何验证RRBS分析结果？

对于重要的甲基化发现，建议使用独立的方法进行验证。常用的验证方法包括：焦磷酸测序，可对特定区域的甲基化进行定量验证，是甲基化验证的金标准；甲基化特异性PCR（MSP），适合定性验证特定区域是否存在甲基化；靶向亚硫酸氢盐测序，使用引物扩增目标区域后测序验证。验证实验应使用独立的样品群体，以确认发现的可靠性。此外，可将RRBS结果与公共数据库中的甲基化数据进行交叉比对，评估结果的一致性。

问题七：RRBS能否用于甲基化定量分析？

RRBS是一种定量分析方法。对于每个CpG位点，RRBS通过统计比对到该位点的测序读段中甲基化读段和未甲基化读段的比例，计算甲基化水平。甲基化水平以β值表示，范围从0到1。例如，某CpG位点的甲基化水平为0.7，表示在该位点有70%的DNA分子处于甲基化状态。需要注意的是，RRBS的定量准确性受测序深度、比对质量、PCR偏好性等因素影响，在解读甲基化水平变化时应综合考虑这些因素。

问题八：RRBS研究需要多少生物学重复？

生物学重复的数量取决于研究目的、样品异质性和统计效力要求。一般来说，进行差异甲基化分析时，每组样品建议不少于3个生物学重复，以获得具有统计学意义的结果。对于临床样品研究，由于个体间异质性较大，建议每组不少于6-10个样品；对于细胞系或模式动物研究，由于遗传背景一致，每组3-5个样品通常可以满足需求。在设计研究方案时，可通过预实验评估样品间甲基化变异程度，据此确定合适的样品数量。

问题九：RRBS数据能否与其他甲基化数据进行整合分析？

RRBS数据可以与其他来源的甲基化数据进行整合分析。例如，可以将RRBS数据与甲基化芯片数据进行比较，评估不同平台的一致性；可以将RRBS数据与MeDIP-seq数据整合，获得更全面的基因组甲基化信息；可以将RRBS数据与RNA-seq数据整合，分析甲基化与基因表达的相关性。整合分析时需要注意不同技术的数据特点，进行适当的标准化和校正。使用生物信息学工具可以实现多源数据的可视化和关联分析。

问题十：RRBS检测周期一般多长？

RRBS CpG岛分析的检测周期通常为4-6周，具体包括：样品接收和质量检测1-2天，文库构建5-7天，上机测序3-5天，生物信息学分析7-14天，报告撰写和质量审核3-5天。如果样品数量较多或需要增加测序深度，周期可能相应延长。对于有特殊时效要求的项目，可通过优化实验流程、增加测序通量等方式缩短周期。建议在项目启动前与检测机构充分沟通，明确项目进度安排。