基于拟时序肿瘤细胞状态演化的七基因预后模型及TYRP1在黑色素瘤中的多队列验证与免疫微环境重塑研究

《International Journal of Genomics》：Pseudotime-Derived Seven-Gene Prognostic Signature and TYRP1 in Melanoma

　　

背景：黑色素瘤（melanoma）具有高度恶性及复杂的免疫微环境异质性。目前仍缺乏整合肿瘤细胞状态多样性与临床结局的稳健预后框架。方法：研究人员对公共黑色素瘤数据集进行整合性单细胞RNA测序（scRNA-seq）分析，采用Harmony批次校正、无监督聚类和Monocle 3拟时序轨迹重建以刻画肿瘤细胞状态转变。对拟时序动态基因在四个独立队列（GSE19234、GSE22153、GSE65904、TCGA）中进行单因素Cox筛选，随后采用LASSO-Cox回归构建多基因预后模型。验证包括免疫浸润估计、通路活性评分、肿瘤突变负荷（TMB）分析及蛋白水平表达评估。结果：scRNA-seq鉴定出14个主要细胞群，其中肿瘤细胞再聚类解析出7种功能状态：黑色素细胞样、神经嵴样、主要组织相容性复合体II类（MHC-II）高、应激/过渡、增殖、干扰素（IFN）刺激和侵袭性。拟时序分析重建了起源于黑色素细胞样状态的方向性演化轨迹，应激/过渡状态作为中心枢纽。基因模块富集进一步刻画各状态的核心转录程序。一个七基因预后模型（GBP4、ITGB4、GBP2、TNFSF13B、TSC22D3、CIITA、TYRP1）在全部四个队列中实现一致的生存分层（均p < 0.01），并按一致性指数（C-index）排名优于76个已发表特征。高风险患者表现为免疫抑制微环境、促肿瘤通路激活及癌症–免疫循环的广泛抑制。风险评分与TMB在统计学上独立，但二者联合分层显著改善预后分辨率（p < 0.001）。TYRP1作为模型中唯一风险基因，经免疫组织化学（IHC）显示肿瘤特异性蛋白上调，并在三个独立转录组队列中一致预测较差总生存（OS）。结论：本研究提出了黑色素瘤整合单细胞与多队列预后框架，揭示肿瘤细胞状态演化与免疫微环境重塑之间的机制联系，并确立TYRP1作为独立预后生物标志物以指导精准免疫治疗决策。

研究背景与问题：黑色素瘤是致死性最高的皮肤恶性肿瘤，尽管发病率在皮肤癌中居第17位，仍造成显著全球疾病负担；GLOBOCAN 2022估计全球每年约331,722例新发和58,667例死亡，澳大利亚、北美和欧洲等高人类发展指数地区发病率持续上升。以免疫检查点抑制剂（ICIs）为核心的全身治疗显著改善晚期黑色素瘤结局，III期CheckMate 067试验中nivolumab联合ipilimumab在6.5年随访时中位总生存（OS）达72.1个月，显著优于单药和化疗，但约50%患者仍因原发或获得性耐药发生疾病进展。治疗抵抗机制复杂及患者间预后异质性显著，亟需跨队列可泛化的精准预后模型以支持个体化分层和治疗决策。肿瘤内异质性是治疗反应和临床结局差异的核心生物学基础；Tirosh et al.将单细胞RNA测序（scRNA-seq）用于黑色素瘤多细胞生态系统解析，揭示恶性细胞在细胞周期状态、空间分布和耐药程序上的广泛转录异质性；Tsoi et al.在低传代黑色素瘤细胞系中鉴定黑色素细胞样、过渡型、神经嵴样和未分化四种连续分化状态，并显示其与铁死亡敏感性和药物反应相关；近期肿瘤微环境（TME，由肿瘤细胞、免疫细胞和基质等构成）单细胞图谱研究显示高、低风险患者免疫浸润组成和功能状态明显不同，肿瘤细胞状态转变与免疫逃逸程序紧密耦合。但现有研究多停留在静态细胞类型刻画，尚缺乏整合拟时序（pseudotime）动态状态转变与多队列预后建模的系统框架。拟时序轨迹分析可重建肿瘤细胞状态动态演化路径，从连续转录转变中提取动态基因集作为预后候选；基于LASSO-Cox正则化回归的预后模型在多种实体瘤中显示优于传统临床分期，但现有黑色素瘤预后模型普遍存在单队列验证、候选基因生物学来源不清及模型间缺乏头对头比较等限制。TYRP1是黑色素生成通路关键酶，在黑色素细胞和黑色素瘤细胞中特异表达，正常非黑色素细胞组织不表达；早期研究提示高TYRP1 mRNA与较短远处转移无生存和OS相关，并与Breslow肿瘤厚度正相关；TYRP1又被重新评估为CAR-T治疗靶点，其在包括PD-1阻断抵抗的免疫排斥肿瘤等黑色素瘤亚型中高表达，凸显治疗靶点和预后标志物双重价值，但尚缺乏在预后模型框架内对TYRP1作为核心风险基因的多队列系统验证。基于此，研究人员整合公共黑色素瘤scRNA-seq数据，应用Monocle 3拟时序轨迹分析识别与肿瘤细胞状态转变相关的动态基因集，在四个独立队列中以单因素Cox回归筛选出24个预后一致基因，再经LASSO-Cox正则化回归精简为七基因预后风险评分模型；该模型在四队列中显示稳健生存分层，并通过与76个已发表黑色素瘤预后特征头对头比较确认竞争优势；研究人员还从免疫微环境重塑、肿瘤突变负荷（TMB）关联和癌症–免疫循环抑制等维度阐释模型生物学基础，并对模型唯一风险基因TYRP1进行蛋白表达和多队列转录组独立验证。该研究为黑色素瘤精准预后分层和免疫治疗候选选择提供新的分子工具和理论依据。

关键技术与方法：研究使用公共scRNA-seq数据集GSE277165（10例原发黑色素瘤患者肿瘤组织样本），经Seurat质控、Harmony批次校正、无监督聚类和UMAP可视化；以Monocle 3将黑色素细胞样群设为根节点重构拟时序，graph_test识别动态基因并行共表达模块分析。多队列验证纳入GSE19234（N = 44）、GSE22153（N = 54）、GSE65904（N = 210）和TCGA-SKCM（N = 469），经limma去批次后行单因素Cox筛选和十折交叉验证LASSO-Cox建模；采用timeROC、C-index比较、IOBR整合七种算法估计免疫浸润、ESTIMATE评分、GSVA、maftools处理体细胞突变，以及HPA免疫组织化学（IHC）验证TYRP1蛋白表达。

研究结果：
3.1 单细胞转录组图谱揭示黑色素瘤TME的细胞组成：通过整合scRNA-seq、质控、Harmony校正、无监督聚类和UMAP，研究人员鉴定14个主要细胞群，并以经典标记基因注释肿瘤细胞（MLANA、S100B、PMEL）、癌症相关成纤维细胞（CAFs，ACTA2）、T细胞（CD3D）、巨噬细胞（LYZ）、B细胞（MS4A1）、自然杀伤（NK）/NKT细胞（NKG7）、内皮细胞（PECAM1）、角质形成细胞（KRT14）、肥大细胞（TPSB2）和增殖细胞（MKI67），标记表达与聚类注释一致，支持细胞类型分类稳健性。
3.2 肿瘤细胞再聚类描绘七种不同功能状态：对肿瘤细胞独立再聚类获得七种功能状态，标记分析显示黑色素细胞样状态高表达PMEL、TYRP1、DCT和CTNNB1；神经嵴样状态特异表达BCAN和CRYAB；MHC-II高状态富集HLA-DRA、CD74和A2M；应激/过渡状态上调FOS、JUN和ATF3；增殖状态由TOP2A、PCNA和MCM4定义；IFN刺激状态显示CXCL9、CXCL10和STAT1升高；侵袭性状态特异上调VGF、ITGA3和CCND1，代表性标记基因表达图进一步验证上述注释。
3.3 拟时序轨迹分析重建肿瘤细胞状态动态演化：Monocle 3以黑色素细胞样群为根节点重建连续状态转变轨迹，黑色素细胞样细胞集中于近零拟时序，神经嵴样细胞拟时序逐渐升高，应激/过渡状态拟时序分布最广并位于多状态连接的中心枢纽，MHC-II高细胞中位拟时序最高并位于主轨迹末端；增殖、IFN刺激和侵袭性状态聚集于低拟时序并在UMAP占据不同位置，原文提示这些状态可能由特定微环境线索驱动，而非主分化轴上的顺序节点。
3.4 基因模块与通路富集分析揭示状态转变的核心转录程序：对拟时序动态基因的共表达分析鉴定三个肿瘤状态相关模块；Module_1在侵袭性和黑色素细胞样状态最活跃，富集TNFα Signaling via NF-κB、上皮–间质转化（EMT）、凋亡、缺氧、p53通路和UV反应；Module_2在IFN刺激状态选择性激活，仅富集干扰素α反应和干扰素γ反应通路，表明免疫微环境持续干扰素信号是该状态关键上游驱动。
3.5 多队列预后基因筛选与LASSO-Cox模型构建：整合四队列并去批次后，PCA显示样本在降维空间均匀分布（Dim1 = 10.9%，Dim2 = 8.6%）；以拟时序动态基因为候选，四队列单因素Cox筛选出24个预后一致基因，森林图显示多数基因为保护性（HR < 1），TYRP1为唯一显著风险基因（HR > 1）；CNV分析显示6号染色体多个HLA家族成员频繁扩增（AMP），9号染色体TYRP1主要缺失（DEL）；LASSO-Cox在最优λ = 0.023时保留7个非零系数基因（GBP4、ITGB4、GBP2、TNFSF13B、TSC22D3、CIITA、TYRP1），十折交叉验证显示该λ对应最小预测误差。
3.6 预后模型在独立队列中显示稳健生存分层与预测性能：Kaplan–Meier分析显示四队列高风险组OS均显著短于低风险组（均p < 0.01）；时间依赖ROC显示TCGA队列1、3、5年AUC为0.62、0.66、0.67，GSE19234队列1、2、3年为0.93、0.81、0.91，GSE22153为0.81、0.82、0.76，GSE65904为0.65、0.68、0.67；与76个已发表模型C-index头对头比较显示本模型在四队列均居前列，七基因表达谱主成分分析（PCA）显示高、低风险患者在PC1–PC2空间明显分离。
3.7 预后风险模型与TIME激活状态密切相关：整合七种算法免疫浸润估计显示，低风险患者几乎所有算法框架下CD4⁺ T细胞、CD8⁺ T细胞、B细胞、NK细胞和树突状细胞（DC）等效应亚群浸润更高，呈免疫炎症型；高风险肿瘤免疫浸润整体耗竭，呈免疫排斥或“冷肿瘤”状态；低风险组MHC I/II类抗原呈递分子、共抑制检查点分子（包括HAVCR2、LAG3、PDCD1、TIGIT、CTLA4）及多数共刺激分子广泛上调，VTCN1和TNFSF9组间差异不显著；风险评分与StromaScore（r = ?0.38，q = 0）、ImmuneScore（r = ?0.55，q = 0）、ESTIMATEScore（r = ?0.52，q = 0）负相关，与TumorPurity正相关（r = 0.53，q = 0），提示高风险患者肿瘤细胞纯度较高、免疫和基质浸润减少并呈免疫抑制微环境。
3.8 风险模型与TMB及免疫通路活性的关联：风险评分与多数癌症相关通路显著正相关，包括Cell_cycle、DNA_replication、Mismatch_repair、Proteasome和Spliceosome，与癌症–免疫循环全部7步显著负相关，尤其T细胞招募、浸润、识别和杀伤步骤；突变图谱显示TTN（67%）、MUC16（62%）和BRAF（51%）为最常见突变基因，高风险患者总体突变频率更高；高、低风险组TMB差异不显著（p = 0.21），风险评分与TMB无有意义线性关联（R = 0.033，p = 0.49），提示风险评分捕获的预后信息独立于TMB；按风险评分和TMB中位数联合分层时，H-TMB + 低风险组OS最佳，L-TMB + 高风险组预后最差，四组差异显著（p < 0.001），联合分层较单变量提高预后分辨率。
3.9 TYRP1作为唯一风险基因的蛋白表达与预后验证：HPA数据库IHC显示TYRP1在正常外阴皮肤表皮细胞无染色，而黑色素瘤肿瘤细胞呈强、弥漫阳性，定位于细胞质和细胞膜，提示肿瘤特异性蛋白上调；TCGA（N = 469）、GSE19234（N = 44）和GSE65904（N = 210）中，TYRP1高表达患者OS均显著较短（log-rank p < 0.0001、p = 0.00014、p = 0.008），支持TYRP1为独立验证的预后生物标志物。

讨论与结论：讨论指出，该研究建立从单细胞转录组异质性到多队列预后验证再到免疫微环境功能解析的整合框架，从拟时序解析的恶性细胞状态出发识别动态基因、经多队列单因素Cox筛选和LASSO-Cox缩减为七基因评分，区别于纯数据驱动的机器学习特征。拟时序主轨迹从黑色素细胞样状态经神经嵴样去分化至MHC-II高免疫调节终点，原文提示分化丢失与免疫逃逸获得之间存在耦合，并与表型可塑性驱动适应性免疫抵抗模型一致；应激/过渡状态分布最广，是连接多状态的中枢节点而非单一线性路径中间站，原文提示靶向该状态可能具有阻断塑性转变的广谱潜力；增殖、IFN刺激和侵袭性状态在拟时序独立聚集，原文提示其出现由特定微环境线索驱动而非主分化轴顺序进展，这一“主干–分支”拓扑可能是患者对免疫和靶向治疗异质反应的细胞基础。七基因模型各成员具有明确免疫生物学依据：GBP2和GBP4属干扰素诱导GTPase家族，其升高不仅提示持续IFN-γ信号输入，还与炎症小体激活和MHC-I呈递增强相关，低风险患者富集反映免疫监视活跃；CIITA为MHC-II主转录调控因子，使肿瘤细胞直接参与抗原呈递，支持肿瘤细胞非经典激活CD4⁺ T细胞塑造适应性免疫；TSC22D3（GILZ）在高风险患者上调，作为糖皮质激素信号核心效应分子，可通过抑制NF-κB和AP-1通路广泛抑制T细胞活化，构成内源性、非检查点依赖的免疫抑制机制；TNFSF13B（BAFF）指向TME中B细胞功能状态，近期研究提示肿瘤引流淋巴结BAFF信号激活与三级淋巴结构形成及抗肿瘤体液免疫相关；ITGB4提示侵袭表型与免疫逃逸的结构关联，整合素β4介导的基底膜锚定信号可激活PI3K-AKT轴并下调肿瘤细胞多种免疫识别分子表面表达；七基因功能汇聚表明模型捕获肿瘤细胞“免疫可见性”的综合状态。模型与TIME高度一致，低风险患者免疫炎症表型伴效应T细胞富集和PD-1、LAG-3、TIGIT、HAVCR2等检查点广泛上调，原文认为这并非简单免疫激活，而是功能性耗竭样且可经PD-1/PD-L1阻断逆转的状态，与ICB反应预测特征重叠；高风险患者高肿瘤纯度、基质和免疫细胞浸润双重耗竭，癌症–免疫循环中T细胞招募、识别和杀伤等多步骤系统抑制，原文提示其抵抗机制可能超出单一检查点分子，需联合策略。风险评分与TMB无线性相关且联合分层预后分辨率更优，原文支持二者分别反映肿瘤固有免疫编程和新生抗原负荷两个不同生物学侧面，整合有助于细化ICB患者选择。TYRP1作为唯一风险基因具有双重转化价值：蛋白和转录组多队列一致验证确立其独立预后标志物地位，其在免疫排斥高风险肿瘤中选择性高表达符合新兴CAR-T靶抗原标准，即高肿瘤表达伴低正常组织毒性，原文提示TYRP1定向过继细胞治疗可能成为本模型定义的高风险、ICB获益有限患者的补充策略。局限性包括依赖回顾性公共数据集、需前瞻临床队列验证；单细胞数据来自特定患者人群，细胞状态图谱在不同种族和疾病分期中的泛化性需系统评估；基因模块与细胞状态因果关系需功能实验阐明；未来应在前瞻队列验证风险评分对ICB疗效的预测价值，探索靶向应激/过渡状态的干预策略，并在体内系统评估TYRP1定向治疗联合免疫检查点阻断的协同潜力。研究结论（译）：本研究提出黑色素瘤整合单细胞与多队列预后框架，揭示肿瘤细胞状态演化与免疫微环境重塑之间的机制联系，并确立TYRP1作为独立预后生物标志物以指导精准免疫治疗决策。