IF31.7｜碱基编辑筛选揭示癌症耐药机制遗传图谱

大规模碱基编辑筛选系统绘制十种抗癌药耐药变异图谱，并揭示四类功能变异。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

肿瘤靶向治疗长期疗效常受限于耐药，而耐药多由基因组单核苷酸变异（SNV）驱动，如药物靶点或信号通路蛋白的点突变。传统方法依赖患者复发后肿瘤活检测序，样本获取困难、周期长，且只能回顾性推断高频变异，难以直接比较不同变异效应。CRISPR碱基编辑技术可在内源基因高效引入C>T或A>G点突变，为前瞻性功能筛选提供可能。已有研究利用碱基编辑解析细胞因子或抑制剂耐药机制，但缺乏对多种药物、多个基因的系统性变异-功能图谱。

科学问题

卡在哪里： 现有耐药研究多依赖临床样本回顾性分析，难以快速、前瞻性地鉴定并比较大量未知意义变异（VUS）的功能，且缺乏统一的变异功能分类框架，限制了对耐药机制的全面理解和临床转化。

本文要回答： 利用碱基编辑筛选在四种癌细胞系中前瞻性鉴定十种肿瘤药物耐药相关变异，建立功能分类体系，并结合单细胞转录组揭示其作用机制。

技术路线

作者主要用了：CRISPR碱基编辑筛选、单细胞转录组（perturb-seq）、竞争性增殖实验、免疫荧光、染色质分级。

作者构建了靶向11个癌症基因的gRNA文库（22,816条），在四种携带不同驱动突变的癌细胞系中表达诱导型CBE或ABE碱基编辑器，进行46个独立混合筛选。通过比较药物处理组与对照组gRNA丰度变化，鉴定影响药物敏感性的变异。随后利用单细胞转录组（perturb-seq）解析耐药变异的转录程序，并通过竞争性增殖实验、免疫荧光和染色质分级等方法验证关键变异的功能。

图 1　Fig. 1。图1展示碱基编辑筛选整体设计及质控。a为筛选流程示意图；b显示CBE和ABE筛选中靶向必需基因剪接位点的gRNA显著耗竭，证明编辑活性；c比较不同细胞系对驱动癌基因编辑的敏感性，PC9对EGFR、HT-29对BRAF和MYC编辑更敏感；d展示HT-29中BRAF错义变异的功能域分布，关键结构域（如RAS结合域、激酶域）的gRNA显著耗竭，并鉴定出已知耐药变异L505F。

核心亮点

亮点 1｜四类功能变异分类框架

整合46个筛选数据，将影响药物敏感性的变异分为四类：药物成瘾变异（仅在药物存在时增殖优势，无药时有害）、经典耐药变异（仅在有药时优势）、驱动变异（有无药均有优势）和药物增敏变异（仅在有药时有害）。以曲美替尼筛选为例，四类变异均被检出，且对照gRNA无一命中，信噪比高。该分类为理解耐药机制提供了系统框架。

图 2　Fig. 2。图2展示曲美替尼和达拉非尼+西妥昔单抗筛选结果及验证。a、b为gRNA z-score散点图，显示药物处理组与质粒库比较，鉴定出四类变异；c为EGFR与西妥昔单抗、MEK1与曲美替尼的共晶结构，标注耐药变异位置；d为增殖实验验证，所有测试变异均导致耐药，药物成瘾变异在无药时生长减慢；e、f显示药物成瘾变异细胞基础MAPK信号和p21升高，β-半乳糖苷酶染色增强，药物处理可逆转。

亮点 2｜药物成瘾变异与间歇给药策略

在BRAF/EGFR抑制剂联合筛选中，鉴定出MAPK通路激活型药物成瘾变异（如KRAS Q61R/E62G、MEK2 Y134H）。这些变异在无药时导致细胞增殖减缓、p21升高、β-半乳糖苷酶染色增强，提示癌基因诱导衰老；而药物处理可逆转这些表型。患者样本中KRAS与BRAF激活突变互斥，支持药物成瘾依赖于细胞原有驱动突变。该发现提示携带此类变异的患者可能受益于药物假期。

图 3　Fig. 3。图3展示驱动变异和耐药变异。a比较AKT1和PIK3CA氨基酸位置的COSMIC突变计数与gRNA z-score，显示已知热点（AKT1 E17K、PIK3CA E545K/E542K）和罕见驱动变异；b为pictilisib筛选，显示AKT1驱动变异导致耐药，而PIK3CA激活变异在无药时有优势但药物仍敏感；c为sotorasib和adagrasib筛选，显示KRAS药物结合口袋附近变异（R68G、D69G）导致交叉耐药，下游MEK1/2激活变异也导致耐药。

亮点 3｜PARP1变异对抑制剂的不同响应

在尤文肉瘤细胞中筛选PARP1/2变异对奥拉帕利和尼拉帕利的敏感性。发现PARP1螺旋亚结构域（HD）的错义变异（如I691T）增敏两种抑制剂，而ART亚结构域的Y889C对尼拉帕利耐药却对奥拉帕利增敏。机制上，Y889C在奥拉帕利存在下增加PARP捕获，而在尼拉帕利存在下减少捕获。该结果提示奥拉帕利可能有效治疗携带Y889C的尼拉帕利耐药肿瘤。

图 4　Fig. 4。图4展示PARP1/2筛选结果及验证。a为PARP1和PARP2的gRNA z-score沿氨基酸位置分布，标注耐药和增敏变异，显示增敏变异集中于催化域HD亚结构域；b为竞争实验验证PARP1变异对奥拉帕利和尼拉帕利的敏感性差异；c为PARP1 Y889C的免疫荧光和染色质分级实验，显示奥拉帕利增加PARP捕获而尼拉帕利减少；d为PARP1催化域结构，标注Y889位置。

亮点 4｜EGFR截短变异增敏EGFR抑制剂

在PC9细胞中，靶向EGFR C端调控区的gRNA引入剪接变异，导致外显子27后内含子保留、移码和提前终止密码子，产生截短蛋白。该截短蛋白表面表达增加，MAPK信号对EGFR抑制剂更敏感，且对多种EGFR抑制剂（吉非替尼、奥希替尼、厄洛替尼、拉帕替尼、西妥昔单抗）增敏，但对化疗药物无影响。临床中携带EGFR C端截短的非小细胞肺癌患者对EGFR抑制剂有响应，提示该变异可作为患者分层标志物。

图 5　Fig. 5。图5展示EGFR抑制剂筛选结果。a、b为PC9细胞中CBE和ABE筛选对吉非替尼和奥希替尼的gRNA z-score，显示T790M对吉非替尼强耐药，奥希替尼耐药变异位于D1006/D1012附近，C797S因碱基编辑限制未检出；c为prime editing筛选验证C797S等变异对奥希替尼的耐药性；d为竞争实验显示D1012N、A1013D/G变异在奥希替尼存在下生长优势增强。

亮点 5｜耐药变异的免疫逃逸特征

通过perturb-seq分析耐药变异的转录组，发现药物成瘾变异和部分驱动变异下调JAK-STAT通路和抗原呈递基因（如B2M、HLA-A）。在HT-29细胞和患者来源类器官中验证了B2M和HLA蛋白表达下调。共培养实验显示，MEK抑制剂预处理可增强T细胞杀伤，但携带MAP2K1 S194P药物成瘾变异的类器官除外。将perturb-seq差异基因与临床试验患者转录组关联，发现耐药变异（尤其药物成瘾变异）与更差的无进展生存期相关，提示耐药变异可能同时驱动免疫逃逸。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： CRISPR碱基编辑筛选、单细胞转录组（perturb-seq）、竞争性增殖实验、免疫荧光、染色质分级
• 公开数据： ERP146490、ERP148732、ERP144241、ERP141719、ERP156437
• 代码： 公开仓库 https://github.com/CansuDincer/BEstimate、https://github.com/MarioniLab/BE_perturb_seq_drug_resistance、https://github.com/MatthewACoelho/Res1_analysis；本地已克隆：BEstimate、Res1_analysis、BE_perturb_seq_drug_resistance、PRIDICT
• 谁该点开原文： 从事肿瘤基因组学、功能基因组学或药物耐药机制研究的科研人员及临床转化研究者。

带走一句

碱基编辑筛选能前瞻性、系统性地绘制耐药变异图谱，结合单细胞转录组可解析变异功能类别和机制，为临床患者分层、联合用药和药物假期策略提供依据。