Ren Conglin, Li Mingshuang, Zheng Yang, Wu Fengqing, Du Weibin, Quan Renfu
The Third Clinical Medical College of Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China.
The First Affiliated Hospital of Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China.
PeerJ. 2021 May 14;9:e11427. doi: 10.7717/peerj.11427. eCollection 2021.
The pathogenesis of rheumatoid arthritis (RA) is complex. This study aimed to identify diagnostic biomarkers and transcriptional regulators that underlie RA based on bioinformatics analysis and experimental verification.
We applied weighted gene co-expression network analysis (WGCNA) to analyze dataset GSE55457 and obtained the key module most relevant to the RA phenotype. We then conducted gene function annotation, gene set enrichment analysis (GSEA) and immunocytes quantitative analysis (CIBERSORT). Moreover, the intersection of differentially expressed genes (DEGs) and genes within the key module were entered into the STRING database to construct an interaction network and to mine hub genes. We predicted microRNA (miRNA) using a web-based tool (miRDB). Finally, hub genes and vital miRNAs were validated with independent GEO datasets, RT-qPCR and Western blot.
A total of 367 DEGs were characterized by differential expression analysis. The WGCNA method divided genes into 14 modules, and we focused on the turquoise module containing 845 genes. Gene function annotation and GSEA suggested that immune response and inflammatory signaling pathways are the molecular mechanisms behind RA. Nine hub genes were screened from the network and seven vital regulators were obtained using miRNA prediction. CIBERSORT analysis identified five cell types enriched in RA samples, which were closely related to the expression of hub genes. Through ROC curve and RT-qPCR validation, we confirmed five genes that were specific for RA, including CCL25, CXCL9, CXCL10, CXCL11, and CXCL13. Moreover, we selected a representative gene (CXCL10) for Western blot validation. Vital miRNAs verification showed that only the differences in has-miR-573 and has-miR-34a were statistically significant.
Our study reveals diagnostic genes and vital microRNAs highly related to RA, which could help improve our understanding of the molecular mechanisms underlying the disorder and provide theoretical support for the future exploration of innovative therapeutic approaches.
类风湿关节炎(RA)的发病机制复杂。本研究旨在基于生物信息学分析和实验验证,确定RA潜在的诊断生物标志物和转录调节因子。
我们应用加权基因共表达网络分析(WGCNA)来分析数据集GSE55457,并获得与RA表型最相关的关键模块。然后我们进行了基因功能注释、基因集富集分析(GSEA)和免疫细胞定量分析(CIBERSORT)。此外,将差异表达基因(DEG)与关键模块内的基因的交集输入STRING数据库,以构建相互作用网络并挖掘枢纽基因。我们使用基于网络的工具(miRDB)预测微小RNA(miRNA)。最后,通过独立的GEO数据集、RT-qPCR和蛋白质印迹法对枢纽基因和重要miRNA进行验证。
通过差异表达分析共鉴定出367个DEG。WGCNA方法将基因分为14个模块,我们重点关注包含845个基因的蓝绿色模块。基因功能注释和GSEA表明免疫反应和炎症信号通路是RA背后的分子机制。从网络中筛选出9个枢纽基因,并通过miRNA预测获得7个重要调节因子。CIBERSORT分析确定了RA样本中富集的5种细胞类型,它们与枢纽基因的表达密切相关。通过ROC曲线和RT-qPCR验证,我们确认了5个RA特异性基因,包括CCL25、CXCL9、CXCL10、CXCL11和CXCL13。此外,我们选择了一个代表性基因(CXCL10)进行蛋白质印迹验证。重要miRNA验证表明,只有has-miR-573和has-miR-34a的差异具有统计学意义。
我们的研究揭示了与RA高度相关的诊断基因和重要微小RNA,这有助于提高我们对该疾病潜在分子机制的理解,并为未来探索创新治疗方法提供理论支持。