甲状腺癌中m6A甲基化调控因子的表达及其预后价值
作者:
通讯作者:
作者单位:

1.中南大学湘雅医院,甲状腺外科,湖南长沙 410013;2.中南大学湘雅医院,国际医疗部,湖南长沙 410013;3.中南大学湘雅医院,中南大学湘雅医学院,湖南长沙 410013

作者简介:

王文龙,中南大学湘雅医院博士研究生,主要从事甲状腺癌表观遗传调控和分子诊断方面的研究。

基金项目:

国家自然科学基金资助项目(81672885);湖南省自然科学基金资助项目(2019JJ40500);中南大学研究生创新型课题基金资助项目(2020zzts259)。


Expressions of m6A methylation regulators and their prognostic value in thyroid cancer
Author:
Affiliation:

1.Department of Thyroid Surgery, Central South University, Changsha 410013, China;2.International Medical Department, Xiangya Hospital, Central South University, Changsha 410008, China, Central South University, Changsha 410013, China;3.Xiangya School of Medicine, Central South University, Changsha 410013, China

Fund Project:

  • 摘要
  • |
  • 图/表
  • |
  • 访问统计
  • |
  • 参考文献
  • |
  • 相似文献
  • |
  • 引证文献
  • |
  • 资源附件
  • |
  • 音频文件
  • |
  • 视频文件
    摘要:

    背景与目的 甲状腺癌发病率逐年升高,但其发病机制仍不清楚,阐述甲状腺癌的发病机制对改善甲状腺癌患者的预后至关重要。研究表明,m6A甲基化调控因子高度参与癌症发生发展,具有良好的潜在预后价值。因此,本研究通过生物信息学方法分析甲状腺癌中m6A甲基化调控因子的表达并构建基于m6A甲基化调控因子的甲状腺癌预后模型。方法 从TCGA数据库下载甲状腺癌m6A甲基化调控因子的表达数据和相应的临床病理资料,通过Wilcoxon检验分析20个m6A甲基化调控因子在肿瘤和正常组织的差异表达;用一致性聚类分析将甲状腺癌患者分为两个聚类,比较两个聚类患者临床病理因素和总体生存率的差异;Lasso Cox回归分析构建风险预测模型并用ROC曲线下面积(AUC)评估模型的预测能力。结果 19个m6A甲基化调控因子在甲状腺癌和正常组织表达具有统计学差异(均P<0.05),其中HNRNPC、IGF2BP2、FMR1在甲状腺癌组织中明显高表达,而其余表达下调。聚类分析示,cluster 1生存期低于cluster 2(P<0.05),颈淋巴结转移发生率明显高于cluster 2(P<0.01)。基于Lasso Cox回归分析筛选的4个基因(IGF2BP2、RBM15、YTHDF1、YTHDF3)构建风险评估模型,相比低风险组患者,高风险患者生存期明显缩短(P=0.007);ROC曲线示,该模型可以预测甲状腺癌患者的预后(AUC=0.731)。结论 m6A甲基化调控因子在甲状腺癌中存在差异表达,基于甲状腺癌中关键m6A甲基化调控因子所构建的预后风险模型具有较好的预测能力,可为临床决策提供一定的依据。

    Abstract:

    Background and Aims The incidence of thyroid carcinoma has been steadily increasing and the molecular mechanism underlying the tumorigenesis is still unclear. Thus, exploring the underlying mechanism is of great importance for improving the prognosis of PTC patients. Studies demonstrated that m6A methylation regulators are widely involved in the occurrence and development of cancers and have superior prognostic value. Therefore, this study was conducted to investigate the expressions of m6A methylation regulators in thyroid cancer and construct a prognostic risk model for thyroid cancer based on m6A methylation regulators using bioinformatics approaches.Methods The gene expression profiles of m6A RNA methylation regulators and the corresponding clinical information were downloaded from The Cancer Genome Atlas (TCGA). The differential expressions of 20 m6A methylation regulators between tumor and normal tissues were analyzed by Wilcoxon test. The thyroid cancer patients were divided into two clusters by consensus clustering, and the differences in clinicopathologic factors and overall survival rate between the two clusters of patients were compared. Subsequently, the risk model was constructed by Lasso Cox regression analysis and its predictive value was evaluated by the area under the ROC (AUC).Results There were 19 m6A methylation regulators that showed significantly different expressions between thyroid cancer and normal tissues (all P<0.05), in which, the HNRNPC, IGF2BP2, FMR1 were remarkably upregulated, while the remaining were down-regulated in thyroid cancer tissue. The clustering analysis showed that the overall survival of cluster 1 was poorer than that of cluster 2 (P<0.05), while the incidence of cervical lymph node metastasis was significantly higher than that of cluster 2 (P<0.01). A prognostic risk model was constructed based on 4 genes (IGF2BP2, RBM15, YTHDF1 and YTHDF3) that were screened by Lasso Cox regression analysis, in which, the high-risk patients had a worse prognosis than that of low-risk patients (P=0.007). The ROC analysis indicated a reliable prediction performance of the model for patients with thyroid cancer (AUC=0.731).Conclusion There are differential expressions in m6A methylation regulators in thyroid cancer, and the constructed prognostic risk model based on the hub m6A methylation regulators in thyroid cancer has better predictive ability, which will provide certain recommendations for clinical decision-making.

    表 1 492例甲状腺癌患者基本临床病理资料Table 1 The clinicopathologic information of the 492 patients with thyroid cancer
    图1 m6A甲基化调控因子的基因表达及相关性 A:20个m6A甲基化调控因子在甲状腺癌和正常组织表达水平;B:Spearman相关分析m6A甲基化调控因子相关性;C:PPI网络评估m6A甲基化调控因子相互作用关系Fig.1 Gene expression and correlations of the m6A methylation regulators A: The mRNA levels of the 20 m6A regulators in thyroid cancer; B: Spearman correlation analysis of co-expressions among the m6A regulators; C: PPI network of the m6A regulators
    图2 一致性聚类分析m6A甲基化调控因子 A:k值为2时,甲状腺癌样本之间不存在交叉,聚类稳定性好;B:主成分分析区分cluster 1和cluster 2;C:不同的聚类与临床病理的关系;D:总体生存率在两组聚类的情况Fig.2 Consensus clustering analysis of the m6A RNA methylation regulators A: The most appropriate selection with clustering stability when k=2 used; B: Principal component analysis to distinguish cluster 1 and cluster 2; C: Heatmap and clinicopathologic features of the two clusters; D: The overall survival of cluster 1 and cluster 2
    图3 m6A甲基化调控因子构建风险预后模型及检验 A:Lasso Cox回归模型的m6A基因及相关系数;B:风险评分和生存状态分布;C:高、低危组的生存曲线;D:ROC曲线评估模型的预测效能Fig.3 Construction od prognostic risk model and validation A: Coefficients calculated by Lasso Cox regression analysis; B: Risk score and survival status distribution; C: Kaplan-Meier analysis of overall survival in high-risk and low-risk groups; D: Assessment of the prediction accuracy by ROC curve
    参考文献
    相似文献
    引证文献
引用本文

王文龙,沈聪,孙博韬,白宁,李新营,陈嘉欣,彭婀敏.甲状腺癌中m6A甲基化调控因子的表达及其预后价值[J].中国普通外科杂志,2021,30(8):934-941.
DOI:10.7659/j. issn.1005-6947.2021.08.008

复制
分享
文章指标
  • 点击次数:
  • 下载次数:
历史
  • 收稿日期:2021-06-05
  • 最后修改日期:2021-07-16
  • 录用日期:
  • 在线发布日期: 2021-09-02