CRISPR文库筛选脱靶gRNA活性验证

来源:CRISPR文库筛选脱靶gRNA活性验证(DNA Lab Space 实验教程)| 编译:DNA Lab Space

📎 相关工具:代谢通路设计

筛完库拿到几十条候选gRNA,一验脱靶就出事:预测高分位点实际切不动,报告基因喊阳性、全基因组测序却找不到影子。钱花了,细胞耗了,数据还不敢写进文章。问题多半不在CRISPR本身,而在验证链条的某个环节漏了。本文按样本准备、切割富集、建库测序、数据分析四段走,把每一步的典型故障、判断方法和解决参数摆出来,你对着查就行。

第一步:先排查这几个低级错误

别急着怀疑生物学。新手翻车,八成是下面几条。

1. gDNA投入量不够。 GUIDE-seq类实验每样本至少500 ng—1 µg gDNA,amplicon验证也要100 ng;低于这个量,覆盖度掉下去,脱靶位点直接假阴性。解决:重新提gDNA,Qubit定量,别拿Nanodrop的RNA污染值凑数。

2. UMI忘了加或者加错位置。 建库时5'端接8—10 nt随机UMI,否则PCR重复全被当成独立位点,脱靶数量虚高3—10倍。解决:接头合成时把UMI写进引物,建库后按UMI去重再计数。

3. 引物二聚体没查。 接头引物和gRNA spacer互补区容易二聚,跑2%琼脂糖胶看到弥散条带就是它。解决:用引物设计工具重设,Tm控制在60—65°C,二聚体ΔG>-9 kcal/mol。

4. 细胞传代太多。 文库筛选用细胞别超过p15,传代多了gRNA representation掉一半,背景噪声压不住。解决:电转后24小时冻一份备份,后续验证全用同一批。

教程配图2
▲ 教程配图2

第二步:分维度系统排查

模块一:样本与文库准备

1. gDNA里残留gRNA质粒

现象:脱靶位点富集在质粒骨架附近,而非基因组真实位点。

原理:质粒DNA污染gDNA,PCR直接把gRNA表达框扩出来,测序读数全落在预期外的“假位点”。

解决:gDNA提取后加DpnI 37°C消化1小时,或者用Hirt法去质粒;质粒残留高的样本直接重提。

验证:qPCR测gRNA表达框拷贝数,正常应<0.1拷贝/细胞。

2. 文库覆盖度不足

现象:对照样本里部分gRNA读数<50,组间比较全是噪声。

原理:电转效率低或MOI过高,多个gRNA进同一个细胞,representation塌缩。

解决:MOI控制在0.1—0.3,电转后24小时取1×10⁶细胞提gDNA,NGS测覆盖度要求≥500×。

验证:看gRNA分布Gini系数,<0.2才算合格。

模块二:切割与富集

3. Cas9活性不够,靶向位点都切不动

现象:靶位点indel<20%,脱靶信号更检测不到。

原理:RNP比例不对或Cas9蛋白活性下降,切割效率低,末端标记自然少。

解决:RNP电转用Cas9:gRNA摩尔比1:1.2,室温孵育10分钟再电转;Cas9蛋白别反复冻融,分装后-80°C保存。

验证:靶位点做amplicon NGS,indel≥30%才继续往下做。切不动靶点,谈脱靶没意义。

4. GUIDE-seq接头连接效率低

现象:双链断裂末端标记少,阳性对照gRNA也检不出几个位点。

原理:末端修复不充分,或者接头浓度不够,连接酶没干活。

解决:37°C末端补平30分钟,接头用10 µM,16°C连接过夜;连接后AMPure XP 1.8×纯化,别省这步。

验证:跑阳性对照gRNA(EMX1或HEK4),应检出≥50个脱靶位点,否则整批重做。

图3:CRISPR文库筛选脱靶gRNA活性验证排查流程示意图
▲ 图3:CRISPR文库筛选脱靶gRNA活性验证排查流程示意图

模块三:PCR与测序建库

5. PCR循环数过多

现象:重复序列、同源基因家族位点假阳性暴增,Sanger验证全阴。

原理:指数扩增把痕量污染放大,背景噪声盖过真信号。

解决:循环数压到18—22,用高保真酶(KAPA HiFi),模板100 ng,每个样本至少3个技术重复。

验证:做无模板对照(NTC),NTC读数应<总读数0.5%。

6. 测序深度不够

现象:脱靶位点reads<10,统计上不显著,但qPCR明明有富集。

原理:PE150测10M reads只够低复杂度分析,基因组范围找位点远远不够。

解决:每样本≥30M reads,双端150 bp;文库复杂度不够就加测,别靠降阈值硬凑。

验证:掺入已知阳性位点做spike-in,回收率应>80%。

模块四:数据分析

7. 阈值定太松

现象:候选脱靶位点几百个,挑10个做Sanger验证全阴性。

原理:没做背景噪声建模,泊松分布没拟合,fold change一刀切。

解决:用对照样本建背景模型,p<0.01且fold change>3才进候选;多重比对位点单独标记。

验证:随机挑10个位点做amplicon NGS,真阳性率>70%才算阈值合理。

第三步:进阶疑难问题

8. 多重比对漏检

脱靶位点落在Alu或同源基因家族,BWA默认只报一条比对,读长又短,直接漏掉。现象是CIRCLE-seq体外找得到、细胞内GUIDE-seq找不到。解决:用minimap2 -x sr或BWA-MEM加-a参数输出全部比对,配合UMI去重。验证:对候选位点做位点特异性PCR加长读长测序,确认是不是同源区域错配。

9. 克隆性富集假阳性

某个脱靶位点在所有重复里都高,但细胞群体水平验证切不动。原理是建库前细胞已发生克隆扩增,一个克隆占了30%群体。解决:电转后72小时内收样,别超过96小时;每次独立电转做3个生物学重复。验证:数字PCR测该位点编辑效率,同时看靶向位点是否也异常富集——两者同高,大概率是克隆效应。

避坑总结+实操建议

三条核心结论:

  • 对照比样本多。 NTC、阳性gRNA、非靶向gRNA三件套缺一不可,少一个都没法判断背景。
  • UMI加去重是刚需。 不是可选项,建库时就设计进去,事后补不了。
  • 先过靶向切割这关。 靶位点indel≥30%再谈脱靶,target切不动,off-target全是噪声。

两条日常好习惯:

  • 每次电转后24小时留一份gDNA冻-80°C,后面复测不用重做细胞,省一周时间。
  • 建库前跑一次Bioanalyzer,片段分布不对直接重做,别硬上测序——省下的测序费不够返工。

工具衔接

验证引物设计别手算,直接用站内引物设计工具,输入位点序列自动出Tm和二聚体评估;做Cas9蛋白表达优化时,密码子优化工具能按HEK293或大肠杆菌偏好改序列,蛋白可溶性优化工具则帮你判断表达后是不是全在包涵体里。三个工具串起来,从gRNA验证到蛋白制备一条线走完,少踩很多坑。

延伸阅读:更多代谢通路设计、多基因组装等合成生物学工具,可访问 DNA Lab Space。

📎 相关文章

qPCR熔解曲线多峰退火温度优化 10-27 细胞复苏存活率低冻存液渗透压 10-27 细胞转染毒性大培养基更换时机 10-27 细胞划痕实验愈合过快血清浓度降低 10-27