HMGA2在胰腺癌组织中表达对患者预后的影响及促进肿瘤进展的机制分析

Influence of HMGA2 expression in pancreatic cancer tissues on prognosis of patients and mechanism analysis of tumor progression promotion

  • 摘要:
    探讨高迁移率族蛋白A2(HMGA2)在胰腺癌组织中表达对患者预后的影响及促进肿瘤进展的机制。
    采用回顾性队列研究和实验研究方法。收集2018年1月至2025年2月复旦大学附属中山医院收治的161胰腺癌患者的临床病理资料;男102例,女59例;年龄为64(56~70)岁。161例患者中,根据不同时期收集患者资料分为3个队列,队列1为2018年9月至2019年10月收治的29例患者,采集肿瘤及癌旁样本分别制作组织芯片;队列2为2020年11月至2021年12月收治的116例患者,采集肿瘤样本制作组织芯片;队列3为2025年1—2月收治的16例患者,采集肿瘤样本制作普通石蜡切片。采用免疫组织化学染色方法评估HMGA2在胰腺癌中的表达情况,并分析其与患者预后的关系。通过生物信息学分析HMGA2在胰腺癌中的生物学功能,并通过实时荧光定量聚合酶链反应(qRT⁃PCR)、蛋白印迹实验、CCK⁃8实验、流式细胞实验进行验证。联合免疫组织化学染色和流式细胞实验分析HMGA2表达水平与胰腺癌微环境特征的相关性,并采用胰腺癌单细胞测序数据分析HMGA2与肿瘤相关中性粒细胞(TANs)的关系。观察指标:(1)HMGA2在胰腺癌中的表达及患者预后情况。(2)HMGA2在胰腺癌中的功能分析。(3)胰腺癌HMGA2表达与免疫微环境的相关性分析。(4)TANs参与HMGA2介导的免疫逃逸分析。正态分布的计量资料两组间比较采用独立样本t检验,多组间比较采用单因素方差分析。偏态分布的计量资料两组间比较采用Mann⁃Whitney U检验。相关性分析采用Pearson相关性分析。采用Kaplan⁃Meier法计算生存率并绘制生存曲线,Log⁃rank检验进行生存分析。
    (1)HMGA2在胰腺癌中的表达及患者预后情况。队列1、2患者免疫组织化学染色结果显示:HMGA2主要表达于肿瘤细胞,其在胰腺癌和癌旁组织中的表达水平分别为1(0~1)和2(1~2),两者比较,差异有统计学意义(Z=-3.77,P<0.05)。队列1中,HMGA2高表达和低表达胰腺癌患者中位总生存时间分别为17个月和36个月,两者比较,差异有统计学意义(风险比=3.11,95%可信区间为1.05~9.23,P<0.05);队列2中,HMGA2高表达和低表达胰腺癌患者中位总生存时间分别为22个月和未达到,两者比较,差异有统计学意义(风险比=3.47,95%可信区间为2.00~6.02,P<0.05)。(2)HMGA2在胰腺癌中的功能分析。癌症基因组图谱数据库分析结果显示:HMGA2高表达患者在GO条目和KEGG通路中基因富集于染色质分离、DNA复制、细胞周期等细胞增殖相关过程;HMGA2高表达患者KRAS突变频率高于低表达患者;经典型和基底型胰腺癌HMGA2表达水平分别为0.95(0.44~1.52)和1.91(0.96~2.87),两者比较,差异有统计学意义(Z=-2.41,P<0.05);HMGA2与基底型胰腺癌KRT5表达呈正相关(r=0.38,P<0.05),与经典型胰腺癌GATA6表达呈负相关(r=-0.18,P<0.05)。在胰腺癌细胞系PANC⁃1中进行功能学实验,采用小干扰RNA(siRNA)抑制HMGA2的表达。qRT⁃PCR分析结果显示:对照组、siRNA#1组、siRNA2#组、siRNA3#组HMGA2表达水平分别为0.99±0.03、0.43±0.05、0.76±0.06、0.31±0.04,4组比较,差异有统计学意义(F=181.60,P<0.05);蛋白印迹实验结果显示:与对照组比较,siRNA#1组和siRNA#3组的HMGA2条带灰度均更低;CCK⁃8实验结果显示:接种96 h后,对照组、siRNA#1组、siRNA#3组细胞吸光度值分别为2.23±0.12、1.83±0.07、1.73±0.06,3组比较,差异有统计学意义(F=46.09,P<0.05)。流式细胞实验结果显示:以对照组为参照,siRNA#1组早、晚期凋亡细胞率共计19.99%±1.99%,siRNA#3组早、晚期凋亡细胞率共计22.94%±1.51%。(3)胰腺癌HMGA2表达与免疫微环境的相关性分析。流式细胞实验结果显示:队列1中,HMGA2高表达和低表达患者CD8+ T细胞浸润水平、TANs浸润水平分别为(76±58)、(118±92)个和(119±42)、(58±49)个,两者上述指标比较,差异均有统计学意义(t=-2.30,2.13,P<0.05);队列3中,HMGA2高表达和低表达患者CD8+ T细胞表面PD⁃1、TIM⁃3表达分别为(1 398±1 214)、(3 930± 2 267)个和(263±254)、(2 030±979)个,两者上述指标比较,差异均有统计学意义(t=2.59,2.18,P<0.05)。(4)TANs参与HMGA2介导的免疫逃逸分析。癌症基因组图谱、国际癌症基因组联合体、基因表达综合数据库、欧洲生物信息学研究所数据库胰腺癌细胞测序结果显示:胰腺癌HMGA2表达与TANs浸润呈正相关,与CD8+ T细胞浸润呈负相关。GSE155698单细胞测序样本结果显示:与HMGA2低表达比较,HMGA2高表达胰腺癌患者TANs浸润水平更高,亚群TANs_ cluster 0、TANs_cluster 1、TANs_cluster 2、TANs_cluster 3浸润水平也均增加,其中TANs_cluster 2高表达血管内皮生长因子。
    HMGA2在胰腺癌组织中高表达,高表达患者预后更差;抑制HMGA2表达可抑制胰腺癌细胞增殖,促进凋亡;与HMGA2低表达比较,HMGA2高表达胰腺癌中CD8+ T细胞浸润水平更低,TANs及其亚群的浸润水平更高。

     

    Abstract:
    Objective To investigate the influence of high mobility group protein AT-Hook-2 (HMGA2) expression in pancreatic cancer tissues on prognosis of patients and mechanism of tumor progression promotion.
    Methods The retrospective cohort and experiment study was conducted. The clinicopathological data of 161 patients with pancreatic cancer who were admitted to Zhongshan Hospital of Fudan University from January 2018 to February 2025 were collected. There were 102 males and 59 females, aged 64 (range, 56-70)years. Among 161 patients, data were collected and divided into 3 cohorts based on different time periods. Cohort 1 included 29 patients admitted from September 2018 to October 2019, of whom tumor and paracancerous tissue samples were collected to prepare tissue microarrays (TMAs). Cohort 2 included 116 patients admitted from November 2020 to December 2021, and their tumor samples were collected for TMAs preparation. Cohort 3 included 16 patients admitted from January to February 2025, and their tumor samples were collected to make conventional paraffin sections. The immunohistochemical staining method was used to evaluate the expression of HMGA2 in pancreatic cancer and analyze its relationship with the prognosis of patients. The biological function of HMGA2 in pancreatic cancer was analyzed by bioinformatics and verified by quantitative real-time polymerase chain reaction (qRT-PCR), Western blot, cell counting kit-8 (CCK8) assay, flow cytometry. The immunohistochemical staining and flow cytometry were used to identify the correlation between HMGA2 expression and microenvironment of pancreatic cancer. The relationship between HMGA2 and tumor-associated neutrophils (TANs) was analyzed using single cell sequencing data of pancreatic cancer. Observation indicators: (1) expression of HMGA2 in pancreatic cancer and prognosis of patients; (2) functional analysis of HMGA2 in pancreatic cancer; (3) correlation between HMGA2 expression and immune microenvironment in pancreatic cancer; (4) analysis of TANs participation in HMGA2 mediated immune escape. Comparison of measurement data with normal distribution between groups was conducted using the independent sample t test, and one-way analysis of variance was used for multiple group comparisons. Comparison of measurement data with skewed distribution between groups was conducted using the Mann-Whitney U test. The Pearson correlation analysis was used for inter group correlation analysis. The Kaplan-Meier method was used to calculate survival rate and plot survival curves, and Log-rank test was used for survival analysis.
    Results (1) Expression of HMGA2 in pancreatic cancer and prognosis of patients. For patients of cohort 1 and 2, immunohistochemical staining analysis showed that HMGA2 mainly localizing in tumor cells, with expression levels of 1 (range, 0-1) in cancer tissues and 2 (range, 1-2) in adjacent tissues, respectively, showing a significant difference between them (Z=-3.77, P<0.05). In cohort 1, the median survival time of pancreatic cancer patients with high and low HMGA2 expression was 17 months and 36 months, respectively, showing a significant difference between them (hazard ratio=3.11, 95% confidence interval as 1.05-9.23, P<0.05). In cohort 2, the median survival time of pancreatic cancer patients with high and low HMGA2 expression was 22 months and not reached, respectively, showing a significant difference between them (hazard ratio=3.47, 95% confidence interval as 1.997-6.016, P<0.05). (2) Functional analysis of HMGA2 in pancreatic cancer. The results of gene set enrichment analysis of patients with high and low HMGA2 expression in pancreatic cancer patients in The Cancer Genome Atlas database showed that the differential genes in GO entry and KEGG pathway were enriched in cell proliferation related processes such as chromatin separation, DNA replication, cell cycle, etc. Patients with high HMGA2 expression had a high KRAS mutation frequency than patients with low HMGA2 expression. In classic-type and basal-type pancreatic cancer, the expression levels of HMGA2 were 0.95 (range, 0.44-1.52) and 1.91 (range, 0.96-2.87), respectively, showing a significant difference between the two groups (Z=-2.41, P<0.05). HMGA2 was positively correlated with the expression of KRT5 in basal-type pancreatic cancer (r=0.38, P<0.05), while it was negatively correlated with the expression of GATA6 in classic-type pancreatic cancer (r=-0.18, P<0.05). Functional experiments were performed in the pancreatic cancer cell line PANC-1, in which the small interfering RNA (siRNA) was used to inhibit HMGA2 expression. Results of qRT-PCR analysis showed that the HMGA2 expression levels in the control group, siRNA#1 group, siRNA#2 group, and siRNA#3 group were 0.99±0.03, 0.43±0.05, 0.76±0.06, and 0.31±0.04, respectively, showing a significant difference among the four groups (F=181.60, P<0.05). Western Blotting revealed that the grayscale values of HMGA2 bands in the siRNA#1 group and siRNA#3 group were lower than those in the control group. Results of the CCK-8 assay showed that at 96 hours after cell seeding, the absorbance values in the control group, siRNA#1 group, and siRNA#3 group were 2.23±0.12, 1.83±0.07, and 1.73±0.06, respectively, showing a significant difference among the three groups (F=46.09, P<0.05). Flow cytometry indicated that taking the control group as a reference, the early and late apoptotic cells rates were 19.99%±1.99% in the siRNA#1 group and 22.94%±1.51% in the siRNA#3 group. (3) Correlation between HMGA2 expression and immune microenvironment in pancreatic cancer. Results of flow cytometry showed that in cohort 1, the infiltration levels of CD8+ T cells and tumor-associated neutrophils (TANs) were 76±58 and 118±92 in patients with high HMGA2 expression, versus 119±42 and 58±49 in patients with low HMGA2 expression, respectively,showing significant differences between the two groups (t=-2.30, 2.13, P<0.05). In cohort 3, the expression levels of PD-1 and TIM-3 on the surface of CD8+ T cells were 1 398±1 214 and 3 930±2 267 in patients with high HMGA2 expression, versus 263±254 and 2 030±979 in patients with low HMGA2 expression, respectively, showing significant differences between the two groups (t=2.59, 2.18, P<0.05). (4) Analysis of TANs participation in HMGA2 mediated immune escape. Sequencing results of pancreatic cancer cells from The Cancer Genome Atlas, International Cancer Genome Consortium, Gene Expression Omnibus, and European Bioinformatics Institute database showed that HMGA2 expression in pancreatic cancer was positively correlated with TANs infiltration and negatively correlated with CD8+ T cell infiltration. Results from single-cell sequencing samples of GSE155698 demonstrated that compared with patients with low HMGA2 expression, those with high HMGA2 expression had a higher TANs infiltration level, the infiltration levels of TANs subsets (TANs_cluster 0, TANs_cluster 1, TANs_cluster 2, and TANs_cluster 3) also increased. Among these subsets, TANs_cluster 2 highly expressed vascular endothelial growth factor.
    Conclusions HMGA2 is highly expressed in pancreatic cancer, and the prognosis of patients with high expression is worse. Inhibition of HMGA2 expression can suppress cell proliferation and promote apoptosis in pancreatic cancer. Compared with low HMGA2 expression, pancreatic cancer with high HMGA2 expression exhibits lower infiltration level of CD8+ T cells, as well as higher infiltration levels of TANs and their subsets.

     

/

返回文章
返回