1q amplification and PHF19 expressing high-risk cells are associated with relapsed/refractory multiple myeloma.

Multiple Myeloma is an incurable plasma cell malignancy with a poor survival rate that is usually treated with immunomodulatory drugs (iMiDs) and proteosome inhibitors (PIs). The malignant plasma cells quickly become resistant to these agents causing relapse and uncontrolled growth of resistant clones. From whole genome sequencing (WGS) and RNA sequencing (RNA-seq) studies, different high-risk translocation, copy number, mutational, and transcriptional markers can be identified. One of these markers, PHF19, epigenetically regulates cell cycle and other processes and is already studied using RNA-seq. In this study, we generate a large (325,025 cells and 49 patients) single cell multi-omic dataset and jointly quantify ATAC- and RNA-seq for each cell and matched genomic profiles for each patient. We identify an association between one plasma cell subtype with myeloma progression that we call relapsed/refractory plasma cells (RRPCs). These cells are associated with chromosome 1q alterations, TP53 mutations, and higher expression of PHF19. We also identify downstream regulation of cell cycle inhibitors in these cells, possible regulation by the transcription factor (TF) PBX1 on chromosome 1q, and determine that PHF19 may be acting primarily through this subset of cells.

Corrigendum: Protein arginine methyltransferase 1 is a therapeutic vulnerability in multiple myeloma.

[This corrects the article DOI: 10.3389/fimmu.2023.1239614.].

Identifying novel mechanisms of biallelic TP53 loss refines poor outcome for patients with multiple myeloma.

Biallelic TP53 inactivation is the most important high-risk factor associated with poor survival in multiple myeloma. Classical biallelic TP53 inactivation has been defined as simultaneous mutation and copy number loss in most studies; however, numerous studies have demonstrated that other factors could lead to the inactivation of TP53. Here, we hypothesized that novel biallelic TP53 inactivated samples existed in the multiple myeloma population. A random forest regression model that exploited an expression signature of 16 differentially expressed genes between classical biallelic TP53 and TP53 wild-type samples was subsequently established and used to identify novel biallelic TP53 samples from monoallelic TP53 groups. The model reflected high accuracy and robust performance in newly diagnosed relapsed and refractory populations. Patient survival of classical and novel biallelic TP53 samples was consistently much worse than those with mono-allelic or wild-type TP53 status. We also demonstrated that some predicted biallelic TP53 samples simultaneously had copy number loss and aberrant splicing, resulting in overexpression of high-risk transcript variants, leading to biallelic inactivation. We discovered that splice site mutation and overexpression of the splicing factor MED18 were reasons for aberrant splicing. Taken together, our study unveiled the complex transcriptome of TP53, some of which might benefit future studies targeting abnormal TP53.

Identifying 1q amplification and PHF19 expressing high-risk cells associated with relapsed/refractory multiple myeloma.

Multiple Myeloma is an incurable plasma cell malignancy with a poor survival rate that is usually treated with immunomodulatory drugs (iMiDs) and proteosome inhibitors (PIs). The malignant plasma cells quickly become resistant to these agents causing relapse and uncontrolled growth of resistant clones. From whole genome sequencing (WGS) and RNA sequencing (RNA-seq) studies, different high-risk translocation, copy number, mutational, and transcriptional markers have been identified. One of these markers, PHF19, epigenetically regulates cell cycle and other processes and has already been studied using RNA-seq. In this study a massive (325,025 cells and 49 patients) single cell multiomic dataset was generated with jointly quantified ATAC- and RNA-seq for each cell and matched genomic profiles for each patient. We identified an association between one plasma cell subtype with myeloma progression that we have called relapsed/refractory plasma cells (RRPCs). These cells are associated with 1q alterations, TP53 mutations, and higher expression of PHF19. We also identified downstream regulation of cell cycle inhibitors in these cells, possible regulation of the transcription factor (TF) PBX1 on 1q, and determined that PHF19 may be acting primarily through this subset of cells.

Protein arginine methyltransferase 1 is a therapeutic vulnerability in multiple myeloma.

Multiple myeloma (MM) is a devastating plasma cell malignancy characterized by the expansion of aberrant monoclonal plasma cells in the bone marrow, leading to severe clinical manifestations and poor prognosis, particularly in relapsed/refractory cases. Identifying novel therapeutic targets is crucial to improve treatment outcomes in these patients. In this study, we investigated the role of the protein arginine methyltransferase 1 (PRMT1) in MM pathogenesis and explored its potential as a therapeutic target. We observed that PRMT1, responsible for most asymmetric di-methylation in cells, exhibited the highest expression among PRMT family members in MM cell lines and primary MM cells. Importantly, PRMT1 expression was significantly elevated in relapsed/refractory patients compared to newly diagnosed patients. High expression of PRMT1 expression was strongly associated with poor prognosis. We found that genetic or enzymatic inhibition of PRMT1 impaired MM cell growth, induced cell cycle arrest, and triggered cell death. Treatment with MS023, a potent PRMT type I inhibitor, demonstrated a robust inhibitory effect on the viability of primary cells isolated from newly diagnosed and proteasome inhibitor-relapsed/refractory patients in a dose-dependent manner. Suppression of PRMT1 downregulated genes related to cell division and upregulated genes associated with apoptosis pathway. We also found that genes related to immune response and lymphocyte activation were significantly upregulated in PRMT1-suppressed cells. Notably, the activation status of T cells was strikingly enhanced upon co-culturing with PRMT1-KO MM cells. In vivo studies using a xenograft model revealed that targeting PRMT1 by either CRISPR/Cas9-mediated knockout or MS023 treatment significantly attenuated MM tumor growth and prolonged the survival of tumor-bearing mice. Histological analysis further confirmed increased apoptotic cell death in MS023-treated tumors. Collectively, our findings establish PRMT1 as an indispensable and novel therapeutic vulnerability in MM. The elevated expression of PRMT1 in relapsed/refractory patients underscores its potential as a target for overcoming treatment resistance. Moreover, our results highlight the efficacy of MS023 as a promising therapeutic agent against MM, offering new avenues for therapeutic approaches in relapsed/refractory MM.

Multiomic Mapping of Acquired Chromosome 1 Copy-Number and Structural Variants to Identify Therapeutic Vulnerabilities in Multiple Myeloma.

PURPOSE: Chromosome 1 (chr1) copy-number abnormalities (CNA) and structural variants (SV) are frequent in newly diagnosed multiple myeloma (NDMM) and are associated with a heterogeneous impact on outcomes, the drivers of which are largely unknown. EXPERIMENTAL DESIGN: A multiomic approach comprising CRISPR, gene mapping of CNAs and SVs, methylation, expression, and mutational analysis was used to document the extent of chr1 molecular variants and their impact on pathway utilization. RESULTS: We identified two distinct groups of gain(1q): focal gains associated with limited gene-expression changes and a neutral prognosis, and whole-arm gains, which are associated with substantial gene-expression changes, complex genetics, and an adverse prognosis. CRISPR identified a number of dependencies on chr1 but only limited variants associated with acquired CNAs. We identified seven regions of deletion, nine of gain, three of chromothripsis (CT), and two of templated insertion (TI), which contain a number of potential drivers. An additional mechanism involving hypomethylation of genes at 1q may contribute to the aberrant gene expression of a number of genes. Expression changes associated with whole-arm gains were substantial and gene set enrichment analysis identified metabolic processes, apoptotic resistance, signaling via the MAPK pathway, and upregulation of transcription factors as being key drivers of the adverse prognosis associated with these variants. CONCLUSIONS: Multiple layers of genetic complexity impact the phenotype associated with CNAs on chr1 to generate its associated clinical phenotype. Whole-arm gains of 1q are the critically important prognostic group that deregulate multiple pathways, which may offer therapeutic vulnerabilities.

Alternative splicing in multiple myeloma is associated with the non-homologous end joining pathway.

Alternative splicing plays a pivotal role in tumorigenesis and proliferation. However, its pattern and pathogenic role has not been systematically analyzed in multiple myeloma or its subtypes. Alternative splicing profiles for 598 newly diagnosed myeloma patients with comprehensive genomic annotation identified primary translocations, 1q amplification, and DIS3 events to have more differentially spliced events than those without. Splicing levels were correlated with expression of splicing factors. Moreover, the non-homologous end joining pathway was an independent factor that was highly associated with splicing frequency as well as an increased number of structural variants. We therefore identify an axis of high-risk disease encompassing expression of the non-homologous end joining pathway, increase structural variants, and increased alternative splicing that are linked together. This indicates a joint pathogenic role for DNA damage response and alternative RNA processing in myeloma.

Molecular characterization and clinical outcome of B-cell precursor acute lymphoblastic leukemia with IG-MYC rearrangement.

Rarely, immunophenotypically immature B-cell precursor acute lymphoblastic leukemia (BCP-ALL) carries an immunoglobulin- MYC rearrangement (IG-MYC-r). This can result in diagnostic confusion with Burkitt lymphoma/leukemia and use of individualized treatment schedules of unproven efficacy. Here we compare the molecular characteristics of these conditions and investigate historic clinical outcome data. We identified 90 cases registered in a national BCP-ALL clinical trial/registry. When present, diagnostic material underwent cytogenetic, exome, methylome and transcriptome analyses. The outcomes analyzed were 3-year event-free survival and overall survival. IG-MYC-r was identified in diverse cytogenetic backgrounds, co-existing with either established BCP-ALL-specific abnormalities (high hyperdiploidy, n=3; KMT2A-rearrangement, n=6; iAMP21, n=1; BCR-ABL1, n=1); BCL2/BCL6-rearrangements (n=15); or, most commonly, as the only defining feature (n=64). Within this final group, precursor-like V(D)J breakpoints predominated (8/9) and KRAS mutations were common (5/11). DNA methylation identified a cluster of V(D)J-rearranged cases, clearly distinct from Burkitt leukemia/lymphoma. Children with IG-MYC-r within that subgroup had a 3-year event-free survival of 47% and overall survival of 60%, representing a high-risk BCP-ALL. To develop effective management strategies this group of patients must be allowed access to contemporary, minimal residual disease-adapted, prospective clinical trial protocols.

The spatio-temporal evolution of multiple myeloma from baseline to relapse-refractory states.

Deciphering Multiple Myeloma evolution in the whole bone marrow is key to inform curative strategies. Here, we perform spatial-longitudinal whole-exome sequencing, including 140 samples collected from 24 Multiple Myeloma patients during up to 14 years. Applying imaging-guided sampling we observe three evolutionary patterns, including relapse driven by a single-cell expansion, competing/co-existing sub-clones, and unique sub-clones at distinct locations. While we do not find the unique relapse sub-clone in the baseline focal lesion(s), we show a close phylogenetic relationship between baseline focal lesions and relapse disease, highlighting focal lesions as hotspots of tumor evolution. In patients with ≥3 focal lesions on positron-emission-tomography at diagnosis, relapse is driven by multiple distinct sub-clones, whereas in other patients, a single-cell expansion is typically seen (p < 0.01). Notably, we observe resistant sub-clones that can be hidden over years, suggesting that a prerequisite for curative therapies would be to overcome not only tumor heterogeneity but also dormancy.

High-resolution simulations of chromatin folding at genomic rearrangements in malignant B cells provide mechanistic insights into proto-oncogene deregulation.

Genomic rearrangements are known to result in proto-oncogene deregulation in many cancers, but the link to 3D genome structure remains poorly understood. Here, we used the highly predictive heteromorphic polymer (HiP-HoP) model to predict chromatin conformations at the proto-oncogene CCND1 in healthy and malignant B cells. After confirming that the model gives good predictions of Hi-C data for the nonmalignant human B cell-derived cell line GM12878, we generated predictions for two cancer cell lines, U266 and Z-138. These possess genome rearrangements involving CCND1 and the immunoglobulin heavy locus (IGH), which we mapped using targeted genome sequencing. Our simulations showed that a rearrangement in U266 cells where a single IGH super-enhancer is inserted next to CCND1 leaves the local topologically associated domain (TAD) structure intact. We also observed extensive changes in enhancer-promoter interactions within the TAD, suggesting that it is the downstream chromatin remodeling which gives rise to the oncogene activation, rather than the presence of the inserted super-enhancer DNA sequence per se. Simulations of the IGH-CCND1 reciprocal translocation in Z-138 cells revealed that an oncogenic fusion TAD is created, encompassing CCND1 and the IGH super-enhancers. We predicted how the structure and expression of CCND1 changes in these different cell lines, validating this using qPCR and fluorescence in situ hybridization microscopy. Our work demonstrates the power of polymer simulations to predict differences in chromatin interactions and gene expression for different translocation breakpoints.

Perspectives on the Risk-Stratified Treatment of Multiple Myeloma.

The multiple myeloma treatment landscape has changed dramatically. This change, paralleled by an increase in scientific knowledge, has resulted in significant improvement in survival. However, heterogeneity remains in clinical outcomes, with a proportion of patients not benefiting from current approaches and continuing to have a poor prognosis. A significant proportion of the variability in outcome can be predicted on the basis of clinical and biochemical parameters and tumor-acquired genetic variants, allowing for risk stratification and a more personalized approach to therapy. This article discusses the principles that can enable the rational and effective development of therapeutic approaches for high-risk multiple myeloma.

Genetic subtypes of smoldering multiple myeloma are associated with distinct pathogenic phenotypes and clinical outcomes.

Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with significant heterogeneity in disease progression. Existing clinical models of progression risk do not fully capture this heterogeneity. Here we integrate 42 genetic alterations from 214 SMM patients using unsupervised binary matrix factorization (BMF) clustering and identify six distinct genetic subtypes. These subtypes are differentially associated with established MM-related RNA signatures, oncogenic and immune transcriptional profiles, and evolving clinical biomarkers. Three genetic subtypes are associated with increased risk of progression to active MM in both the primary and validation cohorts, indicating they can be used to better predict high and low-risk patients within the currently used clinical risk stratification models.

Structural variants shape the genomic landscape and clinical outcome of multiple myeloma.

Deciphering genomic architecture is key to identifying novel disease drivers and understanding the mechanisms underlying myeloma initiation and progression. In this work, using the CoMMpass dataset, we show that structural variants (SV) occur in a nonrandom fashion throughout the genome with an increased frequency in the t(4;14), RB1, or TP53 mutated cases and reduced frequency in t(11;14) cases. By mapping sites of chromosomal rearrangements to topologically associated domains and identifying significantly upregulated genes by RNAseq we identify both predicted and novel putative driver genes. These data highlight the heterogeneity of transcriptional dysregulation occurring as a consequence of both the canonical and novel structural variants. Further, it shows that the complex rearrangements chromoplexy, chromothripsis and templated insertions are common in MM with each variant having its own distinct frequency and impact on clinical outcome. Chromothripsis is associated with a significant independent negative impact on clinical outcome in newly diagnosed cases consistent with its use alongside other clinical and genetic risk factors to identify prognosis.

Myeloma Genome Project Panel is a Comprehensive Targeted Genomics Panel for Molecular Profiling of Patients with Multiple Myeloma.

PURPOSE: We designed a comprehensive multiple myeloma targeted sequencing panel to identify common genomic abnormalities in a single assay and validated it against known standards. EXPERIMENTAL DESIGN: The panel comprised 228 genes/exons for mutations, 6 regions for translocations, and 56 regions for copy number abnormalities (CNA). Toward panel validation, targeted sequencing was conducted on 233 patient samples and further validated using clinical FISH (translocations), multiplex ligation probe analysis (MLPA; CNAs), whole-genome sequencing (WGS; CNAs, mutations, translocations), or droplet digital PCR (ddPCR) of known standards (mutations). RESULTS: Canonical immunoglobulin heavy chain translocations were detected in 43.2% of patients by sequencing, and aligned with FISH except for 1 patient. CNAs determined by sequencing and MLPA for 22 regions were comparable in 103 samples and concordance between platforms was R2 = 0.969. Variant allele frequency (VAF) for 74 mutations were compared between sequencing and ddPCR with concordance of R2 = 0.9849. CONCLUSIONS: In summary, we have developed a targeted sequencing panel that is as robust or superior to FISH and WGS. This molecular panel is cost-effective, comprehensive, clinically actionable, and can be routinely deployed to assist risk stratification at diagnosis or posttreatment to guide sequencing of therapies.

Diagnostic Evidence GAuge of Single cells (DEGAS): a flexible deep transfer learning framework for prioritizing cells in relation to disease.

We propose DEGAS (Diagnostic Evidence GAuge of Single cells), a novel deep transfer learning framework, to transfer disease information from patients to cells. We call such transferrable information "impressions," which allow individual cells to be associated with disease attributes like diagnosis, prognosis, and response to therapy. Using simulated data and ten diverse single-cell and patient bulk tissue transcriptomic datasets from glioblastoma multiforme (GBM), Alzheimer's disease (AD), and multiple myeloma (MM), we demonstrate the feasibility, flexibility, and broad applications of the DEGAS framework. DEGAS analysis on myeloma single-cell transcriptomics identified PHF19high myeloma cells associated with progression. Availability: https://github.com/tsteelejohnson91/DEGAS .

Chromothripsis as a pathogenic driver of multiple myeloma.

Analysis of the genetic basis for multiple myeloma (MM) has informed many of our current concepts of the biology that underlies disease initiation and progression. Studying these events in further detail is predicted to deliver important insights into its pathogenesis, prognosis and treatment. Information from whole genome sequencing of structural variation is revealing the role of these events as drivers of MM. In particular, we discuss how the insights we have gained from studying chromothripsis suggest that it can be used to provide information on disease initiation and that, as a consequence, it can be used for the clinical classification of myeloma precursor diseases allowing for early intervention and prognostic determination. For newly diagnosed MM, the integration of information on the presence of chromothripsis has the potential to significantly enhance current risk prediction strategies and to better characterize patients with high-risk disease biology. In this article we summarize the genetic basis for MM and the role played by chromothripsis as a critical pathogenic factor active at early disease phases.

Plasma cells expression from smouldering myeloma to myeloma reveals the importance of the PRC2 complex, cell cycle progression, and the divergent evolutionary pathways within the different molecular subgroups.

Sequencing studies have shed some light on the pathogenesis of progression from smouldering multiple myeloma (SMM) and symptomatic multiple myeloma (MM). Given the scarcity of smouldering samples, little data are available to determine which translational programmes are dysregulated and whether the mechanisms of progression are uniform across the main molecular subgroups. In this work, we investigated 223 SMM and 1348 MM samples from the University of Arkansas for Medical Sciences (UAMS) for which we had gene expression profiling (GEP). Patients were analysed by TC-7 subgroup for gene expression changes between SMM and MM. Among the commonly dysregulated genes in each subgroup, PHF19 and EZH2 highlight the importance of the PRC2.1 complex. We show that subgroup specific differences exist even at the SMM stage of disease with different biological features driving progression within each TC molecular subgroup. These data suggest that MMSET SMM has already transformed, but that the other precursor diseases are distinct clinical entities from their symptomatic counterpart.

Impact of Etiological Cytogenetic Abnormalities on the Depth of Immunoparesis and Survival in Newly Diagnosed Multiple Myeloma.

INTRODUCTION/BACKGROUND: Immunoparesis, or low polyclonal immunoglobulin levels, is commonly seen in multiple myeloma (MM), and is associated with poor clinical outcomes. MM can be divided into subgroups with distinct biology and outcomes based on etiologic cytogenetic abnormalities. These include hyperdiploidy and translocations of t(11;14), t(4;14), t(14;16), and t(14;20), with the latter 3 associated with high-risk disease. We hypothesized that the different etiologic cytogenetic abnormalities drive bone marrow microenvironmental changes, resulting in different degrees of immunoparesis, and subgroup-dependent effects on clinical outcomes. MATERIALS AND METHODS: We performed a retrospective review of 985 newly diagnosed patients enrolled in the Myeloma IX and XI trials. Immunoglobulin levels, survival outcomes, and infection rates were evaluated for each cytogenetic subgroup. RESULTS: A significant proportion of patients with high-risk t(4;14), t(14;16), or t(14;20) had suppressed polyclonal immunoglobulins compared to standard-risk patients with hyperdiploidy or t(11;14). The clinical impact of immunoparesis depended on the cytogenetic subgroup, with the degree of IgM suppression effecting progression-free and overall survival only in the hyperdiploid subgroup. There was no significant difference in infection rates amongst the etiologic subgroups. CONCLUSION: These findings demonstrate that the etiologic cytogenetic subgroup influences the degree and clinical impact of immunoparesis. This suggests that the underlying cytogenetic abnormality affects remodeling of the bone marrow plasma cell niche, resulting in suppressed normal plasma cell function, and low immunoglobulin levels.