Objective: To determine the expression level of RAG1 and its clinical significance in myelodysplastic syndromes (MDS).
Methods: To explore the candidate genes, the microarray datasets GSE19429, GSE58831, and GSE2779 were downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) in MDS were screened using RStudio, and overlapped DEGs were obtained with Venn Diagrams. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses, and protein-protein interaction network were performed. Quantitative real-time PCR (qRT-PCR) was employed to confirm the microarray results.
Results: This study identified 26 DEGs. Functional enrichment analyses indicated that these DEGs were significantly enriched in the immune response, and hematopoietic cell lineage. Eight core genes, for example, RAG1 and PAX5, were identified with a high degree of connectivity. The result of qRT-PCR showed that RAG1 was significantly down-regulated in MDS patients, which helped in distinguishing MDS patients from normal controls. The area under the curve of the receiver operator characteristic was 0.913 (P < 0.0001). MDS patients with low RAG1 expression level had a poor long-term survival (P = 0.031). What's more, the expression of RAG1 was significantly increased in the patients who received treatment.
Conclusion: The results showed that the expression of RAG1 was down-regulated in MDS patients. Lower RAG1 expression was associated with adverse clinical outcomes. RAG1 may be a potential prognostic biomarker for MDS.
Keywords: Bioinformatics; Clinical outcome; Myelodysplastic syndromes; Prognosis; RAG1.