MELK/TMPO-AS1/hsa-let-7b-5p网络在肺腺癌中的预后意义

《Frontiers in Oncology》：Prognostic significance of the MELK/TMPO-AS1/hsa-let-7b-5p network in lung adenocarcinoma

　　

摘要

引言：母系胚胎亮氨酸拉链激酶（MELK）是G2/M检查点的关键调控因子，也是一种公认的泛癌致癌基因；然而，其在肺腺癌（LUAD）中的调控机制及临床意义仍不完全清楚。本研究旨在探讨MELK在LUAD中的分子、预后、免疫及治疗相关性。

方法：采用整合多组学方法，包括基因表达、生存、转录组、免疫浸润、调控网络、分子对接和分子动力学分析。通过生物信息学方法研究了潜在的MELK相关ceRNA调控轴，并在肺癌细胞系中通过qRT-PCR进行了验证。进一步利用分子对接和分子动力学模拟探索候选MELK结合化合物的治疗潜力。

结果：MELK在LUAD中显著过表达（*log2FC = 4.19），并与较差的总生存期显著相关（HR = 1.63），其中在I期疾病患者（HR = 2.07）和女性吸烟者（HR = 1.50）中预后关联更强。MELK与FOXM1呈强正相关（R = 0.834），支持其在G2/M调控程序中的协同参与。高MELK表达与效应性免疫细胞浸润减少以及耗竭CD8+ T细胞和调节性T细胞富集增加相关。整合调控分析确定了推测的**TMPO-AS1/hsa-let-7b-5p/MELK/FOXM1** ceRNA网络，其特征为TMPO-AS1和MELK表达升高，hsa-let-7b-5p表达降低。这些表达模式在肺癌细胞系中的qRT-PCR分析中得到了进一步支持。分子对接和分子动力学模拟确定*橙皮苷*为候选MELK结合化合物，具有较优的预测结合亲和力和稳定的复合物行为。

结论：这些发现提供了MELK在LUAD中失调的整合视图，将其G2/M相关活性与转录后调控、免疫特征和潜在治疗靶向联系起来。*TMPO-AS1/hsa-let-7b-5p/MELK/FOXM1*轴可能代表了理解MELK介导的LUAD进展以及识别预后生物标志物和治疗机会的一个有前景的分子框架。

论文解读：MELK/TMPO-AS1/hsa-let-7b-5p网络在肺腺癌中的预后意义

一、研究背景与目的

肺癌是全球癌症死亡的主要原因之一，其中非小细胞肺癌（NSCLC）是最常见类型，包含肺腺癌（LUAD）和肺鳞癌（LUSC）。LUAD起源于远端气道黏液分泌细胞，而LUSC起源于支气管上皮。LUAD组织学亚型包括贴壁为主型（40–60%）、腺泡型（40–50%）和乳头型（5–10%），这些亚型通常侵袭性较低且预后较好；而实体为主型（20–30%）、黏液型（2–10%）、微乳头型（2–5%）等则更具侵袭性和转移性。尽管靶向治疗和免疫检查点抑制剂已取得进展，但许多癌症确诊时已处于晚期，治疗疗效受限，且存在毒性、耐药性和肿瘤异质性等问题。因此，深入了解肺癌进展的分子机制至关重要。

细胞周期失调是肺癌发生的一个标志性特征。母系胚胎亮氨酸拉链激酶（MELK）是G2/M检查点的关键调控因子，在胚胎和祖细胞中调节有丝分裂进程和G2/M检查点控制。MELK含有一个N端激酶结构域，介导其磷酸化依赖的调控功能，并与关键有丝分裂调控因子相互作用以协调细胞周期进程。在肿瘤尤其是LUAD中，MELK过表达，通过PLK1-CDC25C-CDK1轴促进增殖，破坏检查点，导致基因组不稳定和抗凋亡。MELK还与上皮-间质转化（EMT）和免疫逃逸相关。功能性MELK–FOXM1轴可能驱动持续增殖信号和肿瘤进展。FOXM1是G2/M程序的核心转录调控因子，控制有丝分裂进程所需的多个基因表达。在转录后水平，竞争性内源RNA（ceRNA）网络通过lncRNA和mRNA竞争microRNA来调控癌症进展。因此，本研究旨在系统探讨MELK在LUAD中的分子、预后、免疫和治疗相关性，并验证其相关ceRNA调控轴。

二、主要技术方法

研究人员采用整合多组学生物信息学方法，结合实验验证。具体包括：利用OncoMX、TIMER 2.0、GSCA、UALCAN、ENCORI、OncoDB、GEPIA2、Lung Cancer Explorer（LCE）等多个数据库进行MELK表达和预后分析；使用Kaplan-Meier plotter（KMP）进行生存分析；采用CancerSEA分析MELK参与的细胞过程；通过UALCAN和TIMER 2.0分析TP53突变相关性；利用STRING构建蛋白-蛋白相互作用（PPI）网络；通过SIGNOR数据库识别上游转录因子；使用TIMER 2.0和GSCA分析免疫浸润；利用CellTracer和Monocle 2进行单细胞转录组分析和拟时间轨迹分析；采用AutoDock Vina进行分子对接，GROMACS进行分子动力学（MD）模拟；通过miRNet、ENCORI、miRWalk、RNA22v2等鉴定miRNA和lncRNA；最后在肺癌细胞系（A549、H460、H23）和正常肺成纤维细胞系（MRC-5）中通过qRT-PCR验证TMPO-AS1、hsa-let-7b-5p和MELK的表达。细胞系来源：A549和H460来自AIIMS Bhopal，H23来自印度浦那国家细胞科学中心（NCCS），MRC-5来自印度新德里AIIMS。

三、研究结果

3.1 MELK功能注释与泛癌分析

MELK基因定位于染色体9p13.2，编码丝氨酸/苏氨酸激酶。GO功能富集分析显示MELK参与G2/M转换、蛋白磷酸化、染色质重塑等多个关键生物学过程。OncoMX数据库分析显示MELK在多种癌症中显著上调，其中肺癌表达最高（log2FC = 4.19，P = 2.25e-267）。TIMER、GSCA和UALCAN数据库进一步验证了MELK在肺肿瘤组织中一致显著上调。

3.2 MELK在肺癌患者中的预后意义

KMP生存分析显示，高MELK表达患者的总生存期（OS）显著缩短（HR = 1.63，95% CI = 1.44-1.84，P = 1.1e-15），中位生存期从95个月降至47.6个月。UALCAN确认了这一趋势（HR = 1.776，P = 0.000526）。高MELK表达也与较短的首进展生存期（FP）相关（HR = 1.82，P = 3.9e-12）。分层分析显示，在LUAD中高MELK表达与较差的OS显著相关（HR = 1.58，P = 1.7e-07），而在LUSC中无显著关联。I期LUAD患者中，高MELK表达与中位生存期降低相关（HR = 2.07，P = 0.00035）。吸烟者中MELK高表达与不良生存相关（HR = 1.54，P = 0.001），其中女性吸烟者关联更强（HR = 1.50，P = 0.042）。

3.3 MELK在LUAD中的表达谱及其临床病理学特征

多个TCGA数据库（UALCAN、ENCORI、OncoDB、GEPIA2、LCE）及R语言验证均显示MELK在LUAD肿瘤组织中显著高表达。TCGAnalyzeR火山图显示MELK转录本为上调mRNA（logFC = 3.38）。MELK表达随肿瘤分期进展而升高，且与淋巴结转移状态正相关。LCE数据库显示女性吸烟者MELK水平高于非吸烟女性（P = 0.003）。

3.4 MELK在LUAD中作为肿瘤进展、EMT和TP53相关增殖的调控因子

CancerSEA数据库显示MELK与细胞周期（R = 0.57）、增殖（R = 0.54）、DNA修复（R = 0.48）、DNA损伤反应（R = 0.43）、侵袭（R = 0.37）和EMT（R = 0.27）显著相关。TP53突变的LUAD样本中MELK表达显著高于非突变样本。MELK与G2/M细胞周期调控因子CDK1（R = 0.865）、CCNB1（R = 0.840）、PLK1（R = 0.836）和CDC25C（R = 0.843）呈强正相关。MELK与EMT标志物TWIST1（R = 0.53）和SNAI1（R = 0.33）正相关。TNMplot分析显示MELK表达从正常到肿瘤到转移样本逐渐升高（P = 2.9e-134）。

3.5 LUAD中MELK异质性模型及关键转录因子的鉴定

通过UALCAN和ENCORI鉴定出与MELK正相关（CKAP2L、KIF4A等）和负相关（NAPSA、TMEM125等）的基因。STRING分析显示MELK与FOXM1、CDC20、BUB1、CCNB1等细胞周期基因存在紧密互作网络。SIGNOR数据库识别出包括FOXM1在内的多个候选调控因子。FOXM1在所有数据库中与MELK均呈高正相关（R > 0.8）。FOXM1高表达与LUAD患者较差的OS显著相关（HR = 2.11，P < 1e-16）。FOXM1表达随肿瘤分期、淋巴结转移和吸烟状态升高。

3.6 MELK相关的肿瘤免疫微环境

Spring viewer显示MELK和FOXM1在多种免疫细胞群中表达。ROC分析显示MELK（AUC = 0.693）和FOXM1（AUC = 0.689）对免疫反应具有中等预测价值。TIMER 2.0分析显示MELK与CD8+ T细胞（R = -0.231）、CD4+ T细胞（R = -0.334）、B细胞（R = -0.337）、巨噬细胞（R = -0.39）和树突状细胞（R = -0.327）呈负相关。单细胞分析显示MELK在CD8+耗竭T细胞、恶性细胞、浆细胞和Treg细胞中富集。拟时间轨迹分析进一步验证了这些细胞群中MELK的高表达。

3.7 MELK潜在治疗靶点

分子对接显示OTSSP167（MELK抑制剂）结合亲和力为-9.672 kcal/mol。多种G2/M检查点靶向药物（Onvansertib、Volasertib、Adavosertib等）和微管靶向药物均与MELK有较好结合。天然黄酮类化合物中，橙皮苷（Hesperidin）结合亲和力最高（-10.050 kcal/mol），其次为Eriocitrin（-10.030 kcal/mol）、Narirutin（-9.866 kcal/mol）等。

3.8 分子动力学模拟证实MELK-配体复合物稳定

100 ns分子动力学模拟显示MELK-OTSSP167和MELK-橙皮苷复合物均保持结构稳定。RMSD、RMSF、回转半径（Rg）、溶剂可及表面积（SASA）和氢键分析表明橙皮苷与MELK结合稳定，与OTSSP167行为相似。

3.9 LUAD中靶向MELK和FOXM1的调控miRNA的鉴定与筛选

miRNet预测31个共同靶向MELK和FOXM1的miRNA，其中7个在LUAD中显著下调。ENCORI相关性分析显示hsa-let-7b-5p与MELK（R = -0.360）和FOXM1（R = -0.393）呈负相关。生存分析显示高hsa-let-7b-5p表达与较好的OS相关（HR = 0.72，P = 0.026）。ROC分析显示hsa-let-7b-5p具有高诊断潜力（AUC = 0.89）。hsa-let-7b-5p在LUAD肿瘤组织中显著下调，且随病理分期进展逐渐降低。

3.10 hsa-let-7b-5p相关lncRNA的鉴定

Enrichr数据库预测97个与hsa-let-7b-5p相互作用的lncRNA，其中TMPO-AS1在LUAD中显著上调，且与较差的OS相关（HR = 2.16，P = 4.1e-10）。ENCORI显示TMPO-AS1与hsa-let-7b-5p呈负相关（R = -0.217，P = 1.20e-05），与MELK（R = 0.601）和FOXM1（R = 0.513）呈正相关。TMPO-AS1表达随肿瘤分期升高，且在女性吸烟者中显著升高。RNA22v2预测了TMPO-AS1与hsa-let-7b-5p之间的稳定结合（折叠能-13.90 kcal/mol）。构建了TMPO-AS1/hsa-let-7b-5p/MELK-FOXM1 ceRNA调控网络。

3.11 细胞系模型中TMPO-AS1/hsa-let-7b-5p/MELK轴的表达谱

qRT-PCR结果显示，与正常肺成纤维细胞MRC-5相比，肺癌细胞系（A549、H460、H23）中TMPO-AS1表达显著上调（H460 >20倍，A549约10-15倍，H23约4-8倍），hsa-let-7b-5p表达显著下调（H23下调最明显，约10倍），MELK表达显著上调（A549和H460约4-5倍，H23约3-4倍）。这些结果验证了生物信息学预测的调控关系。

四、讨论与结论总结

讨论部分指出，LUAD是一种常见且致命的恶性肿瘤，常因诊断较晚导致预后不良。MELK作为LUAD增殖的关键驱动因子，通过调控PLK1-CDC25C-CDK1轴支持G2/M进程。研究发现MELK与FOXM1在LUAD中强关联，尤其是在晚期和吸烟相关肿瘤中。高MELK表达与效应免疫细胞浸润减少及耗竭CD8+ T细胞和Treg富集相关，支持肿瘤免疫逃逸。天然黄酮类化合物（尤其是橙皮苷）与MELK结合稳定，可能为联合治疗提供依据。ceRNA网络分析确定了TMPO-AS1/hsa-let-7b-5p/MELK-FOXM1调控轴，其中TMPO-AS1可能通过海绵吸附hsa-let-7b-5p来解除其对MELK和FOXM1的抑制，从而促进LUAD进展。

研究结论：该研究强调MELK是肺腺癌（LUAD）中一个关键的致癌和免疫调节驱动因子，并表明其受FOXM1-TMPO-AS1/hsa-let-7b-5p ceRNA轴调控。该网络的失调促进细胞周期进程和上皮-间质转化（EMT），形成免疫抑制微环境。MELK既可作为预后生物标志物，也可作为治疗靶点。提出的干预策略包括增强hsa-let-7b-5p活性以抑制MELK和FOXM1，利用siRNA和反义寡核苷酸靶向TMPO-AS1，以及将天然MELK抑制剂（如橙皮苷）与选择性抑制剂（如OTSSP167）联合使用以改善激酶阻断效果。这些策略揭示了ceRNA轴的治疗潜力，并为未来LUAD的研究提供了基础。