HAUS相关预后特征在肝肝细胞癌中的开发与验证：从泛癌分析到单细胞分辨率的多组学整合

《Journal of Hepatocellular Carcinoma》：Development and Validation of a HAUS-Related Prognostic Signature in Liver Hepatocellular Carcinoma: A Multi-Omics Integration from Pan-Cancer Analysis to Single-Cell Resolution

背景：HAUS家族蛋白（HAUS1–HAUS8）对有丝分裂纺锤体微管成核（microtubule nucleation）至关重要。尽管HAUS失调与肿瘤进展相关，但其致癌功能是否具有保守性或组织特异性尚不明确。因此，需要进行泛癌分析以识别通用的HAUS驱动因子和依赖情境的成员，用于精准肿瘤学。方法：研究人员使用TCGA（The Cancer Genome Atlas）数据系统评估了33种癌症类型中全部八种HAUS基因的表达和预后意义。聚焦于肝肝细胞癌（liver hepatocellular carcinoma, LIHC），研究人员基于HAUS共表达模式识别了分子亚型，并表征了其与临床病理特征、肿瘤微环境（tumor microenvironment, TME）、免疫检查点表达及治疗反应的关联。采用LASSO（Least Absolute Shrinkage and Selection Operator）和Cox回归分析开发了HAUS相关预后特征，并在独立的ICGC（International Cancer Genome Consortium）队列中进行了验证。HAUS成员与关键特征基因之间的调控关系通过免疫组织化学（immunohistochemistry, IHC）和蛋白质印迹法（Western blotting）进行了实验验证。结果：大多数HAUS成员在多种癌症中过表达，并与不良临床结局相关。在LIHC中，两种基于HAUS的亚型在生存、临床病理特征、免疫特征、突变负荷、干性指数及预测治疗反应方面存在显著差异。包含DTYMK（deoxythymidylate kinase，脱氧胸苷酸激酶）和SPP1（secreted phosphoprotein 1，分泌性磷蛋白1）的预后特征能有效将LIHC患者分层为不同的风险组，高风险组显示出显著更差的生存。将HAUS相关风险评分与临床病理变量整合的列线图（nomogram）表现出强大的预测性能。重组骨桥蛋白（recombinant osteopontin）在LIHC细胞系中诱导HAUS1和DTYMK表达，支持了功能性的SPP1–HAUS1/DTYMK轴。结论：该研究确立了HAUS基因在癌症中的广泛致癌相关性，并证明基于HAUS的特征能够实现LIHC中的预后分层。

研究背景与问题：肝肝细胞癌（LIHC）是全球常见的恶性肿瘤之一，发病率和死亡率居高不下，尽管靶向治疗和免疫治疗取得了一定进展，但患者预后仍不理想，存在显著的异质性。HAUS家族蛋白作为微管成核和有丝分裂纺锤体形成的关键调控因子，近年来被发现与多种肿瘤进展相关，但HAUS的致癌功能在不同组织中的保守性或特异性尚不清楚，缺乏系统性的泛癌分析。此外，HAUS表达与肿瘤微环境、治疗反应及患者结局的关系，以及哪些HAUS成员最具预后生物标志物或治疗靶点潜力，仍有待阐明。因此，开展本研究以识别泛癌相关的HAUS驱动因子和情境依赖成员，为精准肿瘤学提供依据。

研究人员开展了什么研究：该研究利用TCGA数据库对33种癌症类型进行了全面的泛癌分析，系统评估了全部八种HAUS基因的表达模式和预后意义；随后聚焦于LIHC，基于HAUS共表达模式进行共识聚类分析识别分子亚型，深入分析了亚型与临床病理特征、肿瘤微环境、免疫检查点表达及治疗反应的关联；进一步采用LASSO回归和Cox回归分析筛选关键基因，构建了HAUS相关的预后特征并在独立ICGC队列中验证；最后通过免疫组织化学和蛋白质印迹实验验证了HAUS相关基因的表达调控关系。

研究结论与意义：研究确立了HAUS基因在多种癌症中的广泛致癌相关性，尤其是在LIHC中的预后价值。基于HAUS分子模式构建的DTYMK-SPP1双基因预后特征能够有效分层LIHC患者的生存风险，整合该风险评分的列线图具有良好的预测准确性。研究还发现SPP1可能通过正调控HAUS1和DTYMK参与LIHC进展，为理解HAUS家族的生物学功能和开发LIHC的精准治疗策略提供了新视角。该论文发表于《Journal of Hepatocellular Carcinoma》。

主要关键技术方法：本研究使用TCGA和ICGC-LIRI的LIHC队列以及GEO数据库（GSE149614）的10例原发LIHC肿瘤和8例非肿瘤组织的单细胞测序数据作为样本来源。采用的技术方法包括：基于RNA测序数据的差异表达分析和Kaplan-Meier生存分析；共识聚类分析识别HAUS分子亚型；单样本基因集富集分析（ssGSEA）评估免疫细胞浸润；免疫表型评分（IPS）和肿瘤免疫功能障碍与排除（TIDE）算法预测免疫检查点抑制剂反应；基因集变异分析（GSVA）进行功能富集分析；基于癌症基因组计划（CGP）细胞系药理基因组数据的Ridge回归模型预测药物敏感性（IC50）；单细胞RNA测序分析使用Seurat v3.0进行质量控制、降维、聚类和UMAP可视化；LASSO回归和多变量Cox回归构建预后特征；时间依赖性ROC曲线和AUC评估模型性能；免疫组织化学分析使用组织芯片（HLivH150CS03）和ImageJ-IHC Profiler进行定量；蛋白质印迹法检测重组骨桥蛋白处理后HAUS1和DTYMK蛋白表达变化。

研究结果：

泛癌HAUS家族基因表达谱分析：通过UCSC Xena和GTEx数据库对33种癌症类型的RNA测序数据进行分析，发现大多数HAUS成员在包括胆管癌（CHOL）、结肠腺癌（COAD）、多形性胶质母细胞瘤（GBM）、头颈部鳞状细胞癌（HNSC）、肾透明细胞癌（KIRC）、脑低级别胶质瘤（LGG）和LIHC在内的多种实体瘤中显著上调；而在急性髓系白血病（LAML）、卵巢浆液性囊腺癌（OV）和前列腺腺癌（PRAD）中主要表现为下调。

泛癌HAUS家族基因预后意义分析：Cox回归分析显示，在肾上腺皮质癌（ACC）、肾嫌色细胞癌（KICH）、LGG、LIHC、间皮瘤（MESO）和PRAD中，大多数HAUS家族基因是总生存期（OS）的危险因素；对于疾病特异性生存期（DSS），在ACC、KICH、LGG、LIHC、MESO和PRAD中为危险因素；对于无进展生存期（PFS），在ACC、KICH、LGG、LIHC和PRAD中亦为危险因素。相反，在直肠腺癌（READ）和胸腺瘤（THYM）中，HAUS家族基因大多表现出保护效应。整合分析表明，HAUS高表达与LGG和LIHC的生存不良显著相关。

LIHC中HAUS家族基因高表达与进展及不良预后的关联：TCGA和ICGC-LIRI数据验证显示，全部八种HAUS基因在LIHC组织中均显著上调。人蛋白质图谱数据库的免疫组织化学分析显示，HAUS2和HAUS4在LIHC组织中显著过表达。单变量Cox回归分析确定HAUS1、HAUS2、HAUS5、HAUS6和HAUS8为潜在预后危险因素，其高表达与较短的总生存期相关。HAUS1、HAUS2、HAUS6和HAUS8的表达与晚期病理分期显著相关，且随LIHC进展而明显升高。

基于HAUS家族基因共识聚类识别两种LIHC亚型：Pearson相关分析显示HAUS1、HAUS2、HAUS3、HAUS5、HAUS6和HAUS8存在显著正相关（r > 0.5）。基于这六个HAUS成员进行共识聚类，将370例LIHC患者分为HAUS高表达簇（C1）和HAUS低表达簇（C2）。主成分分析（PCA）确认了两种亚型的 distinct 分离。生存分析表明C1簇患者的中位总生存期显著短于C2簇。临床上，HAUS高表达簇与更具侵袭性的肿瘤特征密切相关，包括更高的组织学分级（G3: 43.23% vs. 25.84%; G4: 4.52% vs. 2.39%）、更晚的临床分期（III期: 31.33% vs. 18.97%）、更高的血清甲胎蛋白（AFP）水平（>400 ng/mL: 36.36% vs. 14.46%）、更年轻的诊断年龄以及更严重的纤维化程度（Ishak评分3–4: 44.55% vs. 22.82%）。

HAUS分子模式与LIHC患者TME的相关性：ssGSEA分析显示C1簇活化CD8⁺ T细胞浸润减少，总CD4⁺ T细胞和免疫抑制性Th2细胞比例增加。C1簇中多种免疫检查点基因显著上调，包括CD274（PD-L1）、PDCD1LG2（PD-L2）、PDCD1（PD-1）、CTLA4、LAG3、TIGIT、ICOS、LGALS9、CSF1R、IL10/IL10RB、VEGFA、TGFB1-3、TGFBR1以及TNFSF/TNFRSF家族成员。IPS和TIDE算法分析显示C1簇IPS评分显著低于C2簇，包括IPS-PD1/PD-L1/PD-L2阻断剂和IPS-CTLA4阻断剂评分，表明C1簇免疫原性降低，免疫检查点抑制剂获益可能性降低。TIDE评分在C1簇中显著升高，提示更高的免疫检查点阻断耐药风险。机制上，C1簇中T细胞排除增加，髓源性抑制细胞（MDSCs）和M2型肿瘤相关巨噬细胞富集，免疫检查点（IFNG、CD274）和Merck18特征上调，而T细胞功能障碍和CD8⁺ T细胞活性减弱。

HAUS分子模式的功能分析：GSVA-HALLMARK分析显示C1簇中细胞周期相关评分（G2M检查点、有丝分裂纺锤体、DNA修复、E2F靶点、MYC靶点）显著更高。GSVA-KEGG分析突出C1簇中DNA修复和细胞周期通路（同源末端连接、碱基切除修复、错配修复、核苷酸切除修复、DNA复制）以及PI3K/AKT/mTOR、mTORC1、Wnt/β-catenin、NOTCH和TGF-β信号通路的显著富集。突变谱分析显示C1簇TP53（40% vs. 15%）和OBSCN（12% vs. 4%）突变频率更高，CTNNB1（19% vs. 31%）和ALB（6% vs. 13%）突变频率更低。C1簇对Nutlin 3a（p53稳定剂）、roscovitine（Cdc2/CDK2/CDK5抑制剂）和PD0332991（CDK4/6抑制剂）的预测IC50值显著更高，表明对这些药物敏感性降低。

HAUS分子模式与干性特征及治疗反应的相关性：C1簇mRNA干性指数（mRNAsi）评分显著高于C2簇，癌症干细胞相关标志物EPCAM、ANPEP、CD24、CD44和CD47表达升高。基于CGP数据库的药物敏感性预测显示，C2簇对多种MEK抑制剂（PD0325901、RDEA119、AZD6244、CI-1040）、EGFR抑制剂（Erlotinib）、VEGFR抑制剂（AMG 706）以及Imatinib和AZD0530的预测IC50值显著更低，提示对这些药物更敏感；而C1簇对顺铂、吉西他滨、多柔比星和甲氨蝶呤的预测IC50值更低，表明对这些化疗药物更敏感。

HAUS家族基因可能是肝癌干细胞样细胞的新型标志物：对10例原发LIHC肿瘤和8例非肿瘤组织的单细胞测序分析（GSE149614）将细胞注释为T细胞、巨噬细胞、内皮细胞、肝癌干细胞样细胞（LCSC-like cells）、NK细胞、纤维母细胞、浆细胞、B细胞、肝细胞和肥大细胞。六个HAUS成员在几乎所有细胞群中广泛表达，其中HAUS1、HAUS5、HAUS6和HAUS8在LCSC样细胞中显示特异性高表达，其表达与LCSC样细胞的存在呈正相关。进一步将LCSC样细胞分为16个亚群，HAUS1、HAUS5、HAUS6和HAUS8在LCSC(6)、LCSC(7)、LCSC(8)、LCSC(15)和LCSC(16)中高表达。功能富集分析显示这些亚群显著富集于细胞周期、G2/M检查点、E2F靶点、DNA修复、氧化磷酸化、上皮-间质转化（EMT）、ROS通路、MTORC1信号、MYC靶点、糖酵解、缺氧、p53通路和PI3K/AKT/MTOR信号等通路。

HAUS分子模式相关预后特征的构建：通过limma包识别HAUS高、低表达簇之间的差异表达基因，获得6440个上调和362个下调基因。经LASSO惩罚Cox回归分析，在最佳λ截断值处筛选出2基因预后特征（DTYMK和SPP1）。两基因在TCGA-LIHC和ICGC-LIRI队列的LIHC组织中均显著上调。风险评分计算公式为：风险评分 = (0.55 × DTYMK mRNA表达水平) + (0.11 × SPP1 mRNA表达水平)。基于中位风险评分将患者分为高风险和低风险组，高风险组死亡率更高、生存时间更短。Kaplan-Meier分析显示高风险组总生存期显著降低。时间依赖性ROC曲线显示该模型具有良好的敏感性和特异性，TCGA队列中1年、3年和5年生存率的AUC值分别为0.755、0.677和0.698。在ICGC-LIRI独立验证队列中，高风险组同样显示出显著更差的生存结局，1年和3年OS的AUC分别为0.697和0.753。单变量和多变量Cox回归分析证实风险评分是LIHC总生存期的独立预后因素。整合病理分期、T分期、M分期和风险评分的列线图，校准曲线与45度对角线接近，表明预测生存概率与实际观察结果具有良好的一致性。

SPP1正调控LIHC中HAUS1和DTYMK表达：对69例LIHC患者组织芯片的免疫组织化学分析显示，HAUS1主要定位于核周区域，在肿瘤组织中表达上调；DTYMK在细胞质和核中均有检测，肿瘤细胞核中表达显著更高。SPP1与HAUS1、SPP1与细胞质/核DTYMK、HAUS1与细胞质/核DTYMK之间均存在正相关。在HepG2和Huh7细胞中，重组骨桥蛋白（200 ng/mL）处理显著诱导HAUS1和DTYMK蛋白表达，提示SPP1可能作为HAUS1和DTYMK的正调控因子参与LIHC进展。

讨论部分总结：HAUS复合物是进化上保守的微管成核调控因子，该研究通过泛癌分析揭示了HAUS成员在多种癌症中的广泛致癌相关性，特别是在LIHC中的预后价值。HAUS高表达与LIHC的免疫逃逸、干细胞样特性及治疗耐药密切相关。研究开发的DTYMK-SPP1双基因预后特征具有独立的预后预测能力，为LIHC患者分层提供了实用工具。SPP1可能通过正调控HAUS1和DTYMK形成功能性调控轴，为理解LIHC进展机制和开发治疗靶点提供了新线索。然而，研究仍存在局限性：生物信息学分析依赖于回顾性TCGA数据，可能受治疗异质性、样本质量和人口统计学失衡的影响；多数发现为关联性而非因果性；尽管免疫组织化学和蛋白质印迹提供了实验支持，但仍需更多体外和体内研究阐明SPP1、HAUS和DTYMK之间的机制联系；在独立队列和前瞻性临床研究验证之前，其临床转化价值尚待确立。

研究结论：该研究表征了LIHC中HAUS基因的表达模式，并开发了HAUS相关预后风险模型以评估其临床相关性。ROC和生存分析表明，HAUS表达特征及相应的风险评分在LIHC中具有预后价值。生物信息学分析进一步显示HAUS表达与免疫检查点表达、T细胞耗竭标志物及肿瘤微环境中免疫抑制细胞群相关。免疫组织化学染色和蛋白质印迹分析为转录组学分析识别的差异表达模式提供了初步实验支持。然而，该HAUS相关特征的临床价值需要在更大规模的独立队列中进行验证后方可考虑转化应用。此外，由于目前结果主要为关联性，需要功能研究来阐明HAUS失调是否直接促进LIHC进展。未来工作应集中于阐明HAUS与免疫逃逸相关的机制，并评估其作为LIHC治疗靶点的潜力。