Identification of Key Genes and Long Non‑coding RNA Expression Profiles in Osteoporosis with Rheumatoid Arthritis Based on Bioinformatics Analysis

Overview

Journal BMC Musculoskelet Disord

Publisher Biomed Central

Specialties Orthopedics
Physiology

Date 2024 Aug 8

PMID 39118036

Authors

Jin-Yu An

Xing-Na Ma

Hui-Long Wen

Hui-Dong Hu

Affiliations

Soon will be listed here.

Abstract

Background: Although rheumatoid arthritis (RA) is a chronic systemic tissue disease often accompanied by osteoporosis (OP), the molecular mechanisms underlying this association remain unclear. This study aimed to elucidate the pathogenesis of RA and OP by identifying differentially expressed mRNAs (DEmRNAs) and long non-coding RNAs (lncRNAs) using a bioinformatics approach.

Methods: Expression profiles of individuals diagnosed with OP and RA were retrieved from the Gene Expression Omnibus database. Differential expression analysis was conducted. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) pathway enrichment analyses were performed to gain insights into the functional categories and molecular/biochemical pathways associated with DEmRNAs. We identified the intersection of common DEmRNAs and lncRNAs and constructed a protein-protein interaction (PPI) network. Correlation analysis between the common DEmRNAs and lncRNAs facilitated the construction of a coding-non-coding network. Lastly, serum peripheral blood mononuclear cells (PBMCs) from patients with RA and OP, as well as healthy controls, were obtained for TRAP staining and qRT-PCR to validate the findings obtained from the online dataset assessments.

Results: A total of 28 DEmRNAs and 2 DElncRNAs were identified in individuals with both RA and OP. Chromosomal distribution analysis of the consensus DEmRNAs revealed that chromosome 1 had the highest number of differential expression genes. GO and KEGG analyses indicated that these DEmRNAs were primarily associated with " platelets (PLTs) degranulation", "platelet alpha granules", "platelet activation", "tight junctions" and "leukocyte transendothelial migration", with many genes functionally related to PLTs. In the PPI network, MT-ATP6 and PTGS1 emerged as potential hub genes, with MT-ATP6 originating from mitochondrial DNA. Co-expression analysis identified two key lncRNA-mRNA pairs: RP11 - 815J21.2 with MT - ATP6 and RP11 - 815J21.2 with PTGS1. Experimental validation confirmed significant differential expression of RP11-815J21.2, MT-ATP6 and PTGS1 between the healthy controls and the RA + OP groups. Notably, knockdown of RP11-815J21.2 attenuated TNF + IL-6-induced osteoclastogenesis.

Conclusions: This study successfully identified shared dysregulated genes and potential therapeutic targets in individuals with RA and OP, highlighting their molecular similarities. These findings provide new insights into the pathogenesis of RA and OP and suggest potential avenues for further research and targeted therapies.

References

Eroglu S, Karatas G . Platelet/lymphocyte ratio is an independent predictor for osteoporosis. Saudi Med J. 2019; 40(4):360-366. PMC: 6506656. DOI: 10.15537/smj.2019.4.24009. View

Miyazaki T, Iwasawa M, Nakashima T, Mori S, Shigemoto K, Nakamura H . Intracellular and extracellular ATP coordinately regulate the inverse correlation between osteoclast survival and bone resorption. J Biol Chem. 2012; 287(45):37808-23. PMC: 3488055. DOI: 10.1074/jbc.M112.385369. View

Zhu D, Chen X, Hu W, Dong S, Lu B, Rong Y . Multiple Functional Variants at 13q14 Risk Locus for Osteoporosis Regulate RANKL Expression Through Long-Range Super-Enhancer. J Bone Miner Res. 2018; 33(7):1335-1346. DOI: 10.1002/jbmr.3419. View

Zhang H, Meltzer P, Davis S . RCircos: an R package for Circos 2D track plots. BMC Bioinformatics. 2013; 14:244. PMC: 3765848. DOI: 10.1186/1471-2105-14-244. View

Moshayedi S, Tasorian B, Almasi-Hashiani A . The prevalence of osteoporosis in rheumatoid arthritis patient: a systematic review and meta-analysis. Sci Rep. 2022; 12(1):15844. PMC: 9508181. DOI: 10.1038/s41598-022-20016-x. View

Ritchie M, Phipson B, Wu D, Hu Y, Law C, Shi W . limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43(7):e47. PMC: 4402510. DOI: 10.1093/nar/gkv007. View

Li X, Zhang J, Meng S, Li Y, Wang R . Mean platelet volume is negatively associated with bone mineral density in postmenopausal women. J Bone Miner Metab. 2012; 30(6):660-5. DOI: 10.1007/s00774-012-0362-4. View

Raterman H, Bultink I, Lems W . Osteoporosis in patients with rheumatoid arthritis: an update in epidemiology, pathogenesis, and fracture prevention. Expert Opin Pharmacother. 2020; 21(14):1725-1737. DOI: 10.1080/14656566.2020.1787381. View

Lopez Herraez D, Martinez-Bueno M, Riba L, Garcia de la Torre I, Sacnun M, Goni M . Rheumatoid arthritis in Latin Americans enriched for Amerindian ancestry is associated with loci in chromosomes 1, 12, and 13, and the HLA class II region. Arthritis Rheum. 2013; 65(6):1457-67. DOI: 10.1002/art.37923. View

10.

Collins L, Hajizadeh S, Holme E, Jonsson I, Tarkowski A . Endogenously oxidized mitochondrial DNA induces in vivo and in vitro inflammatory responses. J Leukoc Biol. 2004; 75(6):995-1000. DOI: 10.1189/jlb.0703328. View

11.

Chen H, Boutros P . VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011; 12:35. PMC: 3041657. DOI: 10.1186/1471-2105-12-35. View

12.

Mattick J, Amaral P, Carninci P, Carpenter S, Chang H, Chen L . Long non-coding RNAs: definitions, functions, challenges and recommendations. Nat Rev Mol Cell Biol. 2023; 24(6):430-447. PMC: 10213152. DOI: 10.1038/s41580-022-00566-8. View

13.

Papaioannou A, Kennedy C, Ioannidis G, Sawka A, Hopman W, Pickard L . The impact of incident fractures on health-related quality of life: 5 years of data from the Canadian Multicentre Osteoporosis Study. Osteoporos Int. 2008; 20(5):703-14. PMC: 5101052. DOI: 10.1007/s00198-008-0743-7. View

14.

Jang H, Lee H, Noh E, Kim J, You Y, Lee G . Aqueous extract of Chrysanthemum morifolium Ramat. inhibits RANKL-induced osteoclast differentiation by suppressing the c-fos/NFATc1 pathway. Arch Oral Biol. 2021; 122:105029. DOI: 10.1016/j.archoralbio.2020.105029. View

15.

Rodan G, Martin T . Therapeutic approaches to bone diseases. Science. 2000; 289(5484):1508-14. DOI: 10.1126/science.289.5484.1508. View

16.

Boilard E, Nigrovic P, Larabee K, Watts G, Coblyn J, Weinblatt M . Platelets amplify inflammation in arthritis via collagen-dependent microparticle production. Science. 2010; 327(5965):580-3. PMC: 2927861. DOI: 10.1126/science.1181928. View

17.

Li M, Ma K, Feng Z, Wang J, Zhou X, Zhou L . Differential long non-coding RNA expression profiles in the peripheral blood and CD4 T cells of patients with active rheumatoid arthritis. Exp Ther Med. 2020; 20(1):461-471. PMC: 7271723. DOI: 10.3892/etm.2020.8681. View

18.

Doncheva N, Morris J, Gorodkin J, Jensen L . Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J Proteome Res. 2018; 18(2):623-632. PMC: 6800166. DOI: 10.1021/acs.jproteome.8b00702. View

19.

Yu G, Wang L, Han Y, He Q . clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16(5):284-7. PMC: 3339379. DOI: 10.1089/omi.2011.0118. View

20.

Kim S, Schneeweiss S, Liu J, Daniel G, Chang C, Garneau K . Risk of osteoporotic fracture in a large population-based cohort of patients with rheumatoid arthritis. Arthritis Res Ther. 2010; 12(4):R154. PMC: 2945054. DOI: 10.1186/ar3107. View