Dataset Information

Identification of immune related cells and crucial genes in the peripheral blood of ankylosing spondylitis by integrated bioinformatics analysis.

ABSTRACT:

Background

Ankylosing spondylitis (AS) is a progressive rheumatic disease and studies reveal that the immune system is critical for the pathogenesis of AS. In the present study, various bioinformatics analysis methods were comprehensively applied, designed to identify potential key genes and inflammation states of AS.

Methods

The transcriptome profiles of GSE25101 and GSE73754 obtained from the Gene Expression Omnibus (GEO) database were merged for subsequent analyses. The differentially expressed genes (DEGs) were identified using the Bioconductor package Limma and threshold values. Functional enrichment and pathway enrichment analyses were performed using the clusterProfiler package and Gene Set Enrichment Analysis (GSEA). Next, protein-protein interaction (PPI) network of the identified DEGs was constructed by the online database, the Search Tool for the Retrieval of Interacting Genes (STRING), visualization and analysis were performed through Cytoscape software. Subsequently, we applied CIBERSORT algorithm to identify subpopulation proportions of immune cells in peripheral blood samples. Finally, we validated the hub genes with the GSE18781 dataset. Samples were collected from patients to validate gene and protein expression using qRT-PCR and ELISA.

Results

A total of 334 DEGs were identified, including 182 upregulated and 152 downregulated DEGs, between AS patients and normal human controls, which were primarily involved in immune response, autophagy, and natural killer cell-mediated cytotoxicity. The most prominent module and candidate biomarkers were identified from the PPI network. Biomarkers were selected for validation and their expressions were significantly decreased in peripheral blood samples which was consistent with transcriptome sequencing results. Nine genes with AUC > 0.70 were considered to be AS hub genes for ROC curve analysis, including GZMA, GZMK, PRF1, GNLY, NKG7, KLRB1, KLRD1, IL2RB and CD247. Furthermore, CIBERSORT results suggest that AS contained a higher proportion of CD8+ T cells, naive CD4+ T cells, neutrophils, and lower levels of gamma delta T cells compared with the normal controls.

Conclusion

In this study, we identified DEGs combined with their closely related biological functions and propose that granule-associated proteins and immune infiltration maybe involved in the progression of ankylosing spondylitis. These validated hub genes may provide new perspectives for understanding the molecular mechanisms of ankylosing spondylitis.

SUBMITTER: Zheng Y

PROVIDER: S-EPMC8432305 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundCircular RNA (circRNA) is a type of closed circular noncoding RNA (ncRNA), mostly formed by back-splicing or alternative splicing of pre-messenger RNA (mRNA). The aim of this study was to explore the expression profile of circRNA in peripheral blood mononuclear cells (PBMCs) of patients with ankylosing spondylitis (AS) and discover potential molecular markers of AS.MethodsThe circRNA microarray technology was used to detect the expression of circRNAs in the peripheral blood of 6 patients with AS and 6 healthy controls (HC). To screen the differentially expressed circRNAs by fold change (FC) and P value, these differentially expressed circRNAs were analyzed by bioinformatics. In 60 cases of AS and 30 cases of HC, 4 circRNAs were subjected to real-time fluorescence quantitative polymerase chain reaction (RT-qPCR), and their correlation with various clinical indicators was analyzed. Finally, the receiver operating characteristic (ROC) curve was used to analyze their potential as AS diagnostic markers.ResultsThe microarray results showed that there were 1369 significantly differently expressed (P < 0.05, FC > 1.5) circRNAs between the AS and HC groups (675 upregulated and 694 downregulated). The results of bioinformatics analysis suggested that they were mainly involved in "enzyme binding," "adenosine ribonucleotide binding," "MAPK signaling pathway", etc. The RT-qPCR results showed that the expressions of hsa_circRNA_001544 (U = 486.5, P < 0.05) and hsa_circRNA_102532 (U = 645, P < 0.05) were significantly different between the AS group and the HC group. The AS group was further divided into two subgroups: active AS (ASA) and stable AS (ASS). After analysis, it was found that compared with the HC group, hsa_circRNA_001544 was significantly increased in both ASA (U = 214, P < 0.05) and ASS groups (U = 273, P < 0.05), while hsa_circRNA_008961 (U = 250, P < 0.05) and hsa_circRNA_102532 (U = 295, P < 0.05) were only significantly increased in the ASA group. Furthermore, hsa_circRNA_012732 was significantly different between the ASA and ASS groups (U = 194, P < 0.05), and there was no statistical significance among the remaining groups. Correlation analysis results showed that hsa_circRNA_012732 was negatively correlated with Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), high-sensitivity C-reactive protein (hsCRP), and globulin (GLOB) and positively correlated with lymphocyte count (LY), mean corpusular volume, and albumin (ALB), and hsa_circRNA_008961 was negatively correlated with platelet (PLT) count. ROC curve analysis showed that hsa_circRNA_001544 (95% CI = 0.610-0.831, P < 0.05) and hsa_circRNA_102532 (95% CI = 0.521-0.762, P < 0.05) were statistically significant, and their area under curve (AUC) values were 0.720 and 0.642, respectively.ConclusionsThere are differentially expressed circRNAs in PBMCs of AS patients, and they may be involved in the occurrence and development of AS. Among these differentially expressed circRNAs, hsa_circRNA_012732 has the potential to become an indicator of disease activity, and hsa_circRNA_001544 has the potential to become a molecular marker for AS diagnosis.

Project description:Background Long non-coding RNAs (lncRNAs) are important regulators in ankylosing spondylitis (AS). Few studies have examined the lncRNA-RNA binding protein (RBP) interaction in AS. This study performed bioinformatics analysis and clinical verification to identify key lncRNAs and propose their RBP interaction. Methods Three GEO datasets of AS were analyzed by differential expression analysis. The differentially expressed lncRNAs between the AS and control groups were screened out, and the intersecting lncRNAs were regarded as target lncRNAs. Functional was performed to identify target lncRNAs by enrichment analysis, co-expressed RNA analysis, and lncRNA-RBP interaction analysis. Finally, this study analyzed the differential expression level and clinical value of lncRNAs between the AS and control groups. Results Linc00304, linc00926, and MIAT were differentially expressed and upregulated. Enrichment analysis indicated that the key KEGG terms were the T-cell receptor signaling pathway and B-cell receptor signaling pathway. The key molecular function term was protein binding, and the key biological process term was adaptive immune response. In qRT-PCR results, 44 samples were validated. linc00304 expression was positively correlated with bath ankylosing spondylitis disease activity index (BASDAI), bath ankylosing spondylitis functional index (BASFI), erythrocyte sedimentation rate (ESR), and c-reactive protein (CRP). linc00926 expression was only positively correlated with ESR, whereas MIAT expression was positively correlated with BASFI, ESR, and CRP. Logistic regression revealed that linc00304, ESR, and CRP were the independent risk factors for BASDAI activation. The area under the curve (AUC) of serum linc00304 level in the diagnosis of AS was 0.687 (cutoff value: 0.413, specificity: 0.423, sensitivity: 0.900). AUC of linc00926 was 0.664 (cutoff value: 0.299, sensitivity: 0.882, specificity: 0.417). AUC of MIAT was 0.623 (cutoff value: 0.432, specificity: 0.443, sensitivity: 0.890) (all P <0.05). Conclusion Overall, this study uncovered three novel lncRNAs, which were upregulated in AS, and proposed a new lncRNA-RBP-mRNA interaction that might regulate adaptive immune response.

Project description:Several previous studies have attempted to investigate the regulatory mechanisms underlying gene expression in ankylosing spondylitis (AS). However, the specific molecular pathways underlying this condition remain unclear. Previous research used next-generation RNA sequencing to identify a series of differentially expressed genes (DEGs) in peripheral blood mononuclear cells (PBMCs) when compared between patients with AS and healthy controls, thus implying that these DEGs may be related to AS. Furthermore, by screening these DEGS, it may be possible to facilitate clinical diagnosis and optimize treatment strategies. In order to test this hypothesis, we recruited 15 patients with AS and 15 healthy controls. We randomly selected five subjects from each group of patients for RNA sequencing analysis. Sequence reads were generated by an Illumina HiSeq2500 platform and mapped on to the human reference genome using HISAT2. We successfully identified 973 significant DEGs (p < 0.05) in PBMCs. When compared with controls, 644 of these genes were upregulated (with a fold?change?(FC) > 2) in AS patients and 329 were downregulated (FC < 0.5). Our analysis identified numerous genes related to immune response. Gene Ontology (GO) analysis indicated that these DEGs were significantly related to the positive regulation of epidermal growth factor-activated receptor activity, the positive regulation of the ERBB (erb-b2 receptor tyrosine kinase) signaling pathway, the differentiation of trophoblast giant cells, oxygen transport, immune-related pathways, and inflammation-related pathways. The DEGs were also closely related to the TNF and NF-?B signaling pathways. Six DEGs were verified by quantitative real-time polymerase chain reaction (qRT-PCR). Receiver operating characteristic (ROC) curve analysis indicated that IL6 may represent a useful biomarker for diagnosing AS. The development of new biomarkers may help us to elucidate the specific mechanisms involved in the development and progression of AS.