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 5263: Whole-genome sequencing as an alternative to analyze copy number abnormalities in acute myeloid leukemia and myelodysplastic syndrome

Whole-genome sequencing as an alternative to analyze copy number abnormalities in acute myeloid leukemia and myelodysplastic syndrome

Whole-genome sequencing as an alternative to analyze copy number abnormalities in AML, MDS and MM.

Whole Genome Sequencing As an Unbiased and Simplified Approach to Molecular Karyotyping in AML/MDS

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

Whole-Genome Sequencing for Copy Number Abnormalities in Multiple Myeloma Supersedes Karyotype Analysis and Fluorescent in Situ Hybridization

Abstract 4628: Validation of an automated, scalable comprehensive genomic profiling assay for hematologic malignancies

Abstract 5069: Retrospective analysis of acute myeloid leukemia using next generation cytogenomics identifies variants missed at diagnosis

Next-generation sequencing-based multigene mutational screening for acute myeloid leukemia using MiSeq: applicability for diagnostics and disease monitoring

Abstract 5891: Plasma-based BDA-NGS Test Holds Potential in Circulating Nucleic Acid Based Oncogenic Aberration Detection

Abstract 2628: Molecular Diagnosis for Pediatric Cancer Through Integrative Analysis of Whole-Genome, Whole-Exome and Transcriptome Sequencing

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

Abstract 5174: Utilization of Cell free messenger RNA Sequencing for non-invasive molecular profiling of multiple myeloma and Acute Myeloid Leukemia

Abstract A22: Complete mutation concordance with amplicon-based NGS between peripheral blood and bone marrow at 0.5% VAF threshold

Abstract 1763: Long-read sequencing of pediatric leukemia identifies clinically relevant genomic rearrangements

Abstract 2739: Low-pass whole genome sequencing detects copy number variations in circulating tumor DNA

Abstract 219: Single-tube NGS profiling allows identification of molecular signature in ALL patients

Abstract 7419: Seq2Karyotype (S2K): A method for deconvoluting heterogeneity of copy number alterations using single-sample whole-genome sequencing data

Targeted Next-Generation Sequencing of Circulating Cell-Free DNA Vs Bone Marrow in Patients with Acute Myeloid Leukemia.

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

Next-generation sequencing reveals gene mutations landscape and clonal evolution in patients with acute myeloid leukemia