Autophagy dysfunction may be involved in the pathogenesis of ankylosing spondylitis

Identification of sex-specific biomarkers related to programmed cell death and analysis of immune cells in ankylosing spondylitis

Ankylosing spondylitis (AS) stands as a persistent inflammatory ailment predominantly impacting the axial skeleton, with the immune system and inflammation intricately entwined in its pathogenesis. This study endeavors to elucidate gender-specific patterns in immune cell infiltration and diverse forms of cell demise within the AS milieu. The aim is to refine the diagnosis and treatment of gender-specific AS patients, thereby advancing patient outcomes. In the pursuit of our investigation, two datasets (GSE25101 and GSE73754) pertinent to ankylosing spondylitis (AS) were meticulously collected and normalized from the GEO database. Employing the CIBERSORT algorithm, we conducted a comprehensive analysis of immune cell infiltration across distinct demographic groups and genders. Subsequently, we discerned differentially expressed genes (DEGs) associated with various cell death modalities in AS patients and their healthy counterparts. Our focus extended specifically to ferroptosis-related DEGs (FRDEGs), cuproptosis-related DEGs (CRDEGs), anoikis-related DEGs (ARDEGs), autophagy-related DEGs (AURDEGs), and pyroptosis-related DEGs (PRDEGs). Further scrutiny involved discerning disparities in these DEGs between AS patients and healthy controls, as well as disparities between male and female patients. Leveraging machine learning (ML) methodologies, we formulated disease prediction models employing cell death-related DEGs (CDRDEGs) and identified biomarkers intertwined with cell death in AS. Relative to healthy controls, a myriad of differentially expressed genes (DEGs) linked to cell death surfaced in AS patients. Among AS patients, 82 FRDEGs, 29 CRDEGs, 54 AURDEGs, 21 ARDEGs, and 74 PRDEGs were identified. In male AS patients, these numbers were 78, 33, 55, 24, and 94, respectively. Female AS patients exhibited 66, 41, 40, 17, and 82 DEGs in the corresponding categories. Additionally, 36 FRDEGs, 14 CRDEGs, 19 AURDEGs, 10 ARDEGs, and 36 PRDEGs exhibited differential expression between male and female AS patients. Employing machine learning techniques, LASSO, RF, and SVM-RFE were employed to discern key DEGs related to cell death (CDRDDEGs). The six pivotal CDRDDEGs in AS patients, healthy controls, were identified as CLIC4, BIRC2, MATK, PKN2, SLC25A5, and EDEM1. For male AS patients, the three crucial CDRDDEGs were EDEM1, MAP3K11, and TRIM21, whereas for female AS patients, COX7B, PEX2, and RHEB took precedence. Furthermore, the trio of DDX3X, CAPNS1, and TMSB4Y emerged as the key CDRDDEGs distinguishing between male and female AS patients. In the realm of immune correlation, the immune infiltration abundance in female patients mirrored that of healthy controls. Notably, key genes exhibited a positive correlation with T-cell CD4 memory activation when comparing male and female patient samples. This study engenders a more profound comprehension of the molecular underpinnings governing immune cell infiltration and cell death in ankylosing spondylitis (AS). Furthermore, the discernment of gender-specific disparities among AS patients underscores the clinical significance of these findings. By identifying DEGs associated with diverse cell death modalities, this study proffers invaluable insights into potential clinical targets for AS patients, taking cognizance of gender-specific nuances. The identification of gender-specific biological targets lays the groundwork for the development of tailored diagnostic and therapeutic strategies, heralding a pivotal step toward personalized care for AS patients.

Identification of sex-specific biomarkers related to programmed cell death and analysis of immune cells in ankylosing spondylitis

Ankylosing spondylitis (AS) stands as a persistent inflammatory ailment predominantly impacting the axial skeleton, with the immune system and inflammation intricately entwined in its pathogenesis. This study endeavors to elucidate gender-specific patterns in immune cell infiltration and diverse forms of cell demise within the AS milieu. The aim is to refine the diagnosis and treatment of gender-specific AS patients, thereby advancing patient outcomes. In the pursuit of our investigation, two datasets (GSE25101 and GSE73754) pertinent to ankylosing spondylitis (AS) were meticulously collected and normalized from the GEO database. Employing the CIBERSORT algorithm, we conducted a comprehensive analysis of immune cell infiltration across distinct demographic groups and genders. Subsequently, we discerned differentially expressed genes (DEGs) associated with various cell death modalities in AS patients and their healthy counterparts. Our focus extended specifically to ferroptosis-related DEGs (FRDEGs), cuproptosis-related DEGs (CRDEGs), anoikis-related DEGs (ARDEGs), autophagy-related DEGs (AURDEGs), and pyroptosis-related DEGs (PRDEGs). Further scrutiny involved discerning disparities in these DEGs between AS patients and healthy controls, as well as disparities between male and female patients. Leveraging machine learning (ML) methodologies, we formulated disease prediction models employing cell death-related DEGs (CDRDEGs) and identified biomarkers intertwined with cell death in AS. Relative to healthy controls, a myriad of differentially expressed genes (DEGs) linked to cell death surfaced in AS patients. Among AS patients, 82 FRDEGs, 29 CRDEGs, 54 AURDEGs, 21 ARDEGs, and 74 PRDEGs were identified. In male AS patients, these numbers were 78, 33, 55, 24, and 94, respectively. Female AS patients exhibited 66, 41, 40, 17, and 82 DEGs in the corresponding categories. Additionally, 36 FRDEGs, 14 CRDEGs, 19 AURDEGs, 10 ARDEGs, and 36 PRDEGs exhibited differential expression between male and female AS patients. Employing machine learning techniques, LASSO, RF, and SVM-RFE were employed to discern key DEGs related to cell death (CDRDDEGs). The six pivotal CDRDDEGs in AS patients, healthy controls, were identified as CLIC4, BIRC2, MATK, PKN2, SLC25A5, and EDEM1. For male AS patients, the three crucial CDRDDEGs were EDEM1, MAP3K11, and TRIM21, whereas for female AS patients, COX7B, PEX2, and RHEB took precedence. Furthermore, the trio of DDX3X, CAPNS1, and TMSB4Y emerged as the key CDRDDEGs distinguishing between male and female AS patients. In the realm of immune correlation, the immune infiltration abundance in female patients mirrored that of healthy controls. Notably, key genes exhibited a positive correlation with T-cell CD4 memory activation when comparing male and female patient samples. This study engenders a more profound comprehension of the molecular underpinnings governing immune cell infiltration and cell death in ankylosing spondylitis (AS). Furthermore, the discernment of gender-specific disparities among AS patients underscores the clinical significance of these findings. By identifying DEGs associated with diverse cell death modalities, this study proffers invaluable insights into potential clinical targets for AS patients, taking cognizance of gender-specific nuances. The identification of gender-specific biological targets lays the groundwork for the development of tailored diagnostic and therapeutic strategies, heralding a pivotal step toward personalized care for AS patients.

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

OBJECTIVES

The aim of this study was to explore the long non-coding RNA (lncRNA) expression profiles in serum of patients with ankylosing spondylitis (AS). The role of these lncRNAs in this complex autoimmune situation needs to be evaluated.

METHODS

We used high-throughput whole-transcriptome sequencing to generate sequencing data from three patients with AS and three normal controls (NC). Then, we performed bioinformatics analyses to identify the functional and biological processes associated with differentially expressed lncRNAs (DElncRNAs). We confirmed the validity of our RNA-seq data by assessing the expression of eight lncRNAs via quantitative reverse transcription polymerase chain reaction (qRT-PCR) in 20 AS and 20 NC samples. We measured the correlation between the expression levels of lncRNAs and patient clinical index values using the Spearman correlation test.

RESULTS

We identified 72 significantly upregulated and 73 significantly downregulated lncRNAs in AS patients compared to NC. qRT-PCR was performed to validate the expression of selected DElncRNAs; the results demonstrated that the expression levels of MALAT1:24, NBR2:9, lnc-DLK1-35:13, lnc-LARP1-1:1, lnc-AIPL1-1:7, and lnc-SLC12A7-1:16 were consistent with the sequencing analysis results. Enrichment analysis showed that DElncRNAs mainly participated in the immune and inflammatory responses pathways, such as regulation of protein ubiquitination, major histocompatibility complex class I-mediated antigen processing and presentation, MAPkinase activation, and interleukin-17 signaling pathways. In addition, a competing endogenous RNA network was constructed to determine the interaction among the lncRNAs, microRNAs, and mRNAs based on the confirmed lncRNAs (MALAT1:24 and NBR2:9). We further found the expression of MALAT1:24 and NBR2:9 to be positively correlated with disease severity.

CONCLUSION

Taken together, our study presents a comprehensive overview of lncRNAs in the serum of AS patients, thereby contributing novel perspectives on the underlying pathogenic mechanisms of this condition. In addition, our study predicted MALAT1 has the potential to be deeply involved in the pathogenesis of AS.

A review of long non-coding RNAs in ankylosing spondylitis: pathogenesis, clinical assessment, and therapeutic targets

Ankylosing spondylitis (AS) is a chronic immune-mediated type of inflammatory arthritis characterized by inflammation, bone erosion, and stiffness of the spine and sacroiliac joints. Despite great efforts put into the investigation of the disease, the pathogenesis of AS remains unclear, posing challenges in identifying ideal targets for diagnosis and treatment. To enhance our understanding of AS, an increasing number of studies have been conducted. Some of these studies reveal that long non-coding RNAs (lncRNAs) play crucial roles in the etiology of AS. Some certain lncRNAs influence the development of AS by regulating inflammatory responses, autophagy, apoptosis, and adipogenesis, as well as the proliferation and differentiation of cells. Additionally, some lncRNAs demonstrate potential as biomarkers, aiding in monitoring disease progression and predicting prognosis. In this review, we summarize recent studies concerning lncRNAs in AS to elucidate the underlying mechanisms in which lncRNAs are involved and their potential values as biomarkers for disease assessment and druggable targets for therapy.

Autophagy and Senescence: The Molecular Mechanisms and Implications in Liver Diseases

The liver is the primary organ accountable for complex physiological functions, including lipid metabolism, toxic chemical degradation, bile acid synthesis, and glucose metabolism. Liver function homeostasis is essential for the stability of bodily functions and is involved in the complex regulation of the balance between cell proliferation and cell death. Cell proliferation-halting mechanisms, including autophagy and senescence, are implicated in the development of several liver diseases, such as cholestasis, viral hepatitis, nonalcoholic fatty liver disease, liver fibrosis, and hepatocellular carcinoma. Among various cell death mechanisms, autophagy is a highly conserved and self-degradative cellular process that recycles damaged organelles, cellular debris, and proteins. This process also provides the substrate for further metabolism. A defect in the autophagy machinery can lead to premature diseases, accelerated aging, inflammatory state, tumorigenesis, and cellular senescence. Senescence, another cell death type, is an active player in eliminating premalignant cells. At the same time, senescent cells can affect the function of neighboring cells by secreting the senescence-associated secretory phenotype and induce paracrine senescence. Autophagy can promote and delay cellular senescence under different contexts. This review decodes the roles of autophagy and senescence in multiple liver diseases to achieve a better understanding of the regulatory mechanisms and implications of autophagy and senescence in various liver diseases.

Acrylamide induces human chondrocyte cell death by initiating autophagy‑dependent ferroptosis

Acrylamide (ACR) is formed during heat treatment of foodstuffs and ACR may serve as a probable malignant neoplastic disease agent in all organs and tissues of the human body. However, it is unknown if ACR is associated with ankylosing spondylitis (AS) pathogenesis. Cell viability and proliferation were determined using CCK-8 assay and EdU staining. Flow cytometry was used to determine cell death and cell cycle arrest. Intracellular lipid reactive oxygen species, Fe²⁺ and mitochondrial membrane potential (MMP) were analyzed using a C11-BODIPY581/591 fluorescent probe, FerroOrange staining and a JC-1 MMP Assay kit, respectively. The present study showed that ACR decreased chondrocyte cell viability in a dose-dependent manner and that ACR significantly promoted chondrocyte senescence. ACR also elevated the expression of cell cycle arrest-associated proteins, including p53, cyclin-dependent kinase inhibitor 1 and cyclin-dependent kinase inhibitor protein, in human chondrocytes. Similarly, DNA damage was also enhanced following ACR treatment in chondrocytes. In addition, the ferroptosis-specific inhibitor ferrostatin-1 (Fer-1) and the autophagy inhibitor 3-methyladenine abolished ACR-induced cell death in chondrocytes. ACR was shown to activate autophagic flux and induce mitochondrial dysfunction by increasing the MMP. Western blot analysis of ferroptosis-related proteins demonstrated that ACR decreased the expression of glutathione peroxidase 4, solute carrier family 7 member 11, transferrin receptor protein 1 and ferritin heavy chain 1 in chondrocytes whereas Fer-1 abolished these effects. ACR treatment significantly elevated the phosphorylation levels of AMP-activated protein kinase (AMPK) and serine/threonine-protein kinase ULK1 in human chondrocytes. Notably, the effect of ACR was diminished by knockdown of AMPK, as evidenced by reduced lipid reactive oxygen species accumulation and Fe²⁺ levels. Hence, ACR inhibited cell proliferation and contributed to cell death by inducing autophagy-dependent ferroptosis while promoting autophagy by activating AMPK-ULK1-mTOR signaling in human chondrocytes. It was hypothesized that the presence of ACR in foodstuffs may increase the risk of AS and that decreasing ACR in food products is of importance.

RNA sequencing and bioinformatics analysis of differentially expressed genes in the peripheral serum of ankylosing spondylitis patients

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.

Joint together: The etiology and pathogenesis of ankylosing spondylitis

Spondyloarthritis (SpA) refers to a group of diseases with inflammation in joints and spines. In this family, ankylosing spondylitis (AS) is a rare but classic form that mainly involves the spine and sacroiliac joint, leading to the loss of flexibility and fusion of the spine. Compared to other diseases in SpA, AS has a very distinct hereditary disposition and pattern of involvement, and several hypotheses about its etiopathogenesis have been proposed. In spite of significant advances made in Th17 dynamics and AS treatment, the underlying mechanism remains concealed. To this end, we covered several topics, including the nature of the immune response, the microenvironment in the articulation that is behind the disease's progression, and the split between the hypotheses and the evidence on how the intestine affects arthritis. In this review, we describe the current findings of AS and SpA, with the aim of providing an integrated view of the initiation of inflammation and the development of the disease.

New Insights into the Regulatory Role of Ferroptosis in Ankylosing Spondylitis via Consensus Clustering of Ferroptosis-Related Genes and Weighted Gene Co-Expression Network Analysis

Background: The pathogenesis of ankylosing spondylitis (AS) remains undetermined. Ferroptosis is a newly discovered form of regulated cell death involved in multiple autoimmune diseases. Currently, there are no reports on the connection between ferroptosis and AS. Methods: AS samples from the Gene Expression Omnibus were divided into two subgroups using consensus clustering of ferroptosis-related genes (FRGs). Weighted gene co-expression network analysis (WGCNA) of the intergroup differentially expressed genes (DEGs) and protein–protein interaction (PPI) analysis of the key module were used to screen out hub genes. A multifactor regulatory network was then constructed based on hub genes. Results: The 52 AS patients in dataset GSE73754 were divided into cluster 1 (n = 24) and cluster 2 (n = 28). DEGs were mainly enriched in pathways related to mitochondria, ubiquitin, and neurodegeneration. Candidate hub genes, screened by PPI and WGCNA, were intersected. Subsequently, 12 overlapping genes were identified as definitive hub genes. A multifactor interaction network with 45 nodes and 150 edges was generated, comprising the 12 hub genes and 32 non-coding RNAs. Conclusions: AS can be divided into two subtypes according to FRG expression. Ferroptosis might play a regulatory role in AS. Tailoring treatment according to the ferroptosis status of AS patients can be a promising direction.

DNA methylation and transcription of the FOXO3a gene are associated with ankylosing spondylitis

To explore the association between methylation level and transcript level of Forkhead box O3a (FOXO3a) gene with ankylosing spondylitis (AS) susceptibility. Methylation levels of the FOXO3a promoter were measured in 84 AS patients and 83 healthy controls. A total of 77 patients and 66 healthy subjects were included in subsequent mRNA level testing. DNA methylation levels of 107 CpG sites on 6 CpG islands in the FOXO3a gene were investigated. This study indicated that CpG-4 and CpG-5 islands were markedly hypomethylated in AS patients. The methylation level of CpG-4 island in AS patients was negatively correlated with erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), and Ankylosing Spondylitis Disease Activity Score (ASDAS). Moreover, FOXO3a mRNA levels were significantly decreased in AS patients and were obviously negatively correlated with the methylation levels of CpG-2 and CpG-5 islands in AS patients without treatment. The sensitivity and specificity of differential methylated CpG sites of FOXO3a were 74.7 and 85.4%, respectively. Besides, FOXO3a mRNA had a sensitivity of 80.0% and a specificity of 68.8%. DNA methylation and transcription of FOXO3a might be related to AS susceptibility and play a crucial role in the diagnosis of AS, but many open questions remain.

Bioinformatics Analysis of the Molecular Mechanism and Potential Treatment Target of Ankylosing Spondylitis

Ankylosing spondylitis (AS) is an autoimmune disease that mainly affects the spinal joints, sacroiliac joints, and adjacent soft tissues. We conducted bioinformatics analysis to explore the molecular mechanism related to AS pathogenesis and uncover novel potential molecular targets for the treatment of AS. The profiles of GSE25101, containing gene expression data extracted from the blood of 16 AS patients and 16 matched controls, were acquired from the Gene Expression Omnibus (GEO) database. The background correction and standardization were carried out utilizing the transcript per million (TPM) method. After analysis of AS patients and the normal groups, we identified 199 differentially expressed genes (DEGs) with upregulation and 121 DEGs with downregulation by the limma R package. The results of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) biological process enrichment analysis revealed that the DEGs with upregulation were mainly associated with spliceosome, ribosome, RNA-catabolic process, electron transport chain, etc. And the DEGs with downregulation primarily participated in T cell-associated pathways and processes. After analysis of the protein-protein interaction (PPI) network, our data revealed that the hub genes, comprising MRPL13, MRPL22, LSM3, COX7A2, COX7C, EP300, PTPRC, and CD4, could be the treatment targets in AS. Our data furnish new hints to uncover the features of AS and explore more promising treatment targets towards AS.

Aberrant antigen processing and presentation: Key pathogenic factors leading to immune activation in Ankylosing spondylitis

The strong association of HLA-B*27 with ankylosing spondylitis (AS) was first reported nearly 50 years ago. However, the mechanistic link between HLA-B*27 and AS has remained an enigma. While 85-90% of AS patients possess HLA-B*27, majority of HLA-B*27 healthy individuals do not develop AS. This suggests that additional genes and genetic regions interplay with HLA-B*27 to cause AS. Previous genome-wide association studies (GWAS) identified key genes that are distinctively expressed in AS, including the Endoplasmic Reticulum Aminopeptidase (ERAP) 1 and ERAP2. As these gene-encoding molecules are primarily implicated in the process of peptide processing and presentation, potential pathological interaction of these molecules with HLA-B*27 may operate to cause AS by activating downstream immune responses. The aberrant peptide processing also gives rise to the accumulation of unstable protein complex in endoplasmic reticulum (ER), which drives endoplasmic reticulum-associated protein degradation (ERAD) and unfolded protein response (UPR) and activates autophagy. In this review, we describe the current hypotheses of AS pathogenesis, focusing on antigen processing and presentation operated by HLA-B*27 and associated molecules that may contribute to the disease initiation and progression of AS.