类风湿关节炎患者滑膜组织lncRNA表达及基于CeRNA网络的生物信息学分析
作者: |
1肖剑伟,
1蔡旭,
1陈新鹏,
1曾苗雨,
1郭粉莲,
1洪易炜,
1叶志中
1 深圳市福田区风湿病专科医院风湿免疫科,广东 深圳 518000 |
通讯: |
叶志中
Email: yezhizhong0823@163.com |
DOI: | 10.3978/j.issn.2095-6959.2020.05.011 |
基金: | 深圳市医疗卫生三名工程(SZSM201602087);深圳市福田区卫生公益性科研项目(FTWS20180)。 |
摘要
目的:通过生物信息学分析来识别类风湿关节炎(rheumatoid arthritis,RA)滑膜组织病变进展相关的差异表达基因。方法:通过NCBI GEO数据库获取GSE55235和GSE55457的基因表达谱。采用Perl语言对下载的数据进行样本数据合并及基因重注释;采用R语言进行批次矫正及差异分析,根据差异长链非编码RNA(lncRNA)和mRNA构建竞争性内源RNA(ceRNA)网络及进行GO富集分析和KEGG通路分析;使用cytoHubba插件筛选Hub基因,分析与差异LncRNA的相关性。结果:分析显示与正常滑膜组织对比,RA患者滑膜组织143个mRNA、3个lncRNA存在明显差异表达。根据差异基因构建lncRNA-miRNA-mRNA互作网络,网络由2个LncRNA节点,16个miRNA节点、17个mRNA节点以及44个边组成。GO功能富集分析主要集中在细胞死亡的正调控、成纤维细胞增殖的调节、免疫应答调节细胞表面受体信号通路等功能。KEGG通路分析显示35条通路被富集,其中涉及IL-17代谢通路、MAPK信号通路、WNT信号通路、TNF信号通路等。其中Hub基因MYC,CDKN1A,JUN,FOS与LncRNA MEG3在RA滑膜组织中均呈低表达,且lncRNA MEG3与MYC,CDKN1A,JUN,FOS表达具有相关性。结论:通过生物信息学网络分析,lncRNA MEG3可能作为ceRNA在RA的疾病发展中发挥着重要作用,为RA提供一些新的候选诊断生物标志物或潜在的治疗靶点。
关键词:
类风湿关节炎;生物信息学分析;lncRNA;ceRNA;Perl语言
Expression of LncRNA in synovial tissues of patients with rheumatoid arthritis and bioinformatics analysis based on CeRNA network
CorrespondingAuthor: YE Zhizhong Email: yezhizhong0823@163.com
DOI: 10.3978/j.issn.2095-6959.2020.05.011
Foundation: This work was supported by the Sanming Project of Medicine in Shenzhen (SZSM201602087) and Shenzhen Futian Medical Science and Technology Research Fund (FTWS2018066), China.
Abstract
Objective: To identify differentially expressed genes related to the progression of synovial tissue lesions in rheumatoid arthritis (RA) by bioinformatics analysis. Methods: Gene expression profiles of GSE55235 and GSE55457 were acquired from the NCBI GEO database. Data were combined and re-annotated in Perl; Batch correction and identify differentially expressed genes were performed using R; the competing endogenous RNA (ceRNA) network was constructed based on the differential long non-coding RNA (lncRNA) and mRNA. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed. The Hub gene was screened using the cytoHubba plugin and analyzed for the correlation between differential lncRNA. Results: Compared with normal synovial tissue, 143 mRNA and 3 LncRNA in synovial tissue of RA patients were significantly different. The LncRNA-miRNA-mRNA interaction network was constructed based on the differential gene. The network consisted of two LncRNA nodes, 16 miRNA nodes, 17 mRNA nodes and 44 edges. GO functional enrichment analysis mainly were significantly enhanced in positive regulation of programmed cell death, regulation of fibroblast proliferation and immune response-regulating cell surface receptor signaling pathway. KEGG pathway analysis showed that pathways associated with IL-17 signaling pathway, MAPK signaling pathway, Wnt signaling pathway, and TNF signaling pathway. Hub gene MYC, CDKN1A, JUN, FOS and LncRNA MEG3 were down-regulated in RA synovial tissue, and LncRNA MEG3 was correlated with the expression of MYC, CDKN1A, JUN and FOS. Conclusion: The bioinformatics network analysis showed that LncRNA MEG3 may play an important role in the development of RA in ceRNA and providing some new candidate diagnostic biomarkers or potential therapeutic targets for RA.
Keywords:
rheumatoid arthritis; bioinformatics analysis; lncRNA; ceRNA; Perl