4. 靶点选择与sgRNA设计实战

在CRISPR动物细胞编辑实验中，sgRNA的质量直接决定切割效率与脱靶风险。一个合格的sgRNA设计流程，通常包含四个关键环节：靶基因外显子定位、候选sgRNA序列枚举、脱靶预测与评分、以及最终筛选与实验验证。本节以小鼠Pten基因（Gene ID: 19211）为例，完整走一遍从序列检索到候选sgRNA排序的实战过程。

一、靶基因外显子定位：先确定"切哪里"

CRISPR-Cas9的切割位点位于PAM（NGG）上游3 bp处。若目标是敲除基因功能，应优先选择位于功能结构域编码区且靠近5'端的组成型外显子，这样移码突变更容易破坏蛋白功能。若目标是敲入或点突变，则应选择靠近目标位点的外显子。

以Pten为例，其编码区共9个外显子，其中外显子5编码磷酸酶催化核心基序（HCXXGXXR）。通过NCBI Gene或Ensembl检索可获得外显子坐标。假设从Ensembl获取到外显子5的基因组坐标为chr19:32,758,000–32,758,180（示例坐标，仅用于教学演示），长度为181 bp。我们需要提取该区间的参考序列，并同时获取其反向互补序列，因为sgRNA可能位于正链或负链。

实际操作中，可用UCSC Genome Browser的"Get DNA"功能导出序列，或使用命令行工具samtools faidx从参考基因组中提取：

samtools faidx GRCm39.fa chr19:32758000-32758180
# 输出正链序列
samtools faidx GRCm39.fa chr19:32758000-32758180 | \
  rev | tr ACGTacgt TGCAtgca | rev
# 输出反向互补序列

将两条序列分别保存为exon5_fwd.txt和exon5_rev.txt，作为后续sgRNA枚举的输入。

二、候选sgRNA枚举：从序列到20 nt

对于每条序列，扫描所有符合"20 nt + NGG"模式的片段。以正链为例，从第1个碱基开始，每滑动1个碱基取20 nt，检查其下游3 nt是否为NGG。若符合，则记录该20 nt序列、PAM类型（NGG/NAG）、链方向及起始位置。

假设在exon5_fwd中枚举出12条候选，在exon5_rev中枚举出9条，共21条候选sgRNA。此时需要先做一轮基础过滤：

• GC含量在40%–70%之间，避免过高导致二级结构、过低导致结合不稳定；
• 排除连续4个以上相同碱基（如TTTT），以免造成转录终止或合成困难；
• 排除含常见限制性酶切位点（如BbsI的GAAGAC）的序列，便于后续载体克隆；
• 若使用U6启动子，5'端第一位碱基优选G，若为其他碱基可考虑加G或更换启动子。

经过滤后，假设剩余14条候选进入脱靶预测。

三、脱靶预测工具使用：从序列到风险评分

脱靶预测的核心是：将候选sgRNA序列与全基因组进行比对，允许一定数量的错配（通常≤4），并评估错配位置与PAM邻近程度对切割活性的影响。常用工具包括CRISPOR、Cas-OFFinder、Benchling和CHOPCHOP。

以CRISPOR为例，操作流程如下：

1. 打开CRISPOR网页，选择基因组版本（如Mouse GRCm39/mm39）；
2. 粘贴候选sgRNA序列（每行一条，5'→3'，不含PAM）；
3. 选择PAM类型为SpCas9的20 bp-NGG；
4. 提交后，工具返回每条sgRNA的MIT特异性评分（0–100，越高越特异）、CFD评分（0–1，越高越特异）以及脱靶位点列表。

MIT评分基于错配位置和数量计算，公式近似为：

MIT = (1 / (1 + Σ 错配权重)) × 100
其中错配权重 = Π (1 - W_i)，W_i为第i位错配的切割概率权重

CFD评分则更细致地考虑每个位置的错配类型。实际筛选中，通常优先保留MIT≥70且CFD≥0.8的候选。若某条sgRNA存在位于外显子区的脱靶位点（尤其错配≤2），应直接排除。

四、设计结果评估与筛选：综合排序

将14条候选的评分汇总，按以下优先级排序：

1. 脱靶安全性：MIT≥70、CFD≥0.8，且无外显子区高相似脱靶；
2. 靶向活性预测：使用Rule Set 2或Doench 2016评分，优先选择活性评分≥0.5的候选；
3. 位置合理性：优先位于外显子5的5'端，且切割位点位于功能域编码区内；
4. 序列质量：GC 45%–65%，无连续4个相同碱基，5'端为G。

假设最终筛选出3条候选：sgPten-1（MIT=82，CFD=0.91，活性=0.68）、sgPten-2（MIT=75，CFD=0.85，活性=0.72）、sgPten-3（MIT=88，CFD=0.94，活性=0.55）。此时建议至少合成2条进行平行实验，因为预测评分与实际切割效率并非完全一致。

五、从设计到实验的衔接

筛选出的sgRNA需要克隆到表达载体中。以px459（SpCas9-2A-Puro）为例，其BbsI位点可用于插入sgRNA。设计时需在20 nt序列两端添加BbsI兼容的粘性末端：

正向寡核苷酸：5'-CACCG + 20nt + -3'
反向寡核苷酸：5'-AAAC + 20nt反向互补 + C-3'

退火后与BbsI酶切的px459连接，转化感受态细胞，挑单克隆测序验证插入序列。验证正确的质粒即可用于细胞转染或电转。

整个流程可概括为：外显子定位 → 序列枚举 → 基础过滤 → 脱靶预测 → 综合评分 → 载体克隆 → 实验验证。每一步的输入输出和判断依据如下图所示。

最后需要强调：任何预测工具都只是辅助。实际编辑效率受细胞类型、转染方式、染色质开放状态等多因素影响。设计完成后，建议先用T7E1酶切或Sanger测序对2–3条候选进行小规模验证，再选择效率最高的sgRNA开展后续实验。