A Comprehensive Analysis of Single-Cell Chromatin Accessibility and Gene Expression Identifies Intra-Tumor Heterogeneity and Molecular Treatment Responses in Relapsed/Refractory Multiple Myeloma

Introduction: Multiple myeloma (MM) is a heterogeneous malignancy of clonal plasma cells that accumulate in the bone marrow (BM). Despite new treatment approaches, in most patients resistant subclones are selected by therapy, resulting in the development of refractory disease. While the subclonal architecture in newly diagnosed patients has been investigated in great detail, intra-tumor heterogeneity in relapsed/refractory (RR) MM is poorly characterized. Recent technological and computational advances provide the opportunity to systematically analyze tumor samples at single-cell (sc) level with high accuracy and througput. Here, we present a pilot study for an integrative analysis of sc Assay for Transposase-Accessible Chromatin with high-throughput sequencing (scATAC-seq) and scRNA-seq with the aim to comprehensively study the regulatory landscape, gene expression, and evolution of individual subclones in RRMM patients. Methods: We have included 20 RRMM patients with longitudinally collected paired BM samples. scATAC- and scRNA-seq data were generated using the 10X Genomics platform. Pre-processing of the sc-seq data was performed with the CellRanger software (reference genome GRCh38). For downstream analyses the R-packages Seurat and Signac (Satija Lab) as well as Cicero (Trapnell Lab) were used. For all patients bulk whole genome sequencing (WGS) data was available, which we used for confirmatory studies of intra-tumor heterogeneity. Results: A comprehensive study at the sc level requires extensive quality controls (QC). All scATAC-seq files passed the QC, including the detected number of cells, number of fragments in peaks or the ratio of mononucleosomal to nucleosome-free fragments. Yet, unsupervised clustering of the differentially accessible regions resulted in two main clusters, strongly associated with sample processing time. Delay of sample processing by 1-2 days, e.g. due to shipment from participating centers, resulted in global change of chromatin accessibility with more than 10,000 regions showing differences compared to directly processed samples. The corresponding scRNA-seq files also consistently failed QC, including detectable genes per cell and the percentage of mitochondrial RNA. We excluded these samples from the study. Analysing scATAC-seq data, we observed distinct clusters before and after treatment of RRMM, indicating clonal adaptation or selection in all samples. Treatment with carfilzomib resulted in highly increased co-accessib