Citations

Abstract Introduction: Stromal tumor-infiltrating lymphocytes (sTILs) have established prognostic and predictive significance in triple-negative breast cancer (TNBC). However, the roles of other immune cells in TNBC are less well-established. We performed high-plex quantitative spatial profiling in a cohort of early-stage TNBC to 1) apply spatial context to tumoral immune landscapes and 2) identify immune proteins associated with clinical outcomes, independently of TILs and other established prognostic clinicopathologic variables, in patients (pts) treated with or without adjuvant chemotherapy (CTX). Methods: The Mayo TNBC cohort comprises pts with centrally-verified, CTX-naive tumors resected from 1985-2012. Using a cohort-based TMA, with Nanostring GeoMX DSP, we quantitated 58 proteins within spatially-distinct intra-epithelial, cytokeratin-positive tumor segments and adjacent cytokeratin-negative/nuclei-positive stromal segments. Differentially-expressed (DE) proteins were identified using a negative binomial generalized linear model (SNR>2, p< 0.05) and a target DE protein set was dichotomized (80th percentile). After adjusting for prognostic clinicopathologic variables, proteins associated with recurrence-free survival (RFS, defined as time from surgery to either local, regional, and distant recurrence, or death by any cause) were identified by performing variable selection using the Akaike Information Criterion (AIC) obtained from fitting all possible Cox proportional hazards regression models (performed separately for intra-epithelial/stromal segments, and in groups +/- adjuvant CTX. Results: From the TNBC TMA, DSP data (N=250 tumors) included 169 pts who received adjuvant CTX+ and 81 who did not (CTX-). Overall, 85/250 developed recurrent disease. In the CTX+ group, intra-epithelial tumor segments from pts without recurrent disease were enriched in 10 immune proteins, including CD8, markers involved in antigen presentation/dendritic cells (CD11c, CD40, HLA-DR) or NK cells (CD56) (FC: 1.4-2.1, p<0.05); CD14 was increased in stroma (FC: 1.5, p<0.05). In contrast, in the CTX- group, both the intra-epithelial tumor and stromal segments from pts without recurrences were enriched in immune proteins (N= 12 and 15 respectively; FC 1.6-5.5, p< 0.05) most markedly CD40, IDO1 and HLA-DR (FC: 3.2-5.5, p< 0.05). Overall, CD3, CD4, CD27, CD44, and ICOS among others were enriched only in the CTX- group; CD14 and CD56 were enriched only in the CTX+ group. Based on these spatial data, biologic function and DSP data from another set of TNBC (FinXX trial), CD11c, CD14, CD27, CD40, CD56, and IDO1 were selected for RFS analysis. After applying our model selection criterion and adjusting for pt age at surgery, tumor size, lymph node status, and sTILs, intra-epithelial CD56 was independently associated with improved RFS in the CTX+ group (HR: 0.31[0.12, 0.81]). In the CTX- group, intra-epithelial CD11c was independently associated with improved RFS (0.10 [0.01, 0.81]). Conclusion: In this early-stage TNBC cohort, spatially-distinct tumor immune landscapes were associated with RFS but differed according to receipt of CTX after surgical resection. In the patients who received CTX, the intra-epithelial compartment, rather than stromal compartment, was immune-enriched in pts without recurrences. Among a targeted protein set, intra-epithelial CD56 remained associated with improved outcomes, independent of sTILs and other clinicopathologic features. In the CTX- group, spatial landscapes were more balanced, and intra-epithelial CD11c was independently associated with improved outcomes. These data provide insight into the spatial context of intrinsic immune landscapes in TNBC, and identify candidate prognostic immune biomarkers which may inform therapeutic strategies. Citation Format: Jodi M Carter, Saranya Chumsri, David W Hillman, David M Zahrieh, Yaohua Ma, Xue Wang, Jennifer M Kachergus, Judy C Boughey, Minetta C Liu, Krishna R Kalari, JC Villasboas, Roberto A Leon Ferre, Fergus J Couch, Matthew P Goetz, E. Aubrey Thompson. Intra-epithelial tumor immune landscapes are associated with clinical outcomes in early-stage triple-negative breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS16-01.

Abstract Introduction: The immune microenvironment of high-risk HR+ breast cancer (BC) is poorly understood. BEAUTY is a prospective neoadjuvant chemotherapy (NAC) study of stage I-III BC patients treated with neoadjuvant weekly taxane followed by anthracycline-based chemotherapy. Among clinical luminal B BC (St. Gallen criteria) from BEAUTY, we used high-plex digital spatial profiling (DSP) to characterize the intra-epithelial tumor and stromal immune landscapes to 1) assess the impact of NAC on these landscapes and 2) identify immune biomarkers predictive of response to NAC. Methods: The tissue set included FFPE sections of resected tumors from 35 patients (median 51y; range: 21-71y) with clinically-defined luminal B BC (ER > 10%/grade 2/Ki67 ≥ 15% or ER > 10%/grade 3), and 16 paired pre-NAC biopsies. Nanostring GeoMX DSP was used to quantitate 58 immune and BC biomarker proteins in intra-epithelial, cytokeratin-positive tumor segments and adjacent stromal (cytokeratin-negative/SYTO13 (nuclear stain)-positive) segments. Data were normalized using the geometric mean of two negative controls (RbIgG and MsIgG1). A general linear model with negative binomial identified differentially-expressed (DE) proteins in pre/post-NAC tumors. Based on DE protein data and biologic function, a targeted protein subset (N=19) was evaluated in pre-NAC biopsies for associations between protein abundance and residual cancer burden (RCB) class (0-II vs. III) (Wilcoxon-rank sum test: p-value < 0.025 considered significant) or RCB ‘breast only’ index (Spearman correlation). Results: Comparing tumor segments in pre-NAC specimens, intra-epithelial segments were predictably enriched in cytokeratin, ER, PR and Ki-67, but also CD127 and NY-ESO-1; whereas stromal segments were enriched in proteins associated with T cell subsets (e.g. CD3, CD4, CD8), macrophages (CD68 and CD163), antigen presentation/dendritic cells (CD11c, HLA-DR) and immune checkpoint proteins (PD-L1, PD-L2, B7-H3, TIM-3, VISTA) among others. While NAC did not alter the spatial distribution of proteins (intra-epithelial vs. stromal segments), it markedly attenuated the immune landscape, with decreased abundance of functionally-diverse immune proteins (e.g. CD45, CD3, CD4, CD127, granzyme B, CTLA4; STING, B7-H3, CD11c, and CD68, log2FC: 0.27-1.52 p < 0.05), including low abundance proteins PD-1, PD-L1, PD-L2, CD20, and OX40L (log 2FC: 0.15-1.16, p < 0.05). Both PR and Ki-67 decreased in post-NAC tumors whereas ER was not significantly altered (p < 0.05). CD8, CD14, CD163, HLA-DR, IDO-1 and TIM-3 were not significantly altered by NAC. In the pre-NAC biopsies (which had subsequent residual tumor burdens of RCB class 0 or I (n=1 each), II (n=6), and III (n=8), and median RCB ‘breast only’ index of 3.39 (range: 0.00-38.02), 19 proteins were used for RCB analysis (CD3, CD4, CD8, CD14, CD34, CD44, CD45, CD68, CD127, CD163, CTLA4, granzyme B, STING, B7-H3, fibronectin, Ki-67, HLA-DR, SMA and TIM3). Among them, stromal CD127 was significantly higher in RCB class III than RCB class 0-II (p=0.021) and RCB ‘breast only’ index was positively correlated with intra-epithelial granzyme B (rho = 0.61; p=0.012), and negatively correlated with intra-epithelial Ki-67 (rho= -0.71; p=0.0022). Conclusion: In this series of clinical luminal B BC, NAC markedly attenuated the tumoral immune landscapes with a small set of “NAC-resistant” immune proteins. Among a targeted set, stromal CD127 was significantly higher in RCB III, and RCB breast-only index was positively correlated with intra-epithelial granzyme B. These data provide insight into the impact of NAC on HR+ BC, and identify potential immune biomarkers to predict response to neoadjuvant chemotherapy. Citation Format: Jodi M Carter, Judy C Boughey, Jun He, Vera J Suman, Xue Wang, Jennifer M Kachergus, Krishna R Kalari, Liewei Wang, Richard Weinshilboum, Ann M Moyer, Sarah A McLaughlin, Alvaro Moreno-Aspitia, Donald W Northfelt, Richard J Gray, James N Ingle, E. Aubrey Thompson, Matthew P Goetz. Neoadjuvant chemotherapy selectively alters spatially-defined immune landscapes in clinical luminal B HR+/HER2- breast cancers: Analysis of the breast cancer genome guided therapy study (BEAUTY) [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD7-05.

AbstractIntron retention (IR) has been implicated in the pathogenesis of complex diseases such as cancers; its association with Alzheimer's disease (AD) remains unexplored. We performed genome‐wide analysis of IR through integrating genetic, transcriptomic, and proteomic data of AD subjects and mouse models from the Accelerating Medicines Partnership‐Alzheimer's Disease project. We identified 4535 and 4086 IR events in 2173 human and 1736 mouse genes, respectively. Quantitation of IR enabled the identification of differentially expressed genes that conventional exon‐level approaches did not reveal. There were significant correlations of intron expression within innate immune genes, like HMBOX1, with AD in humans. Peptides with a high probability of translation from intron‐retained mRNAs were identified using mass spectrometry. Further, we established AD‐specific intron expression Quantitative Trait Loci, and identified splicing‐related genes that may regulate IR. Our analysis provides a novel resource for the search for new AD biomarkers and pathological mechanisms.

AbstractBackgroundUsing an objective mathematical algorithm to assess corticolimbic involvement of neurofibrillary tangle pathology, we identified three AD subtypes: hippocampal sparing (HpSp) AD, typical AD, and limbic predominant AD. Typical AD brains were representative of the expected patterns of hippocampal and cortical involvement as outlined by Braak tangle stage. In contrast, we discovered an extreme phenotype that exists outside of the Braakian‐concept of neurofibrillary tangle patterns. Hippocampal sparing AD cases demonstrate unexpected sparing of the hippocampus relative to severely involved association cortices and limbic predominant AD cases demonstrate inundation of the hippocampus relative to mildly involved association cortices. Using these AD subtypes, our objective was to leverage our understanding of disease spectrum to uncover transcriptomic changes that underlie selective vulnerability of the hippocampus in AD.MethodWe performed RNA‐sequencing in 40 neuropathologically diagnosed AD cases (10 HpSp, 20 Typical and 10 Limbic) and 15 controls to uncover gene expression changes associated with phenotype differences allowing us to prioritize 44 genes for further exploration. Validation with NanoString and quantification with digital pathology were performed in an expanded cohort of 158 AD cases (36 HpSp, 79 Typical and 35 Limbic) and 32 controls. We quantified digital pathology measures of early tangle pathology (CP13), mature tangle pathology (Ab39), amyloid‐β (33.1.1). In addition, we measured cellular markers for microglia (CD68), endothelia (CD34), and astroglia (GFAP). Deep learning based prediction models were employed, which nominated the SLC38A2 gene as a predictor of AD.ResultThe SLC38A2 gene was originally prioritized in our RNA‐sequencing data based upon differential expression between controls and typical AD. We observed upregulation of SLC38A2 in all AD subtypes. RNA‐seq gene expression measures validated well with NanoString for SLC38A2 (R=0.97, p<0.001). Regression analysis identified early tangle (p=0.011) and late tangle measures (p=0.026) as significant predictors of gene expression, but not amyloid measures (p=0.74). Microgliosis approached significance (p=0.052), but with no contribution from endothelial cell burden (p=0.17) or astrogliosis (p=0.58). Age approached significance (p=0.06), but not sex (p=0.88) or APOE‐ε4 status (p=0.10).ConclusionOur data supports consideration of intra‐disease divergence with regard to case stratification and may reveal genes previously masked by heterogeneity of cohorts.

Abstract Background: The intrinsic immune response in triple negative breast cancer (TNBC) has both prognostic and predictive utility, and immune-related biomarkers are now actionable targets in TNBC. Current modalities for immune profiling are limited by lack of spatial context and precise quantitation. We used quantitative digital spatial profiling (DSP) to 1) characterize the immune architecture of tumoral and stromal regions in TNBC and 2) identify prognostic immune signatures. Methods: From the FinXX trial of adjuvant capecitabine therapy, tumor sections (FFPE) matched for patient characteristics, treatment arm and RFS-based outcome (N=44) were assessed with Nanostring GeoMxTM DSP. Briefly, with 3-plex immunofluorescence (pan-cytokeratin, CD45, CD68 with SYTO-13 dye for nuclei), spatially-defined CD45-enriched or CD68-enriched stromal segments and adjacent tumor segments (N=950) were selected for digital quantitation (NanoString nCounter) of 39 immune proteins (e.g. immune cell profiling proteins: CD3, CD4, CD8, CD11c, CD20, CD56, CD68, HLA-DR, PD-1, PD-L1, Granzyme B), immune drug targets (e.g. VISTA, STING, IDO1), and activation status-related proteins (e.g. ICOS, PD-L2, CD40, CD40L, CD44). Normalized, differentially-expressed (DE) proteins were evaluated, and weighted gene co-expression network analysis (WGCNA) was performed to identify modules of highly co-expressing immune proteins, and association with RFS as a categorical variable was determined. Results: 20 DE proteins in tumor segments were associated with RFS. Among them, those with the highest RFS-associated fold changes (FC) included CD20, CD40, CD56, HLA-DR, Granzyme B, and PD-L2 (FC): 2.1-4.4) (p &lt; 0.05). In contrast, only 4 DE proteins were associated with RFS in stromal segments (Granzyme B, IDO-1, PD-L1, and PD-L2 (FC: 1.5-3.3) (p &lt; 0.05). By WGCNA, 2 tumor-based immune modules [IM1 (lymphocyte infiltration) and IM2 (lymphocyte activation)] were identified based on their co-expression. The top three co-expressed proteins in IM1 were CD4, CD40, and CD45 (r &gt;0.9). The three top co-expressed proteins in IM2 were PD-L1, PD-L2, and Granzyme B (r &gt; 0.85). Together with CD56 expression, IM1 and IM2 comprised a set of 3 intratumoral immune features that associated with RFS. For recurring tumors, all 3 immune features clustered with low scores . Among non-recurrent tumors, two distinct, non-overlapping clusters with inverse characteristics were identified: high IM2, low IM1, low CD56 vs. low IM2, high IM1, and high CD56. Overall, across all 3 immune features, PDL2, CD40 and CD56 expression most strongly associated with RFS. Conclusion: With spatially-defined, precise quantification of immune proteins, we identified 3 immune features, including two distinct immune modules, IM1 and IM2 that correlate with durable RFS. These immune modules of highly co-expressing proteins were associated with RFS only when present in the tumor compartment but not tumor-adjacent, immune cell-rich stroma. Citation Format: Jodi M. Carter, Saranya Chumsri, Yaohua Ma, Xue Wang, Jennifer M. Kachergus, Douglas Hinerfeld, Heather Ann Brauer, Sarah Warren, Heikki Joensuu, Edith A. Perez, Aubrey Thompson. Quantitative digital spatial profiling reveals novel prognostic intratumoral immune signatures in triple-negative breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1004.

AbstractIntroductionAbnormal gene expression patterns may contribute to the onset and progression of late‐onset Alzheimer's disease (LOAD).MethodsWe performed transcriptome‐wide meta‐analysis (N = 1440) of blood‐based microarray gene expression profiles as well as neuroimaging and cerebrospinal fluid (CSF) endophenotype analysis.ResultsWe identified and replicated five genes (CREB5, CD46, TMBIM6, IRAK3, and RPAIN) as significantly dysregulated in LOAD. The most significantly altered gene, CREB5, was also associated with brain atrophy and increased amyloid beta (Aβ) accumulation, especially in the entorhinal cortex region. cis‐expression quantitative trait loci mapping analysis of CREB5 detected five significant associations (P < 5 × 10−8), where rs56388170 (most significant) was also significantly associated with global cortical Aβ deposition measured by [18F]Florbetapir positron emission tomography and CSF Aβ1‐42.DiscussionRNA from peripheral blood indicated a differential gene expression pattern in LOAD. Genes identified have been implicated in biological processes relevant to Alzheimer's disease. CREB, in particular, plays a key role in nervous system development, cell survival, plasticity, and learning and memory.

AbstractIntroductionMAPT H1 haplotype is implicated as a risk factor for neurodegenerative diseases including Alzheimer's disease (AD).MethodsUsing Alzheimer's Disease Genetics Consortium (ADGC) genome‐wide association study (GWAS) data (n = 18,841), we conducted a MAPT H1/H2 haplotype–stratified association to discover MAPT haplotype–specific AD risk loci.ResultsWe identified 11 loci—5 in H2‐non‐carriers and 6 in H2‐carriers—although none of the MAPT haplotype–specific associations achieved genome‐wide significance. The most significant H2 non‐carrier–specific association was with a NECTIN2 intronic (P = 1.33E‐07) variant, and that for H2 carriers was near NKX6‐1 (P = 1.99E‐06). The GABRG2 locus had the strongest epistasis with MAPT H1/H2 variant rs8070723 (P = 3.91E‐06). Eight of the 12 genes at these loci had transcriptome‐wide significant differential expression in AD versus control temporal cortex (q < 0.05). Six genes were members of the brain transcriptional co‐expression network implicated in “synaptic transmission” (P = 9.85E‐59), which is also enriched for neuronal genes (P = 1.0E‐164), including MAPT.DiscussionThis stratified GWAS identified loci that may confer AD risk in a MAPT haplotype–specific manner. This approach may preferentially enrich for neuronal genes implicated in synaptic transmission.

Abstract Background: Growing evidence supports the critical role of preexisting immune response in triple negative breast cancer (TNBC). However, there are limitations with current evaluation approaches: inability to functionally assess the type of immune infiltration with traditional pathologic evaluation and loss of spatial distribution in conventional high-plex gene or protein expression analyses from the whole tumor section. Here, we report the initial analysis of immune protein expression as a function of spatial distribution and clinical outcomes in TNBC samples. Methods: NanoString GeoMxTM Digital Spatial Profiling (DSP) was used to quantify 39 immune-related proteins in stromal and tumor segments from 44 TNBC samples from the FinXX trial. Samples were matched for patient characteristics, treatment arm (capecitabine vs. 5-fluorouracil), and outcome based on recurrence-free survival (RFS) with 22 samples from patients who recurred and 22 samples from patients with durable RFS. Regions of interest (ROIs) were selected based on expression of cytokeratin (tumor), CD45 (leukocytes), or CD68 (macrophages). Each ROI was segmented into tumor (pancytokeratin-positive area) and stroma (pancytokeratin-negative/nuclear SYTO13-positive area). The general linear model was used for statistical analysis of differential expression with RFS as a categorical variable (recur yes or no). Results: A total of 950 tumor and stroma segments were included in this initial analysis. In both tumor and stroma segments, over-expression of T cell activation markers (CD137, GITR) was associated with better outcome, whereas T cell markers (CD3, CD4, CD8) were not significantly associated with outcome. In tumor segments alone, improved outcome was significantly associated with increased protein expression [&gt; 2-fold change (FC) at p&lt;0.001] of CD56, PD-L2, HLA-DR, CD137, GITR, and CD40. In CD45-enriched stroma, improved outcome was associated with elevated expression (FC&gt;2.0, p&lt;0.001) of PD-L2; whereas durable RFS was associated with elevated PD-L2 and IDO1 expression in CD68-enriched stroma. In contrast, macrophage/dendritic cell markers CD68 and CD11c were not associated with outcome. In tumor cells adjacent to CD45-enriched stroma, durable RFS was associated with increased abundance (FC&gt;2.0, p&lt;0.001) of PD-L2, CD56, CD27, GITR, CD20, HLA-DR, and IDO1. Similarly, in tumor segment associated with CD68-enriched stroma, proteins associated with RFS included PD-L2, HLA-DR, CD56, GITR, and CD137. Among 39 immune function proteins, only elevated CTLA4 expression in CD68-enriched stroma was associated with recurrence (FC=0.55, p&lt;0.001), whereas tumor-segment CTLA4 was associated with RFS. Conclusions: Using an in-depth analysis to precisely quantify the abundance of multiple immune function proteins in a spatially defined manner, we observed that PD-L2, IDO1, and T cell activation markers were robustly associated with durable RFS in both tumor and stromal segments. In contrast, MHC components (HLA-DR, beta-2-microglobulin), B cell markers (CD20), and NK cell markers (CD56) were strongly associated with favorable outcome in tumor but not in stromal segments. Our study highlights the conclusion that the immune landscape of TNBC is far too complex to be encompassed by any single molecular marker, and more detailed analyses of the DSP data reported here are ongoing with a view towards using quantitative multiplex analyses to refine our understanding of how therapeutic outcome is influenced by interactions among immune cells and between immune and tumor cells. Acknowledgements: Supported by the Breast Cancer Research Foundation (BCRF18-161), Bankhead Coley (6BC05 Florida Department of Health), 26.2 with Donna Foundation, US National Cancer Institute (CA15083), and the Canadian Cancer Center Citation Format: Saranya Chumsri, Douglas Hinerfeld, Jennifer M. Kachergus, Yaohua Ma, Heather A Brauer, Sarah Warren, Xue Wang, Torsten O. Nielsen, Karama Asleh, Heikki Joensuu, Edith A. Perez, E. A. Thompson. Digital spatial mapping of the immune landscape of triple negative breast cancer reveals novel features of immune-tumor cell interaction [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr PD5-06.

Abstract Background: Waldenstršms macroglobulinemia (WM) is an indolent B-lymphoplasmacytic lymphoma, which is highly sensitive to 20S proteasome inhibitors such as bortezomib. However, clinical resistance to the drug can be observed early at onset of therapy (primary or inherent resistance) or develop over time, through repeated exposure to the drug (acquired resistance). Although prior investigations into the biological mechanisms associated with bortezomib resistance (BR) have identified aberrations at the genomic (PSMB5 mutation) and transcriptional level (loss of XPB-1), there is limited data on the epigenetic abnormalities that support resistance to 20S proteasome inhibition. With this in mind, we conducted a comprehensive methylation and mRNA expression analysis in WM cells focusing on the molecular patterns that characterize inherent BR (iBR) as well as acquired BR (aBR). Materials and Methods: Bortezomib sensitive (BS) human WM cell lines, BCWM.1 and MWCL-1, the iBR WM cell line RPCI-WM1 (developed from a bortezomib-refractory WM patient) as well as aBR subclones (BCWM.1/BR, MWCL-1/BR) were used in the study. Methylation status of WM cells was profiled by Illumina Infinium 450k methylation array. Methylation IDAT files were background subtracted using GenomeStudio v1.9. Beta mixture quantile dilation (BMIQ) method in R watermelon package was used to adjust type II probes into distribution characteristic of type I probes. For each gene, probes were partitioned into 6 regions using R IMA package and the ones belonging to the same region were averaged, including transcription start sites (TSS; both 200 bp and 1500 bp upstream), 5' untranslated region (UTR), and 3'UTR, the first exon, and genebody. Gene regions showing a significant adjusted p-value and a beta difference of >0.1 or <-0.1 were selected for differential hypermethylation and hypomethylation, respectively. Correlation with gene expression was conducted only on genes showing >1.5 fold change in expression. Overview of methodology is shown in Figure 1. Results: We first assessed the differential methylation pattern (TSS, 200bp upstream of gene body) of bortezomib-sensitive (BS, BCWM.1, IC50 6nM) and iBR cells, (RPCI-WM1). A total of 3,801 hypermethylated and 1,865 hypomethylated DNA regions were observed in iBR cells relative to BS cells (Figure 2A). In iBR cells, a negative correlation between hypermethylated DNA and corresponding gene expression was found in 1,039 genes (Spearman r = -0.12) (Figure 2B). In parallel, we compared the methylation pattern (TSS 200bp) of BS WM cells (BCWM.1 and MWCL-1) with their aBR derivatives, which display ~20 – 50 fold resistance to bortezomib. aBR subclones showed differential hypermethylation of 1,521 DNA regions and hypomethylation of 651 genomic regions relative to BS WM cells. These findings were not unexpected as aBR cell lines, being derivatives of BS WM cells, are clonally related to the former, whereas iBR cells are derived from a different clone. Meta-analysis of iBR and aBR followed by pathway enrichment (NextBio, Illumina Inc), revealed a significant clustering of commonly hypermethylated genes within the Validated Transcriptional Targets of TP63 isoforms Canonical Pathway (p<1.9-E5, Figure 2C) and the CEBPA binding site Geneset 2 (p<8.2-E9, Figure 2D); pathways whose genes are actively involved in suppression of cell growth and differentiation. Further validation and functional relevance of candidate gene targets from this analysis, which are relevant to WM pathobiology, are currently being investigated. Conclusion: This is the first study to examine and correlate the global methylation and gene expression profiles in WM tumor cells that are inherently wired to resist 20S proteasome inhibition as well as tumor cells, which have acquired resistance through continuous exposure to bortezomib. Although considerable differences were observed in the methylation profiles of WM cells displaying inherent vs. acquired resistance, a number of common genes emerged and pointed towards transcriptional repressors, whose inactivation through methylation has been linked to other cancers. Altogether, these findings uncover cellular systems in WM cells, which are affected by DNA methylation, subsequent mRNA suppression and which may be associated as a whole with both inherent and acquired forms of bortezomib resistance. Disclosures No relevant conflicts of interest to declare.