RRBS甲基化位点测定
CNAS认证
CMA认证
技术概述
RRBS甲基化位点测定是一种高效、精准的DNA甲基化检测技术,全称为Reduced Representation Bisulfite Sequencing,即简化代表性亚硫酸氢盐测序技术。该技术结合了限制性酶切、亚硫酸氢盐转化和高通量测序技术,能够以较低的成本实现对基因组中CpG岛区域的高精度甲基化分析。作为表观遗传学研究的重要手段,RRBS技术已成为癌症研究、发育生物学和精准医学领域不可或缺的检测方法。
DNA甲基化是最重要的表观遗传修饰之一,主要发生在CpG二核苷酸的胞嘧啶5号碳原子上。在人类基因组中,CpG位点分布并不均匀,大约60%的CpG位点聚集在基因组的CpG岛区域,这些区域通常位于基因的启动子区,与基因表达调控密切相关。异常的DNA甲基化模式与多种疾病的发生发展密切相关,特别是肿瘤抑制基因启动子区的超甲基化和全基因组的低甲基化,是癌症发生的重要特征。
RRBS技术的核心优势在于其"简化代表性"设计理念。通过使用甲基化不敏感的限制性内切酶(如MspI)对基因组DNA进行酶切,富集含有CpG岛的DNA片段,从而在不损失关键信息的前提下大幅降低测序成本。与传统全基因组甲基化测序(WGBS)相比,RRBS可以以约1/10的成本覆盖基因组中大部分功能重要的CpG位点,这使得大规模样本的甲基化研究成为可能。
该技术的发展源于科学家对甲基化检测成本效益的追求。2005年,Meissner等人首次提出了RRBS方法,此后经过多次技术优化,包括多重酶切策略、改进的文库构建流程和更精确的生物信息学分析算法,使得RRBS技术的准确性和可重复性不断提升。目前,RRBS已成为表观基因组学研究的主流技术平台之一。
- 高性价比:相比全基因组甲基化测序成本大幅降低
- 高覆盖率:可覆盖基因组中超过80%的CpG岛区域
- 单碱基分辨率:能够精确测定每个CpG位点的甲基化水平
- 可重复性强:技术重复性好,适合大规模队列研究
- 样本需求量低:仅需纳克级DNA即可完成检测
检测样品
RRBS甲基化位点测定对样品类型具有广泛的适用性,可检测多种来源的基因组DNA样品。了解各类样品的具体要求和处理方法,对于确保检测结果的准确性至关重要。
组织样品是RRBS检测中最常用的样品类型之一。新鲜冷冻组织是最佳选择,能够最大程度保留DNA的完整性和甲基化信息。对于临床样本,手术切除或活检获取的组织应在离体后尽快进行冷冻保存,推荐使用液氮速冻后转移至-80°C冰箱长期保存。石蜡包埋组织(FFPE)也可用于RRBS检测,但由于福尔马林固定过程可能导致DNA损伤和甲基化改变,需要特殊的DNA提取和修复流程。
血液样品是临床研究中最易获取的样品类型。外周血单个核细胞(PBMC)常用于表观遗传流行病学研究,采集后应使用抗凝管保存,尽快分离PBMC并冷冻保存。血清或血浆中的游离DNA(cfDNA)也可用于RRBS分析,在肿瘤液体活检领域具有重要应用价值。对于胎儿甲基化研究,羊水细胞和脐带血都是合适的样品来源。
细胞样品包括培养细胞系和原代细胞。贴壁细胞应在胰酶消化后离心收集,悬浮细胞直接离心收集。细胞沉淀应用PBS洗涤后冷冻保存。需注意,细胞培养过程中的应激状态可能影响甲基化水平,应在标准条件下培养并在合适时机收获细胞。
DNA样品可直接送检,但需满足特定质量要求。DNA应具有较高纯度,OD260/280比值应在1.8-2.0范围内,OD260/230比值应大于2.0。DNA完整性可通过琼脂糖凝胶电泳或安捷伦生物分析仪评估,主带应大于1kb。样品应避免反复冻融,推荐在TE缓冲液中-20°C或-80°C保存。
- 新鲜冷冻组织:推荐用量20-50mg,-80°C保存
- 石蜡包埋组织:需提供5-10张切片,厚度5-10μm
- 全血样品:采集量2-5mL,EDTA抗凝管收集
- 血清/血浆:采集量1-2mL,分离后-80°C保存
- 细胞沉淀:细胞数量不低于1×10^6个
- 纯化DNA:浓度不低于20ng/μL,总量不低于500ng
检测项目
RRBS甲基化位点测定涵盖多层次、多角度的甲基化分析内容,从单碱基水平的甲基化定量到基因组区域的甲基化模式分析,为研究者提供全面的表观遗传信息。
单CpG位点甲基化水平测定是最基础也是最重要的检测内容。RRBS技术能够在全基因组范围内测定数百万个CpG位点的甲基化比例,每个位点的甲基化水平以0-1之间的数值表示,0代表完全未甲基化,1代表完全甲基化。这种单碱基分辨率的定量分析能够揭示细微的甲基化变化,对于识别疾病相关的甲基化生物标志物具有重要价值。
启动子区域甲基化分析是功能研究的核心内容。基因启动子区的CpG岛甲基化状态与基因表达呈负相关,启动子高甲基化通常导致基因沉默。RRBS可对全基因组基因启动子区进行系统性甲基化分析,识别差异甲基化启动子,结合基因表达数据进行关联分析,阐明甲基化对基因表达的调控机制。
基因体甲基化分析提供基因功能区域的甲基化信息。包括第一外显子、外显子、内含子区域的甲基化水平测定。研究表明,基因体甲基化与基因表达呈正相关,可能参与转录延伸和剪接调控。RRBS数据可用于分析基因体甲基化与基因结构特征的关系。
差异甲基化区域(DMR)识别是比较研究的核心内容。通过比较不同组别样品的甲基化数据,采用统计学方法识别显著差异甲基化的基因组区域。DMR分析可用于筛选疾病特异性甲基化标志物,研究发育过程中的甲基化重编程,以及环境因素对甲基化的影响。
甲基化单核苷酸多态性(mSNP)分析是容易被忽视的重要内容。亚硫酸氢盐处理后,SNP位点可能发生C到T的转化,影响甲基化定量准确性。专业的RRBS分析流程应包含SNP位点识别和校正,确保甲基化测定的可靠性。
- 全基因组CpG位点甲基化水平测定
- 启动子CpG岛甲基化分析
- 基因体甲基化模式分析
- 增强子区域甲基化检测
- 差异甲基化区域(DMR)识别与注释
- 甲基化水平与基因表达关联分析
- 甲基化单倍型分析
- 重复序列区域甲基化检测
检测方法
RRBS甲基化位点测定采用标准化的实验流程和生物信息学分析流程,确保检测结果的准确性、可重复性和批次间一致性。
基因组DNA提取是检测的第一步。对于组织样品,采用酚-氯仿法或商业化DNA提取试剂盒进行提取,确保DNA纯度和完整性。对于石蜡包埋组织,需使用专门的FFPE DNA提取试剂盒,并进行DNA修复以逆转福尔马林引起的交联损伤。血液样品可采用全血DNA提取试剂盒直接提取。DNA提取后需进行质量检测,包括浓度测定、纯度评估和完整性分析。
限制性酶切是RRBS技术的关键步骤。使用MspI限制性内切酶对基因组DNA进行酶切,该酶识别CCGG序列,切割产生CpG富集的DNA片段。酶切反应在37°C进行,反应时间通常为过夜。酶切完成后进行琼脂糖凝胶电泳验证酶切效率,合格的酶切产物应呈现特征性的片段分布模式。
片段选择用于富集目标区域。酶切产物通过琼脂糖凝胶电泳分离,切取40-220bp范围的DNA片段进行回收纯化。这一步骤有效富集了CpG岛区域的DNA片段,同时去除了非目标区域的大片段DNA,是实现"简化代表性"的关键操作。
末端修复和接头连接为测序做准备。片段选择后的DNA进行末端修复,将粘性末端转化为平末端,随后在3'端添加单个腺嘌呤碱基(A尾),最后连接测序接头。接头连接后的产物进行纯化,去除未连接的接头和短片段。
亚硫酸氢盐转化是甲基化检测的核心步骤。使用亚硫酸氢钠处理DNA,将未甲基化的胞嘧啶转化为尿嘧啶,而甲基化的胞嘧啶保持不变。转化反应通常在变性条件下进行,反应温度约为50-55°C,时间约12-16小时。转化效率是影响检测结果准确性的关键因素,优化后的转化效率应大于99%。转化后进行DNA纯化和脱盐处理。
PCR扩增富集文库。转化后的DNA进行PCR扩增,使用特异性引物对连接了接头的片段进行扩增。PCR循环数需要优化,在保证足够产量的同时避免过度扩增造成的偏差。扩增产物进行纯化和质量检测,合格的文库应具有正确的片段大小分布和足够的浓度。
高通量测序采用Illumina平台进行。文库经质量检验后进行簇生成,然后在测序仪上进行双端测序。测序读长通常为150bp双端测序(PE150),测序深度根据研究目的确定,一般每个样品测序量约为5-10GB数据量。
生物信息学分析流程包括数据质控、比对、甲基化定量和差异分析等步骤。原始数据首先进行质量评估,去除低质量reads和接头序列。然后使用专门的比对软件(如Bismark、BSMAP)将数据比对到参考基因组。比对后的数据进行甲基化定量,计算每个CpG位点的甲基化比例。差异甲基化分析采用统计学方法比较不同组别间的甲基化差异。
- 基因组DNA提取与质量检测
- MspI限制性内切酶酶切
- 琼脂糖凝胶电泳片段选择(40-220bp)
- 末端修复、加A尾、接头连接
- 亚硫酸氢盐转化(未甲基化C→U)
- PCR扩增与文库纯化
- 文库质量评估(片段分布、浓度)
- Illumina高通量测序(PE150)
- 原始数据质控与过滤
- Bismark比对与甲基化定量
- 差异甲基化区域识别与功能注释
检测仪器
RRBS甲基化位点测定的顺利完成依赖于一系列精密仪器设备的支持,从样品制备到数据获取,每个环节都需要专业设备的保障。
DNA定量和质控设备是样品处理阶段的核心装备。NanoDrop超微量分光光度计用于快速测定DNA浓度和纯度,通过OD260、OD280和OD230数值评估DNA质量。Qubit荧光计采用荧光染料法进行DNA精确定量,不受RNA和蛋白质干扰。安捷伦2100生物分析仪或Fragment Analyzer用于DNA完整性分析,可生成数字化的电泳图谱,准确评估DNA片段大小分布。
样品前处理设备包括高速冷冻离心机、精密移液器、涡旋振荡器和恒温混匀仪等。高速冷冻离心机最大转速可达15000rpm以上,用于DNA沉淀和纯化过程中的离心分离。精密移液器覆盖0.1μL-1000μL量程,确保移液准确性。恒温混匀仪用于酶切反应和转化反应中的温控孵育,温度控制精度可达0.1°C。
凝胶电泳系统用于片段分离和质量检测。水平凝胶电泳系统用于常规DNA电泳分析,配备紫外凝胶成像系统进行结果观察。制备型聚丙烯酰胺凝胶电泳系统用于高精度片段选择,可实现10bp级别的片段分离。自动片段回收系统可提高片段选择效率和可重复性。
PCR扩增设备是文库构建的关键设备。梯度PCR仪可用于PCR条件优化,温度梯度范围通常为20-30°C。实时荧光定量PCR仪用于文库定量和质量评估,通过标准曲线法测定文库浓度。普通PCR仪满足常规PCR扩增需求,温度控制均匀性和升降温速率是关键性能指标。
高通量测序平台是数据获取的核心设备。Illumina NovaSeq系列是目前最常用的高通量测序平台,单次运行可产生数TB的测序数据。Illumina NextSeq系列适合中等规模的测序项目,运行周期短、灵活性高。Illumina MiSeq适用于小规模测序和快速验证。测序仪配备自动加样系统和流动池,实现高通量、自动化的测序流程。
生物信息学分服务器用于数据处理和存储。高性能计算服务器配备多核CPU和大容量内存,可并行处理多个样品的测序数据。存储系统提供TB级的数据存储空间,满足测序数据归档和分析需求。专业生物信息学软件包括FastQC、Trimmomatic、Bismark、 methylKit、DSS等,覆盖从数据质控到差异分析的完整流程。
- NanoDrop 2000/8000超微量分光光度计
- Qubit 4荧光计及DNA定量试剂盒
- Agilent 2100 Bioanalyzer生物分析仪
- Eppendorf 5424/5810高速冷冻离心机
- Bio-Rad T100梯度PCR仪
- ABI StepOnePlus实时荧光定量PCR仪
- Illumina NovaSeq 6000高通量测序系统
- Illumina NextSeq 550中通量测序系统
- 高性能生物信息学分析服务器
应用领域
RRBS甲基化位点测定在生命科学研究和医学应用中具有广泛的应用价值,为理解表观遗传调控机制和开发临床应用提供了强有力的技术支撑。
肿瘤表观遗传学研究是RRBS技术应用最广泛的领域。肿瘤发生发展过程中伴随大量异常的DNA甲基化改变,包括抑癌基因启动子超甲基化和全基因组低甲基化。通过RRBS技术系统分析肿瘤组织的甲基化谱,可识别肿瘤特异性甲基化标志物,阐明肿瘤发生机制。此外,循环肿瘤DNA(ctDNA)的甲基化分析为肿瘤早期诊断和疗效监测提供了新策略。
发育与分化研究利用RRBS技术揭示表观遗传重编程过程。哺乳动物发育过程中经历两次全基因组甲基化擦除和重建,即原始生殖细胞发育过程中的甲基化擦除和受精后植入前发育过程中的甲基化重编程。RRBS技术可用于分析不同发育阶段、不同组织器官的甲基化动态变化,理解细胞命运决定的表观遗传机制。
干细胞与再生医学研究依赖甲基化分析评估细胞状态。胚胎干细胞和多能干细胞具有独特的甲基化特征,与多能性维持密切相关。在干细胞定向分化过程中,甲基化谱发生变化以建立细胞类型特异性的表观遗传状态。RRBS可用于监测干细胞培养过程中的甲基化稳定性,评估分化效率和细胞纯度。
复杂疾病研究借助RRBS技术探索疾病的表观遗传基础。除肿瘤外,神经精神疾病、自身免疫病、代谢性疾病等复杂疾病均存在甲基化异常。全基因组关联研究(EWAS)通过比较患者和健康对照的甲基化谱,识别疾病相关的甲基化变异位点,为疾病机制研究和风险预测提供依据。
环境与健康研究利用RRBS评估环境因素的表观遗传效应。环境污染物、营养因素、心理应激等均可影响DNA甲基化状态。通过队列研究和动物模型,分析环境暴露与甲基化改变的关系,揭示环境因素致病的表观遗传机制,为预防医学提供科学依据。
精准医学应用将甲基化标志物转化为临床诊断工具。基于甲基化的肿瘤早筛产品已进入临床应用,通过检测血液中ctDNA的甲基化特征实现肿瘤早期筛查。药物基因组学研究表明,甲基化状态可预测药物疗效和毒副作用,指导个体化用药。甲基化标志物还可用于疾病预后评估和治疗反应监测。
植物与微生物研究也广泛应用RRBS技术。植物甲基化模式与动物存在差异,RRBS可用于分析植物发育、逆境响应的甲基化调控机制。微生物甲基化分析有助于理解表观遗传对基因表达的调控,在合成生物学和发酵工程中具有应用价值。
- 肿瘤甲基化标志物筛选与机制研究
- 肿瘤液体活检与早筛技术开发
- 胚胎发育与细胞分化研究
- 干细胞多能性与分化调控研究
- 神经精神疾病表观遗传机制研究
- 代谢性疾病甲基化关联分析
- 环境暴露与甲基化改变队列研究
- 药物甲基化生物标志物研究
- 植物表观遗传与逆境响应研究
常见问题
在进行RRBS甲基化位点测定过程中,研究人员可能遇到多种技术和应用相关的问题。以下针对常见问题进行详细解答,帮助研究者更好地理解和应用该技术。
问题一:RRBS与全基因组甲基化测序(WGBS)有何区别,应如何选择?
RRBS和WGBS都是基于亚硫酸氢盐转化和测序的甲基化检测技术,但在覆盖范围和成本上存在差异。WGBS覆盖全基因组所有CpG位点,包括基因间区和重复序列,但成本较高。RRBS通过酶切富集策略覆盖CpG富集区域,主要是启动子CpG岛,覆盖约10-15%的全基因组CpG位点,但成本仅为WGBS的1/10左右。如果研究目标是分析基因启动子区域甲基化或筛选功能重要的甲基化标志物,RRBS是性价比更高的选择;如果需要全面分析全基因组甲基化模式或研究重复序列甲基化,则应选择WGBS。
问题二:RRBS检测需要多少DNA量,样品质量有何要求?
RRBS检测通常需要500ng以上的基因组DNA,最低可接受200ng,但DNA量过低可能影响文库构建成功率。DNA纯度要求OD260/280比值在1.8-2.0范围内,OD260/230比值大于2.0。DNA完整性要求主带大于1kb,无明显降解。对于珍贵样品或低浓度样品,可采用优化的文库构建方案,如增加PCR循环数或使用DNA修复试剂。石蜡包埋组织DNA通常存在降解和损伤,需使用专门的修复方案。
问题三:亚硫酸氢盐转化效率如何评估,对结果有何影响?
亚硫酸氢盐转化效率是影响甲基化定量准确性的关键因素。转化效率可通过 spiked-in对照评估,添加已知甲基化状态的DNA片段作为内对照。转化效率低于99%会低估真实的甲基化水平,因为未转化的未甲基化胞嘧啶会被错误识别为甲基化胞嘧啶。优化转化条件包括延长反应时间、增加变性步骤、使用新鲜配制的转化试剂等。现代商业化转化试剂盒已能达到99%以上的转化效率。
问题四:如何确定合适的测序深度?
测序深度取决于研究目的和样品数量。对于差异甲基化分析,每个CpG位点需要足够的覆盖深度才能进行统计学可靠的定量。一般建议每个样品产生5-10GB数据量,对应约10-20×的平均覆盖深度。如果进行精细的等位基因特异性甲基化分析或检测低频甲基化变异,可能需要更高的测序深度。考虑到成本因素,大规模队列研究可在保证统计效力的前提下适当降低测序深度。
问题五:RRBS数据如何进行差异甲基化分析?
差异甲基化分析包括位点水平和区域水平两个层次。位点水平分析使用统计检验(如Fisher精确检验、逻辑回归)比较不同组别在每个CpG位点的甲基化比例差异,并进行多重检验校正。区域水平分析(DMR识别)考虑基因组位置相邻的CpG位点的甲基化模式,采用滑动窗口或分割算法识别差异甲基化区域。常用的分析软件包括methylKit、DSS、bsseq等,建议使用两种以上方法进行交叉验证,提高结果可靠性。
问题六:甲基化数据与基因表达如何关联分析?
甲基化与基因表达关联分析需要同时获取甲基化数据和转录组数据。最常用的策略是将启动子区域甲基化与基因表达进行相关性分析,一般认为启动子高甲基化与基因低表达相关。关联分析方法包括Pearson或Spearman相关性检验、线性回归模型等。需要注意的是,甲基化与基因表达的关系并非简单的线性关系,还受到转录因子结合、染色质状态等多种因素影响。建议结合功能验证实验确认甲基化对基因表达的调控作用。
问题七:不同批次RRBS实验结果如何保证一致性?
批次效应是高通量测序研究中的常见问题。保证批次间一致性的措施包括:使用统一的实验方案和试剂批次;设置批次内对照样品,如商业化对照DNA;采用随机化设计,避免批次与实验组别混淆;数据归一化处理去除批次效应。在数据分析阶段,可使用ComBat等算法进行批次效应校正。对于大规模研究,建议将所有样品分成若干批次完成,每批次包含各类别样品,便于后续校正。
问题八:RRBS检测结果是否可以进行跨研究比较?
不同研究间的RRBS结果比较需要考虑多种因素。首先,实验方案差异可能导致覆盖位点不同,因为酶切效率和片段选择策略会影响最终的CpG覆盖范围。其次,分析流程差异(如比对软件、甲基化定量方法)会影响甲基化数值。建议在跨研究比较时,重点关注共同覆盖的CpG位点,使用标准化的数据处理流程,并报告详细的方法学信息。大型数据库如TCGA提供的甲基化数据可作为参考标准。
- RRBS与WGBS技术选择需根据研究目的和预算确定
- DNA样品质量和数量直接影响检测成功率
- 亚硫酸氢盐转化效率需达到99%以上确保结果准确
- 测序深度根据研究目的和统计效力要求确定
- 差异甲基化分析建议使用多种方法交叉验证
- 甲基化-表达关联分析需结合功能实验验证
- 批次效应可通过实验设计和数据校正控制
- 跨研究比较需考虑实验方案和分析流程差异