Journal of Inflammation Research (Nov 2023)
The m6A/m1A/m5C-Related Methylation Modification Patterns and Immune Landscapes in Rheumatoid Arthritis and Osteoarthritis Revealed by Microarray and Single-Cell Transcriptome
Abstract
Haishi Zheng,1,2 Yirixiati Aihaiti,2 Yongsong Cai,2 Qiling Yuan,2 Mingyi Yang,2 Zheng Li,2 Ke Xu,2 Peng Xu1,2 1Department of Orthopedics, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, People’s Republic of China; 2Department of Joint Surgery, HongHui Hospital, Xi’an Jiaotong University, Xi’an, People’s Republic of ChinaCorrespondence: Peng Xu, Department of Orthopedics, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, People’s Republic of China, Tel +86 13772090019, Email [email protected]: The goal of this study was to explore the expression characteristics of RNA modification-related genes, reveal immune landscapes and identify novel potential diagnostic biomarkers in osteoarthritis (OA) and rheumatoid arthritis (RA) patients.Patients and Methods: RNA microarray and single-cell sequencing (scRNA-seq) data were downloaded from gene expression omnibus (GEO) database. Differentially expressed RNA modification-related genes were identified and then functionally annotated. Univariate logistic regression and lasso regression analysis were used to identify primary disease genes for OA and RA. Validation was done using scRNA-seq analysis and immunohistochemistry (IHC) in human knee synovial tissues and a murine destabilization of the medial meniscus (DMM) model. Through WGCNA analysis, genes associated with cell pyroptosis or autophagy in OA and RA were identified, which were then combined with differentially expressed RNA modification-related genes to construct a PPI interaction network. Furthermore, hub genes were selected for ceRNA interaction network analysis, correlation analysis with OA and RA molecular subtypes, as well as correlation analysis with 22 immune cells.Results: Six RNA modification-related genes (ADAMDEC1, IGHM, OGN, TNFRSF11B, SCARA3 and PTN) were identified as potential OA and RA pathogenesis biomarkers. Their expression was validated in human knee synovial tissues and a murine DMM model. Functional enrichment of differentially expressed RNA modification-related genes between RA and OA was analyzed using GO, KEGG, GSEA, and GSVA. Based on WGCNA and PPI analysis, the six hub genes related to pyroptosis and RNA modification (CXCL10, CXCL9, CCR7, CCL5, CXCL1, and CCR2) were identified as central nodes for ceRNA interaction, correlation with OA and RA molecular subtypes, and association with 22 immune cells.Conclusion: Our research revealed the significance of RNA modification-related genes in the development of OA and RA pathogenesis, thereby providing a novel research direction for understanding the mechanisms, diagnosis, and treatment of OA and RA.Keywords: m6A, m5C, m1A, osteoarthritis, rheumatoid arthritis, single-cell sequencing
