Identification of CXCR4 Upregulation in Diffuse Large B-Cell Lymphoma Associated with Prognostic Significance and Clinicopathological Characteristics

Background Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous malignant lymphoma with distinct characteristics. Patients with treatment failure after the standard immunochemotherapy have worse prognosis, which implies the necessity to uncover novel targets. The C-X-C chemokine receptor 4 (CXCR4) overexpression has been identified in several hematopoietic malignancies. However, the expression signatures and prognostic significance of CXCR4 in DLBCL associated with clinicopathological features remain unclear. Methods Gene expression profiles of DLBCL were obtained from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Then, a meta-analysis with an integrated bioinformatic analysis was performed to assess the relationship between CXCR4 expression and clinicopathological features of DLBCL. Finally, experimental verification including immunohistochemical (IHC) staining and real-time quantitative PCR (qPCR) was carried out using patient samples. In vitro cell line viability tests were conducted using CXCR4 inhibitor WZ811. Results DLBCL patients with activated B-cell-like (ABC) subtype have higher expression level of CXCR4 with worse survival. Differential expressed genes in the CXCR4-upregulation group were enriched in canonical pathways associated with oncogenesis. DLBCL with CXCR4 upregulation had lower degree of CD8+ T cell infiltration. TIMER analysis demonstrated that the CXCR4 expression was positively correlated with the expression of CD5, MYC, NOTCH1, PDCD1, CD274, mTOR, FOXO1, and hnRNPA2B1 in DLBCL. IHC study in patient samples showed the positive correlation between CXCR4 and nongerminal center B-cell (non-GCB) subtype and mTOR expression. Meanwhile, quantitative polymerase chain reaction results revealed that high CXCR4 mRNA level was correlated to double-hit DLBCL. Finally, cell viability test showed that WZ811 exerted antiproliferation effect in DLBCL cell lines in a dose-dependent manner. Conclusion CXCR4 was upregulated in ABC-DLBCL associated with worse prognosis. Our analysis predicted CXCR4 as a potential target for DLBCL treatment, which may serve as an inhibitor both on BCR signaling and nuclear export warranting further investigation in clinical trials.


Introduction
Diffuse large B-cell lymphoma (DLBCL) is one of the most common types of non-Hodgkin lymphoma (NHL), which is considered aggressive but potentially curable. Despite 60% patients can be cured with the standard immunochemotherapy of R-CHOP (rituximab plus cyclophosphamide, doxorubicin, vincristine and prednisolone), those with treat-ment failure after R-CHOP often have worse outcomes [1]. With the advance in second-generation sequencing to uncover novel molecular biomarkers and signaling pathways, there is a trend toward using XR-CHOP treatment, namely, adding new drugs to the standard R-CHOP regimen. For instance, bortezomib, lenalidomide, ibrutinib, BCL2, and EZH2 inhibitors, as well as inhibitors of PI3K/AKT/mTOR, have now being used as first-line therapy in combination with R-CHOP in patients with estimated poor prognosis according to pathological subtypes in several biomarker-driven clinical studies [2]. Nevertheless, there is limited evidence that DLBCL patients benefit more from XR-CHOP than R-CHOP [3,4]. Therefore, further research on the genetics and pathogenesis of DLBCL is needed to identify novel molecular targets and therapeutic vulnerabilities.
Chemokines and their receptors play an important role in oncogenesis, progression, and dissemination of tumor cells [5]. CXCR4 is a chemokine receptor involved in a series of biological processes including hematopoeisis and immune response, which is highly expressed in a variety of hematological malignancies [6]. CXCR4 is a prognostic biomarker in lots of cancer diseases such as gastrointestinal cancer [7], breast cancer [8], and lung cancer [9] [10,11]. Recent studies revealed that the CXCR4 and its ligand CXCL12 are linked to the pathogenesis of lymphohematopoietic malignancies [12,13]. The CXCR4 expression was reported to be associated with poor prognosis, and the inhibition of CXCR4 renders tumor cells more sensitive to chemotherapy [14].
To date, there has been no consensus on the correlation between the prognostic significance of CXCR4 and clinicopathological characteristics of DLBCL. In a large cohort (training/validation cohort of 468/275), the increased expression of CXCR4 was associated with the activated Bcell-like (ABC) subtype in patients with newly diagnosed DLBCL. However, it was also an independent risk factor for progression-free survival (PFS) in germinal center Bcell-like-(GCB-) DLBCL [15]. In another study by et al., CXCR4 was not associated with the prognosis of DLBCL patients [16]. Differed results seem controversial, and an integrated analysis is lacking. In this study, we used a meta-analysis and an integrated bioinformatics analysis to conduct a comprehensive assessment of the relationship between CXCR4 expression and DLBCL. Results were further verified by immunohistochemical (IHC) staining on tissue samples from independent cohorts with in vitro cell function tests.

Data Mining and Clinical Samples.
We searched DLBCL gene expression datasets from the Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/) and the Gene Expression Omnibus (GEO) database (https://www .ncbi.nlm.nih.gov/geo/). By using"DLBCL" or "diffuselar-geBcell lymphoma," "homo sapiens," and "expression profiling" as qualifiers, we searched and screened "microarray" and "RNA sequencing" data from the year 2006-2018 in the GEO database. The inclusion criteria were as follows: detectable CXCR4 expression level, available clinicopathological data, number of datasets larger than 30 patients, and tissue samples. The diagnosis and classification of DLBCL cases were based on the criteria of the 2008 World Health Organization classification of the lymphoid neoplasms. Two independent validation cohorts, one is tissue microarray including forty-four de novo DLBCL patients, another is formalin-fixed paraffin-embedded (FFPE) samples of 59 de novo DLBCL patients, were both obtained from Zhong-shan Hospital, Fudan University. This study was approved by the Institutional Review Board (IRB) of Zhongshan Hospital Affiliated to Fudan University (B2021-025R) and was conducted in accordance with the principles of the Declaration of Helsinki. Written informed consent was signed by each participant.

ICI Level and Tumor Immune Estimate Resource
(TIMER) Analysis. The CIBERSORT algorithm was performed to evaluate ICI in DLBCL tissues between CXCR4 high and low expression group. This algorithm transforms the normalized gene expression matrix into the composition of infiltrating immune cells. Analysis of ICI level between CXCR4 high and low expression groups was conducted by the "CIBERSORT R" package. Then, samples selected with a CIBERSORT P value <0.05 was filtered and reserved for the following analysis. We also used online tool TIMER2.0 (http://timer.cistrome.org/) to conduct Kaplan-Meier analysis to explore the infiltration of various immune cell subsets in patient survival from the TCGA database. Then, we used the TCGA database to explore how the CXCR4 expression was correlated with the expression of biomarkers that showed significance in the GSEA analysis.

Immunohistochemistry and Evaluation
Criterion. An independent cohort of 44 de novo DLBCL tissue microarrays from Zhongshan Hospital was stained with eight antibodies. IHC staining was performed on the Leica platform 2 Disease Markers (Bond Max) according to the antibody manual (Table S1) and manufacturer's instructions. IHC was assessed by two pathologists independently. They reexamined the immunostained slides to reach a consensus when there were inconsistencies. The percentage of positive CXCR4 stained cell area was scored into four grades: 0 for none; 1, <25%; 2, 25% to 50%; 3, 50% to 75%; and 4, >75%. The intensity was as follows: 0, negative; 1, weakly staining; 2, moderately staining; and 3, strongly staining. The final immunoreactivity scores were obtained by multiplying the percentage and the intensity score (0 to 12), and a final score lower than 2 was defined as CXCR4 negative, while a score ≥ 2 represented CXCR4 positivity.

2.5.
Real-Time Quantitative Polymerase Chain Reaction(qPCR). Five 5 μm unstained FFPE tissue sections were freshly cut for gene expression analysis. Total RNA was extracted using an FFPE Total RNA Isolation Kit (Canhelp Genomics Co., Ltd., Hangzhou, China) according to the protocols. For each specimen, reverse transcription was performed on isolated total RNA using the High-Capacity cDNA Reverse Transcription Kit with RNase Inhibitor (Applied Biosystems, Foster City, CA, USA). Next, the real-time PCR (RT-PCR) reaction was applied on the 7500 Real-Time PCR system (Applied Biosystems, Japan) for analyzing gene expression profiles. The RT-PCR program was initiated at 95°C for 10 minutes, followed by 40 cycles at 95°C for 15 seconds, and 60°C for 1 minute. The mRNA quantity of specific genes, calculated using the 2 -△△CT method, was normalized against GAPDH. All the measurements were performed in triplicate. The sequences of the CXCR4 primer pairs were as follows: CXCR4-F: 5′-GGTG GTCTATGTTGGCGTCT-3 ′ , CXCR4-R: 5 ′ -TGGAGT GTGACAGCTTGGAG-3 ′ .
2.6. In Vitro Proliferation Assay. The effects of CXCR4 inhibitor WZ811 (MedChemExpress) were investigated by a Cell Counting Kit-8 (CCK-8; Dojindo, Kumamoto, Japan) assay. Human DLBCL cell line OCI-Ly3 (ABC type) and SU-DHL4 (GCB type) were purchased from Fuheng Cell Center, Shanghai. For time-dependent assay, tumor cells were plated in a 96-well plate, and the cell density was depended on cultured time. For dose-dependent assay, tumor cells were plated at a density of 3000 cells with different concentrations of WZ811. All assays were performed in quadruplicate. At 0 h, 2 h, and 24 h, the cell proliferation assay was performed according to the manufacturer's instructions. Absorbance was measured at a wavelength of 450 nm using a microplate reader (Flexstation III ROM V2.1.28, USA).

Statistical
Analysis. Stata 14.2 (StataCorp, TX, USA) software was applied to conduct meta-analysis. Correlations between clinicopathological features and gene protein expression were determined using a X 2 test. Study heterogeneity was accessed by a X 2 test based on the Q statistical test. We also applied a sensitivity analysis to evaluate the stability of the pooled results. Univariate and multivariate Cox regression analyses were performed using R software (ver-sion 4.0.4). The effects of the CXCR4 expression on OS and PFS in DLBCL were analyzed using Kaplan-Meier analysis. The nonparametric Spearman r correlation coefficient was employed to assess the correlations between CXCR4 and other biomarkers. Cell viability test was analyzed using one-way ANOVA, and P < 0:05 was considered statistically significant.

Correlations between CXCR4 Expression and DLBCL
Subtypes Associated with Prognosis by Meta-Analysis. Information available for 2515 patients derived from 20 datasets were included in meta-analysis. The TCGA database (workflow type: RNASeq-FPKM) included 46 DLBCL cases with CXCR4 expression values and clinical characteristics (including gender, age, clinical stage, and survival). Nineteen datasets in the GEO database (experiment type: expression profiling by array) were included, with 15 datasets available for DLBCL subtypes and 10 datasets available for survival information (Table 1). No heterogeneity was found across the datasets (I 2 = 12:1%, P = 0:321).
To further clarify whether CXCR4 is an independent risk factor for DLBCL, we conducted univariate and multivariate analyses. Multivariate analysis result indicated that CXCR4 acts as an independent factor in the prediction of prognosis in DLBCL patients (Figure 1(d)). In the integrated analysis of survival ( Figure 1(e)), a trend of association between CXCR4 high expression and unfavorable OS was found in DLBCLABC type (P < 0:001) but not GCB type (P = 0:252), and no association between CXCR4 high expression and unfavorable PFS for both types was found (P > 0:05).

Gene Enrichment Analysis Related to CXCR4
Expression by GSEA and PPI Network. Subjects were categorized into low and high expression groups according to median CXCR4 expression level. Gene Ontology (GO) analysis showed that differential expression genes (DEGs) in subjects with higher CXCR4 expression were enriched in biological process (BP): regulation of RNA splicing, epigenetic modification, and mRNA export from the nucleus, etc.; cellular  (Table S2) exhibited that these genes were significantly enriched in signaling pathways such as the JNK-MAPK pathway, PI3K-AKT, mTOR signaling pathway, B-cell receptor (BCR) pathway, B-cell poor survival, IL-2 STAT5, IL-2 PI3k pathway, and NOTCH1 signal pathway. And genes had higher transcriptional activity of smad2/3/4 heterotrimer, adipogenesis, insulin resistance, VH rearrangement, doxorubicin resistance and alkylating resistance, CD5 target up, MYC and TFRC-up, and TP53 expression and degradation ( Figure 2(d)). However, no gene sets were significantly enriched in the CXCR4 low group at FDR < 25%.

ICI and Gene
Correlation Analysis by TIMER2.0. The CIBERSORT algorithm was used to calculate the proportion of immune cells in 816 samples from databases of GPL570 (GSE11318, GSE23501, GSE19246, GSE31312, GSE53786, GSE64555, GSE56313, GSE74266, and GSE93986). We found that high CXCR4 expression tumor tissues harbored a higher level of naive B-cells (P < 0:01), naive CD4+ T cells (P < 0:01), resting memory CD4+ T cells (P < 0:05), and regulatory T cells (Tregs) (P < 0:05), while low CXCR4 expression group harbored a higher level of CD8+ T cells (P < 0:05), macrophage 1 (P < 0:001), macrophage 2 (P < 0:01), and resting dendritic cells (P < 0:05) than high expression ones ( Figure 3(a)). The TIME2.0 website integrates a series of algorithms with TCGA database or the tumor expression data submitted by users themselves to estimate specific immune cell types. The run results are more reliable, and the results are displayed visually after a comprehensive evaluation through complex calculation methods. We used Kaplan-Meier plots to explore the relationship between ICI and DLBCL prognosis (CIBERSORT analytic module). We found that CD8 + T cell (P = 0:0355) infiltration was significantly correlated with DLBCL prognosis. There was no significant correlation between prognosis and macrophage 1 and macrophage 2 infiltration (Figure 3(b)). Based on the above GSEA and PPI results, we assessed the relationship between CXCR4 expression level and the expression of MYC, PDCD1(PD-1), CD274(PD-L1), Notch1, mTOR, CD5, FOXO1, and hnRNPA2 gene markers that showed significance in the GSEA analysis above using TIMER 2.0 (TCGA database). The results showed that the expression of CXCR4 was positively correlated with the expression of these genes (Figure 3(c)).

IHC Validation of CXCR4 Expression in DLBCL Patient
Samples Associated with Clinicopathological Features. Tissue   Table 2. No association was identified between the expression of CXCR4 and age, gender, Ann Arbor stage, and IPI score (P > 0:05). The protein expression level of CXCR4 was found to be positively associated with the cell of origin (COO) of the non-GCB subtype (P = 0:022), and with high levels of the mTOR expression (Figures 4(d)-4(f)). No correlation was found between the expression level of CXCR4 and CD8 + T lymphocyte infiltration, c-MYC, tp53, and PD-1 and PD-L1 expression.
Further, to verify the mRNA level of CXCR4 in DLBCL and its correlation to clinical characteristics, another independent cohort of 59 de novo DLBCL patients, all   7 Disease Markers diagnosed within 2 years and with complete treatment information, was applied for qPCR using their FFPE samples. Based on the qPCR results, high CXCR4 mRNA level was correlated to double-hit DLBCL (P = 0:0175) (Figure 4(g)). However, no correlation was found between the mRNA level of CXCR4 and double-expression DLBCL and Bruton's tyrosine kinase (BTK) inhibitor or lenalidomide treatment response.    (Tables 3   and 4). OCI-ly3 cells exhibited more sensitivity to WZ811 treatment than SU-DHL4 cells within the range <80 μM ( Figure 5). However, according to IC50 results, with the escalation of WZ811 concentration, SU-DHL4 cells exhibited

Discussion
Although the standard immunochemotherapy of R-CHOP has cured a majority of patients with DLBCL, treatment failure after R-CHOP still remains an intractable problem due to worse outcomes, which further underlines the importance of finding novel targets. Alizadeh et al. established the COO classification of DLBCL by using cDNA microarrays for gene expression profiling (GEP): the GCB and ABC subtype [17]. DLBCL patients with ABC subtype showed significantly worse prognosis compared with those with GCB subtype [18]. In our meta-analysis, CXCR4 was highly expressed in DLBCL patients with the ABC subtype and this was related to poor OS (Figures 1(a) and 1(b)). Our GSEA results revealed that the tumorigenicity of CXCR4 in DLBCL is related to the PI3K/AKT/mitogen-activated protein kinase (MAPK) pathways, mTOR, JNK-MAPK signal transducer and activator of transcription (STAT) pathways, Smad2/3/ 4, NOTCH1, and WNT signal pathways (Figure 2(d)) [19,20]. PI3K/AKT/mTOR pathway involved in crucial functions such as cellular proliferation, cell cycle regulation, and cell motility in DLBCL [21]. JAK-STAT signaling pathways were associated with pathogenesis of ABC-DLBCL [22]. DLBCL patients with the NOTCH1 mutations have worse PFS and OS [23]. These results are consistent with   [24][25][26][27][28]. Gene correlation analysis (TIMER) demonstrated that the CXCR4 expression is positively related to the degra-dation of tumor suppressor geneTP53, upregulation of the protooncogene C-MYC, and CD5expression (Figure 3(c)). CD5 + DLBCL has a poor prognosis and is classified as ABC-DLBCL with elusive genetic features [29]. This confirmed the role of CXCR4 in the ABC but not the GCB subtype of DLBCL, suggesting that CXCR4 exerts an adverse effect on the prognosis of the ABC subtype DLBCL through the above molecular signaling pathways. According to the study of Breinholt et al. [30], large B-cell lymphoma patient samples with significantly higher prevalence of MYC-BCL2double-hit showed significantly lower expression of PD-L1 and tumor-associated macrophages. Meanwhile, our qPCR analysis showed that the high CXCR4 expression correlated to molecular cytogenetic abnormality of double-hit DLBCL ( Figure 4(g)). The intrinsic relationship between CXCR4 and double-hit DLBCL deserves worth further research. However, we did not find correlations between the protein expression level of CXCR4 and CD8 + T lymphocyte infiltration, CD5, c-MYC, tp53, and PD-1 and PD-L1 expression. This discrepancy may due to the limited sample size of the IHC cohort.
Immune cells in the tumor microenvironment contribute to tumor progression and the antitumor responses. Their biological significance in DLBCL varies from different studies. CD8 + T lymphocytes and M1 macrophage play an important role in antitumor immune responses. Our results showed that the low CXCR4 expression group harbored a higher level of CD8 + T cells (P < 0:05) (Figure 3(a)) and M1 macrophage (P < 0:001) ( Figure S1), which may explain a better survival result of low CXCR4 expression group. One study assessed the safety, efficacy, and immunobiological effects of the CXCR4 antagonist (BL-8040, motixafortide) combined with a PD-1 inhibitor and chemotherapy in metastatic pancreatic ductal adenocarcinoma (PDAC) [10]. CXCR4 blockade promotes CD8 + effector T cell tumor infiltration and is synergistic with PD-1 inhibitors in PDAC mouse models. This study is in line with our study results. Besides, the high CXCR4 expression in DLBCL harbored a higher level of regulatory T cells (Tregs) (P < 0:05). It is known that Tregs are important modulators for the interaction between lymphoma cells and tumor microenvironment. Zhou [31] reported that a high infiltration of FOXP3/CTLA-4 doublepositive cells was significantly associated with poor prognosis. Chang [32] found that higher infiltration of intr-tumoral CD25 + FOXP3 + lymphocytes correlate with a favorable prognosis in patients with DLBCL. Previous studies on the biological significance of immune cells infiltration in DLBCL [33,34] are inconsistent. This may be due to the sample size and different materials (peripheral blood, tissue, and bone marrow) and methods (immunohistochemistry, flow cytometry, and sequencing). The effects of ICI in DLBCL on the development of DLBCL need to be verified by more fundamental experiments.
According to our GSEA results, CXCR4 coexpression genes are enriched during metabolism biological processes of adipogenesis and insulin resistance (Table S2). Recently, numerous studies have underlined the role of cholesterol metabolism in the tumorigenesis. Cholesterol metabolism regulates oncogenic signaling pathways, including the  12 Disease Markers PTEN/AKT/mTORC1 [35][36][37], MYC [38], Wnt [39], MAPK [40], and Hippo pathway [41]. Rinks et al. [42] reduced cellular cholesterol by blocking the high-density lipoprotein nanoparticles (HDL NP), which can also block the enhanced BCR signaling. In their study, they achieved both cellular cholesterol reduction and apoptosis in resistant ABC DLBCL cell lines by using ibrutinib and the spleen tyrosine kinase inhibitor together with HDL NP. Thus, some new cholesterol metabolites have recently become promising target drugs for cancer treatment [43][44][45]. In another study, metformin, a drug that can enhance insulin sensitivity, was identified to be associated with improved response rate and PFS in diabetic patients and is a possible therapeutic drug against DLBCL [46]. These results illustrate that reduction of cellular cholesterol or restoration of insulin sensitivity in lymphoma is a new method to induce tumor cell apoptosis. Our combined GSEA results indicate that cholesterol metabolism and insulin resistance may regulate complex signaling pathways and antitumor immunity and play important roles in high CXCR4 expression DLBCL. We hypothesize that CXCR4 antagonists may simultaneously inhibit MAPK, PI3k/AKT, and mTOR signaling pathway, inhibit cholesterol metabolism, and correct insulin resistance, which may exert a stronger inhibitory effect on the BCR signaling pathway.
Our GSEA results showed that genes from DLBCL patients with high CXCR4 expression are enriched in the forkhead box O-(FOXO-) mediated transcription of cell cycle genes. KOBAS analysis showed increased hub genes associated with the FOXO signaling pathway. FOXOs are tumor suppressors and can be used as key regulators of cell biology [47]. FOXO3a belongs to this family and functions as a trigger for apoptosis through the expression of genes necessary for cell death, and its inactivation is essential for the proliferation of immune cells. According to the study results of Kapoor et al. [48], Ibrutinib-resistant DLBCL (IB-R) cell lines showed downregulation of FOXO3a and PTEN levels in the nuclei and activation of AKT signaling. Therefore, they used a PI3K inhibitor and AKT inhibitor (idelalisib and MK2206) in IB-R cells, which resulted in increased ibrutinib-induced apoptosis. The drug selinexor (an exportin 1 inhibitor) together with ibrutinib can also increase the nuclear abundance of FOXO3a and PTEN in IB-R cells, leading to the restoration of ibrutinib-induced apoptosis. Our hub gene analysis revealed the coexpression of CXCR4 and hnRNPA2B1, which was confirmed by the (TIMER) gene correlation analysis. HnRNPA2B1 is involved in carcinogenesis by interacting with several proteins. One study found that it promotes the survival of tumor cells by preventing the activation of p53 [49]. In another study by Pan et al., CXCR4 was coimmunoprecipitated with cyclophilin A (CyPA). Stimulation of CXCR4 induces CyPA phosphorylation and nuclear translocation. CyPA formed a complex with heterogeneous nuclear ribonucleoprotein (hnRNP) A2, which underwent nuclear export in response to the activation of CXCR4 and was blocked by RNAi of CyPA. Overall, these findings demonstrate that the CXCR4 antagonist exerts both effects on BCR signaling and nuclear export, restores the function of FOXO3a ( Figure 6) and PTEN, reduces the nuclear export of hnRNPA2, and overcomes acquired resistance to BTK inhibitor in DLBCL patients.
According to our cell viability test, WZ811 suppressed the cell viability of both ABC-DLBCL (OCI-Ly3) and GCB-DLBCL (SU-DHL4), which confirmed the results of  13 Disease Markers our bioinformatics analysis ( Figure 5). Currently, there are a rising number of clinical trials using CXCR4 antagonists for treating different malignant diseases [50]. One phase I study (NCT00903968) enrolled relapse/refractory multiple myeloma (RRMM). The study drugs were bortezomib (1.3 mg/ m 2 ) and CXCR4 antagonist (the maximal-tolerated dose 0.32 mg/kg). Results of this clinical trial indicated that combination of these two drugs is safe and effective in RRMM. A phase Ib/II study also aimed to determine the safety and tolerability of the anti-CXCR4 antibody (ulocuplumab), alone and in combination with lenalidomide and dexamethasone, or with bortezomib and dexamethasone in patients with RRMM. This study showed that ulocuplumab is safe with acceptable adverse events and leads to a high response rate in combination with lenalidomide and dexamethasone in patients with RRMM. These results suggest that CXCR4 inhibitors are promising antimyeloma drugs. However, few clinical trials have been conducted on DLBCL; so, there is a potential for further exploration in clinical trials.
Certain limitations existed in our study. For example, there were a small number of verification samples in this study. Additionally, we simply divided the expression level of CXCR4 into high and low groups according to the median CXCR4 expression level and used negative and positive in the immunohistochemical stains. However, choosing the best immunohistochemical cutoff value may help us more sensitively and specifically identify high-risk DLBCL patients with the positive CXCR4 expression. Also, we used FFPE samples for qPCR due to lacking of fresh biopsy specimen, which may reduce the accuracy of experimental results. Limited variety of cell lines was also the shortcoming of our experiment. We hope to make up for the defects of these experiments in the future.

Conclusion
Overall, our study utilized a comprehensive analysis with experimental validation to imply the therapeutic significance of CXCR4 as a potential target in DLBCL, which may pave the way for precision therapy and translational medicine in the prognosis and treatment of DLBCL in the future.

Data Availability
The data used to support the findings of this study are included within the article.

Ethical Approval
This study was approved by the ethics committee of the Zhongshan Hospital, Fudan University (B2021-025R), and written informed consent was signed by each participant.

Consent
Written consent was obtained from all patients.

Conflicts of Interest
The authors declare no conflicting interests.