《PLOS One》：Identification of prognostic biomarkers for nasopharyngeal carcinoma through microRNA?mRNA network analysis combined with screening of immune?related genes
编辑推荐：
背景：鼻咽癌（nasopharyngeal carcinoma，NPC）是一种常见的头颈部恶性肿瘤，放疗为其主要治疗手段。复发与转移是导致治疗失败的主要原因。尽管免疫治疗可延长患者生存期，但NPC发生发展的分子机制仍不明确。微小RNA（microRNA，miR
背景：鼻咽癌（nasopharyngeal carcinoma，NPC）是一种常见的头颈部恶性肿瘤，放疗为其主要治疗手段。复发与转移是导致治疗失败的主要原因。尽管免疫治疗可延长患者生存期，但NPC发生发展的分子机制仍不明确。微小RNA（microRNA，miRNA）及其靶基因在多种肿瘤中异常表达，并与肿瘤发生密切相关。因此，本研究旨在鉴定NPC中的关键miRNA及免疫相关靶基因，构建miRNA-mRNA调控网络，筛选并验证核心免疫基因，以识别NPC的预后分子生物标志物。
方法：研究人员从基因表达综合数据库（Gene Expression Omnibus，GEO）下载了NPC相关的miRNA数据集（GSE32960）和mRNA数据集（GSE12452、GSE13597、GSE64634），利用R软件鉴定差异表达miRNA（differentially expressed miRNAs，DEMs）和差异表达mRNA（differentially expressed genes，DEGs）。通过TargetScan、miRDB和miRWalk预测DEMs的靶基因，并与DEGs取交集以确定核心基因。在Cytoscape中构建miRNA-mRNA调控网络。将核心基因进一步与免疫相关基因集取交集，随后进行功能富集分析。采用LASSO-Cox回归建立预后风险模型，并通过CIBERSORT评估免疫细胞浸润情况。此外，还进行了模型验证、免疫组织化学分析和药物-基因相互作用分析。体外实验中，在C666-1和NP69细胞系中进行实时定量聚合酶链反应（RT-qPCR）以验证关键基因的表达。
结果：研究人员鉴定了213个关键miRNAs的潜在靶基因。三个miRNA（hsa-miR-30e-5p、hsa-miR-30a-5p和hsa-miR-20a-5p）与NPC预后相关，17个核心免疫相关基因富集于细胞因子介导的信号传导、趋化因子应答和JAK-STAT通路。五个基因（CXCL11、SOCS1、ROBO1、LGR5和RORB）对NPC显示出良好的诊断价值；其表达与CD8+T细胞、巨噬细胞和中性粒细胞浸润呈正相关（P<0.05）。药物-基因分析预测5-氟尿嘧啶靶向LGR5，阿地白介素靶向SOCS1。RT-qPCR结果显示，SOCS1在NPC细胞中上调约3.7倍（P<0.01）。
结论：本研究鉴定了鼻咽癌中潜在的新型关键基因，为该疾病的治疗提供了新的治疗靶点。
一、研究背景与科学问题
鼻咽癌（nasopharyngeal carcinoma，NPC）是头颈部最常见的恶性肿瘤之一，尤其在东亚和东南亚地区发病率显著高于世界其他区域。放疗是NPC的主要治疗手段，但即使放疗技术不断进步，早期至晚期患者单纯放疗后的5年生存率仅为30%–50%。复发和转移仍是治疗失败和预后不良的主要原因。近年来，以程序性死亡受体-1（programmed death receptor-1，PD-1）/程序性死亡配体-1（programmed death ligand-1，PD-L1）免疫检查点抑制剂为代表的免疫治疗虽在一定程度上延长了患者生存期，但仍有一部分患者对该治疗无应答或产生耐药性。肿瘤的恶性表型往往由潜在的分子异常驱动，NPC的发生、进展和远处转移与遗传改变及相关信号通路的失调密切相关。微小RNA（microRNA，miRNA）是一类长度为22–25个核苷酸的内源性非编码RNA，通过与靶mRNA直接结合在转录后水平调控基因表达。大量研究表明，多种miRNA在NPC中异常表达并调节多条信号通路，与患者的临床分期、病理分化程度、放化疗敏感性及长期预后密切相关。然而，miRNA及其靶基因驱动NPC进展和导致不良预后的具体机制仍不清楚。此外，现有的大多数研究集中于长链非编码RNA（long non-coding RNA，lncRNA）-miRNA-mRNA三方网络，而基于miRNA-mRNA二元调控网络的预后模型开发相对不足。更重要的是，虽然已有包含免疫相关基因的NPC预后模型被报道，但这些模型通常仅依赖mRNA表达数据，忽略了上游miRNA对免疫相关基因的转录后调控，难以阐明免疫相关基因异常表达的深层驱动因素，也使具有潜在治疗价值的miRNA被排除在预后评估体系之外。因此，有必要通过整合miRNA和mRNA多组学表达数据，构建miRNA-mRNA调控网络，聚焦核心免疫相关基因，探究其在NPC进展中的作用和潜在分子机制，为NPC的预后评估和免疫治疗提供理论依据和新靶点。
二、研究设计与主要发现
本研究旨在通过生物信息学分析鉴定NPC的关键miRNA和免疫相关靶基因，构建miRNA-mRNA调控网络，筛选并验证核心免疫基因作为预后分子生物标志物。研究得出的主要结论包括：通过整合分析，鉴定出三个与NPC预后相关的关键miRNA（hsa-miR-30e-5p、hsa-miR-30a-5p和hsa-miR-20a-5p），以及17个核心免疫相关差异表达靶基因；这些基因显著富集于细胞因子介导的信号传导、趋化因子应答和JAK-STAT信号通路。基于LASSO-Cox回归构建的五基因（CXCL11、SOCS1、ROBO1、LGR5和RORB）预后风险模型在NPC中展现出良好的诊断效能，其中SOCS1的诊断准确性极高（AUC=0.935）。这五个关键基因的表达均与CD8+T细胞、巨噬细胞和中性粒细胞等多种免疫细胞的浸润水平呈显著正相关。药物预测分析显示5-氟尿嘧啶靶向LGR5，阿地白介素靶向SOCS1。体外RT-qPCR实验证实SOCS1在NPC细胞系中较正常鼻咽上皮细胞系显著上调约3.7倍（P<0.01）。该研究的科学意义在于：首次通过整合miRNA-mRNA调控网络与免疫相关基因分析，揭示了SOCS1等基因在NPC中的潜在重要作用，为NPC的精准预后评估和免疫治疗提供了新的候选生物标志物和治疗靶点。
三、关键技术方法概述
研究人员主要采用了以下关键技术方法：1）从公共数据库获取数据：从基因表达综合数据库（Gene Expression Omnibus，GEO）下载NPC相关的miRNA微阵列数据集GSE32960（含16例正常组织和312例非癌性鼻咽组织样本）以及mRNA微阵列数据集GSE12452（10例正常、31例肿瘤）、GSE13597（3例正常、25例肿瘤）和GSE64634（4例正常、12例肿瘤）；同时利用癌症基因组图谱-头颈鳞状细胞癌（The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma，TCGA-HNSC）数据集进行模型构建和验证。2）差异表达分析与批次校正：使用R软件的limma包和sva包对三个mRNA数据集进行合并和批次效应去除，并分别鉴定差异表达miRNA（differentially expressed miRNAs，DEMs）和差异表达mRNA（differentially expressed genes，DEGs）。3）靶基因预测与网络构建：联合TargetScan、miRDB和miRWalk三个数据库预测DEMs的靶基因，与DEGs取交集获得核心候选基因，并在Cytoscape软件中可视化miRNA-mRNA调控网络。4）免疫相关基因筛选与功能富集：从ImmPort数据库获取免疫相关基因集，与核心候选基因取交集得到核心免疫相关基因，并进行基因本体论（Gene Ontology，GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路富集分析。5）预后模型构建与评估：采用LASSO-Cox回归分析筛选预后相关基因并建立风险评分模型，通过Kaplan-Meier生存曲线和时间依赖性受试者工作特征（receiver operating characteristic，ROC）曲线评估模型性能。6）免疫浸润分析与药物预测：运用CIBERSORT算法量化免疫细胞浸润水平，并利用药物-基因相互作用数据库（Drug-Gene Interaction database，DGIdb）预测潜在靶向药物。
四、研究结果详述
3.1 数据下载结果： mRNA表达谱包含GSE12452、GSE13597和GSE64634三个微阵列数据集，共纳入17例正常组织和68例NPC组织样本；miRNA表达谱基于GSE32960数据集，包含16例正常组织和312例非癌性鼻咽组织。
3.2 miRNA数据集差异表达分析： 在GSE32960数据集中鉴定出49个差异表达miRNA。
3.3 mRNA数据集差异表达分析： 整合三个mRNA数据集后，共鉴定出1050个差异表达mRNA。
3.4 miRNA靶基因鉴定与候选基因筛选： 通过三个数据库的交集获得3704个高置信度靶基因，再与DEGs取交集得到213个共享基因；进一步与免疫相关基因集交叉比对，最终确定17个核心免疫相关差异表达靶基因。调控网络分析显示，hsa-miR-30e-5p、hsa-miR-30a-5p和hsa-miR-20a-5p具有最高的连接度。
3.5 GO和KEGG富集分析： 17个核心免疫相关基因显著富集于细胞因子介导的信号通路、趋化因子应答和JAK激活等生物学过程，以及细胞因子-受体相互作用和JAK-STAT信号通路等免疫相关通路。
3.6 PPI网络与共表达相关性分析： 蛋白-蛋白互作（protein-protein interaction，PPI）网络显示SOCS1、CCL2、PTGS2、CXCL11和RSAD2等基因作为枢纽节点；共表达热图显示大多数核心基因间存在显著正相关。
3.7 预后风险模型的构建与验证： 通过LASSO-Cox回归筛选出CXCL11、SOCS1、ROBO1、LGR5和RORB五个预后相关基因，并建立了风险评分公式。时间依赖性ROC曲线显示模型预测1年、3年和5年总生存期的曲线下面积（area under the curve，AUC）分别为0.574、0.622和0.636。Kaplan-Meier分析表明高风险组的总生存期显著短于低风险组（P<0.001，HR=1.65）。
3.8 免疫细胞浸润分析： NPC肿瘤微环境中存在广泛且异质的免疫细胞浸润，主要包括巨噬细胞、T细胞亚群和中性粒细胞。五个关键基因的表达水平与CD8+T细胞、巨噬细胞和中性粒细胞等多种免疫细胞的浸润水平呈显著正相关。
3.9 关键基因在训练集中的诊断效能评估： SOCS1表现出极高的诊断准确性（AUC=0.935），CXCL11诊断效能较高（AUC=0.809），RORB（AUC=0.673）和LGR5（AUC=0.654）表现中等，而ROBO1诊断准确性较低（AUC=0.473）。
3.10 关键基因在HNSC肿瘤与正常组织中的差异表达： 与正常组织相比，RORB在肿瘤组织中显著下调，而LGR5、ROBO1、SOCS1和CXCL11均显著上调。
3.11 关键基因对1年生存的预测价值： 五个关键基因单独预测1年总生存期的AUC均低于0.6，表明单个基因不宜作为独立的预后预测指标，进一步支持多基因风险评分模型在预后评估中的优势。
3.12 关键基因蛋白表达验证： 人类蛋白质图谱（Human Protein Atlas，HPA）数据库免疫组化结果显示，RORB和SOCS1在肿瘤组织中下调，LGR5上调，ROBO1表达模式出现特异性改变，而CXCL11未见显著差异。
3.13 药物-基因相互作用预测： 预测发现5-氟尿嘧啶靶向LGR5，阿地白介素靶向SOCS1，为NPC靶向治疗的药物筛选提供了初步依据。
3.14 RT-qPCR结果： 实时定量聚合酶链反应（RT-qPCR）结果显示，SOCS1在NPC细胞系中的mRNA表达水平较正常鼻咽上皮细胞系显著升高约3.7倍（P<0.01）。
五、讨论与结论
讨论部分指出，通过整合miRNA-mRNA调控网络与免疫相关基因分析，研究揭示了hsa-miR-30e-5p、hsa-miR-30a-5p和hsa-miR-20a-5p在NPC预后中的潜在作用，并鉴定出17个核心免疫相关差异表达靶基因，这些基因显著富集于细胞因子介导的信号转导、趋化因子应答和JAK-STAT信号通路，表明免疫调控机制在NPC进展中发挥关键作用。基于CXCL11、SOCS1、ROBO1、LGR5和RORB五个关键基因构建的预后模型在NPC中展现出良好的诊断效能。值得注意的是，SOCS1在NPC细胞中显著上调，这与传统认为的抑癌基因角色不同，提示其在NPC中可能具有组织特异性和功能异质性，可能与EB病毒（Epstein-Barr virus，EBV）感染诱导的反馈性上调有关，参与调节肿瘤微环境中的炎症反应和免疫抑制状态，帮助肿瘤细胞逃避免疫监视。研究也承认存在局限性：预后风险模型基于回顾性公共数据库队列，需前瞻性临床队列验证；缺乏体内动物实验验证关键基因的生物学功能；miRNA-mRNA调控网络仅基于生物信息学预测，仍需进一步的功能实验阐明精确的分子相互作用。
研究结论： 通过分析多个基因表达谱及TCGA-HNSC数据集，研究人员筛选出五个可能与NPC发展相关的基因（CXCL11、SOCS1、ROBO1、LGR5和RORB），这些基因在肿瘤进展和转移过程中发挥重要的免疫功能。研究结果还表明，SOCS1是NPC患者一个有前景的生物标志物，可作为预后和免疫治疗应答的预测指标。
