Abstract:Diagnostic bone marrow (BM) DNA and matched plasma-derived circulating tumour DNA (ctDNA) demonstrated excellent correlations in myelodysplastic syndromes (MDS) and acute myeloid leukaemia (AML),1-4 and dynamic ctDNA monitoring contributes to predict relapse in MDS and AML patients undergoing allogeneic hematopoietic stem cell transplantation (allo-HSCT).5 However, few researches are reported on the application of ctDNA dynamic monitoring assessments during the whole disease course in MDS and AML.1 Here, we assessed the feasibility and utility of ctDNA as a novel and minimally invasive biomarker based on targeted next-generation sequencing (NGS) to monitor treatment outcome, track clonal evolution and predict survival. To the best of our knowledge, we firstly demonstrated the application of ctDNA concentration as an effective prognostic biomarker in MDS and AML. The details of clinical data and methods are shown in Tables S1–S3. Thirty-five adult patients were involved in our study. Twenty-seven patients had pre- and post-treatment plasma-derived ctDNA, and paired baseline BM DNA NGS assessments were included for both concordance analysis and dynamic ctDNA analysis. Four patients with only baseline BM DNA and plasma-derived ctDNA NGS assessments were used for concordance analysis. Three patients with only pre- and post-treatment ctDNA assessments and one patient with only post-treatment ctDNA assessments were used for dynamic ctDNA analysis (Figure 1). A total of 46 mutated genes with 135 mutations were detected in BM DNA and plasma-derived ctDNA, involving 42 genes with 100 mutations (74.1%) detected both in BM DNA and plasma-derived ctDNA (Figure 2A, Figure S1). Among the 100 concordant mutations, the concordance of variant allele frequencies (VAFs) assessment based on plasma and BM was high (R = .854 p < .001) (Figure 2B). Comparing pre-treatment mean ctDNA concentrations with blasts in BM at baseline measured by cytological analysis, we found pre-treatment mean ctDNA concentrations were strongly associated with BM blasts at baseline (R = .618, p < .001) (Figure 2C), indicating mean ctDNA concentration were related to tumour burden. We also reached a moderate correlation between pre-treatment mean ctDNA VAF and BM blasts (R = .533, p = .001) (Figure S2). For patients with pre- and post-treatment paired samples, dynamic change of VAFs was analyzed. Patients were segregated into three groups, including those achieved a complete response (CR), those with CR with incomplete hematologic recovery (CRi) (for AML), morphologic CR (mCR) (for MDS)/morphologic leukaemia-free state (mLFS) (for AML) or stable disease (SD) (for MDS) and those with progressive disease (PD). Sixteen patients with CR had dramatically decreased VAFs (median VAF, 16.7% pre-treatment vs. 0% post-treatment, p < .001). Persistent VAFs were observed in CRi, mCR, mLFS and SD patients (median VAF, 3.6% pre-treatment vs. 3.3% post-treatment, p = .804). VAFs from PD patients tended to increase (median VAF, 4.7% pre-treatment vs. 6.7% post-treatment, p = .053) (Figure 2D). The post-treatment VAFs were lower in patients with CR versus those without CR (median VAF, 0% vs. 3.7%, p < .001) (Figure S3A). Similarly, delta VAFs from pre-treatment to post-treatment were greater in patients with CR versus those without CR (median delta VAF, 12.8% vs. −.3%, p < .001) (Figure S3B). As expected, post-treatment ctDNA positivity was associated with both shorter progression-free survival (PFS, median PFS, 5.6 vs. not reached [NR] months, p < .001) (Figure 3A) and overall survival (OS, median OS, 11.0 vs. NR months, p < .001) (Figure 3B). In the meantime, increased mean ctDNA VAF (defined as increasing post-treatment compared with pre-treatment ctDNA VAF) predicted decreased PFS (median PFS, 2.8 vs. 25.0 months, p < .001) (Figure 3C) and OS (median OS, 8.1 vs. 16.8 months, p = .014) (Figure 3D). Furthermore, we discovered patients with increased mean ctDNA concentration have poorer PFS (median PFS, 2.8 vs. NR months, p < .001) (Figure 3E) and OS (median OS, 7.9 vs. NR months, p < .001) (Figure 3F). However, the significant prognostic impact of pre-treatment ctDNA status on PFS and OS with mean ctDNA VAF (Figure S4A,B) or mean ctDNA concentration (Figure S4C and D) was not discovered, indicating pre-treatment ctDNA status was not a good biomarker to stratify the prognosis. To investigate the prognostic differences between mean ctDNA VAF and mean ctDNA concentration change from post- to pre-treatment, we compared the discrimination ability of the prognostic model with the concordance index (C-index). The model based on mean ctDNA concentration change showed better discrimination of PFS with C-index value of .79, compared to mean ctDNA VAF with C-index value of .75. The area under curve (AUC) for predicting 1-year PFS indicated more excellent conformity based on ctDNA concentration change, compared to mean ctDNA VAF (AUC, .84 vs. .82) (Figure 3G). Similarity was also discovered in OS analysis (C-index, .70 vs. .67; AUC, .78 vs. .72) (Figure 3H). Next, we tried to explore clonal evolution patterns in refractory or relapsed MDS and AML. We observed two patterns of progression from MDS to secondary AML (sAML). Patient 1 (P1) showed a linear pattern, where successive clones occurred within the previous parental clone. The acquisition of FLT3 mutation occurred during disease relapse and subsequently transformation to AML (Figure 4A). P7 and P8 showed a branched pattern (Figure 4B, Figure S5A, respectively). For P7, the original clone harboring PTPN11 mutation was suppressed while a new clone carrying FLT3-ITD mutation expanded, leading to AML transformation. However, we only discovered one linear pattern in refractory or relapsed AML. Three representative patients (P11, P26 and P35) achieved CR and relapsed within 6−7 months (Figure 4C, Figure S5B,C, respectively). The founding clone in the primary tumour in P11 contained mutations in WT1, GATA2 and RAD21 that were all recurrent at relapse. A new subclone occurred and evolved to become the dominant clone at relapse by acquiring additional SETD2 mutations. Compared with flow cytometry measurable residue disease results, we found that ctDNA molecular residue detection was earlier about 4 months. In summary, our study revealed that serial plasma-derived ctDNA assessments can reflect treatment response, survival, and clonal evolution in adult MDS and AML as a minimally invasive method, which warrant the prospective use of ctDNA as a biomarker in disease monitoring. This work was supported by the National Natural Science Foundation of China (grant number: 81800121) and Zhejiang Medical Association Clinical Research Fund project (grant number: 2017ZYC-A14). The authors also thank all the patients who participated in this study. The authors declare that they have no competing interests. The National Natural Science Foundation of China Grant Number: 81800121; Zhejiang Medical Association Clinical Research Fund project, Grant Number: 2017ZYC-A14. The data used in the present study are available from the corresponding author upon reasonable request. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

Abstract:Abstract Multiple Myeloma (MM) is the second most common hematologic malignancy, marked by uncontrolled clonal expansion of plasma cells. MM genome is complex and heterogeneous, with a high frequency of structural variants (SVs) and copy-number abnormalities (CNAs). Single-cell sequencing technologies offer advantages over traditional bulk methods in cancer genomics research for evaluating cellular heterogeneity and investigating evolution of cellular subpopulations within the tumor. Thus, there is an impetus to translate the growing understanding of the genomic landscape of MM into dissecting the tumor heterogeneity using scRNA-seq. So far, there are no reported studies systematically comparing tumor and immune populations differences between MM NON-progressors (NPs) and RAPID-progressor (RPs) and investigating how clonal plasma cells contribute to progression in a large cohort. In addition, the bone marrow microenvironment plays an important role in the evolution of premalignant MM and MM progression. A previous study of single-cell transcriptomics analysis of Monoclonal gammopathy of undetermined significance (MGUS) and MM tumor microenvironment (TME) revealed that natural killer (NK) cell abundance is frequently increased in the early stages and associated with altered chemokine receptor expression. This study shed the light on the role of immune cells on disease progression from asymptomatic MGUS to symptomatic MM. However to date, how immune cells influence disease progression within symptomatic MM is still unclear. Here, we subjected 344 CD138-negative Bone Marrow Mononuclear cells (BMMC) samples to scRNA-seq. To control for the first line therapy, 272 patients were treated with RVD and Autologous stem cell transplant (ASCT) and were stratified into 3 groups based on progression: 200 NON-progressors (NPs) (PFS> 5 years), 7 Intermediate Progressors (IPs, PFS 2-4 years) and 38 RAPID-progressors (RPs) (PFS<=18 months). In addition, there are 72 subjects treated with RVD but without ASCT, including 34 NPs, 2 IPs and 36 RPs. From MMRF CoMMpass study (NCT01454297), we also have whole exome sequencing (WES) and bulk RNA-seq from CD138-positive fraction of BMMC samples. Our preliminary analysis of scRNA-seq of 6 samples treated by the same first-line therapies (RVD but without ASCT) revealed that plasma cells (PCs) from two patients (MMRF2550 and MMRF2562) with chromosome 13q deletion clustered together, whereas PCs from a patient (MMRF2187) with t(8;14) Myc translocation formed a distinct cluster. Interestingly, PCs from a patient (MMRF2271) with both chromosome 13q deletion and t(8;14) formed another independent cluster. These observations highlight the important role of genetic drivers in transcriptome profiles of plasma cells. In addition, we further investigated how progression features correlate with the immune microenvironment and identified differentially expressed genes between NK cells from RPs versus those from NPs, which could potentially serve as MM progression markers. Interestingly, several cytotoxicity genes are significantly upregulated in FPs, such as GZMB and KLRB1 (log Fold Change=1.44, and 8.41 respectively). Overall, our preliminary results provide a small glimpse of the interconnected nature of driver genetic alterations and progression features in MM tumor and TME. This study will provide a sufficiently broad, deep, and diverse vast dataset for accurately characterizing MM at single-cell resolution to help interrogate how genetic alterations and disease progression interplay MM tumor and TME. We hope this study could identify novel candidate targets for therapeutic approaches, and ultimately stratify patients by risk of progression for early intervention in the clinic. Disclosures Kumar: Roche-Genentech: Consultancy, Research Funding; Merck: Research Funding; Takeda: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; KITE: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Amgen: Consultancy, Research Funding; Novartis: Research Funding; Abbvie: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Beigene: Consultancy; Carsgen: Research Funding; Bluebird Bio: Consultancy; Astra-Zeneca: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; Oncopeptides: Consultancy; Tenebio: Research Funding; BMS: Consultancy, Research Funding; Antengene: Consultancy, Honoraria; Adaptive: Membership on an entity's Board of Directors or advisory committees, Research Funding; Sanofi: Research Funding. Oh: Abbvie: Membership on an entity's Board of Directors or advisory committees; Blueprint Medicines: Membership on an entity's Board of Directors or advisory committees; Celgene Bristol Myers Squibb: Membership on an entity's Board of Directors or advisory committees; Constellation: Membership on an entity's Board of Directors or advisory committees; CTI Biopharma: Membership on an entity's Board of Directors or advisory committees; Disc Medicine: Membership on an entity's Board of Directors or advisory committees; Geron: Membership on an entity's Board of Directors or advisory committees; Incyte: Membership on an entity's Board of Directors or advisory committees; Kartos Therapeutics: Membership on an entity's Board of Directors or advisory committees; PharamaEssentia: Membership on an entity's Board of Directors or advisory committees; Sierra Oncology: Membership on an entity's Board of Directors or advisory committees. Vij: BMS: Research Funding; Takeda: Honoraria, Research Funding; Sanofi: Honoraria, Research Funding; BMS: Honoraria; GSK: Honoraria; Oncopeptides: Honoraria; Karyopharm: Honoraria; CareDx: Honoraria; Legend: Honoraria; Biegene: Honoraria; Adaptive: Honoraria; Harpoon: Honoraria. Avigan: Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; Pharmacyclics: Research Funding; Kite Pharma: Consultancy, Research Funding; Juno: Membership on an entity's Board of Directors or advisory committees; Partner Tx: Membership on an entity's Board of Directors or advisory committees; Karyopharm: Membership on an entity's Board of Directors or advisory committees; Bristol-Myers Squibb: Membership on an entity's Board of Directors or advisory committees; Aviv MedTech Ltd: Membership on an entity's Board of Directors or advisory committees; Takeda: Membership on an entity's Board of Directors or advisory committees; Legend Biotech: Membership on an entity's Board of Directors or advisory committees; Chugai: Membership on an entity's Board of Directors or advisory committees; Janssen: Consultancy; Parexcel: Consultancy; Takeda: Consultancy; Sanofi: Consultancy.

Cell-Free DNA from Blood or Bone Marrow as Alternatives to Bone Marrow Cells for Molecular Diagnostics and Monitoring of Multiple Myeloma

Serial Monitoring of Circulating Tumour DNA on Clinical Outcome in Myelodysplastic Syndromes and Acute Myeloid Leukaemia.

The Utility of Non‐invasive Liquid Biopsy for Mutational Analysis and Minimal Residual Disease Assessment in Extramedullary Multiple Myeloma

Evaluation of Circulating Tumour DNA for the Mutational Characterisation of Multiple Myeloma

Sensitivity, specificity, and accuracy of molecular profiling on circulating cell-free DNA in refractory or relapsed multiple myeloma patients, results of MM-EP1 study

Peripheral Blood Single-Cell Genotyping and Phenotyping in Multiple Myeloma Reveals Shared Mutations across Multiple Hematopoietic Cell Lineages

Monitoring treatment response and disease progression in myeloma with circulating cell-free DNA

Development of Circulating Tumor DNA (ctDNA) for Molecular Measurable Residual Disease (MRD) in Acute Myeloid Leukemia (AML)

Cell-Free DNA Analysis Demonstrates Comprehensive Genomic Profiling and a Facile Method to Sensitively in Multiple Myeloma

Sequencing of Circulating Cell-Free DNA in Patients with AML Detects Clinically Significant Mutations Not Detected in Bone Marrow: The Role for Complementary Peripheral Blood and Bone Marrow Genomic Analysis

Reliability of Liquid Biopsy and Next Generation Sequencing in Monitoring Residual Disease Post-Hematopoietic Stem Cell Transplant

Correlation between Peripheral Blood and Bone Marrow Somatic Mutations Among Patients with Suspected or Established Myelodysplastic Syndromes from the National MDS Study

Enrichment-Free Single-Cell Detection and Morphogenomic Profiling of Myeloma Patient Samples to Delineate Circulating Rare Plasma Cell Clones

Cell-free DNA measurable residual disease as a predictor of postallogeneic hematopoietic cell transplant outcomes.

High Levels of Circulating Tumor Cells Are Associated with Increased Bone Marrow Proliferation in Newly Diagnosed Multiple Myeloma Patients

Single-Cell RNA-Seq Analysis of CD138-Depleted Bone Marrow Samples Reveals Genetic Alterations and Disease Progression Correlate with Tumor and Bone Marrow Immune Microenvironment in the Mmrf Commpass Study

High level of circulating cell-free tumor DNA at diagnosis correlates with disease spreading and defines multiple myeloma patients with poor prognosis

Mutational and Copy Number Profiling of Circulating Tumor DNA in Acute Myeloid Leukemia Using Targeted Next Generation Sequencing

Single Cell Multi-Omic Correlation of Single Nucleotide Variants, Copy Number Variation and Surface Epitopes for Clonal Profiling of Myeloma

Clinical Significance of Circulating Cell-Free DNA Detection in Multiple Myeloma: A Meta-Analysis

Flow Cytometry Based Detection of MRD in Bone Marrow of Patients with Multiple Myeloma: A Comparison Between Fluorescent-Based Cytometry Versus Cytof