岳蔷薇 徐乐 张东升

    山东第一医科大学附属省立医院口腔科，济南 250021

    头颈部鳞状细胞癌(head and neck squamous cell carcinoma，HNSCC)是全身最常见的恶性肿瘤之一，死亡率为40%～50%[1]。尽管肿瘤研究飞速进展，但HNSCC 总生存率仍未显著提高[2]，因此亟需新的预后评估分子标志物。N6-甲基腺苷(N6-methyladenosine，m6A) 是真核生物 RNA 中最常见的修饰方式[3-4]，异常表达的m6A 调节基因可能与HNSCC 发生发展密切相关[5]。肿瘤微环境是肿瘤生长的温床，其中浸润性免疫细胞如同双刃剑影响癌细胞[6]。研究[7-8]证实m6A RNA 甲基化与免疫功能密切相关。然而，在HNSCC 中关于肿瘤m6A RNA 甲基化、免疫功能和肿瘤临床特征之间的相关性仍不明确。因此，本实验基于癌症基因组图谱(The Cancer Genome Atlas，TCGA)，通过生物信息学分析探索了HNSCC 中m6A 调节基因的表达、功能、临床意义及其与免疫功能之间的关系，为HNSCC预后和治疗提供指导。

    1 材料和方法

    1.1 数据获取和处理

    在TCGA 数据库(https://tcga-data.nci.nih.gov/tcga/)中下载TCGA-HNSC 队列的转录组数据和临床数据，其中包括502个肿瘤样本和44个正常对照组织样本。临床数据主要包括肿瘤患者相应临床病理特征，包括年龄、性别、病理分级、肿瘤分期和生存数据。

    1.2 m6A调节因子差异表达及相互作用分析

    使用limma R包筛选肿瘤和正常样本之间的差异表达基因，筛选阈值设定为调整P值(adjustedPvalue，adj.P) <0.05，且差异倍数|log2FC|≥1。使用R 语言的ggplot2 软件包绘制差异基因提琴图及Pearson相关分析图。

    1.3 m6A调节因子风险预后模型构建及评估

    利用单因素COX 分析评估m6A 调节因子预后相关性，然后将P值小于0.05的预后相关基因进行LASSO 回归分析构建风险预后模型 ......