Multi-omics analysis in IgM monoclonal gammopathies reveals epigenetic influence on oncogenesis via DNA methylation.

Currently, the role of DNA methylation in the IgM-monoclonal gammopathy disease spectrum remains poorly understood. In the present study, a multi-omics prospective analysis was conducted integrating DNA methylation, RNA-seq and WES data in 34 subjects [23 WM, 6 IgM-MGUS, 5 normal controls]. Analysis was focused on defining differences between IgM-gammopathies (WM/IgM-MGUS) compared to controls, and specifically between WM and IgM-MGUS. Between groups, genome-wide DNA methylation analysis demonstrated a significant number of differentially methylated regions which were annotated according to genomic region. Next, integration of RNA-seq data was performed to identify potentially epigenetically deregulated pathways. We found that pathways involved in cell cycle, metabolism, cytokine/immune signaling, cytoskeleton, tumor microenvironment, and intracellular signaling were differentially activated and potentially epigenetically regulated. Importantly, there was a positive enrichment of CXCR4 signaling pathway along with several interleukin (IL-6, IL-8, IL15) signaling pathways in WM compared to IgM-MGUS. Further assessment of known tumor suppressor genes and oncogenes uncovered differential promoter methylation of several targets with concordant change in gene expression, including CCND1 and CD79B. Overall, this report defines how aberrant DNA methylation in IgM-gammopathies may play a critical role in the epigenetic control of oncogenesis and key cellular functions.

MicroRNA and long non-coding RNA analysis in IgM monoclonal gammopathies reveals epigenetic influence on cellular functions and oncogenesis.

Not available.

Molecular classification and identification of an aggressive signature in low-grade B-cell lymphomas.

Non-follicular low-grade B-cell lymphomas (LGBCL) are biologically diverse entities that share clinical and histologic features that make definitive pathologic categorization challenging. While most patients with LGBCL have an indolent course, some experience aggressive disease, highlighting additional heterogeneity across these subtypes. To investigate the potential for shared biology across subtypes, we performed RNA sequencing and applied machine learning approaches that identified five clusters of patients that grouped independently of subtype. One cluster was characterized by inferior outcome, upregulation of cell cycle genes, and increased tumor immune cell content. Integration of whole exome sequencing identified novel LGBCL mutations and enrichment of TNFAIP3 and BCL2 alterations in the poor survival cluster. Building on this, we further refined a transcriptomic signature associated with early clinical failure in two independent cohorts. Taken together, this study identifies unique clusters of LGBCL defined by novel gene expression signatures and immune profiles associated with outcome across diagnostic subtypes.

Targeting glycogen synthase kinase 3 for therapeutic benefit in lymphoma.

Targeting the B-cell receptor and phosphatidylinositol 3-kinase/mTOR signaling pathways has shown meaningful, but incomplete, antitumor activity in lymphoma. Glycogen synthase kinase 3 (GSK3) α and ß are 2 homologous and functionally overlapping serine/threonine kinases that phosphorylate multiple protein substrates in several key signaling pathways. To date, no agent targeting GSK3 has been approved for lymphoma therapy. We show that lymphoma cells abundantly express GSK3α and GSK3ß compared with normal B and T lymphocytes at the messenger RNA and protein levels. Utilizing a new GSK3 inhibitor 9-ING-41 and by genetic deletion of GSK3α and GSK3ß genes using CRISPR/CAS9 knockout, GSK3 was demonstrated to be functionally important to lymphoma cell growth and proliferation. GSK3ß binds to centrosomes and microtubules, and lymphoma cells treated with 9-ING-41 become arrested in mitotic prophase, supporting the notion that GSK3ß is necessary for the progression of mitosis. By analyzing recently published RNA sequencing data on 234 diffuse large B-cell lymphoma patients, we found that higher expression of GSK3α or GSK3ß correlates well with shorter overall survival. These data provide rationale for testing GSK3 inhibitors in lymphoma patient trials.

Genetic overlap between autoimmune diseases and non-Hodgkin lymphoma subtypes.

Epidemiologic studies show an increased risk of non-Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen-stimulation in one disease leads to another. Here we assess shared genetic risk in genome-wide-association-studies (GWAS). Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B-cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta-analysis. Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta-analysis (combing AD with NHL) implicated genes of apoptosis and telomere length. This GWAS-based analysis four NHL subtypes and three ADs revealed few weakly-associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.

Ibrutinib monotherapy outside of clinical trial setting in Waldenström macroglobulinaemia: practice patterns, toxicities and outcomes.

Ibrutinib-related data in Waldenström macroglobulinaemia (WM) remain sparse, particularly outside of trials. We report on 80 patients [previously treated, n = 67 (84%), treatment-naïve, n = 13 (16%)] with WM, evaluated consecutively at Mayo Clinic, who received ibrutinib off-study after its approval in 2015 for WM. Overall response rate (ORR) was 91%; major-response rate (MRR) was 78%. The median time to first response and best response was 2·9 [95% confidence interval (CI): 2-4] and 5·7 (95% CI: 4-12) months, respectively. The median follow-up was 19 (95% CI: 14-21) months; 18-month progression-free survival (PFS) was 82%. The median time on therapy was 12·5 (95% CI: 9·3-16·7) months, and the median duration-of-response was 32 (range: 23-32) months. Twenty-five patients (31%) had discontinued therapy at last follow-up (68% due to treatment-related toxicities) and 18% of patients required dose reduction. Fatigue (12%) and atrial-fibrillation (11%) were common non-haematological toxicities. IgM rebound occurred in 36% of patients who abruptly discontinued ibrutinib. Following ibrutinib discontinuation, 84% of patients received subsequent treatment, achieving an ORR of 57% and MRR of 50%. The median PFS from commencement of subsequent salvage therapy was 18 months. Ibrutinib therapy, outside of clinical trials, is effective in WM, but is associated with toxicities and challenges, including IgM rebound and a high drug discontinuation rate for reasons other than disease progression.

Impact of MYD88L265P mutation status on histological transformation of Waldenström Macroglobulinemia.

Histological transformation in Waldenström macroglobulinemia (WM) is an uncommon complication, with limited data, particularly regarding the impact of MYD88 L265P mutation on transformation. We examined risk factors and outcomes associated with transformation in WM, highlighting the role of MYD88 L265P mutation. Patients with WM seen at Mayo Clinic, Rochester, USA and University Hospital of Reims, France, between 01/01/1996 and December 31, 2017 were included; 50 (4.3%) of 1147 patients transformed to a high-grade lymphoma, with median time-to-transformation of 4.5 (range 0-21) years in the transformed cohort. The MYD88 L265P mutation status was known in 435/1147 (38%) patients (406 with non-transformed WM and 29 patients in transformed cohort). On multivariate analysis, MYD88 WT status alone was an independent predictor of transformation (odds ratio, 7[95%CI: 2.1-23]; P = .003). Additionally, the MYD88 WT status was independently associated with shorter time-to-transformation (HR 7.9 [95%CI: 2.3-27; P = .001]), with a 5-year transformation rate of 16% for MYD88 WT vs 2.8% with MYD88 L265P mutated patients. Patients with transformation demonstrated a significant increase in risk of death compared to patients who did not transform (HR 5.075; 95%CI: 3.8-6.8; P < .001). In conclusion, the MYD88 WT status is an independent predictor of transformation and associated with a shorter time-to-transformation. Additionally, transformation conferred an inferior overall survival in patients with WM.

The utility of prognostic indices, early events, and histological subtypes on predicting outcomes in non-follicular indolent B-cell lymphomas.

Indolent B-cell lymphomas other than follicular lymphoma account for up to 10% of all B-cell neoplasms. While they typically follow a slowly progressive course, some patients experience rapid progression and early mortality. Prognostic scoring systems have not been adopted, hindering the ability of clinicians or researchers to predict outcomes, or risk-stratify patients during clinical trials. To address this, we evaluated the utility of existing prognostic indices and novel, early disease-related outcomes, to predict subsequent long term survival. Baseline characteristics and outcomes data were generated from a longitudinal cohort study that prospectively enrolled 632 patients newly diagnosed with marginal zone lymphoma, lymphoplasmacytic lymphomas, or B-cell lymphomas not otherwise specified, beginning in 2002. The International Prognostic Index (IPI), Follicular Lymphoma International Prognostic Index (FLIPI), and MALT International prognostic index (MALT-IPI) demonstrated c-statistics that ranged from 0.593-0.612 for event-free survival (EFS), and 0.683-0.714 for overall survival (OS). Patients who attained event-free survival at 12 months (EFS12) experienced similar mortality to the US general population (standardized mortality ratio [SMR] 1.19; 95% CI 0.95-1.46). Patients who did not attain EFS12 had subsequent worse morality (SMR 3.14 (95% CI 2.05-4.59). The MALT-IPI demonstrated utility in predicting subsequent long-term outcomes among patients with non-follicular indolent B-cell lymphomas. This index should be used by clinicians giving guidance to patients at the time of initial diagnosis, and risk stratification during clinical studies. The divergent long-term outcomes experienced by patients who do or do not attain EFS12 suggest there exists a subset of patients who harbor high-risk disease. Future research efforts should focus on methods to identify these patients at the time of diagnosis, in order to enable risk-tailored therapy.

Genetically predicted longer telomere length is associated with increased risk of B-cell lymphoma subtypes.

Evidence from a small number of studies suggests that longer telomere length measured in peripheral leukocytes is associated with an increased risk of non-Hodgkin lymphoma (NHL). However, these studies may be biased by reverse causation, confounded by unmeasured environmental exposures and might miss time points for which prospective telomere measurement would best reveal a relationship between telomere length and NHL risk. We performed an analysis of genetically inferred telomere length and NHL risk in a study of 10 102 NHL cases of the four most common B-cell histologic types and 9562 controls using a genetic risk score (GRS) comprising nine telomere length-associated single-nucleotide polymorphisms. This approach uses existing genotype data and estimates telomere length by weighing the number of telomere length-associated variant alleles an individual carries with the published change in kb of telomere length. The analysis of the telomere length GRS resulted in an association between longer telomere length and increased NHL risk [four B-cell histologic types combined; odds ratio (OR) = 1.49, 95% CI 1.22-1.82,P-value = 8.5 × 10(-5)]. Subtype-specific analyses indicated that chronic lymphocytic leukemia or small lymphocytic lymphoma (CLL/SLL) was the principal NHL subtype contributing to this association (OR = 2.60, 95% CI 1.93-3.51,P-value = 4.0 × 10(-10)). Significant interactions were observed across strata of sex for CLL/SLL and marginal zone lymphoma subtypes as well as age for the follicular lymphoma subtype. Our results indicate that a genetic background that favors longer telomere length may increase NHL risk, particularly risk of CLL/SLL, and are consistent with earlier studies relating longer telomere length with increased NHL risk.

Human Pegivirus infection and lymphoma risk and prognosis: a North American study.

We evaluated the association of Human Pegivirus (HPgV) viraemia with risk of developing lymphoma, overall and by major subtypes. Because this virus has also been associated with better prognosis in the setting of co-infection with human immunodeficiency virus, we further assessed the association of HPgV with prognosis. We used risk factor data and banked plasma samples from 2094 lymphoma cases newly diagnosed between 2002 and 2009 and 1572 frequency-matched controls. Plasma samples were tested for HPgV RNA by reverse transcription polymerase chain reaction (RT-PCR), and those with RNA concentrations <5000 genome equivalents/ml were confirmed using nested RT-PCR methods. To assess the role of HPgV in lymphoma prognosis, we used 2948 cases from a cohort study of newly diagnosed lymphoma patients (included all cases from the case-control study). There was a positive association of HPgV viraemia with risk of lymphoma overall (Odds ratio = 2·14; 95% confidence interval [CI] 1·63-2·80; P < 0·0001), and for all major subtypes except Hodgkin lymphoma and chronic lymphocytic leukaemia/small lymphocytic lymphoma, and this was not confounded by other lymphoma risk factors. In contrast, there was no association of HPgV viraemia with event-free survival (Hazard ratio [HR] = 1·00; 95% CI 0·85-1·18) or overall survival (HR = 0·97; 95% CI 0·79-1·20) for lymphoma overall, or any of the subtypes. These data support the hypothesis for a role of HPgV in the aetiology of multiple lymphoma subtypes.

MYD88 mutation status does not impact overall survival in Waldenström macroglobulinemia.

Waldenström macroglobulinemia (WM) is an immunoglobulin M-associated lymphoma, with majority of cases demonstrating MYD88 locus alteration, most commonly, MYD88L265P . Owing to low prevalence of the wild-type (WT) MYD88 genotype in WM, clinically relevant data in this patient population are sparse, with one study showing nearly a 10-fold increased risk of mortality in this subgroup compared to patients with MYD88L265P mutation. We studied a large cohort of patients with MYD88L265P and MYD88WT WM, evaluated at Mayo Clinic, Rochester, between 1995 and 2016, to specifically assess the impact of these genotypes on clinical course. Of 557 patients, MYD88L265P mutation status, as determined by allele-specific polymerase chain reaction, was known in 219, and 174 (79%) of those exhibited MYD88L265P , 157 of 174 patients had active disease. Of 45 (21%) patients with MYD88WT genotype, 44 had active disease. The estimated median follow-up was 7.0 years; median overall survival was 10.2 years (95% CI: 8.4-16.5) for MYD88L265P versus 13.9 years (95% CI: 6.4-29.3) for the MYD88WT (P = 0.86). The time-to-next therapy from frontline treatment and the presenting features were similar in the two patient populations. For patients with smoldering WM at diagnosis, the median time-to-progression to active disease was 2.8 years (95% CI: 2.2-3.8) in the MYD88L265P cohort and 1.9 years (95% CI: 0.7-3.1) in the MYD88WT cohort (P = 0.21). The frequency of transformation to high-grade lymphoma, or the development of therapy-elated myelodysplastic syndrome was higher in the MYD88WT cohort (16% versus 4% in the MYD88L265P , P = 0.009). In conclusion, MYD88L265P mutation does not appear to be a determinant of outcome, and its presence may not be a disease-defining feature in WM. Our findings warrant external validation, preferably through prospective studies.

BCL2 mutations are associated with increased risk of transformation and shortened survival in follicular lymphoma.

Follicular lymphoma (FL), an indolent neoplasm caused by a t(14;18) chromosomal translocation that juxtaposes the BCL2 gene and immunoglobulin locus, has a variable clinical course and frequently undergoes transformation to an aggressive lymphoma. Although BCL2 mutations have been previously described, their relationship to FL progression remains unclear. In this study, we evaluated the frequency and nature of BCL2 mutations in 2 independent cohorts of grade 1 and 2 FLs, along with the correlation between BCL2 mutations, transformation risk, and survival. The prevalence of BCL2 coding sequence mutations was 12% in FL at diagnosis and 53% at transformation (P < .0001). The presence of these BCL2 mutations at diagnosis correlated with an increased risk of transformation (hazard ratio 3.6; 95% CI, 2.0-6.2; P < .0001) and increased risk of death due to lymphoma (median survival of 9.5 years with BCL2 mutations vs 20.4 years without; P = .012). In a multivariate analysis, BCL2 mutations and high FL international prognostic index were independent risk factors for transformation and death due to lymphoma. Some mutant Bcl-2 proteins exhibited enhanced antiapoptotic capacity in vitro. Accordingly, BCL2 mutations can affect antiapoptotic Bcl-2 function, are associated with increased activation-induced cytidine deaminase expression, and correlate with increased risk of transformation and death due to lymphoma.

FCGR3A/2A polymorphisms and diffuse large B-cell lymphoma outcome treated with immunochemotherapy: a meta-analysis on 1134 patients from two prospective cohorts.

Single nucleotide polymorphisms (SNPs) in FCγ-receptor genes FCGR3A (rs396991) and FCGR2A (rs1801274) influence the affinity of the Fc portion of anti-CD20 immunoglobulin G1 monoclonal antibody. Their roles in diffuse large B-cell lymphoma (DLBCL) treated with rituximab in combination with anthracycline-based chemotherapy remain controversial. To address this question, we genotyped FCGR2A and FCGR3A SNPs in two prospective DLBCL cohorts from Lymphoma Study Association trials (N = 554) and Iowa/Mayo Specialized Program Of Research Excellence (N = 580). Correlations with treatment response and hematological toxicity were assessed in Lymphoma Study Association. Correlation with event-free survival (EFS) and overall survival (OS) was performed in both cohorts, followed by a meta-analysis to increase power. Our study shows the absence of correlation between these SNPs and treatment response. Grades 3 and 4 febrile neutropenia during treatment was more frequently observed in FCGR3A VV (39%) than VF (29%) and FF (32%) carriers (p = 0.04). Our analysis for EFS and OS shows that FCGR3A was not associated with outcome. In a meta-analysis using an ordinal model, FCGR2A (per R allele) was associated with a better EFS (hazard ratio = 0.87; 95%CI, 0.76-0.99; p = 0.04) and OS (hazard ratio = 0.86; 95%CI, 0.73-1.00; p = 0.05) which was not altered after adjustment for the International Prognostic Index. Overall, our data demonstrate that patients with DLBCL with the low-affinity FCγRIIA RR had an unexpectedly better outcome than FCγRIIA H carriers. Whether rituximab efficacy is improved in FCγRIIA RR patients due a clearance reduction or other functions of FCγRIIA in DLBCL should be investigated (clinicaltrials.gov identifiers: NCT00135499, NTC00135499 NCT00140595, NCT00144807, NCT00144755, NCT01087424, and NCT00301821). Copyright © 2016 John Wiley & Sons, Ltd.

Constitutive activation of STAT5A and STAT5B regulates IgM secretion in Waldenstrom's macroglobulinemia.

Activation of the Janus kinase family/signal transducer and activator of transcription (JAK/STAT) signaling pathway has been associated with the pathogenesis and progression of both solid and hematologic malignancies. We have detected constitutive activation of STAT5 in malignant B cells derived from patients with Waldenström's macroglobulinemia (WM). Although short hairpin RNA-mediated knockdown of the STAT5A and STAT5B isoforms did not affect cellular proliferation, loss of STAT5 significantly decreased immunoglobulin M (IgM) secretion. A similar dose-dependent inhibition of IgM secretion was observed when WM cell lines were treated with a small molecule inhibitor of STAT5. These data suggest that STAT5 is involved in regulating IgM production in WM and that inhibition of STAT5 may represent a novel therapeutic strategy for lowering IgM levels in WM patients.

Targeted inhibition of the deubiquitinating enzymes, USP14 and UCHL5, induces proteotoxic stress and apoptosis in Waldenström macroglobulinaemia tumour cells.

Deubiquitinase enzymes (DUBs) of the proteasomal 19S regulatory particle are emerging as important therapeutic targets in several malignancies. Here we demonstrate that inhibition of two proteasome-associated DUBs (USP14 and UCHL5) with the small molecule DUB inhibitor b-AP15, results in apoptosis of human Waldenström macroglobulinaemia (WM) cell lines and primary patient-derived WM tumour cells. Importantly, b-AP15 produced proteotoxic stress and apoptosis in WM cells that have acquired resistance to the proteasome inhibitor bortezomib. In silico modelling identified protein residues that were critical for the binding of b-AP15 with USP14 or UCHL5 and proteasome enzyme activity assays confirmed that b-AP15 does not affect the proteolytic capabilities of the 20S proteasome ß-subunits. In vitro toxicity from b-AP15 appeared to result from a build-up of ubiquitinated proteins and activation of the endoplasmic reticulum stress response in WM cells, an effect that also disrupted the mitochondria. Focused transcriptome profiling of b-AP15-treated WM cells revealed modulation of several genes regulating cell stress and NF-κB signalling, the latter whose protein translocation and downstream target activation was reduced by b-AP15 in vitro. This is the first report to define the effects and underlying mechanisms associated with inhibition of USP14 and UCHL5 DUB activity in WM tumour cells.

RVboost: RNA-seq variants prioritization using a boosting method.

MOTIVATION: RNA-seq has become the method of choice to quantify genes and exons, discover novel transcripts and detect fusion genes. However, reliable variant identification from RNA-seq data remains challenging because of the complexities of the transcriptome, the challenges of accurately mapping exon boundary spanning reads and the bias introduced during the sequencing library preparation. METHOD: We developed RVboost, a novel method specific for RNA variant prioritization. RVboost uses several attributes unique in the process of RNA library preparation, sequencing and RNA-seq data analyses. It uses a boosting method to train a model of 'good quality' variants using common variants from HapMap, and prioritizes and calls the RNA variants based on the trained model. We packaged RVboost in a comprehensive workflow, which integrates tools of variant calling, annotation and filtering. RESULTS: RVboost consistently outperforms the variant quality score recalibration from the Genome Analysis Tool Kit and the RNA-seq variant-calling pipeline SNPiR in 12 RNA-seq samples using ground-truth variants from paired exome sequencing data. Several RNA-seq-specific attributes were identified as critical to differentiate true and false variants, including the distance of the variant positions to exon boundaries, and the percent of the reads supporting the variant in the first six base pairs. The latter identifies false variants introduced by the random hexamer priming during the library construction. AVAILABILITY AND IMPLEMENTATION: The RVboost package is implemented to readily run in Mac or Linux environments. The software and user manual are available at http://bioinformaticstools.mayo.edu/research/rvboost/.

PatternCNV: a versatile tool for detecting copy number changes from exome sequencing data.

MOTIVATION: Exome sequencing (exome-seq) data, which are typically used for calling exonic mutations, have also been utilized in detecting DNA copy number variations (CNVs). Despite the existence of several CNV detection tools, there is still a great need for a sensitive and an accurate CNV-calling algorithm with built-in QC steps, and does not require a paired reference for each sample. RESULTS: We developed a novel method named PatternCNV, which (i) accounts for the read coverage variations between exons while leveraging the consistencies of this variability across different samples; (ii) reduces alignment BAM files to WIG format and therefore greatly accelerates computation; (iii) incorporates multiple QC measures designed to identify outlier samples and batch effects; and (iv) provides a variety of visualization options including chromosome, gene and exon-level views of CNVs, along with a tabular summarization of the exon-level CNVs. Compared with other CNV-calling algorithms using data from a lymphoma exome-seq study, PatternCNV has higher sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION: The software for PatternCNV is implemented using Perl and R, and can be used in Mac or Linux environments. Software and user manual are available at http://bioinformaticstools.mayo.edu/research/patterncnv/, and R package at https://github.com/topsoil/patternCNV/.

Discovery and prioritization of somatic mutations in diffuse large B-cell lymphoma (DLBCL) by whole-exome sequencing.

To gain insight into the genomic basis of diffuse large B-cell lymphoma (DLBCL), we performed massively parallel whole-exome sequencing of 55 primary tumor samples from patients with DLBCL and matched normal tissue. We identified recurrent mutations in genes that are well known to be functionally relevant in DLBCL, including MYD88, CARD11, EZH2, and CREBBP. We also identified somatic mutations in genes for which a functional role in DLBCL has not been previously suspected. These genes include MEF2B, MLL2, BTG1, GNA13, ACTB, P2RY8, PCLO, and TNFRSF14. Further, we show that BCL2 mutations commonly occur in patients with BCL2/IgH rearrangements as a result of somatic hypermutation normally occurring at the IgH locus. The BCL2 point mutations are primarily synonymous, and likely caused by activation-induced cytidine deaminase-mediated somatic hypermutation, as shown by comprehensive analysis of enrichment of mutations in WRCY target motifs. Those nonsynonymous mutations that are observed tend to be found outside of the functionally important BH domains of the protein, suggesting that strong negative selection against BCL2 loss-of-function mutations is at play. Last, by using an algorithm designed to identify likely functionally relevant but infrequent mutations, we identify KRAS, BRAF, and NOTCH1 as likely drivers of DLBCL pathogenesis in some patients. Our data provide an unbiased view of the landscape of mutations in DLBCL, and this in turn may point toward new therapeutic strategies for the disease.

IL-21 in the bone marrow microenvironment contributes to IgM secretion and proliferation of malignant cells in Waldenstrom macroglobulinemia.

Cytokines within the tumor microenvironment play an important role in supporting the growth and survival of B-cell malignancies. One such cytokine, IL-21, promotes the growth of myeloma and Hodgkin lymphoma cells while inducing apoptosis in chronic lymphocytic leukemia. However, the biologic significance of IL-21 has not been examined in Waldenstrom macroglobulinemia (WM), a B-cell lymphoma characterized by elevated serum IgM and a lymphoplasmacytic bone marrow infiltrate. We report here on the presence of IL-21 in the bone marrow of patients with WM and have identified activated T cells as the source of this cytokine. We readily detected the IL-21 receptor on malignant WM B cells and show that IL-21 significantly increases both IgM secretion and cellular proliferation of these cells with no effect on viability. IL-21 rapidly induces phosphorylation of STAT3 in WM cells, and treatment of the WM cell line MWCL-1 with a STAT3 inhibitor abolished the IL-21-mediated increases in cellular proliferation and IgM secretion. IL-21 also increased the expression of known STAT3 targets involved in B-cell differentiation, including BLIMP-1, XBP-1, IL-6, and IL-10. Overall, our data indicate that IL-21 in the bone marrow microenvironment significantly affects the biology of WM tumor cells through a STAT3-dependent mechanism.