《International Journal of Genomics》：Integrative Single-Cell and Bulk Transcriptomics Identifies a Proliferative Malignant State and Supports a Functional Contribution of TROAP in Osteosarcoma
编辑推荐：
背景：骨肉瘤含有异质性的恶性细胞状态，其临床相关性难以从小型单细胞队列中推断。
方法：研究人员重新分析了来自GSE162454的六例未经治疗的骨肉瘤，鉴定了经拷贝数变异（CNV）支持的恶性成骨细胞，并利用患者配对伪批量分析和程序选择敏感性测试表征了一种复发性增
背景：骨肉瘤含有异质性的恶性细胞状态，其临床相关性难以从小型单细胞队列中推断。
方法：研究人员重新分析了来自GSE162454的六例未经治疗的骨肉瘤，鉴定了经拷贝数变异（CNV）支持的恶性成骨细胞，并利用患者配对伪批量分析和程序选择敏感性测试表征了一种复发性增殖性/有丝分裂状态。研究人员评估了一个包含50个基因的模块与来自GSE21257的53例治疗前肿瘤中转移发生的关联。候选基因排名在透明加权场景下进行了审核，并在HOS和U2OS细胞中使用两种独立的小干扰RNA（siRNA）对TROAP进行了敲低。
结果：在52,489个输入条形码中，保留了39,063个单态细胞，其中包括10,456个经CNV支持的恶性成骨细胞。M6由725个细胞组成，在所有六个肿瘤中均被检测到。它在八个程序选择场景中的六个场景中排名第一，但当G2M检查点和E2F靶点被共同排除后降至第五位，这支持了其作为增殖性/有丝分裂状态的解释。其加权模块在转移性肿瘤中更高（p = 0.0147；表观AUC = 0.704），但在巨噬细胞调整后该关联减弱（OR = 1.646，95% CI 0.649–4.170；p = 0.2937），而巨噬细胞评分仍与转移呈负相关（OR = 0.251，95% CI 0.081–0.773；p = 0.0161）。TROAP在M6中显著富集，并得到配对伪批量分析的支持，但它并非在所有加权场景下都是排名最高的候选基因。在两种细胞系中，TROAP敲低均减少了集落形成、总体伤口闭合以及Transwell迁移和侵袭。
结论：在六例骨肉瘤中可检测到一种复发性增殖性/有丝分裂恶性状态，其产生的bulk表达程序与转移发生相关，尽管这种关联对巨噬细胞组成敏感。TROAP被支持作为骨肉瘤细胞表型的功能贡献者，而非该状态的既定驱动因素或独立验证的临床标志物。
论文解读文章
一、研究背景与研究问题
骨肉瘤（Osteosarcoma）是儿童和青少年中最常见的原发性恶性骨肿瘤。尽管多药化疗联合手术切除显著改善了局限性疾病的预后，但转移性或复发性骨肉瘤的预后仍然极差，且缺乏有效的治疗进展。该领域的生物学障碍不仅在于缺乏候选基因，更在于骨肉瘤基因组具有结构复杂性、广泛的拷贝数改变（Copy-Number Alteration, CNA）、染色体碎裂以及显著的瘤间和瘤内异质性，这使得难以确定单一的复发性起始病变。因此，一个有临床价值的靶点可能编码于一种复发性恶性细胞状态，而非普遍存在的DNA改变。
近年来，单细胞和空间转录组学研究开始揭示这一问题，发现人骨肉瘤中包含恶性成骨细胞以及与多种免疫、基质细胞相互交织的复杂生态系统。然而，存在两个关键的转化缺口：第一，在小型单细胞队列中识别的恶性亚群可能仅是患者特异性的扩增，而非可重复的状态；第二，一个定义状态的标志物可能在统计上显著，却未必与患者结局相关或对恶性行为有功能贡献。为此，本研究旨在通过整合单细胞发现、跨队列bulk关联分析与体外功能实验，填补这些空白。
二、研究设计与技术方法
研究人员采用了一套严谨的多阶段整合策略，主要关键技术方法包括：
- 1.
单细胞转录组数据分析：对公共数据库GSE162454中的6例未经治疗的原发性骨肉瘤标本进行重新分析，经过严格的质量控制、标准化（SCTransform）、数据整合（Harmony）和聚类，构建了高质量的骨肉瘤细胞图谱。
- 2.
恶性细胞鉴定：利用inferCNV算法，以免疫和内皮细胞为参考，基于转录组衍生的拷贝数变异（CNV）信号，从成骨细胞/间充质谱系中识别出CNV支持的恶性成骨细胞。
- 3.
恶性细胞重聚类与状态定义：对CNV支持的恶性细胞进行独立重聚类，并通过计算六种Hallmark程序（如E2F靶点、G2M检查点等）的综合评分，结合留一法敏感性分析，定义了复发性增殖性/有丝分裂M6状态。
- 4.
跨队列Bulk关联分析：将源自单细胞分析的50个基因模块应用于独立的bulk转录组队列GSE21257（含53例治疗前活检标本），评估其与5年内转移发生的关联，并采用嵌套交叉验证的惩罚逻辑回归（LASSO）评估模型性能。
- 5.
候选基因优先排序与功能验证：通过涵盖单细胞特异性、bulk效应、选择稳定性等多维证据的透明加权敏感性分析，筛选出高优先级候选基因TROAP，并在HOS和U2OS骨肉瘤细胞系中使用两种独立siRNA进行敲低，通过集落形成、划痕愈合和Transwell实验检测其对细胞表型的影响。
三、主要研究结果
1. 高质量图谱解析了未经治疗骨肉瘤的细胞生态系统
通过严格质量控制，最终保留了39,063个高质量单态细胞，成功解析出髓系细胞、成骨细胞/间充质细胞、T/NK细胞、成纤维细胞、B/浆细胞、内皮细胞和周细胞等主要细胞类型，证实了骨肉瘤是一个多细胞生态系统。
2. 表达衍生的CNV模式支持一个大的恶性成骨细胞区室
转录组衍生的CNV负荷主要集中在成骨细胞/间充质细胞中。应用参考细胞第95百分位数CNV阈值并结合成骨细胞谱系身份，共鉴定出10,456个CNV支持的恶性成骨细胞，这些细胞在不同恶性亚簇间显示出广泛的染色体水平偏差。
3. 恶性细胞重聚类识别出一个复发性增殖性/有丝分裂M6状态
对恶性细胞独立重聚类产生11个亚簇（M0-M10）。其中M6具有最高的原始六程序综合评分，包含725个细胞，在所有六例患者中均被检测到。留一程序敏感性分析显示，M6的排名高度依赖于G2M检查点和E2F靶点程序，支持其作为一个协调的增殖性/有丝分裂状态。
4. M6程序围绕有丝分裂和染色体分离组织
患者配对的伪批量差异表达分析鉴定出大量在M6中上调的基因，包括UBE2C、CENPF、TOP2A等。Hallmark富集分析最强的是G2M检查点、E2F靶点和有丝分裂纺锤体，而炎症反应等通路则呈负富集。探索性拟时序分析展示了候选基因沿M2-M8-M6谱系的协调表达变化。
5. M6衍生程序与GSE21257队列中的转移发生相关
在GSE21257队列中，M6衍生的加权模块评分在5年内发生转移的患者中显著更高（Wilcoxon p = 0.0147），表观AUC为0.704。然而，在单基因水平上，只有ASPM在全基因组FDR校正后仍显著，TROAP处于临界值。嵌套交叉验证的LASSO模型产生了适度的AUC（0.659）。
6. 透明的多域证据支持将TROAP用于功能随访
TROAP在M6中高度富集，在伪批量分析中显著上调，且在GSE21257中与转移呈正相关。虽然它在所有统计指标上并非最强（如KIF20A和ASPM排名更高），但其在多种透明加权场景下均保持在高优先级组中，且具有合理的生物学可解释性和有限的骨肉瘤特异性表征，因此被选作功能验证对象。重要的是，M6模块与转移的关联在调整巨噬细胞评分后被削弱，表明该bulk关联对肿瘤微环境组成敏感。
7. TROAP敲低减少克隆形成、伤口闭合、迁移和侵袭
在HOS和U2OS两种细胞系中，两种独立的siRNA均有效降低了TROAP蛋白水平。TROAP敲低显著减少了集落形成数量、24小时总体伤口闭合率以及Transwell迁移和侵袭的细胞数量（所有比较p < 0.001）。
四、讨论与结论总结
本研究的核心贡献在于建立了一条可追溯的证据链：从鉴定可重复的恶性细胞状态，到在独立队列中评估其临床关联，再到聚焦于一个候选基因进行功能测试。研究表明，M6代表一种协调的增殖性/有丝分裂状态，其在所有六例标本中的复发是支持性的，但由于发现队列较小，不能确立其在骨肉瘤中的广泛流行性。跨队列bulk关联分析显示该程序与转移相关，但bulk数据无法定位信号至特定的M6细胞，且该关联在调整巨噬细胞组成后不再显著，这强调了未来研究中应同时测量恶性和微环境轴的重要性。
TROAP被优先考虑是基于其多维证据的一致性，而非作为唯一的统计最优项。功能实验证实了TROAP对骨肉瘤细胞克隆形成和运动相关表型的贡献，但并未证明其驱动了M6状态或介导了与转移的bulk关联。研究明确指出，TROAP更应被视为一种功能贡献者，而非已确定的驱动因子或独立的临床标志物。
研究结论：整合单细胞发现、患者感知复制、跨队列bulk关联和聚焦扰动，在六例人骨肉瘤中鉴定出一种复发性增殖性/有丝分裂恶性状态。其衍生程序在巨噬细胞调整前与GSE21257中的转移发生相关。TROAP被选为高优先级候选之一，其敲低减少了HOS和U2OS细胞的克隆形成、总体伤口闭合以及Transwell迁移和侵袭。这些发现支持TROAP作为骨肉瘤细胞表型的功能贡献者，而其机制和临床作用有待在更大、更直接的模型中检验。
