CRISPR Array 超多重 PCR 引物 + crRNA 设计器

📚 超多重引物＋crRNA 一体化设计方法学

这套工具背后的完整工程方法，源自一条真正跑通 371 重级呼吸道病原面板的设计管线沉淀。方法学要点如下：

①「四均衡」是多重成败的判定标准

效率均衡（各靶标 E 落在 90–110%，多重定量下限由最弱靶标决定）、Tm 兼容（全部引物对 Tm 收敛在 ±2 ℃）、光谱无串扰（NTC 必须验证无跨通道泄漏）、产物长度可分辨（本项目 50–130 bp 走 CE 分辨）。参数上标准 PCR 要向「多重」收紧：Tm 60–65 ℃、GC 45–55%、引物 20–25 nt、超多重引物浓度降到 0.001–0.1 μM。

②交叉二聚体是 N² 问题，全对 ΔG 筛查是工具的核心价值

N 对引物 = 2N 条寡核苷酸，二聚体组合 C(2N,2)——21 重是 861 对、96 重是 18336 对，人工枚举不可行。经验判据：只盯 3′ 端（可延伸二聚体才是灾难）；交叉二聚体 ΔG 优于 −3~−5 kcal/mol 才需警惕；二级结构自由能弱于 −9 kcal/mol 即为风险。371 重实测：742 条引物、27.5 万对组合里 ΔG ≤ −8 占 0.73%，属高风险但非不可行——靠低引物浓度＋热启动酶＋较高退火温度＋替换二聚体 hub 引物来压。

③SADDLE 模拟退火：把二聚体风险显式建模成 Badness

超多重引物组优化的当前最优解（SADDLE，Nat Commun 2022）：候选池按引物-模板结合能 ΔG° ≈ −11.5 kcal/mol 截取；损失函数 = 全体两两 Badness 之和，Badness 只在两引物存在 ≥4 连续互补碱基时计算（连续配对越长、离各自 3′ 端越近、配对区 GC 越多分越高）；再用模拟退火随机换引物逃局部最优。实验把二聚体 reads 从 90.7% 压到 4.9%（AUROC 0.958），但纯算法不能完全消除强二聚体——湿实验确认不可省。

④crRNA 层：PAM 门控种子排他，防止「偶然 8-mer 误判」

跨靶标交叉校验必须做 PAM 门控——只在 PAM 精确匹配窗口内比对 spacer，否则会把无 PAM 的偶然 8-mer 误判为风险。约定：LbCas12a spacer 23 nt、GC 40–65%、与 PCR 引物反向互补 LCS < 10 nt；同管多重须查 PAM 近端 8 nt 种子跨靶标错配 ≥2；保守性以「核心 = PAM＋8 nt」为灵敏度下限。

⑤371 重实测暴露的真问题：鉴别型重设计

371 重面板 in-silico 结果：有可用引物对 371/371、有合规 crRNA 352/371（T1 经典 189 · T2 三核苷酸 132 · T3 非经典 31）。交叉扩增是 crRNA 失败的根本原因：属内近缘簇中，靶标引物在近缘种基因组上会产出与自身扩增子近乎同源（种子错配 0）的产物，任何 crRNA 都无法区分——交换备选扩增子无效，必须 鉴别型重设计（要求 F/R 引物在每个近缘参考双链上无 ≤2 错配匹配，从源头消除交叉产物）。

⑥一条管线沉淀的实战教训（可直接照搬的 checklist）

• 覆盖率/特异性结论必须实测复核，不得引用文献数字——实测发现 5 处文献或坐标错误；
• e-PCR 必须双链（只搜正链会严重低估覆盖率）＋ 3′ 端错配加权，覆盖率报「≤2 / ≤3 错配」两档；
• 分组是硬约束问题：同管同通道长度差 ≥5 bp，理论上界由冲突图最大团给出，不要用启发式搜索失败当作「不可行」；
• primer3 默认 Tm 是 50 mM Na⁺ 条件，报湿实验要附 @PCR 列；
• 扩增子内部二级结构筛查是当前唯一盲点（发夹阈值 −9 kcal/mol）；
• 交付物要加引物浓度建议列与验证 SOP（单重→两两→全重、NTC 串扰/二聚体、E 90–110%）。