基于Illumina高通量测序技术研究肌少-骨质疏松症与骨质疏松症患者骨组织microRNAs差异表达谱及差异分析

目的 采用高通量测序技术筛选“肌少-骨质疏松症”（sarco-osteoporosis，SO）与骨质疏松症（osteoporosis，OP）患者骨组织标本中miRNA差异表达谱，并对差异显著的miRNA进行生物信息分析。方法 运用Illumina Hiseq 2500测序技术筛查2017年至2019年期间福建省老年医院SO组与OP组患者的6对骨组织样品中miRNA表达谱；差异显著表达的miRNA运用层次聚类图和火山图分析，并结合基因本体论（GO）富集分析、京都基因与基因组百科全书 (KEGG)信号通路分析探讨miRNA参与成肌细胞、成骨细胞增殖和分化的生物过程。 结果 SO组与OP 组进行组间对比，12份骨组织中共鉴定到 3 064个（P<0.05；| log2 Fold Change |>0.0）差异表达的miRNAs，其中上调的1 432个，下调的1 632个，筛选出其中差异表达显著的22个miRNAs︱log2 Fold Change︱≥2 且P <0.05包括15个在SO组骨组织中表达上调的 miRNAs 和 7个表达下调的 miRNAs；层次聚类图和火山图分析提示，SO和OP骨组织中的miRNA差异表达谱具有统计学意义；GO 分析结果显示，miRNA的靶基因富集于成肌细胞与成骨细胞的生物过程，还参与细胞的氧化代谢、增殖和分化以及凋亡等生物过程。KEGG通路分析显示，hsa-miR-382-3p、hsa-miR-27a-5p、hsa-miR-1226-5p、hsa-miR-3934-5p和hsa-miR-451a的靶基因集合显著富集于NF-κB、Wnt、MAPK和P53 等信号通路中。结论 Illumina高通量测序技术能有效筛查出肌少-骨质疏松症骨组织中差异表达的miRNA；通过分析5个差异显著基因的下游靶基因，并结合分子生物学相关数据库与查阅大量文献发现这些差异miRNA可能参与肌少-骨质疏松症疾病状态下成骨细胞和成肌细胞多种生物活性过程。这一研究结果为后续深入研究miRNA及其靶基因调控机体骨髓间充质干细胞增殖与成骨分化的机制提供了新思路和理论依据。

The differential expression profile and difference analysis of microRNAs in bone tissue of patients with sarco-osteoporosis and osteoporosis based on Illumina high-throughput sequencing technology