目的采用RNA测序技术分析成纤维样滑膜细胞(FLSs)中与类风湿关节炎(RA)相关的核心基因。方法(1)将30只SD大鼠随机分为模型组和正常组,各15只。给予模型组大鼠注射弗氏完全佐剂以构建佐剂性关节炎大鼠模型(经典RA模型),给予正常组大鼠注射等量无菌生理盐水。分别采用HE染色法、番红固绿染色法观察两组大鼠膝关节组织形态学变化和软骨组织破坏情况。(2)获取两组大鼠滑膜组织以分离培养FLSs。采用RNA测序技术分析两组FLSs的基因表达情况后,利用R软件筛选两组FLSs中的差异表达基因。对差异表达基因进行基因本体论功能富集分析、京都基因与基因组百科全书通路富集分析、基因集富集分析。通过STRING数据库和Cytoscape软件建立差异表达基因的蛋白-蛋白相互作用(PPI)网络并筛选核心基因。采用实时荧光定量PCR对部分差异表达基因进行验证。结果HE染色和番红固绿染色结果表明造模成功。共筛选到227个差异表达基因,包括133个上调基因与94个下调基因。差异表达基因主要富集于平滑肌细胞增殖负向调控等生物过程,脂肪颗粒等细胞组分,钙离子结合等分子功能,以及脂肪细胞中脂肪分解的调节等信号通路;差异表达基因富集的基因集主要涉及防卫反应、G蛋白偶联受体信号通路及活动分子传感器。共筛选出10个核心基因,包括MMP9、AGT、CECR2、FGR、Hist1h3c、APOB、CDH3、MMP15、Hist2h4、NPPB。实时荧光定量PCR结果显示,模型组大鼠滑膜组织中MMP9、FGR、APOB的mRNA表达量高于正常组(P<0.05),与RNA测序结果相符。结论MMP9、FGR、APOB等核心基因可能通过参与调控免疫炎症反应来参与RA的发生和发展。
ObjectiveTo analyze the core genes associated with rheumatoid arthritis (RA) in fibroblast-like synoviocytes (FLSs) by employing RNA sequencing technique. Methods(1) Thirty SD rats were randomly divided into model group or normal group, with 15 rats in each group. Rats in the model group received injection of complete Freund′s adjuvant to establish the rat model of adjuvant arthritis (classic RA model), and rats in the normal group received injection of sterile normal saline with equivalent volume. The histomorphological changes and cartilage tissue destruction of rats′ knee joint in the two groups were observed by the HE staining and saffron O-fast green staining, respectively. (2) Synovial tissues of rats in the two groups were obtained for FLSs isolation and culture. After analyzing the gene expressions of FLSs in the two groups by the RNA sequencing technique, the R software was used to screen the differentially expressed genes in FLSs of the two groups. The Gene Ontology functional enrichment analysis, Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis, and Gene Set Enrichment Analysis were performed on differentially expressed genes. The protein-protein interaction (PPI) network of differentially expressed genes was constructed and core genes were screened by using the STRING database and Cytoscape software. The real-time fluorescent quantitative PCR was used to validate some differentially expressed genes. ResultsThe results of the HE staining and saffron O-fast green staining revealed that the modeling was successful. A total of 227 differentially expressed genes were screened out, including 133 up-regulated genes and 94 down-regulated genes. Differentially expressed genes were mainly concentrated in biological processes such as negative regulation of smooth muscle cell proliferation, etc.,in cellular compsitions such as lipid particles, etc.,in molecular functions such as calcium ion binding, etc., and in signaling pathways such as regulation of lipolysis in adipocytes, etc. The gene sets enriched by differentially expressed genes were mainly involved in defense response, G protein-coupled receptor signaling pathway, and active molecular transducer. A total of 10 core genes were screened out, including MMP9, AGT, CECR2, FGR, Hist1h3c, APOB, CDH3, MMP15, Hist2h4, and NPPB. The results of the real-time fluorescent quantitative PCR indicated that the mRNA expressions of MMP9, FGR, and APOB in synovial tissues of rats in the model group were higher than those in the normal group (P<0.05), which was consistent with the RNA sequencing results. ConclusionCore genes such as MMP9, FGR, APOB, etc.,probably participate in the occurrence and development of RA by participating in the regulation of immune inflammatory response.