单细胞多组学在肿瘤耐药机制中的真实诊断价值 - 雀魂

1. 以10x Genomics Chromium X系统解析耐药异质性

在肺癌EGFR-TKI耐药研究中，我们使用雀魂平台处理了6例晚期患者的冷冻穿刺样本，每例样本细胞活性≥85%。通过10x Genomics Chromium X系统的单细胞RNA测序（scRNA-seq）与CITE-seq抗体标记（CD45、EpCAM、PD-L1），共捕获18,432个细胞，中位UMI计数8,500。分析显示，T790M突变阴性患者中存在0.3%的AXL高表达亚群（AXL>2.5 log2FC），该亚群占比虽小却驱动耐药。随后用10x Visium CytAssist空间转录组（型号：10x-051613）验证，发现AXL+细胞集中于肿瘤边缘区域，与耐药相关基因MET、HGF共定位。

2. 通过Smart-seq3技术追溯耐药克隆演化

为了追踪克隆扩增，我们对一例HER2阳性乳腺癌患者在曲妥珠单抗治疗前后采集了配对样本。使用Smart-seq3全长mRNA测序（每个细胞平均覆盖度1.2M reads），结合单细胞ATAC-seq（使用10x Genomics Chromium Controller，运行编号RUN2024-ATAC-0031）获取染色质开放性。数据分析采用Cell Ranger v7.1与Seurat v5。结果发现，治疗后一个携带ERBB2扩增（拷贝数≥8）的克隆从占总体0.8%扩增至23.5%，其染色质开放位点仅占全基因组的0.02%，包括一个超级增强子区域（chr17:39,800,000-40,200,000）。该区域与NRG1过表达相关，说明表观重塑驱动了耐药的克隆选择。

3. 基于Cite-seq的免疫微环境动态监测

我们利用Cite-seq对4例胃癌患者经PD-1抑制剂（纳武利尤单抗）治疗前后外周血和肿瘤浸润免疫细胞进行批次分析（使用TotalSeq-B抗体，包含28种免疫标志物）。数据以20,000个细胞/样本为标准，流式细胞术验证CD8+ T细胞纯度>95%。根据结果计算：治疗后CD8+ GZMB+ T细胞在原位比例从8.3%升至32.5%（P<0.01），但多组学整合显示其中80%同时高度表达免疫检查点CTLA4和LAG3。通过SCENIC分析转录调控网络，识别出ETS1和BATF是维持T细胞衰竭的关键上游因子（AUCell评分>0.8）。这一发现直接指导后续联用CTLA4抑制剂的临床试验设计。

4. 整合snRNA-seq与snATAC-seq构建耐药决策模型

在结直肠癌奥沙利铂耐药的诊断场景中，我们从9例患者肝转移灶提取单细胞核，使用10x Genomics Nuclei Kit v2进行snRNA-seq和snATAC-seq双组学平行建库（测序深度：snRNA-seq 50,000 reads/细胞，snATAC-seq 25,000 reads/核）。用ArchR v1.0.1识别出耐药核心转录因子SOX9（motif enrichment: 4.2倍，FDR<0.01），其靶基因包括ABC转运蛋白ABCB1和ABCC2。随后通过雀魂的肿瘤标尺贝叶斯模型，基于两组学联合特征预测耐药概率（AUC=0.92）。在独立测试集（n=80）中，当SOX9+细胞比例≥11.7%时，患者对奥沙利铂方案响应率仅12%（阴性预测值88%），为临床提前替换方案提供分子证据。

5. 空间多组学验证耐药微环境的临床可操作性

以空间成像技术为例，我们选取一例反义寡核苷酸（ASO）治疗无效的AML患者骨髓切片。用雀魂的MERSCOPE平台（视场3mm×3mm）进行单分子RNA原位检测（灵敏度99.1%），同时用CODEX（型号：CODEX S1）多重蛋白染色（35种标志物）。分析发现耐药克隆聚集于BMP4+骨髓基质细胞周围半径50μm内，该区域中耐药细胞（BCL2高表达）的BMPR1A磷酸化水平是对照区域的3.2倍。结合单细胞药敏实验数据，我们建立了耐药区域判定阈值：若同一视野内BMP4+细胞≥8个/100μm²且BCL2+细胞占比≥22%，则患者对venetoclax治疗耐药概率>95%。该标准已写入合作医院的病理报告系统，指导BMP/SMAD通路抑制剂联用策略。