Avigan, David

The bone marrow immune microenvironment (BMME) is a critical factor in multiple myeloma (MM) pathogenesis. Despite higher incidence of MM among self-identified African Americans, the relationship between African genetic similarity and the BMME remains unclear. We present the largest analysis to date assessing the impact of Ancestry on the BMME in myeloma (N=320, including 64 with ≥50% African genetic similarity, AFR-High) leveraging our recently generated Immune Atlas. AFR-High was associated with enrichment of late T-effector cells, driven by rare cytotoxic CD4⁺ T cells, suggesting a distinct adaptive BMME. AFR-High exhibited expansion of late-differentiated adaptive NK cells with high cytotoxic GNLY and GZMH expression. AFR-High depicted reduced CD14⁺ monocytes, increased CD16⁺ non-classical monocytes, and expansion of macrophage-like monocyte populations with enrichment of metabolic adaptation, efferocytosis and phagocytic clearance programs, and immunoregulatory pathways. These findings reveal a shift from an innate inflammatory immune landscape toward adaptive cytotoxic and immunoregulatory remodeling associated with African genetic similarity in MM.

Abstract Although the human genome encodes a vast repertoire of noncoding RNAs that regulate gene expression, the noncoding genome remains underexplored due to technical challenges. Specifically, during transcriptomic sequencing data alignment, the overlap between noncoding and coding loci can create ambiguous read alignments that are subsequently discarded from downstream analysis. For this reason, most of the noncoding genome is excluded from standard genomic annotations used for sequencing alignment. To address this challenge and enable concurrent profiling of the coding and noncoding transcriptome, we systematically integrated standard coding (GENCODE) and noncoding (LncBook) genome annotations, preserving coding gene annotations and removing overlapping noncoding regions. The resulting integrated genome annotation expanded the number of annotated noncoding genes from 40,785 to 138,296 while preserving all coding genes and reducing ambiguous read assignment. To evaluate the utility of our integrated genome annotation for uncovering novel, biologically relevant noncoding RNAs (ncRNAs), we realigned CD138-positive bulk RNA-seq (N = 942) and CD138-negative single-cell RNA-seq (N = 478) data from the MMRF CoMMpass study, generating a comprehensive coding-noncoding atlas of the myeloma bone marrow microenvironment with noncoding genes representing 51% of highly variable genes and displaying significant cell type specificity. Tumor expression profiling based on this integrated profiling identified 15 clusters, including two enriched for amp(1q21) or t(4;14) and associated with shorter progression-free survival (PFS). Differential expression and systematic filtering yielded 19 candidate high-risk ncRNAs, including previously uncharacterized ENSG00000310209 , which was associated with poor PFS (HR = 1.141, P = 0.0025), increased IRF4 activity, Wnt pathway activation, CCL5 signaling, and the accumulation of anergic-like CD8+ T cells. These findings establish integrated coding-noncoding analysis as a strategic approach for discovering functional ncRNAs from transcriptomic sequencing data.

Abstract The CoMMpass study (NCT01454297) is a prospective, longitudinal observational study involving 1,141 newly diagnosed multiple myeloma (NDMM) patients. Bone marrow aspirates from the study participants were collected for comprehensive molecular characterization of tumor (using bulk RNA sequencing and whole genome sequencing) and tumor immune microenvironment using 3’ single-cell RNA sequencing (scRNA-seq). To support exploration of this rich multi-omic dataset, we developed CoMMpass Explorer (CE), an intuitive interactive platform that enables real-time analysis of clinical and genomic data from CoMMpass. A central feature of CE is cohort building, where users can filter and stratify patients by clinical, genomic, and survival data elements to create custom cohorts for analysis. CE provides five main views: Overall Summary for clinical feature distribution, Kaplan-Meier curves, and multivariate Cox proportional hazards models; Mutational Profile for visualizing and comparing somatic mutations; Tumor Profile for bulk differential expression and gene set enrichment; Immune Microenvironment for comparing cell type abundance and cell cycle dynamics between cohorts using scRNA-seq data; and Pseudo-Bulk, which aggregates single-cell expression by cell type to enable cohort-level comparisons of transcriptional programs using DEseq2. CE reproduces prior CoMMpass studies, including WEE1 expression being associated with poorer progression-free survival (Simhal, et al.) and different mutation frequencies between African American and European American patients (Manojlovic, et al.). We further demonstrate CE’s utility by interrogating the 2025 Consensus Genomic Staging (CGS) risk classification system (Avet-Loiseau, et al). We compared CGS high-risk (HR, n=249) vs. CGS standard-risk (SR, n=573) patients. HR patients have significantly worse progression-free survival, overall survival, and time to second-line therapy compared with SR patients. Univariate models recapitulated expected differences in the clinical features that define CGS risk. Bulk tumor RNA-seq and ssGSEA highlighted a more proliferative and genomically unstable transcriptional program in HR group, with upregulation of angiogenesis and cell-cycle-related pathways. In the immune microenvironment, monocytes were more abundant in the SR group. Pseudo-bulk RNA-seq of plasma cells showed overexpression of ANXA1, NSD2, IGF1R, MAF and related genes in HR. CE democratizes access to the CoMMpass multi-omic resource, allowing researchers to interrogate clinical and molecular heterogeneity. By enabling hypothesis generation and replication of published results, CE serves as a valuable tool for identifying prognostic biomarkers and meaningful patterns in gene expression, mutation profiles, and the immune microenvironment. Citation Format: Wanxing Zhang, Chaitanya R. Acharya, Steven M. Foltz, David E. Avigan, Samir Parekh, Ravi Vij, Shaji Kunnathu Kumar, Taxiarchis Kourelis, Sagar Lonial, Hearn Cho, Ioannis S. Vlachos, Sacha Gnjatic, Li Ding, Manoj Bhasin, George Mulligan. CoMMpass Explorer: An interactive platform to explore clinical and genomic data from newly diagnosed multiple myeloma patients from the landmark CoMMpass observational study [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 5504.

Abstract Background: Multiple myeloma is the second most common adult blood cancer and remains incurable, despite recent advances. The IMWG/IMS recently published updated risk criteria, Consensus Genomic Staging (CGS), combining genomic events with clinical measures to identify high-risk patients (Avet-Loiseau, et al. JCO 2025). We applied CGS risk criteria to baseline samples from the Multiple Myeloma Research Foundation CoMMpass Study (NCT01454297), a prospective, longitudinal observational study that enrolled more than1100 participants between 2011-2017. Methods: We assessed patient risk following CGS criteria. Criteria 1: del(17p), with a cutoff of >20% clonal fraction, and/or TP53 mutation. Criteria 2: an IgH translocation including t(4;14), t(14;16), or t(14;20) along with 1q+ and/or del(1p32). Criteria 3: monoallelic del(1p32) along with 1q+ or biallelic del(1p32). Criteria 4: beta-2 microglobulin >= 5.5 mg/L with normal creatinine. Patients meeting one or more criteria were classified as High Risk. To test if risk showed association with genetic ancestry, we grouped patients into those with high African genetic ancestry and those with low African genetic ancestry. Genetic ancestry was measured using ADMIXTURE in comparison to 1000 Genomes. Results: We evaluated the CGS risk status of 877 CoMMpass patients, finding 69.2% (573/828) to be Standard Risk (SR) and 30.8% (255/828) to be High Risk (HR). HR patients showed worse progression-free survival (median PFS 26.3 months vs. 42.5 months; log rank test, p < 0.0001) and overall survival (median OS 63.9 months vs. median not reached; p < 0.0001). 10.6% (91/858) of patients met Criteria 1, 10.7% (93/873) met Criteria 2, 5.6% (49/877) met Criteria 3, and 8.7% (73/835) met Criteria 4. We found 14.2% (123/867) of patients to have high African-like genetic ancestry (AFR-high) and 85.8% (744/867) to have low African-like genetic ancestry (AFR-low). We observed significantly fewer HR patients within the AFR-high group (20%, 23/115) than in the AFR-low group (32.4%, 231/712) (Fisher’s exact test, p = 0.0065). The AFR-high group met HR criteria at a lower rate across all criteria, but only with Criteria 3 was the difference significant (p-value = 0.03). Higher levels of creatinine in AFR-high patients may have impacted meeting the normal creatinine requirement of Criteria 4. Conclusions: HR patients comprise 30.8% of CoMMpass subjects and show significantly worse PFS and OS than SR patients. AFR-high patients met HR criteria at a lower rate than AFR-low patients according to the CGS criteria. Clinical trials that utilize CGS criteria to determine patient eligibility may be challenged to enroll a representative patient population, including AFR-high patients. Screening and eligibility criteria should be carefully considered to mitigate potential enrollment disparities. Citation Format: Steven M. Foltz, Chaitanya R. Acharya, Alexander Gout, Yuxin Jin, David E. Avigan, Ravi Vij, Shaji Kunnathu Kumar, Sagar Lonial, Hearn Cho, Katie Wozniak, Jonathan J. Keats, Craig Cole, John D. Carpten, George J. Mulligan. Multiple myeloma risk classification of CoMMpass using Consensus Genomic Staging identifies differences between patients with high and low African genetic ancestry [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 3887.

Abstract Multiple myeloma is a malignancy of clonally expanded plasma cells shaped by complex interactions with the immune microenvironment (IME). To investigate immune correlates of treatment response and disease progression, we conducted multi-omics profiling including CD138neg single-cell RNA sequencing of 243 bone marrow samples from 102 patients (631,226 cells) and CD138pos bulk RNA and whole-genome sequencing from 209 samples. In longitudinal analyses, interferon-γ signaling associated with markers of impaired T-cell memory after autologous stem cell transplant, whereas naïve B-cell abundance and immunoglobulin diversity correlated with improved progression-free survival (HR = 0.48; P = 2.3e−4). At disease progression, multiple myeloma cells upregulated cancer–testis antigens (CTAg) and immune effector genes, with concurrent B-cell depletion, enrichment of myeloid-derived suppressor cell expression, and phenotypic T-cell exhaustion. These findings highlight dynamic immune–tumor interactions, identifying naïve B-cell reconstitution as a biomarker of durable response and CTAgs as potential targets for high-risk disease at progression. Significance: Longitudinal profiling of multiple myeloma and the IME revealed dynamic immune–tumor interactions across the disease course. Altered expression in CD8+ T cells limited memory phenotype after transplant, whereas naïve B-cell recovery associated with sustained treatment response. At progression, CTAg expression associated with immunosuppression, revealing novel mechanisms of immune dysregulation.

Abstract Multiple myeloma (MM) is an incurable malignancy of clonally expanded plasma cells shaped by complex interactions with the immune microenvironment. To investigate immune factors driving treatment response and resistance, we conducted multi-omics profiling including CD138 neg single-cell RNA sequencing of 243 bone marrow samples from 102 patients (631,226 cells) and CD138 pos bulk RNA and whole-genome sequencing from 209 samples. Longitudinal analyses revealed that interferon gamma signaling impairs T cell memory after autologous stem cell transplant, while naïve B cell abundance and immunoglobulin diversity correlated with improved progression-free survival (HR = 0.48, p = 2.3e -4 ). At disease progression, MM cells upregulated cancer-testis antigens and immune effector genes, with concurrent B cell depletion, enrichment of myeloid-derived suppressor cell genes in monocytes, and T cell exhaustion. These findings highlight dynamic immune-tumor interactions, identifying naïve B cell reconstitution as a biomarker of durable response, and cancer-testis antigens as potential targets for high-risk disease at progression. Statement of Significance Longitudinal profiling of multiple myeloma and the immune microenvironment revealed dynamic immune-tumor interactions across the disease course. Dysfunctional CD8⁺ T cells limited memory formation post-transplant, while naïve B recovery associated with sustained treatment response. At progression, cancer-testis antigen expression associated with immunosuppression, revealing novel mechanisms of immune escape.

Introduction: In recent years, multiple myeloma (MM) research has increasingly focused on unraveling the complex interactions between malignant plasma cells and their surrounding bone marrow microenvironment (BMME). To further elucidate these interactions, we developed a comprehensive single-cell atlas of the MM BMME, integrating both coding and noncoding (ncRNA) transcriptional profiles. By correlating these profiles with specific cytogenetic abnormalities in myeloma cells, we aim to uncover how ncRNAs regulate the phenotypic and functional states of myeloma cells and, consequently, shape the cellular landscape of the BMME. A detailed understanding of how the noncoding transcriptome affects the cellular components and immune landscape of the BMME is essential for advancing our knowledge of MM pathogenesis and identifying novel biomarkers and therapeutic targets. Methods: Our atlas was constructed using single-cell RNA sequencing data from 481 samples of CD138neg cells sorted from bone marrow aspirates of MM patients enrolled in the MMRF CoMMpass study. This dataset includes accompanying whole-genome sequencing and survival data, enabling correlative analysis between transcriptional profiles, cytogenetic abnormalities, and patient outcomes. To generate a combined coding and non-coding transcriptome, we developed an expanded human reference genome by systematically merging mRNA and ncRNA transcripts from the LncBook2.0 and GENECODEv42 reference genomes. The sequencing data was aligned to this expanded genome, followed by quality control processing, batch correction, clustering, and supervised analysis. Results: The resulting coding and non-coding atlas comprises over 1.9 million cells, spanning immune, plasma, and stromal compartments. Remarkably, ncRNAs constituted 46% of the associated genes, including 22,102 long ncRNAs (lncRNAs) and 3,112 small ncRNAs. Subclustering analysis of the major immune compartments (B lymphoid, myeloid, and NK and T lymphoid) yielded 70 subclusters, with ncRNAs representing 24-38% of the top 100 differentially expressed genes within each compartment. This underscores the significant role of ncRNAs in distinguishing immune subpopulations. Moreover, examining the association between specific cytogenetic abnormalities and the immune composition of the BMME, 30 subclusters were found to be differentially abundant (P < 0.05) between cytogenetic abnormalities, supporting our hypothesis that myeloma cell genotypes significantly shape the immune microenvironment. For instance, patients with CCND1 amplification displayed an increase CD4+ central memory T cells with a decrease in CD8+ cytotoxic T cells. Similarly, patients with 1q21 gain exhibited an increase in TGFβ-stimulated monocytes, while a decrease in cytotoxic CD8+ T and CD56dim NK cells was observed in patients with 17p13 deletion. Next, examining which ncRNAs are associated with specific cytogenetic abnormalities, we identified 14 ncRNAs significantly differentially expressed across cytogenetic abnormalities. Increased expression of six of these ncRNAs were associated with poor overall survival both independently and in combination (P < 0.05, HR = 3.7). Notably, three of these ncRNAs have documented roles in cancer progression but have not been previously studied in the context of MM. Subsequently stratifying patients based on their enrichment of these outcome-associated ncRNAs revealed even stronger associations with immune composition, including an increase in CD4+ regulatory T cells (P < 0.05). This further supports our hypothesis that myeloma cell genotypes shape the BMME. Lastly, by constructing a gene regulatory network to predict the interactions between coding and noncoding RNAs, we illustrate the regulatory networks of ncRNAs modulating the phenotypes of myeloma and immune cell subpopulations. Conclusion: In summary, we present the first high-resolution transcriptomic atlas of the MM BMME, integrating both coding and noncoding RNAs. This study reveals distinct immune subpopulations associated with specific myeloma cell genotypes and highlights the central role of ncRNAs in modulating the phenotypic and functional states of both malignant and non-malignant cells in the MM BMME. Our findings emphasize the potential of ncRNAs as therapeutic targets and prognostic markers for MM, offering new avenues for research and clinical intervention.

Introduction Newly-diagnosed multiple myeloma (MM) patients often share tumor genetic abnormalities associated with a high risk of progression and poor survival. We hypothesized that single-cell transcriptomic characterization of the tumor immune microenvironment (TIME) may be combined with tumor-based risk stratification to uncover novel disease biology mechanisms associated with clinical outcomes and improve risk stratification. Methods We performed single-cell RNA sequencing of 361 CD138 negative sorted bone marrow mononuclear cells (BMMCs) samples from 263 MM patients enrolled at time of diagnosis in the Multiple Myeloma Research Foundation CoMMpass study (NCT01454297). Specimens were analyzed at baseline (n=263) across a network of five collaborating academic research centers following assay harmonization. Genetic factors that determine risk were defined as one or more of the following gain- or loss-of-function events derived from baseline genomic sequencing data: del17p14, translocations of the Ig locus with the MAF/A/B or WHSC1/MMSET/NSD2 loci, and gain of chromosome 1q21. Statistical modeling strategies were used to assess single-cell populations or RNA signatures associations with known risk factors or disease progression (FDR<0.05). Results The study generated high-quality single-cell profiles of 731,643 cells from 123 high-risk and 107 standard-risk MM patients based on genetics. Based on distinct transcriptome profiles, BMMCs formed 106 clusters corresponding to 5 major compartments: NK and T cells, B cells, erythrocytes, myeloid cells, and plasma cells. Our preliminary differential abundance analysis comparing high- vs standard-risk patients at baseline uncovered significant differences in the TIME associated with genetic risk. The high-risk group was identified as enriched in effector CD8+, and CD4+ T cells. In contrast, the high-risk group exhibited depletion of naive CD8+ and interferon-associated CD4+ and CD8+ T as compared to the standard-risk group. Within subclusters in the T and NK cells compartment, a higher cytotoxicity signature was observed in high-risk patients, characterized by the expression of granzymes and perforin. Finally, our preliminary survival and progression analysis associated a higher frequency of IDO1-expressing macrophages in worst progression patients, while better outcomes were associated with CD8+ T and B cells expressing androgen response and IL6/STAT3/EMT programs, respectively. Conclusions Current approaches for risk stratification are based on bulk measurements of tumor's genetic or epigenetics features, which do not capture the diversity and dynamics of the TIME. Through single-cell transcriptomics, we were able to capture, in high granularity, the TIME and correlate specific immune cell populations and phenotypes with relapse risk and poor prognostic outcomes. These results suggest that immune subpopulations may be an essential novel aspect for improving current risk stratification models.

Introduction: BCMA-directed CAR-T and Bispecific T-cell engager (BiTE) therapies have shown remarkable overall response rates in heavily pretreated patients with relapsed/refractory multiple myeloma (RRMM). This has led to the approval of idecabtagene vicleucel, ciltacabtagene autoleucel and teclistamab by the FDA, EMA and other regulatory agencies. We performed a global survey on accessibility, practice patterns, and toxicity management of commercially available CAR-T and BiTE therapies. Methods: An online survey was completed by clinicians attending the International Myeloma Society's Immune Effector Cell Therapies in Multiple Myeloma Workshop, in-person or virtually, made available for 30 days after the meeting held in March 2023, Boston, MA. Results: 87 clinicians from 16 countries completed the survey, most practicing in an academic practice (86%) and in an urban environment (89%). Sixty-four percent of the respondents practice in North America followed by Europe (25.3%) including Germany, United Kingdom, Ireland, Italy, Netherlands, Spain, France, Greece, and Cyprus, and 6.9% from Asia and 3.4% from Australia. Seventy-seven percent reported CAR-T cell therapy was approved in their country with 44.8% treating 1-5 patients and 14.9% treating >5 patients per month; 67% reported BiTE therapy was approved in their country with 14.9% treating >5 patients per month. Of evaluable responses (n = 40), 62% administered CAR-T therapy inpatient and 15% administered CAR-T therapy outpatient ≥ 50% of the time. Length of time between decision to treat with CAR-T therapy to infusion was 6-12 weeks for 42% of responders and <6 weeks for only 20%. Most (75%) reported inpatient administration of BiTE therapy step-up dosing; 17.5% reported outpatient step-up dosing administered ≥ 50% of the time. Use of institutional or practice SOPs for CAR-T and BiTE therapy patient selection was reported by 85% and 58%, respectively, with the factors influencing SOP outlined in Table 1. Of available CAR-T therapies, 53% prefer cilta-cel and 23% reported CAR-T therapy preference depended most on patients' co-morbidities/frailty (89%), timing of next available apheresis (67%), age (56%), CNS disease (44%) and other factors. Maintenance therapy after CAR-T therapy was reportedly used by 25% and 5% of respondents in 1-10% and 11-50% of their patients, respectively, to “maintain level of response achieved” (50%) or “achieve a deeper level of response” (25%). About half (43%) of the respondents reported using an SOP or particular criteria for secondary Hemophagocytic Lympho-Histiocytosis/Macrophage Activation Syndrome/CAR-related-HLH (sHLH/MAS/CarHLH) management that varied between the use of ASTCT Working Group IEC associated HLH-like syndrome Criteria (47%) (Hines 2023), CARTOX HLH Criteria (24%) (Neelapu 2017) and HLH-2004 criteria (6%). Tocilizumab was used as cytokine release syndrome (CRS) prophylaxis only by one responder. Keppra neurologic prophylaxis was used by 47.5% and 15% of responders for CAR-T and BiTE therapy, respectively. Delayed neurotoxicity with CAR-T was reported in 1-5% of patients by 57.5% of responders and 6-10% of patients by 15% of responders; most common manifestations were other movement/neuromuscular disorder (62%), nerve palsy (55%), and Parkinsonism (45%). Infectious prophylaxis was used to protect against HSV (92.5%), PJP (90%), bacterial pathogens (55%), fungal pathogens (60%), and IVIG (77.5%). However, criteria for infectious prophylaxis use varies widely incorporating IgG levels, infectious history, CD4+ count, B-cell aplasia and cytopenias (48.1 - 88.9%). Most respondents (46.3 - 86.5%) reported completing an infectious work-up including CMV, EBV, Hep B/C, and HIV ( Table 2). Conclusions: Trends of CAR-T and BiTE therapy use and patterns of concurrent prophylactic and supportive care in MM vary widely among academic centers across the globe, demonstrating a need for further understanding of toxicity pathogenesis, toxicity management and consistent guidelines. This survey also demonstrates the great unmet need for access beyond urban centers among the international myeloma community. Prospective studies and claims data analysis should be utilized to evaluate the longitudinal cost of care, quality of life, duration of response and management of infectious disease, neurological or related complications of CAR-T and bispecific TCE therapies.

Abstract BACKGROUND Standard of care (SOC) and patient survival in glioblastoma have changed little in the past 17 years. We evaluated in a phase 3 trial whether adding an autologous tumor lysate-loaded dendritic cell vaccine (murcidencel) to SOC extends survival. Patients and METHODS Newly diagnosed glioblastoma patients were randomized 2:1 to either murcidencel or placebo. Under a crossover design, all patients could receive murcidencel following tumor recurrence. All parties remained blinded regarding treatments before recurrence. Patients thus received murcidencel at new diagnosis (nGBM) or at recurrence (rGBM) following crossover from placebo. The primary and secondary endpoints compare overall survival (OS) with contemporaneous, matched external controls. Four sets of analyses were conducted to ensure rigorous matching of the controls, reduce biases, and confirm the robustness of the results. RESULTS 331 patients were enrolled. With the crossover, 89% received murcidencel. Median OS (mOS) for nGBM patients (n = 232) was 19.3 months from randomization (22.4 months from surgery) with murcidencel vs. 16.5 months from randomization in the controls (HR = 0.80, p = 0.002). Survival at 48 months from randomization was 15.7% vs. 9.9%, and at 60 months was 13% vs. 5.7%. For rGBM (n = 64), mOS was 13.2 months from relapse vs. 7.8 months in the controls (HR = 0.58, p < 0.001). Survival at 24 months post-recurrence was 20.7% vs. 9.6%, and at 30 months post-recurrence was 11.1% vs 5.1%. In nGBM patients with methylated MGMT (n = 90), mOS was 30.2 months from randomization (33 months from surgery) with murcidencel vs. 21.3 months from randomization in the controls (HR = 0.74, p = 0.027). The treatment was well tolerated, with only 5 serious adverse events deemed at least possibly related to the vaccine. CONCLUSION Clinically meaningful and statistically significant survival extension was seen in both nGBM and rGBM patients treated with murcidencel and SOC compared with contemporaneous, matched external controls who received SOC alone.