4. sgRNA设计与脱靶预测

sgRNA是CRISPR敲除实验的“导航序列”,它的质量直接决定Cas9能否在目标位点高效切割,以及是否会在非目标位点造成意外编辑。一个完整的sgRNA设计流程,本质上是在“高切割效率”和“低脱靶风险”之间寻找平衡点。本节围绕三个核心问题展开:如何用在线工具完成sgRNA的批量设计与筛选、脱靶评分算法在评估什么、以及如何通过实验验证和文库构建把设计落到可执行的实验方案上。

一、sgRNA的序列约束与PAM选择

以最常用的SpCas9为例,sgRNA的20 nt引导序列必须紧邻一个NGG型PAM(N为任意碱基)。设计时通常优先选择位于基因编码区前2/3的外显子,且尽量落在保守结构域上游,这样移码突变更容易破坏蛋白功能。需要避开的区域包括:外显子边界附近(可能产生外显子跳跃)、存在常见SNP的位点(不同细胞系结果不一致)、以及高度重复或同聚物区域(易造成测序和比对错误)。

实际操作中,一个基因通常先产生10~30条候选sgRNA,再通过评分和脱靶预测收敛到2~4条进入湿实验。候选序列的GC含量建议控制在40%~70%,避免连续4个以上相同碱基(如TTTT会触发U6启动子终止信号)。

二、三大在线工具的使用定位

CRISPOR、Benchling和CHOPCHOP是当前使用最广的三类设计平台,它们的定位和输出侧重不同,配合使用效果最好。

  • CRISPOR:输入基因名或序列后,会同时给出多个评分算法(MIT、CFD、Doench 2016等)的结果,并列出每个sgRNA在全基因组范围内的脱靶位点、错配数量和位置。它的优势是脱靶信息透明,适合做候选序列的横向比较。
  • Benchling:在序列编辑器中直接标注PAM和候选sgRNA,支持批量设计、自动去重和引物设计,适合已有质粒图谱或需要与克隆流程衔接的场景。
  • CHOPCHOP:界面简洁,支持多物种和多种Cas9变体(如Cas12a、双切口酶),会给出靶向效率评分和脱靶数量,适合快速初筛。

使用时的通用流程是:在CRISPOR中批量输入候选序列,导出包含CFD评分、MIT评分和脱靶位点数的表格;在Benchling中核对序列是否落在正确外显子并检查SNP;最后用CHOPCHOP交叉验证评分趋势。若三个工具对同一条sgRNA的评分方向一致,可信度较高;若分歧明显,应优先选择脱靶位点少且错配集中在PAM远端(第1~12位)的序列。

三、脱靶评分算法在算什么

脱靶评分的核心逻辑是:将候选sgRNA与基因组中所有可能位点比对,统计错配的数量、位置和分布,再映射为一个风险分数。不同算法的差异主要体现在权重设计上。

MIT评分采用位置权重:靠近PAM的错配(第13~20位)惩罚更重,因为该区域对切割活性影响更大;而PAM远端的错配容忍度较高。CFD评分(Doench 2016)则基于大规模实验数据训练,对每个错配位置和碱基类型给出经验权重,更贴近真实切割效率。简单来说,MIT偏向保守,CFD更接近实验观测。

需要特别注意的是,评分低不等于绝对安全。某些sgRNA在评分上表现良好,但仍可能在同源序列或假基因上产生切割。因此,评分只能作为筛选依据,不能替代实验验证。

sgRNA设计与脱靶预测流程 候选序列生成 外显子前2/3 · NGG PAM GC 40%~70% 在线工具评分 CRISPOR · Benchling CHOPCHOP 脱靶评分 MIT · CFD · Doench 错配位置权重 候选筛选 2~4条进入 湿实验 脱靶评分关键维度 错配数量 0~1个错配:高风险 2~3个错配:中风险 ≥4个错配:低风险 错配位置 PAM近端(13~20):高惩罚 PAM远端(1~12):低惩罚 种子区(1~8):中等权重 碱基类型 G:T / G:A 错配容忍度低 C:T / A:G 容忍度较高 CFD基于实验数据训练 实验验证策略 T7E1 / Surveyor 酶切 PCR扩增靶区域 变性退火形成错配 酶切后电泳检测切割带 切割效率≈ indel 比例 Sanger测序 + TIDE PCR产物直接测序 TIDE分析混合峰 定量indel频率 适合单克隆前评估 脱靶位点深度测序 预测前5~10个脱靶位点 设计引物扩增 NGS检测编辑频率 判断是否需换sgRNA

四、实验验证策略

设计完成后,必须通过实验确认sgRNA的实际切割效率。常用方法有三类:

T7E1/Surveyor酶切法是最经济的初筛手段。取靶细胞基因组DNA,PCR扩增包含靶位点的片段(通常400~800 bp),将PCR产物变性后缓慢退火,若存在编辑则形成错配,T7E1酶识别并切割错配位点,电泳后出现两条较小的切割带。切割带灰度占比可粗略估算indel效率,一般切割效率超过20%的sgRNA值得继续推进。

Sanger测序+TIDE分析适合需要定量结果的场景。将PCR产物直接测序,TIDE软件通过解析混合峰形,可计算出indel频率和类型。它的优势是不需要酶切步骤,且能区分不同编辑类型,但要求靶区域测序质量良好。

脱靶位点深度测序是评估特异性的关键。从CRISPOR预测结果中选取评分最高的5~10个脱靶位点,设计引物进行扩增,用NGS检测这些位点的编辑频率。若某个脱靶位点的编辑频率超过1%,且该位点位于基因功能区,则应考虑更换sgRNA或改用高保真Cas9变体(如eSpCas9、HiFi Cas9)。

五、sgRNA文库构建要点

当需要进行全基因组筛选时,需要构建sgRNA文库。文库构建的核心是保证覆盖度和均一性。以人全基因组敲除文库为例,通常每个基因设计4~6条sgRNA,总库容量约10万条。构建流程包括:

  1. 将sgRNA序列合成在oligo pool中,两端包含与载体同源的接头;
  2. 通过Gibson组装或限制性内切酶克隆将sgRNA插入慢病毒载体(如lentiCRISPRv2);
  3. 电转化感受态细胞,获得足够数量的克隆(通常要求覆盖度≥100×,即10万条sgRNA至少需要1000万个转化子);
  4. 提取质粒,通过NGS验证文库的覆盖度和sgRNA分布均一性。

文库质量的关键指标是:未出现明显缺失的sgRNA(缺失率<5%)、各sgRNA丰度分布均匀(Gini指数较低)、以及无过多空载体或接头异常。若文库覆盖度不足,后续筛选的假阴性率会显著升高。

实操建议:对于常规单基因敲除,优先选择CRISPOR中CFD评分≥0.6且脱靶位点数≤3的sgRNA;若三个工具评分分歧较大,以CFD评分为主,因为它基于实验数据训练,更接近真实切割效率。进入湿实验前,至少保留2条独立sgRNA作为备份,避免单条sgRNA效率不佳导致实验停滞。
常见误区:脱靶评分低并不等于没有脱靶。评分算法基于序列比对和统计模型,无法完全预测染色质可及性、DNA甲基化等表观因素对切割的影响。因此,对于关键实验,脱靶位点的NGS验证不可省略。