Bioinformatics Analysis and Verification of Acute B-Lymphocytic Leukemia in Children with Isolated Bone Marrow Relapse

Abstract

To analyze the hub genes affecting the solely bone marrow relapse of childish acute B-cell lymphoblastic leukemia (B-ALL). The high-throughput RNA sequencing data were downloaded from TCGA database, the differentially expressed genes were screened by DESeq2 package of R, and the differentially expressed genes were grouped by GO function enrichment analysis and KEGG pathway enrichment analysis. Further, the data of STRING database and Cytoscape software were used to construct protein interaction network, screen hub genes and highly interaction protein sub network, perform GO and KEGG analysis of the hub genes and protein sub network respectively. JASPAR database was used to screen the upstream transcription factor of the hub gene promoter. Survival analysis based on the expression of hub genes was performed with clinical information attached to TCGA database. The bone marrow samples and clinical data of the patients were collected, the analysis results of hub genes were verified through clinical samples. 847 differentially expressed genes were collected, including 813 up-regulated genes, 34 down-regulated genes, 11 hub genes were screened out. The results of survival analysis showed that RPS5、RPS15、RPL23、RPL35、RPS8、RPS27A、RPS3、RPL9、RPS21、RPS7 and RPL38 showed significant effect on the survival of the children, and ZNF460 might be involved in their regulation. The high expressions of RPS3, RPS15, RPS8, RPS27A, and RPS21 had been verified in clinical samples of solely bone marrow relapsed patients. RPS3, RPS15, RPS8, RPS27A, RPS21 can be used as biomarkers to indicate the malignant event of solely bone marrow relapse, which may be regulated by ZNF460.