面向非编码基因组的高效扰动

Towards efficient perturbation for the noncoding genome.

作者信息Han Zhang, Shijie Luo, Xiaofeng Wang, Liquan Lin, Ruipu Liang, Chunge Zhong, Yunhan Zhang, Wenchang Zhao, Zhisong Chen, Xiaoya Liu, Yunyan Gao, Feng Chen, Ning Sun, Jialiang Huang, Teng Fei
PMID42717197
发布时间2026-08-11
DOI10.1038/s41467-026-76567-4

实验完整度

高

包含体外细胞系筛选、个体gRNA验证、深度测序、RT-qPCR、Western blot及计算模型构建与验证等多个独立证据层级。

主要模型

A549细胞(人肺癌细胞系) K562细胞(人白血病细胞系) PC-9细胞(人肺癌细胞系) A375细胞(人黑色素瘤细胞系)

重点核对

五种扰动系统在A375细胞中靶向CUL3基因已知调控元件后的基因表达变化(RT-qPCR,n=3技术重复) 在A549和K562细胞中针对69个CRE的gRNA文库筛选(6TG处理 vs DMSO对照,测序深度>500倍覆盖) 在PC-9细胞中针对约200个lncRNA的gRNA文库筛选(细胞适应度读出,维持三周) 双切割片段编辑效率检测(PCR扩增后10% TBE-PAGE,n=3生物学重复) DeepDC模型在独立pgRNA数据集上的预测性能评估(AUC、Accuracy、Precision、Recall、F1分数,5折交叉验证)

摘要

破译非编码基因组(包括重要的顺式调控元件(CRE)和转录非编码RNA基因)的功能仍具有技术挑战性。在此,我们利用大规模并行遗传筛选,系统性地基准测试了五种代表性功能缺失扰动工具的性能,包括单导RNA(gRNA)介导的SpCas9切割或CRISPR干扰,以及配对gRNA(pgRNA)涉及的双SpCas9、Big Papi(配对SpCas9和SaCas9)或双enAsCas12a片段删除方法,在解码非编码基因组作用方面的表现。对于靶向CRE(如增强子),双SpCas9在破坏功能基因组区域方面优于其他方法,效率更高。对于扰动非编码RNA基因,除了双SpCas9外,推荐使用其他RNA靶向方法(如RNA干扰)来区分转录依赖性或非依赖性作用。我们构建了一个深度学习模型DeepDC及其相关网络服务器,以促进最佳双SpCas9 pgRNA设计,从而高效删除基因组片段。总之,我们的工作为选择合适的功能缺失工具以解决非编码基因组的功能复杂性提供了实用指导。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
如何高效选择功能缺失工具来扰动非编码基因组中的顺式调控元件和非编码RNA基因,并提高片段删除效率?
核心机制
双切割片段删除方法(尤其是双SpCas9)通过产生较大的基因组片段缺失,比单切割或CRISPRi更有效地破坏非编码元件的功能;DeepDC模型整合片段长度、GC含量、sgRNA效率、表观遗传信号和Hyena序列评分等特征预测pgRNA效率。
主要证据
在A549和K562细胞中进行的CRE筛选显示双SpCas9识别出最多的功能增强子(AUC=0.70和0.72);在PC-9细胞中的lncRNA筛选同样显示双SpCas9性能最佳(AUC=0.73);DeepDC模型在验证实验中高分pgRNA的编辑效率和基因抑制效果显著高于低分pgRNA。
研究意义
为选择合适的功能缺失工具以解析非编码基因组功能提供了实用指导,并提供了DeepDC网络服务器辅助设计高效的双SpCas9 pgRNA。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

构建五种功能缺失扰动平台

建立并验证靶向非编码基因组的五种CRISPR扰动工具。

构建慢病毒表达质粒,包括sgRNA-SpCas9 KO、sgRNA-CRISPRi、pgRNA-双SpCas9、pgRNA-Big Papi和pgRNA-双enAsCas12a;在A375细胞中靶向CUL3基因的已知调控元件,通过RT-qPCR检测基因表达变化;在A549细胞中通过PCR和PAGE检测片段删除效率。

2

CRE筛选基准测试

比较五种工具在靶向增强子等CRE时的性能。

设计针对69个功能CRE的gRNA文库(sgRNA文库3870个gRNA;pgRNA文库分别为双SpCas9 2480对、Papi 2200对、双enAsCas12a 747对),慢病毒转导A549和K562细胞,在DMSO或2.5 μM 6TG处理下进行细胞适应度筛选,通过高通量测序定量gRNA丰度,计算sigmoid评分,定义功能CRE,并进行ROC曲线分析。

3

lncRNA筛选基准测试

比较五种工具在靶向长链非编码RNA基因时的性能。

设计针对约200个lncRNA的gRNA文库(sgRNA文库4217个gRNA;pgRNA文库分别为双SpCas9 3583对、Papi 2803对、双enAsCas12a 2020对),慢病毒转导PC-9细胞,维持三周进行细胞适应度筛选,通过高通量测序定量gRNA丰度,使用MAGeCK RRA算法计算tRRA评分,定义功能lncRNA,并进行ROC曲线分析。

4

DeepDC模型构建与验证

开发预测双SpCas9 pgRNA片段删除效率的深度学习模型,并验证其指导pgRNA设计的能力。

整合154,457个双SpCas9 pgRNA扰动数据,提取片段长度、GC含量、sgRNA效率指数、表观遗传信号和Hyena评分等特征,使用XGBoost架构训练DeepDC模型;在K562细胞中针对HS1和HS2增强子设计高评分(>0.5)和低评分(<0.2)pgRNA,通过TBE-PAGE检测编辑效率,RT-qPCR检测β-珠蛋白基因表达。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
lentiCRISPRv2Addgene52961
pLenti-Cas9-BlastAddgene52962
pLentiGuide-PuroAddgene52963
pX601Addgene61592
pRDA_174Addgene136476
pRDA_052Addgene136474
dCas9-KRAB-MeCP2载体Addgene110821
pLKO.1-puroAddgene8453
XbaI限制酶Takara1634
BamHI限制酶Takara1605
BsmBI限制酶New England BiolabsR0739L
EcoRI限制酶Takara1611
Esp3I限制酶Thermo Fisher ScientificFD0454
BpiI限制酶Thermo Fisher ScientificFD1014
T4 DNA连接酶Thermo Fisher ScientificEL0016
GeneJET PCR纯化试剂盒Thermo Fisher ScientificK0702
GeneJET凝胶提取试剂盒Thermo Fisher ScientificK0692
无内毒素质粒大量提取试剂盒TIANGEN4992194
无内毒素质粒小量提取试剂盒IITIANGENDP118-02
Dulbecco改良Eagle培养基VivaCell BIOSCIENCESC3114-0500
RPMI-1640培养基VivaCell BIOSCIENCESC3001-0500
胎牛血清ExCellFSP500
青霉素-链霉素VivaCell BIOSCIENCESC3420-0100
Opti-MEM减血清培养基Gibco51200038
NeoFect转染试剂----
嘌呤霉素SolarbioB9300
6-硫代鸟嘌呤----
RIPA裂解液BeyotimeP0013C
蛋白酶抑制剂混合物CWBIOCW2200S
FLAG抗体InvitrogenMA1-91878
GAPDH抗体Santa Cruzsc-25778
HA抗体Proteintech51064-2-AP
HPRT1抗体Proteintech67518-1-lg
PCNA抗体BeyotimeAF1363
PMS2抗体BeyotimeAF1339
MLH1抗体BeyotimeAF1327
MHS2抗体BeyotimeAF7509
MSH6抗体BeyotimeAF0237
山羊抗小鼠IgG H&LAbcamab150116
山羊抗兔IgG H&LAbcamab150080
化学发光底物SparkjadeED0025
成像系统Tanon5200
抗FLAG抗体Proteintech80801-2-RR
抗HA抗体Proteintech66006-2-Ig
Alexa Fluor 488标记二抗(抗FLAG)InvitrogenA-11008
Alexa Fluor 488标记二抗(抗HA)InvitrogenA-11001
流式细胞仪BD LSRFortessa--
UNlQ-10柱式Trizol总RNA提取试剂盒Sangon BiotechB511321
MultiScribe逆转录酶Thermo Fisher4311235
2X UltraSYBR混合物SparkjadeAH0105
QuantStudio 5实时PCR系统Thermo FisherA28319
Ex Taq DNA聚合酶TakaraRR001B
Phanta Max超保真DNA聚合酶VazymeP505
T7核酸内切酶ISangon BiotechB610169
T4多核苷酸激酶thermoEK0032

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
CRE筛选
细胞系(A549、K562)、gRNA文库覆盖度(>500倍)、6TG浓度(2.5 μM)、筛选持续时间(三周)、测序模式(PE150)
阅读提示:Methods中“CRE or enhancer screens”小节及Figure 2 legend
lncRNA筛选
细胞系(PC-9)、gRNA文库覆盖度(>500倍)、筛选持续时间(三周)、测序模式(PE150)
阅读提示:Methods中“LncRNA screens”小节及Figure 3 legend
双切割片段编辑效率检测
细胞系、pgRNA转导MOI(~0.3)、感染后时间(7天)、PCR扩增体系、PAGE胶浓度(10% TBE-PAGE)、生物学重复数(n=3)
阅读提示:Methods中“Dual-cut fragment editing assay”小节及Figure 1d legend
DeepDC模型验证
细胞系(K562)、高/低评分pgRNA数量(n=11和n=7)、RT-qPCR检测的基因(HBE1、HBG1/2、HBB)、技术重复数(n=3)
阅读提示:Methods中“Applying the DeepDC model for enhancer characterization”相关段落及Figure 5 legend
T7EI assay
sgRNA数量(n=72)、PCR扩增体系、T7 Endonuclease I用量(1 μL)、消化时间(37°C 30分钟)、凝胶浓度(8% polyacrylamide TBE gels)
阅读提示:Methods中“T7 Endonuclease I assay”小节及Supplementary Fig. 1i, j