Original Article |
1 Laboratory of Experimental Hematology, Faculty of Medicine, Institut Jules Bordet, Université Libre de Bruxelles (ULB), Brussels
2 Functional Genomics and Translational Research Unit, Faculty of Medicine, Institut Jules Bordet, Université Libre de Bruxelles (ULB), Brussels
3 Machine Learning Group, Faculty of Sciences, Université Libre de Bruxelles (ULB), Brussels, Belgium
Correspondence: Basile Stamatopoulos, Université Libre de Bruxelles, Institut Jules Bordet, Laboratoire d'Hématologie Expérimentale, Boulevard de Waterloo 121, 1000 Bruxelles, Belgium., E-mail:bstamato{at}ulb.ac.be
|
|
|---|
Design and Methods: In an attempt to unveil unfavorable cellular properties linked to high ZAP70 expression, we used gene expression profiling to identify genes associated with disparities in B cells from chronic lymphocytic leukemia patients expressing high versus low ZAP70 mRNA, measured by quantitative real-time PCR. Two groups of 7 patients were compared, selected on the basis of either high or low ZAP70 mRNA expression.
Results: Twenty-seven genes were differentially expressed with an FDR<10%, and several genes were significant predictors of treatment-free survival (TFS) and/or overall survival; PDE8A and FCRL family genes (down-regulated in ZAP70+ patients) could predict TFS and overall survival; ITGA4 mRNA (up-regulated in ZAP70+ patients) could significantly predict overall survival. Importantly, gene set enrichment analysis revealed overrepresentation of adhesion/migration genes. We therefore investigated in vitro adhesion/migration capacity of chronic lymphocytic leukemia cells into a stromal microenvironment or in response to conditioned medium. We showed that ZAP70+ cells had better adhesion/migration capacities and only ZAP70+ patient cells responded to microenvironment contact by CXCR4 downregulation.
Conclusions: We concluded that several prognostic factors are the reflection of microenvironment interactions and that the increased adhesion/migratory capacity of ZAP70+ cells in their microenvironment can explain their better survival and thus the aggressiveness of the disease.
Key words: chronic lymphocytic leukemia, ZAP70, microarrays, microenvironment, prognosis.
|
|
|---|
-associated protein 70)1 and LPL (lipoprotein lipase).7 ZAP70 expression in CLL correlates strongly with IgVH mutation status8 and seems to be one of the most promising prognostic factors for future clinical use. In 2003, Wiestner et al. found a 93% correlation between mutational status and ZAP70:1 patients with 20% of their B cells expressing ZAP70 by flow cytometry (FC) (i.e. at the same level as T cells) were generally unmutated; on the contrary, patients with less than 20% of their B cells expressing ZAP70 generally had mutated IgVH. ZAP70 status determination by FC, however, is often inaccurate at the positivity limit because of the low resolution between the positive (+) and negative (–) populations due to the gating procedure and the choice of antibody.9 To offset these drawbacks, ZAP70 mRNA absolute quantification has been proposed. In a previous study, we showed that ZAP70 mRNA expression was highly correlated with IgVH mutational status, in addition to survival and treatment-free time.10 Moreover, we and others suggested that ZAP70 could be an even better prognostic factor than mutational status.10–12 We therefore used real-time RT-PCR quantification (qPCR) to select patients with high and low ZAP70 expression and compared their gene expression profiles in order to assess the biological significance of ZAP70 expression in CLL. We hypothesized that a gene expression profile comparison of these two CLL subsets widely differing ZAP70 expression would reveal gene expression patterns associated not only with prognosis, but also with biology, particularly with regard to microenvironment interactions implicated in the survival of leukemic cells, a relevant aspect of the disease.
|
|
|---|
|
View this table: [in a new window] [Download PPT slide] |
Table 1. Characteristics of patients with low and high ZAP70 expression.
|
Bioinformatic analysis
A comparative gene expression profile was determined in 14 patients (7 ZAP70high and 7 ZAP70low patients). We identified significant differences between sample groups using BRB array tools (Biometric Research Branch, National Cancer Institute). Only genes defined as present by the Affymetrix algorithm in at least 30% of either of the two groups were considered for further analysis, and we calculated two-sample t tests (with a random variance model) of the two groups for each gene. We then addressed the multiple comparison problems by estimating the false discovery rate (FDR) in a simple manner as the ratio of expected number of false positives at this given p value threshold to the number of positives actually found. Using BRB gene set expression comparison tools, overrepresentation of gene ontology (GO) categories, Biocarta, KEGG, and Broad/MIT pathways were investigated by the Hotelling T-square test.
Flow cytometry analysis, IgVH gene mutational and cytogenetic abnormalities determination
Cytoplasmic ZAP70 protein by FC and IgVH mutational status were performed as previously described10 (see also Online Supplementary Appendix). We evaluated the expression of CD38, CXCR4 (chemokine (C-X-C motif) receptor 4 or CD184) ) and CD69 on the cell surface by FC in a CD19+ gate with a panel of fluorochrome-labeled monoclonal antibodies (phycoerythrin-conjugated, Immunosource). Classical cytogenetic abnormalities by standard karyotype analysis were investigated for the 14 patients included in gene expression profile analysis. Additional interphase FISH was performed to screen for most common aberrations using Chromoprobe Multiprobe® - CLL System (Cytocell, Amplitech, Compiegne, France).
Real-time PCR analysis
We used 25 ng of cDNA (produced by a standard reverse transcription) in a qPCR reaction with SYBR® Green PCR Master Mix (Applied Biosystems, Rotterdam, The Netherlands) and 0.32 µM of gene-specific forward and reverse primers (Invitrogen Life technologies, Merelbeke, Belgium). Primer sequences are listed in the Online Supplementary Table S1. Standard real-time PCR was performed as previously described10 (Online Supplementary Appendix).
Co-culture of mesenchymal stromal cells and chronic lymphocytic leukemia cells
MSCs were generated from bone marrow (BM) aspirates of normal volunteers after obtaining informed consent, as previously described.13 Passage one or two MSCs were used for the co-culture experiment. In order to study the influence of soluble factors produced by MSCs, conditioned serum-free medium (CM) from the culture of 1-day-old MSCs was prepared. Peripheral blood mononuclear cells from CLL ZAP70 + and – patients and also from healthy donors were resuspended in RPMI-1640 + 10% FBS (Biowhittaker, Verviers, Belgium) at a final concentration of 106/mL; 2 mL of cell suspension was placed in the wells of a 6-well plate or seeded on MSCs in identical plates. In parallel experiments, cells were seeded onto 24 mm Transwell diffusion chambers (0.4 µm microporous filter; Corning Incorporated, NY, USA) and placed into stroma-coated 6-well plates. ZAP70, CD69, and CXCR4 were all evaluated by FC in CLL cells cultured alone, with, and without contact with stromal cells. Moreover, in the case of contact with stromal cells, we separated cells in suspension in the medium from cells adhering to the stromal microenvironment in order to measure the above described markers. We also measured CXCR4 expression in fresh CLL blood samples from ZAP70 + and – patients.
Adhesion and migration assay
CLL cell suspensions were incubated in fibronectin (BDbioscience, Erembodegem, Belgium) precoated 24-well plates. After four hours of incubation, adherent cells were collected by trypsinization, concentrated, and counted by a Trypan blue exclusion assay (Invitrogen): 106 cells were also resuspended in 200 µL of medium and plated in the upper chamber of a 6.6 mm diameter Transwell culture insert in bare polycarbonate with a 5 µm pore size (Corning Incorporated, NY, USA). The lower chamber of each well contained 500 µL of medium alone, with SDF1
(stromal cell-derived factor 1
also known as CXCL12 -chemokine (C-X-C motif) ligand 12) (200 ng/mL) (R&D Systems, Minneapolis, MN, USA) or MSC-CM. After four hours of incubation, cells were recovered from the lower chamber, concentrated and counted using a Trypan blue exclusion assay. Migration index was calculated as number of cells transmigrating in the presence of the chemoattractant per number of cells transmigrated in absence of the chemoattractant. In the cells found in the lower chamber, ZAP70 was measured by FC in order to compare its expression in migrating versus non-migrating cells in response to MSC-CM.
Statistical analysis
ROC curve analyses were performed to determine the "+" and "–" status cut-offs (for all studied variables) that best distinguished between + and – ZAP70 cases. Treatment-free survival (TFS) and overall survival (OS) distributions were plotted using Kaplan-Meier estimates and were compared using the log-rank test. Significant differences were evaluated using the Wilcoxon matched pairs or Mann Whitney test. All tests were two-sided. An effect was considered to be statistically significant at p<0.05, and all analyses were performed with Prism GraphPad 5.0 software.
|
|
|---|
2(1)=6.132]. Moreover, the median OS was also significantly different [p=0.0018;
2(1)=9.701]. ZAP70 average expression assessed by qPCR was 23.57±8.71 in the low group and 1096±131.50 in the high group (p<0.0001). These results were confirmed in a patient cohort of 85 patients with a median follow-up of 74 months (range 8–299 months) (Figure 1A and B). This cohort was derived from the 108 cohort patient cohort previously published10 for which enough RNA was available. The characteristics of these patients are stated in the Online Supplementary Table S2. ZAP70 expression was expressed as fold change of ZAP70 in the Namalwa cell line and normalized with cyclophilin expression. Patients expressing ZAP70 with a fold change above 114 were deemed "+".
![]() View larger version (28K): [in a new window] [Download PPT slide] |
Figure 1. Disease progression and underlying gene expression in chronic lymphocytic leukemia patients with high versus low ZAP70 expression. ZAP70 expression was expressed as fold change of ZAP70 in the Namalwa cell line and normalized with cyclophilin expression. Patients expressing ZAP70 more than 114 fold were deemed "+" (n=41) while patients expressing ZAP70 with less than 114-fold were deemed "–" (n=44). (A) Treatment-free survival as a function of ZAP70 measured by qPCR. (B) Overall Survival as a function of ZAP70 measured by qPCR plotted using Kaplan-Meier analysis and compared using the log-rank test. (C) Multidimensional scaling of samples using the ZAP70 gene expression signature (n=39, FDR<10%). Red samples represent ZAP70high and black samples represent ZAP70low patients. (D) Hierarchical clustering of 39 probe sets discriminated between ZAP70high and ZAP70low chronic lymphocytic leukemia samples. Color scale indicates units of standard deviation from the mean expression of each gene.
|
After applying an FDR, we found only 39 probe sets representing 27 different genes with an FDR
10% (Figure 1D). Of these, 13 genes were up-regulated and 14 were down-regulated in the ZAP70high group. The two ZAP70 probe sets were, not surprisingly, at the top of this list with an FDR<10–7. As expected, a multidimensional scaling of samples as well as a hierarchical clustering using this ZAP70-specific gene signature exhibited a clear separation of ZAP70high and ZAP70low patients (Figure 1C and D).
The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE12734(http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE12734). In order to validate this pattern, a panel of six genes in the list of 27 (with FDR<10%), five genes among the list of 135 probe sets (with p<0.001) and two genes in the list of 937 probe sets (p<0.05) were selected for qPCR validation on CD19+ purified samples in an extended cohort of 85 patients with a median follow-up of 74 months (range 8–299 months). These 13 randomly chosen genes were all confirmed as differentially expressed (Table 2, Online Supplementary Figure S1).
|
View this table: [in a new window] [Download PPT slide] |
Table 2. Real-time PCR validation of selected genes differentially expressed between ZAP70 + and – patients.
|
Gene set expression comparison reveals interaction differences with the microenvironment
Using BRB array tools, we performed gene set enrichment analysis investigating gene GO categories, in addition to Biocarta, KEGG, and Broad/MIT pathways that had a higher than expected number of genes differentially expressed between our two classes of samples using Hotelling T-square test. In other words, among differentially expressed genes we examined whether an enrichment/abnormal rate of genes involved in a same pathway/function is overrepresented in comparison with what is expected by hazard.
Interestingly, several GO categories/pathways involved in migration, motility, adhesion, cytoskeleton, actin modification, and microenvironment interaction (via the CXCR4 pathway) were significantly over-represented (Table 3). Additional details on involved genes are available in the Online Supplementary Table S4.
|
View this table: [in a new window] [Download PPT slide] |
Table 3. Over-represented gene list between ZAP70high and ZAP70low patients.
|
Migration in response to SDF1
was very low in comparison to MSC-CM (p<0.0001). However, no significant difference in terms of migration was found between ZAP70+ (n=12) and – (n=12) patients (Figure 2A) but there was a clear trend towards greater migration in cells from ZAP70+ patients. We also measured and compared ZAP70 expression using FC in cells from the upper and lower chamber: the number of ZAP70+ cells was significantly increased in the migrating cell population (n=20; p=0.0019) (Figure 2B). We verified that ZAP70 was not induced by CM (data not shown), indicating that this difference was linked to migrating cells.
![]() View larger version (27K): [in a new window] [Download PPT slide] |
Figure 2. Migration of chronic lymphocytic leukemia cells in response to SDF1 and MSC conditioned medium and modulation of ZAP70, CXCR4 and CD69 in response to the stromal microenvironment. (A) Results are expressed as the mean of migration index in response to SDF1 or MSC conditioned medium (CM). (B) ZAP70 measured by FC in the upper and lower migration chamber. Chronic lymphocytic leukemia cells from ZAP70– (C) and + (D) patients were co-plated either with stromal cells, with stromal cells separated by a 0.45 µm transwell not allowing contact between the two cell types, or alone. After 4 h of incubation, ZAP70, CXCR4 and CD69 were measured by FC on adherent and non-adherent cells. Furthermore, ZAP70 was measured on chronic lymphocytic leukemia cells cultured alone and in transwell conditions.
|
0.0223) and + (n=10, p
0.0180) patient cells but also in normal B cells isolated from healthy donors (n=10, p
0.0020; Online Supplementary Figure S2). Moreover, there was no statistical difference between ZAP70 expression among non-adherent cells, cells separated by a transwell, or cells cultured in medium (Figure 2C and D, Online Supplementary Figure S2). Although CXCR4 mRNA was not differentially expressed between our two groups, CXCR4 signaling pathway genes were over-represented. Therefore, we also measured CXCR4, the receptor for SDF1
, on adherent B cells. Very interestingly, CXCR4 cell surface expression between adherent and non-adherent cells was only significantly down-regulated in ZAP70+ patients (% of cells: p=0.0059; Molecules of Equivalent Soluble Fluorochrome (MESF): p=0.0098), but not in ZAP70– patients (% of cells: p=0.1952; MESF: p=0.9219) (Figure 2C and D). Furthermore, when CXCR4 was measured in 54 fresh CLL blood samples (29 ZAP70– and 25 ZAP70+), ZAP70+ patients expressed significantly less CXCR4 molecules per cell (p=0.0183). The characteristics of these patients are stated in the Online Supplementary Table S5. We also observed that CD69 mRNA was significantly increased in ZAP70+ patients (Online Supplementary Table S2) and that CD69 protein was up-regulated after stromal cell contact (Figure 2C and 2D).
|
|
|---|
More interestingly, in addition to their prognostic value, the expression of these factors is linked to microenvironment crosstalk. LPL expression could be induced by B-cell receptor (BCR) stimulation in CLL cells but not in normal B cells.25 Several studies have shown significantly higher LPL expression in IgVH-unmutated patients7,10,26 which is associated with autoreactive BCR activity.27 Taken together, our data strongly suggest an important role of continual stimulation by the tumoral microenvironment in patients with a poor prognosis. Fc receptor-like (FCRL) genes (also known as FcRH, IRTA, IFGP, SPAP) were up-regulated in ZAP70– patients. These genes belong to a large family of lymphocyte receptors with immunoreceptor tyrosine–based inhibition motifs (ITIM).28,29 These ITIMs can be phosphorylated following external stimuli and can subsequently recruit phosphatases. They can therefore have an inhibitory role in BCR30 or in MAPK signaling.31 Stimuli given by the microenvironment could confer growth advantage and extended survival to leukemic cells.32 Moreover, BCR stimulation is linked to cell survival, activation, and G1 progression.33 All these finding are in line with higher expression of FCRL genes in ZAP70– patients, with potential ITIM inhibitory pathways explaining their lower response to external survival stimuli, which would be associated with a good prognosis. Finally, we also showed that CD69 (recently proposed as a strong predictor of CLL prognosis)24 could be up-regulated while CXCR4 could be down-regulated by microenvironment contact. Similarly, Ocana et al. demonstrated that low CXCR3 was associated with Rai disease stages III and IV34 and recently Ding et al. showed that CD38 expression (in CLL cells with evidence for CD38 expression) was up-regulated after two weeks in contact with MSCs. In contrast, CD38 expression remained unchanged in cells with minimal CD38 expression after co-culture with MSC.35 Considering all these data, we can tentatively conclude that expression of markers such as ZAP70, LPL, CD38, CD69, FCRL, CXCR4, and CXCR3 are probably linked to the microenvironment, and classification of patients into poor or good prognosis groups with regard to these factors seems to be a reflection of microenvironment interactions.
Gene set enrichment analysis revealed several pathways and GO categories linked to migration, motility, adhesion, cytoskeleton, actin modification, and microenvironment interaction (via the CXCR4 pathway). These results are in line with findings by Deaglio et al. who reported a genetic signature based simply on migration index and independent of molecular factors.36 This signature is composed of genes involved in the control of cell motility, adhesion, cell-cell contact and cytoskeletal organization and also included ZAP70 indicating that ZAP70 is linked to migration, thus reinforcing our hypothesis. The microenvironment is considered to be crucial in determining expansion37 and survival38 of CLL cells. Therefore, we investigated the effects of microenvironment in a BM-MSC model, which we previously established.38 We showed that ZAP70+ cells whether from ZAP70+ or ZAP70– patients or normal B cells isolated from healthy donors adhered significantly more to MSCs, indicating that ZAP70 played an important role in the chemotaxis and/or the adhesion process, even in normal B cells. Moreover, ITGA4 is significantly over-expressed on cells from ZAP70+ patients. Interestingly, this gene has been implicated in migration into lymph nodes39 and is modulated by SDF1
.40 In CLL, the interaction of ITGA4 with fibronectin enhances BCL2 expression41 and influences chemosusceptibility.42 Furthermore, our qPCR data indicated that ITGA could significantly predict OS (p=0.0075).
Surprisingly, we and others36 did not find significant differences in cell migration in response to SDF1
between cells from ZAP70 + and – patients. However, we observed that ZAP70+ cells migrated significantly more than ZAP70 – cells when we compared the upper and lower chamber cells in our migration assay. Although there were no significant migratory differences between ZAP70 + and –patients, there was a clear trend towards increased migration of cells from ZAP70+ patients in response to MSC-CM. The better responsiveness of ZAP70+ cells to their microenvironment that we observed is in line with the observation of Richardson et al. who showed that ZAP70 identifies cells with an increased propensity to migrate to the lymph node and an increased ability to respond to survival signals.43 Moreover, we showed that only ZAP70+ patients respond to microenvironment contact by a downmodulation of CXCR4 (probably by protein internalization, as previously described)44 indicating once more that only ZAP70+ patients could respond to microenvironment stimuli and suggesting an association between aggressiveness of the disease and chemotaxis towards the microenvironment via the SDF1
/CXCR4 pathway. Furthermore, the role of ZAP70 in migration has already been demonstrated in T cells.45,46 Indeed, expression of wild type ZAP70 in ZAP70-deficient cells (P116 neo) strongly enhanced its migratory capability and this remains true even in absence of CD3 complex surface expression and thus of functional TCR
subunits as is the case in B cells.46 Tacchini et al. also showed that SDF1
induces a rapid and transient tyrosine phosphorylation of ZAP70, suggesting that this chemokine can directly activate ZAP70. We could thus speculate that, in CLL cells, microenvironment stimuli could transitorily activate ZAP70 resulting in downstream cascade phosphorylation. Stimulation by the SDF1
/CXCR4 pathway could activate ZAP70-dependent genes (such as Vav146 which contributes to signal amplification and diversification events)47 leading to actin cytoskeleton changes and thus chemotaxis.
In conclusion, our data demonstrated that in ZAP70+ patients, leukemic cells can better interact and cross-talk with their protective microenvironment, explaining their increased survival and the aggressiveness of the disease.
BS performed research and statistical analysis, analyzed data, made figures and tables, and wrote manuscript. BH-K performed and revised statistical analysis. CE contributed to microarray experiments. AS performed real-time PCR confirmations. CDB contributed to flow cytometry experiments. DH. initiated project and performed research. NM, DB and PM contributed to patient samples and data. LL supervised and designed research, and revised manuscript.
The authors declare no competing financial interests.
Funding: this work was financed by F.R.I.A. grant (Fonds de Recherche pour lIndustrie et lAgriculture) and the Télévie fund, both of which are affiliated with the F.R.S-F.N.R.S. (Fonds de la Recherche Scientifique, FNRS).
Received for publication October 27, 2008. Revision received January 20, 2009. Accepted for publication January 27, 2009.
|
|
|---|
4 integrin are important for migration of chronic lymphocytic leukemia cells into lymph nodes. Blood 2002;99:2977-84.
modulates VLA-4 integrin-mediated multiple myeloma cell adhesion to CS-1/fibronectin and VCAM-1. Blood 2001;97:346-51.
4β1 integrin by fibronectin induces in vitro resistance of B chronic lymphocytic leukemia cells to fludarabine. J Leukoc Biol 2002;71:495-502.This article has been cited by other articles:
![]() |
S. Deaglio and F. Malavasi Chronic lymphocytic leukemia microenvironment: shifting the balance from apoptosis to proliferation Haematologica, June 1, 2009; 94(6): 752 - 756. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||