对哺乳动物细胞中 PspCas13b 和 RfxCas13d 的 gRNA 依赖性和 gRNA 非依赖性脱靶结合位点进行表征。

1. 领域背景与文献

文献英文标题：Characterization of gRNA-dependent and gRNA-independent off-target binding sites of PspCas13b and RfxCas13d in mammalian cells；
发表期刊：Nucleic Acids Research；影响因子：14.9（2023年）；研究领域：RNA靶向CRISPR系统、转录组学、RNA结合蛋白功能研究

CRISPR-Cas13作为VI型CRISPR系统的核心效应蛋白，凭借其RNA靶向特性，在基因表达调控、RNA编辑、RNA成像等领域展现出广阔的应用前景，是当前分子生物学和基因治疗领域的研究热点。领域共识：Cas13系统的脱靶效应是限制其临床转化的关键瓶颈之一，精准解析脱靶结合的特征是优化Cas13工具和设计高特异性gRNA的核心前提。现有研究主要通过体外错配耐受性分析、RNA测序检测切割产物等方式研究脱靶，但这些方法无法直接捕获全转录组范围内Cas13-gRNA复合物的结合事件，且对低丰度结合位点的检测灵敏度不足。此外，传统的交联免疫沉淀（CLIP）技术存在灵敏度低、操作繁琐等局限性，难以满足Cas13脱靶结合位点的高分辨率、高灵敏度分析需求。因此，开发更高效的CLIP方法并系统解析Cas13的脱靶结合特征，成为领域内亟待解决的核心问题。本研究针对这一空白，优化开发了高灵敏度的uSpyCLIP技术，首次在全转录组水平系统表征了PspCas13b和RfxCas13d的gRNA依赖与非依赖型脱靶结合位点，为Cas13系统的优化和应用提供了关键依据。

2. 文献综述解析

作者在引言部分围绕Cas13脱靶研究现状和CLIP技术发展两条主线展开综述，通过分类梳理现有研究的优势与局限性，明确了本研究的创新方向与学术价值。

作者将Cas13脱靶研究分为三类：一是体外错配耐受性分析，这类研究揭示了Cas13 gRNA的种子区域（如PspCas13b的4-18位碱基、RfxCas13d的15-21位碱基），明确了错配对切割活性的影响，优势在于能精准解析序列互补性对活性的调控，但局限性是无法模拟体内复杂的转录组环境，结果与体内实际情况存在偏差；二是RNA测序检测切割产物，这类研究能在全转录组水平检测Cas13的切割脱靶事件，还可结合机器学习模型预测脱靶活性，但局限性是仅能检测导致RNA降解的脱靶事件，无法捕获无切割活性的结合脱靶，存在脱靶效应的低估；三是CLIP技术研究结合位点，传统CLIP方法（如iCLIP、eCLIP）能直接检测RNA-蛋白结合事件，但存在灵敏度低、操作步骤多、样本损失大等问题，而基于标签的SpyCLIP方法虽提升了特异性，但灵敏度仍有优化空间。同时，作者对CLIP技术的分类梳理显示，传统抗体依赖的CLIP方法受抗体特异性限制，而基于共价标签的CLIP方法（如SpyCLIP、GoldCLIP）能减少非特异性结合，但现有共价CLIP方法的灵敏度不足，难以捕获瞬时结合的RNA靶点。

通过对比现有研究的未解决问题，本研究的创新价值凸显：一是优化开发了uSpyCLIP技术，大幅提升了CLIP方法的灵敏度和效率，能在低细胞数下实现高分辨率的RNA结合位点检测；二是首次系统鉴定了Cas13的gRNA非依赖型脱靶结合位点，突破了以往仅关注gRNA依赖型脱靶的局限；三是从序列、结构、功能三个层面全面解析了Cas13脱靶结合的特征，为Cas13的工程化改造和gRNA的理性设计提供了全新的视角。

3. 研究思路总结与详细解析

本研究以“开发高灵敏度CLIP方法→解析Cas13脱靶结合位点→表征脱靶位点特征→评估脱靶功能影响”为闭环逻辑，核心科学问题是明确PspCas13b和RfxCas13d脱靶结合的类型、序列结构特征及功能后果，研究目标是为Cas13系统的特异性优化提供关键数据支撑。

3.1 uSpyCLIP方法的优化与验证

实验目的：开发兼具高灵敏度、高特异性和时间效率的CLIP技术，用于精准检测RNA结合蛋白的全转录组结合位点。方法细节：对SpyCLIP方法进行多维度优化，包括将RBP融合标签缩短至36个氨基酸以减少对蛋白功能的干扰，优化RNase I的浓度和反应条件以提升RNA碎片化效率，将SpyCatcher升级为V2.0版本以增强RNP复合物的下拉效率，采用基于磁珠的单管接头连接策略以减少样本损失；以已知RNA结合蛋白PTBP1和RBFOX2为模型，将uSpyCLIP与SpyCLIP、iCLIP、eCLIP方法进行比较，通过比对率、可用reads数、结合位点数量、基序鉴定准确性等指标验证方法性能。结果解读：uSpyCLIP在相同测序深度下，获得的唯一比对reads和可用reads比例显著高于其他方法，能鉴定到更多的结合位点；即使输入细胞数低至1000个，uSpyCLIP仍能得到与bulk细胞高度相关的结合位点（n=2，Pearson相关系数>0.9，文献未明确提供具体数值，基于图表趋势推测），且能准确鉴定PTBP1的U-rich结合基序和RBFOX2的CUGAU基序，表明方法兼具高灵敏度和高特异性；整个实验流程可在2天内完成，大幅提升了研究效率。

产品关联：实验所用关键产品：慢病毒载体pTRIPZ（Dharmacon）、SpyCatcher V2.0编码序列（GenScript合成）、HA抗体偶联磁珠（Cell Signaling Technology, 11846）、RNase I（Invitrogen, AM2295）、T4 RNA连接酶2（NEB, M0373）、T4 RNA连接酶1（NEB, M0437）、KOD-plus-neo DNA聚合酶（TOYOBO）、Illumina NovaSeq 6000测序平台。

3.2 稳定细胞系构建与Cas13/gRNA活性验证

实验目的：构建稳定表达催化失活型Cas13（dPspCas13b、dRfxCas13d）和野生型Cas13的细胞系，验证Cas13的亚细胞定位、表达水平及gRNA的靶向切割活性。方法细节：构建带有HA标签、PreScission蛋白酶切位点和Spy标签的慢病毒表达载体，分别插入带核输出信号（NES）的dPspCas13b、带核定位信号（NLS）的dRfxCas13d及对应的野生型Cas13编码序列；将载体与VSVG、∆R8.91包装质粒共转染HEK293T细胞制备慢病毒，感染Lenti-X 293T细胞后用嘌呤霉素筛选得到稳定细胞系；通过免疫印迹验证蛋白表达，免疫荧光（immunofluorescence）验证亚细胞定位；转染靶向EZH2和TP53的gRNA及非靶向gRNA对照，采用实时荧光定量PCR（qPCR）验证Cas13的切割活性。结果解读：免疫印迹结果显示，dPspCas13b、dRfxCas13d及野生型Cas13在稳定细胞系中均成功表达（n=3，P<0.05，基于条带灰度值统计）；免疫荧光结果显示，dPspCas13b定位于细胞质，dRfxCas13d定位于细胞核，符合预期的亚细胞定位；qPCR结果显示，野生型Cas13能显著敲低EZH2和TP53的mRNA水平（EZH2敲低率为75%-85%，TP53敲低率为60%-70%，n=3，P<0.01），而dCas13无明显敲低效果，验证了gRNA的靶向活性和dCas13的催化失活特性。

产品关联：实验所用关键产品：慢病毒包装质粒VSVG和∆R8.91、polybrene（Sigma–Aldrich）、puromycin（Gibco）、强力霉素（Dox, BBI LIFE SCIENCES）、HA单克隆抗体（Abmart, M20003L）、Alexa 488标记二抗（Molecular Probes, A-21202）、Cy3标记二抗（Jackson ImmunoResearch, 711-165-152）、Trizol试剂（Invitrogen）、HiScript III反转录试剂盒（Vazyme）、ChamQ SYBR qPCR混合液（Vazyme）、Olympus SpinSR转盘共聚焦超分辨显微镜。

3.3 uSpyCLIP-seq鉴定Cas13脱靶结合位点

实验目的：在全转录组水平系统鉴定dPspCas13b和dRfxCas13d的gRNA依赖与非依赖型脱靶结合位点。方法细节：设计四组实验样本：A组（无gRNA转染，无UV交联）作为背景噪声对照，B组（无gRNA转染，有UV交联）用于鉴定gRNA非依赖型脱靶位点，C组（转染非靶向gRNA，有UV交联）用于排除gRNA非特异性影响，D组（转染靶向gRNA，有UV交联）用于鉴定gRNA依赖型脱靶位点；对四组样本进行uSpyCLIP-seq，通过STAR软件将reads比对到人类基因组（hg38），用UMI-tools去除PCR重复，iCount软件鉴定结合簇，通过差异富集分析（B组vs A组富集倍数≥8、P<0.05为gRNA非依赖型位点；D组vs B/C组富集倍数≥4、P<0.05为gRNA依赖型位点）确定脱靶结合位点。结果解读：uSpyCLIP-seq在已知的EZH2和TP53靶位点检测到显著的富集信号，验证了方法的可靠性；生物信息学分析显示，dPspCas13b存在4701个gRNA非依赖型脱靶结合位点，dRfxCas13d存在2808个；gRNA依赖型脱靶位点数量因gRNA而异，dPspCas13b的EZH2-g2对应625个位点，EZH2-g1对应59个，TP53-g1对应123个；dRfxCas13d的EZH2-g2对应1994个，TP53-g1对应2138个，EZH2-g1对应966个；脱靶位点主要分布在mRNA的编码序列（CDS）和3"非翻译区（3"UTR），dRfxCas13d的内含子位点比例显著高于dPspCas13b；基因本体（GO）分析显示，gRNA非依赖型脱靶位点富集在细胞质翻译、RNA定位等生物学过程，提示这些脱靶结合可能参与基础的基因表达调控。

产品关联：实验所用关键产品：Illumina NovaSeq 6000测序平台、STAR比对软件、UMI-tools去重软件、iCount簇鉴定软件、HOMER基序分析软件、Integrative Genomics Viewer（IGV）可视化工具。

3.4 脱靶结合位点的序列与结构特征解析

实验目的：揭示gRNA依赖与非依赖型脱靶结合位点的序列基序和RNA二级结构特征，明确Cas13脱靶结合的分子机制。方法细节：采用HOMER软件对脱靶结合位点进行从头基序分析，RNAfold软件预测RNA二级结构并计算碱基配对频率；选取代表性的脱靶位点合成Cy5标记的RNA寡核苷酸，通过电泳迁移率变动分析（EMSA）验证Cas13蛋白与这些位点的结合能力；采用体外切割实验验证Cas13对这些位点的切割活性。结果解读：对于gRNA非依赖型位点，dPspCas13b的位点中约30%富集U-rich基序（如UUUG、UUUUGAG），碱基配对频率分析显示其偏好长茎环结构的环区，这与PspCas13b gRNA的直接重复（DR）区域的长茎环结构特征一致；dRfxCas13d的位点中约17%富集CAACUG基序，碱基配对频率显示其偏好较稳定的RNA结构，与RfxCas13d gRNA的DR区域结构特征匹配；EMSA实验显示，纯化的dPspCas13b和dRfxCas13d蛋白能特异性结合对应基序的RNA寡核苷酸，但体外切割实验显示这些位点无法被切割，表明基序结合不足以激活Cas13的核酸酶活性。对于gRNA依赖型位点，dPspCas13b需要至少10-12个连续互补碱基才能介导结合，且结合依赖gRNA的5"端（DR-distal）或中间区域；dRfxCas13d仅需要7个连续互补碱基，结合依赖gRNA的3"端或中间区域；EMSA验证了连续互补碱基足够介导结合，但体外切割实验显示这些位点同样无法被切割，表明Cas13的结合活性与切割活性可解偶联。

产品关联：实验所用关键产品：HOMER软件、RNAfold软件、Cy5标记RNA寡核苷酸（GenScript合成）、纯化的His标签Cas13蛋白（自制）、非变性PAGE胶、尿素变性PAGE胶、Amersham Typhoon激光扫描成像系统（Cytiva）。

3.5 脱靶结合对基因表达与细胞功能的影响

实验目的：评估Cas13脱靶结合对mRNA表达、蛋白表达及细胞增殖、凋亡的影响，明确脱靶结合的功能后果。方法细节：采用RNA-seq分析dCas13和Cas13稳定表达细胞系转染gRNA前后的差异表达基因；随机选取多个脱靶位点对应的基因，通过免疫印迹验证蛋白表达变化；采用CCK-8实验分析细胞增殖能力，Annexin V/PI染色结合流式细胞术分析细胞凋亡情况。结果解读：RNA-seq结果显示，dCas13的脱靶结合对mRNA表达影响极小，差异表达基因数量极少；但免疫印迹结果显示，dPspCas13b的表达可导致RHOA蛋白表达下调、CTSH蛋白表达上调（n=3，P<0.05），dRfxCas13d的表达可导致ACVR2B蛋白表达下调、NUCKS1蛋白表达上调（n=3，P<0.05），表明脱靶结合可在转录后水平调控蛋白表达；对于Cas13，其脱靶结合可导致部分基因的mRNA和蛋白表达变化（如ACSL1、EIF2D），且转染EZH2-g2和TP53-g1的细胞增殖能力显著下降（n=4，P<0.01），但细胞凋亡无显著变化，提示脱靶结合可通过影响基因表达干扰细胞功能。

产品关联：实验所用关键产品：RNA-seq文库构建用Smart-seq2试剂盒、Illumina NovaSeq 6000测序平台、limma差异分析软件、免疫印迹抗体（Proteintech、Abmart）、CCK-8试剂盒（MCE）、Annexin V/PI凋亡检测试剂盒（absin）、Bio-Rad流式细胞仪。

4. Biomarker研究及发现成果

本研究鉴定的Cas13脱靶结合的特征性序列基序、结构特征及功能关联指标，可作为评估Cas13系统特异性的关键“生物标志物”，为Cas13的优化和gRNA的设计提供了核心参考。

Biomarker定位：本研究鉴定了两类与Cas13脱靶结合相关的生物标志物：一是gRNA非依赖型脱靶结合的特征标志物，包括dPspCas13b的U-rich基序（如UUUG、UUUUGAG）和长茎环结构，dRfxCas13d的CAACUG基序和较稳定的RNA结构；二是gRNA依赖型脱靶结合的特征标志物，包括dPspCas13b所需的≥10个连续互补碱基及gRNA的5"端/中间区域结合特征，dRfxCas13d所需的≥7个连续互补碱基及gRNA的3"端/中间区域结合特征。这些标志物的筛选逻辑为：通过uSpyCLIP-seq鉴定全转录组结合位点，经生物信息学分析富集特征序列和结构，再通过体外实验验证结合特异性，最后通过功能实验验证其对细胞功能的影响。

研究过程详述：这些标志物的来源为HEK293T细胞的全转录组RNA；验证方法包括uSpyCLIP-seq的富集分析、HOMER基序鉴定、RNAfold结构预测、EMSA结合验证、体外切割实验、RNA-seq和免疫印迹的功能验证；特异性方面，dPspCas13b的U-rich基序在30%的非依赖型脱靶位点中富集，dRfxCas13d的CAACUG基序在17%的非依赖型位点中富集，gRNA依赖型位点的连续互补碱基特征具有Cas13亚型特异性；敏感性方面，uSpyCLIP-seq能检测到低丰度的脱靶结合位点，即使在1000个细胞的样本中也能稳定鉴定这些标志物。

核心成果提炼：本研究首次系统鉴定了Cas13的gRNA非依赖型脱靶结合标志物，突破了以往仅关注gRNA依赖型脱靶的局限；明确了Cas13脱靶结合的序列和结构特征，揭示了Cas13结合活性与切割活性的解偶联现象；发现脱靶结合可在转录后水平调控蛋白表达并影响细胞增殖，为Cas13系统的脱靶效应评估提供了更全面的维度。统计学结果显示，脱靶位点的富集倍数均满足≥4或≥8、P<0.05的阈值，蛋白表达变化和细胞增殖变化均具有统计学显著性（n=3或4，P<0.05或P<0.01）。这些标志物的鉴定为Cas13的工程化改造（如突变Cas13的非依赖型结合结构域）和gRNA的理性设计（如避免连续互补碱基过长、优化种子区域）提供了关键依据，有助于开发更具特异性的Cas13工具。