Low CCL19 expression is associated with adverse clinical outcomes for follicular lymphoma patients treated with chemoimmunotherapy

Abstract Background This study aimed to recognize the hub genes associated with prognosis in follicular lymphoma (FL) treated with first-line rituximab combined with chemotherapy. Method RNA sequencing data of dataset GSE65135 (n = 24) were included in differentially expressed genes (DEGs) analysis. Weighted gene co-expression network analysis (WGCNA) was applied for exploring the coexpression network and identifying hub genes. Validation of hub genes expression and prognosis were applied in dataset GSE119214 (n = 137) and independent patient cohort from Cancer Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College (n = 32), respectively, by analyzing RNAseq expression data and serum protein concentration quantified by ELISA. The Gene Set Enrichment Analysis (GSEA), gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments analysis were performed. CIBERSORT was applied for tumor-infiltrating immune cells (TIICs) subset analysis. Results A total of 3260 DEGs were obtained, with 1861 genes upregulated and 1399 genes downregulated. Using WGCNA, eight hub genes, PLA2G2D, MMP9, PTGDS, CCL19, NFIB, YAP1, RGL1, and TIMP3 were identified. Kaplan–Meier analysis and multivariate COX regression analysis indicated that CCL19 independently associated with overall survival (OS) for FL patients treated with rituximab and chemotherapy (HR = 0.47, 95% CI [0.25–0.86], p = 0.014). Higher serum CCL19 concentration was associated with longer progression-free survival (PFS, p = 0.014) and OS (p = 0.039). TIICs subset analysis showed that CCL19 expression had a positive correlation with monocytes and macrophages M1, and a negative correlation with naïve B cells and plasma cells. Conclusion CCL19 expression was associated with survival outcomes and might be a potential prognostic biomarker for FL treated with first-line chemoimmunotherapy.