Next Article in Journal
Compensatory Neural Recruitment for Error-Related Cerebral Activity in Patients with Moderate-To-Severe Obstructive Sleep Apnea
Previous Article in Journal
Could Neutrophil-to-Lymphocyte Ratio (NLR) Serve as a Potential Marker for Delirium Prediction in Patients with Acute Ischemic Stroke? A Prospective Observational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phenotypic Characterization of Diffuse Large B-Cell Lymphoma Cells and Prognostic Impact

by
Julie Devin
1,
Alboukadel Kassambara
1,2,
Angélique Bruyer
1,
Jérôme Moreaux
1,2,3,† and
Caroline Bret
1,2,3,*,†
1
CNRS UMR9002, Institute of Human Genetics, 34090 Montpellier, France
2
Department of Biological Hematology, St Eloi Hospital, 34295 Montpellier, France
3
University of Montpellier, Faculty of Medicine, 34090 Montpellier, France
*
Author to whom correspondence should be addressed.
Correspondence: The authors share last authorship.
J. Clin. Med. 2019, 8(7), 1074; https://doi.org/10.3390/jcm8071074
Submission received: 19 June 2019 / Revised: 12 July 2019 / Accepted: 15 July 2019 / Published: 22 July 2019
(This article belongs to the Section Hematology)

Abstract

:
Multiparameter flow cytometry (MFC) is a fast and cost-effective technique to evaluate the expression of many lymphoid markers in mature B-cell neoplasms, including diffuse large B cell lymphoma (DLBCL), which is the most frequent non-Hodgkin lymphoma. In this study, we first characterized by MFC the expression of 27 lymphoid markers in 16 DLBCL-derived cell lines to establish a robust algorithm for their authentication. Then, using the expression profile in DLBCL samples of the genes encoding B lymphoid markers that are routinely investigated by MFC, we built a gene expression-based risk score, based on the expression level of BCL2, BCL6, CD11c, and LAIR1, to predict the outcome of patients with DLBCL. This risk score allowed splitting patients in four risk groups, and was an independent predictor factor of overall survival when compared with the previously published prognostic factors. Lastly, to investigate the potential correlation between BCL2, BCL6, CD11c, and LAIR1 protein level and resistance to treatment, we investigated the response of the 16 DLBCL cell lines to cyclophosphamide, etoposide, doxorubicin, and gemcitabine. We found a correlation between BCL6 overexpression and resistance to etoposide. These results show the interest of MFC for the routine characterization of DLBCL cells and tumors samples for research and diagnostic/prognostic purposes.

Graphical Abstract

1. Introduction

Diffuse large B-cell lymphoma (DLBCL) represents 30% to 40% of all non-Hodgkin lymphoma (NHL) cases in adults. It displays heterogeneous clinical and biological characteristics. Although most patients with DLBCL achieve long-term remission, a third of them relapse after standard treatment with the combination of rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) [1]. In the last decade, gene expression profiling of DLBCL cells allowed the identification of two main molecular subtypes that are related to the cell of origin and are associated with different clinical outcomes: germinal-center B-cell–like (GCB) and activated B-cell–like (ABC). The GCB subgroup (50% of all DLBCL) is associated with a good prognosis, and tumor cells display a healthy germinal-center B cell gene expression profile. The ABC subgroup (30% of cases) has a poorer outcome, and tumor cells display a healthy peripheral blood activated B cell gene expression profile, particularly an NF-kB signature. The remaining 20% of DLBCL correspond to primary mediastinal B cell lymphomas, or are unclassified [2,3]. This DLBCL subdivision is progressively incorporated in the clinical practice, due to the availability of methods to assess the cell of origin status and the development of novel therapeutic agents that target the different DLBCL subtypes [4]. In the clinic, DLBCL are usually classified using immunohistochemistry (IHC) methods and the Hans’ algorithm as well as molecular approaches, including the NanoString technology [5]. According to the Hans’ algorithm, DLBCL are subdivided in GCB tumors that are CD10+ or CD10-/BCL6+/MUM1-/IRF4-, and in non-GCB tumors that are CD10-/MUM1+/IRF4+ (BCL6 can be positive or negative) [6]. The Lymph2X assay is one of the molecular methods that use the NanoString technology and a 20-gene panel for cell of origin determination. It can also be used with formalin-fixed and paraffin-embedded tissue samples [7].
Multiparameter flow cytometry (MFC) is a fast and cost-effective technique that is widely used for the diagnosis and follow-up of lymphoproliferative disorders. In the clinical routine, it allows characterizing DLBCL cells in dissociated lymph node samples and also in blood, bone marrow, and body fluids, in the case of extra-nodal dissemination. In this study, we evaluated the interest of the routine DLBCL characterization by MFC using 16 DLBCL-derived cell lines and gene profiling data from cohorts of patients with DLBCL. First, we assessed whether MFC and a panel of 27 lymphoid markers could be used to develop a rapid method to authenticate DLBCL-derived cell lines for research purposes. Then, based on our previous experience in building powerful risk scores for multiple myeloma [8], acute myeloid leukemia [9], and B cell lymphoma [10], we developed a gene expression-based risk score for DLBCL at diagnosis using four independent cohorts of newly-diagnosed patients and the list of genes encoding the B lymphoid markers tested in our MFC-based routine DLBCL cell analysis. We found that this risk score, based on the expression of BCL2, BCL6, LAIR1/CD305, and CD11c, allows the identification of patients with high-risk DLBCL. Lastly, we analyzed the correlation between BCL2, BCL6, LAIR1/CD305, and CD11c protein level and the response to conventional treatments in the 16 DLBCL-derived cell lines.

2. Materials and Methods

2.1. Cell Culture

The 16 DLBCL cell lines (DOHH2, HT, OCI-LY19, DB, OCI-LY1, SUDHL-4, SUDHL-5, SUDHL-10, NUDHL-1, OCI-LY7, WSU-DLCL2, SUDHL-6, NUDUL-1, U2932, OCI-LY3, and RI-1) were purchased from the American Type Culture Collection or from DSMZ (Leibniz-Institut DSMZ-Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH, Germany). They were grown in RPMI-1640 Glutamax medium (Gibco, Invitrogen, Cergy Pontoise, France), supplemented with 10% fetal bovine serum (FBS) (PAA laboratory GmbH, Pasching, Austria) (U2932, SUDHL-4, HT, DOHH2, SUDHL-10, RI-1, and WSU-DLCL2 cells), 20% FBS (OCI-LY3, DB, SUDHL-5, NUDHL-1, and SUDHL-6 cells), or 15% FBS (NUDUL-1 cells). OCI-LY1, and OCI-LY7 cells were cultured in IMDM Glutamax (Gibco, Invitrogen, Cergy Pontoise, France), supplemented with 20% FBS, and OCI-LY19 cells in MEM alpha modified Glutamax (Gibco, Invitrogen, Cergy Pontoise, France) with 20% FBS. Cultures were maintained at 37 °C in a humidified atmosphere with 5% CO2. Contamination by Mycoplasma species was regularly monitored.

2.2. Flow Cytometry Analysis

Expression of 27 normal and pathological B lymphoid markers (CD19, CD20, FMC7, CD22, CD23, Kappa, Lambda, CD10, CD5, CD38, CD27, CD39, CD43, CD62L, CD81, CD200, BCL2, BCL6, Ki67, IgM, LAIR1, CD123, CD11c, CD25, CD103, CD71, and CD180) was evaluated by labeling the 16 DLBCL cell lines with specific antibodies. Surface staining of the cell suspension was performed at the recommended volume per test in the dark at room temperature (see supplementary Table S1 for references and combinations). Among the markers that were studied, BCL2, BCL6, and Ki67 were evaluated by intra-cytoplasmic staining, using the fix and perm solutions of the kit GAS-002 (Nordic-Mubio©, Susteren, The Netherlands).
Flow cytometry data were acquired with a Canto II flow cytometer (Becton Dickinson©, Le-Pont-de-Claix, France). The instrument setup and calibration were in accordance with the EuroFlow standard operating procedures [11,12].
For each marker, the mean intensity of fluorescence was evaluated after the initial gating of the cells of interest based on the CD45/SSC (side scatter) plot to exclude cell debris. Singlets were included in the FSC-H/FSC-A (forward scatter) plot, and events were analyzed, according to the EuroFlow panel guidelines. With the standard antibody panel, DLBCL-derived cell lines with a B-cell phenotype were selected based on the CD19/CD3 plot (or CD20/CD3 plot, if CD19 was weakly expressed). The total B-cell population represented the “OR” Boolean gate between the kappa-positive or lambda-positive B cell population. Events present on the kappa/lambda diagonal were removed [13]. For the other lymphoid markers, the total B-cell population gated on CD19 (or CD20) represented the “AND” Boolean gate.

2.3. Cell Viability Assay

DLBCL cells were cultured in 96-well flat-bottom microtiter plates in the appropriate medium with FBS in the presence of increasing concentrations of mafosfamide (Santa-Cruz Biotechnology, Dallas, TX, USA), doxorubicin or etoposide (Selleckchem, Houston, TX, USA), or the BCL6 inhibitor 79-6 (Calbiochem, San Diego, CA, USA) for 4 days.
The number of viable cells was determined using the Cell Titer Glo Luminescent Cell Viability Assay from Promega (Promega Corp, Madison, WI, USA). This test is based on the quantification of cellular ATP, as a marker of metabolically active cells, using a Centro LB 960 luminometer (Berthold Technologies, Bad Wildbad, Germany). Data were expressed as the mean percentage of six replicates and were normalized to the untreated control.

2.4. Gene Expression Profiling and Statistical Analyses

Gene expression microarray data were from four independent cohorts of patients with DLBCL. The first cohort (i.e., training cohort) included 233 patients treated with R-CHOP (R-CHOP Lenz cohort). Results were validated using the CHOP Lenz cohort (n = 181 patients treated with CHOP) [14], the Melnick cohort (n = 69 patients) [15], and the FFPE R-CHOP cohort (n = 72 patients) [7]. The patients’ pre-treatment clinical characteristics were previously described in References [7,14,15]. Affymetrix gene expression data are publicly available via the online Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) under the accession numbers GSE10846 (Lenz cohorts), GSE23501 (Melnick cohort) and GSE53786 (FFPE R-CHOP cohort) (Affymetrix HG-U133 plus 2.0 microarrays). Data were normalized with a Microarray Suite version 5.0 (MAS 5.0) using Affymetrix default analysis settings, and global scaling as the normalization method. The trimmed mean target intensity of each array was arbitrarily set to 500. Expression profiling data (accession number GSE56315) of DLBCL samples (n = 73) and normal centrocyte and centro-blast samples (n = 5/each) from human tonsils [16] were also compared.
One or several probe sets were available for 24 of the 27 B-cell markers evaluated by MFC. In the presence of several probe sets for the same gene, the probe set with the highest coefficient of variation was retained (Supplementary Table S1). Overall survival (OS) differences between groups were calculated using the log-rank test. Multivariate analysis was performed using the Cox proportional hazards model and Genomicscape (http://genomicscape.com) [17]. Survival curves were plotted using the Kaplan-Meier method. All analyses were done with R.2.10.1 and Bioconductor, version 2.5.

2.5. Building the B-Cell Marker Risk Score

The prognostic significance of probe sets was evaluated in the training cohort using the Maxstat R function, which allows us to determine the optimal cut point for continuous variables [18,19], and Benjamini Hochberg multiple testing correction, and then was validated in the validation cohorts [18,19].
Multivariate COX analysis was used to identify genes including the independent prognostic marker (i.e., BCL6, BCL2, CD11c, and LAIR1). Then, patients with DLBCL were split in five groups, according to their BCL6, CD11c, BCL2, and LAIR1 expression profile [20]. Kaplan Meier analysis of these five different groups was performed and when two consecutive groups did not show any significant difference, they were merged.

2.6. Interaction Effect Quantification

For each combination, the percentage of expected growing cells in the case of effect independence was calculated, according to the Bliss equation [21,22].
fuc = fuAfuB
where fuc is the expected fraction of cells unaffected by the drug combination in the case of effect independence, and fuA and fuB are the fractions of cells unaffected by treatment A and B, respectively. The difference between the fraction of living cells in the cytotoxicity test and the fuc was considered as an estimation of the interaction effect, with negative values indicating synergism and positive values antagonism.

3. Results

3.1. A Barcode to Identify DLBCL-Derived Cell Lines by MFC

In our laboratory, we routinely use 16 DLBCL-derived cell lines that display a GCB profile (n = 12; DOHH2, HT, OCI-LY19, DB, OCI-LY1, SUDHL-4, SUDHL-5, SUDHL-10, NUDHL-1, OCI-LY7, WSU-DLCL2, and SUDHL-6) or an ABC profile (n = 4, NUDUL-1, U2932, OCI-LY3, and RI-1). To precisely determine their B-cell phenotype and check their identity using flow cytometry after each thawing, we evaluated the expression of 27 markers routinely investigated in mature B cell lymphoma (membrane markers: Kappa, Lambda, CD19, CD20, CD22, CD23, CD5, CD10, CD27, CD38, CD39, CD43, CD62L, CD81, CD200, FMC7, IgM, LAIR1, CD123, CD103, CD25, CD11c, CD71, and CD180, and intra-cytoplasmic markers: Ki67, BCL2, and BCL6) (Supplementary Table S2). We then generated an algorithm that is based on the mutually exclusive expression of the immunoglobulin kappa and lambda light chains, and on the selective expression (absent/present or low/bright expression) of specific markers within these two subgroups (Figure 1).
A first discrimination was done based on light chain expression (seven cell lines kappa-positive versus nine lambda-positive cell lines. In the kappa-positive cell lines (Figure 1A), the lymphoid marker FMC7 allowed discriminating three groups: bright, low, and near absent. Within the FMC7bright group, CD27 expression separated SUDHL-4 (CD27bright) and U2932 (CD27low) cells. Within the FMC7low group, CD62L discriminated between OCI-LY1 (CD62Lpositive) and OCI-LY7 (CD62Lnegative) cells, and, within the FMC7negative group, CD10 was not expressed by RI-1. Among the CD10positive cell lines, CD38 was bright in SUDHL-6 cells and weak in HT cells. In the kappa-positive group (Figure 1A), the lymphoid markers FMC7, CD27, CD62L, CD10, and CD38 allowed for discriminating the different cell lines.
In the lambda-positive group (Figure 1B), the different cell lines were separated using the lymphoid markers CD39, CD10, CD20, CD27, CD19, BCL2, and CD38.
We could separate the nine lambda-positive cell lines in two groups based on CD39 expression (Figure 1B). In the CD39positive group, CD10 staining allowed us to form two other sub-groups (CD10positive and CD10negative cells). Among the CD10positive cells, CD27 discriminated between NUDUL-1 (CD27positive) and DB (CD27negative) cells. Among the CD10negative cells, CD19 allowed us to separate the NUDHL-1 (CD19positive) and OCI-LY3 (CD19negative) cell lines. In CD39negative cells, the lymphoid marker CD20 identified three sub-groups: CD20bright, CD20low, and CD20negative (OCI-LY19 cells). In the CD20bright sub-group, BCL2 intracytoplasmic staining discriminated between WSU-DLCL2 (BCL2positive) and SUDHL-10 (BCL2negative) cells. Lastly, in the CD20low sub-group, DOHH2 (CD38positive) and SUDHL-5 (CD38negative) were separated by CD38 staining.
We validated these results by performing a principal component analysis (PCA) incorporating MFI of the B-cell markers in our algorithm (Supplementary Figure S1). Validation of the algorithm with an independent set of DLBCL cell lines will be of interest.

3.2. CD39 is A Useful Marker to Discriminate Between ABC and GCB DLBCL Tumor Samples

We then performed a significance analysis of microarrays (SAM) analysis of the probe sets that correspond to the 27 markers that we routinely use for the characterization of pathological mature B cells by MFC to compare their expression in ABC and GCB DLBCL samples with the aim to find additional markers to categorize tumor samples. For this analysis, we used the gene expression profiling data of ABC and GCB DLBCL samples (n = 167 and n = 183, respectively) from the CHOP and R-CHOP Lenz cohorts (non-classified tumor samples were excluded). Four of these 27 genes were differentially expressed in these two groups. CD10 and BCL6 were upregulated in GCB samples, whereas CD39 and BCL2 were upregulated in ABC samples (Figure 2A). These results are in agreement with previous studies showing that CD10, BCL6, and BCL2 are differentially expressed between ABC and GCB DLBCL [23] (CD10 and BCL6 are used in the Hans’ algorithm) and that BCL2 is overexpressed in the ABC subtype [24]. Furthermore, we confirmed CD39 gene upregulation in ABC DLBCL samples (Figure 2B, p < 0.0001 ABC vs. GCB group) and also in the ABC DLBCL-derived cell lines compared with the GCB DLBCL-derived cell lines (p < 0.05) (Figure 2C).

3.3. The B Cell Marker Risk Score

First, using the Maxstat R function [25], we investigated the prognostic value of the 27 B-cell markers in two cohorts of patients with DLBCL. High expression of BCL2, LAIR1, CD39, and CD103 was associated with a poor outcome in both the R-CHOP Lenz cohort (training cohort, n = 233) and in the CHOP Lenz cohort (n = 181 patients, one of the validation cohorts) (Figure 3A and Supplementary Figure S2A). Conversely, high expression of BCL6, CD10, CD11c, and CD81 was associated with a significant longer OS (Figure 3B and supplementary Figure S2B). We then investigated the expression of these eight markers in the ABC and GCB molecular subgroups (n = 167 and n = 183 samples, respectively, from the CHOP and R-CHOP Lenz cohorts, as above). BCL2, CD39, and CD103 were significantly overexpressed in ABC DLBCL, whereas BCL6, CD10, CD11c, and CD81 were overexpressed in GCB DLBCL (Supplementary Figure S3). These data were validated at the protein level only for CD39 (Figure 2C) and CD10 in the DLBCL-derived cell lines by MFC (Supplementary Figure S4).
Using a multivariate COX analysis, we found that only BCL6, BCL2, CD11c, and LAIR1 remained independent prognostic factors (Table 1). On the basis of these data, we used BCL6, BCL2, CD11c, and LAIR1 gene expression to create a risk score. We split patients in the training cohort in five groups, according to the tumor BCL6, CD11c, BCL2, and LAIR1 expression (high/low), and then performed Kaplan Meier analyses to determine the overall survival (OS) in the function of the expression level. This allowed us to identify four groups with significantly different OS (Figure 4A). Group 1 (very low risk, 31.3% of all patients) comprised patients with low BCL2 and LAIR1/CD305 expression and high BCL6 and CD11c expression. Group 2 (low risk, 40.8% of patients) comprised patients with high BCL6 or CD11c expression and low BCL2 and LAIR1/CD305 expression, or low BCL2 and CD11c expression and high BCL6 and CD11c expression. Group 3 (medium risk, 18.9% of patients) comprised patients with high BCL2 or LAIR1/CD305 expression and high BCL6 or CD11c expression, or low BCL2 and LAIR1/CD305 expression and high BCL6 and CD11c expression, or high BCL2 and LAIR1/CD305 expression and high BCL6 and CD11c expression. Group 4 (high risk, 9% of patients) included patients with high BCL2 and LAIR1/CD305 expression and low CD11c and BCL6 expression, or high BCL2 and LAIR1/CD305 expression and low CD11c or BCL2 expression. Group 1 and group 2 did not reach the median OS, whereas group 3 and group 4 had a median OS of 46 months and 9 months, respectively (p = 3.2E-11, Figure 4A). The prognostic value of the score was validated in the CHOP Lenz cohort (n = 181) and in the FFPE R-CHOP cohort (n = 72), and showed a trend (p = 0.06) in the Melnick cohort (n = 69) (Figure 4B–D).
We then asked whether the four genes included in the prognostic risk score (BCL2, BCL6, LAIR1, and CD11c) were differentially expressed between malignant and normal B cells. To this aim, we compared their expression in normal centrocytes (n = 5), normal centroblasts (n = 5), and DLBCL samples (n = 73) (GSE12195 dataset) [26]. BCL2, LAIR1, and CD11c were significantly overexpressed in DLBCL samples compared with normal centrocytes (p = 0.002, p = 0.006 and p = 0.00008, respectively) and centro-blasts (p = 0.0006, p = 0.002 and p = 0.0003, respectively) (Figure 5). Conversely, BCL6 was downregulated in DLBCL compared with normal centrocytes (p = 0.006) and centro-blasts (p = 0.002) (Figure 5). Previous studies showed that the BCL2 transcript and protein expression are low in germinal center cells [27], and that BCL6 is strongly expressed in centro-blasts and centrocytes compared with memory B cells and plasma-blasts in a normal population of tonsil B cells [28].
Moreover, CD10 and CD103 were significantly downregulated in DLBCL samples compared with normal centrocytes (p = 0.0004 and p = 0.013, respectively) and centroblasts (p = 0.0002 and p = 0.004, respectively). CD39 was significantly overexpressed in DLBCL samples compared with normal centrocytes (p = 0.0007) and centroblasts (p = 0.001). CD81 expression was comparable in the three groups (Supplementary Figure S5).
These results suggest that BCL2, LAIR1, and CD11c and BCL6 could be used to discriminate DLBCL cells from normal B cells during MFC analysis of dissociated tumor samples.
We then compared, in the R-CHOP Lenz cohort (n = 233 patients), the prognostic value of our risk score and of other previously identified prognostic factors, including the ABC and GCB molecular subgroups, age, International prognostic Index (IPI), and DNA repair score (Table 1A). In two-by-two comparisons, the risk score, age, IPI, and DNA repair score remained independent prognostic factors of OS (Table 1B). When all prognostic factors were tested together in multivariate COX analyses, only the risk score and the DNA repair score maintained their prognostic value (Table 1C).
Moreover, univariate COX analysis showed that the individual lymphoid markers BCL2, CD39, CD103, CD81, LAIR1, BCL6, CD5, and CD11c had a prognostic value in the Lenz R-CHOP cohort (Table 2A). When we tested all parameters together, BCL2, BCL6, CD5, LAIR1, and CD81 expression remained significant (Table 2B).
To verify the applicability of our risk score in routine flow cytometry analysis, we correlated the gene expression profile signal of the markers included in the risk score with their mean fluorescence intensity (MFI) in the 16 DLBCL cell lines, due to the absence of a high enough number of primary tumor cell samples for a statistically robust analysis, and found a significant correlation for all of them (r = 0.68 for BCL2, r = 0.59 for BCL6, r = 0.80 for LAIR1, and r = 0.84 for CD11c, p < 0.05, Supplementary Figure S6).

3.4. BCL6 Protein Expression is Correlated with DLBCL Cell Response to Etoposide

Lastly, we investigated the correlation between BCL2, LAIR1, BCL6, CD5, and CD11c expression in MCF and response to doxorubicin, mafosfamide, and etoposide in the 16 DLBCL cell lines. Only BCL6 expression was significantly and positively correlated with etoposide resistance (i.e., the half maximal inhibitory concentration, IC50) (r = 0.50, p < 0.05) (Figure 6). Since BCL6 is the most frequently involved oncogene in DLBCL, we then determined whether BCL6 inhibition could overcome etoposide drug resistance in DLBCL cell lines. To this aim, we used the low molecular weight compound 79-6 that binds to the corepressor binding groove of the BCL6 BTB domain. This BCL6 inhibitor inhibits tumor growth in mice xenografted with DLBCL cells without any toxicity [29]. A synergy combining 79-6 with low etoposide concentrations was identified (Supplementary Figure S7A). These data suggest that high BCL6 expression could be of interest to identify etoposide-resistant patients. BCL6 targeting could present a therapeutic interest to overcome etoposide resistance (Supplementary Figure S7A).

4. Discussion

DLBCL cell lines are widely used in laboratories to explore the cellular and molecular mechanisms involved in DLBCL pathogenesis. However, routine controls, such as cell authentication and contamination testing, are required to minimize behavioral changes due to cross-contamination with other cell lines, mycoplasma contamination, or phenotypic modifications caused by extensive culturing [30,31]. Therefore, a rapid and low-cost method for cell identity confirmation after each thawing is required, as currently done for mycoplasma detection [32]. A variety of methods are available for authentication testing. Among them, short tandem repeat profiling, which relies on a polymerase chain reaction-based assay to assess polymorphic tetra-nucleotide or penta-nucleotide repeats, has been recommended for human cell lines [30]. Nevertheless, this technique has some limitations, including possible genetic drift with continuous cell passaging [33]. In this case, we describe a simple phenotype-based method for the routine authentication of DLBCL cell lines using MFC. This test can be performed in about half day and does not require DNA extraction and analysis. Moreover, flow cytometry could be a useful method for the detection of very low cross-contamination, the risk of which can increase during extensive culture time and repeated freezing-thawing cycles. However, more attention should be provided to instrument and setting calibration for multiparameter MFC profiling. Validation of the algorithm using an independent set of DLBCL cell line will be important.
Our B-cell biomarkers could be used for different applications, such as discriminating ABC and GCB DLBCL, predicting the outcome of patients with DLBCL, and providing information related to the treatment response. For instance, we found that CD39 could be an additional marker to discriminate between ABC and GCB DLBCL tumors and cell lines, and that its upregulation in DLBCL is associated with a poor outcome (Figure 3). CD39 is encoded by the ENTPD1 gene, and is a typical cell surface enzyme with a catalytic site that faces the extracellular compartment. CD39 is an ectonucleotidase that catalyzes the hydrolysis of γ- and β-phosphate residues of nucleoside tri-phosphates and di-phosphates into mono-phosphate derivatives [34]. Cardoso et al. showed that CD39 expression analysis by flow cytometry is useful to discriminate B-cell lymphomas, particularly Burkitt lymphoma and ABC DLBCL [35].
In the clinical practice, DLBCL samples are classified according to their putative cell of origin by IHC and in the different subgroups using the Hans’s algorithm. By comparing the classifications obtained by IHC and gene expression profiling, Gutierrez-García et al. showed that the Hans’s algorithm misclassified 48% and 15% of GCB and ABC tumors, respectively. CD39 staining could be of interest to help better define the cell of origin in patients with DLBCL [36].
In current practice, patients with DLBCL are classified according the Ann Arbor classification system. IPI, age-adjusted IPI, and cell of origin are associated with a prognostic value [7,37,38,39]. More recently, gene-based risk scores have been proposed [10,40]. In the present study, we assessed the prognostic value of genes that encode lymphoid markers and built a risk score based on the expression of BCL2, BCL6, LAIR1, and CD11c. Our risk score allowed dividing patients with DLBCL in four risk groups (low risk, intermediate risk, and high risk). This approach is an independent predictor factor for OS when compared with the previously published prognostic factors.
These data also suggest that deregulated expression of genes from different signaling pathways is involved in DLBCL pathophysiology. For instance, BCL2 and BCL6 are two markers linked to germinal center B cells [41]. BCL2 is an anti-apoptosis factor located in the mitochondrial membrane and is important in normal B-cell development and differentiation. In DLBCL, the association of MYC and BCL2 rearrangements (double-hit lymphoma, DHL) or of MYC and BCL2 overexpression (double expresser lymphoma, DEL) has a significant prognostic value. These tumors show a more aggressive behavior and more frequent treatment failure after conventional therapy than non-DHL/non-DEL tumors [42,43]. DHL is a B-cell lymphoma that carries rearrangements between MYC and another oncogene, usually BCL2 and more rarely BCL6, BCL3, or CCND1. DHL and triple-hit lymphoma represent about 10% of all DLBCL [44]. The t(14;18)(q21;q32)/IGH-BCL2 translocation is observed in 20% to 30% of DLBCL, usually in the GCB group [45]. DEL, in which both MYC and BCL2 are overexpressed, represent 25% to 35% of all DLBCL, and are more often classified in the ABC/non-GCB rather than in the GCB group. DEL are associated with poor prognosis compared with non-DEL that overexpress either MYC or BCL2 [46,47]. BCL2 overexpression was associated with bad prognosis in our risk score. In the previously published BCL2 IHC score [48], patients with BCL2-overexpressing DLBCL benefitted from treatment with venetoclax. Venetoclax is a highly selective BCL2 inhibitor and is considered a breakthrough therapy for refractory or relapsed chronic lymphocytic leukemia [49]. It is currently in phase I/II development as monotherapy and/or combination therapy for NHL (including DLBCL, mantle cell lymphoma, and follicular lymphoma) and acute myeloid leukemia [50].
BCL6 is a DNA-binding protein and a transcriptional repressor involved in mediating growth suppression, and in blocking the terminal differentiation of B lymphocytes. BCL6 has also been implicated in class switch recombination and somatic hypermutation [41]. In GCB-DLBCL, gain-of-function mutations of EZH2 cooperate with BCL6 overexpression to inhibit centro-blast terminal differentiation [51]. In ABC-DLBCL, constitutive activation of BCL6 by chromosomal translocations and loss of function of PRDM1 block the terminal differentiation of B lymphocytes to plasma cells [52]. The most common chromosomic translocation in DLBCL (30% of cases) involves BCL6 on chromosome 3q27 and IGH on chromosome 14q32. However, BCL6 has many other potential translocation partners [45].
A recent meta-analysis found that, in the absence of a concomitant MYC translocation, the BCL2 and BCL6 translocation status does not provide any additional prognostic information for OS or progression free survival (PFS) [53].
Leukocyte-associated immunoglobulin-like receptor-1 (LAIR1), which is also known as CD305, is a member of the immunoglobulin (Ig) superfamily, and is a type І transmembrane glycoprotein. LAIR-1 includes a C2-type Ig-like domain and two immuno-receptor tyrosine-based inhibition motifs (ITIMs). It is broadly expressed by nearly all immune cells, such as T-cells, B-cells, and natural killer (NK) cells. The known LAIR1 ligands are extracellular matrix collagen and C1q, which is the first complement component. In vivo, LAIR1 inhibits B-cell receptor (BCR)-mediated signaling and controls kinase pathways involved in B-cell proliferation [54,55,56]. LAIR1 overexpression was associated with a bad prognosis in pediatric acute lymphoblastic leukemia (ALL) [57]. Unlike normal pre-B cells, patient-derived ALL cells express LAIR1 at high levels. Genetic studies revealed that LAIR1 critically modulates oncogenic signaling through recruitment of the inhibitory phosphatases PTPN6 and INPP5D. High LAIR1 expression blocks hyperactivation of the BCR pathway by inhibition of SYK, SRC, and the ERK kinase that promotes negative selection of auto-reactive B cells during B-cell development [58]. LAIR1 is also a biomarker of the “host response” cluster in the DLBCL classification described by Morin and colleagues [59]. The presence of tumor-infiltrating lymphocytes and increased expression of multiple components of the T-cell receptor (TCR) complex (TCR alpha and beta and CD3 subunits), CD2, and additional molecules associated with T/NK-cell activation and the complement cascade characterized the “host response” cluster [59]. The “host response” cluster in DLBCL displays an inflammatory/immune response and also high expression of monocyte/macrophage and dendritic cell markers. LAIR1 expression may be related to the tumor environment inflammatory response.
CD11c is an integrin alpha X chain protein (ITGAX gene). Integrins are heterodimeric integral membrane proteins composed of an alpha chain and a beta chain. CD11c combines with the beta 2 chain (ITGB2) to form a leukocyte-specific integrin, referred to as inactivated-C3b (iC3b) receptor 4 (CR4). CD11c is expressed in dendritic cells, macrophages, and monocytes and is a differentiation marker of inflammation. Moreover, CD11c is a diagnostic marker of hairy cell leukemia, and ITGAX expression is associated with aggressive prostate cancer [60]. The prognostic value of CD11c varies depending on the tumor type. For instance, CD11c expression in chronic lymphocytic leukemia cells is associated with a significantly higher incidence of Richter’s transformation [61]. Conversely, in gastric cancer, high CD11c expression is associated with a decreased risk of death and relapse, and may be regarded as an alternative indicator of a favorable prognosis [62]. CD11c is also a marker of dendritic cells. It has been reported that CD11c-positive dendritic cells promote the formation of germinal centers in the steady state and their response after immunization [63]. Furthermore, it has been suggested that, in follicular lymphoma, the presence of CD11c-positive dendritic cells in the tumor micro-environment is involved in tumor progression or recurrence by supporting regulatory T-cell infiltration [64].
Since GEP is not included in the current routine diagnostic work-up, immunohistochemistry or flow cytometry represent interesting substitutes. The validation of our risk score at the protein level using flow cytometry or immunohistochemistry could lead to a simple, fast, and cheaper prognostic tool for routine DLBCL diagnosis and prognosis. Interestingly, BCL2, BCL6, LAIR1, and CD11c expression by flow cytometry has already been reported in studies of normal and malignant B cells including clinical significance analyses [11,65,66,67,68]. Flow cytometry or immunohistochemistry analyses will allow us to analyze, simultaneously, the microenvironment that is recognized as being indirectly a reflection of the tumor clone and presenting a major impact in clinical course and response to treatment. Furthermore, it will be important to better understand the functional role of BCL2, BCL6, LAIR1, and CD11c in DLBCL pathogenesis and pathophysiology.
We found that high BCL6 expression, determined by flow cytometry, was associated with significantly decreased etoposide-induced toxicity. Kurosu et al. showed that, in B cell lymphoma, BCL6 overexpression inhibits the increase in ROS levels and apoptosis in response to etoposide and other chemotherapeutic reagents [69]. Moreover, we observed that the BCL6 inhibitor 79-6 could synergize with low concentrations of etoposide in DLBCL-derived cell lines. A recent study reported that BCL6 depletion induces BCL2 and BCL-XL upregulation in DLBCL cells. According to this observation, in most DLBCL cell lines, effective tumor cell killing requires the concomitant inhibition of BCL6 and BCL2 [70]. We validated the synergy between concomitant inhibition of BCL2 and BCL6 (Supplementary Figure S7B). Furthermore, combination of etoposide with IC20 of 79-6 and IC20 of Venetoclax demonstrated a significantly higher toxicity on DLBCL cells (Supplementary Figure S7D). These findings suggest that the combination of BCL2 and BCL6 inhibitors could be of interest to potentiate etoposide-induced toxicity in BCL6-overexpressing DLBCL cells.

5. Conclusions

We developed a flow cytometry-based algorithm for the authentication of DLBCL-derived cell lines used in research to limit the risk of cross-contamination or misidentification. Moreover, our risk score allows the identification of patients with high-risk DLBCL. Since it is based on the expression of four genes that encode markers routinely analyzed by MFC or IHC, it could represent a powerful tool for simple outcome prediction in DLBCL.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0383/8/7/1074/s1, Figure S1: Principal component analysis representing MFI of FCM biomarkers used for DLBCL cell lines discrimination in our identification algorithm. Figure S2: Survival curves of patients with DLBCL in the CHOP Lenz cohort divided in two groups according to the expression of the prognostic lymphoid markers. Figure S3: Expression of prognostic markers in the ABC and GCB DLBCL samples from the Lenz cohorts (n = 350). Figure S4: CD10, BCL2, BCL6, CD103, LAIR1, CD81 and CD11c expression in GCB and ABC DLBCL-derived cell lines. Figure S5: Comparison of CD10, CD81, CD39 and CD103 gene expression in DLBCL samples and normal centrocytes and centroblasts (GSE12195 dataset). Figure S6: Correlation between gene expression of the prognostic lymphoid markers (Affymetrix microarrays) and their protein expression in 16 DLBCL cell lines by MFC. Figure S7: Effect of the combination of etoposide (increasing concentrations) and the BCL6 inhibitor 79-6 in DLBCL-derived cell lines. Table S1: Antibodies references. Table S2: Mean fluorescence intensity of 27 lymphoid markers in 16 DLBCL-derived cell lines obtained by multiparameter flow cytometry.

Author Contributions

J.M. and C.B. designed and supervised the study. C.B., J.M., and J.D. performed and analyzed the data. J.M. and A.K. performed statistical analysis. C.B., J.M., and J.D. wrote the manuscript. J.M., A.B., and C.B. revised the manuscript for important intellectual content. J.M. and C.B. acquired funding.

Funding

This research was funded by Grants from French INCA (Institut National du Cancer) Institute, grant number: PLBIO18; 2018-160), ANR (TIE-Skip; 2017-CE15-0024-01), France Lymphome Espoir (Prix Jeune Chercheur 2015-2016) and SIRIC Montpellier Cancer. Grant INCa_Inserm_DGOS_12553 supported this work. A grant from the Fondation ARC pour la recherche sur le cancer (DOC20160604280 and DOC20180506917) supported JD.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DLBCL Diffuse large B-cell lymphoma
GEPGene expression profiling
COOCell of origin
GCBGerminal-center B-cell–like subgroup
ABCActivated B cell–like subtype
RRituximab
CHOPCyclophosphamide, doxorubicin, vincristine, and prednisone
EFSEvent-free survival
OSOverall survival
FCMFlow cytometry
FSCForward scatter
SSCSide scatter
NF-kBNuclear factor kB
ERKExtracellular signal-regulated kinases
IC50Inhibitory concentration 50
INPP5DSrc homology 2 domain containing inositol polyphosphate 5-phosphatase 1
IHCImmunohistochemistry
PCRPolymerase chain reaction
PTPN6Protein Tyrosine Phosphatase Non-Receptor Type 6
SAMSignificance analysis of microarray
STRShort tandem repeat
SYKSpleen tyrosine kinase

References

  1. Siegel, R.; Naishadham, D.; Jemal, A. Cancer statistics, 2012. CA. Cancer J. Clin. 2012, 62, 10–29. [Google Scholar] [CrossRef] [PubMed]
  2. Rosenwald, A.; Wright, G.; Chan, W.C.; Connors, J.M.; Campo, E.; Fisher, R.I.; Gascoyne, R.D.; Muller-Hermelink, H.K.; Smeland, E.B.; Giltnane, J.M.; et al. The use of molecular profiling to predict survival after chemotherapy for diffuse large-B-cell lymphoma. N. Engl. J. Med. 2002, 346, 1937–1947. [Google Scholar] [CrossRef] [PubMed]
  3. Alizadeh, A.A.; Eisen, M.B.; Davis, R.E.; Ma, C.; Lossos, I.S.; Rosenwald, A.; Boldrick, J.C.; Sabet, H.; Tran, T.; Yu, X.; et al. Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature 2000, 403, 503–511. [Google Scholar] [CrossRef] [PubMed]
  4. Roschewski, M.; Staudt, L.M.; Wilson, W.H. Diffuse large B-cell lymphoma-treatment approaches in the molecular era. Nat. Rev. Clin. Oncol. 2014, 11, 12–23. [Google Scholar] [CrossRef] [PubMed]
  5. Veldman-Jones, M.H.; Lai, Z.; Wappett, M.; Harbron, C.G.; Barrett, J.C.; Harrington, E.A.; Thress, K.S. Reproducible, Quantitative, and Flexible Molecular Subtyping of Clinical DLBCL Samples Using the NanoString nCounter System. Clin. Cancer Res. 2015, 21, 2367–2378. [Google Scholar] [CrossRef] [PubMed]
  6. Hans, C.P.; Weisenburger, D.D.; Greiner, T.C.; Gascoyne, R.D.; Delabie, J.; Ott, G.; Müller-Hermelink, H.K.; Campo, E.; Braziel, R.M.; Jaffe, E.S.; et al. Confirmation of the molecular classification of diffuse large B-cell lymphoma by immunohistochemistry using a tissue microarray. Blood 2004, 103, 275–282. [Google Scholar] [CrossRef] [PubMed]
  7. Scott, D.W.; Wright, G.W.; Williams, P.M.; Lih, C.-J.; Walsh, W.; Jaffe, E.S.; Rosenwald, A.; Campo, E.; Chan, W.C.; Connors, J.M.; et al. Determining cell-of-origin subtypes of diffuse large B-cell lymphoma using gene expression in formalin-fixed paraffin embedded tissue. Blood 2014, 123, 1214–1217. [Google Scholar] [CrossRef] [PubMed]
  8. Moreaux, J.; Klein, B.; Bataille, R.; Descamps, G.; Maïga, S.; Hose, D.; Goldschmidt, H.; Jauch, A.; Rème, T.; Jourdan, M.; et al. A high-risk signature for patients with multiple myeloma established from the molecular classification of human myeloma cell lines. Haematologica 2011, 96, 574–582. [Google Scholar] [CrossRef] [PubMed]
  9. Bou Samra, E.; Klein, B.; Commes, T.; Moreaux, J. Development of gene expression-based risk score in cytogenetically normal acute myeloid leukemia patients. Oncotarget 2012, 3, 824–832. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Bret, C.; Klein, B.; Moreaux, J. Gene expression-based risk score in diffuse large B-cell lymphoma. Oncotarget 2012, 3, 1700–1710. [Google Scholar] [CrossRef] [Green Version]
  11. van Dongen, J.J.M.; Lhermitte, L.; Böttcher, S.; Almeida, J.; van der Velden, V.H.J.; Flores-Montero, J.; Rawstron, A.; Asnafi, V.; Lécrevisse, Q.; Lucio, P.; et al. EuroFlow antibody panels for standardized n-dimensional flow cytometric immunophenotyping of normal, reactive and malignant leukocytes. Leukemia 2012, 26, 1908–1975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. van Dongen, J.J.M.; Orfao, A. EuroFlow: Resetting leukemia and lymphoma immunophenotyping. Basis for companion diagnostics and personalized medicine. Leukemia 2012, 26, 1899–1907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Alaterre, E.; Raimbault, S.; Garcia, J.-M.; Rème, T.; Requirand, G.; Klein, B.; Moreaux, J. Automated and simplified identification of normal and abnormal plasma cells in Multiple Myeloma by flow cytometry. Cytom. Part B Clin. Cytom. 2017, 94, 484–492. [Google Scholar] [CrossRef] [PubMed]
  14. Lenz, G.; Wright, G.; Dave, S.S.; Xiao, W.; Powell, J.; Zhao, H.; Xu, W.; Tan, B.; Goldschmidt, N.; Iqbal, J.; et al. Stromal gene signatures in large-B-cell lymphomas. N. Engl. J. Med. 2008, 359, 2313–2323. [Google Scholar] [CrossRef] [PubMed]
  15. Shaknovich, R.; Geng, H.; Johnson, N.A.; Tsikitas, L.; Cerchietti, L.; Greally, J.M.; Gascoyne, R.D.; Elemento, O.; Melnick, A. DNA methylation signatures define molecular subtypes of diffuse large B-cell lymphoma. Blood 2010, 116, e81–e89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Dybkær, K.; Bøgsted, M.; Falgreen, S.; Bødker, J.S.; Kjeldsen, M.K.; Schmitz, A.; Bilgrau, A.E.; Xu-Monette, Z.Y.; Li, L.; Bergkvist, K.S.; et al. Diffuse large B-cell lymphoma classification system that associates normal B-cell subset phenotypes with prognosis. J. Clin. Oncol. 2015, 33, 1379–1388. [Google Scholar] [CrossRef]
  17. Kassambara, A.; Rème, T.; Jourdan, M.; Fest, T.; Hose, D.; Tarte, K.; Klein, B. GenomicScape: An easy-to-use web tool for gene expression data analysis. Application to investigate the molecular events in the differentiation of B cells into plasma cells. PLoS Comput. Biol. 2015, 11, e1004077. [Google Scholar] [CrossRef]
  18. Kassambara, A.; Hose, D.; Moreaux, J.; Walker, B.A.; Protopopov, A.; Reme, T.; Pellestor, F.; Pantesco, V.; Jauch, A.; Morgan, G.; et al. Genes with a spike expression are clustered in chromosome (sub)bands and spike (sub)bands have a powerful prognostic value in patients with multiple myeloma. Haematologica 2012, 97, 622–630. [Google Scholar] [CrossRef]
  19. Hothorn, T.; Lausen, B. On the exact distribution of maximally selected rank statistics. Stat. Med. 2012, 3178–3191. [Google Scholar] [CrossRef]
  20. Alaterre, E.; Raimbault, S.; Goldschmidt, H.; Bouhya, S.; Requirand, G.; Robert, N.; Boireau, S.; Seckinger, A.; Hose, D.; Klein, B.; et al. CD24, CD27, CD36 and CD302 gene expression for outcome prediction in patients with multiple myeloma. Oncotarget 2017, 8, 98931–98944. [Google Scholar] [CrossRef]
  21. Greco, W.R.; Bravo, G.; Parsons, J.C. The search for synergy: A critical review from a response surface perspective. Pharmacol. Rev. 1995, 47, 331–385. [Google Scholar] [PubMed]
  22. Combes, E.; Andrade, A.F.; Tosi, D.; Michaud, H.-A.; Coquel, F.; Garambois, V.; Desigaud, D.; Jarlier, M.; Coquelle, A.; Pasero, P.; et al. Inhibition of Ataxia-Telangiectasia Mutated and RAD3-related (ATR) overcomes oxaliplatin resistance and promotes anti-tumor immunity in colorectal cancer. Cancer Res. 2019, 79, 2933–2946. [Google Scholar] [CrossRef] [PubMed]
  23. Shaffer, A.L.; Young, R.M.; Staudt, L.M. Pathogenesis of human B cell lymphomas. Annu. Rev. Immunol. 2012, 30, 565–610. [Google Scholar] [CrossRef] [PubMed]
  24. Deeb, S.J.; D’Souza, R.C.J.; Cox, J.; Schmidt-Supprian, M.; Mann, M. Super-SILAC allows classification of diffuse large B-cell lymphoma subtypes by their protein expression profiles. Mol. Cell. Proteom. MCP 2012, 11, 77–89. [Google Scholar] [CrossRef] [PubMed]
  25. Lausen, B.; Schumacher, M. Maximally Selected Rank Statistics. Biometrics 1992, 48, 73–85. [Google Scholar] [CrossRef]
  26. Pasqualucci, L.; Compagno, M.; Lim, W.K.; Grunn, A.; Nandula, S.V.; Scandurra, M.; Bertoni, F.; Ponzoni, M.; Califano, A.; Bhagat, G.; et al. Mutations in Multiple Genes Cause Deregulation of the NFkB Pathway in Diffuse Large B-Cell Lymphoma. Blood 2008, 112, 801. [Google Scholar]
  27. Shen, Y.; Iqbal, J.; Huang, J.Z.; Zhou, G.; Chan, W.C. BCL2 protein expression parallels its mRNA level in normal and malignant B cells. Blood 2004, 104, 2936–2939. [Google Scholar] [CrossRef]
  28. Kjeldsen, M.K.; Perez-Andres, M.; Schmitz, A.; Johansen, P.; Boegsted, M.; Nyegaard, M.; Gaihede, M.; Bukh, A.; Johnsen, H.E.; Orfao, A.; et al. Multiparametric flow cytometry for identification and fluorescence activated cell sorting of five distinct B-cell subpopulations in normal tonsil tissue. Am. J. Clin. Pathol. 2011, 136, 960–969. [Google Scholar] [CrossRef]
  29. Cerchietti, L.C.; Ghetu, A.F.; Zhu, X.; Da Silva, G.F.; Zhong, S.; Matthews, M.; Bunting, K.L.; Polo, J.M.; Farès, C.; Arrowsmith, C.H.; et al. A small-molecule inhibitor of BCL6 kills DLBCL cells in vitro and in vivo. Cancer Cell 2010, 17, 400–411. [Google Scholar] [CrossRef]
  30. Masters, J.R. Human cancer cell lines: Fact and fantasy. Nat. Rev. Mol. Cell Biol. 2000, 1, 233–236. [Google Scholar] [CrossRef]
  31. Reid, Y.A. Characterization and authentication of cancer cell lines: An overview. Methods Mol. Biol. 2011, 731, 35–43. [Google Scholar] [PubMed]
  32. Maïga, S.; Brosseau, C.; Descamps, G.; Dousset, C.; Gomez-Bougie, P.; Chiron, D.; Ménoret, E.; Kervoelen, C.; Vié, H.; Cesbron, A.; et al. A simple flow cytometry-based barcode for routine authentication of multiple myeloma and mantle cell lymphoma cell lines. Cytom. Part A 2015, 87, 285–288. [Google Scholar] [CrossRef] [PubMed]
  33. Parson, W.; Kirchebner, R.; Mühlmann, R.; Renner, K.; Kofler, A.; Schmidt, S.; Kofler, R. Cancer cell line identification by short tandem repeat profiling: Power and limitations. FASEB J. 2005, 19, 434–436. [Google Scholar] [CrossRef] [PubMed]
  34. Robson, S.C.; Sévigny, J.; Zimmermann, H. The E-NTPDase family of ectonucleotidases: Structure function relationships and pathophysiological significance. Purinergic Signal. 2006, 2, 409–430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Cardoso, C.C.; Auat, M.; Santos-Pirath, I.M.; Rudolf-Oliveira, R.C.M.; da Silva, J.P.; Lange, B.G.; Siegel, D.; de Moraes, A.C.R.; Del Moral, J.A.G.; Santos-Silva, M.C. The importance of CD39, CD43, CD81, and CD95 expression for differentiating B cell lymphoma by flow cytometry. Cytom. Part B Clin. Cytom. 2017, 94, 451–458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gutiérrez-García, G.; Cardesa-Salzmann, T.; Climent, F.; González-Barca, E.; Mercadal, S.; Mate, J.L.; Sancho, J.M.; Arenillas, L.; Serrano, S.; Escoda, L.; et al. Gene-expression profiling and not immunophenotypic algorithms predicts prognosis in patients with diffuse large B-cell lymphoma treated with immunochemotherapy. Blood 2011, 117, 4836–4843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Tilly, H.; da Silva, M.G.; Vitolo, U.; Jack, A.; Meignan, M.; Lopez-Guillermo, A.; Walewski, J.; André, M.; Johnson, P.W.; Pfreundschuh, M.; et al. Diffuse large B-cell lymphoma (DLBCL): ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 2015, 26, v116–v125. [Google Scholar] [CrossRef]
  38. Wright, G.; Tan, B.; Rosenwald, A.; Hurt, E.H.; Wiestner, A.; Staudt, L.M. A gene expression-based method to diagnose clinically distinct subgroups of diffuse large B cell lymphoma. Proc. Natl. Acad. Sci. USA 2003, 100, 9991–9996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Lenz, G.; Wright, G.W.; Emre, N.C.T.; Kohlhammer, H.; Dave, S.S.; Davis, R.E.; Carty, S.; Lam, L.T.; Shaffer, A.L.; Xiao, W.; et al. Molecular subtypes of diffuse large B-cell lymphoma arise by distinct genetic pathways. Proc. Natl. Acad. Sci. USA 2008, 105, 13520–13525. [Google Scholar] [CrossRef] [Green Version]
  40. Bret, C.; Klein, B.; Cartron, G.; Schved, J.-F.; Constantinou, A.; Pasero, P.; Moreaux, J. DNA repair in diffuse large B-cell lymphoma: A molecular portrait. Br. J. Haematol. 2015, 169, 296–299. [Google Scholar] [CrossRef]
  41. Schmidlin, H.; Diehl, S.A.; Blom, B. New insights into the regulation of human B-cell differentiation. Trends Immunol. 2009, 30, 277–285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Hsi, E.D.; Yegappan, S. Lymphoma immunophenotyping: A new era in paraffin-section immunohistochemistry. Adv. Anat. Pathol. 2001, 8, 218–239. [Google Scholar] [CrossRef] [PubMed]
  43. Bilalovic, N.; Blystad, A.K.; Golouh, R.; Nesland, J.M.; Selak, I.; Trinh, D.; Torlakovic, E. Expression of bcl-6 and CD10 protein is associated with longer overall survival and time to treatment failure in follicular lymphoma. Am. J. Clin. Pathol. 2004, 121, 34–42. [Google Scholar] [CrossRef] [PubMed]
  44. Aukema, S.M.; Siebert, R.; Schuuring, E.; van Imhoff, G.W.; Kluin-Nelemans, H.C.; Boerma, E.-J.; Kluin, P.M. Double-hit B-cell lymphomas. Blood 2011, 117, 2319–2331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Li, S.; Young, K.H.; Medeiros, L.J. Diffuse large B-cell lymphoma. Pathology 2018, 50, 74–87. [Google Scholar] [CrossRef]
  46. Johnson, N.A.; Slack, G.W.; Savage, K.J.; Connors, J.M.; Ben-Neriah, S.; Rogic, S.; Scott, D.W.; Tan, K.L.; Steidl, C.; Sehn, L.H.; et al. Concurrent expression of MYC and BCL2 in diffuse large B-cell lymphoma treated with rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisone. J. Clin. Oncol. 2012, 30, 3452–3459. [Google Scholar] [CrossRef]
  47. Hu, S.; Xu-Monette, Z.Y.; Tzankov, A.; Green, T.; Wu, L.; Balasubramanyam, A.; Liu, W.; Visco, C.; Li, Y.; Miranda, R.N.; et al. MYC/BCL2 protein coexpression contributes to the inferior survival of activated B-cell subtype of diffuse large B-cell lymphoma and demonstrates high-risk gene expression signatures: A report from The International DLBCL Rituximab-CHOP Consortium Program. Blood 2013, 121, 4021–4031. [Google Scholar] [CrossRef]
  48. Tsuyama, N.; Sakata, S.; Baba, S.; Mishima, Y.; Nishimura, N.; Ueda, K.; Yokoyama, M.; Terui, Y.; Hatake, K.; Kitagawa, M.; et al. BCL2 expression in DLBCL: Reappraisal of immunohistochemistry with new criteria for therapeutic biomarker evaluation. Blood 2017, 130, 489–500. [Google Scholar] [CrossRef]
  49. Souers, A.J.; Leverson, J.D.; Boghaert, E.R.; Ackler, S.L.; Catron, N.D.; Chen, J.; Dayton, B.D.; Ding, H.; Enschede, S.H.; Fairbrother, W.J.; et al. ABT-199, a potent and selective BCL-2 inhibitor, achieves antitumor activity while sparing platelets. Nat. Med. 2013, 19, 202–208. [Google Scholar] [CrossRef]
  50. Deeks, E.D. Venetoclax: First Global Approval. Drugs 2016, 76, 979–987. [Google Scholar] [CrossRef]
  51. Morin, R.D.; Johnson, N.A.; Severson, T.M.; Mungall, A.J.; An, J.; Goya, R.; Paul, J.E.; Boyle, M.; Woolcock, B.W.; Kuchenbauer, F.; et al. Somatic mutations altering EZH2 (Tyr641) in follicular and diffuse large B-cell lymphomas of germinal-center origin. Nat. Genet. 2010, 42, 181–185. [Google Scholar] [CrossRef] [PubMed]
  52. Mandelbaum, J.; Bhagat, G.; Tang, H.; Mo, T.; Brahmachary, M.; Shen, Q.; Chadburn, A.; Rajewsky, K.; Tarakhovsky, A.; Pasqualucci, L.; et al. BLIMP1 is a tumor suppressor gene frequently disrupted in activated B cell-like diffuse large B cell lymphoma. Cancer Cell 2010, 18, 568–579. [Google Scholar] [CrossRef] [PubMed]
  53. Schmidt-Hansen, M.; Berendse, S.; Marafioti, T.; McNamara, C. Does cell-of-origin or MYC, BCL2 or BCL6 translocation status provide prognostic information beyond the International Prognostic Index score in patients with diffuse large B-cell lymphoma treated with rituximab and chemotherapy? A systematic review. Leuk. Lymphoma 2017, 58, 2403–2418. [Google Scholar] [CrossRef]
  54. Meyaard, L. The inhibitory collagen receptor LAIR-1 (CD305). J. Leukoc. Biol. 2008, 83, 799–803. [Google Scholar] [CrossRef] [PubMed]
  55. Son, M.; Santiago-Schwarz, F.; Al-Abed, Y.; Diamond, B. C1q limits dendritic cell differentiation and activation by engaging LAIR-1. Proc. Natl. Acad. Sci. USA 2012, 109, E3160–E3167. [Google Scholar] [CrossRef] [PubMed]
  56. van der Vuurst de Vries, A.R.; Clevers, H.; Logtenberg, T.; Meyaard, L. Leukocyte-associated immunoglobulin-like receptor-1 (LAIR-1) is differentially expressed during human B cell differentiation and inhibits B cell receptor-mediated signaling. Eur. J. Immunol. 1999, 29, 3160–3167. [Google Scholar] [CrossRef]
  57. Chen, Z.; Shojaee, S.; Buchner, M.; Geng, H.; Lee, J.W.; Klemm, L.; Titz, B.; Graeber, T.G.; Park, E.; Tan, Y.X.; et al. Signalling thresholds and negative B-cell selection in acute lymphoblastic leukaemia. Nature 2015, 521, 357–361. [Google Scholar] [CrossRef] [PubMed]
  58. Limnander, A.; Depeille, P.; Freedman, T.S.; Liou, J.; Leitges, M.; Kurosaki, T.; Roose, J.P.; Weiss, A. STIM1, PKC-δ and RasGRP set a threshold for proapoptotic Erk signaling during B cell development. Nat. Immunol. 2011, 12, 425–433. [Google Scholar] [CrossRef] [PubMed]
  59. Monti, S.; Savage, K.J.; Kutok, J.L.; Feuerhake, F.; Kurtin, P.; Mihm, M.; Wu, B.; Pasqualucci, L.; Neuberg, D.; Aguiar, R.C.T.; et al. Molecular profiling of diffuse large B-cell lymphoma identifies robust subtypes including one characterized by host inflammatory response. Blood 2005, 105, 1851–1861. [Google Scholar] [CrossRef] [Green Version]
  60. Williams, K.A.; Lee, M.; Hu, Y.; Andreas, J.; Patel, S.J.; Zhang, S.; Chines, P.; Elkahloun, A.; Chandrasekharappa, S.; Gutkind, J.S.; et al. A systems genetics approach identifies CXCL14, ITGAX, and LPCAT2 as novel aggressive prostate cancer susceptibility genes. PLoS Genet. 2014, 10, e1004809. [Google Scholar] [CrossRef]
  61. Umit, E.G.; Baysal, M.; Durmus, Y.; Demir, A.M. CD11c expression in chronic lymphocytic leukemia revisited, related with complications and survival. Int. J. Lab. Hematol. 2017, 39, 552–556. [Google Scholar] [CrossRef] [PubMed]
  62. Wang, Y.; Xu, B.; Hu, W.-W.; Chen, L.-J.; Wu, C.-P.; Lu, B.-F.; Shen, Y.-P.; Jiang, J.-T. High expression of CD11c indicates favorable prognosis in patients with gastric cancer. World J. Gastroenterol. 2015, 21, 9403–9412. [Google Scholar] [CrossRef] [PubMed]
  63. Ohl, K.; Schippers, A.; Tenbrock, K. CD11c-Specific Deletion Reveals CREB as a Critical Regulator of DC Function during the Germinal Center Response. J. Immunol. Res. 2018, 2018, 8947230. [Google Scholar] [CrossRef] [PubMed]
  64. Chevalier, N.; Mueller, M.; Mougiakakos, D.; Ihorst, G.; Marks, R.; Schmitt-Graeff, A.; Veelken, H. Analysis of dendritic cell subpopulations in follicular lymphoma with respect to the tumor immune microenvironment. Leuk. Lymphoma 2016, 57, 2150–2160. [Google Scholar] [CrossRef]
  65. Perez-Andres, M.; Paiva, B.; Nieto, W.G.; Caraux, A.; Schmitz, A.; Almeida, J.; Vogt, R.F.; Marti, G.E.; Rawstron, A.C.; Van Zelm, M.C.; et al. Human peripheral blood B-cell compartments: A crossroad in B-cell traffic. Cytom. Part B Clin. Cytom. 2010, 78, S47–S60. [Google Scholar] [CrossRef] [PubMed]
  66. Perbellini, O.; Falisi, E.; Giaretta, I.; Boscaro, E.; Novella, E.; Facco, M.; Fortuna, S.; Finotto, S.; Amati, E.; Maniscalco, F.; et al. Clinical significance of LAIR1 (CD305) as assessed by flow cytometry in a prospective series of patients with chronic lymphocytic leukemia. Haematologica 2014, 99, 881–887. [Google Scholar] [CrossRef] [Green Version]
  67. Menendez, P.; Vargas, A.; Bueno, C.; Barrena, S.; Almeida, J.; De Santiago, M.; López, A.; Roa, S.; San Miguel, J.F.; Orfao, A. Quantitative analysis of bcl-2 expression in normal and leukemic human B-cell differentiation. Leukemia 2004, 18, 491–498. [Google Scholar] [CrossRef]
  68. Kitano, M.; Moriyama, S.; Ando, Y.; Hikida, M.; Mori, Y.; Kurosaki, T.; Okada, T. Bcl6 protein expression shapes pre-germinal center B cell dynamics and follicular helper T cell heterogeneity. Immunity 2011, 34, 961–972. [Google Scholar] [CrossRef]
  69. Kurosu, T.; Fukuda, T.; Miki, T.; Miura, O. BCL6 overexpression prevents increase in reactive oxygen species and inhibits apoptosis induced by chemotherapeutic reagents in B-cell lymphoma cells. Oncogene 2003, 22, 4459–4468. [Google Scholar] [CrossRef] [Green Version]
  70. Dupont, T.; Yang, S.N.; Patel, J.; Hatzi, K.; Malik, A.; Tam, W.; Martin, P.; Leonard, J.; Melnick, A.; Cerchietti, L. Selective targeting of BCL6 induces oncogene addiction switching to BCL2 in B-cell lymphoma. Oncotarget 2016, 7, 3520–3532. [Google Scholar] [CrossRef]
Figure 1. Algorithm for the identification of DLBCL derived-cell lines using multi-parameter flow cytometry. Diagram showing the identification strategy based on the mutually exclusive expression of the kappa (a) or lambda (b) light chain of immunoglobulins.
Figure 1. Algorithm for the identification of DLBCL derived-cell lines using multi-parameter flow cytometry. Diagram showing the identification strategy based on the mutually exclusive expression of the kappa (a) or lambda (b) light chain of immunoglobulins.
Jcm 08 01074 g001
Figure 2. CD39 expression discriminates between GCB and ABC DLBCL primary tumors and cell lines. Cluster gram showing the expression of MME, BCL6, ENTPD1, and BCL2 genes in ABC DLBCL samples (n = 167) and GCB DLBCL samples (n = 183) in the ABC and GCB DLBCL samples from the R-CHOP and CHOP Lenz cohorts (gene expression profile, GEP, data) (a). (B) Boxplots illustrate the expression signal (GEP) of ENPD1 in ABC DLBCL samples and in GCB DLBCL samples in the Lenz cohort (b). Boxplots showing CD39 mean fluorescence intensity (MFI, logarithmic scale) in 16 DLBCL derived-cell lines (c). The boxes indicate the 25th and 75th percentile values. The line in the middle corresponds to the median. The vertical lines, to the 10th and the 90th percentiles. The outliers are identified as the third quartile plus 1.5 IQR (interquartile range); * p-value < 0.05 (Mann-Whitney U-test). GCB: germinal center B cell–like, ABC: Activated B cell–like, DLBCL: diffuse large B cell lymphoma.
Figure 2. CD39 expression discriminates between GCB and ABC DLBCL primary tumors and cell lines. Cluster gram showing the expression of MME, BCL6, ENTPD1, and BCL2 genes in ABC DLBCL samples (n = 167) and GCB DLBCL samples (n = 183) in the ABC and GCB DLBCL samples from the R-CHOP and CHOP Lenz cohorts (gene expression profile, GEP, data) (a). (B) Boxplots illustrate the expression signal (GEP) of ENPD1 in ABC DLBCL samples and in GCB DLBCL samples in the Lenz cohort (b). Boxplots showing CD39 mean fluorescence intensity (MFI, logarithmic scale) in 16 DLBCL derived-cell lines (c). The boxes indicate the 25th and 75th percentile values. The line in the middle corresponds to the median. The vertical lines, to the 10th and the 90th percentiles. The outliers are identified as the third quartile plus 1.5 IQR (interquartile range); * p-value < 0.05 (Mann-Whitney U-test). GCB: germinal center B cell–like, ABC: Activated B cell–like, DLBCL: diffuse large B cell lymphoma.
Jcm 08 01074 g002
Figure 3. Eight lymphoid markers could predict the outcome of patients with DLBCL. Kaplan-Meier curves indicated that, in the training cohort (R-CHOP Lenz cohort, n = 233), high expression of BCL2, LAIR1, CD39, or CD103 is associated with a poor outcome (i.e., overall survival in function of time) (a), whereas high expression of BCL6, CD81, CD10, or Cd11c is associated with better overall survival (b). Red, overexpression, and blue, downregulation. Curves were compared with the log rank test. HR: hazard ratio.
Figure 3. Eight lymphoid markers could predict the outcome of patients with DLBCL. Kaplan-Meier curves indicated that, in the training cohort (R-CHOP Lenz cohort, n = 233), high expression of BCL2, LAIR1, CD39, or CD103 is associated with a poor outcome (i.e., overall survival in function of time) (a), whereas high expression of BCL6, CD81, CD10, or Cd11c is associated with better overall survival (b). Red, overexpression, and blue, downregulation. Curves were compared with the log rank test. HR: hazard ratio.
Jcm 08 01074 g003
Figure 4. The risk score predicts the outcome of patients with DLBCL in four independent cohorts. Kaplan-Meier curves of overall survival in function of time in the R-CHOP Lenz cohort (a), CHOP Lenz cohort (b), R-CHOP Melnick cohort (c), and FFPE RCHOP cohort (d). Data were compared with the log rank test.
Figure 4. The risk score predicts the outcome of patients with DLBCL in four independent cohorts. Kaplan-Meier curves of overall survival in function of time in the R-CHOP Lenz cohort (a), CHOP Lenz cohort (b), R-CHOP Melnick cohort (c), and FFPE RCHOP cohort (d). Data were compared with the log rank test.
Jcm 08 01074 g004
Figure 5. Comparison of BCL2 (a), BCL6 (b), LAIR1 (c), and CD11c (d) gene expression in DLBCL samples and in normal centrocytes and centro-blasts (GSE12195 dataset). The box-plot diagrams include the median value and the interquartile rage (IQR). The error bars represent the minimum values under the median and the outliers are identified as the third quartile plus 1.5 IQR (R I386 3.4.0 software, R Foundation). Results were compared using the Mann-Whitney U-test.
Figure 5. Comparison of BCL2 (a), BCL6 (b), LAIR1 (c), and CD11c (d) gene expression in DLBCL samples and in normal centrocytes and centro-blasts (GSE12195 dataset). The box-plot diagrams include the median value and the interquartile rage (IQR). The error bars represent the minimum values under the median and the outliers are identified as the third quartile plus 1.5 IQR (R I386 3.4.0 software, R Foundation). Results were compared using the Mann-Whitney U-test.
Jcm 08 01074 g005
Figure 6. BCL6 expression correlates with the response to etoposide in DLBCL cell lines. Box-plots illustrate the inhibitory concentration 50 (IC50) of etoposide on a logarithmic scale. The boxes represent the 25th and 75th percentile values. The line in the middle corresponds to the median. The vertical lines correspond to the 10th and the 90th percentiles, and the circles, to the outliers. * p value < 0.05 (Mann-Whitney U-test) (a). Linear regression analysis of BCL6 expression (MFI) in function of the IC50 of etoposide in 16 DLBCL cell lines. r represents the Spearman correlation coefficient, p value <0.05 (Spearman correlation test) (b).
Figure 6. BCL6 expression correlates with the response to etoposide in DLBCL cell lines. Box-plots illustrate the inhibitory concentration 50 (IC50) of etoposide on a logarithmic scale. The boxes represent the 25th and 75th percentile values. The line in the middle corresponds to the median. The vertical lines correspond to the 10th and the 90th percentiles, and the circles, to the outliers. * p value < 0.05 (Mann-Whitney U-test) (a). Linear regression analysis of BCL6 expression (MFI) in function of the IC50 of etoposide in 16 DLBCL cell lines. r represents the Spearman correlation coefficient, p value <0.05 (Spearman correlation test) (b).
Jcm 08 01074 g006
Table 1. Cox univariate and multivariate analyses of prognostic factors of overall survival in patients with diffuse large B-cell lymphoma (R-CHOP Lenz cohort, n = 233).
Table 1. Cox univariate and multivariate analyses of prognostic factors of overall survival in patients with diffuse large B-cell lymphoma (R-CHOP Lenz cohort, n = 233).
A.Overall Survival (n = 233)
Prognostic VariableHRp-Value
Age (>60 years)2.20<0.0001
GCB-ABC molecular subgroups2.75<0.0001
IPI1.79<0.0001
DNA repair score3.87<0.0001
Risk Score2.41<0.0001
B.Overall Survival (n = 233)
Prognostic VariableHRp-Value
Risk Score2.248<0.0001
Age (>60 years)1.8460.03
Risk Score2.269<0.0001
GCB-ABC molecular subgroups1.359NS
Risk Score2.113<0.0001
IPI1.582<0.0001
Risk Score2.263<0.0001
DNA repair score2.909<0.0001
C.Overall Survival (n = 233)
Prognostic VariableHRp-Value
Age (>60 years)1.80NS
GCB-ABC molecular subgroups2.23NS
IPI0.20NS
DNA repair score3.44<0.0001
Risk Score1.610.007
The indicated prognostic factors were tested individually (A), two by two (B), or in multivariate analysis (all variables) (C) using a Cox regression model. p-values and hazard ratios (HR) are shown. NS: not significant at the 5% threshold. The IPI groups were defined as follows: low risk group = IPI score 0 or 1, low-intermediate risk group = IPI score 2, high-intermediate risk group = IPI score 3, and high-risk group = IPI score 4 or 5. IPI, international prognostic index. GCB, germinal-center B-cell–like subgroup. ABC, activated B cell–like subtype.
Table 2. Cox univariate and multivariate analyses of the effect of lymphoid markers on overall survival in patients with diffuse large B-cell lymphoma (gene expression data from the R-CHOP Lenz cohort, n = 233).
Table 2. Cox univariate and multivariate analyses of the effect of lymphoid markers on overall survival in patients with diffuse large B-cell lymphoma (gene expression data from the R-CHOP Lenz cohort, n = 233).
A.Overall Survival (n = 233)
Lymphoid MarkerHRp-Value
BCL2 2.110.006
BCL60.26<0.0001
CD392.10.01
CD1031.820.03
CD11c0.27<0.0001
LAIR12.240.004
CD100.37<0.0001
CD50.590.04
B.Overall Survival (n = 233)
Lymphoid MarkerHRp-Value
BCL21.960.01
BCL60.440.007
CD391.22NS
CD1031.30NS
CD11c0.33<0.0001
LAIR12.820.002
CD100.39NS
CD50.470.01
The indicated prognostic factors were tested as single variables (A) or multi-variables (B) using a Cox-regression model. p-values and hazard ratios (HR) are shown. NS: not significant at the 5% threshold.

Share and Cite

MDPI and ACS Style

Devin, J.; Kassambara, A.; Bruyer, A.; Moreaux, J.; Bret, C. Phenotypic Characterization of Diffuse Large B-Cell Lymphoma Cells and Prognostic Impact. J. Clin. Med. 2019, 8, 1074. https://doi.org/10.3390/jcm8071074

AMA Style

Devin J, Kassambara A, Bruyer A, Moreaux J, Bret C. Phenotypic Characterization of Diffuse Large B-Cell Lymphoma Cells and Prognostic Impact. Journal of Clinical Medicine. 2019; 8(7):1074. https://doi.org/10.3390/jcm8071074

Chicago/Turabian Style

Devin, Julie, Alboukadel Kassambara, Angélique Bruyer, Jérôme Moreaux, and Caroline Bret. 2019. "Phenotypic Characterization of Diffuse Large B-Cell Lymphoma Cells and Prognostic Impact" Journal of Clinical Medicine 8, no. 7: 1074. https://doi.org/10.3390/jcm8071074

APA Style

Devin, J., Kassambara, A., Bruyer, A., Moreaux, J., & Bret, C. (2019). Phenotypic Characterization of Diffuse Large B-Cell Lymphoma Cells and Prognostic Impact. Journal of Clinical Medicine, 8(7), 1074. https://doi.org/10.3390/jcm8071074

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop