切换至 "中华医学电子期刊资源库"

中华关节外科杂志(电子版) ›› 2021, Vol. 15 ›› Issue (03) : 294 -301. doi: 10.3877/cma.j.issn.1674-134X.2021.03.006

基础论著

髋关节滑膜细胞组成及代谢差异生物信息学分析
李志文1, 康焱1, 张紫机1, 盛璞义1, 廖威明1, 毛谷平1,()   
  1. 1. 510080 广州,中山大学附属第一医院关节外科
  • 收稿日期:2021-03-25 出版日期:2021-06-01
  • 通信作者: 毛谷平
  • 基金资助:
    国家自然科学基金(81972051); 广东省医学科研基金(A2020269)

Bioinformatics analysis of differences in synovial cell composition and metabolism of hip joint

Zhiwen Li1, Yan Kang1, Ziji Zhang1, Puyi Sheng1, Weiming Liao1, Guping Mao1,()   

  1. 1. The First Affiliated Hospital of Sun Yat-sen University, Guangzhou 510080, China
  • Received:2021-03-25 Published:2021-06-01
  • Corresponding author: Guping Mao
引用本文:

李志文, 康焱, 张紫机, 盛璞义, 廖威明, 毛谷平. 髋关节滑膜细胞组成及代谢差异生物信息学分析[J]. 中华关节外科杂志(电子版), 2021, 15(03): 294-301.

Zhiwen Li, Yan Kang, Ziji Zhang, Puyi Sheng, Weiming Liao, Guping Mao. Bioinformatics analysis of differences in synovial cell composition and metabolism of hip joint[J]. Chinese Journal of Joint Surgery(Electronic Edition), 2021, 15(03): 294-301.

目的

利用基因表达综合(GEO)数据库中的肥胖及正常体重的髋关节炎滑膜组织的单细胞核糖核酸(RNA)测序结果(GSE152815),分析两者中细胞组成及功能的差异,探讨其代谢重编程及微环境稳态失衡的变化,以期为相关领域的研究者提供参考。

方法

利用GEO数据库中的GSE152815数据,通过质控、无监督聚类、轨迹分析、基因富集分析、基因调控网络分析等方法,采用Wilcoxon秩和检验和Kruskal-Wallis秩和检验等确定基因表达差异之间的统计学意义。

结果

通过质控、无监督聚类,并根据各个亚群的标志分子,可以将滑膜组织分为4群,分别为滑膜成纤维祖细胞(SFPCs),滑膜肥大软骨细胞(SHTCs),滑膜成纤维细胞(SFBs)和滑膜纤维软骨细胞(SFCs)。轨迹分析发现SFPCs可向SHTCs、SFBs及SFCs分化,在分化过程中,脂肪酸代谢相关的基因,载脂蛋白E(APOE)以及ATP结合盒转运体A1(ABCA1)等基因的表达量逐渐上升(Z =25.17、17.89,均为P<0.05),而葡萄糖代谢的相关基因的表达量则逐渐下降(Z =15.32,P<0.05)。单细胞基因调控网络分析发现转录因子SOX4和MITF可以通过联合调控APOE和ABCA1等分子来调控细胞的脂肪酸代谢功能。

结论

肥胖患者及正常体重患者的滑膜的微环境紊乱,可能会导致部分滑膜细胞的代谢紊乱,由原来的葡萄糖代谢重编程为脂肪酸代谢,而这可能与转录因子SRY盒转录因子4(SOX4)和小眼畸形转录因子(MITF)调控的APOE、ABCA1等基因密切相关。明确代谢重编程在这些细胞中的变化,有助于研究不同代谢重编程对细胞功能的影响,为研究骨关节炎等疾病的代谢重编程提供思路以及新途径。

Objective

To use the single-cell RNA sequencing results (GSE152815) of obese and normal-weight hip arthritis synovial tissues in gene expression omnibus(GEO) database to analyze the differences in cell composition and function between the two, and to explore their metabolic reprogramming and the imbalance of microenvironment homeostasis, in order to provide references for researchers in related fields.

Methods

Using the GSE152815 data in the GEO database, through quality control, unsupervised clustering, trajectory analysis, gene enrichment analysis, gene regulatory network analysis and other methods, Wilcoxon rank sum test and Kruskal-Wallis rank sum test were used to determine the statistical significance of differences in gene expression.

Results

Through quality control, unsupervised clustering, and according to the marker molecules of each subgroup, synovial tissue can be divided into four groups, namely synovial fibroblast progenitor cells (SFPCs), synovial hypertrophic chondrocytes (SHTCs), synovial fibroblast cells (SFBs) and synovial fibrochondrocytes (SFCs). Trajectory analysis found that SFPCs can differentiate into SFBs and SFCs. During the differentiation process, the expression of fatty acid metabolism-related genes-apolipoprotein E (APOE) and ATP binding cassette transporter A1 (ABCA1) gradually increase(Z =25.17, 17.89, both P<0.05), while the ability to metabolize glucose gradually declines(Z =15.32, P<0.05). Analysis of single-cell gene regulatory network found that transcription factors [SRY-box transcription factor 4(SOX4) and microphthalmia-associated transcription factor (MITF)] could regulate cell fatty acid metabolism by co-regulating molecules such as APOE and ABCA1.

Conclusions

The results confirmed that the microenvironment disorder of the synovial membrane of obese patients and normal-weight patients will cause the metabolic disorder of some synovial cells, reprogramming from the original glucose metabolism to fatty acid metabolism, which may be closely related to genes such as APOE and ABCA1 regulated by transcription factor SOX4 and MITF. Clarifying the changes of metabolic reprogramming in these cells will help to study the influence of different metabolic reprogramming on cell function, and provide ideas and new ways for studying metabolic reprogramming of osteoarthritis and other diseases.

图1 滑膜组织中4个亚群的UMAP展示图
图2 滑膜组织中4个亚群的标志基因的小提琴展示图
图3 滑膜组织中4个亚群的标志基因的热图展示图
图4 滑膜组织中4个亚群的来源NW(正常体重)和OB(肥胖)占比条形图
图5 Monocle3细胞亚群轨迹分析图
图6 SFBs(滑膜成纤维细胞)和SFCs(滑膜纤维软骨细胞)亚群间的GSEA(基因集合富集分析)分析图
图7 SHTCs(滑膜肥大软骨细胞)和SFCs(滑膜纤维软骨细胞)亚群间的GSEA(基因集合富集分析)分析图
图8 葡萄糖氧化代谢通路与脂肪酸氧化代谢通路相关基因在4个亚群中表达量的热图
图9 OB(肥胖)患者的所有细胞亚群的轨迹分析展示图
图10 基因随细胞分化轨迹的表达变化图
图11 NW(正常体重)患者的所有细胞亚群轨迹分析展示图
图12 基因随分化轨迹的表达变化图
图13 OB(肥胖)患者的细胞分化轨迹1和2上的SFCs(滑膜纤维软骨细胞)细胞亚群的基因表达量差异的火山图展示,显示这两条分化轨迹的SFCs在基因表达方面的差异(P<0.05且基因相对表达量倍数>2)
图14 NW(正常体重)患者的细胞分化轨迹1和2上的SFCs(滑膜纤维软骨细胞)细胞亚群的基因表达量差异的火山图
图15 SOX4(SRY盒转录因子4)和MITF(小眼畸形转录因子)转录因子在滑膜组织4个亚群中的平均表达量热图展示
图16 SOX4(转录因子SRY盒转录因子4)和MITF(小眼畸形转录因子)可联合基因的网络图展示
[1]
Barnett R. Osteoarthritis [J/OL]. Lancet, 2018, 391(1134): 1985. doi: 10.1016/S0140-6736(18)31064-X.
[2]
Rahmati M, Nalesso G, Mobasheri A, et al. Aging and osteoarthritis: Central role of the extracellular matrix[J]. Ageing Res Rev, 2017, 40(11): 20-30.
[3]
文星钊,张志奇.非编码RNA在骨关节炎中的研究进展[J/CD].中华关节外科杂志(电子版),2020,14(2):189-195.
[4]
Chen Z, Ma Y, Li X, et al. The immune cell landscape in different anatomical structures of knee in osteoarthritis: a gene expression-based study[J/OL]. Biomed Res Int, 2020,2020: 9647072. doi: 10.1155/2020/9647072.
[5]
Lopes E, Filiberti A, Husain SA, et al. Immune contributions to osteoarthritis[J]. Curr Osteoporos Rep, 2017, 15(6): 593-600.
[6]
Dicarlo M, Teti G, Cerqueni G, et al. Synovium-derived stromal cell-induced osteoclastogenesis: a potential osteoarthritis trigger[J]. Clin Sci (Lond), 2019, 133(16): 1813-1824.
[7]
孙祥燚,张雷,陈烁,等.微小RNA-564基因下调对滑膜间充质干细胞成软骨分化的影响[J/CD].中华关节外科杂志(电子版),2019,13(1):52-57, 51.
[8]
Al-Khami AA, Rodriguez PC, Ochoa AC. Energy metabolic pathways control the fate and function of myeloid immune cells[J]. J Leukoc Biol, 2017, 102(2): 369-380.
[9]
Kumar V. Inflammation research sails through the sea of immunology to reach immunometabolism[J]. Int Immunopharmacol, 2019, 73(8): 128-145.
[10]
李琼,陈宇婷,接力刚,等.单细胞转录组测序技术在类风湿性关节炎的研究进展[J].分子诊断与治疗杂志,2020,12(6):835-838.
[11]
Nanus DE, Wijesinghe SN, Pearson MJ, et al. Regulation of the inflammatory synovial fibroblast phenotype by metastasis-associated lung adenocarcinoma transcript 1 long noncoding RNA in obese patients with osteoarthritis[J]. Arthritis Rheumatol, 2020, 72(4): 609-619.
[12]
Stuart T, Butler A, Hoffman P, et al. Comprehensive integration of single-cell data[J/OL]. Cell, 2019, 177(7): 1888-1902.e21. doi: 10.1016/j.cell.2019.05.031.
[13]
Korsunsky I, Millard N, Fan J, et al. Fast, sensitive and accurate integration of single-cell data with Harmony[J]. Nat Methods, 2019, 16(12): 1289-1296.
[14]
Trapnell C, Cacchiarelli D, Grimsby J, et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells[J]. Nat Biotechnol, 2014, 32(4): 381-386.
[15]
Qiu X, Hill A, Packer J, et al. Single-cell mRNA quantification and differential analysis with Census[J]. Nat Methods, 2017, 14(3): 309-315.
[16]
Qiu X, Mao Q, Tang Y, et al. Reversed graph embedding resolves complex single-cell trajectories[J]. Nat Methods, 2017, 14(10): 979-982.
[17]
Subramanian A, Tamayo P, Mootha VK, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles [J]. Proc Natl Acad Sci USA, 2005, 102(43): 15545-15550.
[18]
Aibar S, González-Blas CB, Moerman T, et al. SCENIC: single-cell regulatory network inference and clustering[J]. Nat Methods, 2017, 14(11): 1083-1086.
[19]
Ritchie ME, Phipson B, Wu D, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies[J]. Nucleic Acids Res, 2015, 43(7): 47-60.
[20]
Sun H, Wen X, Li H, et al. Single-cell RNA-seq analysis identifies meniscus progenitors and reveals the progression of meniscus degeneration [J]. Ann Rheum Dis, 2020, 79(3): 408-417.
[21]
Enge M, Arda HE, Mignardi M, et al. Single-cell analysis of human pancreas reveals transcriptional signatures of aging and somatic mutation patterns [J/OL]. Cell, 2017, 171(2): 321-330.e14. doi: 10.1016/j.cell.2017.09.004.
[22]
Chen YJ, Chang W, Wu LY, et al. Systematic analysis of transcriptomic profile of chondrocytes in osteoarthritic knee using next-generation sequencing and bioinformatics [J/OL]. J Clin Med, 2018, 7(12): 535. doi: 10.3390/jcm7120535.
[23]
Schelker M, Feau S, Du J, et al. Estimation of immune cell content in tumour tissue using single-cell RNA-seq data[J]. Nat Commun, 2017, 8(1): 2032-2046.
[24]
Ji Q, Zheng Y, Zhang G, et al. Single-cell RNA-seq analysis reveals the progression of human osteoarthritis[J]. Ann Rheum Dis, 2019, 78(1): 100-110.
[25]
Gunja NJ, Athanasiou KA. Passage and reversal effects on gene expression of bovine meniscal fibrochondrocytes[J/OL]. Arthritis Res Ther, 2007, 9(5): R93. doi: 10.1186/ar2293.
[26]
Satizabal CL, Samieri C, Davis-Plourde KL, et al. APOE and the association of fatty acids with the risk of stroke, coronary heart disease, and mortality[J]. Stroke, 2018, 49(12): 2822-2829.
[27]
Qian H, Zhao X, Cao P, et al. Structure of the human lipid exporter ABCA1[J/OL]. Cell, 2017, 169(7): 1228-1239.e10. doi: 10.1016/j.cell.2017.05.020.
[28]
Jiao Y, Zhao J, Zhang Z, et al. SRY-Box containing gene 4 promotes liver steatosis by upregulation of SREBP-1c[J]. Diabetes, 2018, 67(11): 2227-2238.
[29]
Vivas-García Y, Falletta P, Liebing J, et al. Lineage-restricted regulation of SCD and fatty acid saturation by MITF controls melanoma phenotypic plasticity[J/OL]. Mol Cell, 2020, 77(1): 120-137.e9. doi: 10.1016/j.molcel.2019.10.014.
[1] 刘婷婷, 林妍冰, 汪珊, 陈幕荣, 唐子鉴, 代东伶, 夏焙. 超声衰减参数成像评价儿童代谢相关脂肪性肝病的价值[J]. 中华医学超声杂志(电子版), 2023, 20(08): 787-794.
[2] 周钰菡, 肖欢, 唐毅, 杨春江, 周娟, 朱丽容, 徐娟, 牟芳婷. 超声对儿童髋关节暂时性滑膜炎的诊断价值[J]. 中华医学超声杂志(电子版), 2023, 20(08): 795-800.
[3] 闫文, 谢兴文, 顾玉彪, 雷宁波, 马成, 于文霞, 高亚雄, 张磊. 微小RNA与全膝关节置换术后深静脉血栓的研究进展[J]. 中华关节外科杂志(电子版), 2023, 17(06): 842-846.
[4] 樊绪国, 赵永刚, 杨砚伟. 腓骨在膝骨关节炎作用的研究观点[J]. 中华关节外科杂志(电子版), 2023, 17(06): 855-859.
[5] 李善武, 叶永杰, 王兵, 王子呓, 银毅, 孙官军, 张大刚. 胫骨高位截骨与单髁置换的早期疗效比较[J]. 中华关节外科杂志(电子版), 2023, 17(06): 882-888.
[6] 马伟强, 马斌林, 吴中语, 张莹. microRNA在三阴性乳腺癌进展中发挥的作用[J]. 中华普外科手术学杂志(电子版), 2024, 18(01): 111-114.
[7] 李婷, 张琳. 血清脂肪酸代谢物及维生素D水平与结直肠癌发生的关系研究[J]. 中华普外科手术学杂志(电子版), 2023, 17(06): 661-665.
[8] 顾志波, 郝林, 陆明, 陈建刚. 光动力纳米载体联合si-P3H4治疗膀胱癌的初步探索[J]. 中华腔镜泌尿外科杂志(电子版), 2023, 17(06): 633-641.
[9] 何吉鑫, 杨燕妮, 王继伟, 李建国, 谢铭. 肠道菌群及肠道代谢产物参与慢性便秘发生机制的研究进展[J]. 中华结直肠疾病电子杂志, 2023, 12(06): 495-499.
[10] 王旭, 师绍敏, 毛燕, 季上, 刘亚玲. 肝酶代谢与骨关节炎相关性的研究进展[J]. 中华老年骨科与康复电子杂志, 2023, 09(06): 379-384.
[11] 黄岩, 刘晓巍, 杨春玲, 兰烨. 急性胰腺炎合并糖尿病患者的临床特征及血糖代谢与病情严重度的相关性[J]. 中华消化病与影像杂志(电子版), 2023, 13(06): 439-442.
[12] 袁媛, 赵良平, 刘智慧, 张丽萍, 谭丽梅, 閤梦琴. 子宫内膜癌组织中miR-25-3p、PTEN的表达及与病理参数的关系[J]. 中华临床医师杂志(电子版), 2023, 17(9): 1016-1020.
[13] 王磊, 李梦, 孙文利, 刘瑞, 王红春, 卢光泽, 赵颖, 郭进艳, 刘红星. 液相色谱质谱法对急性白血病患者血浆代谢组学的特征分析[J]. 中华临床医师杂志(电子版), 2023, 17(08): 850-857.
[14] 张大涯, 陈世锔, 陈润祥, 张晓冬, 李达, 白飞虎. 肠道微生物群对代谢相关脂肪性肝病发展的影响[J]. 中华临床医师杂志(电子版), 2023, 17(07): 828-833.
[15] 周加军, 余永武, 周涵, 刘勇, 张凌. 甲状旁腺切除对继发性甲状旁腺功能亢进患者骨密度及骨代谢的影响[J]. 中华临床医师杂志(电子版), 2023, 17(06): 706-710.
阅读次数
全文


摘要