RRBS甲基化程度测定

技术概述

RRBS甲基化程度测定是一种基于亚硫酸氢盐转化与限制性内切酶消化的高通量测序技术，全称为Reduced Representation Bisulfate Sequencing，即简化代表性亚硫酸氢盐测序。该技术结合了限制性内切酶富集和高通量测序的优势，能够在全基因组范围内对DNA甲基化状态进行精准、的检测分析。

DNA甲基化作为表观遗传学研究中最为重要的修饰方式之一，在基因表达调控、基因组印记、X染色体失活以及肿瘤发生发展等生物学过程中发挥着关键作用。RRBS技术通过利用甲基化敏感的限制性内切酶（如MspI）对基因组DNA进行酶切，富集CpG密集区域，随后进行亚硫酸氢盐转化，将未甲基化的胞嘧啶转化为尿嘧啶，而甲基化的胞嘧啶则保持不变。通过高通量测序和生物信息学分析，可以实现对目标区域甲基化状态的准确定量分析。

与全基因组亚硫酸氢盐测序（WGBS）相比，RRBS技术具有显著的技术优势。首先，该技术通过酶切富集策略大幅降低了测序成本，测序量仅需WGBS的1/10至1/20即可获得高质量的甲基化数据。其次，RRBS技术专注于CpG岛等甲基化热点区域，数据覆盖度和深度更高，能够更准确地检测低水平甲基化变化。此外，该技术成熟稳定，实验流程标准化程度高，数据质量可靠，已成为表观遗传学研究领域广泛应用的主流技术方案。

RRBS甲基化程度测定技术的检测原理基于以下核心机制：MspI酶（CCGG识别位点）能够识别并切割基因组中的CCGG序列，而无论该序列中的胞嘧啶是否发生甲基化修饰。通过酶切可以获得大量富含CpG位点的DNA片段，这些片段经过接头连接、亚硫酸氢盐转化、PCR扩增后进行测序。在测序数据分析中，通过比对转化后的序列与参考基因组，可以准确识别每个CpG位点的甲基化状态和甲基化比例。

该技术的主要技术特点包括：覆盖全基因组约12-15%的CpG位点，主要集中在CpG岛、启动子区域和基因编码区；单碱基分辨率精准检测；可检测差异甲基化区域（DMR）；适合多样本比较分析；数据重复性好，技术误差低。这些特点使RRBS成为表观遗传学研究和生物标志物开发的理想技术选择。

检测样品

RRBS甲基化程度测定技术对样品类型具有较广泛的适用性，可接收多种来源的生物样品进行检测分析。不同类型的样品在采集、保存和运输过程中需要遵循相应的规范要求，以确保DNA质量和甲基化信息的完整性。

• 动物组织样品：包括各种实验动物（如小鼠、大鼠、斑马鱼等）的各器官组织，如肝脏、肾脏、心脏、脑组织、脾脏、肺脏等。组织样品应在手术分离后立即冷冻保存于-80°C环境中，避免反复冻融。样品量要求一般不低于20mg，建议提供50-100mg以保证充分的DNA提取。
• 植物组织样品：包括各种植物的叶片、根、茎、种子、花器官等组织。植物样品采集后应液氮速冻，-80°C保存。由于植物细胞壁和多糖多酚类物质的存在，植物样品的DNA提取需要特殊处理流程。
• 血液样品：包括全血、血清、血浆等血液成分。全血样品推荐使用EDTA抗凝管采集，避免使用肝素抗凝以防干扰后续实验。血液样品可在4°C短期保存或-80°C长期保存，建议采集后尽快进行DNA提取以保证质量。
• 细胞样品：包括各种培养细胞系、原代细胞、临床样本分离细胞等。细胞样品应收集后离心洗涤去除培养基成分，细胞沉淀液氮速冻后-80°C保存。细胞数量要求一般不低于10^6个细胞。
• 微生物样品：包括细菌、真菌等微生物样品，可用于研究微生物表观遗传调控机制。样品需要保证纯培养状态，避免杂菌污染影响检测结果。
• 福尔马林固定石蜡包埋（FFPE）样品：RRBS技术同样适用于临床病理存档的FFPE样品检测，为临床回顾性研究提供了技术支持。但FFPE样品DNA质量通常有所下降，需要特殊的质量控制和数据处理流程。
• 体液样品：包括尿液、唾液、脑脊液、胸腹水等。此类样品中细胞含量较低，需要特殊处理以获得足够量的DNA用于检测。

样品DNA质量要求方面，RRBS甲基化程度测定对DNA样品有严格的技术标准。DNA总量应不低于500ng，浓度不低于20ng/μL；DNA纯度方面，OD260/OD280比值应在1.8-2.0范围内；DNA完整性方面，需要主条带清晰，无明显降解。样品应使用无DNA酶的离心管保存，干冰运输，避免反复冻融。

检测项目

RRBS甲基化程度测定可提供多层次、多维度的检测分析内容，满足不同研究目标的数据需求。根据研究目的和数据分析深度的不同，检测项目可分为基础分析项目和高级分析项目两大类。

基础检测分析项目包括：

• CpG位点甲基化水平定量：对每个检测到的CpG位点进行甲基化比例计算，以甲基化胞嘧啶占总胞嘧啶的百分比表示。可检测单碱基分辨率的甲基化状态，准确度达到1%以内。
• CpG岛甲基化分析：对位于CpG岛区域内的CpG位点进行汇总分析，计算整体甲基化水平。CpG岛通常定义为GC含量大于50%、CpG观测值与期望值比值大于0.6、长度大于200bp的区域。
• 启动子区域甲基化分析：分析基因启动子区域（通常为转录起始位点上下游各1000bp范围）的甲基化状态，研究甲基化对基因表达调控的影响。
• 基因本体区域甲基化分析：包括基因外显子、内含子区域的甲基化状态分析，探索基因内部甲基化分布规律。
• 甲基化密度图谱绘制：可视化展示全基因组甲基化分布特征，包括染色体水平甲基化密度分布和基因组特征区域甲基化分布热图。
• 样品间相关性分析：评估样品间甲基化谱的一致性和差异性，进行样品聚类分析。

高级检测分析项目包括：

• 差异甲基化区域（DMR）鉴定：通过统计学方法比较不同样品组之间的甲基化差异，鉴定显著差异甲基化的基因组区域。常用的统计方法包括Fisher准确检验、Beta-binomial模型等。
• 差异甲基化位点（DMS）分析：鉴定单个CpG位点水平的甲基化差异，提供更高分辨率的甲基化变化信息。
• DMR相关基因功能注释：对差异甲基化区域相关基因进行GO功能注释和KEGG通路富集分析，揭示甲基化变化可能影响的生物学功能和信号通路。
• 甲基化与表达关联分析：结合转录组数据，分析启动子甲基化与基因表达水平的相关性，构建甲基化调控网络。
• 转录因子结合位点甲基化分析：分析位于转录因子结合位点区域的甲基化状态，研究甲基化对转录因子结合的调控作用。
• 印记基因甲基化分析：鉴定和分析印记基因的甲基化状态，研究基因组印记的表观遗传调控机制。
• 重复序列甲基化分析：分析基因组中重复元件（如Alu、LINE、SINE等）的甲基化状态，研究重复序列甲基化与基因组稳定性之间的关系。

检测报告通常包含以下核心数据内容：原始测序数据、比对结果统计、CpG位点覆盖度统计、全基因组甲基化水平统计、各基因组特征区域甲基化统计、差异甲基化分析结果、功能富集分析结果等。

检测方法

RRBS甲基化程度测定的实验流程包括样品准备、文库构建、测序分析和数据解读四个主要阶段。每个阶段都有严格的操作规范和质量控制标准，以确保检测结果的准确性和可靠性。

第一阶段：DNA样品准备

首先对送检样品进行DNA提取和质量检测。DNA提取采用酚-氯仿法或商业化DNA提取试剂盒，提取过程中注意避免DNA降解和污染。提取后的DNA使用NanoDrop分光光度计检测浓度和纯度，琼脂糖凝胶电泳检测DNA完整性。只有符合质量标准的DNA样品方可进入后续实验流程。

第二阶段：RRBS文库构建

RRBS文库构建是整个检测流程的核心环节，包括以下关键步骤：

• 限制性内切酶消化：使用MspI限制性内切酶对基因组DNA进行酶切，MspI识别并切割CCGG序列，产生富含CpG位点的DNA片段。酶切反应在37°C条件下进行，反应时间通常为2-4小时。
• 末端修复与加A尾：对酶切片段进行末端修复，使其成为平末端，随后在3'末端添加单个腺嘌呤核苷酸，为接头连接做准备。
• 甲基化接头连接：连接预先合成的甲基化测序接头，接头序列经过特殊设计，含有测序引物结合位点和样本索引序列。甲基化接头可以抵抗后续亚硫酸氢盐处理，保证接头序列的完整性。
• 片段大小筛选：使用琼脂糖凝胶电泳或固相可逆固定化（SPRI）磁珠筛选特定长度范围的DNA片段，通常选择40-220bp范围内的片段，该范围能够富集大部分CpG岛区域。
• 亚硫酸氢盐转化：使用亚硫酸氢钠对筛选后的DNA片段进行处理，将未甲基化的胞嘧啶转化为尿嘧啶，而甲基化的胞嘧啶保持不变。转化反应通常在变性条件下进行，温度50-65°C，反应时间10-20分钟，转化效率需达到99%以上。
• PCR扩增：使用特异性引物对转化后的DNA进行PCR扩增，富集文库分子。扩增过程中使用高保真聚合酶，并设置阴性对照以监测污染情况。
• 文库质检：使用Agilent Bioanalyzer或TapeStation对文库进行片段大小分布分析，使用qPCR进行文库浓度定量，确保文库质量符合测序要求。

第三阶段：高通量测序

质检合格的文库使用Illumina测序平台进行高通量测序。测序模式采用双端测序（Paired-end），读取长度通常为150bp。根据研究目的和样本数量，可以选择不同的测序数据量，一般建议每个样品测序数据量不低于5G以保证足够的覆盖深度。测序过程中实时监测测序质量指标，包括Q30比例、簇密度、相位误差等参数。

第四阶段：生物信息学分析

测序数据的生物信息学分析流程包括以下步骤：

• 原始数据质控：使用FastQC软件评估原始测序数据的质量，包括碱基质量分布、GC含量分布、接头污染情况等。
• 数据预处理：使用Trimmomatic或Cutadapt软件去除低质量碱基和接头序列，使用Trim Galore软件识别并去除随机引物序列。
• 序列比对：使用Bismark或BSMAP软件将处理后的序列比对到参考基因组，生成比对结果文件。
• 甲基化calling：提取比对序列中的甲基化信息，计算每个CpG位点的甲基化比例和覆盖深度。
• 下游统计分析：使用R语言和Python进行差异甲基化分析、功能富集分析、可视化图表绘制等。

整个实验流程需要严格的质量控制，包括DNA质量检测、酶切效率验证、转化效率评估、文库质检和测序质量监控等多个质控节点，确保检测数据的准确可靠。

检测仪器

RRBS甲基化程度测定涉及多个实验环节，需要使用多种仪器设备完成样品处理、文库构建、质量检测和测序分析等工作。主要检测仪器设备包括以下几类：

样品制备类仪器：

• NanoDrop超微量分光光度计：用于DNA样品浓度和纯度的快速检测，测量范围覆盖紫外可见光区，检测下限可达2ng/μL。
• Qubit荧光光度计：采用荧光染料结合的方法进行DNA准确定量，检测灵敏度高，特别适合低浓度DNA样品的准确测量。
• Agilent 2100 Bioanalyzer生物分析仪：基于微流控芯片技术的核酸定性定量分析系统，可快速检测DNA片段大小分布和浓度，是文库质检的核心设备。
• Agilent TapeStation系统：自动化核酸电泳分析平台，支持高通量文库质检，可同时处理多个样品。
• 普通PCR仪和梯度PCR仪：用于文库扩增PCR反应，梯度PCR仪可用于条件优化。

文库构建类仪器：

• MspI限制性内切酶反应系统：准确控制酶切反应温度和时间，确保酶切效率和特异性。
• EZ DNA Methylation-Gold Kit亚硫酸氢盐转化系统：商业化转化试剂盒配套设备，转化效率可达99%以上。
• AMPure XP磁珠纯化系统：用于文库构建过程中DNA片段的纯化和大小筛选，自动化程度高，重复性好。
• Thermal Cycler控温系统：准确控制PCR扩增过程中各阶段的温度，确保扩增效率和特异性。

高通量测序类仪器：

• Illumina NovaSeq系列测序系统：新一代高通量测序平台，单次运行可产生数Tb级测序数据，适合大规模样本的高通量测序。
• Illumina NextSeq系列测序系统：中等通量测序平台，适合中小规模研究项目的测序需求。
• Illumina MiSeq测序系统：小规模测序平台，适合验证性实验和小样本量研究。

生物信息分析类设备：

• 高性能计算服务器集群：配备多核CPU和大容量内存，用于海量测序数据的存储和计算分析。
• 生物信息学分析软件套件：包括FastQC、Trimmomatic、Bismark、BSMAP、SeqMonk等分析软件。

所有检测仪器设备均需要定期维护校准，建立完善的仪器使用记录和维护档案，确保仪器处于最佳工作状态，保障检测结果的稳定性和可靠性。

应用领域

RRBS甲基化程度测定技术作为表观遗传学研究的核心技术手段，在生命科学基础研究、临床医学、药物研发、农业科学等多个领域具有广泛的应用价值。

基础生命科学研究领域

在基础生命科学研究中，RRBS技术被广泛应用于探索DNA甲基化在生长发育、细胞分化、组织功能维持等生理过程中的作用机制。研究者利用该技术揭示了胚胎发育过程中全基因组甲基化重编程规律、干细胞分化过程中甲基化动态变化特征、以及各类生理状态下组织特异性甲基化图谱。这些研究为理解表观遗传调控的生物学意义提供了重要数据支撑。

肿瘤研究及临床诊断领域

DNA甲基化异常是肿瘤发生发展的标志性特征之一，RRBS技术在肿瘤研究中具有重要应用价值。该技术可用于鉴定肿瘤特异性甲基化标志物，开发新型肿瘤诊断方法。例如，通过对肿瘤组织和正常组织的比较分析，可以发现肿瘤特异性高甲基化或低甲基化区域，这些区域可作为肿瘤早期诊断的生物标志物。此外，甲基化检测还可用于肿瘤分型、预后评估、疗效监测和复发预测等方面。

• 肿瘤早期筛查：基于血液游离DNA的甲基化检测技术已成功应用于多种肿瘤的早期筛查，如肝癌、结直肠癌、肺癌等。
• 肿瘤分子分型：甲基化谱可用于肿瘤分子分型，指导个体化治疗方案制定。
• 预后评估：特定基因的甲基化状态与肿瘤患者预后密切相关，可作为预后评估指标。
• 疗效监测：动态监测肿瘤患者甲基化状态变化，可评估治疗效果和早期发现耐药。

神经科学研究领域

神经系统疾病如阿尔茨海默病、帕金森病、精神分裂症、抑郁症等与DNA甲基化异常密切相关。RRBS技术可用于研究神经系统疾病的表观遗传机制，探索疾病相关基因的甲基化调控规律，为疾病诊断和治疗靶点发现提供依据。

药物研发领域

在药物研发过程中，RRBS技术可用于药物作用机制研究、药物安全性评价和生物标志物开发等方面。表观遗传药物如DNA甲基转移酶抑制剂的临床研究需要甲基化检测技术支持。此外，该技术还可用于药物毒性评价，检测药物暴露后组织甲基化谱的变化，评估药物潜在的表观遗传毒性。

农业科学研究领域

在农业科学领域，RRBS技术被用于研究作物生长发育、逆境响应、产量品质形成等过程的表观遗传调控机制。通过比较不同品种、不同栽培条件下的甲基化差异，可以筛选与重要农艺性状相关的甲基化位点，为分子育种提供理论依据。在动物育种方面，该技术可用于研究优良性状的表观遗传基础，指导分子标记辅助育种。

环境健康研究领域

环境因素如化学污染物、重金属、辐射等可导致DNA甲基化改变，进而影响健康。RRBS技术可用于研究环境暴露与甲基化改变之间的关系，评估环境因素的健康风险，开发环境暴露的生物标志物。

法医学领域

DNA甲基化具有组织特异性和年龄相关性特点，在法医学领域可用于体液鉴定、年龄推断和表型刻画等应用。RRBS技术可为法医学甲基化标记的开发提供技术支持。

常见问题

问题一：RRBS与全基因组亚硫酸氢盐测序（WGBS）有何区别？

RRBS和WGBS都是基于亚硫酸氢盐转化的甲基化检测技术，主要区别在于检测范围和数据特点。WGBS对全基因组进行测序，理论上可以覆盖所有CpG位点，但测序成本高，部分区域覆盖深度不足。RRBS通过酶切富集策略，专注于CpG密集区域的检测，覆盖约12-15%的CpG位点，主要分布在CpG岛、启动子等功能重要区域。在测序成本方面，RRBS的测序量仅为WGBS的1/10至1/20，大幅降低了研究成本。在数据质量方面，RRBS的目标区域覆盖深度更高，数据重复性更好。研究选择时，如果研究目标是全基因组范围的甲基化图谱构建，建议选择WGBS；如果研究重点关注启动子、CpG岛等特定区域，或样本量较大、预算有限，RRBS是更经济的选择。

问题二：RRBS检测的样品DNA质量有什么要求？

RRBS技术对DNA样品质量有一定要求，主要包括以下几个方面：DNA总量要求不低于500ng，建议提供1μg以上以保证实验顺利进行；DNA浓度不低于20ng/μL；DNA纯度方面，OD260/OD280比值应在1.8-2.0范围内，OD260/OD230比值应大于2.0，表示样品中蛋白质和有机溶剂污染较低；DNA完整性方面，要求主条带清晰，无明显降解拖尾，对于基因组DNA，主条带应大于20kb。样品应避免使用含有EDTA的缓冲液长期保存，建议使用TE缓冲液或纯水保存。样品保存和运输过程中应避免反复冻融，推荐干冰运输。对于质量欠佳的样品，如FFPE样品或部分降解的DNA样品，可以通过优化实验方案进行检测，但数据质量可能有所降低。

问题三：RRBS测序数据量需要多少？

RRBS测序数据量需要根据研究目的和样本特点综合确定。一般情况下，每个样品推荐测序数据量为5-10G（双端150bp测序）。该数据量可确保大部分CpG位点获得足够的覆盖深度（通常大于10x），满足甲基化定量分析的统计学要求。如果研究目的是进行差异甲基化分析，建议适当增加测序数据量以提高统计效力。如果研究涉及复杂样本（如肿瘤组织异质性分析）或需要精细的甲基化定量，可以增加测序深度至15-20G。在实际研究中，可以通过预实验评估目标样品的甲基化特征和覆盖度，再确定正式实验的测序策略。

问题四：RRBS检测的重复性如何？

RRBS技术的技术重复性良好，相同样品的多次检测结果高度一致。在严格的质量控制条件下，技术重复样品之间的甲基化水平相关性通常大于0.95。影响重复性的主要因素包括：DNA质量、文库构建操作规范性、亚硫酸氢盐转化效率、测序质量等。通过标准化的实验流程和严格的质量控制措施，可以确保检测结果的技术重复性。需要注意的是，生物学重复之间的差异主要反映真实的生物学变异，与技术误差无关。建议研究设计时设置足够的生物学重复数量，以获得可靠的研究结论。

问题五：RRBS可以检测哪些类型的甲基化？

RRBS技术主要检测CpG二核苷酸位点的5-甲基胞嘧啶（5mC）。该技术基于亚硫酸氢盐转化原理，能够区分甲基化和非甲基化的胞嘧啶。需要注意的是，RRBS不能区分5-甲基胞嘧啶（5mC）和5-羟甲基胞嘧啶（5hmC），因为两者在亚硫酸氢盐处理条件下均不会被转化。如果需要特异性检测5hmC，需要采用氧化亚硫酸氢盐测序等技术方案。此外，RRBS对CpG位点进行检测，对于非CpG甲基化位点（如植物中的CHG、CHH位点）的检测覆盖度有限。

问题六：如何解读RRBS的甲基化数据？

RRBS甲基化数据解读需要结合研究目的和生物学背景进行综合分析。首先，需要关注数据质量指标，包括CpG位点覆盖度、测序深度分布、转化效率等参数。其次，需要从全基因组水平了解甲基化分布特征，如整体甲基化水平、各染色体甲基化分布、基因组特征区域（启动子、基因本体、CpG岛等）的甲基化模式。在差异分析中，需要关注差异甲基化区域的基因组位置、相关基因功能、通路富集等信息。甲基化水平通常以百分比形式表示，0%表示完全非甲基化，100%表示完全甲基化。在研究中，甲基化水平变化大于10-20%通常被认为具有生物学意义，但具体阈值需要根据研究背景和数据特征综合确定。建议结合基因表达数据或其他表观遗传数据进行多维度分析，深入理解甲基化的生物学功能。

问题七：RRBS检测结果是否可以验证？

RRBS检测结果可以通过多种方法进行验证。常用的验证方法包括：焦磷酸测序，可对特定CpG位点进行准确定量验证；甲基化特异性PCR（MSP），用于定性验证特定位点的甲基化状态；靶向亚硫酸氢盐测序，对目标区域进行深度测序验证。验证实验应选择代表性的差异甲基化位点或区域，建议选择3-5个位点进行验证。验证结果的一致性可以证明RRBS检测结果的可靠性，是研究质量控制的重要环节。

问题八：RRBS适用于植物样品检测吗？

RRBS技术同样适用于植物样品的甲基化检测，但需要注意植物样品的特殊性。植物基因组中CpG位点密度较低，MspI酶切位点数量相对较少，可能导致检测覆盖度降低。此外，植物基因组中含有大量重复序列和多倍体现象，增加了数据分析的复杂性。植物DNA提取过程中需要去除多糖、多酚等干扰物质。针对植物样品的RRBS检测，建议选择针对植物优化的实验方案，并使用适当的植物参考基因组进行数据分析。对于重点研究非CpG甲基化的植物研究，可能需要考虑其他技术方案。