Integrated Epigenetic and Transcriptional Single Cell Analysis of t(11;14) Myeloma and Its BCL2 Dependency
Bibliographic record
Abstract
Abstract Multiple myeloma is characterized by recurrent chromosomal translocations that involve the immunoglobulin gene enhancers and partners such as the cyclin D genes (CCND1, CCND2 or CCND3) or other genes like WHSC1 and MAF. t(11;14) results in upregulation of CCND1 with unique morphological, phenotypic markers, and drug sensitivity profiles with exquisite sensitivity to BCL2 inhibitors. It is evident now that this unique sensitivity profile is driven by the BH3-proapoptotic protein priming of BCL2 with high BCL2/MCL1 or BCL2/BCL2L1 ratios. However, the epigenetic mechanisms associated with t(11;14) and their impact on genes regulation and clinical response to venetoclax remains elusive. In the present study we compared the transcriptomics (scRNA-seq) and chromatin accessibility (scATAC-seq) of single plasma cells of MM patients with and without t(11;14) as well as pre- and post-venetoclax exposure in order to establish the epigenomic signature of t(11;14) and/or BCL2-sensitivity in myeloma. Serial BM aspirates (n=24) were collected from 15 relapsed or refractory myeloma patients (RRMM); harboring t(11;14) (n=6 pairs) and 9 without this translocation prior to initiation of salvage therapy and at time of relapse. All t(11;14) MM patients were treated with venetoclax. Unbiased chromatin accessibility and mRNA profiling of CD138 pos cells were performed using the chromium single cell ATAC and RNA-Seq 3' solution (10x Genomics), respectively. Cell Ranger, Seurat and ArchR were used for sample de-multiplexing, barcode processing, single-cell 3' gene, peaks counting, and data analysis. We first compared the scATAC-seq and scRNA-seq profiles found in CD138 pos MM cell isolated from patients harboring t(11;14) with the one obtained in patients without this translocation. As expected, t(11:14) patients had high chromatin accessibility at the CCND1 locus and high mRNA expression. Differentially accessible chromatin analysis identified 147518 peaks that were specific to t(11;14) patients. Of interest, motifs enrichment analysis of accessible peaks identified a "B cell-like" motifs signature with enriched TFs motifs such as TCF4 and PAX5 in t(11;14) patients compared to non t(11:14) enriched for IRF and STAT family of motifs. The integration of the scATAC-seq and scRNA-seq data confirmed the B cell signature of t(11;14) patients with upregulation of B cell markers such as MS4A1, VPREB3, CD79A, CD19, and down-regulation of plasma cell markers such as TDO2, EFEMP1, CD28, SLAMF7, and IL6R. Additionally, we found PAX1, PAX5, TCF3, TCF5, and SPI1 transcription factors to be highly expressed in t(11;14) while the non t(11:14) were enriched for IRF1-9 transcription factors. Of interest, the clustering analysis performed on scATAC-seq data identified 3 non t(11;14) patients with a chromatin accessibility profile similar to that of t(11;14) patients. They expressed B cell markers (PAX5, VPREB3 or FCRLA), overexpressed BCL2 and we are currently examining whether this B cell-like epigenetic signature determines sensitivity to venetoclax. In order to define the epigenetic contribution to the acquired resistance to venetoclax in t(11;14) myeloma, we compared the chromatin accessibility profiles of t(11;14) patients pre- vs. post-venetoclax treatment. Enriched motifs within accessible peaks differed significantly between pre- and post-venetoclax with RELA, REL, RELB and EGR1 motifs predominantly enrichmed in the pre-samples in contrast to JUN, JUNB, JUND and FOSL1/L2 motifs enrichment in the post-samples. Of note, integration analysis of scRNAseq (differentially expressed genes) and ATACseq data (differentially accessible peaks) identified MCL1 and ENSA (a gene 60 Kb centromeric to MCL1 on chr1q) as the top enriched genes and peaks in resistant samples suggesting that copy number gain at the MCL1 locus (which we confirmed by single cell CNV analysis) rather than epigenetic modifications is likely the main determinant of acquired resistant to venetoclax in t(11;14) MM. In the current study we have defined the epigenetic regulome and transcriptome associated with t(11;14) myeloma and its relatedness to B cell rather than plasma cell biology. Our studies also suggest that acquired resistance to venetoclax is largely driven by copy number gain at the MCL1 locus. Disclosures Bahlis: Karyopharm: Consultancy, Honoraria; Amgen: Consultancy, Honoraria; Takeda: Consultancy, Honoraria; Sanofi: Consultancy, Honoraria; Abbvie: Consultancy, Honoraria; Janssen: Consultancy, Honoraria; BMS/Celgene: Consultancy, Honoraria; GlaxoSmithKline: Consultancy, Honoraria; Genentech: Consultancy; Pfizer: Consultancy, Honoraria. Neri: BMS: Consultancy, Honoraria; Janssen: Consultancy, Honoraria; Sanofi: Consultancy, Honoraria; Amgen: Consultancy, Honoraria.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.000 |
| Insufficient payload (model declined to judge) | 0.002 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".