The lymph node microenvironment promotes B-cell receptor signaling, NF-kappaB activation, and tumor proliferation in chronic lymphocytic leukemia

Abstract

Chronic lymphocytic leukemia (CLL), an incurable malignancy of mature B lymphocytes, involves blood, bone marrow, and secondary lymphoid organs such as the lymph nodes (LN). A role of the tissue microenvironment in the pathogenesis of CLL is hypothesized based on in vitro observations, but its contribution in vivo remains ill-defined. To elucidate the effects of tumor-host interactions in vivo, we purified tumor cells from 24 treatment-naive patients. Samples were obtained concurrently from blood, bone marrow, and/or LN and analyzed by gene expression profiling. We identified the LN as a key site in CLL pathogenesis. CLL cells in the LN showed up-regulation of gene signatures, indicating B-cell receptor (BCR) and nuclear factor-κB activation. Consistent with antigen-dependent BCR signaling and canonical nuclear factor-κB activation, we detected phosphorylation of SYK and IκBα, respectively. Expression of BCR target genes was stronger in clinically more aggressive CLL, indicating more effective BCR signaling in this subtype in vivo. Tumor proliferation, quantified by the expression of the E2F and c-MYC target genes and verified with Ki67 staining by flow cytometry, was highest in the LN and was correlated with clinical disease progression. These data identify the disruption of tumor microenvironment interactions and the inhibition of BCR signaling as promising therapeutic strategies in CLL. This study is registered at https://ichgcp.net/clinical-trials-registry/NCT00019370" title="See in ClinicalTrials.gov">NCT00019370.

Figures

Figure 1

The tissue microenvironment affects the tumor biology of CLL cells in vivo. (A) Gene expression analysis of fresh CLL tumor cells derived from PB, BM, and LN in treatment-naive patients in whom all 3 sites were evaluable (n = 12). Principal component analysis of gene expression normalized for patient effect (see “Microarray hybridization, data processing, and analysis”) is shown. (B) Heat map of 151 genes differentially expressed in purified CLL cells obtained simultaneously from LN and PB in 17 patients (> 2-fold change, FDR < 20%). Patient samples are arranged in columns keeping the same order in the LN and PB groups. Gene expression is median-centered and scaled as indicated. Gene symbols highlight select genes. (C) Heat map of 26 differentially expressed genes in purified CLL cells isolated from BM aspirates compared with cells simultaneously obtained from the PB in 19 patients (> 2-fold change, FDR < 20%).

Figure 2

Dynamic changes in CXCR4 and CD69 expression on CLL cells in different anatomic compartments. (A) CXCR4 was quantified by flow cytometry on CD19+-gated cells and expressed as the percentage of cells expressing CXCR4 above isotype control. Samples obtained simultaneously from the same patient are connected by a line. Mononuclear cells were obtained by density gradient centrifugation, maintained on ice, and analyzed within 24 hours. Normal B cells constituted less than 1% of CD19+ cells. Comparison was by paired t test. (B) CD69 was quantified by flow cytometry on CD19+-gated cells and expressed as the percentage of cells expressing CD69 above isotype control. Sample collection and analysis were as in panel A.

Figure 3

Activation of the B-cell receptor on CLL cells induces a characteristic gene expression signature and phosphorylation of SYK. (A) Gene expression changes in UM-CLL (n = 4) and M-CLL (n = 4) samples 6 hours after in vitro IgM cross-linking. Sixty-one genes showing a > 2-fold change between unstimulated and stimulated cells (P < .001) were identified and constitute a “CLL-BCR gene signature.” Gene expression changes in response to BCR engagement is depicted in a heat map according to the scale shown. Genes are sorted top to bottom for the ratio between the average up-regulation in UM-CLL compared with M-CLL cells. (B) Phosphorylation of SYK (p-SYK) after in vitro IgM cross-linking was assessed by flow cytometry in cells gated on CD3−SYK+IgM+ at the indicated time points. (C) Mean and SD of the percentage of pSYK-positive cells in CLL subtypes (UM-CLL n = 14, M-CLL n = 13) and normal B cells (control n = 11) are shown. Differences between groups was assessed by Student t test (P = .013 for UM-CLL vs M-CLL at 45 minutes; all other comparisons were not significant).

Figure 4

BCR activation in LN-resident CLL cells. (A) Genes in the CLL-BCR signature that were most significantly enriched in the LN as identified by GSEA (leading-edge genes) are depicted in a heat map for the 12 patients having contributed cells from all 3 compartments. (B) The BCR score was computed as the average of the mRNA expression level of the leading edge genes for each sample. Shown is the ratio of the BCR score in BM (n = 19) or LN (n = 17) relative to the score of the matched PB sample. Comparison between BM and LN was by Student t test. The sample with exceptionally high score in BM was obtained from CLL_C10. (C, D) Relative gene expression changes of BCR-regulated genes between UM-CLL and M-CLL. Each gene is represented by a bar, select genes are highlighted. (C) Ratio of relative gene expression change after IgM cross-linking between UM-CLL (n = 4) and M-CLL (n = 4). Bars represent 61 genes, as described in Figure 3A. The fold up-regulation of each gene was averaged for each subtype and the ratio between UM-CLL and M-CLL is shown on a log2 scale. (D) Ratio of relative gene expression change in UM-CLL (n = 12) and M-CLL (n = 5) between the LN- and matched PB-derived cells. Bars represent all 61 genes as in panel C. The fold up-regulation of each gene in the LN-derived cells compared with the matched PB sample was averaged within each subtype, and the ratio between UM-CLL and M-CLL is shown on a log2 scale. (E) Phosphorylation of SYK was assessed by flow cytometry in CD3− cells from the indicated anatomical site. Shown are 2 representative patients of a total of 6 analyzed; in 5 patients, staining for pSYK was increased in LN- compared with PB-derived cells (median 7-fold increase).

Figure 5

NF-κB is activated through the canonical pathway in LN-derived CLL cells. (A) Genes in the NF-κB signature that were most significantly enriched in the LN as identified by GSEA (leading edge genes) are depicted in a heat map for the 12 patients having contributed cells from all 3 compartments. (B) The NF-κB score was computed as the average of the mRNA expression level of the leading edge genes for each sample. Shown is the ratio of the BCR score in BM (n = 19) or LN (n = 17) relative to the score of the matched PB sample. Comparison between BM and LN was by Student t test. The sample with the exceptionally high score in the BM was obtained from CLL_C10. (C) Western blot of nuclear lysates from purified CLL cells; JUNB expression was normalized to TBP and the ratio between LN- and PB-derived cells is shown. Western blots of cytoplasmic protein fractions from purified CLL cells were probed for phosphorylated IκBα (p-IκBα, Ser 32-36; D) and for total IκBα (E). IκBα and p-IκBα expression quantified by densitometry was normalized to γ-tubulin and the ratio between LN- and PB-derived cells is shown. (F) Pearson correlation between IκBα expression and the NF-κB gene expression score.

Figure 6

Tumor proliferation and c-MYC activation in the tissue microenvironment. (A) Genes of the E2F and c-MYC signatures that were most significantly enriched in the LN as identified by GSEA (leading-edge genes) are depicted in a heat map for the 12 patients having contributed cells from all 3 compartments. For E2F only genes that were at least 1.5-fold more highly expressed in LN than PB are shown. (B) The E2F and c-MYC scores were computed individually as the average of the mRNA expression level of the leading-edge genes for each sample. Shown is the ratio of the E2F and c-MYC scores in BM (n = 19) or LN (n = 17) relative to the score of the matched PB sample. Comparison between BM and LN was by Student t test. (C) Western blot of nuclear protein extracts from purified CLL cells were probed with anti-E2F1 or anti–c-MYC antibody, quantified by densitometry, and normalized to TBP. The ratios of LN- to PB-derived cells for E2F1 and c-MYC protein expression are indicated.

Figure 7

Tumor proliferation in the LN correlates with disease progression. (A) Magnification (×200) of LN biopsies of 2 representative CLL samples: top panel, preferential Ki67 staining in proliferation centers; bottom panel, diffuse Ki67 positivity. (B) Ki67 expression was assessed by flow cytometry in CD3− cells from the indicated anatomical sites. Shown are 2 representative patients of a total of 6 analyzed. Cutoffs in each sample were chosen so that PB cells were around 1% Ki67+ cells. (C) Ki67 expression in CD3− cells from the indicated anatomic sites. Statistical comparison was by paired Student t test. (D) Kaplan-Meier analysis of time to treatment in all patients with LN biopsies (n = 17) based on E2F scores: dashed line, scores above the median; solid line, scores below the median.