摘要:目的 基于生物信息学方法分析和探讨肉瘤组织中COL6A1基因的表达水平及其意义。方法 首先使用GEO数据库下载GSE16088、GSE21122和GSE49972数据集并分析COL6A1~3基因的表达水平。通过Kaplan-Meier生存曲线分析COL6A1~3基因的表达水平对肉瘤患者预后的影响。采用GSVA包对肉瘤组织当中COL6A1基因的表达水平与免疫细胞浸润的相关性进行Spearman分析。采用LinkedOmics数据库分析与COL6A1共表达的基因,R软件的ggstatsplot包分析COL6A1和COL6A2基因的相关性。使用R软件的ClusterProfiler包对COL6A1共表达基因进行GO功能富集分析和KEGG通路分析。采用GSEA分析肉瘤组织中COL6A1共表达基因的激酶、miRNA及转录因子靶标。结果 与正常组织相比,COL6A1基因在骨肉瘤、软组织肉瘤和肾透明细胞肉瘤等组织中的表达水平均显著升高。生存分析显示COL6A1基因的高表达与较差的总体生存率和疾病特异生存率有关(χ2=6.37、4.57,P<0.05)。免疫细胞浸润分析显示COL6A1基因与中性粒细胞浸润水平之间存在显著正相关(r=0.280,P<0.05)。共表达分析结果显示与COL6A1表达水平显著正相关的基因有2 444个,与COL6A1表达水平显著负相关的基因3 711个(FDR<0.05)。相关性分析显示COL6A1和COL6A2基因之间存在显著正相关(r=0.830,P<0.05)。GO功能富集分析结果显示肉瘤组织中COL6A1共表达基因主要参与免疫细胞浸润等生物学过程。GSEA分析结果显示与COL6A1共表达基因相互作用的5个最重要的激酶靶标是PRKCA、SRC、LYN、PRKAA1和PAK1;重要的miRNA靶标是miR-506、miR-200B、miR-200C、miR-429和miR-330;重要转录因子靶标是V$ETS2_B、V$AP1_C、V$PU1_Q6、V$RORA1_01和V$AP1_Q2。结论COL6A1基因在肉瘤组织中显著高表达,其高表达与肉瘤患者的不良预后密切相关。COL6A1基因有可能作为肉瘤的一个有价值的预后指标和潜在的治疗靶点。
肉瘤是一种间充质来源的异质性恶性肿瘤,比较罕见,约占所有儿童恶性肿瘤的21%,所有成人恶性肿瘤的1%[1]。目前已鉴定出100多种不同的肉瘤实体,大致分为软组织肉瘤和骨肉瘤[2,3]。目前肉瘤的标准治疗方案主要有手术切除、化学疗法和放射疗法,但近几年来患者的疗效和预后并没有显著提高[4]。Ⅵ型胶原A(COL6A)是蛋白质超家族中的一员,由3条主要的多肽链(α1、α2和α3)组成,是主要的细胞外基质之一[5]。研究发现,在宫颈癌、前列腺癌、肺癌以及胰腺癌等肿瘤组织中COL6A1基因呈显著高表达,并且与肿瘤进展密切相关[6,7,8,9],但COL6A1基因在肉瘤组织中的表达水平及其作用机制尚未见相关研究报道。本研究旨在通过生物信息学方法分析COL6A1基因在肉瘤发生发展及预后中的作用及其机制,为肉瘤患者的治疗提供参考依据。现将结果报告如下。
1、材料与方法
1.1 GEO数据库不同类型肉瘤组织中COL6A1~3基因表达水平
在GEO数据库(https: //www.ncbi.nlm.nih.gov/geo/)选取并下载骨肉瘤数据集GSE16088[10]、软组织肉瘤数据集GSE21122[11]和肾透明细胞肉瘤数据集GSE49972[12]的原始测序数据,用R软件Limma软件包对数据标准化后,分析骨肉瘤、软组织肉瘤和肾透明细胞肉瘤组织中COL6A1~3基因表达水平。肉瘤组织和正常组织中差异基因比较用Wilcoxon检验,以P<0.05为差异有统计学意义。
1.2肉瘤组织中COL6A1~3基因的表达水平与患者的生存、免疫细胞浸润的相关性分析
从TCGA数据库(https: //genome-cancer.ucsc.edu/)中下载肉瘤组织以及正常组织样本的转录组RNAseq数据和患者的生存数据。采用R软件的survival包分析肉瘤组织中COL6A1~3基因的表达水平与患者生存的相关性,并通过survminer包绘制患者的Kaplan-Meier(KM)生存曲线[13],并进行Log-rank检验,以P<0.05为差异具有统计学意义。采用GSVA包[14]对肉瘤组织中COL6A1基因的表达水平与24种免疫细胞[15]浸润之间的相关性进行Spearman相关性分析,以P<0.05为差异具有统计学意义。
1.3肉瘤组织中与COL6A1共表达基因的分析以及COL6A1和COL6A2基因表达的相关性分析
通过LinkedOmics数据库(http: //www.linkedomics.org/login.php)[16]分析肉瘤组织当中与COL6A1存在相关性的其他差异表达的基因,即COL6A1的共表达基因,分析结果以皮尔逊相关系数表示。通过R软件的ggstatsplot包对肉瘤组织中COL6A1和COL6A2基因表达的相关性进行Spearman分析,以P<0.05为差异有统计学意义。
1.4肉瘤组织中COL6A1共表达基因的GO功能富集分析和KEGG通路分析
使用R软件的clusterProfiler包对与COL6A1共表达基因进行GO功能富集分析和KEGG通路分析[17]。
1.5肉瘤组织当中COL6A1共表达基因的激酶、miRNA及转录因子靶标分析
在分子特征数据库当中[18]使用GSEA分析COL6A1共表达基因的激酶、miRNA及转录因子靶标,等级标准参数即错误检测率(FDR)<0.05,且进行了500次模拟。
2、结 果
2.1不同类型肉瘤组织中COL6A1~3基因的表达水平
GEO数据库分析结果显示,与正常组织相比,骨肉瘤和软组织肉瘤组织中COL6A1~3基因的表达水平均显著升高(t=2.68~7.98,P<0.05),肾透明细胞肉瘤组织当中COL6A1以及COL6A2基因的表达水平亦显著升高(t=2.95、2.80,P<0.05),COL6A3表达水平无显著差异(P>0.05)。
2.2肉瘤组织中COL6A1~3基因的表达水平与预后的相关性分析
KM生存曲线分析结果显示,肉瘤组织当中COL6A1基因表达水平与患者的总体生存率(OS)和疾病特异生存率(DSS)有关(χ2=6.37、4.57,P<0.05),详见图1A、B;肉瘤组织当中COL6A2以及COL6A3的高表达与患者的OS和DSS无显著相关性(P>0.05)。
2.3肉瘤组织中COL6A1基因的表达水平与免疫细胞浸润的相关性
GSVA包分析结果显示,肉瘤组织中COL6A1基因的表达水平与中性粒细胞、辅助性T细胞2和成熟树突状细胞浸润存在显著正相关(r=0.240~0.280,P<0.05)。
2.4肉瘤组织中与COL6A1共表达基因分析以及COL6A1和COL6A2基因表达的相关性分析
在LinkedOmics数据库中,分析肉瘤组织当中与COL6A1共表达的基因情况,分析结果显示,与COL6A1表达水平显著正相关的基因2 444个,与COL6A1表达水平显著负相关的基因有3 711个(FDR<0.05),图2A、B分别显示了与COL6A1呈现显著正相关和负相关的前50个基因,其中红色表示了高表达基因,蓝色表示低表达基因。R软件的ggstatsplot包进一步分析结果显示,肉瘤组织当中COL6A1和COL6A2的表达呈显著正相关性(r=0.830,P<0.05)。
2.5肉瘤组织中与COL6A1共表达基因的GO功能富集分析和KEGG通路分析
GO功能富集分析结果显示,与COL6A1相关的共表达基因主要位于突触前、神经元细胞体和细胞-基质连接等成分中,主要参与多细胞生物体内平衡、中性粒细胞介导的免疫反应、中性粒细胞参与的活化免疫反应等生物学过程。KEGG通路分析结果显示,这些基因主要在PI3K-Akt信号通路、人乳头瘤病毒感染、神经元活性配体-受体相互作用和MAPK信号通路中富集。
2.6肉瘤组织当中COL6A1共表达基因的激酶、miRNA及转录因子靶标分析
肉瘤组织中与COL6A1共表达基因的GSEA分析结果显示,与COL6A1共表达基因相互作用的5个最重要的激酶靶标是PRKCA、SRC、LYN、PRKAA1和PAK1;重要的miRNA靶标是miR-506、miR-200B、miR-200C、miR-429和miR-330;重要转录因子靶标是V$ETS2_B、V$AP1_C、V$PU1_Q6、V$RORA1_01和V$AP1_Q2。
图1肉瘤组织中COL6A1基因的表达水平与患者预后的相关性
图2与COL6A1共表达基因热图分析
3、讨 论
肉瘤是相对罕见的异质性恶性肿瘤,患者死亡率极高[19],目前的治疗方法主要还是手术切除、放疗和化疗,患者的预后很差[20]。因此迫切需要深入探讨肉瘤发病的分子机制,挖掘其有效的预后和治疗靶点。
COL6A1是脊椎动物中的保守基因,几乎存在于机体内的所有结缔组织当中,以维持组织的完整性[21]。研究发现,肿瘤组织中COL6A1的高表达与肿瘤进展密切相关[22]。COL6A1还被发现参与多种肿瘤细胞的黏附、侵袭和转移过程[9,23]。此外,前列腺癌细胞系中COL6A1基因敲减以后会通过JK-STAT途径抑制前列腺癌细胞的增殖[7]。目前关于COL6A1在肉瘤组织中的表达情况及其作用机制尚不清楚。
本研究首先通过在GEO数据库的GSE16088、GSE21122和GSE49972数据集中对肉瘤组织和正常组织中的COL6A1~3基因表达进行差异分析,结果显示,与正常组织相比,COL6A1基因在骨肉瘤、软组织肉瘤和肾透明细胞肉瘤等组织中的表达水平均显著升高。进一步KM生存曲线分析显示,肉瘤组织中COL6A1基因高表达的患者,预后更差。以上结果表明,COL6A1可能在肉瘤的发生发展中具有重要作用。
为进一步了解COL6A1共表达基因在肉瘤中的作用机制,本研究对肉瘤组织中COL6A1共表达基因进行了GO功能富集分析和KEGG通路分析。GO功能富集分析结果显示,肉瘤组织中COL6A1共表达基因主要参与多细胞生物体内平衡、中性粒细胞介导的免疫反应和中性粒细胞参与的活化免疫反应等生物学过程。免疫细胞浸润相关性分析同样表明,肉瘤组织中COL6A1基因的表达水平与中性粒细胞的细胞浸润水平之间存在显著相关性。这表明COL6A1表达可能主要通过影响肿瘤微环境中的免疫浸润,尤其是中性粒细胞的细胞浸润,进而影响患者的预后和生存。在过去的几年中,越来越多的证据表明免疫系统在肿瘤的控制和进展中起着重要作用[24]。免疫疗法在某些类型的肿瘤,例如非小细胞肺癌以及黑素瘤的治疗当中获得了喜人的结果[25],本研究结果可能对肉瘤的免疫治疗具有一定的意义。此外,本研究KEGG通路分析显示,肉瘤组织中COL6A1共表达基因主要富集于PI3K-Akt信号通路、人乳头瘤病毒感染、神经元活性配体-受体相互作用以及MAPK信号通路中,因此这些基因表达水平的变化可能与肉瘤的发生发展具有重要的关系。
本研究通过GSEA分析肉瘤组织中COL6A1共表达基因的重要激酶、miRNA及转录因子靶标,以进一步揭示COL6A1在肉瘤中的作用机制。本研究结果显示,与COL6A1共表达基因相互作用的5个最重要的激酶靶标是PRKCA、SRC、LYN、PRKAA1、和PAK1;重要的miRNA靶标是miR-506、miR-200B、miR-200C、miR-429和miR-330;重要转录因子靶标是V$ETS2_B、V$AP1_C、V$PU1_Q6、V$RORA1_01和V$AP1_Q2。然而这些激酶、miRNA及转录因子调控机制比较复杂,需要进一步实验验证,以期为临床治疗提供重要指导。虽然LinkedOmics数据库的共表达分析以及R软件ggstatsplot包的相关性分析均表明COL6A1和COL6A2之间呈显著正相关,COL6A家族成员之间的相互作用及其具体的分子机制也仍有待进一步研究。
综上所述,本研究通过对多个数据库当中肉瘤患者的RNA测序数据进行生物信息学分析提示,COL6A1基因在肉瘤组织中显著高表达,COL6A1基因的高表达也与肉瘤患者的不良预后密切相关,COL6A1可能是肉瘤的一个有价值的预后指标和潜在的治疗靶点。本研究可为肉瘤患者的预后及治疗提供新的思路和理论依据。
作者声明:胡伟全、赵琴飞、陈芳芳参与了研究设计;胡伟全、赵琴飞、朱龙玉、胡素萍、陈芳芳参与了论文的写作和修改。所有作者均阅读并同意发表该论文,且均声明不存在利益冲突。
基金资助:国家自然科学基金项目(82002233);江西省中医药科技计划项目(2022A136);江西省卫生健康委科技计划项目(202311903);赣州市指导性科技计划项目(20222ZDX7783);
文章来源:胡伟全,赵琴飞,朱龙玉等.肉瘤组织中COL6A1基因的表达水平及其意义[J].精准医学杂志,2023,38(06):525-529.
分享:
目前运动疗法也已经作为一种非药物干预手段被应用于肿瘤的预防与治疗[5],当前研究认为运动疗法对肿瘤微环境以及肿瘤细胞存在多项干预机制,可以通过调节机体-肿瘤的相互作用来调节肿瘤的生长进程,而这种调节作用是相互、复杂、多面、重叠的[6]。其中运动针对肿瘤血管的抗癌策略与传统策略有所不同,但最终的目标导向同样都是抑制肿瘤的发生发展。
2024-04-23临床预测模型(clinical prediction models, CPMs),又称临床预测规则(clinical prediction rules)、风险预测模型(risk prediction models)、预测模型(predictive models)或风险评分(risk scores),是指利用数学公式估测某研究对象患有某病或将来发生某结局的概率,其包括参数、半参数或非参数模型。现临床广泛使用的主要为logistic回归或Cox回归,也有一些机器学习的算法应用在临床预测模型的构建上。
2024-04-18富含半胱氨酸的酸性分泌蛋白(secreted protein acidic and rich in cysteine, SPARC)是一种钙结合糖蛋白,它属于细胞外基质(extracellular matrix, ECM)的非结构成分家族,广泛分布于体内[1]。SPARC 的表达集中在骨和内脏上皮细胞等发育和重建活跃的器官和组织中,其在肿瘤组织中的表达和分泌被认为是一种重要的影响肿瘤进展的因素[2]。
2024-04-10喉鳞状细胞癌(laryngeal squamous cell carcinoma, LSCC)是头颈部鳞状细胞癌(head and neck squamous cell carcinoma, HNSCC)最常见的一种。2020年,喉癌新增病例184 615例,死亡病例99 840例,其中95%以上为LSCC[1]。早期LSCC患者通过手术或放疗具有良好的治疗效果[2]。但是,由于复发和治疗耐药,晚期 LSCC 患者的 5 年总生存率仍然较低。已有报道肿瘤干细胞(cancer stem cells, CSC
2024-04-09FGF19作为成纤维细胞因子亚家族成员之一,参与了消化代谢与肠道肿瘤的发生,对葡萄糖和脂质代谢都具有调节作用,并且对细胞内环境稳态维持具有潜在的影响[1]。临床研究一直在探讨把FGF19作为肿瘤的潜在评估靶点,FGF19的基因组位点位于染色体11q13.3,在人类癌症中经常扩增,与肿瘤的发生和进展密切相关[2]。然而FGF19参与消化系统肿瘤的发生和进展的机制尚未见报道,为此我们进行了FGF19基因克隆实验[3]。
2024-04-08霍奇金淋巴瘤(HL)是淋巴系统恶性增殖性肿瘤,发病率占淋巴瘤的5%-10%[1]。HL病理表现为来源于B淋巴细胞的多核细胞(即RS细胞),RS细胞的产生与周围的免疫微环境密切相关。经典型(c HL)包括4种亚型,分别为结节硬化型、混合细胞型、富于淋巴细胞型和淋巴细胞减少型。约90%的HL患者为c HL,且RS细胞表达细胞表面CD30抗原[2]。另有10%的HL患者是结节淋巴细胞为主型,其特征细胞为“爆米花细胞”,是表达B淋巴细胞抗原CD20的大细胞[2]。
2024-04-02缺氧诱导因子(HIFs)家族是参与调节细胞代谢的细胞因子,在肿瘤的发生发展过程中起重要作用。HIFs的主要功能包括促进肿瘤发生、血管生成、生长、转移、调节肿瘤细胞增殖、抑制细胞凋亡、调控能量代谢密切相关[1]。HIFs高表达于肿瘤组织,分布于肿瘤微环境,与多种肿瘤细胞的上皮-间质转化(EMT)过程相关[2]。
2024-04-01胃肠胰神经内分泌肿瘤(neuroendocrine tumors, NETs)根据发病位置不同划分NETs不同类型。NETs发病隐匿,往往恶性程度高。其中,胰腺NETs是胰腺癌中恶性程度非常高的肿瘤。而结肠直肠NETs是胃肠道中仅次于胰腺NETs发病率,第二常见的NETs, 在美国每年的发病率约为0.3~1人/10万人[1],在中国的发病率与之相近,在过去40年中,结肠直肠NETs的发病率有所增加[2]。
2024-03-31脑胶质瘤(Glioma)是中枢神经系统最常见的原发性肿瘤,具有高度侵袭性和较高的死亡率,容易复发且治疗难度高[1,2]。环状RNA(circRNA)是一种具有共价闭合环状结构的新型非编码RNA,circRNA与肿瘤、动脉粥样硬化、糖尿病、神经系统疾病等多种疾病相关[3,4]。CircRNAs也在神经系统疾病(如胶质瘤)中发挥重要作用,如circ_0079593可通过海绵吸附miR-324-5p而促进神经胶质瘤的进展[5,6]。
2024-03-26游离DNA是细胞凋亡、坏死和自噬等过程中释放到血液、腹水、脑脊液、唾液、胸腔积液和尿液等体液中的DNA片段,长度为150~200 bp。游离DNA的产生和清除是一个动态过程,半衰期为5~150 min, 因此游离DNA可以对全身死亡细胞进行“快照”,这一特点使其成为许多疾病的检测目标。
2024-03-16人气:14053
人气:13234
人气:13041
人气:12811
人气:11638
我要评论
期刊名称:精准医学杂志
期刊人气:1763
主管单位:山东省教育厅
主办单位:青岛大学
出版地方:山东
专业分类:医学
国际刊号:2096-529X
国内刊号:37-1515/R
邮发代号:24-130
创刊时间:1986年
发行周期:双月刊
期刊开本:大16开
见刊时间:1-3个月
影响因子:0.000
影响因子:0.000
影响因子:0.000
影响因子:0.000
影响因子:0.000
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!