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.

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

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

The use of cell free DNA (cfDNA) for mutational screening of multiple myeloma

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

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

Minimal Residual Disease (MRD) Testing in Newly Diagnosed Multiple myeloma (MM) Patients: A Prospective Head-to-Head Assessment of Cell-Based, Molecular, and Molecular-Imaging Modalities

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

DNA Sequencing of CD138 Cell Population Reveals TP53 and RAS-MAPK Mutations in Multiple Myeloma at Diagnosis

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

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

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

Simplified in-House Deep Sequencing Method of Inmunoglobulin Genes for Minimal Residual Dissease Quantification and Risk Stratification in Multiple Myeloma

Circulating Tumor DNA in the Peripheral Blood As Early Predictor of Clinical Outcome in Relapsed/ Refractory Multiple Myeloma

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

Novel Flow Cytometry-Based Biomarkers Predict Recurrence in Myeloma Patients through Detection of MRD in the Bone Marrow and CTCs in Peripheral Blood

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

RAS/RAF Mutations in Tumor Samples and Cell-Free DNA from Plasma and Bone Marrow Aspirates in Multiple Myeloma Patients

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

Retrospective Analysis of Cytogenetic and Clinical Characteristics on 65 Patients with Multiple Myeloma

Using Mass Spectrometry to Assess Disease Status in Treated Multiple Myeloma Patients