RRBS甲基化位点分析

技术概述

RRBS甲基化位点分析是一种基于亚硫酸氢盐转化与高通量测序相结合的DNA甲基化检测技术，全称为简化的亚硫酸氢盐甲基化测序技术。该技术通过限制性内切酶酶切基因组DNA，富集CpG密度较高的区域，结合亚硫酸氢盐处理将未甲基化的胞嘧啶转化为尿嘧啶，而甲基化的胞嘧啶保持不变，从而实现对全基因组范围内CpG位点甲基化状态的精准检测。相较于全基因组甲基化测序，RRBS技术以更低的测序数据量实现了对启动子区域、CpG岛等关键区域的高覆盖度检测，是目前表观遗传学研究中应用最为广泛的技术之一。

DNA甲基化作为最重要的表观遗传修饰之一，在基因表达调控、基因组印记、X染色体失活以及转座子沉默等生物学过程中发挥着关键作用。异常的DNA甲基化模式与多种疾病特别是肿瘤的发生发展密切相关。RRBS甲基化位点分析技术凭借其高性价比、高覆盖度和高准确性的特点，为科研人员提供了一种强有力的甲基化研究工具，广泛应用于生物医学研究、疾病标志物筛选、生长发育调控机制研究等领域。

RRBS技术的核心原理在于利用甲基化敏感的限制性内切酶（如MspI）对基因组DNA进行酶切，酶切位点包含CpG二核苷酸，通过选择特定大小的片段进行富集，可以实现CpG密度较高区域的有效覆盖。酶切后的DNA片段经过末端修复、加A尾、接头连接等建库步骤，然后进行亚硫酸氢盐处理，最后通过PCR扩增和高通量测序获得甲基化信息。数据分析时，通过比对亚硫酸氢盐处理后的序列与参考基因组，可以准确判断每个CpG位点的甲基化状态。

该技术的主要优势包括：首先，测序成本大幅降低，仅需全基因组亚硫酸氢盐测序十分之一左右的测序量即可覆盖大部分启动子和CpG岛区域；其次，单碱基分辨率的分析精度可以准确检测每一个CpG位点的甲基化水平；第三，通过优化实验流程和生物信息学分析方法，可以实现多样本的高通量并行检测，适合大规模样本的研究需求。

检测样品

RRBS甲基化位点分析技术对样品类型具有广泛的适用性，可以检测多种来源的DNA样品。为了保证检测结果的准确性和可靠性，对样品的质量和数量有一定的要求。

• 动物组织样品：包括新鲜组织、冷冻组织及石蜡包埋组织等，适用于肿瘤组织与癌旁组织的甲基化差异分析、不同发育阶段组织甲基化谱研究等。组织样品需保证DNA完整性好，无明显降解。
• 植物组织样品：包括叶片、根、茎、花、种子等各类植物组织，可用于植物生长发育、逆境胁迫响应、表观遗传变异等方面的研究。植物样品需特别注意去除多糖、多酚等次生代谢物的干扰。
• 血液样品：包括全血、血清、血浆、外周血单核细胞等，是临床研究和转化医学研究中最常用的样品类型之一，可用于肿瘤早期诊断标志物的筛选、疾病预后评估等研究。
• 细胞样品：包括培养细胞系、原代细胞、流式分选细胞等，适用于细胞分化、发育、凋亡等过程中的甲基化调控机制研究。
• 微生物样品：包括细菌、真菌等微生物基因组DNA，可用于微生物表观遗传学研究。
• 体液样品：包括尿液、唾液、脑脊液、胸腹水等，适合作为无创或微创检测的样品来源。

样品DNA质量要求方面，DNA浓度应大于20ng/μL，总量不少于500ng，DNA完整性要好，主带清晰，无明显降解。DNA纯度方面，OD260/OD280比值应在1.8-2.0之间，OD260/OD230比值应大于2.0，无明显蛋白污染、盐离子污染或其他抑制酶反应的物质残留。对于石蜡包埋组织等特殊样品，由于DNA可能存在一定程度降解，可根据实际情况调整建库方案。

检测项目

RRBS甲基化位点分析提供多层次、多维度的检测项目，涵盖从基础的单点甲基化水平检测到综合的甲基化谱分析，满足不同研究目的的需求。

基础检测项目包括：

• 全基因组CpG位点甲基化水平定量分析：对覆盖范围内所有CpG位点的甲基化比例进行定量检测，以数值形式呈现每个位点的甲基化百分比。
• CpG岛区域甲基化状态分析：重点关注CpG岛区域的甲基化分布特征，分析其与基因表达调控的关系。
• 启动子区域甲基化谱分析：系统分析基因启动子区域的甲基化状态，评估启动子甲基化对基因转录的影响。
• 基因本体区域甲基化分析：检测基因编码区、内含子区等区域的甲基化分布特征。

差异甲基化分析项目包括：

• 组间差异甲基化位点筛选：通过统计学方法筛选不同实验组之间甲基化水平存在显著差异的CpG位点。
• 差异甲基化区域识别：识别组间甲基化差异显著的连续区域，这些区域往往与功能基因的调控相关。
• 差异甲基化区域功能注释：对差异甲基化区域进行基因组位置注释、基因功能注释和通路富集分析。

高级分析项目包括：

• 甲基化与基因表达关联分析：结合转录组测序数据，分析甲基化状态与基因表达水平的相关性。
• 甲基化单倍型分析：分析连续CpG位点的甲基化组合模式，揭示更精细的甲基化调控特征。
• 时序性甲基化变化分析：针对时间序列样品，分析甲基化状态的动态变化规律。
• 肿瘤甲基化特征谱构建：针对肿瘤样品构建特征性甲基化谱，筛选潜在的肿瘤标志物。

数据分析报告内容涵盖测序数据质量控制、比对分析结果、甲基化水平统计、差异甲基化分析结果以及各类可视化图表，包括甲基化水平分布图、染色体分布图、聚类热图、散点图等。

检测方法

RRBS甲基化位点分析采用标准化的实验流程和严格的质控体系，确保检测结果的准确性和可重复性。整个检测流程包括样品检测、建库、测序和数据分析四个主要环节。

第一环节为样品质量检测。采用琼脂糖凝胶电泳检测DNA完整性和纯度，使用紫外分光光度计检测DNA浓度和纯度。对于合格样品，使用Qubit荧光计进行精确定量，确保进入建库反应的DNA量准确可控。

第二环节为文库构建。首先使用甲基化敏感的限制性内切酶MspI对基因组DNA进行酶切，该酶识别CCGG位点，可以富集CpG密度较高的区域。酶切产物经过片段大小筛选，选择40-220bp的片段进行后续建库。然后进行末端修复反应，将酶切产生的粘性末端补平为平末端。接下来在片段3'末端添加单个A碱基，便于后续接头连接。随后进行甲基化接头连接，接头序列预先设计含有甲基化胞嘧啶，可抵抗后续亚硫酸氢盐处理。连接产物经片段筛选后进行亚硫酸氢盐处理，将未甲基化的胞嘧啶转化为尿嘧啶。最后进行PCR扩增，引入测序引物结合位点和样本标签序列。

第三环节为高通量测序。构建完成的文库经过质量检测和定量后，使用Illumina平台进行双端测序，测序读长通常为150bp。测序深度根据研究目的和样本数量进行调整，一般每个样品的测序数据量为5-10G，可覆盖约200-400万个CpG位点。

第四环节为生物信息学分析。原始测序数据经过质量评估和过滤，去除低质量序列和接头序列。使用专门的甲基化比对软件将处理后的序列与参考基因组进行比对，同时保留胞嘧啶到胸腺嘧啶的转换信息。比对结果经过去重处理，去除PCR扩增产生的重复序列。然后进行甲基化 calling，计算每个CpG位点的甲基化比例。后续进行差异甲基化分析和功能注释分析。

整个实验流程设置多重质控节点，包括酶切效率质控、文库构建质控、测序数据质控和分析过程质控，确保每个环节均在可控范围内进行。

检测仪器

RRBS甲基化位点分析依托专业的高通量测序平台和配套的分子生物学实验设备，保证检测的高效性和准确性。

核心测序设备采用Illumina系列高通量测序平台，包括NovaSeq系列、HiSeq系列和MiSeq系列等。其中NovaSeq系列测序仪是目前最先进的高通量测序平台之一，具有超高的数据通量和灵活的运行模式，单次运行可产生数十Tb的测序数据，适合大规模样本的并行检测。该平台采用边合成边测序原理，配合成熟的测序化学和先进的流体学设计，可以提供高质量的测序数据，测序错误率低于0.1%。对于小规模项目，MiSeq系列测序仪可提供快速的周转时间和灵活的通量选择。

文库构建相关设备包括：

• 全自动核酸提取系统：用于从各类样品中自动化提取高质量基因组DNA，确保提取效率和纯度。
• 专业PCR扩增仪：用于文库构建过程中的PCR扩增反应，温度控制精准，升降温速率快。
• 片段筛选系统：用于酶切产物和文库片段的大小筛选，可精确控制目标片段范围。
• 全自动文库质检系统：用于文库浓度检测、片段大小分布分析等质量控制环节。

样品检测相关设备包括：

• 超微量分光光度计：用于DNA样品的浓度和纯度检测，可检测微量样品。
• Qubit荧光定量仪：采用荧光染料法对DNA进行精确定量，灵敏度高，不受RNA和蛋白质干扰。
• 凝胶成像系统：配合琼脂糖凝胶电泳检测DNA完整性和纯度。
• 生物分析仪：采用微流控芯片技术检测DNA片段大小分布，分辨率高。

生物信息学分析依托高性能计算集群，配备大容量存储设备和高性能计算节点，可高效处理海量测序数据，完成从原始数据到分析结果的全部流程。

应用领域

RRBS甲基化位点分析技术在生命科学研究和医学应用领域具有广泛的应用价值，为揭示表观遗传调控机制和筛选疾病相关标志物提供了有力工具。

在肿瘤学研究领域，RRBS技术被广泛应用于肿瘤发生发展机制的探索和肿瘤标志物的筛选。肿瘤细胞普遍存在全基因组低甲基化和特定区域高甲基化的特征性改变，这些甲基化异常与癌基因激活、抑癌基因沉默密切相关。通过对比肿瘤组织与正常组织的甲基化谱，可以筛选差异甲基化区域和关键调控位点，揭示肿瘤发生发展的表观遗传学机制。此外，循环肿瘤DNA的甲基化检测为肿瘤早期诊断、疗效评估和预后预测提供了新的无创检测手段。

在发育生物学研究领域，RRBS技术用于解析胚胎发育、细胞分化和器官形成过程中的甲基化动态变化规律。受精后基因组经历大规模甲基化重编程，建立合子型甲基化模式。研究不同发育阶段和不同组织器官的甲基化谱，有助于理解表观遗传调控在发育过程中的作用机制。体细胞重编程为诱导多能干细胞过程中也伴随着剧烈的甲基化改变，相关研究为再生医学提供理论基础。

在植物科学研究领域，RRBS技术应用于植物生长发育调控、逆境胁迫响应和表观遗传育种研究。植物甲基化模式与动物存在差异，植物基因组甲基化主要发生在基因本体区域而非启动子区域。研究植物在不同逆境条件下的甲基化变化，可以揭示植物适应环境的表观遗传机制，为抗逆品种培育提供新思路。

在表观遗传流行病学研究领域，RRBS技术用于分析环境因素与甲基化的关联，探索环境暴露导致的表观遗传改变及其健康效应。包括营养因素、环境污染、心理压力等均可能影响DNA甲基化状态，这些改变可能是环境因素致病的重要中介机制。

在转化医学和临床检测领域，RRBS技术用于疾病诊断标志物的筛选和验证。甲基化标志物具有稳定性好、检测方便等优势，已有多款基于甲基化检测的肿瘤早筛产品获批上市。此外，甲基化检测还可用于疾病分型、药物疗效预测等精准医学应用。

常见问题

问：RRBS与全基因组亚硫酸氢盐测序有什么区别？

答：两种技术均采用亚硫酸氢盐转化和高通量测序进行甲基化检测，主要区别在于覆盖范围和成本。全基因组亚硫酸氢盐测序理论上可以覆盖全基因组所有CpG位点，测序数据量大，成本较高；RRBS通过限制性内切酶酶切富集CpG密度较高的区域，以约十分之一的测序数据量覆盖大部分启动子和CpG岛区域，性价比更高。研究者可根据研究目的和预算选择合适的技术方案。

问：RRBS能检测多少个CpG位点？

答：RRBS覆盖的CpG位点数量与物种基因组特征、测序深度和样本质量相关。以人类样品为例，通常每个样品可覆盖200-400万个CpG位点，约占全基因组CpG位点的10-15%，但这些位点主要集中在启动子、CpG岛等功能重要区域，覆盖了约80%以上的CpG岛和启动子区域。

问：RRBS对样品DNA有什么要求？

答：RRBS建库需要高质量基因组DNA，建议DNA总量不少于500ng，浓度大于20ng/μL，完整性要好，无明显降解。DNA纯度要求OD260/OD280比值在1.8-2.0之间，无蛋白、盐离子或其他抑制物污染。对于DNA量较少的样品，可通过优化建库方案或增加PCR循环数进行处理，但可能对数据质量产生一定影响。

问：亚硫酸氢盐处理后DNA损失大怎么办？

答：亚硫酸氢盐处理确实会导致DNA损伤和损失，这是该技术的固有特性。实验中可通过优化转化条件、控制反应时间和温度来减少DNA损失。同时，建库过程中在亚硫酸氢盐处理前完成接头连接，可以有效保护目的片段。对于起始量较低的样品，可通过增加PCR循环数进行补偿，但需注意扩增偏好性对结果的影响。

问：如何判断RRBS数据的可靠性？

答：评价RRBS数据可靠性可从多个维度进行：测序数据质量方面，Q30比例应大于85%，GC含量分布正常；比对效率方面，比对率应大于70%；甲基化转化效率方面，可利用添加的未甲基化对照DNA评估转化效率，理想情况下应大于99%；覆盖度方面，CpG位点覆盖数量和覆盖深度应符合预期；重复性方面，技术重复样品的相关系数应大于0.95。

问：RRBS能否检测非CpG位点的甲基化？

答：RRBS主要针对CpG位点进行甲基化检测，因为在脊椎动物中，DNA甲基化主要发生在CpG二核苷酸位点。限制性内切酶MspI识别的CCGG序列中含有CpG位点，因此RRBS技术的设计目标即是富集CpG区域。对于植物中存在的CHG和CHH甲基化，需要采用全基因组亚硫酸氢盐测序或其他方法进行检测。

问：如何进行多组学联合分析？

答：RRBS数据可与转录组测序、组蛋白修饰测序、染色质开放区域检测等其他组学数据进行联合分析，从多个层面解析基因表达调控机制。例如，将启动子甲基化数据与基因表达数据进行关联分析，筛选甲基化调控的关键基因；将甲基化数据与组蛋白修饰数据进行整合，分析表观遗传修饰之间的协同或拮抗关系。

问：不同批次实验数据是否可比？

答：规范的实验流程和质量控制可以保证批次间数据的一致性。建议在不同批次实验中设置共同的对照样品，用于评估批次效应。数据分析时可通过标准化方法校正批次效应。对于大规模研究，建议在相同学条件下完成全部样品的建库和测序，减少技术变异对结果的影响。

问：数据分析周期需要多长时间？

答：RRBS项目周期通常为2-4周，具体时间取决于样品数量、测序方案和分析深度。实验建库约需5-7个工作日，上机测序约需2-3天，生物信息学分析约需3-5个工作日。对于大规模样品或深度分析需求，周期可能相应延长。