葡萄膜黑色素瘤(Uveal Melanoma, UVM)是成人中最常見的原發(fā)性眼內(nèi)惡性腫瘤，起源于虹膜、睫狀體和脈絡(luò)膜的黑色素細胞，年發(fā)病率約為每百萬人5-7例，在白人、中老年群體中尤為高發(fā)。該病具有極強的肝轉(zhuǎn)移傾向，約50%的患者在確診后10年內(nèi)發(fā)生轉(zhuǎn)移，其中約95%為單純肝轉(zhuǎn)移，導(dǎo)致患者一旦發(fā)生轉(zhuǎn)移，中位生存期僅6-12個月。UVM的治療面臨多重困境：常規(guī)化療療效有限；因其腫瘤突變負荷(Tumor Mutational Burden, TMB)低、缺乏新抗原、腫瘤微環(huán)境呈“冷”免疫狀態(tài)以及高頻發(fā)生BAP1功能缺失突變等原因，導(dǎo)致其對免疫檢查點抑制劑(Immune Checkpoint Inhibitors, ICIs)如抗PD-1/抗CTLA-4抗體的反應(yīng)極差，遠遜于皮膚黑色素瘤。同時，疾病的罕見性也嚴重限制了臨床試驗的開展和藥物研發(fā)投入。因此，尋找能夠早期預(yù)測轉(zhuǎn)移、有效評估預(yù)后并指導(dǎo)治療的分子標(biāo)志物成為UVM臨床管理的迫切需求。

為深入探究UVM的分子特征，研究者設(shè)計了一套整合分析方法。首先，利用GEO數(shù)據(jù)庫中的單細胞RNA測序(Single-Cell RNA Sequencing, scRNA-seq)數(shù)據(jù)集GSE139829，納入了來自原發(fā)和轉(zhuǎn)移性UVM樣本的50,933個高質(zhì)量細胞。經(jīng)過嚴格的質(zhì)量控制，使用Harmony算法校正批次效應(yīng)，并通過主成分分析(Principal Component Analysis, PCA)、均勻流形逼近與投影(Uniform Manifold Approximation and Projection, UMAP)進行降維可視化。基于經(jīng)典標(biāo)記基因的表達譜，研究者最終將細胞注釋為7種主要類型：惡性細胞（標(biāo)記基因：MITF, MLANA）、T/NK細胞（CD3E, CD3D）、單核/巨噬細胞（CD14, CD68, CD163, C1QA）、B細胞/漿細胞（MS4A1, CD38, SDC1）、內(nèi)皮細胞（PECAM1, CD34, CLDN5）、成纖維細胞（DCN, C1R, C1S）以及光感受器細胞（RCVRN）。

在單細胞層面，研究者提取了惡性細胞并進行差異表達分析，篩選出原發(fā)與轉(zhuǎn)移性惡性細胞之間的差異表達基因(Differentially Expressed Genes, DEGs)。功能富集分析（KEGG和GO）顯示，轉(zhuǎn)移性細胞在免疫調(diào)節(jié)、代謝重編程和蛋白質(zhì)穩(wěn)態(tài)相關(guān)通路中顯著富集，例如白細胞跨內(nèi)皮遷移、蛋白酶體通路、抗原加工與呈遞、谷胱甘肽代謝，以及線粒體內(nèi)膜、氧化磷酸化、ATP合成偶聯(lián)的電子傳遞等過程。

為了從海量數(shù)據(jù)中提煉出關(guān)鍵的預(yù)后信息，研究者采用了整合策略。首先，基于癌癥基因組圖譜(The Cancer Genome Atlas, TCGA)的UVM隊列，通過單因素Cox回歸分析，篩選出900個與總生存期(Overall Survival, OS)顯著相關(guān)的預(yù)后基因。接著，將這900個基因與單細胞分析中得到的惡性細胞DEGs取交集，獲得候選基因集。為了避免過擬合并篩選出最穩(wěn)健的預(yù)測因子，研究者應(yīng)用了最小絕對收縮和選擇算子(Least Absolute Shrinkage and Selection Operator, LASSO)回歸分析。通過10折交叉驗證選擇最優(yōu)正則化參數(shù)λ，最終鎖定了11個UVM特異性預(yù)后基因，構(gòu)建了一個多基因風(fēng)險預(yù)測模型。

這11個關(guān)鍵基因包括：絲氨酸蛋白酶抑制劑B家族成員9(SERPINB9)、聚腺苷二磷酸核糖聚合酶家族成員8(PARP8)、程序性細胞死亡4(PDCD4)、癲癇相關(guān)6同源物樣2(SEZ6L2)、KDEL內(nèi)質(zhì)網(wǎng)蛋白滯留受體3(KDELR3)、脂滴包被蛋白2(PLIN2)、干擾素刺激外切核酸酶20(ISG20)、S100鈣結(jié)合蛋白A13(S100A13)、溶質(zhì)載體家族45成員2(SLC45A2)、γ-氨基丁酸A型受體相關(guān)蛋白樣1(GABARAPL1)和組織蛋白酶F(CTSF)。

研究者為TCGA隊列中的每位患者計算了風(fēng)險評分，并以中位風(fēng)險評分為界，將患者分為高風(fēng)險組和低風(fēng)險組。模型的預(yù)測性能在多個獨立數(shù)據(jù)集中得到了卓越驗證。在TCGA訓(xùn)練集中，該模型對1年、3年和5年生存率的曲線下面積(Area Under the Curve, AUC)分別高達0.897、0.973和1.000，Kaplan-Meier生存分析顯示高風(fēng)險組患者的生存結(jié)局顯著更差。模型在外部驗證集GSE22138和GSE84976中同樣表現(xiàn)出強大的預(yù)測能力和穩(wěn)定性，證實了其良好的泛化性。

為了評估風(fēng)險模型是否是不依賴于其他臨床因素的獨立預(yù)后因子，研究者進行了單因素和多因素Cox回歸分析。結(jié)果顯示，在多因素模型中，排除了年齡等混雜因素后，風(fēng)險評分仍然是UVM患者OS的獨立預(yù)后因素，其風(fēng)險比(Hazard Ratio, HR)高達6.983。為了將這一模型更直觀地應(yīng)用于臨床，研究者構(gòu)建了一個整合了年齡和風(fēng)險評分的列線圖(Nomogram)，可以便捷地預(yù)測個體患者在1年、3年和5年的生存概率。校準(zhǔn)曲線顯示，該列線圖的預(yù)測結(jié)果與實際觀察到的生存率高度一致，展現(xiàn)了其良好的臨床實用性。

該研究進一步揭示了風(fēng)險模型背后的生物學(xué)內(nèi)涵，特別是腫瘤免疫微環(huán)境的差異。通過CIBERSORT算法對腫瘤浸潤免疫細胞進行分析，發(fā)現(xiàn)高風(fēng)險組患者的腫瘤微環(huán)境中，活化的CD4⁺T細胞、CD8⁺T細胞、自然殺傷(Natural Killer, NK)細胞、輔助性T細胞1(T helper 1, Th1)、調(diào)節(jié)性T細胞(Regulatory T cells, Tregs)、髓源性抑制細胞(Myeloid-Derived Suppressor Cells, MDSCs)以及多種抗原呈遞細胞的浸潤水平顯著更高。這種看似“活躍”的免疫狀態(tài)，實際上可能伴隨著更強的免疫抑制和調(diào)節(jié)特性。

更關(guān)鍵的是，免疫檢查點基因的表達譜分析顯示，高風(fēng)險組中多個關(guān)鍵的免疫抑制分子，如CD276(B7-H3)、HAVCR2(TIM-3)、LAG3、PDCD1LG2(PD-L2)和TIGIT的表達均顯著上調(diào)。這強烈提示，高風(fēng)險UVM患者腫瘤細胞的免疫逃逸潛力更強，可能對現(xiàn)有的免疫檢查點阻斷治療產(chǎn)生不同的反應(yīng)，為個性化免疫治療策略的選擇提供了線索。

研究還利用pRRophetic R包，基于基因表達譜預(yù)測了患者對常見化療藥物的敏感性，即半數(shù)抑制濃度(Half Maximal Inhibitory Concentration, IC₅₀)。分析發(fā)現(xiàn)，高風(fēng)險組患者對阿糖胞苷(cytarabine)和甲氨蝶呤(methotrexate)的預(yù)測IC₅₀值顯著更低，意味著他們可能對這些藥物更敏感。相反，高風(fēng)險組對絲裂霉素C(mitomycin C)的預(yù)測IC₅₀值顯著更高，表明其敏感性可能較低。這些發(fā)現(xiàn)提示，該風(fēng)險模型不僅能預(yù)測預(yù)后，還可能為指導(dǎo)不同風(fēng)險亞組患者的化療藥物選擇提供有價值的信息。

研究分別在單細胞水平和群體轉(zhuǎn)錄組水平驗證了這11個預(yù)后基因的表達模式。單細胞數(shù)據(jù)顯示，這些基因在腫瘤微環(huán)境中呈現(xiàn)細胞類型特異性的分布。例如，SERPINB9主要在惡性細胞中高表達，PLIN2富集于巨噬細胞，而S100A13和ISG20則在內(nèi)皮細胞和B/漿細胞中表達較高。在群體水平，研究者分析了部分基因表達與臨床分期的關(guān)系，發(fā)現(xiàn)SERPINB9的表達在III期腫瘤中升高，在IV期又有所下降，顯示出與疾病進展的潛在關(guān)聯(lián)；而PLIN2和S100A13在不同分期間的表達差異未達到統(tǒng)計學(xué)顯著性。

在討論部分，文章深入探討了部分關(guān)鍵基因的已知功能及其在UVM中的潛在意義。例如，SERPINB9作為顆粒酶B的抑制劑，可通過保護腫瘤細胞免受細胞毒性T淋巴細胞(Cytotoxic T Lymphocytes, CTLs)和NK細胞的殺傷來促進免疫逃逸，是潛在的免疫治療靶點。SLC45A2是黑色素細胞相關(guān)抗原，與黑色素合成相關(guān)，已有研究探索以其為靶點的過繼性T細胞療法治療轉(zhuǎn)移性UVM。PLIN2參與脂滴代謝，與腫瘤的免疫抑制和代謝重編程微環(huán)境相關(guān)。S100A13作為S100蛋白家族成員，參與炎癥和細胞遷移等過程，其異常表達可能通過免疫相關(guān)機制影響疾病進程。

綜上所述，本研究通過整合單細胞與群體轉(zhuǎn)錄組數(shù)據(jù)，系統(tǒng)描繪了UVM原發(fā)與轉(zhuǎn)移狀態(tài)的分子圖譜及免疫微環(huán)境特征，并成功構(gòu)建并驗證了一個包含11個基因的穩(wěn)健預(yù)后風(fēng)險模型。該模型不僅能有效區(qū)分UVM患者的預(yù)后，還揭示了高風(fēng)險患者所具有的免疫抑制微環(huán)境、特定的免疫檢查點表達譜以及差異化的化療藥物敏感性。研究特別指出，PLIN2、S100A13和SERPINB9等基因的表達變化與疾病進展階段相關(guān)，是具有潛力的早期預(yù)測生物標(biāo)志物和治療靶點。這項工作為深入理解UVM的轉(zhuǎn)移生物學(xué)、開發(fā)新的預(yù)后工具以及制定更精準(zhǔn)的個體化治療策略奠定了重要的理論基礎(chǔ)。當(dāng)然，研究的結(jié)論主要基于生物信息學(xué)分析，未來需要通過體外實驗、動物模型以及前瞻性臨床研究進一步驗證這些基因的生物學(xué)功能及其臨床轉(zhuǎn)化價值。