关键词:
激素性股骨头坏死
m^(6)A甲基化
微小RNA
转录因子
生物信息学
差异基因
基因调控网络
核心基因
摘要:
背景:m^(6)A修饰与股骨头坏死的发生发展相关,但在激素性股骨头坏死中的作用尚不清楚。目的:基于GEO数据库,采用生物信息学方法分析激素性股骨头坏死中表达差异的m^(6)A基因及互作miRNAs,探寻其潜在发病机制。方法:在GEO数据库中检索并下载与激素性股骨头坏死相关的mRNA表达谱数据集(GSE123568),通过R软件对数据集进行差异基因筛选及GO功能、KEGG通路富集分析。识别差异基因中的m^(6)A差异表达基因(m^(6)A-DEGs)并对其进行GO功能与KEGG通路富集分析,比较m^(6)A-DEGs的表达量并分析它们之间的相关性。最后通过Cytoscape构建m^(6)A-DEGs的PPI互作网络及筛选核心基因。使用TargetScan,miRTarBase和miRBD数据库预测m^(6)A-DEGs相关的潜在miRNAs,同时使用ChIPBase及hTFtarget数据库预测7个核心基因潜在转录因子,然后分别构建m^(6)A-miRNA与转录因子m^(6)A调控网络。最后使用数据集GSE74089验证7个核心m^(6)A-DEGs的表达水平。结果与结论:①从数据集中共筛选出2460个差异表达的基因,其中1455个上调,1005个下调。②从数据集中筛选出了14个m^(6)A-DEGs,包括3个下调和11个上调基因,m^(6)A-DEGs在激素性股骨头坏死中的表达具有显著差异(P<0.05),Spearman分析表明它们之间具有一定相关性。③m^(6)A-DEGs的GO和KEGG富集分析主要集中在骨髓细胞分化与发育、免疫受体与细胞因子受体活性、破骨细胞分化、AMPK与白细胞介素17信号通路。④m^(6)A-DEGs前7个核心基因包括YTHDF3,YTHDF1,YTHDF2,ALKBH5,METTL3,HNRNPA2B1及HNRNPC,它们在miRTarBase,miRDB和TargetScan数据库中共有44个miRNA重叠,在ChIPBase及hTFtarget数据库中共有79个重叠转录因子。⑤在GSE74089数据集中有6个核心m^(6)A-DEGs的表达水平与GSE123568数据集一致。⑥结果证实,根据生物信息学方法筛选的7个m^(6)A-DEGs可能通过调控多个miRNA、转录因子和AMPK及白细胞介素17信号通路表达,进而影响激素性股骨头坏死中骨髓细胞分化发育与破骨细胞分化,为进一步深入研究激素性股骨头坏死的发病机制和靶向治疗提供了数据支持和研究方向。