1. 领域背景与文献
文献英文标题：Clonal dynamics shaped by diverse drug-tolerant persister states in melanoma resistance；
发表期刊：未公开；影响因子：未公开；研究领域：肿瘤学-黑色素瘤耐药机制
黑色素瘤是恶性程度最高的皮肤肿瘤之一，BRAF V600突变是其核心驱动突变，BRAF/MEK抑制剂联合靶向治疗已成为晚期患者的标准治疗方案，但超过70%的患者会在1-2年内出现耐药，其中药物耐受持留细胞（Drug-Tolerant Persister, DTP）是介导非遗传耐药、驱动肿瘤复发的关键细胞亚群。领域发展关键节点包括：2002年BRAF突变靶点的发现奠定了靶向治疗基础；2011年首个BRAF抑制剂维莫非尼获批上市；2014年BRAF/MEK抑制剂联合治疗方案获批，进一步提升疗效；2014年左右DTP概念被提出，开启了非遗传耐药机制的研究热潮；近年谱系追踪技术的发展为解析肿瘤克隆异质性提供了核心工具。当前研究热点聚焦于DTP的克隆起源、转录可塑性、空间异质性及非遗传耐药调控网络，但未解决的核心问题包括：DTP克隆动态的时间演化规律、不同持留状态的转录调控机制、空间组织模式及细胞间通讯的功能关联。
现有研究多局限于单一维度的分析，如仅通过单细胞转录组鉴定DTP的转录状态，或仅通过谱系追踪解析克隆动态，缺乏整合克隆追踪、单细胞转录组、空间成像的技术平台来系统解析DTP的克隆起源与转录可塑性的动态关联，以及不同持留状态的空间分布与互作机制。针对这一研究空白，本研究开发了MeRLin（Melanoma Resistance Lineage tracing）谱系追踪平台，整合细胞条形码标记、单细胞转录组分析及RNA荧光原位杂交（RNA-FISH）空间验证，旨在解析黑色素瘤BRAF/MEK抑制剂耐药过程中DTP的克隆动态、转录状态演化及空间组织特征，为耐药干预提供新的靶点与机制依据。
2. 文献综述解析
作者从谱系追踪技术发展、体外/体内模型优劣、DTP耐药机制研究现状三个维度对领域内现有研究进行分类评述。
谱系追踪技术方面，Confetti、LeGO等系统实现了克隆可视化追踪，但难以结合单细胞转录组解析克隆的转录状态；LARRY、CellTag等静态条形码技术可与单细胞转录组结合，但缺乏便捷的空间验证手段；FateMap、intMEMOIR等技术整合了空间信息，但操作复杂度高、通量有限。模型研究方面，2D细胞培养适合高通量药物筛选，但缺乏体内肿瘤的异质性与微环境复杂性；患者来源异种移植（Patient-Derived Xenograft, PDX）模型更贴近临床肿瘤的基因组特征与空间结构，但谱系追踪技术的应用仍有限。DTP耐药机制方面，已知DTP可通过转录重编程、代谢重塑、表观遗传调控等非遗传机制逃避治疗压力，但DTP的克隆起源时间动态、不同持留状态的调控网络及空间互作机制尚不清楚，且多数研究未整合克隆、转录、空间多维度数据进行系统分析。
本研究的创新价值在于：首次开发MeRLin整合平台，实现了高复杂度克隆标记、单细胞转录组分析与便捷RNA-FISH空间验证的多维度联合分析；首次在PDX模型中动态解析了DTP克隆起源于治疗前的微小亚群，且通过适应性转录重编程而非预存耐药克隆主导肿瘤复发；首次鉴定了四种DTP持留状态（应激样、脂质代谢、PI3K信号、细胞外基质（ECM）重塑），并系统解析了其转录调控网络、空间分布模式及细胞间通讯机制，填补了领域内多维度整合解析DTP耐药机制的空白。
3. 研究思路总结与详细解析
本研究的核心目标是开发MeRLin谱系追踪平台，解析黑色素瘤BRAF/MEK抑制剂耐药过程中DTP的克隆动态、转录可塑性及空间组织特征；核心科学问题包括DTP的克隆起源与演化规律、不同持留状态的转录调控机制、空间分布与细胞间通讯的功能关联；技术路线遵循“平台构建→体外/体内模型验证→多组学数据采集→整合分析机制→空间验证”的闭环逻辑。
3.1 MeRLin谱系追踪平台构建
实验目的：开发可同时实现克隆标记、单细胞转录组分析、空间原位验证的整合式谱系追踪工具。
方法细节：以pCDH-EF1a-eFFly-mCherry为载体骨架，通过限制性酶切替换mCherry为远红荧光蛋白mNeptune2.5，优化载体KpnI酶切位点以适配条形码插入；合成包含半随机重复序列的条形码寡核苷酸，通过Gibson组装构建复杂度达2.89×10^6的条形码文库；进行慢病毒包装后，以感染复数（MOI）<0.4转染WM4237-1黑色素瘤细胞，流式分选mNeptune2.5阳性细胞以确保单个细胞携带唯一条形码。
结果解读：成功构建的MeRLin载体包含荧光素酶、mNeptune2.5及半随机条形码序列，可实现体内生物发光成像、单细胞转录组条形码读取及RNA-FISH空间克隆验证；转染后细胞的条形码多样性满足克隆追踪需求，为后续实验奠定基础。
产品关联：实验所用关键产品：Addgene的pCDH-EF1a-eFFly-mCherry（货号104833）、pcDNA3-mNeptune2.5（货号51310），NEB的Q5定点突变试剂盒（E0552）、NEBuilder HiFi DNA组装试剂盒（E2621），Lucigen的Endura电转化感受态细胞（60242-1），Bio-Rad的Gene Pulser II电转仪，Origene的慢病毒浓缩试剂（TR30025）。
3.2 体外BRAF/MEK抑制剂处理的克隆与转录组分析
实验目的：解析体外培养条件下黑色素瘤耐药克隆的起源与转录特征。
方法细节：将MeRLin标记的WM4237-1细胞分为对照组与BRAF/MEK抑制剂处理组（0.3μM dabrafenib+30nM trametinib），处理1周后收集细胞进行单细胞转录组测序；从测序数据中提取条形码序列，分析克隆多样性、细胞周期分布及差异基因表达。
结果解读：处理组细胞出现G1期阻滞（P<2.2×10^-16，双尾Fisher精确检验），黑色素分化标志物EDNRB、DCT、TYRP1表达上调；条形码多样性显著降低（Shannon指数，P<2.2×10^-16，双尾Hutcheson t检验），耐药克隆起源于对照组中的微小亚群；耐药克隆相对富集氧化磷酸化（OXPHOS）通路活性，下调TGFβ通路活性，且高表达BCL2L1、PIK3CB、STAT2（P<0.05，双尾Wilcoxon秩和检验）。
产品关联：实验所用关键产品：NCI的dabrafenib（NSC764134）、trametinib（NSC758246），10x Genomics的Chromium Next GEM单细胞3"试剂盒v3.1，Illumina的NovaSeq 6000测序平台。
3.3 PDX模型中耐药克隆的动态追踪
实验目的：解析体内黑色素瘤耐药过程中克隆的时间动态与转录状态演化。
方法细节：将MeRLin标记的WM4237-1细胞接种于免疫缺陷NSG小鼠右侧 flank，待肿瘤体积达500mm³时给予含BRAF/MEK抑制剂的饲料；分别在预处理（Day0）、早期微小残留病（Early Minimal Residual Disease, MRD，Day21）、晚期MRD（Day57）、耐药终点（Endpoint）四个时间点采集肿瘤样本，进行bulk RNA-seq与单细胞转录组测序，分析克隆丰度变化、转录状态演化及遗传变异特征。
结果解读：肿瘤在治疗初期体积缩小，随后逐渐复发；克隆动态分为三类：持留亚群（初始响应治疗后持续存活并扩增）、敏感亚群（治疗后丰度显著下降）、低丰度亚群（丰度波动且不主导耐药肿瘤）；持留克隆起源于预处理的微小亚群，且无基因或通路水平的遗传收敛（置换检验，P>0.05）；单细胞转录组显示，持留克隆在耐药终点形成四个转录状态：应激样、脂质代谢、PI3K信号、ECM重塑。
产品关联：实验所用关键产品：Bio-Serv的BRAF/MEK抑制剂饲料（S7581），Miltenyi Biotec的肿瘤解离试剂盒，10x Genomics的Visium空间转录组平台。
3.4 持留细胞转录状态与调控网络解析
实验目的：鉴定持留细胞的转录特征、调控因子及表观遗传/遗传调控机制。
方法细节：用ClonoCluster算法整合克隆条形码与转录组数据进行聚类分析，通过Enrichr进行功能富集分析，SCENIC分析转录调控网络，MAAPER分析可变多聚腺苷酸化（Alternative Polyadenylation, APA），inferCNV分析拷贝数变异（Copy Number Variation, CNV）。
结果解读：四个持留状态分别富集应激响应、脂质代谢、PI3K信号、ECM重塑通路；转录调控因子分别为ATF4（应激样）、ETV5（脂质代谢）、LEF1（PI3K信号）、ETS1（ECM重塑）；APA分析显示，早期MRD出现3"UTR长短变化，耐药终点呈全局3"UTR缩短；CNV分析显示，脂质代谢和PI3K信号组存在BRAF扩增，且BRAF拷贝数与表达水平正相关（rs>0.44，P<0.001）。
产品关联：实验所用关键产品：R包Seurat v5、ClonoCluster v0.0.1、SCENIC v0.12.1、MAAPER、inferCNV v1.20.0。
3.5 持留细胞的空间组织与细胞间通讯分析
实验目的：解析持留细胞的空间分布模式与细胞间通讯机制。
方法细节：对耐药终点肿瘤进行空间转录组分析，用AUCell注释持留状态，Moran指数分析空间自相关，CellChat推断细胞间通讯网络；设计条形码特异性RNA-FISH探针，验证主导持留克隆与应激样标志物SLC2A1的共定位。
结果解读：脂质代谢和PI3K信号状态共定位于肿瘤边界，ECM重塑状态呈马赛克分布；细胞间通讯分析显示，应激样细胞分泌TGFβ，脂质代谢与PI3K信号状态存在BMP自分泌/旁分泌通讯，ECM重塑细胞激活胶原信号；RNA-FISH显示，主导持留克隆（条形码后缀GTTGAACGACCACAA）与SLC2A1显著共定位（所有斑块尺寸，Fisher精确检验，P<0.014）。
产品关联：实验所用关键产品：ACD的RNAscope Multiplex荧光试剂套装v2（323270）、Hs-SLC2A1探针（423141-C2），Leica的TCS SP8 X WLL共聚焦显微镜。
4. Biomarker研究及发现成果解析
Biomarker定位
本研究涉及的Biomarker包括持留细胞的转录标志物、细胞表面标志物及转录调控因子，筛选与验证逻辑为：基于单细胞转录组差异分析筛选候选标志物→通过bulk RNA-seq验证其在多时间点的表达动态→通过RNA-FISH验证空间共定位→通过TCGA数据库验证临床相关性。
研究过程详述
- 转录标志物：应激样状态标志物SLC2A1、BNIP3，脂质代谢状态标志物FASN、APOE，PI3K信号状态标志物AKT3、FGFR1，ECM重塑状态标志物ECM1、MET。验证方法：单细胞转录组差异表达分析（fold change≥2，细胞占比≥50%，FDR<0.05），bulk RNA-seq动态分析显示标志物在耐药终点表达显著上调（P<0.05，双尾t检验），RNA-FISH验证SLC2A1与持留克隆的空间共定位。特异性与敏感性：SLC2A1在应激样持留细胞中的表达阳性率达80%以上（文献未明确提供具体数值，基于图表趋势推测）。
- 细胞表面标志物：泛持留细胞标志物FXYD3，在所有四个持留状态中均高表达，通过单细胞转录组差异分析验证（fold change≥2，细胞占比≥50%，FDR<0.05）。
- 转录调控因子标志物：ECM重塑状态的调控因子ETS1，通过SCENIC分析鉴定，TCGA数据库生存分析显示其高表达与黑色素瘤患者不良预后相关（P=0.018，Kaplan-Meier分析）。
核心成果提炼
功能关联方面，SLC2A1作为应激样持留细胞的标志物，与克隆扩增及耐药形成直接相关；ETS1作为ECM重塑状态的调控因子，其高表达与黑色素瘤患者总生存期缩短相关（P=0.018，TCGA样本量n=470）；FXYD3作为泛持留细胞表面标志物，有望成为DTP靶向治疗的候选靶点。
创新性方面，本研究首次鉴定了黑色素瘤DTP的四种转录状态及对应标志物，系统解析了其克隆起源与动态演化规律；首次揭示了ECM重塑状态的调控因子ETS1与患者预后的关联；首次通过RNA-FISH验证了持留克隆与标志物的空间共定位，为DTP的临床检测与靶向干预提供了新的生物标志物与理论依据。