单细胞蛋白质活性分析揭示了MYOD1突变横纹肌肉瘤中异常的肌生成和IGF2-PI3K通路依赖性。
Single-cell protein activity analysis reveals aberrant myogenesis and IGF2-PI3K pathway dependencies in MYOD1-mutant rhabdomyosarcoma.
1. 领域背景与文献
发表期刊：Science Advances；影响因子：14.9（2023年）；研究领域：肿瘤学，罕见软组织肉瘤分子机制。
领域共识：横纹肌肉瘤是儿童和青年人群高发的软组织肉瘤，MYOD1突变梭形细胞/硬化性横纹肌肉瘤（SRMS）是一种超罕见的临床侵袭性亚型，其定义性特征为MYOD1基因L122R热点突变，该突变会破坏骨骼肌分化过程中MYOD1的正常转录因子功能。现有临床数据显示该疾病对常规化疗响应差，5年生存率不足20%，复发率高，缺乏有效的靶向治疗手段，长期生存极为罕见。近年来单细胞RNA测序研究已经明确，融合阳性横纹肌肉瘤和融合阴性横纹肌肉瘤均存在反映肌发育层级的复杂肿瘤内异质性，这种异质性与治疗敏感性密切相关，但由于MYOD1突变SRMS发病率极低，此前仅2例被纳入单细胞研究，其细胞状态组成、转录调控程序和治疗靶点始终未得到系统解析。
目前仅30%~40%的MYOD1突变SRMS携带可靶向的PI3K/AKT/mTOR通路共突变，仅通过基因组学无法为所有病例识别治疗靶点，同时MYOD1 L122R突变如何重编程转录调控网络、驱动肿瘤进展和耐药的机制也不明确。该领域亟需针对该罕见亚型系统解析其分子特征，发掘机制明确的治疗策略。本研究针对这一研究空白，开展了目前最大样本量的MYOD1突变SRMS单细胞转录组研究，系统解析了其细胞状态架构和可靶向调控依赖性，为该难治疾病提供了新的治疗方向。
2. 文献综述解析
本文作者围绕横纹肌肉瘤领域现有研究进展和局限性展开梳理，按疾病亚型和研究技术两个维度对现有成果进行分类评述。现有研究已经证实，调控网络结合单细胞蛋白活性分析的技术框架可以有效从稀疏的单细胞转录组数据中识别驱动肿瘤细胞状态的核心调控因子，这类调控因子富集于肿瘤必需基因和可靶向的非癌基因依赖，即使在没有可检测转录本或可靶向基因组突变的情况下也能识别潜在治疗靶点，该技术已经在多种肿瘤类型中得到验证。针对横纹肌肉瘤的现有研究已经明确不同亚型存在反映肌发育层级的肿瘤内异质性，这种异质性对治疗选择有重要指导意义，MYOD1突变SRMS的定义性突变已经得到明确，基因组层面的突变谱也有初步报道，证实PI3K/AKT/mTOR通路改变是常见的共突变。
但现有研究存在明确的局限性：首先，由于MYOD1突变SRMS发病率极低，现有大样本单细胞研究几乎未纳入该亚型，其细胞状态架构和转录调控程序始终未得到系统解析；其次，仅基于基因组学识别靶点的策略无法覆盖不存在通路突变的病例，相当比例的患者没有可靶向的基因组改变，现有方法无法为这类患者提供治疗选择；最后，现有研究仅证实MYOD1 L122R突变会干扰正常肌分化，但该突变如何重编程转录调控网络、驱动肿瘤进展的机制尚不明确。
本研究的创新价值十分明确，这是首次针对MYOD1突变SRMS完成大样本单细胞转录组和蛋白活性分析，构建了目前该领域最大的分子数据集，首次系统解析了该亚型独特的细胞状态架构，发现了不依赖突变状态的IGF2-PI3K通路可靶向依赖性，明确了突变MYOD1重编程肌分化调控网络的具体机制，弥补了该领域的核心研究空白，为该超罕见难治疾病提供了经过临床前验证的新治疗策略，也为其他罕见肿瘤的靶点发掘提供了可复制的研究范式。
3. 研究思路总结与详细解析
本研究的核心研究目标是解析MYOD1 L122R突变SRMS的肿瘤内异质性、转录调控程序和可靶向治疗依赖，核心科学问题是突变MYOD1如何重编程肌分化程序并驱动肿瘤进展，存在哪些可靶向的调控依赖性。整体技术路线为：收集患者来源肿瘤样本→单核RNA测序结合调控网络分析识别保守肿瘤细胞状态→通过配受体分析和比较组学识别核心调控通路→在独立bulk队列、患者来源异种移植（PDX）模型、类器官中验证保守性→体内外功能验证靶向治疗效果，形成了从发现到验证的完整研究闭环。
3.1 单细胞转录组测序与恶性细胞鉴定
实验目的：获得MYOD1 L122R突变SRMS的单细胞基因表达谱，区分恶性肿瘤细胞与非肿瘤细胞，明确肿瘤内细胞组成特征。
方法细节：采集6例经新辅助化疗后的MYOD1 L122R突变SRMS患者冰冻肿瘤样本，采用10x Genomics平台进行单核RNA测序，通过SingleR算法注释细胞类型，以样本匹配的免疫细胞作为正常参考，使用InferCNV识别存在大规模拷贝数变异的恶性肿瘤细胞，过滤非恶性细胞后开展后续分析，同时对配对样本进行苏木精-伊红染色和MYOD1免疫组化（IHC）分析明确病理特征。
结果解读：最终获得177051个高质量恶性细胞核，占总数据集的92%，非恶性基质和免疫细胞占比为6%~27%（n=6），这是目前为止该罕见亚型最大的单细胞转录组数据集。病理检测结果显示，所有样本MYOD1均呈弥漫强阳性，肌细胞生成素表达罕见至阴性，符合该亚型的病理特征；治疗后获得部分缓解的样本可见治疗诱导的细胞分化，伴随MYOD1表达丢失，与临床预后特征一致。该分析确认了研究样本的病理和分子特征，为后续分析提供了高质量的基础数据。文献未提及具体测序实验产品，领域常规使用10x Genomics单细胞建库试剂盒。
3.2 单细胞蛋白活性推断与细胞状态鉴定
实验目的：基于调控网络推断单细胞水平的调控蛋白活性，识别MYOD1 L122R突变SRMS中保守的肿瘤细胞状态。
方法细节：首先使用ARACNe-AP算法针对每个样本从头构建样本特异性的单细胞基因调控网络，获得每个调控因子的靶基因调控子；然后使用metaVIPER算法，基于靶基因的差异表达推断每个细胞中2336个调控蛋白的活性；随后采用分辨率优化的Louven算法对蛋白活性进行聚类，通过viperSimilarity计算亚簇之间的调控蛋白活性相似性，结合CytoTRACE分化潜能评分、细胞周期评分和肌源性调控因子活性，将亚簇整合为保守的细胞状态，进一步与其他横纹肌肉瘤亚型的调控网络进行比较分析。
结果解读：聚类分析共得到20个患者特异性亚簇，最终整合为3个在所有患者中均存在的保守细胞状态：祖细胞状态（68265个细胞核）、过渡状态（82053个细胞核）和分化状态（26292个细胞核），三个细胞状态沿分化连续谱分布。祖细胞状态的CytoTRACE分化潜能评分最高，提示分化程度最低，显著高于过渡和分化状态（****P<0.0001，n=177051），该状态间充质干细胞基因特征富集，缺乏经典肌源性转录因子活性；过渡状态分化潜能居中，增殖细胞比例最高，部分激活肌源性调控因子；分化状态CytoTRACE评分最低，分化程度最高，激活肌细胞基因特征和晚期肌源性效应因子，但MYOD1等关键调控因子活性仍然较低，提示仅能达到部分分化，无法完成正常终末分化。与融合阳性/阴性横纹肌肉瘤的野生型MYOD1不同，MYOD1活性在MYOD1突变SRMS的祖细胞状态中富集，随分化进程降低，这种反转的活性模式提示突变MYOD1改变了正常的调控轨迹。比较调控网络分析显示，突变MYOD1与野生型MYOD1的靶基因仅存在11个重叠，几乎完全不同，突变MYOD1的靶基因富集于PI3K/AKT/mTOR信号和增殖、生物合成相关通路，而野生型MYOD1靶基因富集于经典肌生成通路，证实突变MYOD1将调控网络从肌分化重编程为增殖相关通路。文献未提及具体分析产品，领域常规使用ARACNe、metaVIPER开源分析工具。
3.3 核心调控通路识别与保守性验证
实验目的：识别细胞状态特异性的核心调控通路，验证其在独立队列和PDX模型中的保守性。
方法细节：通过Pearson相关性分析识别与分化潜能显著相关的细胞状态特异性核心调控因子；使用CellChat分析配受体相互作用，推断细胞间信号交流；在24例MYOD1突变SRMS的bulk RNA-seq队列、2例PDX模型的snRNA-seq数据中，使用相同的metaVIPER分析流程验证细胞状态和核心调控通路的保守性；对20例MYOD1突变SRMS进行MSK-IMPACT靶向测序，分析基因组改变与调控通路的关联。
结果解读：IGF2是祖细胞状态中活性最高的调控因子，且与分化潜能显著正相关；祖细胞高表达IGF2但不表达IGF1R和IGF2R，而过渡和分化状态IGF2活性低但高表达IGF1R和IGF2R。CellChat分析证实，祖细胞是IGF2配体的主要发送者，过渡和分化细胞是主要接收者，形成了祖细胞向分化细胞的旁分泌IGF2-IGF1R信号轴。与正常骨骼肌相比，所有细胞状态的IGF1R/PI3K/AKT/mTOR通路组分均显著激活（Bonferroni校正P≤10−5），在过渡和分化状态富集更显著，祖细胞还存在PTK2、SYK等独特的可激活调控因子。PDX模型的snRNA-seq分析证实，PDX保留了患者肿瘤的细胞状态架构和调控程序，IGF2-PI3K通路激活特征与患者肿瘤完全一致。24例bulk队列的分析显示，bulk样本可分为祖细胞样和分化样亚组，均存在IGF1R-PI3K/AKT/mTOR通路的显著激活。20例患者的基因组分析显示，IGF2扩增和PI3K/AKT/mTOR通路改变是最常见的共突变，这些突变几乎完全互斥，且在治疗后样本中富集，提示不同的基因组改变最终汇聚到相同的IGF2-PI3K通路激活，与转录调控分析结果一致。文献未提及具体实验产品，领域常规使用CellChat分析工具、MSK-IMPACT靶向测序试剂盒。
3.4 体内外功能验证治疗敏感性
实验目的：功能验证IGF2-PI3K/AKT/mTOR通路是MYOD1突变SRMS的可靶向治疗依赖，评估靶向联合化疗的体内外疗效。
方法细节：在PDX来源的细胞系中，通过外源性IGF2刺激后免疫印迹检测AKT和ERK磷酸化，验证通路激活；在PDX来源细胞和患者来源三维类器官（PDXO）中进行药物筛选，评估不同靶点抑制剂的杀伤效果和联合用药的协同效应；在MYOD1突变SRMS PDX模型中开展体内疗效实验，分为安慰剂组、双PI3K/mTOR抑制剂LY3023414单药组、伊立替康+替莫唑胺化疗组、LY3023414联合化疗组，每组8只小鼠，评估肿瘤体积控制和无进展生存。
结果解读：外源性IGF2刺激可诱导PDX来源细胞中AKT和ERK的磷酸化，证实IGF2可激活下游PI3K和MAPK通路。药物筛选显示，MYOD1突变SRMS细胞对IGF1R、PI3K、AKT、mTOR抑制剂敏感，其中AZD8055（mTOR抑制剂）半数抑制浓度（IC50）为0.11~0.42 μM（n=3），LY3023414（双PI3K/mTOR抑制剂）IC50为0.27~0.72 μM（n=3），而对常规化疗药物伊立替康和替莫唑胺耐药，IC50分别大于29 μM和100 μM（n=3）；联合用药分析显示LY3023414与化疗具有协同效应，三药联合的协同效应最显著。PDXO药物筛选显示，PI3K/AKT/mTOR抑制剂和IGF1R抑制剂的中位响应百分位均低于预设的25%灵敏度阈值，证实该通路是明确的治疗依赖性。体内实验显示，所有治疗组的肿瘤体积控制均显著优于安慰剂组（P<0.01，n=8），LY3023414单药即可获得显著的抗肿瘤活性，联合化疗组达到客观肿瘤退缩，肿瘤体积控制显著优于化疗单药（P=0.005，n=8），所有治疗组无进展生存均显著优于安慰剂（P=7.7×10−5，n=8）。实验所用关键产品：重组IGF2（Sigma-Aldrich，#SRP3070）；磷酸化-AKT抗体（Cell Signaling Technology，#4060T）；磷酸化-ERK1/2抗体（Cell Signaling Technology，#9102）；LY3023414（MedChemExpress，#HY-12513）。
4. Biomarker研究及发现成果解析
本研究识别的核心生物标志物为肿瘤IGF2-PI3K/AKT/mTOR通路活性，属于信号通路活性类预测生物标志物，可用于预测MYOD1突变SRMS对IGF1R-PI3K/mTOR抑制剂的敏感性，筛选验证逻辑为：单细胞调控网络分析差异筛选→独立bulk转录组队列验证→PDX模型和类器官验证→体内功能验证，逻辑链条完整。
该生物标志物的来源为MYOD1突变SRMS肿瘤细胞的转录调控蛋白活性，研究通过metaVIPER算法基于肿瘤细胞基因表达推断该通路核心蛋白的活性，在24例bulk样本中验证了该通路的普遍激活，无论样本是否存在IGF2扩增或PI3K/AKT通路突变，该通路均存在转录水平的激活。功能验证结果显示，该通路激活的肿瘤对IGF1R-PI3K/mTOR抑制剂的敏感性显著高于其他类型药物，体内联合化疗可使肿瘤获得客观退缩，预测价值明确。
核心成果方面，本研究首次在MYOD1突变SRMS中证实，IGF2-PI3K/AKT/mTOR通路是不依赖突变状态的普遍可靶向依赖性，祖细胞通过旁分泌IGF2驱动过渡和分化细胞的通路激活，该通路激活是MYOD1突变SRMS化疗耐药的重要机制，抑制该通路可显著增强化疗的抗肿瘤效果。在PDX模型中，联合治疗对比单纯化疗的无进展生存风险比为0.32（P=0.01，n=8），证实该策略的显著疗效。该发现为超罕见难治的MYOD1突变SRMS提供了经过临床前验证的新治疗策略，也证明了调控网络分析在发掘罕见肿瘤非突变依赖治疗靶点中的价值，为其他罕见肿瘤的研究提供了范式参考。