Mulligan, George

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 Introduction: Multiple Myeloma (MM) is a complex hematologic malignancy, where patients often display one or more potentially high risk cytogenetic alterations associated with poor outcomes. Traditional high risk cytogenetics events in myeloma, such as NSD2 t(4;14) translocations, gain of 1q21, or deletion of 17p13, have often been studied in isolation, but emerging evidence suggests the combination of these high risk alterations has a greater negative prognostic impact than any on their own. This has been recently reflected in risk stratification guidelines from IMS and IMWG, which consider certain events, such as t(4;14), as high risk only in combination with gain of 1q21 or deletion of 1p32. Still, the mechanisms through which multi-hit high risk events disrupt the microenvironment and promote poor outcomes remains unclear. We hypothesize that myelomas with multiple high risk alterations remodel the immune environment to promote tumor progression and immune escape. Methods: To study the role of multi-hit high risk cytogenetics on the immune microenvironment, we analyzed CD138-bone marrow aspirate single-cell RNA sequencing (scRNAseq) data from 337 newly diagnosed MM patients from the MMRF CoMMpass study, generated as part of the Immune Atlas initiative. Cytogenetic alterations for these patients were derived from whole-genome and whole-exome sequencing of CD138+ bone marrow aspirates from the MMRF CoMMpass cohort. Patients were classified as multi-hit high-risk MM if they exhibited either (i) two high-risk cytogenetic events, including NSD2 t(4;14), MAF t(14;16), del(17p13), or (ii) one high-risk event with gain of 1q21. Results: Based on cytogenetic information, 33 patients in the Immune Atlas cohort were classified as 'multi-hit’ high-risk, while 34 patients had only a single high-risk alteration. These multi-hit high risk patents displayed significantly worse progression free survival compared to standard risk patients, adjusting for age, therapy, and ASCT (HR=2.34, p<7.2e-5), while single-hit patients showed a non-significant trend towards worse outcomes (HR=1.51, p=0.071). The most frequent multi-hit risk event in the cohort was NSD2 t(4;14) in combination with 1q21, observed in 23 patients. Compositional analysis of the immune populations in the scRNAseq data revealed that relative to both standard and single-hit high risk patients, multi-hit high risk patients displayed significant reductions in IFN-I stimulated CD4+ and CD8+ T cell populations, marked by decreased expression for antiviral response markers such as ISG15, XAF1, and IFI44L. Enrichment of inflammatory CD16+ non-classical monocytes was also observed. Further analysis of the most frequent multi-hit risk event, NSD2 with 1q21, revealed distinct enrichment of late-activated cytotoxic HLA+ CD8+T cells and cytotoxic CD4+.T cells. Neither alteration was observed in patients with single-hit NSD2 or 1q21. Further analysis on the interactions between high risk events revealed distinct transcriptomic profiles dependent on the high risk events involved, particularly in multi-hit events involving mutations or deletions of TP53 or gain of 1q21. Reduction in IFN-I signaling across various immune populations was primarily restricted to multi-hit risk events involving gain of 1q21. Conversely, patients with biallelic loss of TP53 displayed significant enrichment of these same IFN-I populations. Given this contrast in IFN-I signaling spanning multiple immune subtypes, we hypothesized that the malignant myeloma cells might similarly display altered IFN-I activity. Utilizing CD138+ myeloma bulk RNA sequencing data from CoMMpass, we found that IFN-I response scores strongly correlated between paired scRNA immune and bulk RNAseq tumor datasets (R2: 0.75, p=1.63e-32). Across 660 patients with bulk RNAseq and cytogenetic data, we confirmed significant reductions in IFN-I response scores in tumor cells in multi-hit high risk involving 1q21, and significant enrichment in those with complete loss of TP53 relative to standard-risk MM. Conclusions: Overall, we do observe that MM patients with multiple high risk cytogenetic alterations display distinct immunologic traits relative to standard risk or single-hit high risk patients in addition to significantly worse outcomes. Furthermore, we find variability in the immune alterations observed based on the cytogenetic events involved, further illustrating the heterogeneity in how MM shapes the immune microenvironment.

Abstract Background: The Multiple Myeloma Research Foundation (MMRF) and its partners launched the MyDRUG platform clinical trial (NCT03732703) to evaluate the safety and efficacy of genomically-guided treatments for functional high-risk relapsed/refractory myeloma patients (fHRMM). fHRMM patients with no discernible activating mutations and who had received 1-3 prior therapies (exposed to at least one proteasome inhibitor and an IMiD) were given a combination of daratumumab (D), ixazomib (I), pomalidomide (P) and dexamethasone (d) (D-IPd) quadruplet therapy. The goal of this study is to evaluate the immunomodulatory effects of daratumumab-(anti-CD38)-based therapy on fHRMM patients enrolled in the MMRF-sponsored MyDRUG umbrella clinical trial over the duration of therapy and to correlate immune changes with patient response immediately following the first two cycles of therapy. Methods Bone marrow (BM) aspirates were collected at baseline (BL) prior to therapy, at the end of 2 cycles (EOC2), 4 cycles (EOC4) of the quadruplet therapy and at the end of treatment (EOT) or progressive disease (PD). CD138-negative bone marrow mononuclear cells were isolated for profiling the immune microenvironment using scRNAseq, scTCRseq and CyTOF. Results Single-cell RNA and TCR profiles were generated from ~90,000 cells across 41 samples from 16 subjects. scRNAseq identified CD3+ T cells, NK cells, classical and non-classical monocytes and B cells. CD8+ T cells were further resolved into their subsets including precursor (TPEX) and terminally exhausted T cells (TEX) within the bone marrow microenvironment. Within the CD8+ T cell compartment, an expansion of the effector memory (TEM) cell population was observed following 2-cycles of therapy (EOC2). However, this TEM expansion was not sustained at timepoint EOC4. Towards the end of therapy, enrichment of precursors to putative exhausted and terminal dysfunctional effector memory cells (TEMRA) was observed. Subjects with partial (PR), very good partial (VGPR) and complete (CR) best overall responses showed an initial increase in TEM cells followed by a slow expansion of TEX and TEMRAcells over time. One patient with progressive disease (PD) maintained the highest LAG3 expression across the cohort among CD8 T cell subsets (TEMRA, EM) and CD4 Tregs at baseline while observing a significant upregulation of PDCD1 in CD4 Tregs and CD8 Naive cells by EOC2. Cell-cell communication network analysis identified secreted ligands between responders (CR/VGPR) and partial responders (PR) at EOC2. Samples with PR exhibited upregulation of signaling of TNF from the CD8 Naïve compartment (p=0.048, Wilcoxon), IFNG from the CD8 EM (p=0.073, Wilcoxon) and TEMRA (p=0.048, Wilcoxon) populations along with IFNG (p=0.073, Wilcoxon) from NK cells. CD14+ and CD14+16+ Monocytes on the other hand showed enrichment of LGALS9 (p=0.03 (CD14+), p=0.048 (CD14+CD16), Wilcoxon) in the PR group at EOC2 relative to responders. Finally, longitudinal tracking of T cell clonality showed preferential progressive expansion of TEMRA cells after the therapy. We mapped hyperexpanded (more than 16 copies) clonotypes, presumed to be tumor-reactive, to various CD8+ T cell subsets. These hyperexpanded clonotypes are mapped to the TEMRA, TEX and TEM CD8+ T cells. Analysis of TCR repertoires clonality and diversity also demonstrated a reduced richness and increased clonality of CD8+ TCRs in responders. Summary/Conclusion These preliminary analyses of D-IPd quadruplet therapy of RRMM indicate that paired single cell RNA and TCR sequencing of BM T cells suggests a selective expansion of clonotypic CD8+ T cells. When considering patient response, the only PD patient profiled in this study maintained an exhaustive signature (LAG3, PDCD1) in multiple T cell compartments at BL. The differential expression of the receptor ligand pairs after treatment between different clinical response groups suggests an altered TME response that may inform future therapeutic stratification. These results support further explorations of rational combinations of targeted and immune therapies for greater efficacy in this disease.

Abstract Background: The goal of this study is to evaluate the immunomodulatory effects of daratumumab-(anti-CD38)-based therapy on fHRMM patients enrolled in the MMRF sponsored MyDRUG umbrella clinical trial (NCT03732703) over the duration of therapy and to correlate immune changes with patient response. Methods: fHRMM who had received 1-3 prior therapies (exposed to at least one proteasome inhibitor and an IMiD) were given a combination of daratumumab (D), ixazomib (I), pomalidomide (P) and dexamethasone (d) (D-IPd) quadruplet therapy. Bone marrow (BM) aspirates were then collected from patients prior to therapy (BL), after 2 cycles (EOC2), after 4 cycles (EOC4) of the quadruplet therapy and at the end of treatment or disease progression (EOT). Immune and tumor fractions were enriched by CD138 selection. Here, CD138- bone marrow mononuclear cells were isolated and transcriptionally profiled by scRNAseq and scTCRseq. Results: Single-cell RNA and TCR profiles were generated from ∼90,000 cells across 41 samples from 16 subjects. scRNAseq identified CD3+ T cells, NK cells, classical and non-classical monocytes and B cells. CD8+ T cells were further resolved into their subsets including precursor (TPEX) and terminally exhausted T cells (TEX) within the bone marrow microenvironment. Within the CD8+ T cell compartment, an expansion of the effector memory (TEM) cell population was observed following 2-cycles of therapy (EOC2). However, this TEM expansion was not sustained at timepoint EOC4. Towards the end of therapy, enrichment of precursors to exhausted and terminal dysfunctional effector memory cells (TEMRA) was observed. Subjects with partial (PR), very good partial (VGPR) and complete (CR) best overall responses showed an initial increase in TEM cells followed by a slow expansion of TEX and TEMRA cells over time. Longitudinal tracking of T cell clonality showed preferential progressive expansion of TEMRA cells after the therapy. We mapped hyperexpanded (more than 16 copies) clonotypes, presumed to be tumor-reactive, to various CD8+ T cell subsets. These hyperexpanded clonotypes mapped to the TEMRA, TEX and TEM CD8+ T cells. Analysis of TCR repertoires clonality and diversity also demonstrated a reduced richness and increased clonality of CD8+ TCRs in responders. Conclusions: Paired single cell RNA and TCR sequencing of BM T cells from RRMM patients receiving D-IPd quadruplet therapy revealed a selective expansion of clonotypic CD8+ T cells with a clonal replacement of BM exhausted clonotypes by dysfunctional circulating CD8+ T cells in responders over time. Citation Format: Chaitanya R. Acharya, Darwin D'souza, Reyka Jayasinghe, Mark Hamilton, Zhihong Zhihong, Seunghee Kim-Schulze, Hearn J. Cho, Shaji Kumar, George Mulligan. The immunomodulatory effects of daratumumab-based therapy on functional high-risk relapse/refractory multiple myeloma (fHRMM) patients enrolled in the MyDRUG trial [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 5258.

Introduction As part of the Multiple Myeloma Research Foundation (MMRF) MyDRUG umbrella clinical trial, we compared immune cell transcriptomic and proteomic signatures of multiple myeloma bone marrow samples from patients assessed by single-cell RNA sequencing (scRNA-seq) and mass cytometry (CyTOF). The goal of our study is to evaluate the effects of cobimetinib on the immune repertoire over the treatment course at single-cell resolution. Methods Bone marrow samples were collected from patients prior to therapy with cobimetinib (baseline/BL), after 2 cycles of cobimetinib therapy (EOC2), after 2 cycles of combination cobimetinib with IPd (EOC4) and at the end of treatment or at disease progression (EOT). All samples were processed into CD138+ myeloma cell and CD138- “immune cell” fractions for single-cell analysis. In this study, the CD138- “immune cell” fractions were transcriptionally profiled by single-cell RNA-sequencing (scRNA) at up to four time points. Results Transcriptome profiles of approximately 50,000 cells were analyzed across 22 samples from 9 patients with tumor MAPK Pathway mutations including: BRAF (p.K483Q, p.V600E, n=2); KRAS (p.Q61H, p.G12V, p.G12R, n=3); NRAS (p.Q61H, p.Q61R, p.Q61K, n=4). A subset of paired samples were subject to CyTOF with two immune related panels to supplement the single-cell transcriptomic findings. Focusing on one immune panel, 4 samples were profiled at baseline and end of cycle 2 while two samples were profiled only at baseline capturing over 600,000 cells spanning both timepoints. Our immune focused CyTOF panel contained 41 markers including T and NK cell markers (CD4,CD8A,CD3, GZMB, CD56), Myeloid lineage markers (CD14), B/Plasma markers (CD19, CD20, SLAMF7, CD38) among others. The largest populations captured by CyTOF included 260,000 neutrophils, 139,000 monocytes, 130,000 CD4/CD8 T cells, 43,000 NK cells, 40,000 dendritic cells and 22,000 B/plasmablasts/myeloma cells. Exhaustion markers, including LAG3 and TIGIT, are increased in CD8 and CD4 T cells after MEK inhibition, while PD1 is downregulated. PD-L1 is upregulated in both T cell compartments (CD4/CD8), monocytes and in residual myeloma cells profiled by both CyTOF and scRNA-seq. Both the IFNG and STAT1 pathways are concurrently upregulated post MEK inhibition in transcriptionally distinct CD8 T cell subsets profiled by scRNA-seq, suggesting a mechanism influencing PD-L1 expression. Conclusions These data indicate IFNG pathway activation by CD8+ T cells post MEK inhibition promotes PD-L1 expression on residual myeloma cells, suggesting the PD-L1 and IFNG signaling axes is activated in relapsed refractory myeloma patients after MEKi treatment. Here we present immunophenotypic and transcriptional profiling of longitudinal samples pre/post-MEK inhibition with multiple myeloma patients which can inform precision guided treatment options.

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.