目的:构建miRNA-mRNA网络,识别调控骨肉瘤发生的关键miRNA。方法:下载GEO数据库中GSE28423和GSE28424两个骨肉瘤相关的mRNA和miRNA表达谱数据集,差异表达分析得到失调的mRNA和miRNA,并利用超几何检验获得失调mRNA富集的miRNA。从TargetScan数据库和miRTarBase数据库中获得miRNA与mRNA关系对,结合皮尔斯相关系数构建miRNA-mRNA调控网络。最后结合GSE16091和GSE39040数据集中的临床信息,对网络的全部节点进行生存分析;对hub节点进行功能注释,探索与骨肉瘤发病相关的miRNA和mRNA的调控关系。结果:通过对GSE28423及GSE28424表达谱的差异表达分析,我们得到了174个差异的miRNA及234个差异的mRNA。结合超几何检验、网络节点的预后分析和GO、KEGG功能注释,最终发现miR-371-5p、miR-139-5p、miR-151-3p与骨肉瘤患者预后相关,miR-371-5p可能通过靶向诱癌基因mRNA启动免疫相关途径进而影响骨肉瘤患者的生存,miR-139-5p则可能是通过影响细胞代谢途径来影响疾病发生,而miR-151-3p则可能通过影响免疫功能进而影响骨肉瘤的发生发展。结论:通过生物信息学分析,我们揭示了骨肉瘤发病过程的潜在的关键miRNA。随着免疫治疗的进一步成熟,miR-371-5p、miR-139-5p、miR-151-3p可能成为骨肉瘤免疫治疗的有效靶点。这些关键分子或可成为骨肉瘤预后判断的新指标。