Comprehensive Integration of Whole Genome, Transcriptome and Methylation Profiling Reveals Novel Gene Dysregulation Including IL15, SOCS6 and CARD11 Associated with MYD88 and CXCR4 Genotype Status in WM

Abstract Background: Waldenstrom's Macroglobulinemia/IgM lyphoplasmacytic lymphoma is a B-cell lymphoma defined by highly recurring mutations in MYD88 (95-97%) and CXCR4 (30-40%) by whole genome, PCR and Sanger sequencing (Treon et al NEJM 2012, Hunter et al Blood 2013, Xu et al, BJH 2016). Using next generation RNA sequencing studies (RNASeq) of these same group of patients, we subsequently observed that MYD88 and CXCR4 mutation status were the primary determinants of differential gene expression in WM (Hunter et al, Blood 2016). We have now integrated enhanced reduced representation bisulfite sequencing (ERRBS) data into the existing WGS and RNASeq data to present of more complete picture of genomic regulation in WM. Methods: CD19+ selected bone marrow samples from 54 patients with WM, and CD19- peripheral blood mononuclear samples were used for tumor and germline controls, respectively. WGS and RNASeq were conducted as previously described (Hunter et al, Blood 2014; 2016) and methylation profiling was conducted using ERRBS at the Weil Cornell Medical Epigenomics Core and analyzed using Bismark. Results were filtered for methylation sites supported by at least 10 reads across all of the samples. Differential methylation analysis was conducted at the individual site level and aggregated promoter level using an established edgeR protocol. Differential methylation data was then compared with Salmon derived gene and isoform expression data which was available for 49/54 (91%) of the samples. Results: Clinical characteristics of patients were as follows: Median age 60 (range 40-87 years); BM involvement 43% (range 4-95%); serum IgM 3590 (range 416-8320 ug/dl). Most patients (67%) were untreated. WGS data revealed 33 (61%), 18 (33%), and 3 (6%) were MYD88Mutant/CXCR4Wild-Type, MYD88Mutant/CXCR4Mutant and MYD88Wild-Type/CXCR4Wild-Type, respectively. As was observed in our RNASeq analysis, pairwise multidimensional scaling of the methylation status of the top 2,000 high variance promoters revealed segregation of the MYD88Mutant/CXCR4Wild-Type from both MYD88Mutant/CXCR4Mutant and MYD88Wild-Type/CXCR4Wild-Type (Figure 1A). Log fold change in methylation in differentially methylated promoters from MYD88Mutant/CXCR4Mutant and MYD88Wild-Type/CXCR4Wild-Type samples relative to MYD88Mutant/CXCR4Wild-Type were negatively correlated with log-fold change in gene expression (Rho = -0.463, p Conclusions: The studies provide the first comprehensive insights into epigenomic regulation of WM and show that MYD88 and CXCR4 mutation status drives methylation and confers a distinct transcriptomic profile that includes genes such as IL15, SOCS6, CARD11, HIF1A, and PIK3R5 with important pro-survival and immune regulatory roles. Download : Download high-res image (327KB) Download : Download full-size image Disclosures Hunter: Pharmacyclics: Consultancy. Castillo: Abbvie: Consultancy, Research Funding; Beigene: Consultancy, Research Funding; Genentech: Consultancy; Janssen: Consultancy, Research Funding; Millennium: Research Funding; Pharmacyclics: Consultancy, Research Funding. Treon: Pharmacyclics: Consultancy, Other: Travel, Accommodations, Expenses, Research Funding; BMS: Research Funding; Johnson & Johnson: Consultancy; Janssen: Consultancy, Other: Travel, Accommodations, Expenses.