Abstract OBJECTIVE Our goal is to identify new pathways perturbed in and which may lead to GRS. METHODS We collected brain transcriptome levels from healthy controls (n=5), patients with GRS (n=9), glioma without seizures (non-GRS, n=8), and idiopathic temporal lobe epilepsy (iTLE, n=7) using Illumina-based poly-A capture high-throughput RNA sequencing approach. Differential gene expression analysis was conducted to ascertain expression differences between the four groups. We identified molecular pathways enriched for differentially expressed genes (DEG) using Ingenuity Pathway Analysis. DEGs were defined as those with log2 fold change ≥ 2 and q-value < 0.01 when compared to GRS. RESULTS Principal component analysis revealed that diagnostic grouping was the main determinant of the transcriptome variance. We identified 212 DEGs in the GRS vs. non-GRS analysis. Among these, 108 genes were up-regulated and 104 genes were down-regulated. Among the genes that were down‐regulated, there was significant enrichment for those involved in acute phase response signaling and immune cell trafficking (e.g. H19, BMP5, CXCL8). Compared to iTLE samples, there were 1015 DEGs in the GRS samples. Among these 877 genes were up-regulated and 138 genes were down-regulated. There was significant enrichment for genes involved in cell-to-cell signaling among those that were either uniquely up‐regulated (e.g. NRK, RGS10, SLC4A4, and SELPLG) or down‐regulated in the GRS samples (e.g. SLC9A3 –Na+/H+ exchanger; SLC5A5-sodium-glucose co-transporter family). Further, two of the most down-regulated genes in samples with GRS, SCLC17A6- vesicular glutamate transporter 2 and FOSB, have previously been shown to be involved in glutamatergic transmission. CONCLUSIONS Our findings suggest that inflammation-related processes and alterations in cell-to cell signaling may underlie the pathogenesis of GRS. Our results also provide support for perturbations in glutamatergic transmission in GRS. Our findings require replication in larger human cohorts and functional experiments to explore the function of the genes perturbed in GRS.
O5‐02‐04: INTEGRATIVE GENOMIC PROFILING TO IDENTIFY GENES AND PATHWAYS ASSOCIATED WITH VASCULAR RISK IN AD
Using Whole Exome Sequencing to Identify Genetic Variation and Polymorphisms Associated with Graft Versus Host Disease in Allogeneic Stem Cell Transplant Recipients
Abstract Objective: Allogeneic stem cell transplant (ASCT) is a potentially curative treatment for many hematologic malignancies. HLA match is the principal factor in donor selection but the total antigenic landscape as seen by the donor immune cells extends beyond the HLA and minor HLA antigens. Thus, despite "perfect" HLA match graft versus host disease (GVHD) caused by alloreactive donor immune cells reacting to disparate antigens in the recipient, remains a significant cause of morbidity and mortality. Whole exome sequencing (WES) which assays all the protein coding nucleotide sequences provides a powerful means to assess the entire antigenic makeup of the individual. We examined the WES of sibling donor-recipient pairs undergoing BMT to assess the influence of the extent of genetic disparity and specific gene polymorphisms on GVHD. Methods: Sixteen matched sibling donor-recipient pairs receiving allogeneic stem cell transplants between June 2008 and April 2013 had samples available for retrospective WES. Transplant data including conditioning regimen, age of donor and recipient at the time of transplant, development and type of GVHD and cause of death were recorded. The sequencing was performed using Illumina HiSeq2000 and the single nucleotide variants (SNVs) and small insertions and deletions (INDELs) were identified using the haplotype caller within the Genome Analysis Toolkit GATK v2.7. Functional variants including non-synonymous, frame-shift, truncating, and splice site SNVs and INDELs identified were compared between donor and recipient and among those with GVHD and those without GVHD (non-GVHD). One-sided Fisher's exact test was used to determine the statistical significance of the SNPs in GVHD and non-GVHD. Candidate genes with a p value of <0.05 were further evaluated for their biochemical pathway and possible role in inflammatory or immune responses. Results: Nine of 16 patients had a diagnosis of AML or high grade myelodysplastic syndrome, 3 had a diagnosis of non- Hodgkin lymphoma, 1 had Hodgkin lymphoma, 1 had multiple myeloma, and 2 had CML. Six patients were conditioned with a myeloablative regimen, others with a reduced intensity conditioning regimen. Eight of 16 patients were female and 7 of 16 donors were female. Nine of 16 patients developed stage II-IV acute and/or chronic GVHD with 7 of them developing GVHD of the skin. A total of 12,623 genes with 43,670 functional variants were identified using WES. 134 genes had a p value < 0.05 and enrichment analysis was performed on these genes. 12 out of 134 genes were associated with gene ontology (GO) term "immune system development" with an unadjusted enrichment p-value 0.0015. Three (IL31RA, OSMR, SIRPA) out of 134 genes were associated with "inflammation_Jak-STAT pathway" with an adjusted enrichment p-value 0.07. Variants of IL31RA, OSMR, and PTX3 were identified as genes with polymorphisms associated with GVHD (p= 0.003, 0.01, and 0.01, respectively) with 5 variants identified in IL31RA, 7 variants in OSMR, and 2 variants in PTX3. Further investigation of these three genes revealed their role in inflammatory or immune responses. Discussion: GVHD is a cytokine mediated immune-inflammatory process that can occur despite "perfect" HLA matching. WES identified the extensive genetic variability between donors and recipients implying large allo-reactivity potential that can lead to GVHD. Quantification of genetic variation may allow for more precise determination of GVHD risk and better selection of donors. In our cohort, WES identified polymorphisms in IL31RA, OSMR, and PTX3 as statistically significantly associated with development of GVHD. The protein products of these three genes have each been shown to be involved in inflammatory processes and are influenced by other cytokines. IL31RA and OSMR are found in many tissues including the skin and lungs. They function both independently and as a heterodimer receptor capable of binding multiple cytokines including IL-31, IL-6, and TNF-α. PTX3 is a reactive protein that increases in response to cytokine release related to systemic inflammation. Levels of this protein have been used as a prognostic marker in many disease processes with high levels correlating with poor outcomes. The small sample size of our cohort precludes any firm conclusions but provides promising hypotheses for further evaluation in larger studies. Disclosures No relevant conflicts of interest to declare.
Therapeutic Sensitivity of CD20- Waldenströms Macroglobulinemia Cells Is Determined By Underlying Genomic and Epigenetic Events
Abstract Waldenströms Macroglobulinemia (WM) is a unique form of lymphoplasmacytic lymphoma, in which the tumor compartment is comprised of B-lymphocytes, plasmacytoid lymphocytes and plasma cells. The majority of WM cases are immunophenotypically characterized by expression of surface antigens that are present on both B-cells (CD19/20) and plasma cells (CD38/138). However, in advanced stage disease, which has been treated with various chemoimmunotherapeutics the molecular profile of WM cells can drastically shift resulting in a majority of the WM tumor compartment exhibiting loss of CD20 and concurrent upregulation of atypical surface antigens (Barakat et al, Am J Clin Pathol 2011; Morice et al Mod Pathol 2009). These external immunophenotypic changes are associated with underlying epigenetic and genomic alterations, whose integrated study may provide further insight into the biology of advanced stage WM. Moreover, such an analysis holds clinical potential by uncovering of therapeutic vulnerabilities, which can direct the selection of appropriate therapeutics for patients with CD20- WM. With this in mind, we conducted comprehensive molecular testing coupled with functional studies to understand if underlying genomic/epigenetic changes correlated with drug sensitivity using appropriate preclinical models of CD20+ WM (BCWM.1 cell line) and CD20- WM (RPCI-WM1 cell line). Importantly, the index patient from whom the CD20- cells were developed had received a total of 6 lines of treatment, two of which contained rituximab and was documented as being refractory to rituximab. Immunophenotyping of CD20- cells showed differential downregulation of 19 and upregulation of 7 CD antigens relative to CD20+ WM cells (cutoff at 1.5 fold; >20% gated expression). Upregulated antigens were associated with a phenotype reminiscent of a memory B-cell with plasmacytoid features. Next, we conducted RNA-Seq analysis (Illumina HiSeq), which revealed downregulation of 6,673 and upregulation of 4,594 genes (>1.5 fold cutoff) in CD20- WM cells compared with CD20+ cells. Gene pathway enrichment analysis (NextBio, Illumina Inc.) of mRNA differentially expressed in CD20- cells showed significant enrichment of upregulated genes corresponding to Protein Transport (p<1.8E-98) and Response to Unfolded Protein Gene Ontology groups (p<5.9E-20). Interestingly, pathway analysis of downregulated genes in CD20- WM cells showed enrichment within the B-cell receptor-signaling (BCR) pathway (Broad MSigDB) (p<6.1E-14) relative to CD20+ WM cells. BCR signaling has been linked to dysregulated expression of several micro-RNAs (miRs). Consistent with pathway analysis results, we observed decreased levels of numerous miRs (Exiqon qPCR) in CD20- vs. CD20+ WM cells, which target and regulate components involved in BCR activation. To examine further epigenetic changes that may account for decreased BCR signaling, we examined the global methylation profile (Infinium HumanMethylation 450K array) of CD20- and CD20+ WM cells and noted hypermethylation of the promoter sites (200 base pairs upstream) for a number of transcription factors involved in the BCR pathway. We therefore examined sensitivity of CD20- WM cells to the BTK inhibitor ibrutinib, which abrogates BCR signaling. As anticipated, CD20- WM cells retained >90% viability and demonstrated significantly less apoptosis in presence of ibrutinib (4uM) compared to CD20+ cells, suggesting decreased reliance on canonical BCR signaling for survival. In parallel and guided by pathway analysis results, we exposed both CD20- and CD20+ WM cells to a novel 19S proteasome inhibitor (VLX1570), which caused profound proteotoxicity and ER stress at concentrations ranging between 250 – 500nM. We noted CD20- cells were significantly more sensitive (p<0.01) to 19S proteasome inhibition than CD20+ WM cells indicating the formers critical dependence on tightly regulated protein homeostasis. In summary we describe for the first time, a comprehensive molecular and immunophenotypic investigation of CD20- WM cells compared to CD20+ WM cells, using validated preclinical models of the two disease phenotypes. Overall our data suggests significant architectural differences between CD20+ and CD20- WM, which can potentially impact the behavior of the cells and their sensitivity towards novel therapeutics being developed for treatment of advanced stage WM. Disclosures Martin: Janssen: Honoraria.
O3‐06‐02: POLYGENIC SCORE ANALYSIS OF EXONIC VARIANTS IN AN IMMUNE CO‐REGULATORY NETWORK IDENTIFIES NOVEL PROTEIN‐ALTERING VARIANTS THAT ASSOCIATE WITH ALZHEIMER'S DISEASE