WGCNA聯(lián)合機(jī)器學(xué)習(xí)法探究缺氧與慢性髓系白血病潛在靶點(diǎn)與作用機(jī)制研究

摘要：目的" 通過WGCNA聯(lián)合機(jī)器學(xué)習(xí)法鑒定缺氧與慢性髓系白血?。–ML）潛在關(guān)鍵基因，并探討關(guān)鍵基因與CML中22種免疫細(xì)胞的組成和成分分布。方法" 將高通量基因表達(dá)數(shù)據(jù)庫（GEO）中的兩個CML數(shù)據(jù)集（GSE5550和GSE24739）進(jìn)行整合，消除批次效應(yīng)。采用微陣列數(shù)據(jù)線性模型（LIMMA）在25例CML和16例對照中對200個缺氧基因進(jìn)行差異分析。采用加權(quán)基因共表達(dá)網(wǎng)絡(luò)分析（WGCNA）篩選出相關(guān)模塊，聯(lián)合SVM-RFE分析和LASSO回歸模型篩選關(guān)鍵基因。三種算法做Venn分析。免疫浸潤中的CIBERSORT算法評估CML中22種免疫細(xì)胞組分的組成模式以及關(guān)鍵基因與22種免疫細(xì)胞相關(guān)性。結(jié)果" 共鑒定出34個與CML相關(guān)缺氧基因，其中14個基因上調(diào)，20個基因下調(diào)。GO富集顯示，與長鏈脂肪酸轉(zhuǎn)運(yùn)途徑、血管傷口愈合、B細(xì)胞穩(wěn)態(tài)、負(fù)調(diào)節(jié)淋巴細(xì)胞凋亡以及蛋白酶結(jié)合途徑相關(guān)。KEGG分析顯示，缺氧誘導(dǎo)因子1（HIF-1）信號通路、剪切應(yīng)力和動脈粥樣硬化信號通路明顯富集。ANXA2被鑒定為CML進(jìn)展中的關(guān)鍵基因，表達(dá)水平較正常血液標(biāo)本明顯上調(diào)（Plt;0.01），ANXA2在CML診斷模型中表現(xiàn)出較強(qiáng)的預(yù)測能力（AUC=0.820，95%CI=0.692～0.935）。免疫細(xì)胞浸潤分析顯示，ANXA2與漿細(xì)胞、細(xì)胞毒性T細(xì)胞、巨噬細(xì)胞M0型和肥大細(xì)胞活化等相關(guān)。結(jié)論" 缺氧基因ANXA2可作為CML的免疫檢查點(diǎn)，可能在CML的進(jìn)展中起到生物學(xué)標(biāo)志物的作用。

關(guān)鍵詞：慢性髓系白血病;缺氧;WGCNA;膜聯(lián)蛋白A2

中圖分類號：R733.72" " " " " " " " " " " " " " " " " " "文獻(xiàn)標(biāo)識碼：A" " " " " " " " " " " " " "DOI：10.3969/j.issn.1006-1959.2024.19.001

文章編號：1006-1959（2024）19-0001-10

Study on the Potential Targets and Mechanism of Hypoxia and Chronic Myeloid Leukemia

by WGCNA Combined with Machine Learning Method

FENG Xiaoyun1，QIN Yufeng2，LUO Qingyuan1，MO Jun1，HUANG Xiaobin1，ZHANG Peng1

（1.Center of Tissue Engineering and Stem Cell，Guizhou Medical University，Guiyang 550004，Guizhou，China;

2.Department of Prosthetic Implantology，Affiliated Stomatological Hospital of Guizhou Medical University，Guiyang 550004，Guizhou，China）

Abstract：Objective" To identify the potential key genes of hypoxia and chronic myeloid leukemia （CML） by WGCNA combined with machine learning method， and to explore the composition and distribution of key genes and 22 immune cells in CML.Methods" Two CML datasets （GSE5550 and GSE24739） in the high-throughput gene expression database （GEO） were integrated to eliminate batch effects. The microarray data linear model （LIMMA） was used to analyze the differences of 200 hypoxia genes in 25 CML patients and 16 controls. The related modules were screened by weighted gene co-expression network analysis （WGCNA）， and the key genes were screened by SVM-RFE analysis and LASSO regression model. Three algorithms do Venn analysis. The CIBERSORT algorithm in immune infiltration was used to evaluate the composition pattern of 22 immune cell components in CML and the correlation between key genes and 22 immune cells.Results" A total of 34 CML-related hypoxia genes were identified， of which 14 genes were up-regulated and 20 genes were down-regulated. GO enrichment showed that it was related to long-chain fatty acid transport pathway， vascular wound healing， B cell homeostasis， negative regulation of lymphocyte apoptosis and protease binding pathway. KEGG analysis showed that hypoxia-inducible factor 1 （HIF-1） signaling pathway， shear stress and atherosclerosis signaling pathway were significantly enriched. ANXA2 was identified as a key gene in the progression of CML， and its expression level was significantly higher than that of normal blood samples （Plt;0.01）. ANXA2 showed strong predictive ability in CML diagnostic model （AUC=0.820， 95%CI=0.692-0.935）. Immune cell infiltration analysis showed that ANXA2 was associated with plasma cells， cytotoxic T cells， macrophage M0 type and mast cell activation.Conclusion" Hypoxia gene ANXA2 can be used as an immune checkpoint of CML， which may play a role as a biological marker in the progression of CML.

Key words：Chronic myeloid leukemia;Hypoxia;WGCNA;Annexin A2

慢性髓系白血?。╟hronic myeloid leukemia， CML）是骨髓中髓細(xì)胞惡變產(chǎn)生大量異常粒細(xì)胞所導(dǎo)致的一種疾病。生理學(xué)上，骨髓顯示平均氧分壓在5%～7%[1]。在CML進(jìn)展過程中，由于腫瘤細(xì)胞的快速增殖減少了O2的擴(kuò)散，以及腫瘤內(nèi)異常血管的灌流受損，導(dǎo)致了低氧條件的形成[2]。研究表明[3]，腫瘤組織中的氧分壓低于1.3%，遠(yuǎn)低于生理性氧分壓。骨髓組織結(jié)構(gòu)至少由兩個不同的龕結(jié)構(gòu)組成，其中骨髓內(nèi)膜龕是骨髓中最缺氧的區(qū)域，靠近骨骼。而血管龕氧合較高，靠近血管[4]。在骨髓造血干細(xì)胞（hemopoietic stem cell， HSC）微環(huán)境中，氧分壓是重要環(huán)境因素，缺氧導(dǎo)致了骨髓HSC不能獲得充分氧或者不能利用氧時的增殖抑制現(xiàn)象[4]。缺氧誘導(dǎo)因子1亞基α（hypoxia inducible factor 1 subunit alpha， HIF1α）是缺氧相關(guān)的關(guān)鍵分子，在內(nèi)膜邊緣高表達(dá)，而在血管龕中低表達(dá)。HSC主要定位于具有骨髓內(nèi)膜龕的骨髓，表明缺氧信號傳導(dǎo)在干細(xì)胞代謝和調(diào)節(jié)干細(xì)胞生物學(xué)方面發(fā)揮著重要作用，有助于維持白血病中干細(xì)胞的存活和正常造血功能[5]。CML缺氧特征與腫瘤免疫抑制和免疫逃逸之間存在關(guān)聯(lián)，缺氧會降低NK細(xì)胞的自然殺傷能力，免疫系統(tǒng)的改變會影響CML患者的預(yù)后[6]。CML患者免疫細(xì)胞亞群存在顯著異質(zhì)性，與正常血液相比，CML患者血液的免疫特征是CD3+T細(xì)胞、輔助性T細(xì)胞和CD45+記憶T細(xì)胞增加。但CML骨髓微環(huán)境中免疫細(xì)胞浸潤景觀特點(diǎn)尚不清楚[7]。ANXA2是一種鈣依賴性磷脂結(jié)合蛋白，在各種真核細(xì)胞中普遍表達(dá)。研究表明[8]，缺氧誘導(dǎo)因子-1（hypoxia inducible factor， HIF-1）可誘導(dǎo)ANXA2的表達(dá)，促進(jìn)細(xì)胞的信號轉(zhuǎn)導(dǎo)和新生血管生成，但ANXA2在CML中缺氧環(huán)境下的作用并不清楚。本研究旨在通過生物信息學(xué)篩選CML患者骨髓的缺氧關(guān)鍵基因ANXA2，并進(jìn)一步分析關(guān)鍵基因在CML患者骨髓的免疫細(xì)胞浸潤特點(diǎn)，以期為CML的免疫治療提供潛在的可能靶點(diǎn)。

1資料與方法

1.1 CML微陣列數(shù)據(jù)" 在GEO（https：//www.ncbi.nlm.nih.gov/geo/）下載關(guān)于CML的GSE5550和GSE24739 mRNA表達(dá)譜數(shù)據(jù)。200個缺氧基因來源于分子特征數(shù)據(jù)庫v7.0（MSigDB7.0，www.gsea-msigdb.org）。GSE5550包含9例CML患者和8例正常對照的血液標(biāo)本，GSE24739包含16例CML患者和8例正常對照組的血液標(biāo)本。

1.2 CML缺氧差異表達(dá)基因的數(shù)據(jù)處理與鑒定" 利用R對兩個原始數(shù)據(jù)集進(jìn)行預(yù)處理，包括背景校正、歸一化和log2轉(zhuǎn)換。當(dāng)多個探針對應(yīng)一個共同基因時，取其平均值作為表達(dá)值。此外，Bioconductor的代理變量分析“SVA”包用于去除兩個數(shù)據(jù)集之間的批次效應(yīng)和其他不必要的變異。最后，利用“LIMMA”包篩選缺氧與CML差異基因（differentially expressed genes， GEGs）。其中Plt;0.05，|log2 Fold change （FC）|gt;0.5作為篩選DEGs的條件。

1.3 CML中缺氧DEGs功能富集分析" 根據(jù)基因本體論（gene ontology， GO）和京都基因與基因組百科全書通路（Kyoto encyclopedia of genes and genomes， KEGG）分析CML中缺氧DEGs的功能富集情況，通過R包中的“clusterProfiler”和“DOSE”包完成。Plt;0.05作為顯著性的富集通路。

1.4潛在生物學(xué)標(biāo)志物的篩選" LASSO是一種利用正則化方法提高預(yù)測精度的回歸分析算法，以確定候選基因集。通過R中的“glmnet”包完成LASSO回歸算法，確定與CML與骨髓健康樣本鑒別能力顯著相關(guān)的基因集。支持向量機(jī)（support vector machine， SVM）是一種監(jiān)督式的機(jī)器學(xué)習(xí)技術(shù)，廣泛用于分類和回歸分析，為了避免過擬合，利用RFE算法從數(shù)據(jù)隊(duì)列中篩選最優(yōu)基因。采用特征選擇之支持向量機(jī)遞歸特征消除（SVM recursive feature elimination， SVM-RFE）算法篩選合適的特征基因。根據(jù)基因之間的相互作用，應(yīng)用加權(quán)基因共表達(dá)網(wǎng)絡(luò)分析（weighted gene co-expression network analysis， WGCNA）構(gòu)建基因共表達(dá)網(wǎng)絡(luò)，以確定潛在的候選基因。樣本中超過25%變異基因?qū)⒓傻臄?shù)據(jù)集導(dǎo)入R語言中的“WGCNA”模塊。為了保證網(wǎng)絡(luò)構(gòu)建結(jié)果的可靠性，對離群樣本進(jìn)行剔除。采用選取-軟閾值函數(shù)，由共表達(dá)式相似度得到的軟閾值冪β計(jì)算鄰接度，然后將鄰接變換為拓?fù)渲丿B矩陣（topological overlap matrix， TOM），并計(jì)算相應(yīng)的不相似度（1-TOM）。通過層次聚類和動態(tài)樹切函數(shù)進(jìn)行模塊檢測。為了將表達(dá)譜相似的基因分類到基因模塊中，根據(jù)基于TOM的相異度度量進(jìn)行平均連鎖層次聚類，基因樹狀圖的最小大小基因組為50。對與臨床屬性相關(guān)的模塊，計(jì)算模塊隸屬度（module membership， MM）和基因顯著性（gene significance， GS）。將特征基因網(wǎng)絡(luò)可視化，模塊中的基因信息用于進(jìn)一步分析。本研究的篩選方式在于將LASSO分析、SVM-RFE算法和WGCNA分析做進(jìn)一步Venn分析，得到CML缺氧相關(guān)基因ANXA2。

1.5 ANXA2的表達(dá)及診斷意義" 通過R中的“ggplot2”和“ggpubr”可視化ANXA2差異表達(dá)。為了進(jìn)一步了解ANXA2對CML的臨床診斷意義，通過“pROC”包計(jì)算受試者工作特征曲線（receiver operating characteristic curve， ROC）曲線下面積（area under curve， AUC）。當(dāng)AUCgt;0.7時，表示臨床診斷準(zhǔn)確性較高。

1.6免疫細(xì)胞浸潤分析" CIBERSORT分析的計(jì)算方法是一種建立在遺傳表達(dá)式基礎(chǔ)上的反卷積算法，用于評估一個基因組相對于樣本中其他基因的變化。借助CIBERSORT算法分析22個免疫細(xì)胞的免疫反應(yīng)，并評估正常樣本和CML樣本中這些免疫細(xì)胞與ANXA2的相關(guān)性。為了量化CML中浸潤免疫細(xì)胞的相對比例，采用組織中不同免疫細(xì)胞類型的CIBERSORT方法分析合并表達(dá)數(shù)據(jù)和計(jì)算免疫細(xì)胞浸潤。以Plt;0.05為標(biāo)準(zhǔn)對樣品進(jìn)行過濾。并計(jì)算每種免疫細(xì)胞在樣本中的百分比，并用柱狀圖顯示。22個免疫細(xì)胞的熱圖使用“pheatmap”包制作?！皏ioplot”包用于比較和可視化CML和對照樣本之間的22個免疫細(xì)胞的水平。使用“corrplot”軟件包進(jìn)行了相關(guān)熱圖，揭示ANXA2與22種浸潤免疫細(xì)胞的相關(guān)性。

1.7統(tǒng)計(jì)學(xué)方法" 本研究采用R語言4.2.0進(jìn)行分析，相關(guān)語言包直接從Bioconductor（http：//www.bioconductor.org/）獲取，包括“SVA”“LIMMA”clusterProfiler”“DOSE”“glmnet”“WGCNA”“ggplot2”“ggpubr”“pROC”“pheatmap”“vioplot”和“corrplot”包。當(dāng)Plt;0.05表明差異有統(tǒng)計(jì)學(xué)意義。

2結(jié)果

2.1 CML中缺氧DEGs的鑒定" 數(shù)據(jù)集GSE55235和GSE55457包含25個CML數(shù)據(jù)和16個對照數(shù)據(jù)。去除批處理效應(yīng)后，通過Limma包共鑒定出34個缺氧DEGs，其中14個基因顯著向上調(diào)（CHST2、LXN、TGM2、ANXA2、PFKP、ADORA2B、PLIN2、VEGFA、SLC37A4、CASP6、HK1、MIF、SDC4、SLC2A1），20個基因顯著下調(diào)（NDST2、COL5A1、TGFBI、RS2、BTG1、FOS、ZFP36、ISG20、PNRC1、NAGK、CXCR4、DUSP1、AKAP12、KLF7、EFNA1、BCL2、TNFAIP3、CCNG2、NR3C1、RRAGD）。缺氧DEGs的熱圖和火山圖見圖1。

2.2 CML功能富集分析" 通過R軟件包“ClusterProfil”對CML中34個缺氧DEGs的GO和KEGG進(jìn)行分析，結(jié)果顯示34種缺氧DEGs主要參與化學(xué)引誘活性、趨化因子結(jié)合、纖連蛋白結(jié)合、14-3-3結(jié)合蛋白、C-C趨化因子結(jié)合、糖胺聚糖結(jié)合、蛋白酶結(jié)合、碳水化合物激酶活性、單糖代謝過程、白細(xì)胞增殖的調(diào)節(jié)、細(xì)胞對細(xì)胞外刺激的反應(yīng)、B細(xì)胞增殖的正向調(diào)節(jié)、脂肪酸轉(zhuǎn)運(yùn)、淋巴細(xì)胞凋亡過程的負(fù)調(diào)節(jié)、血管生成參與傷口愈合、B細(xì)胞穩(wěn)態(tài)、血管傷口愈合以及長鏈脂肪酸轉(zhuǎn)運(yùn)途徑。同時，KEGG表明HIF-1信號通路、剪切應(yīng)力和動脈粥樣硬化通路明顯富集，見圖2。

2.3關(guān)鍵特征基因的篩選" 首先，針對可能的診斷因素，采用3種機(jī)器學(xué)習(xí)算法來預(yù)測CML診斷標(biāo)志物。通過LASSO回歸算法對缺氧DEGs進(jìn)行篩選鑒定，確定了10個變量作為CML的診斷標(biāo)志（圖3A、3B）。進(jìn)一步通過SVM-RFE算法（圖3C）在GEGs中識別出16個特征的子集。使用WGCNA算法識別基因共表達(dá)網(wǎng)絡(luò)模塊。首先，將基因從大到小進(jìn)行排序，選擇方差前25%的基因進(jìn)行分析。其次，應(yīng)“flashClust”R軟件包進(jìn)行聚類分析，將閾值設(shè)置為80個，然后檢測并剔除離群值樣本。第三，利用WGCNA軟件包中的“選擇軟閾值”功能，篩選出1～20之間的功率參數(shù)。本研究中選擇β=8（無尺度R2=0.85）的功率作為軟閾值，以確保一個尺度網(wǎng)絡(luò)（圖3D）。將閾值設(shè)置為0.25，以合并集群樹中的類似模塊，共獲得13個模塊，其中基因具有相似的共表達(dá)特征。隨機(jī)選擇所有模塊的顏色來區(qū)分不同的模塊。與其他模塊相比，特征基因（ME）表現(xiàn)出最高的正相關(guān)（R=0.64，P=7E-06，圖3E）。因此，包含357個基因的MEbrown模塊被確定為CML的關(guān)鍵模塊，并保留以供進(jìn)一步分析。此外，對MEbrown模塊中基因模塊成員關(guān)系特定基因（module memberships， MMs）和模塊特征基因（gene significances， GSs）進(jìn)行顯著性計(jì)算，結(jié)果顯示MM和GS之間的相關(guān)系數(shù)較高（Cor=0.5，P=4.8E-24，圖3F）。從集成數(shù)據(jù)集中篩選出的缺氧DEGs與MEbrown模塊中的基因交集為共同基因（圖3G），且MEbrown模塊的基因與CML相關(guān)。最終對WGCNA、LASSO和RFE三種算法進(jìn)行取交集，獲得重疊特征基因?yàn)锳NXA2（圖3H）。

2.4 ANXA2差異表達(dá)分析和ROC檢測" 根據(jù)ANXA2差異表達(dá)結(jié)果分析發(fā)現(xiàn)，與健康樣本相比，CML樣本中ANXA2的表達(dá)水平明顯上調(diào)，同時也表明ANXA2參與了CML的發(fā)生。ROC進(jìn)一步分析發(fā)現(xiàn)，ANXA2在CML診斷模型中表現(xiàn)出較強(qiáng)的能力（AUC=0.820，95%CI=0.692～0.935），見圖4。

2.5免疫細(xì)胞浸潤分析" R語言中的CIBERSORT算法分析了CML每個樣本中22種免疫細(xì)胞的百分比，其中B祖細(xì)胞和B記憶細(xì)胞占比最高，中性粒細(xì)胞和嗜酸性粒細(xì)胞占比最低（圖5）。CML每個樣本中22種免疫細(xì)胞的免疫細(xì)胞浸潤熱圖顯示，GSM129066、GSM609348、GSM609359、GSM609361、GSM609364、GSM609366和GSM609367樣本中肥大細(xì)胞靜息比例顯著上調(diào)，GSM609348樣本中漿細(xì)胞與GSM609362中嗜酸性粒細(xì)胞也顯著上調(diào)，而其他樣本中免疫細(xì)胞則多數(shù)顯示下調(diào)或者沒有發(fā)揮作用（圖6）。免疫細(xì)胞浸潤差異的小提琴圖表明，與對照組相比，CML患者的濾泡輔助性T細(xì)胞水、初始CD4T細(xì)胞和初始B細(xì)胞顯著上調(diào)，而巨噬細(xì)胞M2顯著下調(diào)（圖7）。ANXA2與22種免疫細(xì)胞的相關(guān)性分析顯示，ANXA2與巨噬細(xì)胞M2型、泡輔助性T細(xì)胞呈負(fù)相關(guān)（R=-0.38，Plt;0.05，圖8）。且ANXA2與漿細(xì)胞、CD4記憶T細(xì)胞靜息、自然殺傷細(xì)胞靜息、巨噬細(xì)胞、樹突狀細(xì)胞靜息、樹突狀細(xì)胞活化、肥大細(xì)胞靜息和肥大細(xì)胞活化相關(guān)（圖9）。

3討論

盡管CML患者的十年生存有效期顯著提高，但仍然有部分患者表現(xiàn)出免疫抑制、耐藥性、骨髓抑制和胃腸道反應(yīng)[9-11]。在CML患者中，腫瘤細(xì)胞的生長優(yōu)勢導(dǎo)致產(chǎn)生比正常造血更成熟的細(xì)胞群，并且臨床進(jìn)化更長。迄今為止，盡管大多數(shù)CML和小兒急性淋巴白血病（acute lymphoblastic leukemia， ALL）病例在治療下已治愈或得到良好控制，但仍有很高的復(fù)發(fā)風(fēng)險[12-14]。缺氧誘導(dǎo)腫瘤微環(huán)境中多種促血管生成因子的表達(dá)能夠影響類似于HSC的白血病干細(xì)胞（leukemic stem cells， LSCs）的循環(huán)、靜止、分化、代謝和治療耐藥性。LSCs能夠與骨髓微環(huán)境相互作用，影響周圍細(xì)胞的表型，從而獲得更強(qiáng)的增殖潛力，使腫瘤繼續(xù)進(jìn)展[15]。缺氧與LSCs之間的相關(guān)性提示，缺氧可能在調(diào)節(jié)腫瘤內(nèi)異質(zhì)性和免疫逃逸中發(fā)揮核心作用，信號通路可能和HIF1α相關(guān)[16]。缺氧還可以下調(diào)大多數(shù)NK細(xì)胞受體的表達(dá)和功能，這些受體直接負(fù)責(zé)對腫瘤細(xì)胞發(fā)揮細(xì)胞溶解活性[17]。

本研究通過WGCNA聯(lián)合機(jī)器學(xué)習(xí)法發(fā)現(xiàn)，ANXA2為缺氧與CML的共同關(guān)鍵基因。ANXA2是一種鈣離子依賴性的磷脂結(jié)合蛋白，在內(nèi)皮細(xì)胞、巨噬細(xì)胞、單核細(xì)胞和各種類型的癌細(xì)胞表面表達(dá)。在C-末端結(jié)構(gòu)域包含陰離子磷脂和鈣離子結(jié)合特性，而N-末端是翻譯后修飾和與其他蛋白質(zhì)相互作用的部位，其功能涉及肌動蛋白重組、細(xì)胞信號轉(zhuǎn)導(dǎo)、胞吞、胞吐、轉(zhuǎn)錄、細(xì)胞膜構(gòu)成、新生血管生成以及DNA復(fù)制與修復(fù)[18，19]。據(jù)研究報(bào)道[20，21]，ANXA2被認(rèn)為是膀胱癌、透明細(xì)胞腎細(xì)胞癌、肝細(xì)胞癌和前列腺癌的生物標(biāo)志物。在細(xì)胞表面，ANXA2異四聚體為纖溶酶原和組織型纖溶酶原激活物提供結(jié)合部位，將纖溶酶原轉(zhuǎn)化為纖溶酶，在腫瘤的增殖、侵襲和轉(zhuǎn)移中發(fā)揮重要作用[22]。本研究通過生物信息學(xué)分析發(fā)現(xiàn)，ANXA2在CML中表達(dá)上調(diào);進(jìn)一步通過ROC驗(yàn)證發(fā)現(xiàn)，ANXA2有較高的臨床診斷價值，表明ANXA2參與了CML的發(fā)生和發(fā)展。目前，臨床關(guān)于ANXA2在CML免疫浸潤中的作用研究有限。通過CIBERSORT算法分析發(fā)現(xiàn)，ANXA2與漿細(xì)胞、CD4記憶T細(xì)胞靜息、自然殺傷細(xì)胞靜息、巨噬細(xì)胞、樹突狀細(xì)胞靜息、樹突狀細(xì)胞活化、肥大細(xì)胞靜息和肥大細(xì)胞活化相關(guān)，表明ANXA2的上調(diào)可能引發(fā)CML產(chǎn)生免疫反應(yīng)，揭示ANXA2作為HIF1α下游效應(yīng)分子具有潛在免疫靶向作用。

綜上所述，通過構(gòu)建CML與缺氧相關(guān)的預(yù)后特征基因及免疫相關(guān)性，預(yù)測了ANXA2在CML中的潛在免疫作用，提示ANXA2可能是CML診斷的生物學(xué)標(biāo)志物。但本研究仍存在一定的局限性，這項(xiàng)研究是通過計(jì)算機(jī)生物學(xué)回顧性研究的，仍需要更多的臨床樣本驗(yàn)證和分子生物學(xué)實(shí)驗(yàn)驗(yàn)證。

參考文獻(xiàn)：

[1]Cutica I，Riva S，Orlandi EM，et al.Psychological Factors Affecting the Willingness to Accept a Possible Tyrosine Kinase Inhibitor （TKI） Discontinuation in Chronic Myeloid Leukaemia （CML） Patients[J].Patient Prefer Adherence，2022，16：2963-2975.

[2]Copland M.Treatment of blast phase chronic myeloid leukaemia：A rare and challenging entity[J].Br J Haematol，2022，199（5）：665-678.

[3]Krishnan V，Schmidt F，Nawaz Z，et al.A single-cell atlas identifies pretreatment features ofprimary imatinib resistance in chronic myeloid leukemia[J].Blood，2023，141（22）：2738-2755.

[4]Abruzzese E，Bocchia M，Trawinska MM，et al.Minimal Residual Disease Detection at RNA and Leukemic Stem Cell （LSC） Levels： Comparison of RT-qPCR，d-PCR and CD26+ Stem Cell Measurements in Chronic Myeloid Leukemia （CML） Patients in Deep Molecular Response （DMR） [J].Cancers （Basel），2023，15（16）：4112.

[5]Qiu S，Sheth V，Yan C，et al.Metabolic adaptation to tyrosine kinase inhibition in leukemia stem cells[J].Blood，2023，142（6）：574-588.

[6]Saussele S，Richter J，Guilhot J，et al.Discontinuation of tyrosine kinase inhibitor therapy in chronic myeloid leukaemia （EURO-SKI）： a prespecified interim analysis of a prospective，multicentre，non-randomised，trial[J].Lancet Oncol，2018，19（6）：747-757.

[7]Yao F，Zhao C，Zhong F，et al.Bioinformatics analysis and identification of hub genes and immune-related molecular mechanisms in chronic myeloid leukemia[J].Peer J，2022，10：e12616.

[8]Panuzzo C，Pironi L，Maglione A，et al.mTORC2 Is Activated under Hypoxia and Could Support Chronic Myeloid Leukemia Stem Cells[J].Int J Mol Sci，2023，24（2）：1234.

[9]Amir M，Javed S.A Review on the Therapeutic Role of TKIs in Case of CML in Combination With Epigenetic Drugs[J].Front Genet，2021，12：742802.

[10]Gleixner KV，F(xiàn)ilik Y，Berger D，et al.Asciminib and ponatinib exert synergistic anti-neoplastic effects on CML cells expressing BCR-ABL1T315I-compound mutations[J].Am J Cancer Res，2021，11（9）：4470-4484.

[11]Patterson SD，Copland M.The Bone Marrow Immune Microenvironment in CML： Treatment Responses，Treatment-Free Remission，and Therapeutic Vulnerabilities[J].Curr Hematol Malig Rep，2023，18（2）：19-32.

[12]Madabhavi I，Patel A，Anand A，et al.Severe persistent marrow aplasia in chronic myeloid leukemia （CML）： An unusual complication of imatinib[J].J Cancer Res Ther，2023，19（2）：335-339.

[13]Bührer ED，Amrein MA，F(xiàn)orster S，et al.Splenic red pulp macrophages provide a niche for CML stem cells and induce therapy resistance[J].Leukemia，2022，36（11）：2634-2646.

[14]Langhammer M，Sch?pf J，Jaquet T，et al.Mast cell deficiency prevents BCR：：ABL1 induced splenomegaly and cytokine elevation in a CML mouse model[J].Leukemia，2023，37（7）：1474-1484.

[15]Shalhout SZ，Yang PY，Grzelak EM，et al.YAP-dependent proliferation by a small molecule targeting annexin A2[J].Nat Chem Biol，2021，17（7）：767-775.

[16]Zhao S，Li B，Zhao R，et al.Hypoxia-induced circADAMTS6 in a TDP43-dependent manner accelerates glioblastoma progression via ANXA2/NF-KB pathway[J].Oncogene，2023，42（2）：138-153.

[17]Hossen MB，Islam MA，Reza MS，et al.Robust identification of common genomic biomarkers from multiple gene expression profiles for the prognosis， diagnosis， and therapies of pancreatic cancer[J].Comput Biol Med，2023，152：106411.

[18]Zhang H，Zhu X，F(xiàn)riesen TJ，et al.Annexin A2/TLR2/MYD88 pathway induces arginase 1 expression in tumor-associated neutrophils[J].J Clin Invest，2022，132（22）：e153643.

[19]Yang W，Mei FC，Cheng X.EPAC1 regulates endothelial annexin A2 cell surface translocation and plasminogen activation[J].FASEB J，2018，32（4）：2212-2222.

[20]Ning Y，Li Y，Wang H.ANXA2 is a potential biomarker for cancer prognosis and immune infiltration： A systematic pan-cancer analysis[J].Front Genet，2023，14：1108167.

[21]Murkovic M，Babarovic E，Marijic B，et al.Association of pre-treatment bone marrow morphology and achievement of BCR-ABL1 transcript milestones in CML[J].Pathol Res Pract，2023，246：154517.

[22]Li Z，Pan Y，Yao J，et al.ANXA2 as a novel substrate of FBXW7 promoting esophageal squamous cell carcinoma via ERK phosphorylation[J].Biochem Biophys Res Commun，2023，649：93-100.

收稿日期：2023-09-25;修回日期：2023-10-11

編輯/杜帆

基金項(xiàng)目：1.國家自然科學(xué)基金地區(qū)科學(xué)基金項(xiàng)目（編號：81960036）;2.貴州省科技計(jì)劃項(xiàng)目（編號：黔科合基礎(chǔ)-ZK[2021]重點(diǎn)005）

作者簡介：馮小云（1998.11-），男，貴州貴陽人，碩士研究生，主要從事干細(xì)胞表觀遺傳學(xué)研究

通訊作者：張鵬（1987.4-），男，甘肅張掖人，博士，教授，博士生導(dǎo)師，主要從事干細(xì)胞表觀遺傳學(xué)研究