《Acta Neurologica Scandinavica》：Integrative MicroRNA Expression Profiling Identifies hsa-miR-210-3p and hsa-miR-8075 as Candidate miRNAs Associated With Meningioma Biology
编辑推荐：
背景
脑膜瘤是最常见的原发性颅内肿瘤，临床结局具有异质性。由于当前诊断工具缺乏足够的分子分辨率，亟需新的生物标志物以改进疾病表征。微小RNA（miRNA）作为关键的转录后调控因子，已成为肿瘤生物学中的重要参与者。然而，其在脑膜瘤中的作用仍缺乏充分研究。
方
背景
脑膜瘤是最常见的原发性颅内肿瘤，临床结局具有异质性。由于当前诊断工具缺乏足够的分子分辨率，亟需新的生物标志物以改进疾病表征。微小RNA（miRNA）作为关键的转录后调控因子，已成为肿瘤生物学中的重要参与者。然而，其在脑膜瘤中的作用仍缺乏充分研究。
方法
研究人员对两个公开可用的miRNA表达数据集（GSE126563和GSE88721）进行了整合分析，以鉴定原发性脑膜瘤中持续失调的miRNA。利用4个经典数据库——miRWalk、miRTarBase、TargetScan和miRDB——预测变化最显著miRNA的高置信度靶基因。随后使用GSE43290数据集在转录组水平对预测靶基因进行进一步评估，并采用g:Profiler实施功能富集分析。
结果
在两个数据集中共鉴定出30个持续失调的miRNA。其中，hsa-miR-210-3p是显著性最高的上调miRNA，而hsa-miR-8075是唯一持续下调的miRNA。研究分别为hsa-miR-210-3p和hsa-miR-8075预测到29个和125个靶基因，且验证结果显示其在GSE43290中呈现相反的表达模式。功能富集分析提示，hsa-miR-210-3p的靶基因参与神经元配体转录以及NF-κB信号传导，而hsa-miR-8075的靶基因则富集于突触信号传导和蛋白结合相关通路。
结论
研究结果鉴定出hsa-miR-210-3p和hsa-miR-8075为具有相反表达模式、并可能在脑膜瘤中具有调控相关性的候选miRNA。这些miRNA是值得在脑膜瘤中进一步开展机制学与临床验证的有前景候选分子。
本文发表于《Acta Neurologica Scandinavica》，聚焦脑膜瘤（meningioma）分子标志物不足这一关键问题。脑膜瘤是最常见的原发性颅内肿瘤，虽然多数生长缓慢且组织学上偏良性，但部分病例可表现出复发、侵袭性生长以及向更高世界卫生组织分级进展等不良生物学行为。现有分类体系主要依赖组织病理学标准，尚不能充分解释患者间显著的生物学异质性与临床异质性，因此亟需能够提升诊断精度、辅助预后判断并服务个体化治疗的分子标志物。微小RNA（miRNA，长度较短的非编码RNA，通过与靶mRNA互补结合介导转录后抑制）已被证实在多种肿瘤发生发展中发挥重要作用，并参与增殖、凋亡、血管生成和免疫逃逸等过程，但脑膜瘤相关miRNA研究仍较为有限，尤其缺乏跨队列、可重复的整合性分析。基于此，研究人员开展本研究，旨在通过整合公开miRNA表达谱和外部转录组数据，筛选在原发性脑膜瘤中稳定失调的候选miRNA，并解析其潜在靶基因及功能通路。
研究人员首先从GEO数据库纳入两个miRNA表达谱队列：GSE126563与GSE88721，用于比较脑膜瘤样本与对应对照样本间的差异表达miRNA；同时纳入mRNA转录组队列GSE43290，用于对候选miRNA靶基因进行外部转录组支持分析。该设计的核心价值在于，不依赖单一数据集或单个候选分子，而是通过跨数据集交叉验证提高结果的稳健性。研究的总体结论是：在两个独立队列中共发现30个持续失调的miRNA，其中hsa-miR-210-3p为最显著上调miRNA，hsa-miR-8075为唯一持续下调miRNA；二者的预测靶基因在独立mRNA数据集中呈现与miRNA表达方向相反的变化趋势，提示其具有潜在生物学调控意义；功能富集结果进一步显示，hsa-miR-210-3p相关靶基因主要涉及神经元配体转录调控、NF-κB信号传导及线粒体相关调节，而hsa-miR-8075相关靶基因则涉及突触信号传导、化学突触传递和蛋白结合网络。研究因此提出，hsa-miR-210-3p与hsa-miR-8075是脑膜瘤中值得进一步机制学和临床验证的候选miRNA。
在主要技术方法方面，作者基于GEO公开队列开展二次生物信息学分析：使用GEOquery获取GSE126563、GSE88721和GSE43290表达矩阵及平台注释；以“limma”进行差异表达分析，阈值设定为调整后p值<0.01且|log2FC|≥1；对两个miRNA队列取交集筛选稳定失调miRNA；针对重点miRNA，整合miRWalk、miRTarBase、TargetScan和miRDB四库进行靶基因预测，仅保留至少3个数据库共同支持的高置信度靶标；再利用GSE43290中47例脑膜瘤与4例正常脑膜对照的mRNA数据进行反向表达验证，并通过g:Profiler完成GO、KEGG和Reactome富集分析。
以下为论文主体结果与内容解读。
3.1. Differential Expression of miRNAs in Meningioma
在这一部分，研究人员首先对两个独立脑膜瘤miRNA队列进行差异表达分析。结果显示，GSE126563中共有1102个显著差异表达miRNA，其中上调856个、下调246个；GSE88721中共有230个显著差异表达miRNA，其中上调132个、下调98个。尽管两个队列的样本来源、平台和规模不同，交叉分析后仍识别出30个在两队列中方向一致且统计学显著的重叠差异miRNA，说明这些分子具有较好的跨平台和跨样本重复性。
在这30个共享miRNA中，hsa-miR-210-3p是最显著的上调分子，在GSE126563中的log2FC为5.45，调整后p值为1.26 × 10
?3；在GSE88721中的log2FC为7.09，调整后p值为1.30 × 10
?3。此外，hsa-miR-4443、hsa-miR-937-5p、hsa-miR-324-5p和hsa-miR-744-5p也呈现强烈上调。与之相对，hsa-miR-8075是唯一一个在两个队列中均持续下调的miRNA，在GSE126563中的log2FC为?2.35，在GSE88721中的log2FC为?3.62。研究人员据此将hsa-miR-210-3p和hsa-miR-8075选定为后续整合分析的核心对象。该结果表明，脑膜瘤中存在可重复识别的miRNA失调谱，而这两种表达方向相反的miRNA可能代表不同的调控轴。
3.2. Integrative Prediction of High-Confidence Target Genes for hsa-miR-210-3p and hsa-miR-8075
为提高靶基因预测的可靠性，研究人员对hsa-miR-210-3p和hsa-miR-8075分别整合4个经典miRNA靶基因数据库进行交叉筛选。对于hsa-miR-210-3p，共有172个基因被至少2个数据库预测到，而其中29个基因被至少3个数据库共同支持，构成高置信度靶基因集合；其中9个基因被4个数据库同时识别，属于证据最集中的候选靶标。对于hsa-miR-8075，共筛得125个至少被3个数据库支持的高置信度靶基因，其中BNIP2是唯一被4个数据库一致预测的基因。该部分研究说明，借助多数据库整合能够显著减少单一算法造成的假阳性，提高候选miRNA—靶基因关联的可信度，也为后续外部转录组验证奠定基础。
3.3. Validation of the Predicted miRNA Targets Using Independent Transcriptomic Data
在独立mRNA表达队列GSE43290中，研究人员进一步检验预测靶基因是否符合miRNA与mRNA之间预期的负向调控关系。该数据集包含47例原发性脑膜瘤样本与4例正常脑膜对照。分析22,283个基因后发现，11,185个基因在脑膜瘤中下调，1,089个基因上调，提示脑膜瘤存在广泛的转录组失衡。
对hsa-miR-210-3p而言，研究人员在其预测靶基因中发现12个基因在GSE43290中显著下调，尤其包括ISCU、SDHD、MID1IP1和DIMT1。这一结果与hsa-miR-210-3p在miRNA层面显著上调相一致，支持其可能对这些基因产生抑制作用。对hsa-miR-8075而言，其预测靶基因中有8个在脑膜瘤中显著上调，尤其包括UNC5C、CANX和DAZAP1。这一结果与hsa-miR-8075持续下调相一致，提示其下调后可能导致对应mRNA去抑制（derepression）。这一部分的核心意义在于，研究不仅停留在计算预测阶段，而是通过独立转录组证据增强了候选调控关系的生物学可信度。
3.4. Functional Enrichment Analysis of the Target Genes
在功能层面，研究人员对两种候选miRNA的高置信度靶基因分别开展富集分析。hsa-miR-210-3p的靶基因显著富集于神经元调控相关通路，尤其是Reactome中的“MECP2 regulates transcription of neuronal ligands”，涉及BDNF和SIN3A等基因；转录因子分析识别出RelA，涉及BDNF、HIF3A和E2F3等关键基因，提示其可能与NF-κB信号轴存在联系；共靶标分析还支持hsa-miR-210-3p本身及hsa-miR-147b相关通路可能共同调节ISCU、NDUFA4和SDHD等线粒体相关基因。由此可见，hsa-miR-210-3p可能参与脑膜瘤中的缺氧适应、线粒体代谢调控及神经元样信号调节。
相较之下，hsa-miR-8075的靶基因主要富集于“protein binding”“synaptic signaling”和“chemical synaptic transmission”等条目，涉及KIT、ADCY1、SQSTM1和PDYN等基因，提示该miRNA相关网络可能更多影响神经元通信样过程及蛋白互作网络。该发现表明，脑膜瘤相关miRNA异常不仅涉及经典肿瘤信号，还可能与神经样调控程序和脑组织微环境相关分子网络发生交叉。
在讨论部分，作者强调本研究的创新性主要体现在整合设计：并非基于单一队列或个别候选miRNA，而是通过独立公开数据集的交叉筛选、外部转录组验证和功能通路富集，提升了候选miRNA识别的可重复性与生物学相关性。作者指出，hsa-miR-210-3p作为缺氧诱导型miRNA，在脑膜瘤中的显著升高提示其可能参与肿瘤微环境下的低氧适应，并通过抑制ISCU、SDHD、DIMT1等基因影响线粒体代谢和应激反应。hsa-miR-8075虽然文献报道有限，但其在两个队列中稳定下调，并伴随多个靶基因上调，提示其可能具有肿瘤抑制样作用。作者还指出，miRNA可在脑脊液和血液中检测，因此具有发展为微创诊断标志物的潜力。
同时，作者也客观说明了研究局限性：本研究属于公开微阵列数据的二次分析，可能受到平台差异、样本异质性和数据质量影响；缺乏同一样本配对的miRNA—mRNA数据及实验验证，因此尚不能直接证明调控因果关系；缺少基于体液的临床样本，也限制了结果的即时转化价值。尽管如此，研究仍为脑膜瘤miRNA调控网络的解析提供了较强的数据基础。
研究结论部分可译为：总之，本研究在独立数据集中鉴定出hsa-miR-210-3p和hsa-miR-8075这两个与脑膜瘤具有可重复关联的候选miRNA。整合性靶基因预测与转录组评估支持了它们潜在的生物学相关性，并凸显其作为后续深入研究对象的前景。这些发现为未来在脑膜瘤领域开展功能学研究和临床验证提供了有价值的基础。
