Braggio, Esteban

Abstract In multiple myeloma (MM), despite well-characterized precursor states such as monoclonal gammopathy of undetermined significance (MGUS) and smoldering multiple myeloma (SMM), there is a lack of sufficient biomarkers to predict disease progression. Most genomic analyses have studied the malignant plasma cells, however, cancers form a complex ecosystem with the immune and stromal microenvironment. To characterize the cellular composition and transcriptional programs of each component of the tumor and microenvironment at different stages of MM progression, we employed single-cell RNA sequencing on 48K plasma and 40.8K immune microenvironmental cells from a cohort of 22 patients with varying stages of disease progression and 9 healthy donors. Expression profiles of plasma cells revealed clear tumor-specific differences in known oncogenic drivers in MM (MMSET/FGFR3, CCND1 and MAFB) as well as other clonally expressed genes (LAMP5, HIST1H1C, and AREG), distinguishing them from healthy plasma cells. We identified a subset of cycling plasma cells in malignant samples, observing a range of proliferative capacity across disease stages. Furthermore, our approach allowed a unique head-to-head comparison of gene expression changes in normal and malignant plasma cells from the same individual, revealing early alterations in genes related to immune modulation (NKBIA) or controlling transcription and differentiation (EID1). Some alterations were patient-specific, while others, such as MHC I overexpression and CD27 loss, were recurrently observed across subsets of the cohort. Analysis of the BM microenvironment demonstrated significant infiltration of natural killer cells, non-classical monocytes/macrophages, and T cells, even in the earliest stages of the disease. Further investigation revealed upregulation of MHC II expression at the mRNA level in CD14+ monocytes/macrophages and yet, intriguingly, analysis by CyTOF and immunohistochemistry revealed a shift towards intracellular localization of MHC II in these cells. Co-culture with MM cell lines was sufficient to induce the decrease of extracellular MHC II, providing strong evidence for MM-induced compromised antigen presentation by macrophages, and hinting at a mechanism of immune evasion. Together, our results provide a comprehensive view at the complex interplay of the immune and malignant cells in different stages of the disease. We demonstrate the immune response beginning in premalignant conditions to be heterogeneous, including compromised antigen presentation as well as alterations in cellular composition and signaling. Consideration of the type of immunological response may prove valuable in determination of progression risk, as well as open up potential strategies for therapy. Citation Format: Nicholas J. Haradhvala, Oksana Zavidij, Tarek H. Mouhieddine, Romanos Sklavenitis-Pistofidis, Jihye Park, Mairead Reidy, Abdallah Flaifel, Benjamin Ferland, Salomon Manier, Mark Bustoros, Daisy Huynh, Marzia Capelletti, Brianna Berrios, Mahshid Rahmat, Chia-Jen Liu, Meng Xiao He, Esteban Braggio, Rafael Fonseca, Yosef Maruvka, Jennifer Guerriero, Melissa Goldman, Eliezer Van Allen, Steven McCarroll, Jamil Azzi, Gad Getz, Irene M. Ghobrial. Single-cell RNA sequencing reveals compromised immune microenvironment in precursor stages of multiple myeloma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 139.

293 Background: Cholangiocarcinoma (CCA) is a deadly cancer with ~20% five year survival. Obtaining tissue to confirm diagnosis of CCA is often challenging. This also limits the use of tumor genotyping that could guide potential molecularly targeted treatment (such as FGFR inhibitors for a subset of cases). We investigated whether circulating tumor DNA (ctDNA) analysis in plasma is feasible for noninvasive tumor genotyping in patients with advanced CCA. Methods: We collected and processed plasma samples for ctDNA analysis from 17 patients with CCA at presentation and from 11 healthy volunteers. Following cell-free DNA extraction and quality assessment using digital PCR, we performed low-pass, shallow whole genome sequencing (sWGS). We analyzed sequencing data to identify copy number aberrations and quantify tumor fraction in plasma using recently published computational approaches. Results: Median total cfDNA concentration in patients with CCA was 7,003 genome equivalents (GEs)/mL plasma (range: 1,512- 107,707 GEs/mL) and significantly higher than healthy volunteers (two tailed t-test p < 0.04, median: 1,658 GEs/mL, range: 661-4,938 GEs/mL). We successfully sequenced all tumor libraries, generating a median of 173 million paired-end reads. ctDNA was detectable in 10/17 samples using sWGS. In these 10 samples, median tumor fraction in plasma DNA was 7.1% (range: 3.9%-44.6%). We found recurrent copy number aberrations expected in CCA including 1q gain in 80%, 7p gain in 50%, 1p loss in 70% and 3p loss in 50% of cases. Conclusions: ctDNA analysis for noninvasive tumor genotyping is feasible in patients with advanced cholangiocarcinoma. On-going work is focused on expanding these results using targeted and whole-exome sequencing.

The prevalence, clinicopathologic correlates, and outcomes of previously untreated chronic lymphocytic leukemia (CLL) patients with IGH‐BCL2 and IGH‐BCL3 translocations are not well known. Using the Mayo Clinic CLL database, we identified patients seen between March 1, 2002 and September 30, 2016 who had FISH testing performed within 3 years of CLL diagnosis. The prognostic profile, time to first therapy (TTT), and overall survival (OS) of patients with IGH‐BCL2 and IGH‐BCL3 translocation were compared to patients without these abnormalities (non‐IGH group). Of 1684 patients who met the inclusion criteria, 38 (2.2%) had IGH‐BCL2, and 16 (0.9%) had IGH‐BCL3 translocation at diagnosis. Patients with IGH‐BCL3 translocation were more likely to have high and very‐high CLL‐International Prognostic Index, compared to patients with IGH‐BCL2 translocation and the non‐IGH group. The 5‐year probability of requiring therapy was significantly higher for IGH‐BCL3 compared to IGH‐BCL2 and non‐IGH groups (84% vs 33% vs 29%, respectively, P < 0.0001). The 5‐year OS was significantly shorter for IGH‐BCL3 compared to IGH‐BCL2 and non‐IGH groups (45% vs 89% vs 86%, respectively, P < 0.0001). On multivariable analyses, IGH‐BCL3 translocation was associated with a shorter TTT (hazard ratio [HR] = 2.7; P = 0.005) and shorter OS (HR = 5.5; P < 0.0001); IGH‐BCL2 translocation did not impact TTT and OS. In conclusion, approximately 3% of all newly diagnosed CLL patients have either an IGH‐BCL2 or IGH‐BCL3 translocation. Patients with IGH‐BCL3 translocations have a distinct prognostic profile and outcome. These results support the inclusion of an IGH probe during the routine evaluation of FISH abnormalities in newly diagnosed CLL.

Abstract CLL is a highly heritable cancer with first degree relatives of CLL cases having a 7.5-fold increased CLL risk. Genome-wide association studies (GWAS) and linkage studies have been performed to study inherited predisposition; however a larger proportion of heritability to CLL remains unexplained. Rare coding variants might account for the missing heritability information. Inherited loss of function variants in shelterin complex genes (POT1, ACD, TERF1, TINF2, TERF2, TERF2IP- involved in telomere regulation), CDK1 (critical for cell division) and ATM (tumor suppressor gene) have been found to co-segregate in CLL families and be enriched in CLL cases using exome-wide sequencing data. Our study evaluates rare germline variants from these suspect genes segregating in CLL families who are followed at the Mayo Clinic. Using whole exome sequencing (WES), we sequenced 93 CLL families with at least 2 reported CLL cases consisting of 443 individuals: 160 with CLL, 73 with monoclonal B-cell lymphocytosis (MBL), and 210 relatives. DNA was extracted from buccal cells, coding exons were selectively captured using Agilent 50Mb and SureSelect Human All Exon V4 capture kits; sequencing was performed using Illumina HiSeq 2000. Mayo Clinic's DNASeq pipeline uses Novoalign (initial read alignment), Picard (marking duplicate reads), and the Genome Analysis Toolkit (GATK) for local realignment, recalibration, and variant calling. The variant discovery step leverages GATK's HaplotypeCaller in per sample mode and all of the samples across the cohort are jointly genotyped together. All called variants are evaluated with GATK's Variant Quality Score Recalibration tool and annotated for biological relevance (BIOR). Quality control included removing variants that had &lt;75% call rate across the two capture kits, &lt;8x coverage, or phred score&lt;10, resulting in 317,666 remaining variants. Of these, over 80% of the coding sequence had a median read depth of 23 reads. In our pedigrees, we searched for rare variants within the genes described above. We identified suspect variants with the following criteria: 1) enriched in CLL and MBL samples compared to unaffected samples; 2) multiple affected members with the variant within a family; 3) variants present in all sequenced affecteds within the family; 4) rarely seen in an in-house database of non-cancer controls or 1K Genomes; and 5) predicted to have a functional damaging effect (using SIFT). We identified three novel rare missense variants, defined as functionally deleterious, which each co-segregated within a CLL family. Specifically, these variants from shelterin complex genes; POT1 (rs116916706), TERF2IP (rs138458227), and TERF2 (rs749171225), met the criteria. This study further highlights telomere dysregulation as a key process in CLL development. Investigating rare variants within CLL pedigrees with WES can help identify germline variants impacting predisposition to familial CLL. Citation Format: Alyssa I. Clay-Gilmour, Daniel R. O'Brien, Sara J. Achenbach, Celine M. Vachon, Kari G. Chaffee, Timothy G. Call, Jose F. Leis, Aaron D. Norman, Brian F. Kabat, Sameer A. Parikh, Neil E. Kay, Esteban Braggio, James R. Cerhan, Susan L. Slager. Rare germline variants segregating in chronic lymphocytic leukemia (CLL) families [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1226.

Abstract Richter syndrome represents the evolution of chronic lymphocytic leukemia into an aggressive tumor, most commonly diffuse large B-cell lymphoma. The lack of in vitro and in vivo models has severely hampered drug testing in a disease that is poorly responsive to common chemoimmunotherapeutic combinations as well as to novel kinase inhibitors. Here we report for the first time the establishment and genomic characterization of two patient-derived tumor xenograft (PDX) models of Richter syndrome, RS9737 and RS1316. Richter syndrome xenografts were genetically, morphologically, and phenotypically stable and similar to the corresponding primary tumor. RS9737 was characterized by biallelic inactivation of CDKN2A and TP53, monoallelic deletion of 11q23 (ATM), and mutations of BTK, KRAS, EGR2, and NOTCH1. RS1316 carried trisomy 12 and showed mutations in BTK, KRAS, MED12, and NOTCH2. RNA sequencing confirmed that in both cases &gt;80% of the transcriptome was shared between primary tumor and PDX. In line with the mutational profile, pathway analyses revealed overactivation of the B-cell receptor, NFκB, and NOTCH pathways in both tumors, potentially providing novel tumor targets. In conclusion, these two novel models of Richter syndrome represent useful tools to study biology and response to therapies of this highly aggressive and still incurable tumor. Significance: Two patient-derived xenograft models of Richter syndrome represent the first ex vivo model to study biology of the disease and to test novel therapeutic strategies. Cancer Res; 78(13); 3413–20. ©2018 AACR.

Progressive multifocal leukoencephalopathy rarely occurs in patients with multiple myeloma. Intracranial central nervous system invasion is also an uncommon event in multiple myeloma, occurring in less than 1% of cases. We describe herein an exceptional case of coexisting progressive multifocal leukoencephalopathy and intraparenchymal central nervous system myeloma infiltration. A 73-year-old woman with relapsed multiple myeloma was treated with 15 cycles of lenalidomide and dexamethasone, but therapy had to be stopped because of a hip fracture after a fall. During hospitalization, the patient developed progressive multifocal leukoencephalopathy caused by John Cunningham virus, and a prominent intra-parenchymal CD138-positive infiltrate was detected. VDJ rearrangements of the immunoglobulin heavy chain gene and the mutational profile of plasma cells in bone marrow at the time of diagnosis and in brain biopsy after progression were analyzed by next generation sequencing, showing genetic differences between medullary and extramedullary myeloma cells. The role of long-term treatment with lenalidomide and dexamethasone in the development progressive multifocal leukoencephalopathy or intraparenchymal central nervous system myeloma infiltration remains unknown. However, our results suggest that both events may have arisen as a consequence of treatment-related immunosuppression. Thus, an appropriate clinical approach compatible with the simultaneous treatment of progressive multifocal leukoencephalopathy and multiple myeloma should be developed.

Abstract Bortezomib is highly effective in the treatment of multiple myeloma; however, emergent drug resistance is common. Consequently, we employed CRISPR targeting 19,052 human genes to identify unbiased targets that contribute to bortezomib resistance. Specifically, we engineered an RPMI8226 multiple myeloma cell line to express Cas9 infected by lentiviral vector CRISPR library and cultured derived cells in doses of bortezomib lethal to parental cells. Sequencing was performed on surviving cells to identify inactivated genes responsible for drug resistance. From two independent whole-genome screens, we selected 31 candidate genes and constructed a second CRISPR sgRNA library, specifically targeting each of these 31 genes with four sgRNAs. After secondary screening for bortezomib resistance, the top 20 “resistance” genes were selected for individual validation. Of these 20 targets, the proteasome regulatory subunit PSMC6 was the only gene validated to reproducibly confer bortezomib resistance. We confirmed that inhibition of chymotrypsin-like proteasome activity by bortezomib was significantly reduced in cells lacking PSMC6. We individually investigated other members of the PSMC group (PSMC1 to 5) and found that deficiency in each of those subunits also imparts bortezomib resistance. We found 36 mutations in 19S proteasome subunits out of 895 patients in the IA10 release of the CoMMpass study (https://themmrf.org). Our findings demonstrate that the PSMC6 subunit is the most prominent target required for bortezomib sensitivity in multiple myeloma cells and should be examined in drug-refractory populations. Mol Cancer Ther; 16(12); 2862–70. ©2017 AACR.

Abstract FAM46C is one of the most recurrently mutated genes in multiple myeloma; however its role in disease pathogenesis has not been determined. Here we demonstrate that wild-type (WT) FAM46C overexpression induces substantial cytotoxicity in multiple myeloma cells. In contrast, FAM46C mutations found in multiple myeloma patients abrogate this cytotoxicity, indicating a survival advantage conferred by the FAM46C mutant phenotype. WT FAM46C overexpression downregulated IRF4, CEBPB, and MYC and upregulated immunoglobulin (Ig) light chain and HSPA5/BIP. Furthermore, pathway analysis suggests that enforced FAM46C expression activated the unfolded protein response pathway and induced mitochondrial dysfunction. CRISPR-mediated depletion of endogenous FAM46C enhanced multiple myeloma cell growth, decreased Ig light chain and HSPA5/BIP expression, activated ERK and antiapoptotic signaling, and conferred relative resistance to dexamethasone and lenalidomide treatments. Genes altered in FAM46C-depleted cells were enriched for signaling pathways regulating estrogen, glucocorticoid, B-cell receptor signaling, and ATM signaling. Together these results implicate FAM46C in myeloma cell growth and survival and identify FAM46C mutation as a contributor to myeloma pathogenesis and disease progression via perturbation in plasma cell differentiation and endoplasmic reticulum homeostasis. Cancer Res; 77(16); 4317–27. ©2017 AACR.