- Open Access
Dysregulation of ribosome-related genes in ankylosing spondylitis: a systems biology approach and experimental method
BMC Musculoskeletal Disorders volume 22, Article number: 789 (2021)
Ankylosing spondylitis (AS) is an autoimmune rheumatic disease. Few candidate gene associations have been reported for AS and the current understanding of its pathogenesis remains still poor. Thus, the exact mechanism of AS is needed to urgently be disclosed. The purpose of this study was to identify candidate genes involving in AS disease.
Methods and results
GSE25101 publicly available microarray and GSE117769 RNA-seq datasets of AS patients were obtained for bioinformatics analyses. Gene set enrichment analysis showed that in the microarray dataset, the ribosome pathway was significantly up-regulated in AS compared with controls. Furthermore, some ribosomal components demonstrated overexpression in patients in the RNA-seq dataset. To confirm the findings, 20 AS patients and 20 matching controls were selected from the Rheumatology Research Center clinic, Shariati Hospital. PBMCs were separated from whole blood and RNA contents were extracted. Following the results of datasets analysis, the expression level of rRNA5.8S pseudogene, rRNA18S pseudogene, RPL23, RPL7, and RPL17 genes were measured through real-time PCR. Our findings showed dysregulation of rRNA5.8S and rRNA18S pseudogenes, and also the RPL17 gene in patients.
Considering that genes involved in ribosome biogenesis contributed to some AS-associated biological processes as well as diseases that have comorbidities with AS, our results might advance our understanding of the pathological mechanisms of ankylosing spondylitis.
Ankylosing spondylitis (AS) is the major subtype of spondyloarthropathies, which is one of chronic inflammatory arthritis [1, 2]. The characteristic symptoms of AS particularly are back pain, spinal stiffness, and loss of spinal mobility [3, 4]. Human leukocyte antigen (HLA) B27 is one of the most convincing genetic biomarkers associated with AS. However genetic association studies reveal that HLA-B27 only contributes to approximately 20% of AS heritability and 30% of the overall risk for AS [5, 6]. Previous studies suggested that AS is a multifactorial disease and the susceptibility to this disorder may be own to genetic and environmental factors . Nevertheless, few genes associated with the disease were identified and the actual cause of AS has remained unclear . Therefore, there is a crucial necessity to investigate the molecular mechanisms of AS in order to find more relevant biomarkers (s). Recent studies investigated other genes and inflammatory biomarkers related to the AS pathogenesis. Environmental factors through altering epigenetics mechanisms alongside with genetics polymorphisms control different gene expressions. Hence, altered gene expression profiles that differentiate disease from the healthy condition can be used as a basis of immuno-pathological mechanisms involved in the pathogenesis of AS [9, 10]. High-throughput genomics technologies such as microarrays and RNA sequencing allow analyzing gene expression profiles of thousand genes and provide deep insight into the interaction and network of the genes during complex diseases .
Eukaryotic ribosomes are the organelles consist of four ribosomal RNAs (rRNAs) and around 80 ribosomal proteins (RPs), whose function is to organize spatiotemporal control of gene expression. The RPs organize in small (40S) and large (60S) subunits and mostly are called with ‘S’ and ‘L’ letters, respectively [12, 13]. Ribonucleoprotein (RNP) components such as ribosomes are often targeted in various autoimmune disorders and recognized by B cells derived auto-antibodies [14, 15]. It was also demonstrated that perturbed expression of the RPs and rRNAs occurs in numerous human disorders, notably, in many cancers, ribosomopathies, and autoimmune diseases [16, 17].
In the current study, by expression profiling of AS-transcriptomes datasets and experimental analysis, we confirmed the dysregulation of some ribosomal genes in the disease. Ribosomes have traditionally been viewed as invariant , and their relationship to AS is not yet known, therefore, the finding of different expressions of ribosomal components in AS patients compared to the control group was interesting.
Methods and donors
Gene expression analysis
Microarray data analysis
The expression profile of GSE25101 which included 16 AS and 16 healthy samples was obtained from Gene Expression Omnibus (GEO) database. Characteristics of subjects involved in the microarray dataset are represented in Supplementary file 1 . Following quantile normalization of microarray dataset, the limma package was used to determine differentially expressed genes between the AS and control groups with empirical Bayes analysis [19, 20]. Gene set enrichment analysis was conducted on normalized and log2 transformed data with the Gene Set Enrichment Analysis (GSEA) tool using curated genes from KEGG dataset (category C2) of MSigDB database .
RNA sequencing data analysis
For ribosomal RNA expression analysis, we used a publicly available dataset consisting of 120 samples with four phenotypes, namely; Rheumatoid arthritis, Ankylosing spondylitis, Psoriatic arthritis, and Healthy (GEO Accession: GSE117769 ). We analyzed all 8 ankylosing spondylitis samples in the dataset, along with an equal number of sex and age-matched control samples (Samples characteristics of RNA-seq dataset are shown in Table 1 and detailed information of patient samples are represented in Supplementary file 2). The RAW data of all selected samples were downloaded from EMBL-EBI (The European Bioinformatics Institute, PRJNA483133 accession number). The obtained RNA-seq libraries were processed by Trimmomatic version 0.38 to removed adapters, low-quality bases and read smaller than 50 bp (ILLUMINACLIP: LEADING:3 TRAILING:3 TruSeq3-SE:2:30:10 SLIDINGWINDOW:4:15 MINLEN:50), and the quality of all libraries was checked with FastQC, before and after trimming . Then, the Salmon index file was built based on Ensembl (version 95) reference transcriptome annotation using k-mers of length 31, and the trimmed reads quantified using quasi-mappings of Salmon tool [24, 25]. We used tximport (version 1.10.1) to convert transcript-level to gene-level count estimates . All differential expression analyses were performed with DESeq2 1.22.2 based on the negative binomial distribution, considering only the transcripts which had a TPM (transcripts per million) ≥1 in all samples .
In this study, 20 AS patients were selected from the outpatient rheumatology Research Center (RRC) clinic, Shariati Hospital, Tehran, Iran. All diagnoses of AS by experienced Rheumatologists were made based on the modified New York criteria . The patient group includes 15 men and 5 women with an average age of 35 years. The average duration since diagnosis is 8.7 years (more demographic characteristics are represented in Table 2). The patients were assessed for functional capacities and disease severity by a protocol based on the Assessment of Spondyloarthritis International Society (ASAS) core set ; Disease activity was evaluated by Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) and only active AS patients with BASDAI ≥4 were selected ; damage or deformity of the spine was assessed by Bath Ankylosing Spondylitis Metrology Index (BASMI) ; Disease function was measured by Bath Ankylosing Spondylitis Functional Index (BASFI) . No biologic-treated patients were included. For the control group, 20 healthy persons who matched in age and gender with the patients were selected and they should have no clinical evidence of any type of autoimmune disorders in themselves and their family. Signed informed consent was obtained from all the participants. The study procedure is confirmed by the Ethics Committee of Tehran University of Medical Sciences.
Sample preparation and quantitative real time PCR
Five-milliliter blood was collected from all subjects in test tubes containing ethylenediamine tetra-acetic acid (EDTA). Samples were transferred to the laboratory on ice and possessed in about 2–3 h from the collection. PBMCs were separated from the whole blood by utilizing Ficoll-Hypaque gradient (Innotrain, Germany). RNA extraction was performed using the High Pure RNA Isolation Kit (Progen lab, Germany). Before cDNA synthesis, the concentration of extracted RNA was quantified by spectrophotometry (NanoDrop, Thermo Fisher Scientific, USA). According to Roche cDNA Synthesis Kit’s protocols, isolated RNA was converted to cDNA. Gene expression was measured by using the SYBR Green method for quantitative Real-time polymerase chain reaction (PCR) (Applied Biosystems, Foster City, USA). The specific primer sequences were shown in Table 3. Comparative Ct method was used for the relative gene expression analysis and Hypoxanthine phosphoribosyl transferase 1 (HPRT1) was used as a housekeeping gene for normalization. The expression data did not pass the normality test, so the Mann–Whitney test was used for mRNA expression analysis. The result with p-value of lower than 0.05 was considered significant.
Bioinformatics analysis outcomes
The analysis of the microarray dataset identified that ribosome pathway was significantly enriched by utilizing the GSEA method (normalized enrichment score (NES) = 2.18, false discovery rate (FDR p-value) = 0.000, and p-value = 0.000) (Fig. 1b). Figure 1a is a volcano plot of all genes (represented by dots) involving in the ribosomal pathway. The interest genes which is statistically significant as evidenced by the p-values and the fold change threshold in the microarray dataset are denoted by hexagon dots. The plot demonstrates higher expression profile changes of the ribosome pathway’s genes in the AS group. Three genes with the highest differentially expression included RPL23 (Log2Fc = 1.04, p-value = 0.001), RPL17 (Log2Fc = 0.9, p-value = 0.003), and RPL7 (Log2Fc = 0.8, p-value = 0.003) were selected for the experimental section.
RNA-seq analysis with DESeq2 confirmed that all of the significant genes (more than a thousand genes with p-value < 0.01) in the dataset can distinct very clearly between two groups (Fig. 2a; samples dendrogram and trait heatmap). In line with the notable expression trend of ribosomal genes in the microarray dataset, it found that ENST00000515896.1 which is a pseudogene of rRNA5.8S is one of the most significant genes in the RNA-seq dataset with Log2FC = 3.67 and p-value = 0.0001. Interestingly the gene has the highest correlation with another gene called AC010970.1(Log2Fc = 2.57, p-value = 0.001) which is an 18S ribosomal RNA pseudogene (Fig. 2b).
Evaluation of mRNA expressions of selected genes
Using RT-PCR assays for the certain genes, we confirmed that the expression level of rRNA5.8S pseudogene (fold change = 4.36, p-value = 0.002; Fig. 3a), rRNA18S pseudogene (fold change = 5.13, p-value = 0.000; Fig. 3b) and RPL17 (fold change = 4.56, p-value = 0.000; Fig. 3c) in AS patients were significantly higher than controls, and will be discussed in the next section. mRNA expression levels of RPL7, and RPL23 were not different between two groups (Fig. 3d and e).
The eukaryote ribosome is macromolecular complex, consisted of 2 subunits, the small 40S ribosomal subunit contains one 18S rRNA and 33 ribosomal proteins (RPS), whereas the large 60S subunit including the 28S, 5.8S, and 5S rRNAs together with 47 ribosomal proteins (RPL) . A variety of extra-ribosomal functions for different ribosome’s components were recently confirmed, including the regulation of immune signaling, cellular development, cell cycle progression, and axial skeleton formation. Therefore, dysregulated expression of different ribosomal components elicits a wide spectrum of phenotypes, from developmental defects to diseases [12, 34]. Indeed, an entire class of diseases originates from defects in ribosome biogenesis called ribosomopathies, like Shwachman Diamond syndrome, Diamond-Blackfan anemia, and Dyskeratosis congenita. Despite the unique nature of their biochemical process, the clinical manifestations are extremely variable [35, 36].
Moreover, cancer cells are constantly associated with an increase in protein synthesis and dysregulation of ribosome biogenesis [13, 37]. RP altered expression pattern has been found in various human cancers, including the pancreas, breast, brain, bladder, and many others which could serve as prognostic or diagnostic markers [12, 16, 38, 39]. The increased expression of ribosomal protein mRNAs has been reported in human colon cancer and hepatocellular carcinoma . Noteworthy, some ribosomal proteins like S13 and L23, provoke multidrug resistance in gastric cancer cells . Furthermore, previous studies have linked increased rRNA expression with cancer development in cervical and prostate cancer [42, 43]. Overexpression of the pre-45S rRNA encourages G1/S cell-cycle transition in colorectal cancer and it is related to poor prognosis .
Ribosome-related networks have also been identified to be associated with different autoimmune rheumatic disorders. In a study by Luan et al. it was shown that ribosome-related pathway is associated with susceptibility to systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and multiple sclerosis (MS), suggesting that these diseases are associated to ribosomal genes . Many of the auto-antibodies and auto-reactive T cells in different autoimmune rheumatic diseases recognized ribosomal proteins as an auto-antigens. It was reported that despite the similar frequency of anti-ribosomal antibodies in RA and SLE, they have different specific targets. Po, L7, L5, Sb, S19, S13, and L2 proteins in SLE, and L35a and L37a in RA identified by immunoblotting after two-dimensional gel electrophoresis with anti-ribosome antibodies [46, 47]. Some studies recognized differential ribosomal gene expression profiles in AS patients. In a recent study, Gracey et al. performed sex-dependent functional network analysis in AS patients and healthy controls and found translation/ribosome–related pathways dysregulated especially in male AS patients . In a meta-analysis of differentially expressed genes and associated biological pathways in AS, was found that S11, L27, and L40 mitochondrial ribosomal proteins are among the most significant downregulated genes in AS patients .
In this study, RPL17, RPL23A, and RPL7 were the three highly differentially expressed genes in our GSE analysis of the microarray dataset. Among the genes, RPL17 was also significantly up-regulated in PBMCs from AS patients in the laboratory. Also in a previous bioinformatics analysis on GSE 25101 by Zhao et al. to predict the related genes to AS, RPL17 was among the ribosomal protein which significantly enriched in the selected module of the up-regulated network . RPL17 gene encodes a protein in humans called 60S ribosomal protein L17. In 2012, an extra-ribosomal function in vascular smooth muscle cells (VSMC) growth was reported for RPL17. Its expression was oppositely correlated with VSMC growth by increased cells in G0/G1 and decreased cells in the S phase. This study indicates that RPL17 could play a part in angiogenesis suppression, however, the exact mechanism remains unclear . It was also involved in Diamond-Blackfan anemia and was reported to have higher expression in tumors across one cancer cohort .
RPL23A is one of the components of the large 60S ribosomal subunit and has been expressed in various organs at a high level. RPL23A can stimulate arthritogenic T cells differentiation to effector T lymphocytes and the production of IL-2 and inflammatory cytokines from mice. RPL23A has been expressed in fibroblast-like synoviocytes (FLS) in healthy and arthritis synovial tissues and is recognized by autoantibodies from patients with auto-inflammatory RA disease. The proportion of anti-RPL23A antibody-positive individuals is significantly higher in the RA patients group compared to healthy controls. Besides, a subset of T cells from RA patients produce interferon-γ upon RPL23A exposure [53, 54]. RPL7 is a Protein Coding gene which is associated with systemic autoimmune diseases, such as SLE and other connective tissue diseases . In contrast with our bioinformatics analysis, we did not find significant differences in RPL23A and RPL7 expression in PBMCs from AS patients in the experimental section, and further investigations are needed for determining their role in the pathogenesis of AS.
In the current study, we also discovered overexpression of rRNA5.8S and rRNA18S pseudogenes in the bioinformatics section and also in the lab in PBMCs of AS patients. Due to the high similarity of pseudogenes to their genes, there is the possibility of binding of the pseudogenes designed primers to the 5.8S rRNA and 18S rRNA genes. Both of the ribosomal RNA is part of the 45S precursor which is transcribed by RNA polymerase I and also contains 28S rRNA . In the eukaryotic ribosome, 5.8S rRNA is a non-coding component of the large subunit (60S) and previous studies showed its function in cell growth and protein synthesis . Additionally, the covalent binding of the 5.8S rRNA to p53 was previously revealed and it was suggested that p53 and the RNA could be two integral parts of a single functional entity . 18S rRNA is the structural RNA in small subunit (40S) of the eukaryotic cytoplasmic ribosomes . Bowen-Conradi syndrome is a disease caused by a single mutation in EMG1. The gene is required for the assembly of 40S ribosomal subunit and is involved in the modification of the 18S rRNA . In that patients, skeletal dysmorphology is observed and serious prenatal and postnatal growth defect usually leads to death by 1 year of age [61, 62].
Pseudogenes are non-coding alleles of normal genes that have long been considered as nonfunctional inactive elements, however, recent evidence demonstrated that many of them are transcribed into RNA, functionally active, and can regulate gene expression as long noncoding RNAs (lncRNAs) . Pseudogenes have tissue-specific expression pattern and altered expression of them are involved in various diseases and physiological condition . They play a crucial role in the immune response and are associated with signaling pathways involved in inflammation. The transcription of diverse pseudogenes is induced by pro-inflammatory cytokines such as tumor necrosis factor- α (TNF-α) . There are not any previous reports, to our knowledge, investigating the expression of ribosomal RNA in AS, while some evidence indicates the ribosomal pseudogenes function in immune responses. In this line, it is demonstrated that 5S ribosomal RNA pseudogene 141 plays an essential role in antiviral innate immunity and induction of cytokine response by binding to retinoic acid-inducible gene-I (RIG-1) receptors .
The relation between epigenetic changes, such as DNA methylation and autoimmunity has been well established in the literature . Moreover, global DNA hypo-methylation was reported in rheumatic diseases such as RA and lupus-like disease compared to the healthy group . Furthermore, in another study, the association between AS and DNA methyltransferase 3A (DNMT3A), DNMT3B and DNMT3L, reinforce the notion that methylation variation is involved in AS pathogenesis . In this regard, our finding of the dysregulation of ribosomal genes might be a confirmation for methylation changes that influence the expression of affected genes and consequently protein translation; although, further studies are needed in order to assess this hypothesis.
Strong evidence disclosed the integral role of ribosomal RNA and protein in various physiological and pathological processes associated with AS like innate immune response that regulate inflammatory cytokine induction (NF-kB, interferon-g (IFN-g)), inflammatory signaling, chondrocyte growth, and skeletal development [34, 70,71,72]. There are also observations regarding the role of the ribosomal components in other autoimmune diseases like MS, RA, and SLE and also diseases with high incidence in AS patients like cancers and atherosclerosis [12, 45, 51, 73, 74]. However little is known about the link between ribosome biogenesis with AS and further studies will be required to determine the relevance of ribosomal genes in the disease.
This study focuses on selected ribosomal genes of the microarray and RNA-seq datasets and evaluated their expression in community samples of Iranian AS patients and healthy controls. However, some limitations should be noted. First, the datasets were different in terms of therapeutic criteria, whereas medical treatments might affect the RNA expression . Second, ankylosing spondylitis mainly has comorbidity with other diseases like bowel inflammation, uveitis, psoriasis, and heart disease . These diseases might also affect the biological pathways and gene expression. In this study, unfortunately, due to a lack of information, the possible role of these comorbid diseases in gene expression has not been considered. Third, the number of samples in the experimental section was limited and the results of the study may be affected by the experimental conditions. Therefore, it is suggested to evaluate the RPL23 and RPL7 genes expression once again with more samples and by different experimental methods. These two genes showed significant dysregulation in patients compared to controls in the microarray dataset, while they did not confirm in laboratory analysis.
In summary, our results identify the altered expression of some ribosomal components in bioinformatics and experimental analysis in AS patients. However, it is not clear whether the alternations in gene expression are the cause or effect of the disease. Knowing the significant contribution of ribosomal gene expression variation in diseases and biological processes associated with AS, no wonder that ribosome dysregulation could in part account for AS pathogenesis.
Availability of data and materials
All data generated or analyzed during this study are available upon request.
Assessment of Spondyloarthritis International Society
Bath Ankylosing Spondylitis Disease Activity Index
Bath Ankylosing Spondylitis Functional Index
Bath Ankylosing Spondylitis Metrology Index
Ethylendiamine tetraacetic acid
Gene Expression Omnibus
Gene Set Enrichment Analysis
Human leukocyte antigen
Hypoxanthine phosphoribosyl transferase 1
Long noncoding RNAs
Molecular Signatures Database
Normalized enrichment score
Peripheral blood mononuclear cells
Polymerase chain reaction
Retinoic acid-inducible gene-I
Ribosomal proteins of large subunit
Ribosomal proteins of small subunit
Rheumatology research center
Systemic lupus erythematosus
Tumor necrosis factor alpha
Transcript per million
Vascular smooth muscle cells
Dibra D, Xia X, Gagea M, Lozano G, Li S. A spontaneous model of spondyloarthropathies that develops bone loss and pathological bone formation: a process regulated by IL27RA−/−and mutant-p53. PLoS One. 2018;13.
Zhu W, He X, Cheng K, Zhang L, Chen D, Wang X, et al. Ankylosing spondylitis: etiology, pathogenesis, and treatments. Bone Res. 2019;7:1–16.
Cao H, Huang P, Cui G, Zhang X, Mao K, Lu N, et al. Integrated analyses of gene profile digs out suitable biomarkers for ankylosing spondylitis using three different methods. Int J Hum Genet. 2018;18:87–94.
Masoumi M, Mahmoudi M, Malekshahi Z, Karami J, Vojdanian M, Jamshidi A. Determination of ETS1 gene single nucleotide polymorphism in Iranian patients with ankylosing spondylitis. Rheumatol Res. 2017;2:133–8.
Brophy S, Hickey S, Menon A, Taylor G, Bradbury L, Hamersma J, et al. Concordance of disease severity among family members with ankylosing spondylitis? J Rheumatol. 2004;31:1775–8.
Brown MA, Kennedy LG, MacGregor AJ, Darke C, Duncan E, Shatford JL, et al. Susceptibility to ankylosing spondylitis in twins: the role of genes, HLA, and the environment. Arthritis Rheum. 1997;40:1823–8.
Reveille JD. The genetic basis of spondyloarthritis. Ann Rheum Dis. 2011;70:i44–50.
Wang Z-H, Xiang D, Dong J-J, He S-X, Guo L-M, Lv J, et al. Identification of dysregulated pathways associated with ankylosing spondylitis using pathway interaction network. Iran Red Crescent Med J. 2017;19.
Li Z, Haynes K, Pennisi DJ, Anderson LK, Song X, Thomas GP, et al. Epigenetic and gene expression analysis of ankylosing spondylitis-associated loci implicate immune cells and the gut in the disease pathogenesis. Genes Immun. 2017;18:135–43.
Hammaker D, Firestein GS. Epigenetics of inflammatory arthritis. Curr Opin Rheumatol. 2018;30:188–96.
Rodriguez-Esteban R, Jiang X. Differential gene expression in disease: a comparison between high-throughput studies and the literature. BMC Med Genet. 2017;10:59.
Guimaraes JC, Zavolan M. Patterns of ribosomal protein expression specify normal and malignant human cells. Genome Biol. 2016;17:236.
Bastide A, David A. The ribosome,(slow) beating heart of cancer (stem) cell. Oncogenesis. 2018;7:1–13.
Calich AL, Viana VS, Cancado E, Tustumi F, Terrabuio DR, Leon EP, et al. Anti-ribosomal P protein: a novel antibody in autoimmune hepatitis. Liver Int. 2013;33:909–13.
Carmona-Fernandes D, Santos MJ, Canhao H, Fonseca JE. Anti-ribosomal P protein IgG autoantibodies in patients with systemic lupus erythematosus: diagnostic performance and clinical profile. BMC Med. 2013;11:98.
Dolezal JM, Dash AP, Prochownik EV. Diagnostic and prognostic implications of ribosomal protein transcript expression patterns in human cancers. BMC Cancer. 2018;18:275.
Loging WT, Reisman D. Elevated expression of ribosomal protein genes L37, RPP-1, and S2 in the presence of mutant p53. Cancer Epidemiol Prev Biomarkers. 1999;8:1011–6.
Pimentel-Santos FM, Ligeiro D, Matos M, Mourão AF, Costa J, Santos H, et al. Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects. Arthritis Res Therapy. 2011;13:R57.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47-e.
Mozhgani SH, Zarei-Ghobadi M, Teymoori-Rad M, Mokhtari-Azad T, Mirzaie M, Sheikhi M, et al. Human T-lymphotropic virus 1 (HTLV-1) pathogenesis: a systems virology study. J Cell Biochem. 2018;119:3968–79.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102:15545–50.
Goldberg G NM, Chatfield S, Murphy J, et al. White blood cells from rheumatoid arthritis patients and matched healthy donors 2018. Available from: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE117769.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14:417.
Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Research. 2015;4.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Linden SVD, Valkenburg HA, Cats A. Evaluation of diagnostic criteria for ankylosing spondylitis. Arthritis Rheum. 1984;27:361–8.
Van der Heijde D, Van der Linden S, Bellamy N, Calin A, Dougados M, Khan M, et al. Which domains should be included in a core set for endpoints in ankylosing spondylitis? Introduction to the ankylosing spondylitis module of OMERACT IV. J Rheumatol. 1999;26:945–7.
Garrett S, Jenkinson T, Kennedy LG, Whitelock H, Gaisford P, Calin A. A new approach to defining disease status in ankylosing spondylitis: the Bath ankylosing spondylitis disease activity index. J Rheumatol. 1994;21:2286–91.
Jones S, Porter J, Garrett S, Kennedy L, Whitelock H, Calin A. A new scoring system for the Bath ankylosing spondylitis metrology index (BASMI). J Rheumatol. 1995;22:1609.
Calin A, Garrett S, Whitelock H, Kennedy L, O'hea J, Mallorie P, et al. A new approach to defining functional ability in ankylosing spondylitis: the development of the Bath ankylosing spondylitis functional index. J Rheumatol. 1994;21:2281–5.
Penzo M, Montanaro L, Treré D, Derenzini M. The ribosome biogenesis—Cancer connection. Cells. 2019;8:55.
Zhou X, Liao W-J, Liao J-M, Liao P, Lu H. Ribosomal proteins: functions beyond the ribosome. J Mol Cell Biol. 2015;7:92–104.
Danilova N, Gazda HT. Ribosomopathies: how a common root can cause a tree of pathologies. Dis Model Mech. 2015;8:1013–26.
Armistead J, Triggs-Raine B. Diverse diseases from a ubiquitous process: the ribosomopathy paradox. FEBS Lett. 2014;588:1491–500.
Ruggero D, Pandolfi PP. Does the ribosome translate cancer? Nat Rev Cancer. 2003;3:179–92.
Hong M, Kim H, Kim I. Ribosomal protein L19 overexpression activates the unfolded protein response and sensitizes MCF7 breast cancer cells to endoplasmic reticulum stress-induced cell death. Biochem Biophys Res Commun. 2014;450:673–8.
Yang J, Chen Z, Liu N, Chen Y. Ribosomal protein L10 in mitochondria serves as a regulator for ROS level in pancreatic cancer cells. Redox Biol. 2018;19:158–65.
Barnard GF, Staniunas RJ, Mori M, Puder M, Jessup MJ, Steele GD, et al. Gastric and hepatocellular carcinomas do not overexpress the same ribosomal protein messenger RNAs as colonic carcinoma. Cancer Res. 1993;53:4048–52.
Shi Y, Zhai H, Wang X, Han Z, Liu C, Lan M, et al. Ribosomal proteins S13 and L23 promote multidrug resistance in gastric cancer cells by suppressing drug-induced apoptosis. Exp Cell Res. 2004;296:337–46.
Arthurs C, Murtaza BN, Thomson C, Dickens K, Henrique R, Patel HR, et al. Expression of ribosomal proteins in normal and cancerous human prostate tissue. PLoS One. 2017;12.
Zhou H, Wang Y, Lv Q, Zhang J, Wang Q, Gao F, et al. Overexpression of ribosomal RNA in the development of human cervical cancer is associated with rDNA promoter hypomethylation. PLoS One. 2016;11.
Tsoi H, Lam K, Dong Y, Zhang X, Lee C, Zhang J, et al. Pre-45s rRNA promotes colon cancer and is associated with poor survival of CRC patients. Oncogene. 2017;36:6109–18.
Luan M, Shang Z, Teng Y, Chen X, Zhang M, Lv H, et al. The shared and specific mechanism of four autoimmune diseases. Oncotarget. 2017;8:108355–74.
Absi M, La Vergne JP, Marzouki A, Giraud F, Rigal D, Reboud AM, et al. Heterogeneity of ribosomal autoantibodies from human, murine and canine connective tissue diseases. Immunol Lett. 1989;23:35–41.
Philip M, Schietinger A, Schreiber H. Ribosomal versus non-ribosomal cellular antigens: factors determining efficiency of indirect presentation to CD4+ T cells. Immunology. 2010;130:494–503.
Gracey E, Yao Y, Green B, Qaiyum Z, Baglaenko Y, Lin A, et al. Sexual dimorphism in the Th17 signature of ankylosing spondylitis. Arthritis Rheum. 2016;68:679–89.
Lee YH, Song GG. Meta-analysis of differentially expressed genes in ankylosing spondylitis. Genet Mol Res. 2015;14:5161–70.
Zhao H, Wang D, Fu D, Xue L. Predicting the potential ankylosing spondylitis-related genes utilizing bioinformatics approaches. Rheumatol Int. 2015;35:973–9.
Smolock EM, Korshunov VA, Glazko G, Qiu X, Gerloff J, Berk BC. Ribosomal protein L17, RpL17, is an inhibitor of vascular smooth muscle growth and carotid intima formation. Circulation. 2012;126:2418–27.
Bartels M, Bierings M. How I manage children with diamond-Blackfan anaemia. Br J Haematol. 2019;184:123–33.
Ito Y, Hashimoto M, Hirota K, Ohkura N, Morikawa H, Nishikawa H, et al. Detection of T cell responses to a ubiquitous cellular protein in autoimmune disease. Science. 2014;346:363–8.
Shipman L. Autoimmunity: troublemaker T cells target ribosomal protein. Nat Rev Immunol. 2014;14:776.
von Mikecz A, Hemmerich P, Peter H-H, Krawinkel U. Characterization of eukaryotic protein L7 as a novel autoantigen which frequently elicits an immune response in patients suffering from systemic autoimmune disease. Immunobiology. 1994;192:137–54.
Yu S, Lemos B. A portrait of ribosomal DNA contacts with hi-C reveals 5S and 45S rDNA anchoring points in the folded human genome. Genome Biol Evol. 2016;8:3545–58.
Elela SA, Nazar RN. Role of the 5.8 S rRNA in ribosome translocation. Nucleic Acids Res. 1997;25:1788–94.
Fontoura B, Sorokina EA, David E, Carroll RB. p53 is covalently linked to 5.8 S rRNA. Mol Cell Biol. 1992;12:5145–51.
Meyer A, Todt C, Mikkelsen NT, Lieb B. Fast evolving 18S rRNA sequences from Solenogastres (Mollusca) resist standard PCR amplification and give new insights into mollusk substitution rate heterogeneity. BMC Evol Biol. 2010;10:70.
Warda AS, Freytag B, Haag S, Sloan KE, Görlich D, Bohnsack MT. Effects of the Bowen-Conradi syndrome mutation in EMG1 on its nuclear import, stability and nucleolar recruitment. Hum Mol Genet. 2016;25:5353–64.
Armistead J, Khatkar S, Meyer B, Mark BL, Patel N, Coghlan G, et al. Mutation of a gene essential for ribosome biogenesis, EMG1, causes Bowen-Conradi syndrome. Am J Hum Genet. 2009;84:728–39.
Armistead J, Patel N, Wu X, Hemming R, Chowdhury B, Basra GS, et al. Growth arrest in the ribosomopathy, Bowen–Conradi syndrome, is due to dramatically reduced cell proliferation and a defect in mitotic progression. Biochim Biophys Acta (BBA) - Mol Basis Dis. 2015;1852:1029–37.
Delpu Y, Larrieu D, Gayral M, Arvanitis D, Dufresne M, Cordelier P, et al. Chapter 12 - noncoding RNAs: clinical and therapeutic applications. In: Egger G, Arimondo P, editors. Drug discovery in Cancer epigenetics. Boston: Academic Press; 2016. p. 305–26.
Pink RC, Wicks K, Caley DP, Punch EK, Jacobs L, Carter DR. Pseudogenes: pseudo-functional or key regulators in health and disease? RNA. 2011;17:792–8.
Rapicavoli NA, Qu K, Zhang J, Mikhail M, Laberge R-M, Chang HY. A mammalian pseudogene lncRNA at the interface of inflammation and anti-inflammatory therapeutics. eLife. 2013;2:e00762-e.
Chiang JJ, Sparrer KMJ, van Gent M, Lassig C, Huang T, Osterrieder N, et al. Viral unmasking of cellular 5S rRNA pseudogene transcripts induces RIG-I-mediated immunity. Nat Immunol. 2018;19:53–62.
Mazzone R, Zwergel C, Artico M, Taurone S, Ralli M, Greco A, et al. The emerging role of epigenetics in human autoimmune disorders. Clin Epigenetics. 2019;11:1–15.
Javierre BM, Hernando H, Ballestar E. Environmental triggers and epigenetic deregulation in autoimmune disease. Discov Med. 2011;12:535–45.
Whyte JM, Ellis JJ, Brown MA, Kenna TJ. Best practices in DNA methylation: lessons from inflammatory bowel disease, psoriasis and ankylosing spondylitis. Arthritis Res Therapy. 2019;21:133.
Bianco C, Mohr I. Ribosome biogenesis restricts innate immune responses to virus infection and DNA. Elife. 2019;8:e49551.
Miller SA, Brown AJ, Farach-Carson MC, Kirn-Safran CB. HIP/RPL29 down-regulation accompanies terminal chondrocyte differentiation. Differentiation. 2003;71:322–36.
Oristian DS, Sloofman LG, Zhou X, Wang L, Farach-Carson MC, Kirn-Safran CB. Ribosomal protein L29/HIP deficiency delays osteogenesis and increases fragility of adult bone in mice. J Orthop Res. 2009;27:28–35.
Chang CC, Chang CW, Nguyen PAA, Chang TH, Shih YL, Chang WY, et al. Ankylosing spondylitis and the risk of cancer. Oncol Lett. 2017;14:1315–22.
Nurmohamed MT, van der Horst-Bruinsma I, Maksymowych WP. Cardiovascular and cerebrovascular diseases in ankylosing spondylitis: current insights. Curr Rheumatol Rep. 2012;14:415–21.
Nakamura S, Suzuki K, Iijima H, Hata Y, Lim CR, Ishizawa Y, et al. Identification of baseline gene expression signatures predicting therapeutic responses to three biologic agents in rheumatoid arthritis: a retrospective observational study. Arthritis Res Ther. 2016;18:159.
Rosenbaum J, Chandran V. Management of comorbidities in ankylosing spondylitis. Am J Med Sci. 2012;343:364–6.
Role of the funding source
The funder of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.
This study was supported by a grant from the Pasteur Institute of Iran (Grant No. ZS-9503), and a grant from the Deputy of Research, Tehran University of Medical Sciences (Grant No. 98–3–96-45657).
Ethics approval and consent to participate
This study was performed based on the Declaration of Helsinki guidelines and was approved by the Ethics Committee of Tehran University of Medical Sciences (Approval No: IR.TUMS.VCR.REC.1398.1000. The written informed consent was signed by all participants before enrolling in the study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lari, A., Pourbadie, H.G., Sharifi-Zarchi, A. et al. Dysregulation of ribosome-related genes in ankylosing spondylitis: a systems biology approach and experimental method. BMC Musculoskelet Disord 22, 789 (2021). https://doi.org/10.1186/s12891-021-04662-2
- Ankylosing spondylitis
- Microarray analysis