-
摘要:目的
探讨头颈部区域脂肪肉瘤的发病以及预后影响因素,同时分析不同治疗方案的疗效。
方法回顾性分析2008年1月至2024年1月本院接治的头颈部脂肪肉瘤患者。所有患者均接受随访,利用SPSS软件分析患者预后。
结果共纳入患者30例,其中脂肪肉瘤在眼眶内可达60%,其余脂肪肉瘤主要分布在颈深部各个间隙中。去分化脂肪肉瘤(33%)最常见,黏液多形性脂肪肉瘤(4%)最少见。肿瘤病理类型(P<0.001)和Ki67(P=0.014)对肿瘤控制率有显著影响。针对肿瘤的疾病特异性生存率进行分析,显示各项因素间差异无统计学意义(均P>0.05)。
结论头颈部脂肪肉瘤的预后相对于其他部位的脂肪肉瘤预后更好,但黏液多形性脂肪肉瘤或多形性脂肪肉瘤以及Ki67过高是预后差的信号,而术后的辅助放射治疗并不能显著提高疾病特异性生存率。
Abstract:ObjectiveTo explore the pathogenesis and prognostic factors of liposarcoma in the head and neck region, and simultaneously analyze the efficacy of different treatment regimens.
MethodsA retrospective analysis was performed on all patients with primary untreated head and neck liposarcoma who were diagnosed and underwent surgical treatment at our hospital from January 2008 to January 2024. All patients were monitored during follow-up, and their prognoses were analyzed using SPSS software.
ResultsA total of 30 patients were included in the study. Liposarcoma accounted for up to 60% of the cases in the orbit, while the remaining liposarcomas were primarily located in various interspaces of the neck. Dedifferentiated liposarcoma was the most common type, comprising 33%, while myxoid pleomorphic liposarcoma was the rarest at 4%. The tumor pathological type (P<0.001) and Ki67 (P=0.014) significantly affected the tumor control rate. However, an analysis of disease-specific survival rates revealed no significant differences across various factors (all P>0.05).
ConclusionThe prognosis of head and neck liposarcoma is better compared to that of liposarcomas in other parts of the body. However, myxoid pleomorphic liposarcoma, pleomorphic fat sarcoma, and high Ki67 levels are indicators of poor prognosis. Additionally, postoperative adjuvant radiotherapy does not significantly enhance disease-specific survival rates.
-
Key words:
- Head and neck /
- Liposarcoma /
- Prognosis analysis /
- Treatment
-
0 引言
胶质母细胞瘤(Glioblastoma, GBM)被认为起源于神经胶质干细胞或祖细胞[1],是级别较高的恶性胶质瘤(Ⅳ级),也是中枢神经系统(Central nervous system, CNS)最常见的原发性恶性肿瘤,占所有CNS肿瘤的14.5%,占CNS恶性肿瘤的48.6%。GBM患者的中位总生存期(Overall survival, OS)仅为15个月[2-3]。GBM是高度血管化的肿瘤,其生长依赖于新生血管的形成。血管新生是血管内皮细胞(Endothelial cells, ECs)在特定信号刺激下增殖、迁移和分化的复杂过程[4]。内皮细胞除了参与血管生成,进一步为肿瘤提供氧气和营养物质外,还参与肿瘤血管内侵袭,使肿瘤细胞转移到血管腔内,促进肿瘤的侵袭和转移[5]。因此,内皮细胞在GBM的发生、发展中发挥着重要的作用。相较于肿瘤细胞标志物,内皮细胞相关标志物或可以更全面地评估肿瘤进展、反映肿瘤的预后。
寻找可靠的生物标志物来预测GBM患者的预后,并开发新的分子靶向治疗策略,对于目前GBM的治疗至关重要。然而,GBM目前尚缺乏有效的预后标志物。本研究旨在通过多维度生物信息学分析方法,构建GBM内皮细胞风险评分,探讨并验证内皮细胞在GBM中的预后价值。
1 资料与方法
1.1 GBM转录组数据获取与标准化
GBM转录组数据以及相关的临床资料获取于癌症基因组图谱(The Cancer Genome Atlas, TCGA)数据库(TCGA-GBM,包含174名GBM患者),其中转录组数据被统一标准化为每百万条reads的转录本(Transcripts per million reads, TPM)格式,并进行log2转化。中国脑胶质瘤基因组图谱(Chinese Glioma Genome Atlas, CGGA)中GBM相关数据(CGGA-GBM,包含388名GBM患者)获取于GlioVis数据库(http://gliovis.bioinfo.cnio.es)。
1.2 GBM单细胞测序分析
GBM单细胞数据(GSE103224)获取于开源单细胞数据库TSCH2数据库(http://tisch.comp-genomics.org/),相关细胞分群结果直接获取于该数据库。内皮细胞相关基因标志物按照如下标准筛选:在满足P.adj<0.05条件下,按照log2FC的绝对值排序在上调(Up-regulated)和下调(Down-regulated)基因中各取前100进行后续分析。
1.3 单因素Cox回归分析
单因素Cox回归分析获取TCGA-GBM预后相关基因,相关筛选标准:风险比(Hazard ratio, HR)>1且P.adj<0.05或HR<1且P.adj<0.05。
1.4 内皮细胞预后相关标志物获取
使用韦恩图对上述四种标志物取交集,获取的内皮细胞预后相关标志物满足条件:HR>1且为上调基因或HR<1且为下调基因。使用森林图展示相关标志物及其预后相关信息。
1.5 LASSO回归分析
基于TCGA-GBM转录组数据,使用LASSO回归分析进一步筛选内皮细胞预后相关标志物,并构建内皮细胞相关预后风险评分,选取交叉验证误差均值最小对应的lambda值(lambdamin),并获取风险评分中各基因的系数,构成风险评分(Risk score)公式:
$$ Risk\;score = \sum\nolimits_{i = 1}^n \beta i*Expi $$ 其中Exp为基因的表达量,β为基因所对应的系数。
1.6 qPCR实验验证预后相关标志物表达差异
选取在解放军总医院第一医学中心收集的GBM样本组织及瘤周正常脑组织样本各3例进行定量PCR实验,验证用于构建风险评分的4个目的基因(DUSP6、STC1、VWA1和TM4SF1)在GBM组织及瘤周正常脑组织中的表达差异,qPCR引物见表1。收集患者临床病理信息。所有患者及家属均签署知情同意书,本研究符合2013年修订的《赫尔辛基宣言》要求。根据公式整理和计算qPCR实验数据,采用两独立样本t检验进行统计{stats[4.2.1]包以及car[3.1-0]包},用ggplot2[3.3.6]包对数据进行可视化。
表 1 qPCR实验所用引物Table 1 Primers used in qPCR experimentsGene name Gene ID Forward primer Reverse primer DUSP6 1 848 5′-GAACTGTGGTGTCTTGGTACATT-3′ 5′-GTTCATCGACAGATTGAGCTTCT-3′ STC1 6781 5′-CACGAGCTGACTTCAACAGGA-3′ 5′-GGATGTGCGTTTGATGTGGG-3′ VWA1 64856 5′-GCAGACTCGGGCTACTATGTG-3′ 5′-CACGTTGGACTCAGGCACTA-3′ TM4SF1 4071 5′-TGTGGCAAACGATGTGCGA-3′ 5′-TGACACAGTAGCCAGATCCTG-3′ 1.7 生存分析
基于内皮细胞相关预后风险评分,获取每例GBM患者的风险评分数值,根据该数值的中位数,将患者分为高风险和低风险组。使用Kaplan-Meier法构建生存曲线以鉴定该预后风险评分的预后效能。上述分析及数据可视化基于软件R(4.2.1),相关R包见表2。
表 2 分析及数据可视化所用R包Table 2 R packages for analysis and data visualizationAnalysis item R package used LASSO regression
analysisglmnet (4.1.7) Single factor Cox
regression analysissurvival (3.3.1) and rms (6.3-0) Survival analysis survival (3.3.1), survminer and
ggplot2 (3.3.6)Venn diagram ggplot2 (3.3.6), VennDiagram (1.7.3) Forest plot ggplot2 (3.3.6) Co-expression heat map ggplot2 (3.3.6) 2 结果
2.1 GBM预后相关基因
单因素Cox回归分析共筛选获得2 115个GBM预后相关基因,其中1 494个满足HR>1且P.adj<0.05,621个满足HR<1且P.adj<0.05。
2.2 GBM单细胞测序分析
基于TSCH2数据库对GBM单细胞数据(GSE103224)进行解析,降维后共获得7群细胞,见图1A,分别为AC样肿瘤细胞(AC-like malignant)、内皮细胞(Endothelial)、单核/巨噬细胞(Mono/Macro)、NB样肿瘤细胞(NB-like malignant)、神经元(Neuron)、OC样肿瘤细胞(OC-like malignant)以及OPC样肿瘤细胞(OPC-like malignant)。相关细胞群标志物如图1B所示。根据内皮细胞差异基因及相应筛选标准,共获取200个内皮细胞相关基因标志物,其中100个为上调基因(CLND5、ITM2A、IGFBP7、IFI27以及ESAM等),100个为下调基因(CRYAB、CRYGS、CDK4、S100A6以及APOE等)。
图 1 GBM单细胞测序分析及内皮细胞预后相关标志物分析Figure 1 GBM single-cell sequencing analysis and endothelial cell prognostic marker analysisA: seven groups of cells obtained by analyzing GBM single-cell data (GSE103224); B: expression of seven cell-population marker genes; C: six endothelial cell prognostic markers obtained after Venn diagram analysis; D: forest plot showing the prognostic information related to the prognostic markers of endothelial cells; E: quantitative analysis of the violin plot showing that the six markers were significantly highly expressed in endothelial cells. GBM: glioblastoma.2.3 内皮细胞预后相关标志物
韦恩图分析共获取6个满足条件的交集基因,即内皮细胞预后相关标志物,分别为PLXND1、DUSP6、STC1、ESM1、TM4SF1以及VWA1基因,见图1C。森林图显示了内皮细胞预后相关标志物相关(PLXND1、DUSP6、STC1、ESM1、TM4SF以及VWA1)的预后信息,见图1D。定量分析小提琴图显示上述标志物均在内皮细胞中显著高表达,见图1E。
2.4 内皮细胞相关预后风险评分
LASSO回归分析共筛选出4个基因(DUSP6、STC1、VWA1和TM4SF1)用于风险评分的构建,变量系数谱见图2A。风险评分=0.171*DUSP6+0.144*STC1+0.041*VWA1−0.004*TM4SF1。变量轨迹图见图2B。TCGA-GBM发现集生存分析显示,风险评分是GBM患者预后不佳的影响因素(HR=1.56,95%CI:1.10~2.21,P=0.007),见图2C,CGGA-GBM验证集生存分析,风险评分亦是GBM患者预后不佳的影响因素(HR=1.40,95%CI:1.12~1.75,P=0.003)。qPCR实验结果显示DUSP6、STC1、VWA1和TM4SF1在GBM组织与瘤周正常脑组织中的表达均显著上调(P<0.001),见图3。
图 2 内皮细胞相关预后风险评分Figure 2 Endothelial cell-related prognostic risk scoreA: variable coefficient spectrum of four genes screened by LASSO regression analysis; B: six endothelial cell prognostic markers corresponding with the variable trajectory diagram; C: TCGA-GBM discovery-set survival analysis; D: CGGA-GBM validation-set survival analysis.3 讨论
GBM是成人中脑部最常见、侵袭性最强的原发性恶性肿瘤。尽管GBM患者在接受最大限度的手术切除后可以继续联合放化疗,但预后仍很差,其中位生存期仅约15个月,5年生存率<10%[6-8]。因此我们亟需探索更多的治疗方法,以提高疗效。目前尚缺乏十分可靠的生物标志物来预测当前或新疗法的预后。Sareen等[9]进行的一项系统评价和Meta分析结果显示:O6-甲基鸟嘌呤-DNA甲基转移酶(O6 -methylguanine-DNA methyltransferase, MGMT)启动子甲基化和异柠檬酸脱氢酶1(Isocitrate dehydrogenase 1, IDH1)突变与GBM患者较长的总生存期显著相关,合并风险比分别为1.66(95%CI:1.32~2.09;P<
0.0001 )和2.37(95%CI:1.81~3.12;P<0.00001 )。另有研究显示MGMT启动子甲基化水平在预测替莫唑胺(Temozolomide, TMZ)治疗反应中具有重大的临床意义[10]。1p/19q共缺失和第10染色体丢失也是TMZ反应的阳性预测指标。另一方面,错配修复系统(Mismatch repair, MMR)缺陷导致的高突变表型可预测GBM对TMZ具有耐药性。编码H3.3组蛋白的基因H3F3A的突变是儿童患者预后不良的重要标志物。MYB、MN1和MAPK通路的改变则是积极的预后因素[11]。新的研究通过“液体活检”技术在血液或脑脊液中识别GBM特异性的外泌体,从而实现对GBM患者预后的微创预测[12]。另外,细胞因子信号转导抑制因子蛋白通过JAK/STAT通路和NF-κB信号通路调控GBM的生物合成,从而发挥抑癌作用。SOCS蛋白表达的上调与GBM更好的预后相关[13]。理论上,有许多长链非编码RNA(Long non-coding RNA, lncRNA),如HOTAIR、H19和NEAT1可作为GBM的诊断和预后指标。遗憾的是,迄今为止尚未有lncRNAs成功应用于临床[14-15]。Jarmuzek等[16]对炎性反应和免疫标志物作为GBM患者预后因素的最新研究进行了全面的综述和荟萃分析,结果发现中性粒细胞/淋巴细胞比值(NLR)或血小板/淋巴细胞比值(PLR)较高的患者预后更差。同样,细胞死亡过程中释放入血的游离DNA(Cell-free DNA, cfDNA)(HR=2.35,95%CI:1.27~4.36,P<0.01)较高的患者预后更差[16]。在GBM中,一些miRNA的水平降低。研究表明,这些miRNA的过度表达可增强细胞凋亡和抑制肿瘤进展,提示预后较好[17]。国内有学者研究发现固有免疫分子CD58在高级别胶质瘤中的表达高于低级别胶质瘤和非瘤脑组织,其表达差异与胶质瘤生存期相关,故可作为判定胶质瘤的恶性程度及预后的一项指标[18]。但上述生物标志物在GBM预后的临床预测中仍存在一定的局限性。目前还需要探索更有效、更实用的GBM预后标志物。鉴于血管生成在肿瘤生长和转移中的关键作用,以及GBM的高度血管化特性[19],探究GBM血管内皮细胞相关预后标志物用于构建预后体系已成为新的研究方向。同时,抗血管生成药物也已被广泛探索作为GBM治疗的新选择[5]。GBM细胞分泌大量的促血管生成因子,营造出强烈的促血管新生的肿瘤微环境(Tumor microenvironment, TME)。GBM内皮细胞是血肿瘤屏障(Blood-tumour barrier, BTB)的重要结构基础,使化疗药物难以通过BTB到达病灶[20]。这成为GBM药物治疗的难题。通过应用单细胞测序对BTB进行更深入的了解,并通过生物信息分析方法开发新的内皮细胞生物标志物,从而实现靶向调节GBM内皮细胞的生物学行为,调控GBM的血管新生及BTB的通透性,抑制肿瘤的生长和侵袭,将为GBM提供新的治疗策略[21]。
TME中的免疫细胞在肿瘤生长过程中发挥着重要作用。近年来的多项研究表明,TME中的B细胞与改善临床预后相关。浆细胞是终末分化的B细胞,可产生高度特异性的抗体。关于临床结局与抗体产生之间关系以及Stephanie等的研究数据显示,在不分泌抗体的小鼠中,双重免疫检查点抑制剂治疗对乳腺癌无效,这强烈表明抗体的产生在抗肿瘤免疫反应中具有重要意义[22]。越来越多的证据表明,自然杀伤细胞(Natural killer, NK)细胞的频率、浸润和功能提高了GBM患者的生存率。但GBM细胞可表达独特的MHC-Ⅰ类分子与NK细胞表面的受体结合,从而抑制其功能[23],且NK细胞受体NCR2高表达可能与GBM预后较差有关[24]。T细胞在机体对恶性肿瘤的适应性免疫应答中发挥着重要作用。调节性T细胞(Regulatory T cells, Tregs)是一类独特的T细胞群,通过免疫抑制措施调节机体的整体免疫稳态。其中表达叉头状转录因子3(Forkhead Box P3, FOXP3)转录因子的Tregs尤为重要。该转录因子可下调NFAT和NFκB信号通路,从而下调IL2等重要效应细胞因子的表达。GBM预后较差与Tregs/T效应细胞比例较高有关[23]。
综上所述,内皮细胞基因相关标志物可在一定程度上反映GBM的恶性程度及预后。本研究通过qPCR实验进一步验证了GBM中内皮细胞基因相关标志物(DUSP6、STC1、VWA1 和TM4SF1)在GBM组织中的表达均较瘤周正常脑组织显著上调,并通过采用多基因构建风险评分,且验证集生存分析显示风险评分是GBM患者预后不佳的影响因素,进一步完善了GBM预后分析体系。未来需进一步进行大规模临床分析研究。
Competing interests: The authors declare that they have no competing interests.利益冲突声明:所有作者均声明不存在利益冲突。作者贡献:丁 硕:文章撰写与修改黄志刚、张洋:提出研究选题、设计论文框架侯丽珍、郭伟:整理数据尹高菲:调研整理文献房居高、钟琦:文章审阅 -
表 1 脂肪肉瘤发病部位及病理类型
Table 1 Location and pathological type of liposarcoma
Site Pathological type Total Well-differentiated Undifferentiated Myxoid Pleomorphic Myxo-pleomorphic Orbital cavity 4 7 6 1 0 18 Retropharyngeal space 2 1 0 0 0 3 Carotid space 1 1 0 2 0 4 Parapharyngeal space 0 1 0 0 0 1 Submental space 1 0 0 0 0 1 Larynx and hypopharynx 1 1 0 0 0 2 Temporal fossa 0 0 0 0 1 1 Total 9 11 6 3 1 30 表 2 脂肪肉瘤的治疗及病理类型
Table 2 Treatment and pathologic types of liposarcoma
Treatment Pathological type Total Well-
differentiatedUndiffer-
entiatedMyxoid Pleom-
orphicMyxo-
pleomorphicSurgery 6 7 0 3 0 16 Surgery+
radiation2 3 4 1 1 11 Surgery+
particle
therapy1 0 2 0 0 3 Total 9 10 6 4 1 30 表 3 接受不同治疗方式的脂肪肉瘤患者随访结果
Table 3 Follow-up of patients with liposarcoma treated with different modalities
Treatment No. Survival Recurrence Metastasis Surgery 16 15 6 0 Surgery+radiation 11 8 6 1 Surgery+particle 3 3 0 0 Total 30 2 12 1 表 4 脂肪肉瘤患者临床病理特征与预后的相关性
Table 4 Correlation between clinical pathological features and prognosis in patients with liposarcoma
Pathological
featuresn Percentage
(%)Control
rateP Disease-free
survivalP Gender 0.532 0.314 Male 12 40 38.1 60.6 Female 18 60 54.9 83.7 Site 0.651 0.940 Orbital set 18 60 45.6 77.2 Nonorbital
set12 40 47.6 68.8 Age (years) 0.208 0.073 ≥18-45 13 43 27.8 90.0 ≥45-60 10 33 42.9 68.6 ≥60 7 24 75.0 62.5 History of
benign
mesenchymal
origin tumors0.153 0.276 Yes 8 27 18.8 100.0 No 22 73 55.8 83.0 Tumor
size(cm)0.065 0.851 <3 14 47 61.0 92.0 3-5 6 20 60.0 75.0 ≥5 10 33 10.0 90.0 Ki67 0.014 0.169 <10% 16 54 56.4 100.0 10%-30% 7 23 56.6 87.5 ≥30%-50% 3 10 0.0 50.0 ≥50% 4 13 0.0 66.7 Surgery option 0.311 0.349 Extended resection 5 17 49.6 80.0 Tumor excision 25 83 30.0 89.7 Tumor capsule 0.945 0.12 Yes 17 57 50.3 100.0 No 13 43 46.3 78.8 Postoperative treatment 0.116 0.087 Yes 14 47 57.6 77.1 No 16 53 41.6 100.0 Pathological type <0.001 0.187 Well- differentiated 9 30 48.6 100.0 Dedifferentiated liposarcoma 10 33 68.8 100.0 Myxoid 6 20 40.0 66.7 Pleomorphic 4 13 0.0 66.7 Myxo- pleomorphic 1 4 0.0 0.0 -
[1] Sbaraglia M, Bellan E, Dei Tos AP. The 2020 WHO Classification of Soft Tissue Tumours: news and perspectives[J]. Pathologica, 2021, 113(2): 70-84.
[2] de Bree R, van der Waal I, de Bree E, et al. Management of adult soft tissue sarcomas of the head and neck[J]. Oral Oncol, 2010, 46(11): 786-790. doi: 10.1016/j.oraloncology.2010.09.001
[3] Ramu EM, Houdek MT, Isaac CE, et al. Management of soft-tissue sarcomas; treatment strategies, staging, and outcomes[J]. SICOT-J, 2017, 3: 20. doi: 10.1051/sicotj/2017010
[4] Mastrangelo G, Coindre JM, Ducimetière F, et al. Incidence of soft tissue sarcoma and beyond: a population-based prospective study in 3 European regions[J]. Cancer, 2012, 118(21): 5339-5348. doi: 10.1002/cncr.27555
[5] Luo P, Cai W, Yang L, et al. Retroperitoneal dedifferentiated liposarcoma: Analysis of 61 cases from a large institution[J]. J Cancer, 2018, 9(21): 3831-3838. doi: 10.7150/jca.25715
[6] Furukawa H, Tanemura M, Matsuda H, et al. A giant esophageal liposarcoma radically resected by the cervical approach: a case report[J]. Clin J Gastroenterol, 2022, 15(1): 71-76. doi: 10.1007/s12328-021-01550-z
[7] Barisella M, Giannini L, Piazza C. From head and neck lipoma to liposarcoma: a wide spectrum of differential diagnoses and their therapeutic implications[J]. Curr Opin Otolaryngol Head Neck Surg, 2020, 28(2): 136-143. doi: 10.1097/MOO.0000000000000608
[8] Lawrence W Jr, Anderson JR, Gehan EA, et al. Pretreatment TNM staging of childhood rhabdomyosarcoma: A report of the Intergroup Rhabdomyosarcoma Study Group. Children's Cancer Study Group. Pediatric Oncology Group[J]. Cancer, 1997, 80(6): 1165-1170. doi: 10.1002/(SICI)1097-0142(19970915)80:6<1165::AID-CNCR21>3.0.CO;2-5
[9] Choi JH, Ro JY. The 2020 WHO Classification of Tumors of Soft Tissue: Selected Changes and New Entities[J]. Adv Anat Pathol, 2021, 28(1): 44-58. doi: 10.1097/PAP.0000000000000284
[10] Yang T, Xiao L, Ren H. A well-differentiated liposarcoma of the prevertebral space: a case report[J]. Transl Cancer Res, 2021, 10(7): 3600-3604. doi: 10.21037/tcr-21-143
[11] Gerry D, Fox NF, Spruill LS, et al. Liposarcoma of the head and neck: Analysis of 318 cases with comparison to non–head and neck sites[J]. Head Neck, 2014, 36(3): 393-400. doi: 10.1002/hed.23311
[12] Nili F, Baghai F, Aghai A, et al. Well-differentiated liposarcoma of the floor of the mouth: Report of a rare case and review of the literature[J]. J Oral Maxillofac Pathol, 2016, 20(2): 312-315. doi: 10.4103/0973-029X.185984
[13] Mouret P. Liposarcoma of the hypopharynx. A case report and review of the literature[J]. Rev Laryngol Otol Rhinol (Bord), 1999, 120(1): 39-42.
[14] Breakey RW, Crowley TP, Anderson IB, et al. The surgical management of head and neck sarcoma: The Newcastle experience[J]. J Plast Reconstr Aesthet Surg, 2017, 70(1): 78-84. doi: 10.1016/j.bjps.2016.09.026
[15] Gritli S, Khamassi K, Lachkhem A, et al. Head and neck liposarcomas: A 32 years experience[J]. Auris Nasus Larynx, 2010, 37(3): 347-351. doi: 10.1016/j.anl.2009.08.003
[16] de Bree D, van der Valk P, Kuik DJ, et al. Prognostic factors in adult soft tissue sarcomas of the head and neck: A single-centre experience[J]. Oral Oncol, 2006, 42(7): 703-709. doi: 10.1016/j.oraloncology.2005.11.009