技术概述

RRBS（Reduced Representation Bisulfite Sequencing）即简化代表性亚硫酸氢盐测序，是一种结合了限制性酶切消化和亚硫酸氢盐转换的高通量甲基化检测技术。该技术由Meissner等人于2005年首次开发，旨在以较低的成本实现全基因组范围内CpG岛区域DNA甲基化状态的高精度分析。与全基因组亚硫酸氢盐测序（WGBS）相比，RRBS通过限制性内切酶富集CpG富集区域，在保证数据质量的前提下大幅降低了测序成本，是目前研究DNA甲基化差异的主流技术手段之一。

DNA甲基化是表观遗传学研究中最为重要的修饰形式之一，主要发生在胞嘧啶的5号碳原子上，形成5-甲基胞嘧啶（5mC）。在哺乳动物基因组中，DNA甲基化主要出现在CpG二核苷酸位点。CpG岛通常位于基因的启动子区域，其甲基化状态与基因表达密切相关。当启动子区域的CpG岛发生高甲基化时，通常会抑制基因转录；而低甲基化则有助于基因的正常表达。因此，精准检测CpG岛的甲基化水平对于理解基因表达调控机制具有重要意义。

RRBS技术的核心优势在于其能够针对基因组中甲基化密度较高的区域进行定向富集。通过使用MspI等限制性内切酶（识别序列为CCGG），RRBS可以特异性地切割富含CpG位点的DNA片段。经过酶切、末端修复、加A尾、连接接头等步骤后，利用亚硫酸氢盐处理将未甲基化的胞嘧啶转化为尿嘧啶，而甲基化的胞嘧啶保持不变。最终通过高通量测序，可以精确检测每个CpG位点的甲基化比例。

甲基化差异分析是RRBS技术的核心应用场景。通过对不同样本组别（如疾病组与对照组、处理组与空白组等）进行RRBS测序，可以获得各CpG位点的甲基化水平数据。随后利用生物信息学方法进行差异甲基化区域（DMR）和差异甲基化位点（DMP）的鉴定，从而揭示表观遗传层面的生物学差异。这种分析策略已被广泛应用于肿瘤研究、发育生物学、环境暴露评估等多个领域。

检测样品

RRBS甲基化差异分析对样品类型具有较广的适用性，可满足多种研究需求。以下是常见的检测样品类型：

  • 动物组织样品：包括正常组织、病变组织、肿瘤组织等。取样后应立即液氮速冻，-80°C保存。建议取样量不低于50mg，以保证足够的DNA提取量。
  • 植物组织样品：如叶片、根、茎、种子等。植物细胞具有细胞壁结构，DNA提取需采用专门的植物基因组DNA提取方法。建议取样量不低于100mg。
  • 血液样品：包括全血、血清、血浆、外周血单核细胞（PBMC）等。全血样品建议使用抗凝管采集，-80°C保存；血浆和血清样品建议分离后-80°C保存。
  • 培养细胞样品：贴壁细胞或悬浮细胞均可。收集细胞后用PBS洗涤，离心收集细胞沉淀，液氮速冻后-80°C保存。建议细胞数量不低于1×10^6个。
  • 微生物样品：细菌、真菌等微生物样品也可用于RRBS分析，但需注意微生物基因组中CpG分布特征与真核生物的差异。
  • 甲醛固定石蜡包埋（FFPE）样品：临床病理标本可通过FFPE形式保存。由于甲醛固定可能导致DNA交联和降解，需采用专门的FFPE DNA提取试剂盒。

样品质量是影响RRBS实验成功的关键因素。送检样品应满足以下基本要求：DNA浓度不低于50ng/μL，DNA总量不低于1μg；DNA完整性良好，主条带清晰，无明显降解；OD260/OD280比值在1.8-2.0之间；样品无RNA污染、无蛋白质残留。对于DNA质量较差的样品，建议先进行纯化处理或重新提取。

检测项目

RRBS甲基化差异分析涵盖多个层面的检测内容，可全面解析样品的表观遗传特征：

  • 全基因组CpG岛甲基化检测：覆盖基因组中绝大部分CpG岛区域，可检测启动子区域、基因间区、重复序列等区域的甲基化状态。
  • 差异甲基化位点（DMP）分析：鉴定不同组别间甲基化水平存在显著差异的单个CpG位点。通常设定阈值如|Δβ|>0.2和FDR<0.05进行筛选。
  • 差异甲基化区域（DMR）分析：将相邻的差异甲基化位点合并为区域进行分析。DMR分析可以更准确地反映功能区域的甲基化变化，常用的鉴定工具包括DSS、methylKit、Metilene等。
  • 启动子甲基化分析：聚焦于基因启动子区域（通常定义为转录起始位点上游2kb至下游500bp范围），分析启动子甲基化与基因表达的关系。
  • 基因本体（GO）功能富集分析：对差异甲基化基因进行GO分类注释，分析其在生物学过程、细胞组分、分子功能等层面的分布特征。
  • KEGG通路富集分析：将差异甲基化基因映射到KEGG代谢通路中，挖掘受甲基化调控的关键信号通路。
  • 转录因子结合位点甲基化分析：分析转录因子结合位点区域的甲基化变化，探讨甲基化对转录因子结合能力的影响。
  • 重复序列甲基化分析：检测LINE、SINE、LTR等重复元件的甲基化状态，评估基因组稳定性。

检测项目的选择应根据研究目的进行合理搭配。对于探索性研究，建议进行全面的甲基化差异分析；对于目标基因验证，可针对性地进行特定位点的深度分析。

检测方法

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

第一阶段：DNA提取与质量检测

采用适当的DNA提取方法从样品中分离基因组DNA。动物组织常用酚-氯仿法或硅胶膜吸附柱法；植物组织采用CTAB法或商用植物DNA提取试剂盒；血液样品可采用离心柱法或磁珠法。提取后的DNA需进行浓度测定（Qubit荧光定量）、纯度检测（NanoDrop分光光度计）和完整性评估（琼脂糖凝胶电泳或Agilent Bioanalyzer）。合格的DNA样品进入下游实验流程。

第二阶段：RRBS文库构建

  • 限制性酶切：使用MspI限制性内切酶消化基因组DNA，该酶识别CCGG序列，能够富集CpG密集区域。酶切产物经琼脂糖凝胶电泳分离，回收40-220bp的目标片段。
  • 末端修复与加A尾：使用T4 DNA聚合酶和Klenow片段对酶切片段进行末端修复，生成平末端；随后利用Klenow片段（3'-5'外切酶活性缺失）在3'末端添加单个腺嘌呤碱基。
  • 接头连接：使用T4 DNA连接酶将甲基化接头连接至DNA片段。接头经甲基化修饰，可抵抗后续的亚硫酸氢盐处理。
  • 亚硫酸氢盐转换：使用亚硫酸氢盐处理DNA样品，将未甲基化的胞嘧啶脱氨转化为尿嘧啶，而甲基化的胞嘧啶保持不变。此步骤是甲基化检测的关键，需优化反应条件以提高转换效率并降低DNA降解。
  • PCR扩增：使用高保真聚合酶对转换后的文库进行PCR扩增，引入测序接头和样品标签。扩增循环数需优化以避免过度扩增造成的偏倚。
  • 文库质检：使用Agilent Bioanalyzer或TapeStation检测文库片段大小分布，使用qPCR定量文库浓度。

第三阶段：高通量测序

质检合格的文库经 Pooling后进行上机测序。根据样品数量和检测深度需求，选择适当的测序平台。常用平台包括Illumina NovaSeq、HiSeq、MiSeq等。推荐采用双端测序模式（PE150），以获得更高的比对效率。测序数据量根据研究目的确定，通常每个样品5-10M reads可满足大部分分析需求。

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

  • 数据质控：使用FastQC评估原始测序数据的质量，使用Trimmomatic或Cutadapt去除低质量碱基和接头序列。
  • 序列比对：使用Bismark、BSMAP或BS-Seeker等软件将处理后的序列比对至参考基因组。比对前需对参考基因组进行亚硫酸氢盐转换模拟。
  • 甲基化 calling：根据比对结果，统计每个CpG位点的甲基化与非甲基化read数，计算甲基化比例。
  • 差异分析：使用methylKit、DSS、DSS2等软件包进行差异甲基化位点（DMP）和差异甲基化区域（DMR）的鉴定。
  • 功能注释：对DMR相关基因进行GO和KEGG功能富集分析，挖掘生物学意义。
  • 可视化展示：生成甲基化水平热图、Circos图、曼哈顿图、火山图等多种可视化图形。

检测仪器

RRBS甲基化差异分析涉及多类仪器设备，确保实验的精确性和数据的可靠性：

核酸分析设备：

  • Qubit荧光定量仪：用于DNA样品的精确定量，采用荧光染料结合双链DNA进行检测，灵敏度高，准确性好。
  • NanoDrop分光光度计：快速测定DNA浓度和纯度，OD260/OD280和OD260/OD230比值反映样品纯度。
  • Agilent Bioanalyzer 2100：微流控芯片电泳系统，用于DNA片段大小分布分析和文库质检。
  • Agilent TapeStation：自动化电泳系统，高通量检测文库质量和片段分布。

文库构建设备：

  • PCR扩增仪：用于文库扩增、接头连接等温控反应环节。常用型号包括Bio-Rad T100、Applied Biosystems Veriti等。
  • QIAcube全自动核酸提取仪：自动化完成DNA提取和纯化流程，减少人工操作引入的误差。
  • MAGBIND提取系统：基于磁珠法的自动化提取平台，适用于高通量样品处理。

测序平台：

  • Illumina NovaSeq 6000：高通量测序平台，单次运行可产生数十 Tb数据，适用于大规模样本检测。
  • Illumina HiSeq X Ten：高性能测序系统，适用于批量样本的全基因组甲基化分析。
  • Illumina MiSeq：中小通量测序平台，适用于小规模样本或目标区域甲基化检测。

辅助设备：

  • 台式高速离心机：用于样品离心、沉淀收集等操作。
  • 超低温冰箱： -80°C保存DNA样品和文库，确保样品稳定性。
  • 液氮罐：用于组织样品的长期保存。
  • 生物安全柜：提供洁净操作环境，防止外源污染。

应用领域

RRBS甲基化差异分析在生命科学和医学研究中具有广泛的应用价值：

肿瘤学研究：

DNA甲基化异常是肿瘤发生发展的重要特征之一。RRBS可用于鉴定肿瘤特异性甲基化标志物，分析肿瘤组织与正常组织的表观遗传差异，探索甲基化在肿瘤诊断、分型、预后评估中的应用价值。例如，通过比较肿瘤组织和癌旁组织的甲基化谱，可以发现潜在的肿瘤抑制基因启动子高甲基化事件，揭示肿瘤发生的表观遗传机制。

发育生物学研究：

DNA甲基化在胚胎发育、细胞分化过程中发挥关键调控作用。RRBS可用于分析不同发育阶段、不同组织器官的甲基化动态变化，揭示表观遗传重编程的规律。在干细胞研究中，通过比较诱导多能干细胞与胚胎干细胞的甲基化特征，可以评估重编程的效率和完整性。

复杂疾病研究：

除肿瘤外，DNA甲基化异常与心血管疾病、糖尿病、神经系统疾病等多种复杂疾病密切相关。RRBS可用于分析疾病状态下全基因组甲基化变化，鉴定疾病相关的功能性甲基化位点，为疾病的表观遗传治疗提供靶点。

环境暴露与表观遗传毒性评估：

环境污染物、化学毒物等外源因素可引起DNA甲基化改变。RRBS可用于评估环境暴露后的表观遗传效应，筛选敏感的甲基化标志物，建立暴露-效应关系。在职业卫生领域，可用于评估工人长期暴露于有害因素后的表观遗传损伤。

植物表观遗传学研究：

植物DNA甲基化在基因印记、转基因沉默、胁迫响应等方面发挥重要作用。RRBS可用于分析不同品种、不同组织、不同生长条件下的植物甲基化谱，揭示甲基化在植物发育和环境适应中的调控作用。

药物表观遗传毒性评价：

药物研发过程中需评估化合物对表观遗传的潜在影响。RRBS可用于检测药物处理后的甲基化变化，评价药物的表观遗传毒性，为药物安全性评估提供科学依据。

衰老研究：

衰老过程中DNA甲基化呈现特定变化模式，被称为"表观遗传时钟"。RRBS可用于分析不同年龄个体的甲基化特征，研究衰老相关的甲基化变化，探索延缓衰老的干预策略。

常见问题

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

RRBS和WGBS都是基于亚硫酸氢盐转换的甲基化检测技术，主要区别在于覆盖范围和成本。WGBS覆盖全基因组范围的所有CpG位点，测序深度均匀，但成本较高；RRBS通过酶切富集CpG密集区域（主要是CpG岛），覆盖约10-15%的基因组区域和60-80%的CpG岛，成本约为WGBS的1/5-1/10。对于聚焦启动子甲基化的研究，RRBS是性价比更高的选择。

问：RRBS检测的准确性如何？

RRBS检测准确性主要取决于亚硫酸氢盐转换效率和测序深度。现代试剂盒的亚硫酸氢盐转换效率可达99%以上，通过λDNA或pUC19质粒作为转换对照可监控转换效率。测序深度方面，每个CpG位点10-30×覆盖度即可获得稳定的甲基化定量结果。整体而言，RRBS的检测准确性与WGBS相当，单碱基分辨率可精确反映甲基化状态。

问：哪些因素会影响RRBS实验结果？

主要影响因素包括：（1）DNA质量：降解、污染、低浓度都会影响文库构建效率；（2）酶切效率：不完全酶切会导致覆盖度降低和偏倚；（3）亚硫酸氢盐处理条件：温度、时间、pH值等参数需优化；（4）PCR扩增：过度扩增会导致重复率升高和偏倚；（5）测序深度：过低深度会影响定量准确性。

问：RRBS可用于哪些物种的甲基化分析？

RRBS适用于所有具有参考基因组序列的真核生物。常用模式生物包括人、小鼠、大鼠、斑马鱼、线虫、拟南芥、玉米、水稻等。对于基因组CpG密度较低的物种（如某些昆虫和真菌），RRBS的覆盖效果可能受限，需评估后选择合适的技术策略。

问：样品重复数量有何要求？

为保证统计学效力，建议每组样品设置至少3-5个生物学重复。对于人类临床样品，考虑到个体差异，建议每组不少于10个样品。样品数量过少会导致假阴性率升高，难以鉴定出真实的差异甲基化位点。

问：如何解决批次效应问题？

批次效应是高通量测序中的常见问题。建议：（1）实验设计时将不同组别样品随机分配至不同批次；（2）文库构建时保持操作一致性；（3）数据分析时使用ComBat、removeBatchEffect等方法进行批次效应校正；（4）设置批次对照样品监控技术变异。

问：RRBS数据可以与其他组学数据进行整合分析吗？

可以。RRBS甲基化数据常与转录组测序（RNA-seq）数据整合，分析启动子甲基化与基因表达的相关性；也可与ChIP-seq、ATAC-seq等表观基因组数据进行整合，构建多层次的基因调控网络。整合分析有助于全面理解表观遗传调控机制。

问：差异甲基化结果如何验证？

常用的验证方法包括：焦磷酸测序、甲基化特异性PCR（MSP）、甲基化特异性定量PCR（MS-qPCR）、靶向亚硫酸氢盐测序（Targeted BS-seq）等。焦磷酸测序是金标准验证方法，可定量检测特定区域的甲基化水平。