[摘" "要]" "目的：尋找骨肉瘤抗腫瘤免疫及預后的新的生物標志物及其作用機制。方法：采用生物信息學方法對GEO數(shù)據(jù)庫和Target的骨肉瘤數(shù)據(jù)進行差異分析、功能富集、免疫細胞浸潤分析、生存分析和風險模型構(gòu)建等。常規(guī)H-E和免疫組織化學染色分析冠蛋白1A（coronin 1A， CORO1A）在骨肉瘤組織中的表達情況，并轉(zhuǎn)染si-CORO1A和HA-CORO1A，驗證CORO1A低表達和高表達對程序性細胞死亡-配體1（programmed cell death ligand 1， PD-L1）的影響。結(jié)果：CORO1A基因在骨肉瘤患者中高表達，CORO1A基因的表達與骨肉瘤中抗腫瘤免疫細胞的比例呈正相關(guān)，并預示著較好的生存與預后。CORO1A可促進PD-L1蛋白表達，兩者呈正相關(guān)。結(jié)論：CORO1A基因促進抗腫瘤免疫以利于生存，是骨肉瘤患者獨立預后的生物標志物。

[關(guān)鍵詞]" "骨肉瘤；冠蛋白1A；生物標志物；預后；免疫細胞

[中圖分類號]" "R738.1" " " " " " " "[文獻標志碼]" "A" " " " " " " "[文章編號]" "1674-7887（2025）02-0120-07

Role of CORO1A in anti-tumor immunity and prognosis in osteosarcoma*

TANG Tao**， CHEN Xi， DA Wenyue， WANG Hui***" " " " （Department of Pathology， First People？蒺s Hospital of Changzhou， Jiangsu 213003）

[Abstract]" "Objective： To find new biomarkers for anti-tumor immunity and prognosis of osteosarcoma and mechanisms. Methods： Bioinformatics methods were used to perform differential analysis， functional enrichment， immune cell infiltration analysis， survival analysis and risk model construction on osteosarcoma data from GEO and Target， and H-E and immun-ohistochemical staining were performed to analyse coronin 1A（CORO1A） osteosarcoma tissues. Si-CORO1A and HA-CORO1A were transfected to verify the effect of low and high expression of CORO1A on programmed cell death ligand 1（PD-L1）. Results： The CORO1A gene was highly expressed in osteosarcoma patients. The expression of the CORO1A gene was positively correlated with the proportion of anti-tumor immune cells in osteosarcoma and predicted better survival and prognosis. CORO1A can promote the expression of PD-L1 protein， and the two are positively correlated. Conclusion： The CORO1A gene promotes antitumor immunity to facilitate survival and is an independent prognostic biomarker for osteosarcoma patients.

[Key words]" "osteosarcoma; coronin 1A; biomarker; prognosis; immune cell

骨肉瘤是骨的原發(fā)惡性腫瘤，其特征是惡性腫瘤細胞產(chǎn)生類骨質(zhì)或不成熟骨，發(fā)病率為每年每百萬人2～3例，多發(fā)于青少年[1-2]。近年來隨著醫(yī)療技術(shù)的發(fā)展，骨肉瘤患者的治療效果得到了一定的提高[3-5]。目前的標準治療是以包含多柔比星、順鉑、大劑量甲氨蝶呤（MAP方案）為基礎(chǔ)組成的新輔助化療-手術(shù)-輔助化療綜合治療模式。大多數(shù)骨肉瘤患者明確診斷時已存在亞臨床的轉(zhuǎn)移性疾病[6-8]，患者預后無顯著改善。因此，探尋骨肉瘤發(fā)生發(fā)展的分子機制，尋找評估預后新的生物標志物對患者的治療至關(guān)重要。有證據(jù)[9]表明，除了基因突變外，腫瘤微環(huán)境，特別是免疫細胞，也參與了骨肉瘤的發(fā)展。

本課題組前期對公共數(shù)據(jù)進行進一步的挖掘，采用生物信息學的方法對GEO數(shù)據(jù)庫和Target的骨肉瘤數(shù)據(jù)進行差異分析，找到了眾多影響骨肉瘤患者預后的差異基因[10]，如冠蛋白1A（coronin 1A， CORO1A）、骨髓細胞瘤癌基因（myelocytomatosis oncogene， MYC）、過氧化氫酶（catalase， CAT）等。CORO1A基因與細胞的凋亡和信號轉(zhuǎn)導等有關(guān)，與骨肉瘤相關(guān)的研究較為少見，為進一步探討CORO1A基因在骨肉瘤中發(fā)揮的作用。本研究采用生物信息學技術(shù)分析CORO1A基因與腫瘤免疫微環(huán)境中腫瘤浸潤淋巴細胞的關(guān)系，并通過免疫組織化學染色（immunohistochemistry， IHC）和蛋白印跡探討該基因在骨肉瘤中的表達情況及其與程序性細胞死亡-配體1（programmed cell death ligand 1， PD-L1）的關(guān)系，旨在為骨肉瘤治療提供新的策略。

1" "材料與方法

1.1" "樣本來源" "GEO數(shù)據(jù)庫下載GSE42352數(shù)據(jù)集（測試集）：正常樣本數(shù)15例，腫瘤樣本數(shù)103例[11-12]。從Target網(wǎng)站下載Target-OS數(shù)據(jù)集（驗證集）：腫瘤樣本數(shù)101例，包含臨床信息患者數(shù)96例。臨床樣本來源于2021年5月—2024年5月在常州市第一人民醫(yī)院手術(shù)的骨肉瘤患者，經(jīng)醫(yī)院倫理委員會批準[批文號：（2021）科第069號]，并獲得每例患者的知情同意。U2OS骨肉瘤細胞株（貨號：HTX1634）購自深圳華拓細胞庫。

1.2" "H-E染色和IHC" "標本經(jīng)4%多聚甲醛固定、脫水、透明、浸蠟，石蠟包埋、常規(guī)制片，H-E染色。IHC采用DAKO Omins全自動免疫組化儀，一抗PD-L1（克隆號：22C3）購自Dako公司，一抗CORO1A（克隆號：OTI1A5）、CD8A（克隆號：SP16）和CD56（克隆號：UMAB83）購自北京中杉生物技術(shù)開發(fā)有限公司。操作步驟嚴格按照產(chǎn)品說明書進行。根據(jù)顯微鏡下200倍的染色程度和染色強度評估蛋白質(zhì)表達水平。染色程度評分為0～4分：0分，0%；1分，≤10%；2分，gt;10%～50%；3分，gt;50%～80%；4分，gt;80%～100%。染色強度分為弱染色、中染色、強染色3個等級，相應(yīng)的強度分數(shù)范圍為1～3分。最終的IHC分數(shù)為染色程度評分乘以染色強度分數(shù)。

1.3" "免疫細胞浸潤分析" "通過R語言的CIBERSORT[13]和TISIDB[14]計算不同CORO1A mRNA表達狀態(tài)的骨肉瘤組織中腫瘤浸潤淋巴細胞（tumor-infiltrating lymphocyte， TIL）的相對豐度。

1.4" "差異表達分析" "通過R語言的limma包[15]進行差異分析，根據(jù)組間差異變化倍數(shù)（fold change）及差異顯著性（P value）來篩選差異（共表達）基因，設(shè)定|logFC|＞1和adj.P＜0.05為閾值。

1.5" "GO和KEGG分析" "采用clusterProfiler包[16]進行差異基因的GO和KEGG分析，并使用ggplot2包進行可視化。GO分析包括生物學過程（biological process， BP）、分子功能（molecular function， MF）、細胞組分（cellular component， CC）。

1.6" "ROC曲線評估診斷性能" "通過pROC包[17]繪制ROC曲線，分析15例正常樣本和103例骨肉瘤樣本中CORO1A表達量的差異，并確定產(chǎn)生最高似然比的最佳截斷值以判斷CORO1A對骨肉瘤的識別閾值。并使用survivalROC包[18]驗證風險模型預測骨肉瘤生存預后的準確性。AUCgt;0.5則能夠直觀反映生物標志物對疾病的診斷價值。同時，AUC越接近1，離（0， 1）點越近，證明診斷試驗的真實性越好。

1.7" "生存分析" "使用survival包進行Target-OS數(shù)據(jù)集生存分析，包括總生存時間（overall survival， OS）、無事件（疾?。┥鏁r間（disease free survival， DFS）、疾病進展時間（time to progress， TTP）。并通過驗證風險模型對骨肉瘤患者的總生存預后進行分析。

1.8" "siRNA干擾實驗和HA-pCDNA實驗" "傳代骨肉瘤細胞U2OS至6孔板，待生長至70%左右時，加抗生素無血清培養(yǎng)基饑餓細胞8 h；Lipo3000分別轉(zhuǎn)染CORO1A-siRNA（正義鏈5？蒺-CCAUGACAGUGCC-UAGAAAdTdT-3？蒺，反義鏈5？蒺-UUUCUAGGCACUGUCAUGGdTdT-3？蒺）和HA-CORO1A（正義鏈5？蒺-ATGA-GCCGGCAGGTGGTC-3？蒺，反義鏈5？蒺-CTACTTGGCC-TGGACTGTCTC-3？蒺）至U2OS細胞，6 h后換成完全培養(yǎng)基繼續(xù)培養(yǎng)12～16 h，空載體組作為空白對照。實驗重復5次。

1.9" "統(tǒng)計學方法" "采用Excel（Microsoft）和R統(tǒng)計軟件（版本4.0.2）進行數(shù)據(jù)分析。CORO1A和TIL之間的相關(guān)性采用Spearman相關(guān)分析。Plt;0.05為差異具有統(tǒng)計學意義。

2" "結(jié)" " " 果

2.1" "CORO1A在骨肉瘤中的異常表達及預后價值" "首先，評估分析TIMER中所有TCGA癌癥組織中CORO1A的分布，發(fā)現(xiàn)在大多數(shù)腫瘤中CORO1A表達顯著上調(diào)，包括肉瘤組織（sarcoma， SARC），見圖1A。為進一步驗證CORO1A基因在骨肉瘤中的表達情況，將TCGA的骨肉瘤數(shù)據(jù)進行統(tǒng)計分析，結(jié)果顯示CORO1A基因在骨肉瘤中高表達，與正常組織相比差異有統(tǒng)計學意義（Plt;0.01），見圖1B。通過生存分析發(fā)現(xiàn)，CORO1A基因高表達患者的DFS和OS較低表達患者明顯延長（Plt;0.01），提示CORO1A基因高表達是腫瘤患者較好的預后指標，見圖1C～D。

2.2" "CORO1A基因在骨肉瘤中的診斷價值" "為分析CORO1A基因在骨肉瘤中的診斷價值，用TCGA的數(shù)據(jù)構(gòu)建ROC曲線和nomogram預測模型，結(jié)果顯示CORO1A基因在骨肉瘤中有較高的診斷價值，且骨肉瘤中CORO1A基因的表達情況可以預測患者1、3、5年的生存率，見圖2。

2.3" "CORO1A與腫瘤免疫微環(huán)境的關(guān)聯(lián)" "為了解CORO1A基因在骨肉瘤預后中的潛在機制，采用TISIDB分析所有TCGA腫瘤中哪些類型的TIL可能受到CORO1A的調(diào)節(jié)，以此來探討CORO1A與腫瘤免疫微環(huán)境的關(guān)聯(lián)，見圖3A。在肉瘤中CORO1A與B細胞（r=0.742， Plt;0.001）、激活的CD8+T細胞（active CD8， Act CD8）（r=0.855， Plt;0.001）和自然殺傷細胞（natural killer cell， NK細胞）（r=0.684， Plt;0.001）呈正相關(guān)，見圖3B～D。同時發(fā)現(xiàn)，CORO1A與免疫相關(guān)的C-C基序趨化因子配體（C-C motif chemokine ligand， CCL）2（r=0.660， Plt;0.001）、CCL3（r=0.730， Plt;0.001）和CCL4（r=0.817， Plt;0.001）呈正相關(guān)，見圖3E～G。此外，還采用CIBERSORT驗證CORO1A mRNA差異表達對骨肉瘤患者腫瘤免疫微環(huán)境中狀態(tài)的影響。CORO1A mRNA高表達組（紅色為高表達組，藍色為低表達組）的抗腫瘤免疫細胞比例較高，如CD8+T細胞（CD8T）（Plt;0.05）、記憶B細胞（Plt;0.05），而促腫瘤免疫細胞比例較低，如單核細胞（Mono）（Plt;0.05）、M2巨噬細胞（M2）（Plt;0.05），見圖3H。結(jié)果表明，CORO1A與骨肉瘤的抗腫瘤免疫相關(guān)，進一步驗證CORO1A與骨肉瘤患者較好預后的相關(guān)性。并且提示通過預測這些免疫細胞的靶點小分子藥物，通過改善免疫細胞分布比例，進而達到治療骨肉瘤的目的。

2.4" "CORO1A不同表達組的差異基因分析和富集的生物學過程" "以CORO1A在測試集的103例腫瘤樣本中表達中位值進行分組，其中低表達組51例，高表達組52例；同樣對驗證集Target數(shù)據(jù)庫的101例骨肉瘤患者以CORO1A表達中位值分組，高表達組50例和低表達組51例。使用limma包進行差異分析，設(shè)定|logFC|＞1和adj.P＜0.05為閾值，對兩組差異基因取交集，進行GO和KEGG富集分析。發(fā)現(xiàn)富集到的生物學過程有抗原加工提呈、免疫效應(yīng)過程的調(diào)控及MHC Ⅱ類分子的抗原遞呈等，見圖4，進一步驗證了CORO1A有助于骨肉瘤中抗腫瘤免疫微環(huán)境。

2.5" "骨肉瘤組織中CORO1A高表達與部分腫瘤浸潤淋巴細胞的關(guān)聯(lián)分析" "為進一步驗證骨肉瘤組織中CORO1A高表達對腫瘤浸潤淋巴細胞的影響，對骨肉瘤患者的組織進行IHC。結(jié)果顯示，具有高表達豐度的CORO1A的腫瘤組織中，與抗腫瘤免疫相關(guān)的免疫細胞CD8+T細胞（CD8A標記）和CD56+NK細胞也呈現(xiàn)出高表達豐度，見圖5A。此外，高表達豐度的CORO1A腫瘤組織中，PD-L1有呈現(xiàn)出高表達狀態(tài)。相關(guān)性分析顯示，骨肉瘤組織中CORO1A的表達與PD-L1、CD8A和CD56均呈正相關(guān)（r=0.95、0.74、0.62，均Plt;0.001），見圖5B。結(jié)果顯示，在骨肉瘤組織中，CORO1A高表達組中PD-L1、CD8A和CD56的表達明顯高于CORO1A低表達組（Plt;0.001），見圖5C。

2.6" "CORO1A基因?qū)D-L1蛋白表達的影響" "為進一步驗證CORO1A與PD-L1兩者的關(guān)系，采用體外細胞實驗，分別檢測CORO1A不同表達狀態(tài)下，PD-L1蛋白表達情況。結(jié)果顯示：轉(zhuǎn)染si-CORO1A抑制CORO1A表達后，PD-L1蛋白表達也隨之降低，CORO1A高表達后，PD-L1蛋白表達也升高（Plt;0.001），見圖6。

3" "討" " " 論

CORO1A與細胞信號傳導、細胞黏附和細胞遷移等生物學過程密切相關(guān)。X.F.FEI等[19]研究發(fā)現(xiàn)CORO1A可通過擬素化修飾直接結(jié)合RAS癌基因家族成員7（Ras-related GTP binding protein 7， Rab7）促進多泡體募集和Rab7的激活，隨后多泡體以依賴Rab7的方式降解，導致細胞外囊泡分泌減少，從而在腫瘤發(fā)展中扮演著重要角色。近年來，研究[20]發(fā)現(xiàn)CORO1A在一些癌癥中起到監(jiān)控疾病進展的作用，但在骨肉瘤中的作用鮮有報道。本研究通過對公共數(shù)據(jù)庫的骨肉瘤數(shù)據(jù)進行生物信息學分析發(fā)現(xiàn)，CORO1A在骨肉瘤患者中高表達，且高表達的COR-O1A在骨肉瘤中具有獨立的預后價值，并通過生物信息學分析證實CORO1A的高表達在骨肉瘤的發(fā)展中起到防御作用。

腫瘤浸潤免疫細胞與腫瘤發(fā)生、血管生成以及腫瘤細胞的生長和轉(zhuǎn)移密切相關(guān)，并能調(diào)控免疫細胞的數(shù)量和分化[21]。有證據(jù)[22]表明，腫瘤進展可能是癌細胞逃避宿主免疫監(jiān)視的結(jié)果。明確腫瘤微環(huán)境中浸潤的免疫細胞可能有助于闡明CORO1A在骨肉瘤中的作用機制。本研究發(fā)現(xiàn)CORO1A在骨肉瘤中高表達，CORO1A高表達的患者，其DFS和OS明顯延長，且其抗腫瘤免疫細胞比例較高，而促腫瘤免疫細胞比例較低。此外，還發(fā)現(xiàn)CORO1A的差異表達基因富集到免疫相關(guān)信號通路網(wǎng)絡(luò)，且與CD8A、CD56以及免疫相關(guān)趨化因子CCL2、CCL3和CCL4呈正相關(guān)，這些基因都與NK細胞和T細胞募集有關(guān)，在抑制腫瘤生長和改善預后方面發(fā)揮重要作用[23]。結(jié)果提示CORO1A是骨肉瘤獨立的預后指標，CORO1A是通過介導免疫細胞的浸潤來促進抗腫瘤免疫。

有趣的是，本研究發(fā)現(xiàn)GEO和Target中CORO1A的表達還與PD-L1呈正相關(guān)，抑制CORO1A基因PD-L1蛋白表達明顯被降低。浙江大學蔡志堅團隊[24]研究表明敲除CORO1A可抑制內(nèi)源性細胞外囊泡的產(chǎn)生，逆轉(zhuǎn)抗PD-L1抗體治療耐藥。多項研究[25-26]表明，腫瘤內(nèi)大量CD8+T細胞浸潤能夠在某種程度上緩解對免疫檢查點抑制劑治療的耐藥性。因此驗證CORO1A是否可以預測骨肉瘤對免疫檢查點抑制劑的反應(yīng)是很有希望。本研究結(jié)果可為患者免疫治療，特別是骨肉瘤患者提供新的見解，幫助臨床醫(yī)師為骨肉瘤患者選擇適當?shù)拇胧┎⒏纳祁A后。

[參考文獻]

[1]" "CHEN C L， XIE L， REN T T， et al. Immunotherapy for osteosarcoma： Fundamental mechanism， rationale， and recent breakthroughs[J]. Cancer Lett， 2021， 500：1-10.

[2]" "SAKAMOTO A， IWAMOTO Y. Current status and perspectives regarding the treatment of osteo-sarcoma： chemotherapy[J]. Rev Recent Clin Trials， 2008， 3（3）：228-231.

[3]" "THANINDRATARN P， DEAN D C， NELSON S D， et al. Advances in immune checkpoint inhibitors for bone sarcoma therapy[J]. J Bone Oncol， 2019， 15：100221.

[4]" "MASON N J， GNANANDARAJAH J S， ENGILES J B， et al. Immunotherapy with a HER2-targeting Listeria induces HER2-specific immunity and demonstrates potential therapeutic effects in a phase I trial in canine osteosarcoma[J]. Clin Cancer Res， 2016， 22（17）：4380-4390.

[5]" "YOSHIDA K， OKAMOTO M， SASAKI J， et al. Anti-PD-1 antibody decreases tumour-infiltrating regulatory T cells[J]. BMC Cancer， 2020， 20（1）：25.

[6]" "BERHE S， DANZER E， MEYERS P， et al. Unusual abdominal metastases in osteosarcoma[J]. J Pediatr Surg Case Rep， 2018， 28：13-16.

[7]" "WEDEKIND M F， WAGNER L M， CRIPE T P. Immunoth-erapy for osteosarcoma： where do we go from here？[J]. Pediatr Blood Cancer， 2018， 65（9）：e27227.

[8]" "XU J， XIE L， GUO W. PDGF/PDGFR effects in osteosarcoma and the “add-on” strategy[J]. Clin Sarcoma Res， 2018， 8：15.

[9]" "HE X Y， QUE L Y， YANG F， et al. Single-cell transcriptional profiling in osteosarcoma and the effect of neoadjuvant chemotherapy on the tumor microenvironment[J]. J Bone Oncol， 2024， 46：100604.

[10]" "LV Y G， WU L Y， JIAN H， et al. Identification and characterization of aging/senescence-induced genes in ost-eosarcoma and predicting clinical prognosis[J]. Front Imm-unol， 2022， 13：997765.

[11]" "KUIJJER M L， PETERSE E F， VAN DEN AKKER B E， et al. IR/IGF1R signaling as potential target for treatment of high-grade osteosarcoma[J]. BMC Cancer， 2013， 13：245.

[12]" "KUIJJER M L， VAN DEN AKKER B E W M， HILHORST R， et al. Kinome and mRNA expression profiling of high-grade osteosarcoma cell lines implies Akt signaling as possible target for therapy[J]. BMC Med Genomics， 2014， 7：4.

[13]" "SZKLARCZYK D， GABLE A L， LYON D， et al. STRING v11： protein-protein association networks with increased coverage， supporting functional discovery in genome-wide experimental datasets[J]. Nucleic Acids Res， 2019， 47（D1）：D607-D613.

[14]" "RU B， WONG C N， TONG Y， et al. TISIDB： an integrated repository portal for tumor-immune system interactions[J]. Bioinformatics， 2019， 35（20）：4200-4202.

[15]" "RITCHIE M E， PHIPSON B， WU D， et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies[J]. Nucleic Acids Res， 2015， 43（7）：e47.

[16]" "NEWMAN A M， LIU C L， GREEN M R， et al. Robust enumeration of cell subsets from tissue expression profiles[J]. Nat Methods， 2015， 12（5）：453-457.

[17]" "ROBIN X， TURCK N， HAINARD A， et al. pROC： an open-source package for R and S+ to analyze and compare ROC curves[J]. BMC Bioinformatics， 2011， 12：77.

[18]" "WANG N， GUO H， DONG Z， et al. Establishment and validation of a 7-microRNA prognostic signature for non-small cell lung cancer[J]. Cancer Manag Res， 2018， 10：3463-3471.

[19]" "FEI X F， LI Z J， YANG D Y， et al. Neddylation of Coro1a determines the fate of multivesicular bodies and biogenesis of extracellular vesicles[J]. J Extracell Vesicles， 2021， 10（12）：e12153.

[20]" "DARABI S， ARIAEI A， RUSTAMZADEH A， et al. Cerebrospinal fluid and blood exosomes as biomarkers for amyotrophic lateral sclerosis; a systematic review[J]. Diagn Pathol， 2024， 19（1）：47.

[21]" "DI CARO G， MARCHESI F， LAGHI L， et al. Immune cells： plastic players along colorectal cancer progression[J]. J Cell Mol Med， 2013， 17（9）：1088-1095.

[22]" "SPINNER C A， LAMSOUL I， M？魪TAIS A， et al. The E3 ubiquitin ligase Asb2α in T helper 2 cells negatively regulates antitumor immunity in colorectal cancer[J]. Cancer Immunol Res， 2019， 7（8）：1332-1344.

[23]" "B？魻TTCHER J P， BONAVITA E， CHAKRAVARTY P， et al. NK cells stimulate recruitment of cDC1 into the tumor microenvironment promoting cancer immune control[J]. Cell， 2018， 172（5）：1022-1037.e14.

[24]" "CHEN J， YANG J， WANG W， et al. Tumor extracellular vesicles mediate anti-PD-L1 therapy resistance by decoying anti-PD-L1[J]. Cell Mol Immunol， 2022， 19（11）：1290-1301.

[25]" "TUMEH P C， HARVIEW C L， YEARLEY J H， et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance[J]. Nature， 2014， 515（7528）：568-571.

[26]" "HEGDE P S， KARANIKAS V， EVERS S. The where， the when， and the how of immune monitoring for cancer immunotherapies in the era of checkpoint inhibition[J]. Clin Cancer Res， 2016， 22（8）：1865-1874.

[收稿日期] 2024-11-08