Clinical Features and Prognostic Impact of Coexpression Modules Constructed by WGCNA for Diffuse Large B-Cell Lymphoma

Abstract

Objective. Diffuse large B-cell lymphoma (DLBCL) is a highly aggressive malignant tumor, accounting for 30-40% of non-Hodgkin’s lymphoma. Our aim was to construct novel prognostic models of candidate genes based on clinical features. Methods. RNA-seq and clinical data of DLBCL were retrieved from TCGA database. Coexpression modules were constructed by WGCNA. Then, we investigated the interactions between modules and clinical features. By overall survival analysis, prognostic candidate genes from modules of interest were identified. A coexpression network of prognostic candidate genes was then constructed through WGCNA. GEPIA was used to analyze the expression of a candidate gene between DLBCL and normal samples. Results. 19 coexpression modules were constructed by 12813 genes from 52 DLBCL samples. The number of genes in modules ranged from 34 to 5457. We found that the purple module was significantly related with histological type (p value = 1e-04). Overall survival analysis revealed that MAFA-AS1, hsa-mir-338, and hsa-mir-891a were related with prognosis of DLBCL (p value = 0.027, 0.039, and 0.022, respectively). A coexpression network was constructed for the three prognostic genes. MAFA-AS1 was interacted with 36 genes, hsa-mir-891a was interacted with 11 genes, while no gene showed interaction with hsa-mir-338. Using GEPIA, we found that MAFA-AS1 showed low expression in DLBCL samples (p<0.01). Conclusion. We constructed a coexpression module related with histological type and identified three candidate genes (MAFA-AS1, hsa-mir-338, and hsa-mir-891a) that possessed potential value as prognostic biomarkers and therapeutic targets of DLBCL.