1. 引言

胃癌作为全球性的健康问题，作为最为常见的消化道恶性肿瘤之一，在2020年全球新增病例超过100万例，发病率排名第五，死亡率排名第四，其中东亚地区的发病率显著高于其他地区 [1] 。近年来胃镜检查的推广，使得胃癌的早诊早治率较前有明显提升，我国胃癌的发病率较前明显降低 [2] 。但尽管在内镜诊治、手术治疗、化学药物治疗、靶向药物治疗等方面有了很大进步，但当前胃癌患者的5年生存率仍较低 [3] ，在我国，胃癌的发病率和因胃癌死亡的人数仍位于各种肿瘤的前列。并且因为胃癌患者在肿瘤早期并没有显著的症状，当因为上消化道出血、胃部不适或幽门梗阻等症状就诊时，患者往往已处于胃癌进展期，部分患者已经失去手术根治机会。目前对于不能接受手术的患者进行抗肿瘤治疗时首选铂类及氟尿嘧啶联合序贯化疗方案进行抗肿瘤治疗 [4] ；对于特定靶点阳性的患者，可以同时联合一些分子靶向药物，如曲妥珠单抗、雷莫芦单抗等进行治疗。根据上一次欧洲医学肿瘤学会(ESMO) 2020年虚拟会议上发表的CheckMate 649研究结果，在晚期GC中，化疗和免疫疗法的结合似乎是最有前途的方法 [5] ，但即使是经过上述药物联合治疗，胃癌患者的总体中位生存时间仍然是难以超过12个月 [3] [6] 。近几年来免疫检查点抑制剂(immune checkpoint inhibitors, ICIs)应用为主的免疫治疗为胃癌患者带来了新的希望。目前已有的ICIs如派姆单抗、纳武单抗等在临床试验中表现出了较好的治疗效果 [7] [8] ，ICIs的应用在中国临床肿瘤医学会(CSCO)的指南中的地位越来越高，覆盖范围也越来越广，目前ICIs的应用已从三线治疗升级为一线治疗。但尽管当前晚期胃癌行ICIs联合化疗的客观缓解率(overall response rate, ORR)较高，仍有一部分患者效果欠佳，且治疗相关不良事件(treatment-related adverse event, TRAE)发生率较高。目前许多研究表明，不同患者肿瘤中免疫细胞浸润程度以及免疫检查点表达的水平差异可能是造成患者对于免疫治疗耐药性和不良反应的原因。因此，探究与肿瘤免疫细胞浸润与肿瘤细胞免疫检查点表达的新的治疗靶点以及准确的免疫浸润预后标志物，是提高胃癌患者免疫治疗效率所急需的。

Xα1型胶原蛋白(Collagen type X alpha 1, COL10A1)是属于胶原蛋白家族的一员，该基因编码的是参与软骨内成骨过程的一种小链状胶原蛋白 [9] ，过去该蛋白主要认为与软骨发育不良引起的Schimid型干骺端软骨发育异常(Schmid type metaphyseal chondrodysplasia, SMCD)以及骨关节炎有关 [10] [11] ，近期研究发现该蛋白是细胞外基质中的一种主要基质成分。根据当前研究，细胞外基质已被确定在肿瘤细胞的生长、分化、凋亡、转移等各个方面发挥重要作用 [12] 。近几年的研究显示，COL10A1在各种实体肿瘤中的高表达可能与肿瘤的血管生成有关 [13] ，而且COL10A1在结肠癌 [13] [14] 、乳腺癌 [15] [16] 、肺癌 [17] 等多种肿瘤中与患者的预后有关，是影响总生存率的独立威胁因素。但当前对于COL10A1在胃癌的作用尚不明确。

肿瘤微环境(tumor microenvironment, TME)指的是肿瘤的发生、生长及转移等过程中肿瘤细胞所处的内外环境，它不仅包括肿瘤所在的组织结构、功能和代谢，而且于肿瘤细胞自身的内在环境有关。肿瘤微环境的组成不仅仅与肿瘤细胞自身有关，同时与免疫细胞、基质细胞、内皮细胞以及细胞外基质分子等多种成分有关，他们在调节肿瘤发生、发展的过程中以及抗肿瘤治疗的效果等多种方面发挥着关键的作用 [18] 。依托于基因测序技术的发展，目前生物信息学分析已经用于众多肿瘤的研究之中，近年来肿瘤微环境的研究是当前肿瘤研究的热点之一，使用生物信息学分析来研究肿瘤微环境也已在多种肿瘤中得到应用，如乳腺癌 [19] [20] 、结直肠癌 [21] 、胆囊癌 [22] [23] 、肾细胞癌 [24] 等肿瘤的研究中。为明确COL10A1是否影响了肿瘤微环境的构成，本研究通过使用CIBERSORT研究目标基因差异表达与肿瘤微环境之间的相关性，通过比较高低表达组之间免疫细胞浸润情况，来比较两者之间肿瘤微环境差异；同时还可以通过使用Estimate算法来分别评价目标基因高低表达组之间得分，来比较两者肿瘤微环境的差异 [25] ：Estimate算法不仅将免疫细胞的数量纳入评价标准当中，同时也将基质细胞以及肿瘤细胞纯度等因素纳入对肿瘤微环境的评价之中。通过两种方法，来比较目标基因高低表达对于肿瘤微环境构成的影响。

本研究通过分析公共数据库TCGA、GTEx等数据库中数据以及青岛市市立医院行手术治疗的胃癌患者的肿瘤组织测序数据，探究了COL10A1在胃癌中可能参与的进程，COL10A1高表达与患者预后的相关性以及COL10A1与免疫细胞浸润相关的关系。

2. 资料和方法

2.1. 数据下载及处理

本研究中所用到的数据主要是由TCGA中的STAD项目中的可公开使用的转录组数据，该数据由GDC官网下载，下载时将数据类型选择为“transcriptomeprofiling”，数据种类选为“Gene Expression Quantification”，使用GDC官方下载工具下载RNA-Seq测序数据，共407例，其中测序数据中共有胃癌标本375例，正常癌旁组织32例。然后下载对应的临床数据，同样在官网将数据类型选择为“Clinical”，数据种类选为“Clinical Supplement”，同样通过GDC官网下载工具下载。下载完成后进行数据的初步处理，对于基因表达量部分，首先将RNA-Seq数据由FPKM (Fregments Per Kilobase per Million)格式转换为TPM (transcripts per millionreads)格式，然后根据TPM值将所有展示的基因表达量转换为经log2 (TPM + 1)转换后的数值；对于临床数据部分，将所下载数据中所需要用到的部分进行提取、合并，包括样本对应患者的基础信息，包括性别、年龄、生存状态等，以及疾病相关数据，包括肿瘤分期、生存时间等。GTEx中正常胃组织测序数据由UCSCXENA官网下载，该数据是经过Toil流程统一处理过的TCGA项目中的测序数据以及GTEx中的正常组织测序数据，已经过TPM转化，用于表达差异分析。

2.1.1. 差异表达分析

表达差异分析部分主要分为两个部分进行，一个是通过使用TIMER2.0在线工具，对COL10A1在多种肿瘤及肿瘤旁正常组织的表达量进行对比，大致了解COL10A1在各种肿瘤中的表达情况。作为验证，同时将使用R软件分析合并了GTEx与TCGA的数据，从而得出COL10A1在各肿瘤及正常组织中的表达差异。

然后具体到胃癌中，使用R软件中的limma包来比较正常癌旁组织与胃癌组织中的COL10A1表达量的不同，然后使用ggpubr包将结果进行可视化转换，得到差异表达图形。由于TCGA中的STAD项目里正常组织均为癌旁正常胃组织，因此将对应的癌旁组织及胃癌组织进行配对比较，分析COL10A1在胃癌组织与癌旁组织中的表达差异。

2.1.2. GO、KEGG富集分析

为探究COL10A1在胃癌的发生、发展中所起到的作用，将对分析与COL10A1表达差异相关的差异表达基因(Differentially Expression Genes)进行GO、KEGG富集分析：根据COL10A1的表达量的中位数值作为分界点，按照表达量高低分将胃癌样品分为高表达组及低表达组，然后比较高低表达组之间差异表达的基因，然后将这些基因与COL10A1一起进行富集分析。富集分析过程主要是通过使用cluster Profiler包进行，通过GO分析得到COL10A1可能相关的分子功能、生物过程及细胞成分信息，以校正后p值 < 0.05作为过滤条件，选取前十位进行可视化转换；通过KEGG分析得到COL10A1可能参与的代谢通路，以p值小于0.05作为过滤条件，选取通路并进行可视化转换。

2.1.3. GSEA富集分析

GSEA分析过程同样是建立在根据COL10A1表达量高低分组比较得出的差异表达基因的结果之上的 。应用clusterProfiler包，以由MSigDB数据库下载的C2基因集作为参考基因集，以p值小于0.05作为过滤条件，对根据COL10A1高低表达分成的高表达组与低表达组数据进行GSEA分析。得出在COL10A1高表达与低表达时活跃的前5条通路，并输出可视化图形结果。

2.1.4. 生存分析

为得到COL10A1的表达高低与胃癌患者总生存期(Overall Survival, OS)、无进展生存期(Progression-Free-Survival, PFS)的关系，通过使用R软件中的survival包对由TCGA中下载得到的测序数据以及临床数据进行合并整理，绘制生存曲线图。同时使用Kaplan-Meierplotter在线工具分析COL10A1表达高低与胃癌患者OS与PFS之间的关系，并绘制曲线图。

2.1.5. 生存预测模型构建

首先通过使用survival包对COL10A1进行COX单因素及多因素回归分析，明确COL10A1表达量的高低对于患者的生存时间是否有影响。如有确认了COL10A1表达高低对于患者生存有影响，可通过绘制ROC曲线分析COL10A1表达对于患者生存预测的能力。ROC曲线分析是用于评估一个因素预测能力的手段，是可以用于连续型变量分组的方法。在生存分析中，疾病状态和因素取值均会随时间发生变化。而标准的ROC曲线分析将个体的疾病状态和因素取值视作固定值，未将时间因素考虑在分析之中。因此使用时间依赖性ROC无疑是更好的选择。为进行时间依赖性ROC分析，可使用R软件中的timeROC包分析，分析了在不同时间截点时，COL10A1对患者的生存预测能力；进一步使用survival包、rms包处理基因表达数据及临床数据构建COX风险比例模型，同时将模型图形化，输出列线图(Nomogram)。

2.1.6. 对于肿瘤微环境及免疫细胞浸润影响

为分析COL10A1表达对于胃癌肿瘤微环境的影响，首先使用estimate包对由TCGA下载的肿瘤样本进行分析，得出各个样本的estimate评分，然后使用limma包将COL10A1高表达组与低表达组中的estimate评分作比较来评估COL10A1对于肿瘤微环境的影响。

为分析COL10A1表达对于免疫细胞浸润的影响，首先使用CIBERSORT官网提供的CIBERSORT包参照同样由官网下载的参考文件对所下载的TCGA胃癌样本进行分析处理，得出各个样本中免疫细胞浸润情况，根据p值 < 0.05为条件筛选准确度较高的样本，然后合并基因表达数据与免疫细胞浸润数据，分析COL10A1高表达组与低表达组之间各种免疫细胞浸润情况的差异，并进行可视化，输出合并表示的箱型图和浸润与COL10A1表达量相关的免疫细胞的相关性图。

2.1.7. 对于免疫治疗效果的影响

为分析COL10A1表达量高低对于接受免疫检查点抑制剂治疗的胃癌患者的治疗效果的影响，可以通过分析目标基因高低表达组之间免疫治疗效果的差异、目标基因表达量与常见免疫检查点相关基因表达量、目标基因表达量与样本肿瘤突变负荷(Tumor Mutation Burden, TMB)高低等三个方面来分析。

TCIA数据库是基于TCGA数据库中的数据建立的，其中归纳整理了TCGA中与免疫相关的数据，其中包括TCGA-STAD项目中样本对应的患者接受免疫治疗的效果评价。根据TCIA数据库中胃癌患者接受免疫治疗效果的评分数据，分析比较COL10A1高表达组与低表达组之间免疫治疗效果评分的差异，从而推测目标基因表达量对于免疫治疗效果的影响。

免疫检查点(Immune Checkpoint)，指程序性死亡受体及其配体。当前胃癌免疫治疗中免疫检查点阻断疗法的基本原理是基于免疫细胞T细胞的激活机制。程序性死亡受体表达在T细胞表面，其配体表达在肿瘤细胞和髓源性抑制细胞表面。程序性死亡受体与其配体结合可使T细胞衰竭而无法正常杀伤肿瘤细胞，肿瘤细胞便可逃脱宿主的免疫监视。因此，基于程序性死亡受体及其配体的免疫检查点阻断疗法通过抑制二者的结合，从而提高宿主免疫系统对肿瘤细胞的攻击性。因此通过分析目标基因与常见免疫检查点相关基因的表达之间的关系，可以推测目标基因与患者预后之间的关系。

肿瘤突变负荷(TMB)被定义为每百万碱基中被检测出的，体细胞基因编码错误、碱基替换、基因插入或缺失错误的总数。肿瘤细胞基因突变频率高，则细胞表面，所携带的肿瘤抗原就多，容易受到机体免疫系统的攻击，而被杀死。因此，检测恶性肿瘤的肿瘤突变负荷，就可以预测患者对免疫治疗的效果，目前已成为临床上，预测癌症免疫治疗疗效的、具有良好前景的全新生物标志物，在现有的研究中可证实免疫治疗效果与TMB水平有关。因此分析TCGA下载的胃癌患者的测序数据，得出各样本的TMB值，然后比较目标基因高低表达组之间的TMB水平差异，分析COL10A1表达量与TMB值之间的关系。

2.1.8. 对于化疗药物治疗效果的影响

分析常见化疗药物与目标基因表达量之间的关系主要通过R软件中的pRRophetic包实现。pRRophetic包是基于其制作者(明尼苏达大学Paul Geeleher, Nancy Cox, R. Stephanie Huang)所提出的算法：从基因表达数据预测表型，预测外部细胞系(CCLE)的药物敏感性，也可用于临床数据的预测。为分析COL10A1表达量高低对于胃癌组织对各种常见抗肿瘤药物敏感性之间的关系，通过使用R软件中的pRRophetic包处理分析由TCGA下载的胃癌患者测序数据，评估COL10A1高表达组与低表达组之间药物敏感性差异的分析。

3. 结果

3.1. COL10A1在多肿瘤及胃癌中表达差异分析

为明确COL10A1在多种肿瘤中的表达差异，首先在UCSC数据库下载经过Toil流程统一处理过的统一标准化的GTEx的测序数据，剔除了单个癌种中样本个数小于3个的癌种，最终获得了34个癌种的表达数据。使用非配对的Wilcoxon秩和检验进行差异显著性分析，我们在25种肿瘤中观察到了显著上调，在3种肿瘤中观察到了显著下调(图1)。作为验证，使用在线数据库TIMER2.0中分析COL10A1在多种肿瘤中的表达差异(图2)。两者结果基本相同。其中胃癌组织与正常胃组织表达具有明显差异。

为进一步分析COL10A1在胃癌组织和正常胃组织之间的表达差异，从TCGA数据库的STAD项目中下载了共407例样本的测序数据，其中共有375例样本为胃癌组织测序数据，32例为癌旁正常胃组织的测序数据。整理提取出各个样本中COL10A1的表达量，使用limma包对正常组与肿瘤组之间表达量差异(图3)，将这些正常组织与其对应的肿瘤组织配对比较，并在图中连线表达出来(图4)。将GTEx数据库中的正常胃组织测序数据加入正常组后再行差异分析比较两组中的COL10A1，结果如图所示(图5)。综上，正常组与肿瘤组两组中COL10A1表达差异明显，在胃癌组织中的表达明显高于正常组织。

3.2. GO、KEGG、GSEA富集分析结果

为将TCGA中下载的STAD相关测序数据以表达量中位值为筛选条件，分为高表达组和低表达组，对高表达组与低表达组之间的表达差异分析，筛选出可能与COL10A1的表达差异相关的基因(DEGs)。将数据分组后对高低表达组之间进行差异分析，其中以|logFC| > 1且校正后p值 < 0.05作为筛选条件，共得到差异基因，其中上调的基因共242个、下调的基因共200个，将其中最为差异最为明显的上调基因及下调基因各50个的表达差异分析结果可视化，绘制成环形图(图6)及热图(图7)。

将所获得的基因使用R软件中的cluster Profiler包进行GO分析，其中经GO富集分析后共鉴定出可能参与的284项基因集，其中生物过程(biological process, BP)基因集共213个，细胞组分(cellular component, CC)基因集共21个，生物过程(biological process, BP)基因集50个，分别选取三类中可信度最高的10中基因集进行可视化转换，绘制环形图(图8)及气泡图(图9)。如图所示，结果提示主要与细胞外基质、细胞外机构组织、外部封装结构组织等细胞外结构的构成以及骨化、骨骼系统开发、内胚层的发育等功能。

使用clusterProfiler包对DEGs进行KEGG富集分析，依据筛选条件共选出可信度最高的18个通路信息，将富集分析结果进行可视化转换，绘制条形图(图10)。根据结果所示可信度排名较前的几条通路(p < 0.01)主要与以下几个功能有关：1、细胞间信息传递；2、IL-17信号通路；3、TGF-β信号通路；4、PI3K-Akt信号通路。

使用clusterProfiler包，以Msigdb数据库下载的C2集中的c2.cp.kegg.v7.4.symbols基因集作为参考基因集，对筛选出的DEGs进行GSEA分析，共得到可能参与的通路共37条，其中满足筛选条件的通路共14条，选择其中可信度最高的5条通路进行可视化转换，结果如图(图11)。如图所示随目标基因高表达激活的通路是与细胞间信息传递功能相关通路与TGF-β信号通路，受到抑制的为嗅觉传导相关通路。

3.3. 生存分析结果

将胃癌样本的测序数据与临床数据整合，去除掉测序数据中正常组织部分，然后根据COL10A1的表达量中位值，将样本分为高表达组及低表达组。然后使用limma、survival和survminer包来分析整合数据：比较高、低表达组之间的生存差异得到显著性p值并绘制总生存期(Overall Survival, OS)曲线(图12)。以同样流程得出无进展生存期(Progression-Free-Survival, PFS) (图13)。同时作为验证，使用Kaplan-Meierplotter数据库中的在线分析工具，绘制OS曲线(图14)及PFS曲线(由于K-M plotter数据库没有PFS选项，因此使用FirstProgression曲线代替) (图15)。如图TCGA数据库STAD数据所得出的OS曲线图和K-Mplotter数据库所得到的OS曲线图所示，高、低表达组之间生存曲线之间p值均小于0.05，而PFS图中p值均大于0.05。说明COL10A1表达高低对于胃癌患者的总生存时间有影响，高表达患者总生存时间较短，而COL10A1表达量高低表达组之间无进展生存期长短无差异，对于患者的无进展生存期无影响。

3.4. 生存预测模型构建

为探索COL10A1表达量与患者的生存状态之间的关系，首先处理由TCGA下载的临床数据，将数据根据目标基因的表达量高低将样本分为两组，分别分析两组之间各临床数据之间的差别，结果如图所示(图16、图17)，两组样本之间年龄、性别、肿瘤Grade分期、M和N分期等信息无明显差别，Stage分期有明显差异，T分期有显著差异。

使用survival包对临床数据进行分析，从而明确COL10A1的表达量以及年龄、TNM分期、Grade分级和Stage分期等其他临床数据对生存的影响。首先对上述的各种因素行单因素COX分析，将结果绘制成森林图表示(图18)，如图所示COL10A1表达量高低在单因素COX分析中p值小于0.05，HR大于1；同样步骤将上述各影响因素行多因素COX分析(图19)，如图所示COL10A1表达量高低在分析中的p值仍小于0.05且HR大于1。综上可初步判定COL10A1的表达量高低与患者的预后有关。

使用timeROC包分析临床数据绘制时间依赖性ROC曲线，评估1年、3年、5年时COL10A1表达量对于患者生存预测能力，曲线如图所示(图20)，5年时曲线下面积最大，为0.713，而1年、3年时曲线下面积分别为0.570、0.584，提示COL10A1在1年、3年时具有一定的生存预测能力，而5年时生存预测能力最佳。

进一步使用R软件中survival、rms包使用COL10A1、年龄、性别、肿瘤分期等因素行多因素COX回归分析，构建COX比例风险模型，使用R软件的regplot包将该模型以列线图(nomogram) (图21)的形式表现出来，同时输出预测模型的矫正曲线。该图的内容主要有：1) 预测模型中的变量(每一个变量上都有其可取值范围)；2) 得分(每个影响因素在不同的取值下都具有其对应的单项得分，以及所有变量对应的总得分)；3) 预测概率(总得分对应到1年、3年、5年上的生存概率)。如图所示，假如患者总分为475分，则该患者1年生存率为75.5%、3年生存率为41.2%、5年生存率为27.7%。该模型共纳入了年龄、性别、TNM分期、Grade分期、Stage分期等因素，其中年龄的影响度较高，其次为COL10A1表达量。与列线图同时得出的校准曲线(calibrationcurve) (图22)，如图所示三条线分别为1年、3年、5年预测值的曲线，基本沿对角线分布，预测效果尚可。

3.5. COL10A1对于胃癌肿瘤微环境的影响

使用estimate包对由TCGA下载的肿瘤样本进行分析，得出每个样本的基质得分(Stromal Score)和免疫得分(Immune Score)，由两者联合得出该样本的Estimate评分(ESTIMATE Score)，并得出每个样本的肿瘤纯度。比较高、低表达组之间基质得分、免疫得分和Estimate评分的差异，如图所示，将差异分析结果绘制成小提琴图(图23)，如图所示，基质得分和Estimate评分在高、低表达组之间具有显著的差异，p值小于0.001；免疫评分在高、低表达组之间也有明显差异，其p值小于0.05。综上，COL10A1的高表达者胃癌组织的Estimate评分较高，也意味着该患者胃癌组织纯度更低。

3.6. COL10A1对于胃癌组织免疫细胞浸润影响

使用R软件中由Cibersort官方提供的Cibersort包对由TCGA下载的胃癌测序数据进行整理、分析。根据目标基因COL10A1的表达量高低将该样本分为高、低表达组，然后使用limma包比较高、低表达组之间免疫细胞的差异，并将结果转化为箱型图(图24)。对目标基因表达量以及免疫细胞表达量之前做相关性分析，结果以Lolipop图(图25)表达出来，如图所示，如图所示，M2型巨噬细胞、B记忆细胞在两组中的表达差异显著，其p值小于0.001；激活的树突状细胞在两组中表达差异较明显，p值小于0.01；M0、M1型巨噬细胞，调节性T细胞在两组中有一定差异。可推断COL10A1的高表达可能与免疫细胞大量浸润有关。

3.7. COL10A1对于胃癌患者免疫治疗效果的影响

将TCIA数据库下载的免疫治疗打分文件与下载的TCGA中胃癌测序数据整合，根据样本ID将TCIA中关于免疫治疗效果打分文件与整理后的测序数据整合起来。根据COL10A1的表达量分为高、低表达两组，使用limma包将高、低表达组之间免疫治疗评分进行比较，评估高低表达组在CTLA4阳性、PD-1阴性，CTLA4阴性、PD-1阳性，CTLA4阴性、PD-1阴性，CTLA4阳性、PD-1阳性等四种情况下免疫治疗效果，并使用ggpurb包将最终差异比较结果转换为小提琴图表达出来(图26)，如图所示，在CTLA4阳性、PD-1阴性，以及CTLA-4阴性、PD-1阴性这两组中，高、低表达组之间的免疫治疗评分具有显著差距。

为探究目标基因表达量与常见免疫检查点相关基因之间的关系，提取肿瘤样本测序数据中的免疫检查点相关基因的表达量，使用皮尔森相关系数对这些免疫检查点基因表达量与目标基因表达量进行相关性检验，相关性检验过滤条件设置为p值 < 0.001，将分析结果转化为热力图(图27)，如图所示，目标基因表达与HAVCR2、NRP1、CD70、TNFSF4、PDCD1LG2、CD86、LAIR1、TNFSF9、CD276、TNFSF18、TNFRSF4、CD200等多种免疫检查点之间为正调控关系。

4. 讨论

在本研究中，为了确定COL10A1在胃癌中的作用，首先是通过TIMER2.0及联合TCGA联合GTEx数据在多种瘤中进行差异表达分析，根据分析结果在除了在膀胱癌、乳腺癌、头颈部鳞癌、肾乳头状细胞癌、肺鳞癌、肺腺癌、前列腺癌等肿瘤中肿瘤组织与正常组织中的目标基因表达量有显著差异，同时在胃癌、直肠癌、胆管癌、结肠癌、食管癌等消化道相关的肿瘤中，肿瘤与正常组织中具有显著的表达差异，其中结直肠癌 [14] 、食管癌、乳腺癌 [15] 、肺癌 [17] [26] 中的差异已有研究证明。通过绘制根据目标基因表达量高低的分组中的胃癌患者的总生存期和无进展生存期曲线的差异，可以得知目标基因高表达组的总生存期相较于低表达组的生存期更低，而在无进展生存期中高低表达组之间无明显差异。进一步通过单因素、多因素COX回归分析发现COL10A1的高表达、年龄、高Stage分期等是胃癌患者的预后独立因素。通过汇总上述预后相关因素，构建了对于患者预后的预测模型，通过绘制时间依赖性的ROC曲线及矫正曲线进行评估，该模型对于胃癌患者生存具有一定的预测效果。同时根据已有研究显示，胃癌患者血液中的COL10A1表达量亦有显著差距，且与胃癌患者的预后有关，具有成为胃癌预后的标志物潜力 [27] ，该结果与本研究结果可以互相支持验证。

对于免疫治疗，当前已证实的是免疫微环境特征可以作为评估肿瘤细胞对于免疫治疗的反应及影响的可靠的生物标志物之一。Estimate算法通过对基因表达数据分析，来预测恶性肿瘤组织中的间质细胞和免疫细胞的含量，推测肿瘤纯度，从而评估肿瘤微环境情况。本研究中通过使用R软件中的estimate包为胃癌样本进行打分，然后比较目标基因高低表达组之间的差异，结果显示高表达组中的样本三项得分均显著高于低表达组，预示着目标基因的表达与肿瘤微环境中的基质和免疫细胞呈正相关。进一步使用Cibersort算法进一步探究COL10A1表达与各种免疫细胞的浸润水平之间的关系，结果提示COL10A1高表达组具有更多的免疫细胞浸润，其中M0、M1、M2巨噬细胞，未活化的树突状细胞的浸润水品成正相关。此外，结果表明COL10A1与HAVCR2、NRP1、CD70、TNFSF、PDCD1LG2、CD86、LAIR1、TNFSF9、CD276、TNFSF18、TNFRSF4及CD200常见的免疫检查点基因的表达量成正相关。同时通过分析目标基因COL10A1高低表达组之间患者对于免疫治疗的效果评分，表明高表达组的患者应用免疫治疗效果较差。目前根据现有研究表明，细胞外基质与功能失调的免疫逃逸与免疫治疗的失败密切相关。综上，我们猜测COL10A1的高表达可能通过影响肿瘤微环境的构成及免疫检查点的表达导致胃癌患者接受免疫治疗效果较差，提示COL10A1具有成为潜在的免疫治疗靶点以及成为免疫检查点阻断治疗的预测指标。

通过对目标基因高低表达组之间的差异表达基因的GO、KEGG和GSEA分析提示COL10A1可能参与了TGF-β信号通路 [28] 、PI3K-Akt信号通路 [29] 。目前两种信号通路已经经过长时间研究：TGF-β信号通路在成熟有机体和发育中的胚胎中都参与了许多细胞过程，这些过程包括细胞生长，细胞分化，细胞凋亡，细胞动态平衡等其它细胞功能，目前大量研究表明该通路在肿瘤的进展和转移中起到关键作用；PI3K-AKT是细胞中一个经典的信号通路，它是一条与磷脂酰肌醇有关的信号通路，也是RTK介导衍生的信号通路。在胰岛素的刺激以及细胞生存当中发挥着重要作用；AKT作为一种原癌基因，在多种细胞过程中发挥关键作用，如葡萄糖代谢，细胞凋亡，细胞增殖，转录和细胞迁移；PI3K/AKT信号通路失调见于多种人类疾病，包括癌症、糖尿病、心血管疾病和神经疾病。这些潜在的信号通路为后续继续对目标基因COL10A1的基质研究提供了参考。

综上所述，在本研究中通过分析目标基因COL10A1在胃癌细胞中的高表达，并与胃癌患者的预后较差存在显著相关，通过该基因以及年龄、Stage分期等其他因素所构成的临床预测模型可较好的预测患者1年、3年、5年的总生存期。对于目标基因可能参与的机制的探索，该基因在胃癌中与免疫细胞的浸润和免疫检查点的表达具有显著性相关，并且可能通过TGF-β信号通路、PI3K-Akt信号通路等促进胃癌的发生与发展。结合分析表明目标基因COL10A1具有作为胃癌患者免疫治疗靶点的潜在能力。

5. 结论

COL10A1在胃癌中高表达，当COL10A1表达上调时胃癌患者的总体生存期显著缩短，但对无进展生存期无明显影响，提示其与患者的预后不良有关。COL10A1与肿瘤微环境的构成有关，与免疫细胞浸润和免疫检查点的表达存在显著性相关。通过分析TCIA数据库中数据，COL10A1高表达会使部分患者的免疫治疗效果差，提示COL10A1可能是一个潜在的免疫治疗靶点。

NOTES

^*通讯作者。

参考文献