Next Generation Sequencing Identifies a Distinct Transcriptional Profile, Including Isoform Dysregulation That Segue with Genomic Alterations in Waldenstrom's Macoglobulinemia
Authors: Hunter Zachary , Xu Lian , Yang Guang , Tsakmaklis Nicholas , Liu Xia , Chen Jie , Manning Robert R , Chen Jiaji , Patterson Christopher J , Castillo Jorge J , Anderson Kenneth C , Munshi Nikhil C. , Treon Steven P 
Published: October 2019

Abstract
Background: Whole genome sequencing has identified highly prevalent somatic mutations including MYD88, CXCR4 and ARID1A, as well as gene losses in Waldenström's Macroglobulinemia (WM) (NEJM 367(9):826-33; Blood 123(11):1637-46). At least three genomic defined subpopulations of WM have been identified based on MYD88 and CXCR4 mutation status (Blood 123(18):2791-6).The regulatory impact for these genomic alterations remain to be clarified.
Methods: Next generation transcriptome profiling was performed using Illumina HiSeq. RNA was isolated from CD19-selected bone marrow cells from 57 WM patients, as well as memory (CD19+CD27+) and non-memory (CD19+CD27-) B-cells from 9 healthy donors. Differential gene expression analysis was performed based on MYD88, CXCR4 and ARID1A mutation status, as well as common cytogenetic abnormalities encountered in WM including amplifications in chromosomes 3q, 4, 6p, and deletions in 6q. Reads were aligned to KnownGene HG19/GRCh37 reference using STAR. Read counts per gene were obtained using featureCounts from Rsubread, and analyzed using voom from the edgeR/limma Bioconductor packages in R. Differential isoform expression was assessed using the Cufflinks software suite. Functional enrichment analysis was conducted using Ingenuity Pathway Analysis.
Results: Transcriptional profilingof WM cells showed a stronger correlation with healthy donor memory B-cells, their putative cell of origin. Differential gene expression analysis between WM patients and healthy donors generated a distinct transcriptional profile for WM that included dysregulation of PIM1, RGS7, IL11RA and EPHB4, as well as differential isoform usage in TP53, PRDM1 and XBP1. MYD88mutatedCXCR4wild-type(WT) WM patients consistently over-expressed genes unique to this cohort including, IL17RB, GPER, WNT5A, WNK2, CABELS1, and PRDM5. MYD88mutated patients who harbored nonsense (NS) or frame-shift (FS) mutations in CXCR4 showed functional enrichment of genes that indicated inhibition of Toll-like receptor inflammatory pathways. CXCR4 mutated patients showed strong up-regulation of CXCR7 and TSPAN33, and exhibited isoform level dysregulation of MEF2B, FOXO3, KDM2A and PRKAG2. Unsupervised clustering differentiated MYD88mutatedCXCR4WT from MYD88mutatedCXCR4NS/FS patients, while samples from MYD88WTCXCR4WT clustered with CXCR4NS/FS group. Clusters based on clonality and disease burden were also observed (Figure 1). CXCR4NS/FS transcripts were preferentially expressed over CXCR4WT transcripts despite the subclonal presence of CXCR4mutations. These findings were further validated by comparative DNA versus cDNA Sanger sequencing indicating allelic dysregulation within CXCR4NS/FS subclones. Controlling for MYD88 and CXCR4 mutation status, the presence of ARID1A mutations and cytogenetic abnormalities generated distinct transcriptional profiles. Random forest regression analysis identified subsets of genes strongly associated with bone marrow disease involvement, serum IgM and hemoglobin levels. Notably, increased bone marrow disease burden was associated with increased CXCL13, and decreased TP53 and RBL1 expression. Likewise, higher levels of serum IgM corresponded with increased IL27RA expression.
Conclusion: Using next generation sequencing, we have identified a distinct transcriptional profile, including isoform dysregulation that segue with highly prevalent genomic mutations in WM. The findings provide valuable insights into the molecular pathogenesis and clinical presentation of WM.
Figure 1. Figure 1.

Disclosures
No relevant conflicts of interest to declare.