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.

Genomic Classification and Individualized Prognosis in Multiple Myeloma.

PURPOSE: Outcomes for patients with newly diagnosed multiple myeloma (NDMM) are heterogenous, with overall survival (OS) ranging from months to over 10 years. METHODS: To decipher and predict the molecular and clinical heterogeneity of NDMM, we assembled a series of 1,933 patients with available clinical, genomic, and therapeutic data. RESULTS: Leveraging a comprehensive catalog of genomic drivers, we identified 12 groups, expanding on previous gene expression-based molecular classifications. To build a model predicting individualized risk in NDMM (IRMMa), we integrated clinical, genomic, and treatment variables. To correct for time-dependent variables, including high-dose melphalan followed by autologous stem-cell transplantation (HDM-ASCT), and maintenance therapy, a multi-state model was designed. The IRMMa model accuracy was significantly higher than all comparator prognostic models, with a c-index for OS of 0.726, compared with International Staging System (ISS; 0.61), revised-ISS (0.572), and R2-ISS (0.625). Integral to model accuracy was 20 genomic features, including 1q21 gain/amp, del 1p, TP53 loss, NSD2 translocations, APOBEC mutational signatures, and copy-number signatures (reflecting the complex structural variant chromothripsis). IRMMa accuracy and superiority compared with other prognostic models were validated on 256 patients enrolled in the GMMG-HD6 (ClinicalTrials.gov identifier: NCT02495922) clinical trial. Individualized patient risks were significantly affected across the 12 genomic groups by different treatment strategies (ie, treatment variance), which was used to identify patients for whom HDM-ASCT is particularly effective versus patients for whom the impact is limited. CONCLUSION: Integrating clinical, demographic, genomic, and therapeutic data, to our knowledge, we have developed the first individualized risk-prediction model enabling personally tailored therapeutic decisions for patients with NDMM.

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.

Cystatin C is glucocorticoid responsive, directs recruitment of Trem2+ macrophages, and predicts failure of cancer immunotherapy.

Cystatin C (CyC), a secreted cysteine protease inhibitor, has unclear biological functions. Many patients exhibit elevated plasma CyC levels, particularly during glucocorticoid (GC) treatment. This study links GCs with CyC's systemic regulation by utilizing genome-wide association and structural equation modeling to determine CyC production genetics in the UK Biobank. Both CyC production and a polygenic score (PGS) capturing predisposition to CyC production were associated with increased all-cause and cancer-specific mortality. We found that the GC receptor directly targets CyC, leading to GC-responsive CyC secretion in macrophages and cancer cells. CyC-knockout tumors displayed significantly reduced growth and diminished recruitment of TREM2+ macrophages, which have been connected to cancer immunotherapy failure. Furthermore, the CyC-production PGS predicted checkpoint immunotherapy failure in 685 patients with metastatic cancer from combined clinical trial cohorts. In conclusion, CyC may act as a GC effector pathway via TREM2+ macrophage recruitment and may be a potential target for combination cancer immunotherapy.

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.

11ß-HSD1 inhibition does not affect murine tumour angiogenesis but may exert a selective effect on tumour growth by modulating inflammation and fibrosis.

Glucocorticoids inhibit angiogenesis by activating the glucocorticoid receptor. Inhibition of the glucocorticoid-activating enzyme 11ß-hydroxysteroid dehydrogenase type 1 (11ß-HSD1) reduces tissue-specific glucocorticoid action and promotes angiogenesis in murine models of myocardial infarction. Angiogenesis is important in the growth of some solid tumours. This study used murine models of squamous cell carcinoma (SCC) and pancreatic ductal adenocarcinoma (PDAC) to test the hypothesis that 11ß-HSD1 inhibition promotes angiogenesis and subsequent tumour growth. SCC or PDAC cells were injected into female FVB/N or C57BL6/J mice fed either standard diet, or diet containing the 11ß-HSD1 inhibitor UE2316. SCC tumours grew more rapidly in UE2316-treated mice, reaching a larger (P<0.01) final volume (0.158 ± 0.037 cm3) than in control mice (0.051 ± 0.007 cm3). However, PDAC tumour growth was unaffected. Immunofluorescent analysis of SCC tumours did not show differences in vessel density (CD31/alpha-smooth muscle actin) or cell proliferation (Ki67) after 11ß-HSD1 inhibition, and immunohistochemistry of SCC tumours did not show changes in inflammatory cell (CD3- or F4/80-positive) infiltration. In culture, the growth/viability (assessed by live cell imaging) of SCC cells was not affected by UE2316 or corticosterone. Second Harmonic Generation microscopy showed that UE2316 reduced Type I collagen (P<0.001), whilst RNA-sequencing revealed that multiple factors involved in the innate immune/inflammatory response were reduced in UE2316-treated SCC tumours. 11ß-HSD1 inhibition increases SCC tumour growth, likely via suppression of inflammatory/immune cell signalling and extracellular matrix deposition, but does not promote tumour angiogenesis or growth of all solid tumours.

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 ATP-binding cassette proteins ABCB1 and ABCC1 as modulators of glucocorticoid action.

Responses to hormones that act through nuclear receptors are controlled by modulating hormone concentrations not only in the circulation but also within target tissues. The role of enzymes that amplify or reduce local hormone concentrations is well established for glucocorticoid and other lipophilic hormones; moreover, transmembrane transporters have proven critical in determining tissue responses to thyroid hormones. However, there has been less consideration of the role of transmembrane transport for steroid hormones. ATP-binding cassette (ABC) proteins were first shown to influence the accumulation of glucocorticoids in cells almost three decades ago, but observations over the past 10 years suggest that differential transport propensities of both exogenous and endogenous glucocorticoids by ABCB1 and ABCC1 transporters provide a mechanism whereby different tissues are preferentially sensitive to different steroids. This Review summarizes this evidence and the new insights provided for the physiology and pharmacology of glucocorticoid action, including new approaches to glucocorticoid replacement.

Whole-genome analysis identifies novel drivers and high-risk double-hit events in relapsed/refractory myeloma.

Large-scale analyses of genomic data from patients with newly diagnosed multiple myeloma (ndMM) have been undertaken, however, large-scale analysis of relapsed/refractory MM (rrMM) has not been performed. We hypothesize that somatic variants chronicle the therapeutic exposures and clonal structure of myeloma from ndMM to rrMM stages. We generated whole-genome sequencing (WGS) data from 418 tumors (386 patients) derived from 6 rrMM clinical trials and compared them with WGS from 198 unrelated patients with ndMM in a population-based case-control fashion. We identified significantly enriched events at the rrMM stage, including drivers (DUOX2, EZH2, TP53), biallelic inactivation (TP53), noncoding mutations in bona fide drivers (TP53BP1, BLM), copy number aberrations (CNAs; 1qGain, 17pLOH), and double-hit events (Amp1q-ISS3, 1qGain-17p loss-of-heterozygosity). Mutational signature analysis identified a subclonal defective mismatch repair signature enriched in rrMM and highly active in high mutation burden tumors, a likely feature of therapy-associated expanding subclones. Further analysis focused on the association of genomic aberrations enriched at different stages of resistance to immunomodulatory agent (IMiD)-based therapy. This analysis revealed that TP53, DUOX2, 1qGain, and 17p loss-of-heterozygosity increased in prevalence from ndMM to lenalidomide resistant (LENR) to pomalidomide resistant (POMR) stages, whereas enrichment of MAML3 along with immunoglobulin lambda (IGL) and MYC translocations distinguished POM from the LEN subgroup. Genomic drivers associated with rrMM are those that confer clonal selective advantage under therapeutic pressure. Their role in therapy evasion should be further evaluated in longitudinal patient samples, to confirm these associations with the evolution of clinical resistance and to identify molecular subsets of rrMM for the development of targeted therapies.

Novel inhibitors and activity-based probes targeting serine proteases.

Serine proteases play varied and manifold roles in important biological, physiological, and pathological processes. These include viral, bacterial, and parasitic infection, allergic sensitization, tumor invasion, and metastasis. The use of activity-based profiling has been foundational in pinpointing the precise roles of serine proteases across this myriad of processes. A broad range of serine protease-targeted activity-based probe (ABP) chemotypes have been developed and we have recently introduced biotinylated and "clickable" peptides containing P1 N-alkyl glycine arginine N-hydroxy succinimidyl (NHS) carbamates as ABPs for detection/profiling of trypsin-like serine proteases. This present study provides synthetic details for the preparation of additional examples of this ABP chemotype, which function as potent irreversible inhibitors of their respective target serine protease. We describe their use for the activity-based profiling of a broad range of serine proteases including trypsin, the trypsin-like protease plasmin, chymotrypsin, cathepsin G, and neutrophil elastase (NE), including the profiling of the latter protease in clinical samples obtained from patients with cystic fibrosis.

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.

Corticotroph isolation from Pomc-eGFP mice reveals sustained transcriptional dysregulation characterising a mouse model of glucocorticoid-induced suppression of the hypothalamus-pituitary-adrenal axis.

Glucocorticoids (GC) are prescribed for periods > 3 months to 1%-3% of the UK population; 10%-50% of these patients develop hypothalamus-pituitary-adrenal (HPA) axis suppression, which may last over 6 months and is associated with morbidity and mortality. Recovery of the pituitary and hypothalamus is necessary for recovery of adrenal function. We developed a mouse model of dexamethasone (DEX)-induced HPA axis dysfunction aiming to further explore recovery in the pituitary. Adult male wild-type C57BL6/J or Pomc-eGFP transgenic mice were randomly assigned to receive DEX (approximately 0.4 mg kg-1 bodyweight day-1 ) or vehicle via drinking water for 4 weeks following which treatment was withdrawn and tissues were harvested after another 0, 1, and 4 weeks. Corticotrophs were isolated from Pomc-eGFP pituitaries using fluorescence-activated cell sorting, and RNA extracted for RNA-sequencing. DEX treatment suppressed corticosterone production, which remained partially suppressed at least 1 week following DEX withdrawal. In the adrenal, Hsd3b2, Cyp11a1, and Mc2r mRNA levels were significantly reduced at time 0, with Mc2r and Cyp11a1 remaining reduced 1 week following DEX withdrawal. The corticotroph transcriptome was modified by DEX treatment, with some differences between groups persisting 4 weeks following withdrawal. No genes supressed by DEX exhibited ongoing attenuation 1 and 4 weeks following withdrawal, whereas only two genes were upregulated and remained so following withdrawal. A pattern of rebound at 1 and 4 weeks was observed in 14 genes that increased following suppression, and in six genes that were reduced by DEX and then increased. Chronic GC treatment may induce persistent changes in the pituitary that may influence future response to GC treatment or stress.

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.