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.

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.

Follicular lymphoma (FL) is an indolent B cell malignancy characterized by an extensive but poorly functional T cell infiltrate in the tumor microenvironment. Using mass cytometry, we identified at least 12 subsets of intratumoral CD4+ T cells, 3 of which were unique to FL biopsies versus control tissues. Of these subsets, the frequency of naive T cells correlated with improved patient survival. Although total PD-1+ T cell numbers were not associated with patient outcome, specific PD-1+ T cell subpopulations were associated with poor survival. Intratumoral T cells lacking CD27 and CD28 co-stimulatory receptor expression were enriched in FL and correlated with inferior patient outcomes. In vitro models revealed that CD70+ lymphoma cells played an important role in expanding this population. Taken together, our mass cytometry results identified CD4+ memory T cell populations that are poorly functional due to loss of co-stimulatory receptor expression and are associated with an inferior survival in FL.

Exhausted T-cells in follicular lymphoma (FL) typically express PD-1, but expression of PD-1 is not limited to exhausted cells. Although expected to be functionally suppressed, we found that the population of intratumoral PD-1+ T cells were predominantly responsible for production of cytokines and granules. This surprising finding prompted us to explore the involvement of LAG-3 to specifically identify functionally exhausted T cells. We found that LAG-3 was expressed on a subset of intratumoral T cells from FL and LAG-3+ T cells almost exclusively came from PD-1+ population. CyTOF analysis revealed that intratumoral LAG-3+ T cells were phenotypically heterogeneous as LAG-3 was expressed on a variety of T cell subsets. In contrast to PD-1+LAG-3- cells, intratumoral PD-1+LAG-3+ T cells exhibited reduced capacity to produce cytokines and granules. LAG-3 expression could be substantially upregulated on CD4+ or CD8+ T cells by IL-12, a cytokine that has been shown to induce T-cell exhaustion and be increased in the serum of lymphoma patients. Furthermore, we found that blockade of both PD-1 and LAG-3 signaling enhanced the function of intratumoral CD8+ T cells resulting in increased IFN-γ and IL-2 production. Clinically, LAG-3 expression on intratumoral T cells correlated with a poor outcome in FL patients. Taken together, we find that LAG-3 expression is necessary to identify the population of intratumoral PD-1+ T cells that are functionally exhausted and, in contrast, find that PD-1+LAG-3- T cells are simply activated cells that are immunologically functional. These findings may have important implications for immune checkpoint therapy in FL.

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.

Copy number alterations (CNAs) of 9p24.1 occur frequently in Hodgkin lymphoma, primary mediastinal large B-cell lymphoma (PMBCL), primary central nervous system lymphoma, and primary testicular lymphoma, resulting in overexpression of PD-L1 and sensitivity to PD-1 blockade-based immunotherapy. While 9p24.1 CNA was also reported in diffuse large B-cell lymphoma (DLBCL), little is known about its molecular or clinical significance. In this study, we analyzed the prevalence of 9p24.1 CNA in newly diagnosed DLBCL and examined its association with PD-L1, PD-L2, and JAK2 expression, clinical characteristics, and outcome. We found that 10% of DLBCL cases had CNA of 9p24.1, with 6.5% gains, and 3.5% amplifications. Only the cases with a 9p24.1 amplification had high levels of PD-L1, PD-L2, and JAK2 expression. Gains or amplifications of 9p24.1 were associated with a younger age and the ABC/non-GCB subtype. Compared with DLBCL cases without 9p24.1 CNA, the cases with a 9p24.1 amplification had a trend of better event-free survival. Furthermore, the amplification cases had a gene expression and mutation profile similar to those of PMBCL. Our data suggest that amplification of 9p24.1 identifies a unique subset of DLBCL with clinical and molecular features resembling PMBCL that may be amenable to PD-1 blockade-based immunotherapy.

Marginal zone lymphoma (MZL) is the third most common subtype of B-cell non-Hodgkin lymphoma. Here we perform a two-stage GWAS of 1,281 MZL cases and 7,127 controls of European ancestry and identify two independent loci near BTNL2 (rs9461741, P=3.95 × 10(-15)) and HLA-B (rs2922994, P=2.43 × 10(-9)) in the HLA region significantly associated with MZL risk. This is the first evidence that genetic variation in the major histocompatibility complex influences MZL susceptibility.

Non-Hodgkin lymphoma (NHL) represents a diverse group of hematological malignancies, of which follicular lymphoma (FL) is a prevalent subtype. A previous genome-wide association study has established a marker, rs10484561 in the human leukocyte antigen (HLA) class II region on 6p21.32 associated with increased FL risk. Here, in a three-stage genome-wide association study, starting with a genome-wide scan of 379 FL cases and 791 controls followed by validation in 1,049 cases and 5,790 controls, we identified a second independent FL-associated locus on 6p21.32, rs2647012 (OR(combined)  = 0.64, P(combined)  = 2 × 10(-21)) located 962 bp away from rs10484561 (r(2)<0.1 in controls). After mutual adjustment, the associations at the two SNPs remained genome-wide significant (rs2647012:OR(adjusted)  = 0.70, P(adjusted)  =  4 × 10(-12); rs10484561:OR(adjusted)  = 1.64, P(adjusted)  = 5 × 10(-15)). Haplotype and coalescence analyses indicated that rs2647012 arose on an evolutionarily distinct haplotype from that of rs10484561 and tags a novel allele with an opposite (protective) effect on FL risk. Moreover, in a follow-up analysis of the top 6 FL-associated SNPs in 4,449 cases of other NHL subtypes, rs10484561 was associated with risk of diffuse large B-cell lymphoma (OR(combined)  = 1.36, P(combined)  =  1.4 × 10(-7)). Our results reveal the presence of allelic heterogeneity within the HLA class II region influencing FL susceptibility and indicate a possible shared genetic etiology with diffuse large B-cell lymphoma. These findings suggest that the HLA class II region plays a complex yet important role in NHL.

CD8+ T-lymphocyte subsets defined by killer lectin-like receptor G1 (KLRG1) and CD127 expression have been reported to have an important role in infection, but their role in the setting of lymphoid malignancies, specifically follicular lymphoma (FL), has not been studied.

Treatment with programmed death-1 (PD-1) blocking antibodies results in high overall response rates in refractory and relapsed classical Hodgkin lymphoma (cHL) patients, indicating that PD-1/PD-1 ligand interactions are integral to progression of this disease. Given the genetically driven increased PD-L1/2 expression in HL, we hypothesized that reverse signaling through PD-1 ligands may be a potential mechanism contributing to the growth and survival of Hodgkin Reed-Sternberg (HRS) cells in cHL. Our data show that engagement of PD-L1 using an agonistic monoclonal antibody increases cell survival and proliferation and reduces apoptosis in HL cell lines. We show that HL patients have significantly higher serum levels of soluble PD-1 than healthy controls, and find that both membrane-bound and soluble forms of PD-1 are able to induce PD-L1 reverse signaling in HL cell lines. PD-L1 signaling, which is associated with activation of the MAPK pathway and increased mitochondrial oxygen consumption, is reversed by PD-1 blockade. In summary, our data identify inhibition of reverse signaling through PD-L1 as an additional mechanism that accounts for clinical responses to PD-1 blockade in cHL.

Diffuse large B cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of 3 new genome-wide association studies (GWAS) and 1 previous scan, totaling 3,857 cases and 7,666 controls of European ancestry, with additional genotyping of 9 promising SNPs in 1,359 cases and 4,557 controls. In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; P = 2.33 × 10(-21)), rs2523607 at 6p21.33 (HLA-B; P = 2.40 × 10(-10)), rs79480871 at 2p23.3 (NCOA1; P = 4.23 × 10(-8)) and two independent SNPs, rs13255292 and rs4733601, at 8q24.21 (PVT1; P = 9.98 × 10(-13) and 3.63 × 10(-11), respectively). These data provide substantial new evidence for genetic susceptibility to this B cell malignancy and point to pathways involved in immune recognition and immune function in the pathogenesis of DLBCL.

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.

Chronic lymphocytic leukemia (CLL) is a common lymphoid malignancy with strong heritability. To further understand the genetic susceptibility for CLL and identify common loci associated with risk, we conducted a meta-analysis of four genome-wide association studies (GWAS) composed of 3,100 cases and 7,667 controls with follow-up replication in 1,958 cases and 5,530 controls. Here we report three new loci at 3p24.1 (rs9880772, EOMES, P=2.55 × 10(-11)), 6p25.2 (rs73718779, SERPINB6, P=1.97 × 10(-8)) and 3q28 (rs9815073, LPP, P=3.62 × 10(-8)), as well as a new independent SNP at the known 2q13 locus (rs9308731, BCL2L11, P=1.00 × 10(-11)) in the combined analysis. We find suggestive evidence (P<5 × 10(-7)) for two additional new loci at 4q24 (rs10028805, BANK1, P=7.19 × 10(-8)) and 3p22.2 (rs1274963, CSRNP1, P=2.12 × 10(-7)). Pathway analyses of new and known CLL loci consistently show a strong role for apoptosis, providing further evidence for the importance of this biological pathway in CLL susceptibility.

Metabolic reprogramming is a hallmark of cancer cells. In Waldenstrom Macroglobulinemia (WM), the infiltration of IgM-secreting lymphoplamacytic cells into the bone marrow (BM) could shift the homeostasis of proteins and metabolites towards a permissive niche for tumor growth. Here, we investigated whether alerted metabolic pathways contribute to the pathobiology of WM and whether the cytokine composition of the BM promotes such changes. Metabolomics analysis on WM patients and normal donors' serum samples revealed a total of 75 metabolites that were significantly altered between two groups. While these metabolites belonged to amino acids, glucose, glutathione and lipid metabolism pathways, the highest number of the differentially expressed metabolites belonged to glutathione metabolism. Proteomics analysis and immunohistochemical staining both confirmed the increased protein levels mediating glutathione metabolism, including GCLC, MT1X, QPCT and GPX3. Moreover, treatment with IL-6 and IL-21, cytokines that induce WM cell proliferation and IgM secretion, increased gene expression of the amino acid transporters mediating glutathione metabolism, including ASCT2, SLC7A11 and 4F2HC, indicating that cytokines in the WM BM could modulate glutathione metabolism. Glutathione synthesis inhibition using Buthionine sulphoximine (BSO) significantly reduced WM cells proliferation in vitro, accompanied with decreased NFκB-p65 and MAPK-p38 phosphorylation. Moreover, BSO treatment significantly reduced the tumor growth rate in a WM xenograft model, further highlighting the role of glutathione metabolism in promoting tumor growth and proliferation. In summary, our data highlight a central role for glutathione metabolism in WM pathobiology and indicate that intervening with the metabolic processes could be a potential therapy for WM patients.

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma, and front line therapies have not improved overall outcomes since the advent of immunochemotherapy. By pairing DNA and gene expression data with clinical response data, we identified a high-risk subset of non-GCB DLBCL patients characterized by genomic alterations and expression signatures capable of sustaining an inflammatory environment. These mutational alterations (PIM1, SPEN, and MYD88 [L265P]) and expression signatures (NF-κB, IRF4, and JAK-STAT engagement) were associated with proliferative signaling, and were found to be enriched in patients treated with RCHOP that experienced unfavorable outcomes. However, patients with these high-risk mutations had more favorable outcomes when the immunomodulatory agent lenalidomide was added to RCHOP (R2CHOP). We are the first to report the genomic validation of a high-risk phenotype with a preferential response towards R2CHOP therapy in non-GCB DLBCL patients. These conclusions could be translated to a clinical setting to identify the ~38% of non-GCB patients that could be considered high-risk, and would benefit from alternative therapies to standard RCHOP based on personalized genomic data.

How to identify follicular lymphoma (FL) patients with low disease burden but high risk for early progression is unclear. Building on a prior study demonstrating the early transformation of FLs with high variant allele frequency (VAF) BCL2 mutations at activation-induced cytidine deaminase (AICDA) sites, we examined 11 AICDA mutational targets, including BCL2, BCL6, PAX5, PIM1, RHOH, SOCS, and MYC, in 199 newly diagnosed grade 1 and 2 FLs. BCL2 mutations with VAF ≥20% occurred in 52% of cases. Among 97 FL patients who did not initially receive rituximab-containing therapy, nonsynonymous BCL2 mutations at VAF ≥20% were associated with increased transformation risk (HR 3.01, 95% CI 1.04-8.78, p = 0.043) and a trend toward shorter event-free survival (EFS, median 20 months with mutations versus 54 months without, p = 0.052). Other sequenced genes were less frequently mutated and did not increase the prognostic value of the panel. Across the entire population, nonsynonymous BCL2 mutations at VAF ≥20% were associated with decreased EFS (HR 1.55, 95% CI 1.02-2.35, p = 0.043 after correction for FLIPI and treatment) and decreased overall survival after median 14-year follow-up (HR 1.82, 95% CI 1.05-3.17, p = 0.034). Thus, high VAF nonsynonymous BCL2 mutations remain prognostic even in the chemoimmunotherapy era.

Waldenström macroglobulinemia (WM)/lymphoplasmacytic lymphoma (LPL) is a rare, chronic B-cell lymphoma with high heritability. We conduct a two-stage genome-wide association study of WM/LPL in 530 unrelated cases and 4362 controls of European ancestry and identify two high-risk loci associated with WM/LPL at 6p25.3 (rs116446171, near EXOC2 and IRF4; OR = 21.14, 95% CI: 14.40-31.03, P = 1.36 × 10-54) and 14q32.13 (rs117410836, near TCL1; OR = 4.90, 95% CI: 3.45-6.96, P = 8.75 × 10-19). Both risk alleles are observed at a low frequency among controls (~2-3%) and occur in excess in affected cases within families. In silico data suggest that rs116446171 may have functional importance, and in functional studies, we demonstrate increased reporter transcription and proliferation in cells transduced with the 6p25.3 risk allele. Although further studies are needed to fully elucidate underlying biological mechanisms, together these loci explain 4% of the familial risk and provide insights into genetic susceptibility to this malignancy.

MYD88 mutations are one of the most recurrent mutations in hematologic malignancies. However, recent mouse models suggest that MYD88L265P alone may not be sufficient to induce tumor formation. Interplay between MYD88L265P and other genetic events is further supported by the fact that TNFAIP3 (A20) inactivation often accompanies MYD88L265P. However, we are still lacking information about the consequence of MYD88L265P in combination with TNFAIP3 loss in human B cell lymphoma. Review of our genetic data on diffuse large B cell lymphoma (DLBCL) and Waldenstrom macroglobulinemia (WM), found that a large percentage of DLBCL and WM cases that have a MYD88 mutation also harbor a TNFAIP3 loss, 55% DLBCL and 28% of WM, respectively. To mimic this combination of genetic events, we used genomic editing technology to knock out TNFAIP3 in MYD88L265P non-Hodgkin's lymphoma (NHL) cell lines. Loss of A20 expression resulted in increased NF-κB and p38 activity leading to upregulation of the NF-κB target genes BCL2 and MYC. Furthermore, we detected the increased production of IL-6 and CXCL10 which led to an upregulation of the JAK/STAT pathway. Overall, these results suggest that MYD88L265P signaling can be enhanced by a second genetic alteration in TNFAIP3 and highlights a potential opportunity for therapeutic targeting.

While the effect of TGF-β on malignant B cells in non-Hodgkin lymphoma (NHL) has been previously evaluated, studies to specifically define the role of TGF-β in tumor immunity in B-cell NHL are limited. We found that soluble TGF-β, secreted by both lymphoma cells and intratumoral T cells, is present in the serum of patients with B-cell NHL. Soluble TGF-β promoted regulatory T (T(reg)) cells by enhancing expression of Foxp3 in CD4(+) T cells and suppressed effector helper T (T(H)) cells by inhibiting expression of IFN-γ and IL-17. Blockade of the IL-2 signaling pathway diminished the effect of soluble TGF-β on T cell differentiation. Furthermore, we found that membrane-bound TGF-β is expressed specifically on the surface of malignant B cells in B-cell NHL. TGF-β was able to bind to the surface of lymphoma B cells through an interaction with heparan sulfate (HS) but not through the TGF-β receptor. We showed that pretreatment of lymphoma B cells with TGF-β significantly inhibits the proliferation and cytokine production of intratumoral T cells. Taken together, these results suggest that tumor-associated soluble and membrane-bound TGF-β are involved in the regulation of intratumoral T cell differentiation and function in B-cell NHL.

Genome-wide association studies (GWAS) have previously identified 13 loci associated with risk of chronic lymphocytic leukemia or small lymphocytic lymphoma (CLL). To identify additional CLL susceptibility loci, we conducted the largest meta-analysis for CLL thus far, including four GWAS with a total of 3,100 individuals with CLL (cases) and 7,667 controls. In the meta-analysis, we identified ten independent associated SNPs in nine new loci at 10q23.31 (ACTA2 or FAS (ACTA2/FAS), P=1.22×10(-14)), 18q21.33 (BCL2, P=7.76×10(-11)), 11p15.5 (C11orf21, P=2.15×10(-10)), 4q25 (LEF1, P=4.24×10(-10)), 2q33.1 (CASP10 or CASP8 (CASP10/CASP8), P=2.50×10(-9)), 9p21.3 (CDKN2B-AS1, P=1.27×10(-8)), 18q21.32 (PMAIP1, P=2.51×10(-8)), 15q15.1 (BMF, P=2.71×10(-10)) and 2p22.2 (QPCT, P=1.68×10(-8)), as well as an independent signal at an established locus (2q13, ACOXL, P=2.08×10(-18)). We also found evidence for two additional promising loci below genome-wide significance at 8q22.3 (ODF1, P=5.40×10(-8)) and 5p15.33 (TERT, P=1.92×10(-7)). Although further studies are required, the proximity of several of these loci to genes involved in apoptosis suggests a plausible underlying biological mechanism.