Analysis on the Expression of Basement Membrane Protein in Endometrial Carcinoma and Its Relationship with Immune Invasion
-
摘要:
目的 探讨基底膜蛋白(Hub基因)与子宫内膜癌的预后、免疫浸润及泛癌的关系,建立临床预后模型。 方法 利用癌症基因图谱(the cancer genome atlas,TCGA)子宫内膜癌数据库中的子宫内膜癌组织及癌旁组织表达谱和临床信息,筛选出差异表达基因,对LncRNA共表达的mRNA进行GO富集分析及KEGG信号通路分析,将子宫内膜癌与癌旁组织的差异基因与基底膜蛋白基因进行交集,将筛选出的差异表达基因与生存状态、生存时间进行合并通过Lasso-cox回归分析筛选出Hub基因,利用多因素Cox回归分析建立预后模型,进一步进行泛癌分析以及免疫浸润相关性分析。 结果 确定了6个与子宫内膜癌患者预后相关的Hub基因,ADAMTS5、EVA1C、THBS4、CTSD、ITGAV和LAMA1,发现患者随着风险评分增加,患者生存率明显下降。进行泛癌分析发现这6个基因在大多数癌症类型中显著差异(P < 0.05),且在多形性胶质母细胞瘤(glioblastoma multiform,GBMLGG)、脑低级别胶质瘤(low grade glioma,LGG)、急性髓细胞样白血病(acute myeloid leukemia,LAML)、葡萄膜黑色素瘤(uveal melanom,UVM)、肾上腺皮质癌(adrenocortical carcinom,ACC)等癌症中高表达预后差。并探究了这些基因在肿瘤免疫治疗中的潜在作用发现与Th17 cells、Th2 cells 、NK CD56bright cells等免疫细胞存在显著负相关(P < 0.01),与Tcm、iDC、Eosinophils、aDC等免疫细胞存在显著正相关(P < 0.01)。 结论 基底膜蛋白基因在子宫内膜癌的诊断和预后方面有较高的临床价值,可作为子宫内膜癌患者的预后标志物和潜在治疗靶点。 Abstract:Objective Immune infiltration, as well as their implications across various cancers, and to establish a clinical prognostic model based on the findings. Methods The expression profile and clinical information of endometrial cancer tissue and adjacent tissues in TCGA endometrial cancer database were used, the differential expression genes were screened out for analysis, and the mRNA co-expressed by LncRNA was analyzed by GO enrichment and KEGG signaling pathway.The differential genes between endometrial cancer and adjacent tissues were intersected with the basement membrane protein genes. The selected differentially expressed genes were combined with survival status and survival time to screen out Hub genes by Lasso-cox regression analysis. Multivariate Cox regression analysis was used to establish a prognostic model, and pan-cancer analysis and immunoinfiltration correlation analysis were further performed. Results Six basal membrane protein Hub genes, ADAMTS5, EVA1C, THBS4, CTSD, ITGAV and LAMA1, were identified as associated with the prognosis of patients with endometrial cancer, and found that the survival rate of patients decreased significantly with the increase of risk score. Pan-cancer analysis found that these 6 genes were significantly different in most cancer types, and high expression in GBMLGG (glioma), LGG (low-grade glioma of the brain), LAML (acute myeloid leukemia), UVM (Uveal melanoma), ACC (adrenal cortical cancer) and other cancers had poor prognosis. The potential role of these genes in tumor immunotherapy was also explored, and significant negative correlation was found with Th17 cells, Th2 cells, NK CD56 bright cells and other immune cells (P < 0.01), and significant positive correlation was found with Tcm, iDC, Eosinophils, aDC and other immune cells (P < 0.01). Conclusion Basal membrane protein gene has high clinical value in the diagnosis and prognosis of endometrial cancer, and can be used as a prognostic marker and potential therapeutic target for patients with endometrial cancer. -
Key words:
- Basement membrane protein gene /
- Endometrial carcinoma /
- Prognosis /
- Immune infiltration
-
子宫内膜癌(endometrial carcinoma,EC)是1种源自子宫内膜腺体的恶性肿瘤,属于女性生殖器3大恶性肿瘤之一[1−2],其发病率约占女性癌症的7%[3],占生殖道恶性肿瘤的20%~30%[4]。其发病机制尚不清楚,可能与长期雌激素刺激、肥胖、高血压、糖尿病、不孕或绝经、遗传等因素有关[5−7]。近年来,子宫内膜癌的患病率有所增加,同时随着经济状况改善,人们的平均预期寿命也较前增加了。尽管子宫内膜癌患者的治疗效果似乎很乐观,但用于检测子宫内膜癌生物标志物验证仍然存在挑战[8−9]。基底膜(basallamina,BM)是细胞外基质特化而形成的1种柔软、坚韧的网膜结构,研究表明,BM相关成分的遗传缺陷可能导致患者的疾病表型,表现在视网膜[10−11]、肾脏[12]、血管[13]、骨骼系统[14]和肌肉[15]等各个方面。基底膜蛋白是免疫疾病中自身抗体的靶标,基底膜蛋白表达和周转缺陷是癌症、糖尿病和纤维化的关键致病方面。作为1种屏障结构,BM在癌细胞迁移和浸润过程中起着重要作用,防止其进入血液和淋巴管系统[16]。癌细胞在传播时,BM会限制其对血液和淋巴管的入侵(腔内)和迁移(外渗)方式,这一过程与90%的癌症死亡相关[17]。细胞通过BMS进入周围区域后,5年存活率显著降低[18]。BM的作用在于防止癌细胞侵入、迁移和蔓延,对于癌症的增长和传播具有至关重要的意义[19−20]。
本研究利用公共数据库TCGA中提供的子宫内膜癌及癌旁组织RNA-seq数据,筛选潜在的子宫内膜癌和癌旁组织的相关显著差异基因,随后进行生存预后分析,为临床医生治疗子宫内膜癌提供新的参考。
1. 资料与方法
1.1 资源与处理
研究团队从TCGA数据库中获取了子宫内膜癌(uterine corpus endometrial carcinoma,UCEC)及其癌旁组织的RNA表达谱和临床信息。在去除重复样本后,他们收集了544个子宫内膜癌基因样本和35个癌旁组织样本。将以Fregments Per Kilibase per Million(FRKM)格式存储的RNAseq转换为(transcripts per million reads,TPM)格式,并进行了对数转换,以便进行后续分析。使用R语言对基因表达数据进行处理,将基因探针ID映射为对应的基因符号,生成含有标准国际基因名称的矩阵文件,清除无对应基因名和探针对应多个基因的数据,对多个探针映射到基因的数据取平均值 [21]。从BMBASE数据库[22](https://bmbase.manchester.ac.uk)提取了222个基底膜基因的表达。
1.2 研究方法
1.2.1 差异表达基因的筛选
使用R语言包中的“edgeR”对子宫内膜癌组织和正常组织的CRGs表达数据进行差异性分析,筛选条件为P < 0.05 且| log2 差异倍数( fold change,FC) | 大于1.5。
1.2.2 差异表达基因的 GO和KEGG 富集分析
使用R软件包org.Hs.eg.db(version 3.1.0)中的基因的GO注释及从KEGG Pathway注释,将基因映射到注释集合中,使用R软件包ClusterProfiler (version 3.14.3)进行富集分析并得到结果。设置最小基因集为5,最大基因集为
5000 ,P < 0.05和 FDR<0.25被认为具有统计学意义。利用韦恩图的方式,将子宫内膜癌差异基因与BM的222个基底膜基因进行比较。1.2.3 建立与评估UCEC预后模型
合并筛选出的差异表达基因与子宫内膜癌患者的生存状态和生存时间。经过Lasso-cox回归分析以及多因素Cox回归分析,确定与子宫内膜癌预后相关的差异基因以及计算它们的回归系数,建立子宫内膜癌预后风险预测模型。依据风险评分的中位数,将子宫内膜癌患者分成2组,分别为高危及低危组,并应用“survminer”R包对2组的总生存期(overall survival,OS)进行生存分析。另外,评估子宫内膜癌患者预后的风险评分预测能力,绘制了C-index曲线和ROC曲线下面积(area under curve,AUC)。
1.2.4 泛癌分析
从TCGA数据库中下载了经统一标准化的泛癌数据集:TCGA TARGET GTEx (PANCAN,N=
19131 ,G=60499 ),从中提取了基因在各个样本中的表达数据,筛选了样本来源为:正常实体组织、原发性实体肿瘤、原发性肿瘤、正常组织、原发性血癌-骨髓、原发性血癌-外周血的样本,更进一步的对每1个表达值进行了log2(x+0.001)变换,还剔除了单个癌种中样本个数小于3个的癌种,最终获得了34个癌种的表达数据。1.2.5 免疫浸润相关性分析
采用GSVA软件包中的ssGSEA算法,根据Immunity杂志提供的免疫细胞标志物,分析TCGA数据库中TCGA-UCEC (子宫内膜癌)项目的TPM格式数据和临床数据,评估相应数据的免疫浸润水平,P < 0.05为差异有统计学意义。
1.2.6 统计学处理
使用R语言制作图表和统计学分析(3.6.3版)。采用Wilcoxon秩和检验分析非配对样本中ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1的表达情况。进行单因素logistic回归分析或单因素/多因素cox回归分析对显著差异进行统计学检验。P < 0.05为差异有统计学意义。
2. 结果
2.1 筛选差异表达基因
从TCGA-UCEC资料库中提取相关数据比较子宫内膜癌基因样本和癌旁组织样本,结果显示
5882 个差异基因,其中上调基因2786 个,下调基因3096 个,绘制火山图,见图1A和图1B。2.2 KEGG、GO富集分析不同基因
对差异基因进行GO富集分析,结果显示这些基因在生物过程中主要集中在系统开发、细胞外区域、泡 、动物器官发育、胞外区部分、解剖结构形态发生、细胞外小泡、细胞外细胞器、细胞外外来体、循环系统发展等,见图1C。KEGG通路分析提示,差异基因主要富集于以下途径中,见图1D,癌症的途径、癌症中的蛋白聚糖、焦点粘连、轴突导向、胃癌、血管平滑肌收缩、细胞凋亡、细胞周期、AGE-RAGE信号通路与糖尿病并发症、ECM-受体相互作用等。
2.3 与基底膜蛋白基因交集
从基底膜数据系统提取222个基底膜基因的表达,与子宫内膜癌及癌旁差异基因进行交集,得出108个交集基因,见图2A。
2.4 建立预后模型
对108个交集基因整合生存时间、生存状态和基因表达数据,利用lasso-cox方法进行回归分析,设置Lambda值为
0.0565552971333623 ,最终获得了11个基因,见图2B和2C。进一步进行多因素Cox回归分析结果,检出回归系数和显著性明显的6个基因:ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1,见图2D。2.5 风险评分模型的预后
对于544位子宫内膜癌患者,利用风险评分公式进行计算,得出其相应的风险值,并按照中位值将患者分为高风险和低风险2组。随着风险评分的不断增加,病人的生存率呈现显著下降的风险曲线和散点图,见图3A。通过对风险评分和临床特征进行计算,可以预测子宫内膜癌患者365 d(1 a)、
1095 d(3 a)、1825 d(5 a)生存率的AUC值,研究结果表明,风险评分可以有效预测子宫内膜癌患者1 a、3 a、5 a的生存率的AUC值分别为0.7、0.72、0.76,对临床实践有重要意义,见图3B。K-M生存曲线比较2组之间的OS(总生存期),发现子宫内膜癌患者高风险组患者的OS低于低风险组(P < 0.001),见图3C。2.6 泛癌分析
利用TCGA数据库,研究了ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1在34种癌症患者和正常健康对照人群的表达水平的差异,ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1在大多数癌症类型中显著差异(P < 0.05),见图4。其中ADAMTS5在12种肿瘤中观察到了显著上调(P < 0.05),在19种肿瘤中观察到了显著下调(P < 0.05),如UCEC等。EVA1C在9种肿瘤中观察到了显著上调(P < 0.05),在24种肿瘤中观察到了显著下调(P < 0.05),如UCEC。THBS4 15种肿瘤中观察到了显著上调(P < 0.05),在15种肿瘤中观察到了显著下调(P < 0.05),如UCEC。CTSD在26种肿瘤中观察到了显著上调(P < 0.05),如UCEC,在5种肿瘤中观察到了显著下调(P < 0.05)。ITGAV在18种肿瘤中观察到了显著上调(P < 0.05),9种肿瘤中观察到了显著下调(P < 0.05),如UCEC。LAMA1在13种肿瘤中观察到了显著上调(P < 0.05),如UCEC,在13种肿瘤中观察到了显著下调(P < 0.05)。ADAMTS5、EVA1C、THBS4、ITGAV在USEC中显著下调(P < 0.05),CTSD、LAMA1在UCEC中显著上调(P < 0.05)。分析基因表达在肿瘤中的预后关系,结果表示ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1在GBMLGG、LGG、LAML、UVM、ACC等癌症中高表达预后差,见图5和图6。
2.7 免疫浸润分析
为了更准确地探究ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1表达与免疫浸润的相关性,采用ssGSEA算法分析了SERPINE1表达与24个免疫细胞表达的相关性,见图7。结果表明,ADAMTS5表达与Tcm、Tgd、NK cells、T helper cells、Mast cells中存在显著的正相关(P < 0.01),与Th17 cells、NKCD56bright cells、Cytotoxic cells、TReg 、T cells、iDC存在显著负相关。EVA1C表达与NK cells、Neutrophils、Eosinophils、NK CD56bright cells 、Mast cells 、iDC存在显著正相关,与Th2 cells存在显著负相关。THBS4与NK cells、 Eosinophils存在显著正相关,与Th17 cells、 Cytotoxic cells、pDC、Treg存在显著负相关。CTSD与Treg、T cells、NK CD56dim cells、B cells、iDC、 DC、Neutrophils、Cytotoxic cells、Th1 cells、Macrophages、aDC存在显著正相关,与Tcm、T helper cells存在显著负相关。ITGAV表达与Macrophages、Tgd、Th2 cells、T helper cells、Tcm存在显著正相关,与NK CD56bright cells、pDC、 TReg、Th17 cells、NK cells、Cytotoxic cells存在显著负相关。LAMA1表达与NK CD56dim cells、 DC 、Cytotoxic cells 、pDC、T cells、CD8 T cells、TReg、TFH、Neutrophils、Tem存在显著负相关,未见显著正相关。
3. 讨论
EC作为常见的妇科恶性肿瘤[23],其与基底膜蛋白之间的关系是1个重要的研究领域。基底膜蛋白是参与细胞极性、信号传导、细胞外基质组织和细胞间相互作用等生物学过程的1类在基底膜结构中发挥重要作用的蛋白质。据研究显示,基底膜蛋白在调控细胞增殖、迁移和侵袭中发挥着重要作用[24]。基底膜蛋白的异常表达可能与肿瘤的侵袭性和转移性有关[25],也可能与肿瘤的恶性程度和预后密切相关。此外,通过调节肿瘤微环境中的细胞相互作用和信号传导,基底膜蛋白可能在子宫内膜癌的发生和发展过程中扮演重要角色,基底膜蛋白的异常表达和功能常常与肿瘤浸润、转移以及治疗反应息息相关。研究子宫内膜癌与基底膜蛋白之间的关系有助于深入解析子宫内膜癌的发病机制,并为开发相关的治疗策略提供重要的理论基础。这次研究基底膜蛋白基因在子宫内膜癌组织的基因表达水平及预后,也研究了其在不同癌症类型中显著差异表达。结果提示基底膜蛋白中的ADAMTS5、EVA1C、THBS4、CTSD、ITGAV 和 LAMA1这些基因在子宫内膜癌中显著表达,且对该癌症的预后有重要影响。这些基因在多种癌症类型中显示出表达差异,这种差异可能与它们在肿瘤的生长、侵袭和转移中的作用相关。需要进一步深入研究这些基因的功能以及它们在癌症中的作用机制。
对基底膜蛋白基因中的ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1基因表达与免疫细胞浸润之间的相关性也进行了深入的探讨。从结果来看,这些基因的表达与多种免疫细胞类型之间存在着显著的相关性,与Th17 cells、Th2 cells、NK CD56bright cells等免疫细胞存在显著负相关,与Tcm、iDC、Eosinophils、aDC等免疫细胞存在显著正相关,这些结果可能为癌症免疫治疗和免疫监测提供了重要线索,为更深入地理解癌症的免疫微环境提供了有益的信息。基底膜蛋白基因的表达与多种免疫细胞类型的正负相关性的发现,对于进一步深入深入探讨它们在癌症免疫治疗中的潜在影响至关重要。这些关联性的发现可能为制定个体化的免疫治疗策略提供了线索,也为相关免疫监测方法的开发提供了参考。
综上所述,研究结果初步表明ADAMTS5、EVA1C、THBS4、CTSD、ITGAV、LAMA1与子宫内膜癌的发生和进展密切相关,对于深入理解子宫内膜癌的发病机制、探讨靶向治疗的新策略以及预测患者预后具有重要意义,但鉴于研究方法仅限于生物信息学分析,仍有许多问题有待进一步研究,希望通过持续的深入研究,为临床治疗提供更有效的依据。
-
-
[1] Wang N,Zhang J,Fan X,et al. Identification of risk factors for the prognosis of Chinese patients with endometrial carcinoma[J]. Medicine (Baltimore),2021,100(38):e27305. [2] Li Y T,Yang S T,Wang P H. Molecular pathology and prognosis of endometrial cancer[J]. Taiwan J Obstet Gynecol,2022,61(6):921-922. doi: 10.1016/j.tjog.2022.06.011 [3] Maheshwari E,Nougaret S,Stein E B,et al. Update on MRI in evaluation and treatment of endometrial cancer[J]. Radiographics,2022,42(7):2112-2130. [4] Liu Y,Li H,Zhao C,et al. MicroRNA-101 inhibits angiogenesis via COX-2 in endometrial carcinoma[J]. Mol Cell Biochem,2018,448(2):61-69. doi: 10.1007/s11010-018-3313-0 [5] Nees L K,Heublein S,Steinmacher S,et al. Endometrial hyperplasia as a risk factor of endometrial cancer[J]. Arch Gynecol Obstet,2022,306(2):407-421. [6] Crosbie E J,Kitson S J,Mcalpine J N,et al. Endometrial cancer[J]. Lancet,2022,399(10333):1412-1428. doi: 10.1016/S0140-6736(22)00323-3 [7] Raglan O,Kalliala I,Markozannes G,et al. Risk factors for endometrial cancer: An umbrella review of the literature[J]. Int J Cancer,2019,145(7):1719-1730. doi: 10.1002/ijc.31961 [8] Marin-jimenez J A,Garcia-mulero S,Matias-guiu X,et al. Facts and hopes in immunotherapy of endometrial cancer[J]. Clin Cancer Res,2022,28(22):4849-4860. doi: 10.1158/1078-0432.CCR-21-1564 [9] Karpel H,Slomovitz B,Coleman R L,et al. Biomarker-driven therapy in endometrial cancer[J]. Int J Gynecol Cancer,2023,33(3):343-350. doi: 10.1136/ijgc-2022-003676 [10] Wilson S E. Defective perlecan-associated basement membrane regeneration and altered modulation of transforming growth factor beta in corneal fibrosis[J]. Cell Mol Life Sci,2022,79(3):144-159. [11] Wilson S E. Fibrosis Is a Basement membrane-related disease in the cornea: Injury and defective regeneration of basement membranes may underlie fibrosis in other organs[J]. Cells,2022,11(2):309-330. doi: 10.3390/cells11020309 [12] Kuang H,Liu J,Jia X Y,et al. Autoimmunity in anti-glomerular basement Membrane Disease: A Review of Mechanisms and Prospects for Immunotherapy[J]. Am J Kidney Dis,2023,81(1):90-99. doi: 10.1053/j.ajkd.2022.07.006 [13] Mukwaya A,Jensen L,Lagali N. Relapse of pathological angiogenesis: Functional role of the basement membrane and potential treatment strategies[J]. Exp Mol Med,2021,53(2):189-201. [14] Kanazawa Y,Nagano M,Koinuma S,et al. Effects of aging on basement membrane of tibialis anterior muscle during recovery following muscle injury in rats[J]. Microscopy (Oxf),2022,71(4):245-248. doi: 10.1093/jmicro/dfac016 [15] Fitzgerald J. WARP: A unique extracellular matrix component of cartilage,muscle,and endothelial cell basement membranes[J]. Anat Rec (Hoboken),2020,303(6):1619-1623. [16] Nissen N I,Karsdal M,Willumsen N. Collagens and cancer associated fibroblasts in the reactive stroma and its relation to Cancer biology[J]. J Exp Clin Cancer Res,2019,38(1):115-127. doi: 10.1186/s13046-019-1110-6 [17] Banerjee S,Lo W C,Majumder P,et al. Multiple roles for basement membrane proteins in cancer progression and EMT[J]. Eur J Cell Biol,2022,101(2):151220. doi: 10.1016/j.ejcb.2022.151220 [18] Siegel R L,Miller K D,Fuchs H E,et al. Cancer statistics,2022[J]. CA Cancer J Clin,2022,72(1):7-33. doi: 10.3322/caac.21708 [19] Liu S J,Dang H X,Lim D A,et al. Long noncoding RNAs in cancer metastasis[J]. Nat Rev Cancer,2021,21(7):446-460. [20] Li X,Kuang Q,Peng M,et al. Basement membrane-associated lncRNA risk model predicts Prognosis and Guides Clinical Treatment in Clear Cell Renal Cell Carcinoma[J]. Biomedicines,2023,11(10):2635-2655. doi: 10.3390/biomedicines11102635 [21] 安梦霞,王萍玉. 基于铜死亡相关LncRNAs构建肺鳞癌预后预测模型[J]. 肿瘤防治研究,2023,50(11):1084-1090. [22] Jayadev R,Morais M,Ellingford J M,et al. A basement membrane discovery pipeline uncovers network complexity,regulators,and human disease associations[J]. Sci Adv,2022,8(20):eabn2265. doi: 10.1126/sciadv.abn2265 [23] Siegel L,Miller K D,Wagle N S,et al. Cancer statistics,2023[J]. CA Cancer J Clin,2023,73(1):17-48. doi: 10.3322/caac.21763 [24] Li X,Fliegel L. Permissive role of Na(+)/H(+) exchanger isoform 1 in migration and invasion of triple-negative basal-like breast cancer cells[J]. Mol Cell Biochem,2022,477(4):1207-1216. [25] Deutsch J S,Lai J,Schenk K M,et al. Immune microenvironment of basal cell carcinoma and tumor regression following combined PD-1/LAG-3 blockade[J]. J Immunother Cancer,2023,11(12):e007463. doi: 10.1136/jitc-2023-007463 -