Abstract Currently, the role of DNA methylation in the immunoglobulin M (IgM) monoclonal gammopathy disease spectrum remains poorly understood. In the present study, a multiomics prospective analysis was conducted integrating DNA methylation, RNA sequencing (RNA-seq), and whole-exome sequencing data in 34 subjects (23 with Waldenström macroglobulinemia [WM], 6 with IgM monoclonal gammopathy of undetermined significance [MGUS], and 5 normal controls). Analysis was focused on defining differences between IgM gammopathies (WM/IgM-MGUS) compared with controls, and specifically between WM and IgM-MGUS. Between groups, genome-wide DNA methylation analysis demonstrated a significant number of differentially methylated regions that 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 the CXCR4 signaling pathway along with several interleukin (interleukin 6 [IL-6], IL-8, and IL-15) signaling pathways in WM compared with 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.
Molecular classification and identification of an aggressive signature in low‐grade B‐cell lymphomas
Abstract 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.
Sonicate Fluid Cellularity Predicted by Transcriptomic Deconvolution Differentiates Infectious from Non-Infectious Arthroplasty Failure
Background: Although cellularity is traditionally assessed morphologically, deep sequencing approaches being used for microorganism detection may be able to provide information about cellularity. We hypothesized that cellularity predicted using CIBERSORTx (Stanford University), a transcriptomic-based cellular deconvolution tool, would differentiate between infectious and non-infectious arthroplasty failure. Methods: CIBERSORTx-derived cellularity profiles of 93 sonicate fluid samples, including 53 from subjects who underwent failed arthroplasties due to periprosthetic joint infection (PJI) (abbreviated for the purpose of this study as PJIF) and 40 from subjects who had undergone non-infectious arthroplasty failure (abbreviated NIAF) that had been subjected to bulk RNA sequencing were evaluated. Results: Samples from PJIF and NIAF subjects were differentially clustered by principal component analysis based on the cellularity profile. Twelve of the 22 individual predicted cellular fractions were differentially expressed in the PJIF cases compared with the NIAF cases, including increased predicted neutrophils (mean and standard error, 9.73% ± 1.06% and 0.81% ± 0.60%), activated mast cells (17.12% ± 1.51% and 4.11% ± 0.44%), and eosinophils (1.96% ± 0.37% and 0.42% ± 0.21%), and decreased predicted M0 macrophages (21.33% ± 1.51% and 39.75% ± 2.45%), M2 macrophages (3.56% ± 0.52% and 8.70% ± 1.08%), and regulatory T cells (1.57% ± 0.23% and 3.20% ± 0.34%). The predicted total granulocyte fraction was elevated in the PJIF cases (32.97% ± 2.13% and 11.76% ± 1.61%), and the samples from the NIAF cases had elevated predicted total macrophage and monocyte (34.71% ± 1.71% and 55.34% ± 2.37%) and total B cell fractions (5.89% ± 0.30% and 8.62% ± 0.86%). Receiver operating characteristic curve analysis identified predicted total granulocytes, neutrophils, and activated mast cells as highly able to differentiate between the PJIF cases and the NIAF cases. Within the PJIF cases, the total granulocyte, total macrophage and monocyte, M0 macrophage, and M2 macrophage fractions were differentially expressed in Staphylococcus aureus compared with Staphylococcus epidermidis -associated samples. Within the NIAF cases, the predicted total B cell, naïve B cell, plasma cell, and M2 macrophage fractions were differentially expressed among different causes of failure. Conclusions: CIBERSORTx can predict the cellularity of sonicate fluid using transcriptomic data, allowing for the evaluation of the underlying immune response during the PJIF and NIAF cases, without a need to phenotypically assess cell composition.
Molecular Clusters and Tumor-Immune Drivers of IgM Monoclonal Gammopathies
Abstract Purpose: IgM monoclonal gammopathy of undetermined significance (MGUS) and Waldenström macroglobulinemia (WM) represent a disease spectrum with highly varied therapeutic management, ranging from observation to chemoimmunotherapy. The current classification relies solely on clinical features and does not explain the heterogeneity that exists within each of these conditions. Further investigation is warranted to shed light on the biology that may account for the clinical differences. Experimental Design: We used bone marrow (BM) clonal CD19+ and/or CD138+ sorted cells, matched BM supernatant, and peripheral blood serum from 32 patients (7 MGUS, 25 WM) to perform the first multi-omics approach including whole-exome sequencing, RNA sequencing, proteomics, metabolomics, and mass cytometry. Results: We identified three clusters with distinct pathway activation, immune content, metabolomic, and clinical features. Cluster 1 included only patients with WM and was characterized by transcriptional silencing of genes involved in cell cycle and immune response, enrichment of mitochondrial metabolism, infiltration of senescent T effector memory cells, and aggressive clinical behavior. Genetic/structural alterations of TNFAIP3 were distinct events of this cluster. Cluster 2 comprised both MGUS and WM patients with upregulation of inflammatory response, senescence and glycolysis signatures, increased activated T follicular helper and T regulatory cells, and indolent clinical behavior. Cluster 3 also included both MGUS and WM patients and exhibited intermediate features, including proliferative and inflammatory signaling, as well as glycolysis and mitochondrial metabolism. Conclusions: We have identified three distinct molecular clusters, suggesting a potential biologic classification that may have therapeutic implications.
Phenotype, Function, and Clinical Significance of CD26+ and CD161+Tregs in Splenic Marginal Zone Lymphoma
Abstract Purpose: Regulatory T-cells (Treg) are essential to Tregs homeostasis and modulate the antitumor immune response in patients with lymphoma. However, the biology and prognostic impact of Tregs in splenic marginal zone lymphoma (SMZL) have not been studied. Experimental Design: Biopsy specimens from 24 patients with SMZL and 12 reactive spleens (rSP) from individuals without lymphoma were analyzed by using CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing), CyTOF (mass cytometry) analysis, and flow cytometry to explore the phenotype, transcriptomic profile, and clinical significance of intratumoral Tregs and their subsets. The biological characteristics and cell signaling pathways of intratumoral Treg subsets were confirmed by in vitro functional assays. Results: We found that Tregs are more abundant in SMZL patients' spleens than rSP, and Tregs from patients with SMZL and rSP can be separated into CD161+Treg and CD26+Treg subsets. CD161+Tregs are increased in SMZL but have dysregulated immune function. We found that CD161+Treg and CD26+Tregs have unique gene expression and phenotypic profiles and are differentially correlated with patient outcomes. Specifically, increased CD161+Tregs are significantly associated with a favorable prognosis in patients with SMZL, whereas CD26+Tregs are associated with a poor prognosis. Furthermore, activation of the IL2/STAT5 pathway contributes to the induction of CD26+Tregs and can be reversed by STAT5 inhibition. Conclusions: IL2/STAT5-mediated expansion of CD26+Tregs contributes to a poor clinical outcome in SMZL and may represent a therapeutic opportunity in this disease.
TIGIT Expression Is Associated with T-cell Suppression and Exhaustion and Predicts Clinical Outcome and Anti–PD-1 Response in Follicular Lymphoma
Abstract Purpose: T-cell immunoglobulin and ITIM domain (TIGIT), a member of the immune checkpoint family, is important in normal T-cell biology. However, the phenotypical profile and clinical relevance of TIGIT in follicular lymphoma is largely unknown. Experimental Design: Biopsy specimens from a cohort of 82 patients with follicular lymphoma were analyzed using mass cytometry to explore the phenotype and biological and clinical significance of TIGIT+ T cells. Results: TIGIT is highly expressed on intratumoral T cells and its expression alters T-cell phenotype in follicular lymphoma. TIGIT is abundantly expressed on Treg cells, resulting in an enhanced suppressive property. TIGIT expression on non-Treg/TFH T cells defines a population that exhibits an exhausted phenotype. Clinically, increased numbers of TIGIT+ T cells are associated with inferior patient outcomes and poor survival. We observe that anti–PD-1 therapy with pembrolizumab alters the phenotype of TIGIT+ T subsets and identifies a role for CD28 expression on TIGIT+ T cells in treatment response. Conclusions: The current study provides a comprehensive analysis of the phenotypic profile of intratumoral TIGIT+ T subsets and their prognostic relevance in follicular lymphoma. Inhibition of TIGIT signaling may be an additional mechanism to prevent T-cell suppression and exhaustion in B-cell lymphoma.
Targeting glycogen synthase kinase 3 for therapeutic benefit in lymphoma
Abstract 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.