RNA Sequencing and Bioinformatics Analysis of Differentially Expressed Genes in the Peripheral Serum of Ankylosing Spondylitis Patients

Overview

Journal J Orthop Surg Res

Publisher Biomed Central

Specialty Orthopedics

Date 2023 May 30

PMID 37254181

Authors

Yongchen Bie

Xiujun Zheng

Xiaojiong Chen

Xiangyun Liu

Liqin Wang

Yuanliang Sun

Jianqiang Kou

Affiliations

Soon will be listed here.

Abstract

Background: Ankylosing spondylitis (AS) is a chronic progressive autoimmune disease characterized by spinal and sacroiliac arthritis, but its pathogenesis and genetic basis are largely unclear.

Methods: We randomly selected three serum samples each from an AS and a normal control (NC) group for high-throughput sequencing followed by using edgeR to find differentially expressed genes (DEGs). Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes, Reactome pathway analyses, and Gene Set Enrichment Analysis were used to comprehensively analyze the possible functions and pathways involved with these DEGs. Protein-protein interaction (PPI) networks were constructed using the STRING database and Cytoscape. The modules and hub genes of these DEGs were identified using MCODE and CytoHubba plugins. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to validate the expression levels of candidate genes in serum samples from AS patients and healthy controls.

Results: We successfully identified 100 significant DEGs in serum. When we compared them with the NC group, 49 of these genes were upregulated in AS patients and 51 were downregulated. GO function and pathway enrichment analysis indicated that these DEGs were mainly enriched in several signaling pathways associated with endoplasmic reticulum stress, including protein processing in the endoplasmic reticulum, unfolded protein response, and ubiquitin-mediated proteolysis. We also constructed a PPI network and identified the highly connected top 10 hub genes. The expression levels of the candidate hub genes PPARG, MDM2, DNA2, STUB1, UBTF, and SLC25A37 were then validated by RT-qPCR analysis. Finally, receiver operating characteristic curve analysis suggested that PPARG and MDM2 may be the potential biomarkers of AS.

Conclusions: These findings may help to further elucidate the pathogenesis of AS and provide valuable potential gene biomarkers or targets for the diagnosis and treatment of AS.

Citing Articles

attenuates the intestinal permeability, oxidative stress and endoplasmic reticulum stress: transcriptome and microbiome analyses in weaned piglets.

Yuan J, Meng H, Liu Y, Wang L, Zhu Q, Wang Z Front Microbiol. 2024; 15:1362487.

PMID: 38808274 PMC: 11131103. DOI: 10.3389/fmicb.2024.1362487.

Identification and bioinformatics analysis of lncRNAs in serum of patients with ankylosing spondylitis.

Kou J, Bie Y, Liu M, Wang L, Liu X, Sun Y BMC Musculoskelet Disord. 2024; 25(1):291.

PMID: 38622662 PMC: 11017588. DOI: 10.1186/s12891-024-07396-z.

References

Dean L, Macfarlane G, Jones G . Five Potentially Modifiable Factors Predict Poor Quality of Life in Ankylosing Spondylitis: Results from the Scotland Registry for Ankylosing Spondylitis. J Rheumatol. 2017; 45(1):62-69. DOI: 10.3899/jrheum.160411. View

Fattahi F, Kiani J, Khosravi M, Vafaei S, Mohammadi A, Madjd Z . Enrichment of Up-regulated and Down-regulated Gene Clusters Using Gene Ontology, miRNAs and lncRNAs in Colorectal Cancer. Comb Chem High Throughput Screen. 2019; 22(8):534-545. DOI: 10.2174/1386207321666191010114149. View

Fougerat A, Montagner A, Loiseau N, Guillou H, Wahli W . Peroxisome Proliferator-Activated Receptors and Their Novel Ligands as Candidates for the Treatment of Non-Alcoholic Fatty Liver Disease. Cells. 2020; 9(7). PMC: 7408116. DOI: 10.3390/cells9071638. View

Rashid T, Ebringer A . Ankylosing spondylitis is linked to Klebsiella--the evidence. Clin Rheumatol. 2006; 26(6):858-64. DOI: 10.1007/s10067-006-0488-7. View

Zhang C, Chen J, Cai L, Wu J, Wang J, Cao L . DNA induction of MDM2 promotes proliferation of human renal mesangial cells and alters peripheral B cells subsets in pediatric systemic lupus erythematosus. Mol Immunol. 2018; 94:166-175. DOI: 10.1016/j.molimm.2018.01.003. View

Ryno L, Wiseman R, Kelly J . Targeting unfolded protein response signaling pathways to ameliorate protein misfolding diseases. Curr Opin Chem Biol. 2013; 17(3):346-52. PMC: 5859939. DOI: 10.1016/j.cbpa.2013.04.009. View

Shannon P, Markiel A, Ozier O, Baliga N, Wang J, Ramage D . Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13(11):2498-504. PMC: 403769. DOI: 10.1101/gr.1239303. View

Zak D, Penn-Nicholson A, Scriba T, Thompson E, Suliman S, Amon L . A blood RNA signature for tuberculosis disease risk: a prospective cohort study. Lancet. 2016; 387(10035):2312-2322. PMC: 5392204. DOI: 10.1016/S0140-6736(15)01316-1. View

Pimentel-Santos F, Ligeiro D, Matos M, Mourao A, Costa J, Santos H . Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects. Arthritis Res Ther. 2011; 13(2):R57. PMC: 3132052. DOI: 10.1186/ar3309. View

10.

Wang X, Liotta L . Clinical bioinformatics: a new emerging science. J Clin Bioinforma. 2011; 1(1):1. PMC: 3097094. DOI: 10.1186/2043-9113-1-1. View

11.

Layh-Schmitt G, Colbert R . The interleukin-23/interleukin-17 axis in spondyloarthritis. Curr Opin Rheumatol. 2008; 20(4):392-7. PMC: 3494971. DOI: 10.1097/BOR.0b013e328303204b. View

12.

Zeng L, Lindstrom M, Smith J . Ankylosing spondylitis macrophage production of higher levels of interleukin-23 in response to lipopolysaccharide without induction of a significant unfolded protein response. Arthritis Rheum. 2011; 63(12):3807-17. PMC: 3228355. DOI: 10.1002/art.30593. View

13.

Chin C, Chen S, Wu H, Ho C, Ko M, Lin C . cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014; 8 Suppl 4:S11. PMC: 4290687. DOI: 10.1186/1752-0509-8-S4-S11. View

14.

Li Z, Sun Y, He M, Liu J . Differentially-expressed mRNAs, microRNAs and long noncoding RNAs in intervertebral disc degeneration identified by RNA-sequencing. Bioengineered. 2021; 12(1):1026-1039. PMC: 8806213. DOI: 10.1080/21655979.2021.1899533. View

15.

Lin A, Inman R, Streutker C, Zhang Z, Pritzker K, Tsui H . Lipocalin 2 links inflammation and ankylosis in the clinical overlap of inflammatory bowel disease (IBD) and ankylosing spondylitis (AS). Arthritis Res Ther. 2020; 22(1):51. PMC: 7081573. DOI: 10.1186/s13075-020-02149-4. View

16.

Liu Y, Luo S, Zhan Y, Wang J, Zhao R, Li Y . Increased Expression of PPAR-γ Modulates Monocytes Into a M2-Like Phenotype in SLE Patients: An Implicative Protective Mechanism and Potential Therapeutic Strategy of Systemic Lupus Erythematosus. Front Immunol. 2021; 11:579372. PMC: 7873911. DOI: 10.3389/fimmu.2020.579372. View

17.

Ermann J . Pathogenesis of Axial Spondyloarthritis - Sources and Current State of Knowledge. Rheum Dis Clin North Am. 2020; 46(2):193-206. DOI: 10.1016/j.rdc.2020.01.016. View

18.

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

19.

Burton P, Clayton D, Cardon L, Craddock N, Deloukas P, Duncanson A . Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat Genet. 2007; 39(11):1329-37. PMC: 2680141. DOI: 10.1038/ng.2007.17. View

20.

Lu W, Lu Y, Wang C, Chen T . Expression profiling and bioinformatics analysis of exosomal long noncoding RNAs in patients with myasthenia gravis by RNA sequencing. J Clin Lab Anal. 2021; 35(4):e23722. PMC: 8059713. DOI: 10.1002/jcla.23722. View