Purpose: Whole-genome sequencing has revealed MYD88 L265P and CXCR4 mutations (CXCR4mut) as the most prevalent somatic mutations in Waldenström macroglobulinemia. CXCR4 mutation has proved to be of critical importance in Waldenström macroglobulinemia, in part due to its role as a mechanism of resistance to several agents. We have therefore sought to unravel the different aspects of CXCR4 mutations in Waldenström macroglobulinemia.

Experimental Design: We have scanned the two coding exons of CXCR4 in Waldenström macroglobulinemia using deep next-generation sequencing and Sanger sequencing in 98 patients with Waldenström macroglobulinemia and correlated with SNP array landscape and mutational spectrum of eight candidate genes involved in TLR, RAS, and BCR pathway in an integrative study.

Results: We found all mutations to be heterozygous, somatic, and located in the C-terminal domain of CXCR4 in 25% of the Waldenström macroglobulinemia. CXCR4 mutations led to a truncated receptor protein associated with a higher expression of CXCR4. CXCR4 mutations pertain to the same clone as to MYD88 L265P mutations but were mutually exclusive to CD79A/CD79B mutations (BCR pathway). We identified a genomic signature in CXCR4mut Waldenström macroglobulinemia traducing a more complex genome. CXCR4 mutations were also associated with gain of chromosome 4, gain of Xq, and deletion 6q.

Conclusions: Our study panned out new CXCR4 mutations in Waldenström macroglobulinemia and identified a specific signature associated to CXCR4mut, characterized with complex genomic aberrations among MYD88L265P Waldenström macroglobulinemia. Our results suggest the existence of various genomic subgroups in Waldenström macroglobulinemia. Clin Cancer Res; 22(6); 1480–8. ©2015 AACR.

This article is featured in Highlights of This Issue, p. 1299

Translational Relevance

CXCR4 mutation has proved to be of critical importance in Waldenström macroglobulinemia, in part due to its role as a mechanism of resistance to several agents. We have therefore sought to unravel the different aspects of CXCR4 mutations in Waldenström macroglobulinemia and to characterize the genetic background using targeted next-generation sequencing (NGS) and SNP arrays. Mutational spectrum of eight candidate genes involved in TLR, RAS, and BCR pathway was studied in an integrative study. About 25% of Waldenström macroglobulinemia displays CXCR4 mutation in C-terminal domain responsible for higher CXCR4 expression. CXCR4 mutations pertain to the same clone as to MYD88 L265P mutations, with a clear-cut genomic signature traducing a more complex genome using SNP arrays but were mutually exclusive to BCR pathway mutations (CD79A/B mutations). CXCR4 mutation identified a genomic subgroup of MYD88L265P Waldenström macroglobulinemia.

Whole-genome sequencing has revealed CXCR4 as the second most frequent somatic mutation, identified in approximately 30% of Waldenström macroglobulinemia, similar to germline mutations found in the WHIM (warts, hypogammaglobulinemia, infections, and myelokathexis) syndrome (1–3). MYD88 L265P mutations remain the most frequent mutation reported in Waldenström macroglobulinemia to date, in nearly 90% of Waldenström macroglobulinemia (4, 5). Interestingly, MYD88 L265P may act as a founder mutation because of its high frequency in Waldenström macroglobulinemia; and thus it is suspected that a second hit may accelerate development and progression of the malignant clones leading to full-blown in Waldenström macroglobulinemia. Although the mechanism of deregulation of CXCR4 axis is not fully understood in Waldenström macroglobulinemia, it is possible that CXCR4 mutation might represent one of these secondary events.

C-X-C chemokine receptor type 4 (CXCR4) is a G-protein–coupled receptor that plays an important role in lymphopoiesis and cell trafficking (2, 6), along with its ligand, the stromal cell–derived factor-1 (CXCL12/SDF-1). The SDF1/CXCR4 axis promotes activation of several pathways including RAS, Akt, and NF-κB and interplays with BCR pathway (6–8). The CXCR4 gene is located on the long arm of chromosome 2 at position 21 and codes for a chemokine receptor that promotes migration and survival of various B lymphoid malignancies (8–10).

CXCR4 mutations were identified using Sanger sequencing by Treon and colleagues in nearly 25% of Waldenström macroglobulinemia, and the CXCR4 C1013G mutation was described as the most frequent recurrent CXCR4 mutation in 7% (11). However, the occurrence of potential subclonal of CXCR4 mutation and the proportion of CXCR4 mutation over the other known mutations was not available; in other words, how this mutation lies in the architecture of the MYD88-mutated clone. On contrary, Roccaro and colleagues have reported 30% incidence rate of the CXCR4 C1013G mutation using qPCR; however, they have solely studied this mutation (12).

In this study, we characterized the genomic landscape of CXCR4 somatic mutation in a large cohort of 98 patients with Waldenström macroglobulinemia, combining deep sequencing of target genes allowing study of the clonal architecture in the Waldenström macroglobulinemia cells, including subclonal analysis, and to confirm the incidence rate of these mutations and single polymorphism nucleotide array (SNPa) to decipher the genomic landscape of CXCR4-mutated Waldenström macroglobulinemia.

Patients

Ninety-eight patients diagnosed with Waldenström macroglobulinemia were included in this study (63 males, 35 females). The diagnosis and treatment initiation criteria of Waldenström macroglobulinemia were as published (13, 14). Among the cohort with available clinical follow up (n = 93), 68 patients were treated. First-line treatment was initiated at diagnosis in 41% of the patients. Front-line therapy included chloraminophene alone in 12 cases, fludarabine alone in 3 cases, rituximab (R) alone in 3 cases or in association with chemotherapy in 51 cases (dexamethasone-R- cyclophosphamide in 19 cases, bortezomib-R-dexamethasone in 7, R-bendamustine in 5 cases). Patients were untreated at time of bone marrow collection and gave informed consent prior to research sampling. No familial form of Waldenström macroglobulinemia was included.

Cell selection

Research sampling consisted of bone marrow samples and blood cells for all Waldenström macroglobulinemia. All samples were enriched using immunomagnetic beads (B cell isolation kit for tumoral cells from bone marrow sample, and Pan T cell isolation kit, Miltenyi-Biotec), following Ficoll-paque gradient centrifugation. The purity of samples was confirmed by multiparametre flow cytometry using a marker combination (including CD19, κ, λ, CD2, CD38, CD138, and CD27). More than 90% a light-chain isotype–positive B cells was detected in all cases (mean, 95.8% of tumoral cells; ref. 15).

DNA sequencing

DNA was extracted from isolated cells using Qui Amp kit (Sigma-Aldrich Co.). MYD88 L265P, CD79A, and CD79B mutations were analyzed as previously described (16) in 98 patients. C-terminal CXCR4 gene was amplified from gDNA by PCR as previously described (n = 98 patients; ref. 5). All the exons of CXCR4 gene were analyzed using targeted next-generation sequencing (NGS) in a cohort of 53 patients. NGS was analyzed using the Ion Torrent PGM platform (Life Biotechnologies). Amplicons covering the regions of interest were designed with an amplicon length of 150 to 250 bp. Libraries were sequenced with 200-bp read length on a 318 chip. Mutation calling was performed using the torrent suite variant caller under the low stringency somatic settings (TS4.0). We have studied CXCR4 mutations along with MYD88 L265P (exon 5), CD79A (ITAM domain), CD79B (ITAM domain), CARD11 (exons 5–9), N-RAS (exons 2 and 3), K-RAS (exons 2 and 3), BRAF6 (exon 15), and PTEN (exon 5+7) using NGS (n = 53).

To quantify the mutated clone, we have analyzed the VAF (variant allele frequency) defined as the number of reads that mapped each studied to this position, cover the variant base and show the reference allele, divided by all fragments covering the site. VAF was corrected by the percentage of tumoral cells in B cell selected sample. Mutation calls were considered positive when called by at least 20 variants reads. Complete sequence data of the coding and splice site regions of CXCR4 were generated at a mean depth coverage of 2,000× per nucleotide. NGS assay also allowed a better assessment of the percentage of CXCR4 mutant allele whose sensitivity is of 1%.

Gene expression profiling

Gene expression profiling (GEP) was performed using U133A arrays (Affymetrix) for seven of the patients with Waldenström macroglobulinemia. Total RNA was extracted from purified B tumoral cells population isolated from bone marrow using the TRIzol method. Expression data were normalized using the Robust Multi-array Average (RMA) algorithm. Differential gene expression was analyzed with the Bioconductor R package “limma,” which is well-known to improve variance modeling when the number of replicates is small. All probes with an adjusted P value under 0.05 were considered differentially expressed. Their normalized expression values were then represented in a heatmap, whose colors vary according to the standardized by row values of the input data. In parallel, a gene set enrichment analysis (GSEA) was performed for several gene sets among them the C2 collection of the MSigDB (Supplementary Tables S1 and S2). The PGSEA package (17) was used to calculate z scores for all gene sets, and P values were adjusted for multiple testing (18).

Cytogenetic analysis, FISH, SNPa, and immunophenotypic studies

Genome-wide detection of copy number alteration (CNA) and LOH was performed using the Genome-Wide Human SNP Array 6.0 (Affymetrix). Conventional cytogenetic analysis was performed on DSP30 + IL2 stimulated bone marrow cells. FISH was performed to detect chromosomal aberrations: 6q23 deletion, 17p12, 11q22 deletions, 13q14 deletion, trisomy 12, trisomy 4 (n = 90; ref. 15). Immunophenotypic expression of CXCR4 (R et D systems, UK), CD49d, CD27, CD80, CD86, CD138 gating on CD19+ cells was performed along with CD38 (Immunotech) and determination of the Matutes's score by flow cytometry. The ratio of mean fluorescence intensity (MFIR) was calculated as the ratio of specific fluorescence with an isotype control.

The relationships between the clinical, biologic, and molecular parameters were determined using a nonparametric test (Mann–Whitney), a t test, a χ2, or Fisher exact tests when appropriate. Correlations were tested through Spearman correlation coefficient. The differences between the results of the comparative tests were considered statistically significant if P < 0.05. All statistical analyses were done with the SPSS 15.0 software. Survival was studied using Kaplan–Meier test, comparisons were made with log-rank test

New CXCR4 mutations identified in Waldenström macroglobulinemia

We have first used NGS (deep sequencing) to screen the entire sequence of CXCR4 (n = 53 Waldenström macroglobulinemia samples) with a greater sensitivity and to quantify the allelic frequency of the variant of CXCR4 mutations in B selected cells. We have identified 14 of 53 (26.4%) Waldenström macroglobulinemia with CXCR4 mutations. We then confirmed the observed mutations using NGS using Sanger sequencing (SaS) and extend to the entire cohort of 98 patients with Waldenström macroglobulinemia (including the 53 done by NGS). We identified CXCR4 mutations (CXCR4mut) for a total of 24 of 98 (24.5%), and we have characterized several new mutations across the entire coding exons of CXCR4 (Fig. 1).

Figure 1.

Somatic mutations in C terminus of CXCR4 identified by Sanger sequencing and targeted NGS in patients with Waldenström macroglobulinemia (WM). List of CXCR4 mutations observed in our study. n, number of patients. The mutations previously reported in WM by Treon et al. are noted with *. In the schematic structure of CXCR4, sites of mutation are highlighted in blue (frameshift mutation) or gray (nonsense mutation). The S338 hotspot was represented in black.

Figure 1.

Somatic mutations in C terminus of CXCR4 identified by Sanger sequencing and targeted NGS in patients with Waldenström macroglobulinemia (WM). List of CXCR4 mutations observed in our study. n, number of patients. The mutations previously reported in WM by Treon et al. are noted with *. In the schematic structure of CXCR4, sites of mutation are highlighted in blue (frameshift mutation) or gray (nonsense mutation). The S338 hotspot was represented in black.

Close modal

Overall, 17 different mutations were identified in CXCR4, including 12 that we have described for the first time in Waldenström macroglobulinemia. Interestingly, only one type of CXCR4 mutation was observed in either patient, and all mutations identified were located in the 45 amino acid intracytoplasmic carboxy-terminal tail. These mutations were never observed in the paired T lymphocytes and thus confirmed their somatic feature. All mutations were heterozygous, and no UPDa (acquired uniparental disomy, LOH without variation of copy number) was observed at CXCR4 locus nor variation of copy gene number (gain or deletion) in our cohort, using SNP array.

High allelic frequency of CXCR4 mutation in Waldenström macroglobulinemia

NGS also allowed the quantification the allelic frequency of the variant of CXCR4 mutations. This quantification allows to conclude for the percentage of mutation in the clonal population for a given patient and thus to differentiate between clonal (dominant clonal) and subclonal. The mutation load of CXCR4 varied from 13.5% to 47.58% (mean, 35.2%) using the VAF (n = 14, CXCR4mut Waldenström macroglobulinemia in our series). Interestingly, for 12 patients with Waldenström macroglobulinemia, the high VAF of CXCR4 mutation was suggested to be present in the dominant clone. However, we were also able to identify the presence of CXCR4 mutations at the subclonal level in 4 patients using NGS, and for these 4 patients, the VAF varied from 13.5% to 26% (Fig. 2).

Figure 2.

Representation of VAF of CXCR4 and MYD88 mutation in each patient using NGS. Each row represents a patient (x-axis). VAF of MYD88 L265P mutation is represented in light gray, CXCR4 VAF in black. Percentage indicates allele frequencies determined by NGS corrected by the percentage of tumoral cells in the B selected samples in each 53 patients (y-axis). *, subclonal CXCR4 mutation in patients with Waldenström macroglobulinemia.

Figure 2.

Representation of VAF of CXCR4 and MYD88 mutation in each patient using NGS. Each row represents a patient (x-axis). VAF of MYD88 L265P mutation is represented in light gray, CXCR4 VAF in black. Percentage indicates allele frequencies determined by NGS corrected by the percentage of tumoral cells in the B selected samples in each 53 patients (y-axis). *, subclonal CXCR4 mutation in patients with Waldenström macroglobulinemia.

Close modal

Nonsense (CXCR4NS) and frameshift (CXCR4FS) CXCR4 mutations lead to truncation of CXCR4 C-terminal

Among the CXCR4 mutations observed, 11 of 24 (45.8%) and 13 of 24 (54.2%) were nonsense (CXCR4NS) and frameshift (CXCR4FS) mutations, respectively. The most frequent mutation was the C1013G (S338X) mutation (5 of 98, 5.1%) followed by C1013A (S338X; 3 of 98, 3%). The modelization of C1013G/CXCR4 and C1013A/CXCR4 variants predicts a stop codon in place of a serine at amino acid position 338 (S338X). We also found several other mutations to S338, suggesting a hotspot locus in CXCR4 gene. The S338 hotspot mutation was targeted in 50% of mutated cases in Waldenström macroglobulinemia. Other CXCR4NS mutations were here described. Interestingly, the exact same acquired R334X mutation in Waldenström macroglobulinemia was also previously described in the constitutional WHIM syndrome (19).

Increased expression of CXCR4 in Waldenström macroglobulinemia with CXCR4mut

Because we observed that all mutations were located in the C-terminal tail of CXCR4 and that this C-terminal domain plays a critical role in the desensitization process that contributes to the regulation of the signal transduction and CXCR4 expression (2); we then sought to study the impact of CXCR4 mutation on the expression of CXCR4 on the membrane of Waldenström macroglobulinemia tumor cells using flow cytometry. CXCR4 was expressed in all cases (n = 53), including the 12 CXCR4mut patients. The pattern varied widely using the MFIR to evaluate the CXCR4 expression (MFI varied from 4.6 to 148.7). However, there was a significant greater expression of CXCR4 protein in Waldenström macroglobulinemia with CXCR4mut (P = 0.003; Fig. 3). This increased expression was observed independently of the type of mutation (data not shown).

Figure 3.

Expression of CXCR4 and CD49d according to MYD88 and CXCR4 mutation status. Boxplots are shown with an overlay of the individual data points. ns, not significant.

Figure 3.

Expression of CXCR4 and CD49d according to MYD88 and CXCR4 mutation status. Boxplots are shown with an overlay of the individual data points. ns, not significant.

Close modal

Phenotypic signature associated to CXCR4mut genotype

CXCR4 directly interacts with CD49d in response to CXCL12 signaling in regulating migration and adhesion of Waldenström macroglobulinemia cells to the bone marrow microenvironment (10). CD49d favors lymphocyte homing by cooperating with CXCR4 in stromal cell adhesion and extracellular matrix (9). We then studied the impact of CXCR4 mutation on the expression of the integrin VLA4 (CD49d) using flow cytometry (n = 53). CD49d was expressed in all Waldenström macroglobulinemia cases; but no difference in CD49d expression was observed according to CXCR4 mutation. In contrast, a higher expression of CD49d was observed in Waldenström macroglobulinemia with MYD88 L265P mutation (P = 0.048; Fig. 3). Similarly, we found no difference in CD38 expression according to CXCR4 mutation, despite of the known functional link between CXCR4, CD49d, and CD38 (20). CD138 mediates cell adhesion to the extracellular matrix and is a marker of plasmacytic differentiation (21). Interestingly, we found a significant lower expression of CD138 (syndecan1) in the CXCR4mut group (P = 0.034). We found no difference in the CD27, CD23, CD10, K/L expression and the Matutes' score between CXCR4mut and CXCR4wild subgroups.

Genomic signature of CXCR4mut Waldenström macroglobulinemia

We then sought to study whether a specific high-throughput SNPa signature was associated to CXCR4mut Waldenström macroglobulinemia in a series of 53 Waldenström macroglobulinemia including 12 CXCR4mut. For this analysis, we have considered “presence of one” as any genomic abnormality including gain, loss, and/or copy number without LOH (CN-LOH), identified using SNPa. In this first analysis, we have compared presence versus absence of genomic abnormality according to presence of CXCR4 mutation. We found a relationship between CXCR4mut and a greater frequency of genomic aberrations in Waldenström macroglobulinemia compared with CXCR4wild (91% vs. 68%, P = 0.010).

We then evaluated the number of genomic abnormality as a reflection of the genomic complexity on the basis of the number of SNP abnormalities. We found that CXCR4-mutated Waldenström macroglobulinemia had a greater mean number of abnormality (5.8 vs. 2.8 per patient, P = 0.046). No significant difference was observed between Waldenström macroglobulinemia with clonal or subclonal CXCR4 mutation.

Furthermore, we analyzed the relationship between CXCR4 mutation and certain type of SNP aberration. We have observed that Waldenström macroglobulinemia with CXCR4mut had a greater incidence rate of trisomy 4 (complete or partial), 58% versus 12% respectively (P = 0.002); a greater frequency of gain of Xq, including MYD88 pathway–based IRAK1 gene, and of deletion 8p (P = 0.002 and P = 0.007, respectively; Fig. 4). Finally, we also found a greater frequency of deletion 6q, including the NF-κB pathway–based TNFAIP3 gene (P = 0.038). We found no relationship to deletion 7q, deletion 11q, deletion 17p, deletion 13q14, trisomy 18 in our study. No difference of SNPa feature was observed according to the type of CXCR4 mutation (nonsense (CXCR4NS) versus frameshift (CXCR4FS) mutations.

Figure 4.

Association network of genomic alterations in Waldenström macroglobulinemia (WM). The heatmap represents co-occurrence of CXCR4 mutation with other genetic alterations analyzed in 53 patients with WM using SNPa, cytogenetic, and FISH. The co-occurrence of CXCR4 mutation with other genetic alterations was represented using MeV software (TM4 microarray software suite; refs. 31, 32). Each row corresponds to key locus or regulator genes of the canonical and noncanonical NF-κB pathways targeted by CNA, UPD, or mutation in WM. The columns represent individual patients color coded on the base of gene status. Each patient is represented by a virtual column (black, mutation; CN LOH, deletion or gain of gene; gray, wild-type).

Figure 4.

Association network of genomic alterations in Waldenström macroglobulinemia (WM). The heatmap represents co-occurrence of CXCR4 mutation with other genetic alterations analyzed in 53 patients with WM using SNPa, cytogenetic, and FISH. The co-occurrence of CXCR4 mutation with other genetic alterations was represented using MeV software (TM4 microarray software suite; refs. 31, 32). Each row corresponds to key locus or regulator genes of the canonical and noncanonical NF-κB pathways targeted by CNA, UPD, or mutation in WM. The columns represent individual patients color coded on the base of gene status. Each patient is represented by a virtual column (black, mutation; CN LOH, deletion or gain of gene; gray, wild-type).

Close modal

Mutational landscape of CXCR4mut Waldenström macroglobulinemia

We next attempted to delineate the mutational spectrum and intratumoral heterogeneity related to presence of CXCR4mut for genes involved in the RAS (this pathway may be activated by CXCR4/protein G proteins), BCR and TLR pathways (MYD88 exon 5), known to be interconnected with CXCR4 pathway, using targeted NGS (6).

When looking at the incidence rate of the CXCR4 mutations with presence of MYD88 L265P mutation in the cohort of 98 patients, it appeared that 23 of 24 (95.8%) of the CXCR4mut were MYD88 L265P mutated, except one Waldenström macroglobulinemia with a CXCR4FS mutation. The high correlation between MYD88 L265P and CXCR4 mutational loads tends to suggest the coexistence of the 2 mutations in the same clone (Fig. 2). However, in four patients with Waldenström macroglobulinemia, a lower VAF of CXCR4 in comparison to MYD88 L265P suggested a subclonal CXCR4 mutation in the dominant clone harboring MYD88 mutation.

Regarding the presence of mutations in the RAS pathway in Waldenström macroglobulinemia, including PTEN, K-Ras, and N-Ras, crucial regulators of RAS pathway, we have not identified any mutation in our cohort (n = 53).

BCR pathway can be constitutively activated, either by signals from the microenvironment or by genetic aberration (22). Several mutations of key signaling molecules and regulators of the BCR pathway were reported overtime, particularly in the activated B-cell–like (ABC) diffuse large B-cell lymphoma (DLBCL) and in Waldenström macroglobulinemia (15, 23), including the immunoreceptor tyrosine–based activation motif (ITAM) signaling modules of the CD79A and CD79B BCR coreceptors and the coiled-coil domain of CARD11 (22). Mutations of CD79B and CD79A were observed in 12 of 98 (12.2%) of Waldenström macroglobulinemia. Interestingly, we found no coexistence of CXCR4mut and CD79A/CD79B mutations, but in one patient who had CXCR4mut at the subclonal level. No mutation of CARD11 was observed in our studied group (n = 53). One might conclude that there is no unique somatic mutation in the BCR pathway, but instead a wide variety of molecules that could be altered through molecular alterations in the BCR pathway, likely mutually exclusive in the same clone, such as CXCR4mut and CD79A/CD79B mutations in the present situation in our study.

Transcriptional signature related to presence of CXCR4 mutation

We finally studied the gene expression profile of Waldenström macroglobulinemia on 7 patients with Waldenström macroglobulinemia to understand whether the CXCR4 mutation impacted the gene expression and to identify a potential gene signature of Waldenström macroglobulinemia with CXCR4mut. The differential analysis allowed to identify a gene signature attached to the CXCR4mut molecular profile, corresponding to 32 probes and 27 genes. We found several genes among these, such as MMP8, CRISP3, or BRCC3 involved in cell death and survival or cell movement (ARG1, LYST). These probes were all underexpressed in the five CXCR4wild Waldenström macroglobulinemia compared with the two CXCR4mut Waldenström macroglobulinemia (Fig. 5 and Supplementary Tables S1 and S2). We then sought to identify deregulated pathways using Ingenuity analysis software. We found CXCR4mut to be related to several functional network networks such as cell cycle and DNA replication or cell-to-cell signaling and interaction, cellular growth, and proliferation. Using GSEA for targeted gene sets previously associated with CXCR4 mutation (12), we also identified deregulation of proteasome pathway and BCR/B lymphocyte pathway in CXCR4mut Waldenström macroglobulinemia (P = 0.01 and P = 0.04, respectively; Supplementary Table S1).

Figure 5.

Gene expression profiling. Heatmap of normalized expression values for differentially expressed probes. Colors vary according to the standardized by row values of the input data. The scale of the gene expression data is indicated. Some genes of interests are highlighted.

Figure 5.

Gene expression profiling. Heatmap of normalized expression values for differentially expressed probes. Colors vary according to the standardized by row values of the input data. The scale of the gene expression data is indicated. Some genes of interests are highlighted.

Close modal

Clinical and biologic features of CXCR4mut

We thought to identify clinical–biologic characteristics of Waldenström macroglobulinemia according to CXCR4mut features. The median age was 67 years (range, 36–92 years) in our series. We found that CXCR4mut was associated to a higher IgM monoclonal component (r = 0.309; P = 0.006) and thrombocytopenia (r = −0.226; P = 0.048), markers of adverse prognosis in Waldenström macroglobulinemia. Indeed, 35% and 12% of patients had serum IgM level greater than 30 g/L (P = 0.023) and 27% and 6.5% of patients had low platelet count lower than 100 g/L (P = 0.018), respectively, in patients with CXCR4mut versus CXCR4wild. There was a trend for CXCR4mut patients to have more often IPSS Waldenström macroglobulinemia stage 3 (P = 0.086). No association was found with regard to other markers of tumoral syndrome, such as adenopathy and splenomegaly, extramedullary localization, or lymphocytosis (one would consider a marker of egression of tumoral cells from the bone marrow; Supplementary Table S3). With a median follow-up of 6 years, 21 (22%) patients in our series had died with a median OS (95% confidence interval) not reached and estimated at 63% at 10 years. No significant difference in type of treatment or number of treatment line was observed (data not shown). The estimated overall survival (OS) at 10 years according to CXCR4mut status was lower for CXCR4mut patients, 50% versus 65.5% for patients with CXCR4wild (P = ns, Fig. 6). Interestingly, the difference in OS between CXCR4mut versus CXCR4wild became significant in patients with indolent Waldenström macroglobulinemia (estimated 5-year OS rate at 50% and 93%, respectively, P = 0.019); the CXCR4mut patients required therapy more often and much earlier than the CXCR4wild.

Figure 6.

OS according to presence or absence of CXCR4 mutation in Waldenström macroglobulinemia. The solid line indicates the CXCR4wild patients, the dashed line CXCR4mut patients. ns, not significant.

Figure 6.

OS according to presence or absence of CXCR4 mutation in Waldenström macroglobulinemia. The solid line indicates the CXCR4wild patients, the dashed line CXCR4mut patients. ns, not significant.

Close modal

We have described 12 new CXCR4 mutations for a total of 17 patients with Waldenström macroglobulinemia so far, observed in 25% of patients with Waldenström macroglobulinemia, using targeted deep sequencing. These data suggest a highly heterogeneous pattern of CXCR4 mutations in Waldenström macroglobulinemia and confirms using a more sensitive technique that was reported by Treon and colleagues using Sanger (11). We have also confirmed CXCR4 1013C>G, followed by 1013C>A, mutations to be the most frequent CXCR4 mutations in Waldenström macroglobulinemia, although in a lesser frequency as previously reported, 12.5% (for the 2 mutations using Sanger sequencing in CD19 selected cells) and 30% (only 1013C>A using allele-specific PCR in unselected tumoral samples) in Treon and colleagues and Roccaro and colleagues, respectively (11–12). Overall the high frequency of CXCR4 Whim-like mutation is a new hallmark of Waldenström macroglobulinemia (12, 24).

Interestingly, all mutations, either nonsense or frameshift, occurred in the C-terminus of CXCR4. It was suggested that the truncation of the distal amino acid region of the C-terminus known to regulate signaling of CXCR4 by CXCL12, might deregulate the CXCR4/SDF1 axis signaling pathway, and thus participate in the pathogenesis of Waldenström macroglobulinemia (2, 11, 12, 25). One might propose that the type of mutation does not matter much in Waldenström macroglobulinemia, on the contrary to the loss of C-terminus domain of CXCR4 protein which modified function might play a role in Waldenström macroglobulinemia pathogenesis (2). Functional studies may confirm this hypothesis in the near future.

We found a greater CXCR4 expression on tumor cells of Waldenström macroglobulinemia carrying CXCR4 mutation, independently of the type of mutation, similar to that was described in WHIM syndrome. It is suggested that the increase of CXCR4 receptor expression would induce an altered migration profile in response to SDF1 (25), resulting in an increased egression of tumor cells as demonstrated in mouse model (12). In our study, we also found a decreased expression of CD138, a cell surface syndecan-1 that mediates cancer cell adhesion to the extracellular matrix, and alteration of expression of several genes involved in matrix degradation such as MMP8 or cell migration on gene expression profiling. Taking together, CXCR4 mutation might alter the ability of the tumor cells to interact with the bone marrow microenvironment in Waldenström macroglobulinemia, as shown in CXCR4 S338X cells (12).

Nearly all patients with Waldenström macroglobulinemia with CXCR4 mutations harbored the MYD88 L265P mutation, as initially described by Treon and colleagues suggesting a potential cooperation between these pathways (11). An important dataset obtained in our series came from the quantification of the allelic frequency of the variant (VAF) of CXCR4 mutations into tumoral cells using NGS that may identify minor subclones of less than 1% in our study. Targeted deep sequencing allows to conclude for the percentage of mutations in the clonal population for a given patient and thus to differentiate between clonal (dominant clonal) and subclonal mutations. Combining analysis of VAF of several genes, CXCR4 mutation was expressed in the same clone as MYD88 L265P, as we have identified a high allelic frequency of CXCR4 mutation MYD88 L265P in Waldenström macroglobulinemia tumor cells; but for 4 patients who had CXCR4 mutation present only at the subclonal level. Overall, these data may suggest the following hypothesis, where MYD88 L265P mutation is a founder mutation in Waldenström macroglobulinemia, a first genetic hit that would promote NF-κB and JAK/STAT3 signaling (4, 26). Direct inhibition of MYD88 L265P signaling overcomes CXCL12-trigered survival effects in CXCR4-mutated cells, supporting a primary role for MYD88 signaling in Waldenström macroglobulinemia (25). CXCR4 could then be one of the secondary events in some patients who showed subclonal mutation of CXCR4 compared with MYD88 L265P mutation present in the main clone. However, longitudinal analysis studies are needed to explore the dynamic of clonal architecture to identify driver or additional mutations in Waldenström macroglobulinemia that could contribute to clinical progression or chemoresistance. Indeed, it was shown that a diminished clinical activity of ibrutinib was observed in patients with Waldenström macroglobulinemia with CXCR4 mutations (25;30).

In addition, we also found absence of co-occurrence of CD79A/CD79B and CXCR4 mutations in Waldenström macroglobulinemia. These data may suggest that these types of activating mutations of BCR and CXCR4 pathway were mutually exclusive in MYD88 L265P Waldenström macroglobulinemia (6). A cooperation between CXCR4 and TLR signaling activated by MYD88L265P mutation was described in B cells (6). The functional role of mutations in the BCR pathway is not fully described in Waldenström macroglobulinemia, but CD79B mutation was shown to alter the BCR response in diffuse large B lymphoma (15, 23). We might thus suggest the existence of two subgroups of MYD88 L265P Waldenström macroglobulinemia, one with CXCR4 mutations and one with BCR pathway mutations that may involve in cooperation with the aforementioned driver mutation which primarily affect TLR signaling. These data further may advance the concept that a complex rather than a unique deregulation of the TLR/NF-κB pathway characterizes MYD88 L265P Waldenström macroglobulinemia with others interconnected pathways.

We have observed that CXCR4 mutation was more frequently associated to a complex genomic signature, including gain of chromosome 4, deletion 6q, and gain of chromosome X using SNPa. High-throughput genomic studies have identified multiple mechanisms of genetic changes in Waldenström macroglobulinemia including several recurrent CNA such as deletion 6q, deletion 13q, or gain of chromosome 4 and CN-LOH (15, 27). This genomic pattern segregates further Waldenström macroglobulinemia among the others B- cell tumors. The prognostic role of these genomic alterations is not fully understood, but the presence of more than 3 SNPa abnormalities was associated with symptomatic status (15, 27, 28). Taking together, this finding suggests an inter- and intraclonal heterogeneity and the potential selective advantage of specific combinations of genetic lesions in Waldenström macroglobulinemia, in particular in the subset of CXCR4mut Waldenström macroglobulinemia. The understanding of these genomic abnormalities, including gene mutations and CNA, described herein might help deciphering the pathogenesis of Waldenström macroglobulinemia to determine associated risk of progression, relapse, and drug resistance.

Overall, Waldenström macroglobulinemia with CXCR4 mutation had a specific clinicobiologic and genomic signature associated to features characterized with adverse prognosis in Waldenström macroglobulinemia, including high IgM M-spike and thrombocytopenia (part of the IPSS Waldenström macroglobulinemia score adverse features) and greater frequency of complex genomic aberrations pattern using SNP array. Our data therefore may suggest a worse prognosis for CXCR4 mutation subgroup of asymptomatic Waldenström macroglobulinemia, but we did not observe significant impact on overall survival in our cohort as previously described (11). Further studies are needed to explore the prognosis value of CXCR4 mutation in Waldenström macroglobulinemia in clinical trials. Interestingly, the potential clinical impact of Waldenström macroglobulinemia with CXCR4 mutation has already been suggested by Roccaro and colleagues and Treon and colleagues on the basis of their in vitro studies showing drug resistance in CXCR4-mutated Waldenström macroglobulinemia cells exposed to BTK, mTOR, and PI3K inhibitors but not proteasome inhibitors (29, 30). This study, along with other, thus led us to propose that it is primetime for a systematic evaluation of CXCR4 mutation as part as the pretreatment package of Waldenström macroglobulinemia, to identify patients with CXCR4 mutation that might benefit more from proteasome inhibitors as compared with other options, such a BTK inhibitors.

In conclusion, approximately 25% of Waldenström macroglobulinemia harbor CXCR4 mutation, leading to altered CXCR4 protein at the C-terminal end, irrespective of the pleiotropic pattern of mutations. The study of CXCR4 mutations showed existence of intraclonal (variation in co-expression of MYD88 and CXCR4 mutations) and interclonal (BCR and CXCR4 mutations in MYD88L265P Waldenström macroglobulinemia) heterogeneity when analyzing the molecular landscape of CXCR4 mutation in Waldenström macroglobulinemia. Patients with CXCR4 mutation displayed a specific clinical–biologic–genomic and transcriptomic signature, possibly related to specific physiopathologic mechanism related to CXCR4 mutation. This dataset further suggests the genomic heterogeneity of Waldenström macroglobulinemia, beyond the current indolent versus symptomatic classification.

No potential conflicts of interest were disclosed.

Conception and design: S. Poulain, M. Figeac, F. Lepretre, X. Leleu

Development of methodology: S. Poulain, A. Venet-Caillault, M. Figeac, S. Geffroy, X. Leleu

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): S. Poulain, C. Roumier, M. Figeac, C. Herbaux, G. Marot, E. Bertrand, A. Decambron, E.M. Boyle, A. Renneville, S. Tricot, B. Quesnel, X. Leleu

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): S. Poulain, C. Roumier, M. Figeac, G. Marot, E. Doye, B. Quesnel, X. Leleu

Writing, review, and/or revision of the manuscript: S. Poulain, C. Roumier, C. Preudhomme, X. Leleu

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): S. Poulain, C. Roumier, E. Doye, F. Lepretre, O. Nibourel, B. Quesnel, X. Leleu

Study supervision: S. Poulain, E. Doye, P. Duthilleul, X. Leleu

The authors thank Pauline Vandycke, Valérie Grandières, Claudine Delsaut, Axelle Séghir, and Sabine Ranwez for their excellent technical assistance.

This work was supported by the Comité Septentrion de la Ligue contre le Cancer et la Fondation Française pour la Recherche contre le Myélome et les Gammapathies (FFRMG).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Hunter
ZR
,
Xu
L
,
Yang
G
,
Zhou
Y
,
Liu
X
,
Cao
Y
, et al
The genomic landscape of Waldenström macroglobulinemia is characterized by highly recurring MYD88 and WHIM-like CXCR4 mutations, and small somatic deletions associated with B-cell lymphomagenesis
.
Blood
2014
;
123
:
1637
46
.
2.
Al Ustwani
O
,
Kurzrock
R
,
Wetzler
M
. 
Genetics on a WHIM
.
Br J Haematol
2014
Jan;
164
:
15
23
.
3.
Hernandez
PA
,
Gorlin
RJ
,
Lukens
JN
,
Taniuchi
S
,
Bohinjec
J
,
Francois
F
, et al
Mutations in the chemokine receptor gene CXCR4 are associated with WHIM syndrome, a combined immunodeficiency disease
.
Nat Genet
2003
;
34
:
70
4
.
4.
Treon
SP
,
Xu
L
,
Yang
G
,
Zhou
Y
,
Liu
X
,
Cao
Y
, et al
MYD88 L265P somatic mutation in Waldenstrom's macroglobulinemia
.
N Engl J Med
2012
;
367
:
826
33
.
5.
Poulain
S
,
Boyle
EM
,
Roumier
C
,
Demarquette
H
,
Wemeau
M
,
Geffroy
S
, et al
MYD88 L265P mutation contributes to the diagnosis of Bing Neel syndrome
.
Br J Haematol
2014
;
167
:
506
13
.
6.
Chatterjee
S
,
Behnam Azad
B
,
Nimmagadda
S
. 
The intricate role of CXCR4 in cancer
.
Adv Cancer Re
2014
;
124
:
31
82
.
7.
Conley-LaComb
MK
,
Saliganan
A
,
Kandagatla
P
,
Chen
YQ
,
Cher
ML
,
Chinni
SR
. 
PTEN loss mediated Akt activation promotes prostate tumor growth and metastasis via CXCL12/CXCR4 signaling
.
Mol Cancer
2013
;
12
:
85
.
8.
Nie
Y
,
Waite
J
,
Brewer
F
,
Sunshine
MJ
,
Littman
DR
,
Zou
YR
. 
The role of CXCR4 in maintaining peripheral B cell compartments and humoral immunity
.
J Exp Med
2004
;
200
:
1145
56
.
9.
Burger
JA
,
Gribben
JG
. 
The microenvironment in chronic lymphocytic leukemia (CLL) and other B cell malignancies: insight into disease biology and new targeted therapies
.
Semin Cancer Bio
2014
;
24
:
71
81
.
10.
Ngo
HT
,
Leleu
X
,
Lee
J
,
Jia
X
,
Melhem
M
,
Runnels
J
, et al
SDF-1/CXCR4 and VLA-4 interaction regulates homing in Waldenstrom macroglobulinemia
.
Blood
2008
;
112
:
150
8
.
11.
Treon
SP
,
Cao
Y
,
Xu
L
,
Yang
G
,
Liu
X
,
Hunter
ZR
. 
Somatic mutations in MYD88 and CXCR4 are determinants of clinical presentation and overall survival in Waldenstrom macroglobulinemia
.
Blood
2014
;
123
:
2791
6
.
12.
Roccaro
AM
,
Sacco
A
,
Jimenez
C
,
Maiso
P
,
Moschetta
M
,
Mishima
Y
, et al
C1013G/CXCR4 acts as a driver mutation of tumor progression and modulator of drug resistance in lymphoplasmacytic lymphoma
.
Blood
2014
,
123
:
4120
31
13.
Owen
RG
,
Pratt
G
,
Auer
RL
,
Flatley
R
,
Kyriakou
C
,
Lunn
MP
, et al
Guidelines on the diagnosis and management of Waldenstrom macroglobulinaemia
.
Br J Haemato
2014
;
165
:
316
33
.
14.
Dimopoulos
MA
,
Kastritis
E
,
Owen
RG
,
Kyle
RA
,
Landgren
O
,
Morra
E
, et al
Treatment recommendations for patients with Waldenstrom macroglobulinemia (WM) and related disorders: IWWM-7 consensus
.
Blood
2014
;
124
:
1404
11
.
15.
Poulain
S
,
Roumier
C
,
Galiegue-Zouitina
S
,
Daudignon
A
,
Herbaux
C
,
Aiijou
R
, et al
Genome wide SNP array identified multiple mechanisms of genetic changes in Waldenstrom macroglobulinemia
.
Am J Hematol
2013
;
88
:
948
54
.
16.
Poulain
S
,
Roumier
C
,
Decambron
A
,
Renneville
A
,
Herbaux
C
,
Bertrand
E
, et al
MYD88 L265P mutation in Waldenstrom macroglobulinemia
.
Blood
2013
;
121
:
4504
11
.
17.
Kim
SY
,
Volsky
DJ
. 
PAGE: parametric analysis of gene set enrichment
.
BMC Bioinformatics
2005
;
6
:
144
.
18.
Hochberg
Ba
. 
Controlling the false discovery rate: a practical and powerful approach to multiple testing
.
J Royal Stat Soc, Series B 57
1995
:
289
300
19.
Taniuchi
S
,
Masuda
M
,
Fujii
Y
,
Izawa
K
,
Kanegane
H
,
Kobayashi
Y
. 
The role of a mutation of the CXCR4 gene in WHIM syndrome
.
Haematologica
2005
;
90
:
1271
2
.
20.
Dal Bo
M
,
Tissino
E
,
Benedetti
D
,
Caldana
C
,
Bomben
R
,
Del Poeta
G
, et al
Microenvironmental interactions in chronic lymphocytic leukemia: the master role of CD49d
.
Semin Hematol
2014
;
51
:
168
76
.
21.
Teng
YH
,
Aquino
RS
,
Park
PW
. 
Molecular functions of syndecan-1 in disease
.
Matrix Biol
2012
;
31
:
3
16
.
22.
Lim
KH
,
Yang
Y
,
Staudt
LM
. 
Pathogenetic importance and therapeutic implications of NF-kappaB in lymphoid malignancies
.
Immunol Rev
2012
;
246
:
359
78
.
23.
Davis
RE
,
Ngo
VN
,
Lenz
G
,
Tolar
P
,
Young
RM
,
Romesser
PB
, et al
Chronic active B-cell-receptor signalling in diffuse large B-cell lymphoma
.
Nature
2010
;
463
:
88
92
.
24.
Martinez
N
,
Almaraz
C
,
Vaque
JP
,
Varela
I
,
Derdak
S
,
Beltran
S
, et al
Whole-exome sequencing in splenic marginal zone lymphoma reveals mutations in genes involved in marginal zone differentiation
.
Leukemia
2014
;
28
:
1334
40
.
25.
Cao
Y
,
Hunter
ZR
,
Liu
X
,
Xu
L
,
Yang
G
,
Chen
J
, et al
CXCR4 WHIM-like frameshift and nonsense mutations promote ibrutinib resistance but do not supplant MYD88(L265P) -directed survival signalling in Waldenstrom macroglobulinaemia cells
.
Br J Haematol
2015
;
168
:
701
7
.
26.
Wang
JQ
,
Jeelall
YS
,
Beutler
B
,
Horikawa
K
,
Goodnow
CC
. 
Consequences of the recurrent MYD88(L265P) somatic mutation for B cell tolerance
.
J Exp Med
2014
;
211
:
413
26
.
27.
Braggio
E
,
Philipsborn
C
,
Novak
A
,
Hodge
L
,
Ansell
S
,
Fonseca
R
. 
Molecular pathogenesis of Waldenstrom's macroglobulinemia
.
Haematologica
2012
;
97
:
1281
90
.
28.
Nguyen-Khac
F
,
Lambert
J
,
Chapiro
E
,
Grelier
A
,
Mould
S
,
Barin
C
, et al
Chromosomal aberrations and their prognostic value in a series of 174 untreated patients with Waldenstrom's macroglobulinemia
.
Haematologica
2013
;
98
:
649
54
.
29.
Roccaro
AM
,
Sacco
A
,
Jimenez
C
,
Maiso
P
,
Moschetta
M
,
Mishima
Y
, et al
C1013G/CXCR4 acts as a driver mutation of tumor progression and modulator of drug resistance in lymphoplasmacytic lymphoma
.
Blood
2014
;
123
:
4120
31
.
30.
Cao
Y
,
Hunter
ZR
,
Liu
X
,
Xu
L
,
Yang
G
,
Chen
J
, et al
The WHIM-like CXCR4(S338X) somatic mutation activates AKT and ERK, and promotes resistance to ibrutinib and other agents used in the treatment of Waldenstrom's Macroglobulinemia
.
Leukemia
2015
;
29
:
169
76
.
31.
Saeed
AI
,
Bhagabati
NK
,
Braisted
JC
,
Liang
W
,
Sharov
V
,
Howe
EA
, et al
TM4 microarray software suite
.
Methods Enzymol
2006
;
411
:
134
93
.
32.
Saeed
AI
,
Sharov
V
,
White
J
,
Li
J
,
Liang
W
,
Bhagabati
N
, et al
TM4: a free, open-source system for microarray data management and analysis
.
Biotechniques
2003
;
34
:
374
8
.