Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 14 November 2023
Sec. Cancer Immunity and Immunotherapy
This article is part of the Research Topic Diffuse Large B Cell Lymphoma: Updates on Diagnosis and Treatments View all 12 articles

EBF1, MYO6 and CALR expression levels predict therapeutic response in diffuse large B-cell lymphomas

Alice Turdo&#x;Alice Turdo1†Miriam Gaggianesi&#x;Miriam Gaggianesi2†Caterina D&#x;AccardoCaterina D’Accardo1Gaetana PorcelliGaetana Porcelli1Sebastiano Di BellaSebastiano Di Bella2Dario CricchioDario Cricchio2Irene PillitteriIrene Pillitteri1Rossana PorcasiRossana Porcasi1Melania Lo IaconoMelania Lo Iacono1Francesco VeronaFrancesco Verona1Chiara ModicaChiara Modica2Narges RoozafzayNarges Roozafzay1Ada Maria FlorenaAda Maria Florena1Giorgio Stassi*Giorgio Stassi2*Salvatrice Mancuso*&#x;Salvatrice Mancuso1*‡Matilde Todaro,&#x;Matilde Todaro1,3‡
  • 1Department of Health Promotion, Mother and Child Care, Internal Medicine and Medical Specialties (PROMISE), University of Palermo, Palermo, Italy
  • 2Department of Surgical, Oncological and Stomatological Sciences (DICHIRONS), University of Palermo, Palermo, Italy
  • 3A.O.U.P. “Paolo Giaccone”, University of Palermo, Palermo, Italy

Background: Diffuse large B-cell lymphoma (DLBCL) is a hematological malignancy representing one-third of non-Hodgkin’s lymphoma cases. Notwithstanding immunotherapy in combination with chemotherapy (R-CHOP) is an effective therapeutic approach for DLBCL, a subset of patients encounters treatment resistance, leading to low survival rates. Thus, there is an urgent need to identify predictive biomarkers for DLBCL including the elderly population, which represents the fastest-growing segment of the population in Western countries.

Methods: Gene expression profiles of n=414 DLBCL biopsies were retrieved from the public dataset GSE10846. Differentially expressed genes (DEGs) (fold change >1.4, p-value <0.05, n=387) have been clustered in responder and non-responder patient cohorts. An enrichment analysis has been performed on the top 30 up-regulated genes of responder and non-responder patients to identify the signatures involved in gene ontology (MSigDB). The more significantly up-regulated DEGs have been validated in our independent collection of formalin-fixed paraffin-embedded (FFPE) biopsy samples of elderly DLBCL patients, treated with R-CHOP as first-line therapy.

Results: From the analysis of two independent cohorts of DLBCL patients emerged a gene signature able to predict the response to R-CHOP therapy. In detail, expression levels of EBF1, MYO6, CALR are associated with a significant worse overall survival.

Conclusions: These results pave the way for a novel characterization of DLBCL biomarkers, aiding the stratification of responder versus non-responder patients.

1 Introduction

Diffuse large B-cell lymphoma (DLBCL) is the most common subtype of non-Hodgkin’s lymphoma and one of the highest mortality rates for all countries in the world within the elderly subjects (1). DLBCL is a heterogeneous disease at molecular and genetic level, characterized by a different biological behavior.

Although more than 50% of patients affected by DLBCL successfully respond to standard therapy, approximately 40% experience a relapse, making this neoplasia the leading cause of morbidity due to limited treatment options (2). Moreover, DLBCL commonly occurs in patients with comorbidities or in very elderly patients who warrants geriatric assessment prior treatment. Thus, a comprehensive examination of treatment efficacy versus the occurrence of side effects is required in order to predict tolerability, cardiotoxicity and the broad quality of life in frail patients (3).

Several studies have shown that the magnitude of clinical benefit rate in therapies for the treatment of DLBCL, which is mainly based on the use of immunotherapy in combination with chemotherapy (R-CHOP), reflects the molecular heterogeneity, including gene copy-number alterations and mutations (2, 4). Of note, in the last decades, the addition of rituximab to the standard CHOP therapy, significantly improved, by 10-15%, the overall survival of DLBCL patients (2).

Nonetheless comprehensive mechanisms underlying the refractoriness to R-CHOP have not been determined, several clinical parameters have been associated with treatment resistance and worse outcomes. The main prognostic model applied to DLBCL is based on the International Prognostic Index (IPI). The scoring system allows to stratify patients from low risk (0/1 score) to high risk (4/5 score) groups, depending on age, serum lactate dehydrogenase (LDH) levels, the eastern cooperative oncology group (ECOG) performance status, number of extranodal sites and Ann Arbor stage disease (5). Newly diagnosed DLBCL patients, treated with R-CHOP, are categorized according to the revised IPI, which facilitates the prognostic classification of patients (6).

Apart from the scoring system incorporating clinical parameters, advances in molecular characterization led to distinguish two different molecular subtypes of DLBCLs with a different biological behavior, the germinal center B-cell (GCB) lymphoma and the activated B-cell (ABC) lymphoma, this last associated with a poorer prognosis. These molecular subtypes of DLBCL are likewise arising from distinct cell of origin at diverse stages of lymphoid differentiation and specifically GCB from normal germinal-center B cells, while ABC from a post-germinal B cell (7).

Although several integrative approaches and models to detect patients at increased risk of relapse have been proposed, the identification of decisive driver biomarkers that can predict therapy response is still an unmet need. In the present study, in order to identify the gene expression profile of elderly (≥65-year-old) DLBCL patient’s responders and non-responders to the therapy with CHOP and R-CHOP, we benefited from a publicly available dataset (GSE10846) (7). Using a multiplexed gene expression analysis, furtherly validated by immunohistochemical evaluation, it has been identified a gene signature predictive of therapeutic response, in an independent cohort of DLBCL patients. From the molecular analysis emerged that expression of EBF1, MYO6 and CALR is able to select patients with distinct outcomes. Here, we provided biomarkers that could be of clinical interest to stratify elderly DLBCL patients, predicting the response to standard therapy, and develop novel therapeutic strategies based on the knowledge acquired, regarding validated molecular targets.

2 Materials and methods

2.1 Study populations

DLBCL tumor specimens and patients’ clinical data were obtained at the Hematology Unit, “P. Giaccone” Hospital of Palermo. Elderly patients (≥ 65-year-old) have been selected for the study and further classified in two cohorts of responder (n=13) and non-responder (n=6) to first-line R-CHOP therapy (validation cohort). A panel of hematologist and pathologist at the “P. Giaccone” Hospital followed the ESMO Clinical Practice Guidelines and Italian Society of Hematology guidelines for diagnosis, treatment and follow-up of DLBCL patients.

2.2 Statistical analysis

The training cohort (GSE10846, n=414) has been divided in two groups: patients treated with CHOP (n=181) and patients treated with R-CHOP (n=233). Patients’ cohort has been filtered by age (≥ 65-year-old) (n=188) and subsequently divided in responder (n=94) and non-responder (n=94) (7).Differential expressed genes (DEGs) (fold change >1.4, p-value <0.05, κ=387) have been clustered in the responder and non-responder groups and according to LDH levels. The activated B-cell (ABC) and germinal center B-cell (GCB) molecular subtypes have been reported as annotations.

Finally, the top 30 up-regulated genes of responder and non-responder patients have been used to perform enrichment analysis in order to identify the main signatures involved in gene ontology (MSigDB), considering molecular function, biological process and cellular component (p-value<10-7). The signatures associated with the first 30 upregulated genes were also computed by the QIAGEN Ingenuity Pathway Analysis software.

The association between features and patients’ overall survival was assessed by using Cox proportional- hazards model. Specifically, in the univariate analysis Cell of origin (COO), ECOG performance status, Extranodal Sites, IPI, LDH, sex and stage parameters were dichotomized according to (8, 9). The dichotomization of the identified signature was defined by using the median expression of each gene (MYO6, EBF1 and CALR). In the multivariate analysis, we combined our signature with each previously described feature.

To generate the Kaplan-Meier curves of overall survival by using the GSE10846 dataset, the initial population was filtered by age (≥ 65-year-old). “High” and “Low” groups were defined by using the median expression of each gene (MYO6, EBF1 and CALR) in the patient cohort.

All analyses were performed with R survival, survminer, and coxph libraries. Graphs were created by using the ggplot2 library.

2.3 RNA extraction and droplet digital PCR

Total RNA from FFPE tumor tissue specimens was isolated by using RNeasy FFPE Kit (Qiagen). 300 ng of total RNA was retro-transcribed with the high-capacity c-DNA reverse transcription kit (Applied Biosystem). In order to perform a four-gene multiplex assay, we used specific Droplet digital PCR (ddPCR- QX200 Droplet Reader) gene expression assays with FAM (n=2) and HEX (n=2) fluorophores. To optimize the multiplex reactions, from 100 to 300nM gene-specific primers have been used in combination with ddPCR supermix for probes (No-dUTP) and 25 ng of cDNA samples. Droplets were generated using the QX200 Droplet Generator (Bio-Rad) and dispensed into a 96 well-PCR plate. PCRs were performed in a ProFlex PCR System (Applied Biosystem) with the following protocol: 1x (95°C for 10 min), 50x (94°C for 30 sec, 56°C for 1 min), 1x (98°C for 10 min). After gene target amplification, samples were analyzed using QX200 Digital Droplet Reader (Bio-Rad). Gene expression analyses (copies/µl) were performed using QX Manager Software (1.2 Standard Edition) and normalized by using GAPDH.

2.4 Immunohistochemistry

FFPE lymphoma tissue specimens, stratified by age ≥ 65-year-old, were obtained from 11 responder patients and 4 non responder patients treated with R-CHOP.

Antigen retrieval was performed using the PT link system (Dako, Agilent Technologies, Santa Clara, CA, USA). Thereafter, sections were permeabilized with the 0.1% TRITON X-100 PBS for 10 min on ice, followed by 3% H202 and 10% human serum blocking incubation.

All slides were exposed overnight at 4°C to primary antibodies against Calreticulin (CAL-R) (ab22683; mouse IgG1; Abcam, Cambridge Science Park, UK), Myosin VI (MYO6) (MUD-19; mouse IgG1, Sigma-Aldrich), and EBF-1 (HPA061169; rabbit; Sigma-Aldrich). Staining was revealed using a biotin-streptavidin system (Dako LSAB2 System-HRP) and detected with the DAB substrate chromogen system (Dako). Nuclei were counterstained with Mayer’s Hematoxylin (Lillie’s Modification) Histological Staining Reagent (Dako). Staining was analyzed using an Olympus BX60 microscope. Immunohistochemical analysis were quantified with Image J.

3 Results

3.1 The analysis of a large cohort of elderly DLBCL patients revealed a gene signature associated to prognosis

Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous disease, causing high mortality in elderly patients. Despite the adverse effects, CHOP- and R-CHOP-based therapies result effective in the two-thirds of DLBCL patients, of which the rest portion experiences disease recurrence. Being DLBCL elderly patients the more susceptible to standard therapy side effects, in order to identify the genes predictive of therapy response, a gene expression analysis of 188 pretreatment biopsies of patients with an age ≥ 65 was retrieved from a publicly available dataset (GSE10846) (Figure 1A) (7).

FIGURE 1
www.frontiersin.org

Figure 1 Gene expression analysis of a large cohort of DLBCL patients reveal a gene signature associated to prognosis. (A) Workflow chart indicating the process to select DLBCL responder and non-responder gene signatures in the training cohort of the GSE10846 database. (B) Heatmap of differential expressed genes (DEGs) (fold change >1.4, p-value <0.05, κ=387) in responder versus non responders DLBCL patients. The LDH levels and the molecular subtype classification are shown. (C) Top ten up-regulated genes (log2) in responder (blue) and non-responder (orange) DLBCL patient cohort. (D) Enrichment analysis in gene ontology (MSigDB) in responder and non-responder DLBCL patient cohort. (E) Protein network analysis, generated with DEGs listed in Supplementary Table 1 The networks were generated through the use of QIAGEN IPA (QIAGEN Inc., https://digitalinsights.qiagen.com/IPA), in responder and non-responder DLBCL patient cohort.

Following unsupervised hierarchical clustering, the dichotomization of training cohort patients in responder and non-responder to CHOP and R-CHOP therapy allowed the identification of differentially expressed genes (DEGs) associated to a poor outcome (Figure 1B). The median age of the responder and non-responder cohorts of patients, to standard therapy, was comparable (74.04 versus 75.12-year-old) thus allowing the exclusion of age-related deaths. Of note, responder patients were mainly characterized by the GCB-like (62,5%) molecular subtype associated with a favorable outcome (Figure 1B). In accordance with well-established negative prognostic LDH parameter, responder DLBCL patients harbored lower LDH levels (1,0887 versus 1,8112) with respect to non-responder patients (Figure 1B; Supplementary Figures S1A, B). Analysis of gene expression profile, including 387 genes, of responder versus non-responder patients showed ten most differentially expressed genes (p-value ≤ 0.001) (Figure 1C; Supplementary Table 1). Specifically, ATAD3A, CALR, CWF19L1, GALT, MAGEA9, MAPKI8IP3, PSLNR, SEPTIN7P13, SLC19A1 and SLC38A5 resulted up-regulated in non-responder patients, while high expression levels of EBF1, EDNRA, CCDC18, CCDC186, FYB, MALAT1, MIS18BP1, MYO6, THRAP3 and TOP1, characterized responder DLBCL patients (Figure 1C).

The enrichment analysis of top 30 upregulated genes computed with Molecular Signatures Database (MSigDB), revealed six signatures associated with cell cycle, cell division and cytoskeleton organization, which are related to B cell malignant neoplasia (Figures 1D, E; Supplementary Table 2). Together these data provide evidence that these gene signatures may select responder from non-responder patients, identifying patients with a better life expectancy.

3.2 The expression levels of three genes predicted the response of DLBCL patients to R-CHOP therapy

In order to validate the expression levels of previously identified genes in dictating the dichotomization in life expectancy, we analyzed a cohort of naïve DLBCL patients, diagnosed and in follow-up at the Hematology/Oncology Unit of the “P. Giaccone” Hospital in Palermo, treated with R-CHOP as first -line therapy (validation cohort), selecting the frail cohort of DLBCL patients (≥ 65-year-old) (Table 1).

TABLE 1
www.frontiersin.org

Table 1 Clinical parameters of patients with Diffuse Large B Cell Lymphoma (DLBCL) treated with rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) therapy.

To overcome the low abundance and integrity of RNA content on formalin-fixed paraffin embedded (FFPE) biopsy samples, we adopted an implemented quantitative multiplex droplet digital PCR-based assay (Figure 1C; Figure 2A), from which emerged that two out ten DEGs, previously identified, EBF1 and MYO6 resulted up-regulated, at both mRNA and protein levels, in responder DLBCL patients (Figures 2B–D; Supplementary Figures 2A, B; Supplementary Table 3). Furthermore, the gene expression analysis of the ten upregulated genes, arisen from the non-responder included in the dataset (GSE10846), displayed an increasing trend of CALR mRNA levels, although not reaching statistical significance, which was paralleled by high protein expression levels (Figures 2B–D; Supplementary Table 4). These data indicate that EBF1, MYO6 and CALR could predict DLBCL patients’ response to R-CHOP therapy and aid in the stratification of responder versus non-responder patients.

FIGURE 2
www.frontiersin.org

Figure 2 EBF1 and MYO6 result highly expressed in responder DLBCL patients whilst CALR high expression characterize non responder DLBCL patients. (A) Workflow chart indicating the validation of the identified signature in our cohort of FFPE DLBCL samples. (B) Absolute mRNA levels (copies/µl) of EBF1, MYO6 and CALR in responder and non-responder DLBCL patients (n=19). Data are represented as mean  ±  SEM of three independent experiments. (C) Representative droplet digital PCR (ddPCR) scatter plots showing single, double, triple or quadruple positive droplets for EBF1, MYO6, CALR and CWF19L1 of FFPE samples of responder (pt#1) and non-responder DLBCL patient (pt#14). (D) Representative IHC analysis for EBF1, MYO6 and CALR of patients as in (C) Scale bar is 100 µm.

3.3 EBF1, MYO6 and CALR signature is associated with survival probability in DLBCL patients

To investigate the clinical significance of the identified signature, the magnitude of EBF1, MYO6, and CARL expression levels has been correlated to the survival data of ≥ 65-year-old DLBCL patients of the training cohort.

Transcriptome microarray analysis of a cohort of 154 DLBCL revealed a significant negative correlation between the signature EBF1low, MYO6low and CALRhigh expression and survival probability of patients, assuming a more pronounced significance than single gene expression (Figure 3A). Univariate analysis denoted that EBF1high, MYO6high and CALRlow signature expression is an independent positive prognostic factor of overall survival showing a higher statistical significance over several important clinical parameters, such as ECOG performance status, extranodal sites and stage (Supplementary Figure 3A). Importantly, EBF1high, MYO6high and CALRlow signature expression significantly increased the prognostic value of ECOG performance status, extranodal sites, IPI, LDH and stage (Supplementary Figure 3B).

FIGURE 3
www.frontiersin.org

Figure 3 Expression of EBF1, MYO6 and CALR correlate with prognosis in DLBCL. (A) Kaplan Meier overall survival (OS) curves of elderly DLBCL patients (GSE10846) stratified by high or low EBF1, MYO6 and CALR expression levels. (B) Functional protein association network of EBF1, MYO6 and CALR based on STRING database.

In line with the signature prognostic value, from the STRING network analysis emerged an implication of EBF1, MYO6, and CALR in the regulation of tumor progression together with the control of B lymphocyte gene transcription, intracellular vesicle transport and protein folding (Figure 3B; Supplementary Figure S3C). Being restricted the signature expression to DLBCLs as compared to non-tumoral lymphoid cells, a specific targeting of EBF1, MYO6 and CALR could be exploited for therapeutic intervention (Supplementary Figure S4A).

To investigate whether the EBF1, MYO6 and CALR differential expression could be an age-independent signature to predict R-CHOP response, we assessed their expression levels in a heterogeneous age group, which revealed a comparable expression level in adults as well as elderly DLBCL patients (Figure 4A). Of note, beside in the over 65-year-old, EBF1low, MYO6low and CALRhigh signature expression is able to dichotomize the response to R-CHOP treatment in under 65-year-old DLBCL patients (Figure 4B; Supplementary Figure S4B). Although EBF1high, MYO6high and CALRlow signature expression was an independent positive prognostic factor in under 65-year-old DLBCL patients, it did not increase the prognostic value of the selected clinical features in a multivariate analysis (Supplementary Figures 4C, D).

FIGURE 4
www.frontiersin.org

Figure 4 EBF1, MYO6 and CALR characterize DLBCL as compared to normal lymphoid tissues. (A) RNA seq expression data of EBF1, MYO6 and CALR in DLBCL patients at different age, retrieved from the TCGA database and analyzed by UALCAN. (B) Kaplan Meier overall survival (OS) curves of adult DLBCL patients (<65-year-old) (GSE10846) stratified by high or low EBF1, MYO6 and CALR expression levels. N/A stands for not available.

Our study unveils a robust predictive and prognostic signature able to determine the response to R-CHOP treatment in both under and over 65-year-old DLBCL patients. Notably, EBF1high, MYO6high and CALRlow signature expression ameliorates the prognostic power of the most important clinical parameters and, in particular, of the IPI clinical risk scoring system in the elderly patients.

4 Discussion

The clinical use of the so far identified gene signatures, mostly associated with tumor microenvironmental components, led to unsatisfied clinical outcomes of several DLBCL patients that remains poor (7, 10). Thus, it is becoming increasingly clear that the heterogeneity of patients, affected by DLBCL, has been underestimated, posing an urgent need to identify novel specific biomarkers for predicting the response to standard therapy.

Here, we found a new signature that significantly associates with the progression of disease that may be exploited for curative therapies in advanced DLBCLs. Droplet digital PCR and immunohistochemistry analysis of a cohort of ≥65-year-old naïve DLBCL patients revealed that EBF1, MYO6 and CALR expression levels stratify patients for the response to the standard R-CHOP therapy, regardless the IPI score, currently used in clinical settings. MYO6 is a motor protein, classified as unconventional myosin protein due to its reverse direction movement towards the actin filaments. MYO6 is involved in vesicular and macromolecules transport, cell migration and signaling (11). Albeit being implicated in prostate and breast cancer progression, here we uncover a novel role as favorable predictive biomarker in DLBCLs. CALR controls the protein folding by regulating the protein glucosylation-deglucosylation cycle and calcium homeostasis in the endoplasmic reticulum (12). CALR genetic alterations have been observed in several cancer types and correlated to a worse outcome. Specifically, CALR driver mutations has been described in myeloproliferative disorders but have not yet been associated to other hematological neoplasia (13, 14). Interestingly, the main driver mutations in CALR exon 9 change protein localization, from a cytoplasmic form to a membrane bound homodimers (15). Secretion of mutated CALR has also been observed in liquid biopsies samples as urine in bladder urothelial cancer patients (16). Our pioneer findings prospectively pose CALR as an extremely powerful biomarker in DLBCL patients and further analysis are necessary to clarify the mutational status of CALR also in DLBCL patients.

Being EBF1, together with E2A and Pax5, involved in the B-cell lineage commitment by regulating cell transcription, it has been also implicated in the development of B-cell-acute lymphoblastic leukemias (B-ALL) (17). EBF1 protein levels have been otherwise associated to better prognosis in colorectal cancer and cholangiocarcinoma and to a worse outcome in triple-negative breast cancer (1820). The herein reported findings indicate that the EBF1low/MYOlow/CALRhigh gene expression foresees the failure in therapeutic response of patients associated with a low-intermediate IPI risk, airing this identified signature as significant prognostic power. While further studies are needed to investigate MYO6 regulation of lymphoma cell dynamics and to design novel therapeutic approaches, a monoclonal antibody targeting the neoepitope generated by mutated CALR in myeloproliferative disorders has been generated. This approach could be prospectively applied to DLBCL patients harboring alterations in CALR (21). To date, specific agonists of EBF1 are not clinically available, however, evidence showed that the inhibition of EBF1 may be influenced by Notch and IL-7 signaling, whose modulation by already accessible compounds could indirectly interfere with EBF1 (17). Our findings provide evidence that EBF1, likely by preserving B cell entity, is required together with MYO and CALR for R-CHOP response and significantly associates with patient overall survival.

High-grade DLBCLs requires more intensive therapies mainly characterized by the addition of other chemotherapeutic treatments. Although the major limitation of this study relies in the lack of functional validation of the newly identified genes, EBF1, MYO6 and CALR could be considered targetable candidates to aid advanced DLBCL management.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

Ethical approval was not required for the studies involving humans because the study was performed on already published dataset and retrospective samples. The studies were conducted in accordance with the local legislation and institutional requirements. For the human samples used in this study, a written informed consent was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and the institutional requirements.

Author contributions

AT: Conceptualization, Funding acquisition, Supervision, Writing – original draft, Writing – review & editing. MG: Conceptualization, Funding acquisition, Supervision, Writing – original draft, Writing – review & editing. CD’A: Methodology, Writing – review & editing. GP: Methodology, Writing – review & editing. SB: Methodology, Writing – review & editing. DC: Methodology, Writing – review & editing. IP: Resources, Writing – review & editing. RP: Resources, Writing – review & editing. MLI: Methodology, Writing – review & editing. FV: Methodology, Writing – review & editing. CM: Methodology, Writing – review & editing. NR: Methodology, Writing – review & editing. AF: Resources, Writing – review & editing. GS: Funding acquisition, Writing – review & editing, Writing – original draft. SM: Funding acquisition, Writing – review & editing, Writing – original draft. MT: Conceptualization, Funding acquisition, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The research leading to these results has received funding from AIRC under IG 2018 - ID. 21492 project and from PSN2015, 6.2, CUP176J17000470001 project to MT; from Department PROMISE - “Angelo Ferrante” project to AT, IP, AMF and SM; from Department DICHIRONS - PJ_PROGRIC2022_MIS1A project to MG; from Italian Ministry for University - PNRR M4C2I1.3-HEAL ITALIA project PE 00000019 CUP B73C22001310006 to GS and MT.

Acknowledgments

AT is a research fellow funded by “Programma Operativo Complementare” 2014-2020 from Regione Siciliana. AT and MG have been supported by FONDO FINALIZZATO ALLA RICERCA DI ATENEO (FFR) – 2023.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1266265/full#supplementary-material

References

1. Hounsome L, Eyre TA, Ireland R, Hodson A, Walewska R, Ardeshna K, et al. Diffuse large B cell lymphoma (DLBCL) in patients older than 65 years: analysis of 3 year Real World data of practice patterns and outcomes in England. Br J Cancer (2022) 126(1):134–43. doi: 10.1038/s41416-021-01525-4

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Coiffier B, Lepage E, Briere J, Herbrecht R, Tilly H, Bouabdallah R, et al. CHOP chemotherapy plus rituximab compared with CHOP alone in elderly patients with diffuse large-B-cell lymphoma. N Engl J Med (2002) 346(4):235–42. doi: 10.1056/NEJMoa011795

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Handforth C, Clegg A, Young C, Simpkins S, Seymour MT, Selby PJ, et al. The prevalence and outcomes of frailty in older cancer patients: a systematic review. Ann Oncol (2015) 26(6):1091–101. doi: 10.1093/annonc/mdu540

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Schmitz R, Wright GW, Huang DW, Johnson CA, Phelan JD, Wang JQ, et al. Genetics and pathogenesis of diffuse large B-cell lymphoma. N Engl J Med (2018) 378(15):1396–407. doi: 10.1056/NEJMoa1801445

PubMed Abstract | CrossRef Full Text | Google Scholar

5. International Non-Hodgkin's Lymphoma Prognostic Factors P. A predictive model for aggressive non-Hodgkin's lymphoma. N Engl J Med (1993) 329(14):987–94. doi: 10.1056/NEJM199309303291402

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sehn LH, Berry B, Chhanabhai M, Fitzgerald C, Gill K, Hoskins P, et al. The revised International Prognostic Index (R-IPI) is a better predictor of outcome than the standard IPI for patients with diffuse large B-cell lymphoma treated with R-CHOP. Blood (2007) 109(5):1857–61. doi: 10.1182/blood-2006-08-038257

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Lenz G, Wright G, Dave SS, Xiao W, Powell J, Zhao H, et al. Stromal gene signatures in large-B-cell lymphomas. N Engl J Med (2008) 359(22):2313–23. doi: 10.1056/NEJMoa0802885

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Risueno A, Hagner PR, Towfic F, Fontanillo C, Djebbari A, Parker JS, et al. Leveraging gene expression subgroups to classify DLBCL patients and select for clinical benefit from a novel agent. Blood (2020) 135(13):1008–18. doi: 10.1182/blood.2019002414

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Ruppert AS, Dixon JG, Salles G, Wall A, Cunningham D, Poeschel V, et al. International prognostic indices in diffuse large B-cell lymphoma: a comparison of IPI, R-IPI, and NCCN-IPI. Blood (2020) 135(23):2041–8. doi: 10.1182/blood.2019002729

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Feng P, Li H, Pei J, Huang Y, Li G. Identification of a 14-gene prognostic signature for diffuse large B cell lymphoma (DLBCL). Front Genet (2021) 12:625414. doi: 10.3389/fgene.2021.625414

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Wang D, Zhu L, Liao M, Zeng T, Zhuo W, Yang S, et al. MYO6 knockdown inhibits the growth and induces the apoptosis of prostate cancer cells by decreasing the phosphorylation of ERK1/2 and PRAS40. Oncol Rep (2016) 36(3):1285–92. doi: 10.3892/or.2016.4910

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Fucikova J, Spisek R, Kroemer G, Galluzzi L. Calreticulin and cancer. Cell Res (2021) 31(1):5–16. doi: 10.1038/s41422-020-0383-9

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Mozes R, Gango A, Sulak A, Vida L, Reiniger L, Timar B, et al. Calreticulin mutation specific CAL2 immunohistochemistry accurately identifies rare calreticulin mutations in myeloproliferative neoplasms. Pathology (2019) 51(3):301–7. doi: 10.1016/j.pathol.2018.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Jutzi JS, Marneth AE, Jimenez-Santos MJ, Hem J, Guerra-Moreno A, Rolles B, et al. CALR-mutated cells are vulnerable to combined inhibition of the proteasome and the endoplasmic reticulum stress response. Leukemia (2023) 37(2):359–69. doi: 10.1038/s41375-022-01781-0

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Varricchio L, Falchi M, Dall'Ora M, De Benedittis C, Ruggeri A, Uversky VN, et al. Calreticulin: challenges posed by the intrinsically disordered nature of calreticulin to the study of its function. Front Cell Dev Biol (2017) 5:96. doi: 10.3389/fcell.2017.00096

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kageyama S, Isono T, Matsuda S, Ushio Y, Satomura S, Terai A, et al. Urinary calreticulin in the diagnosis of bladder urothelial carcinoma. Int J Urol (2009) 16(5):481–6. doi: 10.1111/j.1442-2042.2009.02287.x

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Lukin K, Fields S, Hartley J, Hagman J. Early B cell factor: Regulator of B lineage specification and commitment. Semin Immunol (2008) 20(4):221–7. doi: 10.1016/j.smim.2008.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Qiu Z, Guo W, Dong B, Wang Y, Deng P, Wang C, et al. EBF1 promotes triple-negative breast cancer progression by surveillance of the HIF1alpha pathway. Proc Natl Acad Sci U.S.A. (2022) 119(28):e2119518119. doi: 10.1073/pnas.2119518119

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Shen Z, Chen Y, Li L, Liu L, Peng M, Chen X, et al. Transcription Factor EBF1 Over-Expression Suppresses Tumor Growth in vivo and in vitro via Modulation of the PNO1/p53 Pathway in Colorectal Cancer. Front Oncol (2020) 10:1035. doi: 10.3389/fonc.2020.01035

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Armartmuntree N, Murata M, Techasen A, Yongvanit P, Loilome W, Namwat N, et al. Prolonged oxidative stress down-regulates Early B cell factor 1 with inhibition of its tumor suppressive function against cholangiocarcinoma genesis. Redox Biol (2018) 14:637–44. doi: 10.1016/j.redox.2017.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Tvorogov D, Thompson-Peach CAL, Fosselteder J, Dottore M, Stomski F, Onnesha SA, et al. Targeting human CALR-mutated MPN progenitors with a neoepitope-directed monoclonal antibody. EMBO Rep (2022) 23(4):e52904. doi: 10.15252/embr.202152904

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: diffuse large B-cell lymphoma, R-CHOP, therapy resistance, elderly patients, gene expression signature, biomarkers of response

Citation: Turdo A, Gaggianesi M, D’Accardo C, Porcelli G, Bella SD, Cricchio D, Pillitteri I, Porcasi R, Lo Iacono M, Verona F, Modica C, Roozafzay N, Florena AM, Stassi G, Mancuso S and Todaro M (2023) EBF1, MYO6 and CALR expression levels predict therapeutic response in diffuse large B-cell lymphomas. Front. Immunol. 14:1266265. doi: 10.3389/fimmu.2023.1266265

Received: 24 July 2023; Accepted: 17 October 2023;
Published: 14 November 2023.

Edited by:

Shimin Hu, University of Texas MD Anderson Cancer Center, United States

Reviewed by:

Fabrizio Vianello, University Hospital of Padua, Italy
Irma Olarte, General Hospital of Mexico, Mexico
Zarir Karanjawala, UC Davis Medical Center, United States

Copyright © 2023 Turdo, Gaggianesi, D’Accardo, Porcelli, Bella, Cricchio, Pillitteri, Porcasi, Lo Iacono, Verona, Modica, Roozafzay, Florena, Stassi, Mancuso and Todaro. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Giorgio Stassi, Z2lvcmdpby5zdGFzc2lAdW5pcGEuaXQ=; Salvatrice Mancuso, c2FsdmF0cmljZS5tYW5jdXNvQHVuaXBhLml0

These authors have contributed equally to this work and share first authorship

These authors have contributed equally to this work and share last authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.