Hub genes and their key effects on prognosis of Burkitt lymphoma

World J Clin Oncol. 2023 Oct 24;14(10):357-372. doi: 10.5306/wjco.v14.i10.357.

Abstract

Background: Burkitt lymphoma (BL) is an exceptionally aggressive malignant neoplasm that arises from either the germinal center or post-germinal center B cells. Patients with BL often present with rapid tumor growth and require high-intensity multi-drug therapy combined with adequate intrathecal chemotherapy prophylaxis, however, a standard treatment program for BL has not yet been established. It is important to identify biomarkers for predicting the prognosis of BLs and discriminating patients who might benefit from the therapy. Microarray data and sequencing information from public databases could offer opportunities for the discovery of new diagnostic or therapeutic targets.

Aim: To identify hub genes and perform gene ontology (GO) and survival analysis in BL.

Methods: Gene expression profiles and clinical traits of BL patients were collected from the Gene Expression Omnibus database. Weighted gene co-expression network analysis (WGCNA) was applied to construct gene co-expression modules, and the cytoHubba tool was used to find the hub genes. Then, the hub genes were analyzed using GO and Kyoto Encyclopedia of Genes and Genomes analysis. Additionally, a Protein-Protein Interaction network and a Genetic Interaction network were constructed. Prognostic candidate genes were identified through overall survival analysis. Finally, a nomogram was established to assess the predictive value of hub genes, and drug-gene interactions were also constructed.

Results: In this study, we obtained 8 modules through WGCNA analysis, and there was a significant correlation between the yellow module and age. Then we identified 10 hub genes (SRC, TLR4, CD40, STAT3, SELL, CXCL10, IL2RA, IL10RA, CCR7 and FCGR2B) by cytoHubba tool. Within these hubs, two genes were found to be associated with OS (CXCL10, P = 0.029 and IL2RA, P = 0.0066) by survival analysis. Additionally, we combined these two hub genes and age to build a nomogram. Moreover, the drugs related to IL2RA and CXCL10 might have a potential therapeutic role in relapsed and refractory BL.

Conclusion: From WGCNA and survival analysis, we identified CXCL10 and IL2RA that might be prognostic markers for BL.

Keywords: Burkitt lymphoma; Functional enrichment analysis; Microarray data; Prognosis; Therapeutic target; Weighted gene co-expression network analysis.