Undetectable circulating tumor DNA (ctDNA) levels correlate with favorable outcome in metastatic melanoma patients treated with anti-PD1 therapy
Journal of Translational Medicine volume 17, Article number: 303 (2019)
Treatment with anti-PD1 monoclonal antibodies improves the survival of metastatic melanoma patients but only a subgroup of patients benefits from durable disease control. Predictive biomarkers for durable benefit could improve the clinical management of patients.
Plasma samples were collected from patients receiving anti-PD1 therapy for ctDNA quantitative assessment of BRAFV600 and NRASQ61/G12/G13 mutations.
After a median follow-up of 84 weeks 457 samples from 85 patients were analyzed. Patients with undetectable ctDNA at baseline had a better PFS (Hazard ratio (HR) = 0.47, median 26 weeks versus 9 weeks, p = 0.01) and OS (HR = 0.37, median not reached versus 21.3 weeks, p = 0.005) than patients with detectable ctDNA. Additionally, the HR for death was lower after the ctDNA level became undetectable during follow-up (adjusted HR: 0.16 (95% CI 0.07–0.36), p-value < 0.001). ctDNA levels > 500 copies/ml at baseline or week 3 were associated with poor clinical outcome. Patients progressive exclusively in the central nervous system (CNS) had undetectable ctDNA at baseline and at subsequent assessments. In multivariate analysis adjusted for LDH, CRP, ECOG and number of metastatic sites, the ctDNA remained significant for PFS and OS. A positive correlation was observed between ctDNA levels and total metabolic tumor volume (TMTV), number of metastatic sites and total tumor burden.
Assessment of ctDNA baseline and during therapy was predictive for tumor response and clinical outcome in metastatic melanoma patients and reflected the tumor burden. ctDNA evaluation provided reliable complementary information during anti-PD1 antibody therapy.
The treatment of metastatic melanoma has changed significantly during the last decade. Approved medical treatments available are targeted therapy (BRAF and MEK inhibitors) and immunotherapy (immune checkpoint inhibitors—ICI) . Only patients with BRAFV600—mutant disease (45–50%) can benefit from both targeted and immunotherapy, while all metastatic melanoma patients can receive immunotherapy because predictive biomarkers to select patients for immunotherapy are not available yet in the clinical setting. The first approved immunotherapy showing improved overall survival (OS) for metastatic melanoma patients was ipilimumab, an anti-CTLA-4 blocking monoclonal antibody . The anti-PD-1 monoclonal antibodies, nivolumab and pembrolizumab are the newer generation of ICI providing a higher response rate, longer progression free (PFS) and OS with a lower frequency of adverse events [3, 4]. The ICI enhance the immune system response to melanoma by blocking inhibitory receptors on T lymphocytes. Even though ICI have a significant impact on survival, half of the patients do not respond and will progress on therapy within 1 year. Currently no predictive biomarkers are available to identify patients who benefit most from immunotherapy. Across different tumor types, PD-L1 staining is used to predict response; however melanoma patients with PD-L1 negative tumors can achieve durable objective response [5,6,7]. Therefore, PD-L1 expression is not used as a predictive biomarker. In this context one technique amply investigated is the analysis of the tumor circulating DNA (ctDNA) present in the plasma of cancer patients . ctDNA has the advantage to be minimally invasive compared to invasive tumor biopsies. Additionally serial blood samples can be easily assessed for successive analyses during therapy to document and potentially predict the course of the disease. Our group and others showed that longitudinal quantitative measurement of BRAFV600mut in ctDNA is a useful tool to monitor responses to targeted therapy [9,10,11,12]. Furthermore we showed that monitoring during anti-PD1 therapy was informative in a series of seven representative patients with clinical benefit from anti-PD1 treatment . In a prospective trial Lee et al.  showed that undetectable ctDNA levels at baseline or within the first 8 weeks is an independent predictor biomarker for response and better OS in a population of patients receiving pembrolizumab or nivolumab ± ipilimumab. Furthermore Cabel et al.  showed the same trend in other types of tumors: non-small cell lung cancer, uveal melanoma and microsatellite-instable colorectal cancer. In our prospective translational research study we evaluated metastatic melanoma patients receiving pembrolizumab in monotherapy over a follow up period of almost 2 years with longitudinal BRAFV600mut and NRASQ61/G12/G13mut ctDNA monitoring during treatment.
Materials and methods
The study population was identified among 141 metastatic melanoma patients who participated in an institutional therapeutically non-interventional clinical trial protocol on the outcome of immune-checkpoint inhibition for patients with advanced melanoma at the Universitair Ziekenhuis Brussel (UZ Brussel, Brussels Belgium). Patients were eligible for this sub-study on ctNDA if they had been treated with pembrolizumab and at the condition that BRAF V600 or NRAS Q61/G12/G13 mutations were detected in tumor tissue and/or plasma samples. Outcome data were collected prospectively between September 2014 and April 2017 in patients who had initiated treatment with the anti-PD1 antibody therapy (pembrolizumab 2 mg/kg 3 weeks) in the early access program prior to reimbursement of the drug in Belgium. Blood samples were prospectively collected after obtaining informed consent with an Ethical Committee (UZ Brussel) approved document. Plasma samples were collected together with routine blood collections (every 2 or 3 weeks during the entire period of immunotherapy).
ctDNA extraction and PCR analysis by ddPCR and Idylla platform
Silica-based extraction of DNA and subsequent allele-specific quantitative PCR (qPCR) to detect the BRAF wild-type gene and the G1798 > A and T1799 > A changes in the BRAF gene were performed with Idylla™ (Biocartis) on 1 ml of the stored plasma. The G1798 > A change is present in patients with V600K, V600R, and V600M mutations, whereas the T1799 > A change is present in patients with V600E, V600K, V600E2, and V600D mutations. A linear correlation between Cq values reported by prototype Idylla assays and digital droplet PCR was previously established, allowing precise and sensitive quantification of mutant ctDNA fragments in plasma down to 3 mutant copies per PCR reaction with an analytical sensitivity of 0.01%.
For the detection of NRASQ61/G12/G13 mutations DNA was extracted using the QIAamp Circulating Nucleic Acid Kit (Qiagen), followed by droplet digital PCR (ddPCR, Bio-rad). ddPCR was performed using the Bio-Rad QX-200 system (Biorad, Hercules, USA). Assays were purchased from Bio-Rad at 20× concentration (see list below).
ddPCR reaction mixtures contained a final concentration of 250 nM for each of the probes, 450 nM for the forward and reverse primers, 1× ddPCR™ Supermix for Probes (No dUTP) (Bio-Rad #186-3024 USA) and 8.8 µl of DNA in a final volume of 22 µl. Twenty microliter of this ddPCR reaction volume were loaded in appropriate wells of a DG8 cartridge (Bio-Rad #186-4008, USA) with 70 µl of generator oil (Bio-Rad #186-3030, USA) into the oil well. Samples are partitioned into approximately 20,000 water–oil emulsion droplets, each 0.85 nl in volume, using the QX200™ Droplet generator™ (Bio-Rad). Forty microliter of this water–oil emulsion were used for the ddPCR assay by transferring it into a 96-wells plate sealed with a PX1™ PCR plate Sealer (Bio-Rad, USA). ddPCR were performed with a T100™ thermal cycler (Bio-Rad, USA) under the following conditions: 1 cycle of 95 °C for 10 min, 40 cycles of 94 °C for 30 s and 55 °C for 1 min, and 1 cycle of 98 °C for 10 min. Cycled droplets were read individually with the QX200TM droplet reader (Bio-Rad). No template control wells (NTC) were included in the assays. Data was analyzed using QuantaSoft TM software version 1.6.6.0320.
Clinical outcome and statistical analysis
The ctDNA was dichotomized as detectable and undetectable to discriminate between responders and non-responders. Tumor responses were evaluated using the immune-related response criteria in solid tumors (irRC) at 2- or 3-monthly intervals. Patients (n = 7) who had stopped treatment and/or had died before the first response evaluation (all due to clinical progressive disease) were categorized as progressive disease (PD). Patients were grouped in three categories based on their best objective tumor response (BOR): objective response [complete remission (CR) and partial response (PR), CR/PR group], stable disease (SD), and progressive disease (PD). The CR/PR were considered responders, while SD/PD were considered non-responders. PFS and OS were described via the Kaplan–Meier method from the start of therapy to the date of first reported PD or death for PFS and death for OS. Probability of PFS or OS were compared between subgroup using the Log-rank test. In order to determine the adjusted relative risk of death after the undetectable ctDNA occurrence during the follow-up, we developed an extended multivariable Cox model in the overall population and in the defined subgroups. Undetectable ctDNA during the follow-up was introduced in the model as a time-dependent variable using the time to the first occurrence. The nonparametric Mann–Whitney test, one-tailed was used to compare the medians for the ctDNA values between the groups of responders versus non-responders. For additional comparisons of the means between the two groups the unpaired t test with Welch’s correction was used.
In order to identify factors associated with the occurrence of OS and PFS, univariate and multivariable regression analyses were performed using cox regression model. ctDNA, LDH, CRP, number of tumor sites and ECOG variables were considered in the analyses. The values of ctDNA at any time were considered as time-dependent variable. The results are presented as crude and adjusted hazard ratios (HR) with 95% confidence intervals (95% CI). Total metabolic tumor volume (TMTV) was assessed by 18-fluorodeoxyglucose positron emission tomography (18FDG-PET/CT) using MIM Encore SoftwareVR. TMTV was defined as the sum of all tumor-associated voxels with a standardized uptake value (SUV) higher than the mean SUV measured in a reference region in normal liver tissue + 3 standard deviations. Baseline disease burden was determined by the sum of the product of bi-dimensional diameters (SPOD) for every metastasis. Based on the median the patients were divided in three categories: low, intermediate and high burden disease.
Statistical analyses were carried out using Graphpad Prism 7 and IBM® SPSS® statistics software v24. All reported p values are two-sided. The p values lower than 0.05 were considered statistically significant.
Patient characteristics and ctDNA analysis
The baseline clinical characteristics of the 85 metastatic melanoma patients with a BRAFV600E/K or NRASQ61/G12/G13 mutation who were included in the study are summarized in Table 1. The median age was 57 years (range 27–82). The majority of patients had an ECOG of 0 (56 patients [65.9%]) and had stage IV-M1c (AJCC staging 7th edition) at baseline (67 patients [78.8%]). Brain metastases were present in 31 patients (36.5%) and LDH was higher than the upper limit of normal in 37 patients (43.5%).
Tumor BRAFV600 mutations were present in 63 patients (57 BRAFV600E and 6 BRAFV600K) and 22 patients had an NRASQ61/G12/G13 mutations (Fig. 1a). These mutations were quantitatively assessed in plasma samples collected every 3 weeks at the occasion of pembrolizumab administration. A total of 495 plasma samples from these 85 patients were analyzed (327 samples for BRAFV600 mutations and 168 samples for NRASQ61/G12/G13 mutations); 63 samples were obtained at baseline (in 35 samples the level of ctDNA was undetectable; 22 samples tested positive for BRAFV600 mutations and 6 samples for NRASQ61/G12/G13 mutations). Distribution of ctDNA BRAFV600 and NRASQ61/G12/G13 copy numbers (medians and means) for baseline, week 3, 6 and 9 are illustrated in Fig. 1b. In total 457 samples were obtained during treatment (with a median of five samples per patient [range 1 to 19 samples per patient], and a median of five samples for BRAFV600 mutations and seven samples for NRASQ61/G12/G13 mutations).
Correlation between baseline ctDNA analysis and clinical outcome
After a median follow-up of 84 weeks (95% CI 63.2–104.8), 60 (70.6%) patients had progressed and 45 (56.2%) had died (2 deaths were not melanoma related). The median PFS and OS of the total study population (n = 85 patients) was 13 weeks (95% CI 9.5–16.5) and 56 weeks (95% CI 7.2–104.8), respectively. The median PFS and OS for the subpopulation of 63 patients who had evaluable plasma samples for baseline ctDNA analysis was 9 weeks (95% CI 4.8–13.1) and 56 weeks (95% CI 7–104.9), respectively. The median OS in the subgroup with undetectable baseline ctDNA (n = 35) was not reached, while in the group of patients with detectable baseline ctDNA the median OS was 21 weeks (95% CI 0–43; p = 0.003) (Fig. 2a). The median PFS in the group with baseline undetectable ctDNA was 26 weeks (95% CI 0–71.1) compared to 9 weeks (95% CI 6.9–11; p = 0.008) for the group with detectable baseline ctDNA (Fig. 2b). The 1-year survival rates were respectively 70% (95% CI 54.3–85.7) in the subgroup with undetectable base line ctDNA levels versus 32% (95% CI 14.36–49.64) in the subgroup with detectable baseline ctDNA levels. Furthermore the 2.5-years survival was 54% (95% CI 34.4–73.6) versus 16% (95% CI − 7.52 to 39.52). Comparing the ctDNA levels in the subgroup of responders (CR/PR) versus non-responders (SD/PD), the median pretreatment ctDNA level in CR/PR-patients was 0 copies/ml of plasma (IQR 0–69) versus 31 copies/ml of plasma (IQR 0–647) for SD/PD patients (Additional file 1: Figure S1A) with a statistically significant difference between the two groups (p = 0.0345, Mann–Whitney test). The baseline ctDNA absolute levels for the two groups (63 patients) are shown in Additional file 1: Figure S1B. The difference for the mean value of BRAF and NRAS mutant copy number between the two group of patients was statistically significant (p = 0.0052, unpaired t test with Welch’s correction). The mean copy number of BRAF and NRAS mutant in the responders group (CR/PR) was 42.69 copies (± 19.09) compared to 1199 copies (± 392.7) in the non-responders group. In the CR/PR group the majority of patients had undetectable ctDNA (11 patients—69%), while 5 patients (31%) had detectable ctDNA ranging from 25 to 233 copies/ml of plasma (Fig. 3). In the SD/PD group 22 patients (49%) had undetectable baseline ctDNA and 23 patients (51%) had detectable ctDNA ranging from 31 to 11.641 copies/ml of plasma.
ctDNA monitoring during the anti-PD1 therapy
In an extended multivariable Cox model, undetectable ctDNA during the follow-up was introduced in the model as a time-dependent variable using the time to the first occurrence. The HR for death was lower after the ctDNA level became undetectable during follow-up (adjusted HR: 0.16 (95% CI 0.07–0.36), p-value < 0.001). The difference in OS was statistically significant for the group of patients who had at least one undetectable ctDNA value during follow up compared with the patients who had ctDNA detectable in the plasma at each time point (64.2% alive in the first group versus 15.6% alive in the second group, p < 0.001) (Fig. 4). Moreover in the subpopulation of 75 patients (excluding patients with exclusively brain PD) 76.7% of those who had at least one undetectable ctDNA value during follow up were alive compared to 15.6% alive in the group that had detectable ctDNA in the plasma at each time point (p < 0.001) (Fig. 4). A comparative analysis was performed for BRAF/NRAS ctDNA levels at second, third and fourth pembrolizumab administration. The median ctDNA at later time points (week 3, 6 and 9) remained 0 (IQR 0–0) in the responders group compared to detectable ctDNA median value in the non-responders group. The difference between the two groups was statistically significant at these three time points (Additional file 2: Figure S2A, D, G). The significant difference in terms of PFS and OS between the patients with undetectable ctDNA and detectable ctDNA was observed as well at these subsequent three time points: at week 3 median PFS was 89 weeks (95% CI 20.5–157.4) versus 6 weeks (95% CI 4.9–7.1) (p < 0.001) and median OS not reached versus 25 weeks (95% CI 9.2–40.7) (p < 0.001); at week 6 median PFS was estimated at 86 weeks (95% CI not available) versus 3 weeks (95% CI 2.3–3.6) (p < 0.001) and median OS not reached versus 42 weeks (95% CI 15.2–68.7) (p = 0.002); at week 9 median PFS was estimated at 83 weeks (95% CI not available; 75th percentiles 9) versus 0 (95% CI not available; 75th percentiles 0, at the median all patients had PD) weeks (p = 0.001) and median OS not reached (75th percentiles 81 weeks) versus 26 weeks (95% CI 0.7–51.2) (p = 0.002). The Kaplan–Meier curves are shown in Additional file 2: Figure S2B, C, E, F, H, I.
ctDNA levels in the subgroup of patients progressive exclusively in the CNS
Within the subgroup of “PD patients”, we identified 10 patients who presented with PD exclusively in the central nervous system (CNS). The majority of the patients (9/10) had CNS involvement when starting anti-PD-1 therapy and one patient developed brain metastases during treatment (MEL40). These patients had undetectable ctDNA at baseline and at subsequent assessments (36 samples analyzed in total). In this subgroup 9 patients had BRAFV600E mutations and 1 patient had a BRAFV600K mutation. The median number of samples tested per patient was 3 (ranging from 1 to 8 samples per patient) with a median of 7 pembrolizumab administrations (ranging from 1 to 16). All samples showed undetectable ctDNA even at later time points after week 9 (Additional file 3: Table S1). The median PFS in this subgroup was 3 weeks (95% CI 0–6.8) and the median OS was 33 weeks (95% CI 0–85.7) (Additional file 4: Figure S3A). An interesting finding was observed for patient MEL36 that showed PD with multiple brain metastases under targeted therapy (dabrafenib 150 mg bid/trametinib 2 mg qd) and received pembrolizumab as second line treatment. This patient was presenting a PR in the visceral metastases under targeted therapy with no metabolic activity on the PET-CT. The analysis of the plasma ctDNA at baseline showed undetectable levels despite the PD in the brain. Early on during pembrolizumab treatment the ctDNA levels became detectable in parallel with PD observed as well in the visceral metastases in addition to the brain PD (Additional file 4: Figure S3B).
High baseline or on-treatment levels of BRAF/NRAS mutations in ctDNA correlate with poor clinical outcome
In the PD group 11 out of 50 patients (22%) had baseline ctDNA levels higher than 500 copies/ml. For this subgroup the median PFS was 5 weeks (95% CI 3–6.9) and median OS was 8 weeks (95% CI 4.8–11.1) (Additional file 4: Figure S3C). Eight additional patients in the PD group had more than 500 copies/ml at week 3 assessment (no baseline available for 6 patients and 2 patients had increased ctDNA levels from baseline to week 3: 31 to 1294 copies/ml and 177 to 986 copies/ml) (Additional file 3: Table S2). For this additional subgroup of 8 patients, the median PFS was 6 weeks (95% CI 3.2–8.7) and the median OS was 11 weeks (95% CI 5.5–16.5) (Additional file 4: Figure S3D). The median value for ctDNA at baseline in the group of 11 patients was 3 652 (IQR 1241–8322) copies/ml and for the group of 8 patients was 1 973 (IQR 1140–27,192) copies/ml. In these two groups the majority of the patients (17/19; 89%) had baseline LDH higher than the upper limit of normal, and baseline CRP was higher than the upper limit of normal in 14/19 (74%) patients. Additionally, the majority [14/18 (78%); 1 was NE] had more than three metastatic sites. Moreover, an important increase in the ctDNA levels was observed (from 104.5 copies at baseline to 1980 copies/ml only 4 days after the first pembrolizumab administration) in one patient (MEL 67) who died 3 weeks later (data not shown). This patient had normal baseline LDH, elevated CRP (> 10xN) and more than three metastatic sites.
The multivariable Cox regression analysis including LDH, CRP, number of metastatic sites (> 3) and ECOG variables showed that undetectable ctDNA at baseline remained significantly correlated with PFS and OS (Table 2). Additionally, the ECOG status and baseline CRP remained significant for PFS, while for OS only the number of metastatic sites remained an additional significant co-variate. This analysis was performed in the subpopulation of 75 patients (excluding CNS exclusively PD patients).
Correlation between baseline ctDNA levels and disease burden
Baseline ctDNA levels correlated with total metabolic tumor volume (TMTV) assessed by 18FDG-PET/CT (Fig. 5a). We found a positive correlation between ctDNA levels and TMTV (Pearson correlation r = 0.321, p = 0.04). The TMTV evaluation was available in 40 patients. Additionally the total tumor burden (volume of the metastatic disease) and the number of metastatic sites were assessed. This assessment was available for 62 patients. A positive correlation was observed between the number of metastatic sites and ctDNA levels (Spearman correlation, r = 0.5374, p < 0.0001) (Fig. 5b). The patients were divided in three groups after tumor burden assessment: low, intermediate and high tumor burden. The comparison between the group showed a significant difference between the low and high tumor burden group (p = 0.0022; median 206.5 [95% CI 144–3346] versus 0 [95% CI 0–0]) and between the high and intermediate tumor burden group (p = 0.048; median 206.5 [95% CI 144–3346] versus 0 [95% CI 0–233]) (Fig. 5c).
PD-1 immune checkpoint inhibitors have become a standard of care for patients diagnosed with metastatic melanoma. However, since the majority of patients remains in need of alternative treatment options there is a high need for the identification of reliable novel biomarkers that can predict response and probability for long-term PFS, detect early disease progression, and assist in timely clinical decision making. The recent research is focused on the characterization of tumor immune microenvironment together with the analysis of mutational tumor burden. Current findings showed that a T cell inflamed gene expression and the mutational burden of the tumor are independent prognostic factors for the response to anti-PD1 therapy across tumor types . Additionally the MHC I expression predicted resistance to anti-CTLA4 therapy, while MHC class II membrane expression and interferon-γ (IFN-γ) gene signatures predicted response to anti-PD-1 . Moreover although the PD-L1 expression does not predict response in all patients, it can give an information regarding a higher probability of response but being only one factor among many that affect how patients respond [18, 19]. Another additional promising candidate biomarkers is the use of the real-time ctDNA analysis . A growing number of findings are available in different type of cancers regarding the ability of ctDNA to be informative for response to therapies, relapse detection or resistance mechanism uncovering [21, 22]. The first ctDNA liquid biopsy approved for use in clinical settings was in lung cancer patients for the identification of EGFR mutations for first-line therapy or identifying resistance mutations that will allow for treatment with third generation EGFR inhibitors [23, 24]. In melanoma there is no liquid biopsy test approved in the clinical settings at present. Nevertheless several studies showed the utility of ctDNA as a diagnostic, predictive and prognostic biomarker [9,10,11, 13, 14]. Thus the landscape for ctDNA utilization in clinical setting may change in the near future for melanoma patients especially for diagnostic purposes in advanced metastatic emergency situations if tumor biopsy is invasive and risky or the tissue results are not informative or will be available after a significant delay. Moreover, ctDNA can also ease the decision in the daily clinical practice when radiological evaluation is problematic especially for patients receiving immunotherapy. In this context, an important advantage of ctDNA is the possibility of non-invasive serial testing for monitoring treatment response and resistance to therapy. Regarding the predictive value for ctDNA in melanoma there is one study showing that ctDNA level is an independent predictive biomarker for response, PFS and OS in melanoma patients receiving anti-PD1 antibody-based therapy . Our results confirm these observations. In their study the patients received pembrolizumab or nivolumab ± ipilimumab, while our study population was more homogeneous all patients receiving pembrolizumab monotherapy. In our population similar findings were observed: undetectable baseline or early decline to undetectable ctDNA levels (during the first 6–9 weeks) correlated with prolonged PFS and OS. This pattern was observed in the CR/PR group. An additional study showed ctDNA correlation with response and PFS but again in a heterogeneous population where patients received either targeted therapy or immunotherapy . In line with previous observations in melanoma patients [14, 25], the subgroup of patients among the PD patients that had exclusive PD in the CNS presented undetectable ctDNA levels at any time point during treatment. This was also observed in breast and lung cancer patients with restricted CNS disease thus it seems to be a general observation [26,27,28,29]. One hypothesis is that the blood–brain barrier plays a major role for this event because ctDNA can be easily detected in cerebrospinal fluid (CSF) of patients with brain lesions. Moreover, it can be the modality of choice to confirm metastatic disease in patients with leptomeninges carcinomatosis or brain lesions with limited accessibility or difficult to sample lesions due to the high invasiveness of the biopsy procedure.
Additionally, high ctDNA levels (> 500 copies/ml) at baseline or week 3 observed only in the group of PD patients were indicative of a very poor clinical outcome. Moreover, in the PD group the patients with ctDNA levels < 500 copies/ml at baseline never presented a decline to undetectable ctDNA during the first 9 weeks. A cut-off for OS lower or higher than 2 years for patients on targeted therapy was shown to be 216 mutated copies/ml in the manuscript of Sanmamed et al. . Finding a precise cut-off informative for response to therapy is still a challenge at this point in time. We hypothesize that high ctDNA corresponds mainly to high proliferative activity/status of tumor cells (aggressive tumor phenotype with poor clinical outcome), which in the majority of the patients correlates as well with high tumor burden.
An added value for our approach is that using the cartridge system on Idylla platform for mutation analysis, the results can be available the same day as plasma collection rendering possible a clinical decision taking also into account ctDNA level. This can be extremely important when the interpretation of the radiological evaluation is problematic for patients receiving immunotherapy. With the development of cartridges testing multiple genes and implicitly multiple mutations, this technique can be used to identify, besides the dynamics of the disease, the mutational spectrum from disease sites. However, advancements are needed to overcome the difficulties related to undetectable ctDNA levels when progression occurs in cerebral sites. This should be considered if the technique will be introduced as routine clinical practice. This observation reinforces the importance of brain MRI or CT use during follow up of metastatic melanoma patients. Nevertheless, for the majority of the patients ctDNA monitoring was informative for tumor response or disease progression.
In clinics, ctDNA provides reliable complementary information to imaging during anti-PD-1 antibody therapy improving the clinical follow-up and evaluation of the efficacy of anti-PD-1 therapy. ctDNA assessment was predictive for tumor response and clinical outcome in metastatic melanoma patients except for patients with disease progression exclusively in the CNS. An added value for our approach is that using an automatized cartridge system, the results can be available the same day as plasma collection rendering possible in the future a clinical decision taking into account the ctDNA level.
Availability of data and materials
All raw data are available for review if necessary.
anti-programmed cell death 1
- BRAF :
B-Raf proto-oncogene, serine/threonine kinase
- NRAS :
NRAS proto-oncogene, GTPase
circulating tumor DNA
Droplet Digital polymerase chain reaction
progression free survival
(adjusted) hazard ratio
central nervous system
Eastern Cooperative Oncology Group
mitogen-activated protein kinase kinase 1
immune checkpoint inhibition
cytotoxic T-lymphocyte associated protein 4
programmed cell death 1 ligand 1
quantitative polymerase chain reaction
Expanded Access Program
American Joint Committee on Cancer
best overall response
positron emission tomography/computed tomography
- EGFR :
epidermal growth factor receptor
magnetic resonance imaging
Luke JJ, Flaherty KT, Ribas A, Long GV. Targeted agents and immunotherapies: optimizing outcomes in melanoma. Nat Rev Clin Oncol. 2017;14:463–82. https://doi.org/10.1038/nrclinonc.2017.43.
Hodi FS, O’Day SJ, McDermott DF, Weber RW, Sosman JA, Haanen JB, Gonzalez R, Robert C, Schadendorf D, Hassel JC, Akerley W, van den Eertwegh AJM, Lutzky J, Lorigan P, Vaubel JM, Linette GP, Hogg D, Ottensmeier CH, Lebbé C, Peschel C, Quirt I, Clark JI, Wolchok JD, Weber JS, Tian J, Yellin MJ, Nichol GM, Hoos A, Urba WJ. Improved survival with ipilimumab in patients with metastatic melanoma. N Engl J Med. 2010;363:711–23. https://doi.org/10.1056/NEJMoa1003466.
Robert C, Long GV, Brady B, Dutriaux C, Maio M, Mortier L, Hassel JC, Rutkowski P, McNeil C, Kalinka-Warzocha E, Savage KJ, Hernberg MM, Lebbé C, Charles J, Mihalcioiu C, Chiarion-Sileni V, Mauch C, Cognetti F, Arance A, Schmidt H, Schadendorf D, Gogas H, Lundgren-Eriksson L, Horak C, Sharkey B, Waxman IM, Atkinson V, Ascierto PA. Nivolumab in previously untreated melanoma without BRAF mutation. N Engl J Med. 2015;372:320–30. https://doi.org/10.1056/NEJMoa1412082.
Robert C, Schachter J, Long GV, Arance A, Grob JJ, Mortier L, Daud A, Carlino MS, McNeil C, Lotem M, Larkin J, Lorigan P, Neyns B, Blank CU, Hamid O, Mateus C, Shapira-Frommer R, Kosh M, Zhou H, Ibrahim N, Ebbinghaus S, Ribas A, KEYNOTE-006 investigators. Pembrolizumab versus ipilimumab in advanced melanoma. N Engl J Med. 2015;372:2521–32. https://doi.org/10.1056/NEJMoa1503093.
Daud AI, Wolchok JD, Robert C, Hwu W-J, Weber JS, Ribas A, Hodi FS, Joshua AM, Kefford R, Hersey P, Joseph R, Gangadhar TC, Dronca R, Patnaik A, Zarour H, Roach C, Toland G, Lunceford JK, Li XN, Emancipator K, Dolled-Filhart M, Kang SP, Ebbinghaus S, Hamid O. Programmed death-ligand 1 expression and response to the anti-programmed death 1 antibody pembrolizumab in melanoma. J Clin Oncol. 2016;34:4102–9. https://doi.org/10.1200/JCO.2016.67.2477.
Larkin J, Chiarion-Sileni V, Gonzalez R, Grob JJ, Cowey CL, Lao CD, Schadendorf D, Dummer R, Smylie M, Rutkowski P, Ferrucci PF, Hill A, Wagstaff J, Carlino MS, Haanen JB, Maio M, Marquez-Rodas I, McArthur GA, Ascierto PA, Long GV, Callahan MK, Postow MA, Grossmann K, Sznol M, Dreno B, Bastholt L, Yang A, Rollin LM, Horak C, Hodi FS, Wolchok JD. Combined nivolumab and ipilimumab or monotherapy in untreated melanoma. N Engl J Med. 2015. https://doi.org/10.1056/nejmoa1504030.
Morrison C, Pabla S, Conroy JM, Nesline MK, Glenn ST, Dressman D, Papanicolau-Sengos A, Burgher B, Andreas J, Giamo V, Qin M, Wang Y, Lenzo FL, Omilian A, Bshara W, Zibelman M, Ghatalia P, Dragnev K, Shirai K, Madden KG, Tafe LJ, Shah N, Kasuganti D, de la Cruz-Merino L, Araujo I, Saenger Y, Bogardus M, Villalona-Calero M, Diaz Z, Day R, Eisenberg M, Anderson SM, Puzanov I, Galluzzi L, Gardner M, Ernstoff MS. Predicting response to checkpoint inhibitors in melanoma beyond PD-L1 and mutational burden. J Immunother Cancer. 2018;6:32. https://doi.org/10.1186/s40425-018-0344-8.
Bettegowda C, Sausen M, Leary RJ, Kinde I, Wang Y, Agrawal N, Bartlett BR, Wang H, Luber B, Alani RM, Antonarakis ES, Azad NS, Bardelli A, Brem H, Cameron JL, Lee CC, Fecher LA, Gallia GL, Gibbs P, Le D, Giuntoli RL, Goggins M, Hogarty MD, Holdhoff M, Hong SM, Jiao Y, Juhl HH, Kim JJ, Siravegna G, Laheru DA, Lauricella C, Lim M, Lipson EJ, Marie SKN, Netto GJ, Oliner KS, Olivi A, Olsson L, Riggins GJ, Sartore-Bianchi A, Schmidt K, Shih IM, Oba-Shinjo SM, Siena S, Theodorescu D, Tie J, Harkins TT, Veronese S, Wang T-L, Weingart JD, Wolfgang CL, Wood LD, Xing D, Hruban RH, Wu J, Allen PJ, Schmidt CM, Choti MA, Velculescu VE, Kinzler KW, Vogelstein B, Papadopoulos N, Diaz LA. Detection of circulating tumor DNA in early- and late-stage human malignancies. Sci Transl Med. 2014;6:224ra24. https://doi.org/10.1126/scitranslmed.3007094.
Schreuer M, Meersseman G, van Den Herrewegen S, Jansen Y, Seremet T, Bott A, Chevolet I, Wilgenhof S, Maertens G, Neyns B. Applications for quantitative measurement of BRAF V600 mutant cell-free tumor DNA in the plasma of patients with metastatic melanoma. Melanoma Res. 2016;26:157–63. https://doi.org/10.1097/CMR.0000000000000224.
Schreuer M, Meersseman G, Van Den Herrewegen S, Jansen Y, Chevolet I, Bott A, Wilgenhof S, Seremet T, Jacobs B, Buyl R, Maertens G, Neyns B. Quantitative assessment of BRAF V600 mutant circulating cell-free tumor DNA as a tool for therapeutic monitoring in metastatic melanoma patients treated with BRAF/MEK inhibitors. J Transl Med. 2016;14:95. https://doi.org/10.1186/s12967-016-0852-6.
Haselmann V, Gebhardt C, Brechtel I, Duda A, Czerwinski C, Sucker A, Holland-Letz T, Utikal J, Schadendorf D, Neumaier M. Liquid profiling of circulating tumor DNA in plasma of melanoma patients for companion diagnostics and monitoring of BRAF inhibitor therapy. Clin Chem. 2018;64:830–42. https://doi.org/10.1373/clinchem.2017.281543.
Sanmamed MF, Fernández-Landázuri S, Rodríguez C, Zárate R, Lozano MD, Zubiri L, Perez-Gracia JL, Martín-Algarra S, González A. Quantitative cell-free circulating BRAFV600E mutation analysis by use of droplet digital PCR in the follow-up of patients with melanoma being treated with BRAF inhibitors. Clin Chem. 2015;61:297–304. https://doi.org/10.1373/clinchem.2014.230235.
Seremet T, Planken S, Schreuer M, Jansen Y, Delaunoy M, El Housni H, Lienard D, Del Marmol V, Heimann P, Neyns B. Illustrative cases for monitoring by quantitative analysis of BRAF/NRAS ctDNA mutations in liquid biopsies of metastatic melanoma patients who gained clinical benefits from anti-PD1 antibody therapy. Melanoma Res. 2018;28:65–70. https://doi.org/10.1097/CMR.0000000000000415.
Lee JH, Long GV, Boyd S, Lo S, Menzies AM, Tembe V, Guminski A, Jakrot V, Scolyer RA, Mann GJ, Kefford RF, Carlino MS, Rizos H. Circulating tumour DNA predicts response to anti-PD1 antibodies in metastatic melanoma. Ann Oncol. 2017;28:1130–6. https://doi.org/10.1093/annonc/mdx026.
Cabel L, Riva F, Servois V, Livartowski A, Daniel C, Rampanou A, Lantz O, Romano E, Milder M, Buecher B, Piperno-Neumann S, Bernard V, Baulande S, Bieche I, Pierga JY, Proudhon C, Bidard FC. Circulating tumor DNA changes for early monitoring of anti-PD1 immunotherapy: a proof-of-concept study. Ann Oncol. 2017. https://doi.org/10.1093/annonc/mdx212.
Cristescu R, Mogg R, Ayers M, Albright A, Murphy E, Yearley J, Sher X, Liu XQ, Lu H, Nebozhyn M, Zhang C, Lunceford JK, Joe A, Cheng J, Webber AL, Ibrahim N, Plimack ER, Ott PA, Seiwert TY, Ribas A, McClanahan TK, Tomassini JE, Loboda A, Kaufman D. Pan-tumor genomic biomarkers for PD-1 checkpoint blockade-based immunotherapy. Science. 2018. https://doi.org/10.1126/science.aar3593.
Rodig SJ, Gusenleitner D, Jackson DG, Gjini E, Giobbie-Hurder A, Jin C, Chang H, Lovitch SB, Horak C, Weber JS, Weirather JL, Wolchok JD, Postow MA, Pavlick AC, Chesney J, Hodi FS. MHC proteins confer differential sensitivity to CTLA-4 and PD-1 blockade in untreated metastatic melanoma. Sci Transl Med. 2018. https://doi.org/10.1126/scitranslmed.aar3342.
Ayers M, Lunceford J, Nebozhyn M, Murphy E, Loboda A, Kaufman DR, Albright A, Cheng JD, Kang SP, Shankaran V, Piha-Paul SA, Yearley J, Seiwert TY, Ribas A, McClanahan TK. IFN-γ-related mRNA profile predicts clinical response to PD-1 blockade. J Clin Invest. 2017;127:2930–40. https://doi.org/10.1172/JCI91190.
Ayers M, Lunceford J, Nebozhyn M, Murphy E, Loboda A, Albright A, Cheng J, Kang SP, Ebbinghaus S, Yearley J, Shankaran V, Seiwert T, Ribas A, McClanahan T. Relationship between immune gene signatures and clinical response to PD-1 blockade with pembrolizumab (MK-3475) in patients with advanced solid tumors. J Immunother Cancer. 2015;3:P80. https://doi.org/10.1186/2051-1426-3-S2-P80.
Husain H, Velculescu VE. Cancer DNA in the circulation: the liquid biopsy. JAMA. 2017;318:1272–4. https://doi.org/10.1001/jama.2017.12131.
Heitzer E, Perakis S, Geigl JB, Speicher MR. The potential of liquid biopsies for the early detection of cancer. NPJ Precis Oncol. 2017;1:36. https://doi.org/10.1038/s41698-017-0039-5.
Hench IB, Hench J, Tolnay M. Liquid biopsy in clinical management of breast, lung, and colorectal cancer. Front Med. 2018;5:9. https://doi.org/10.3389/fmed.2018.00009.
Bernabé R, Hickson N, Wallace A, Blackhall FH. What do we need to make circulating tumour DNA (ctDNA) a routine diagnostic test in lung cancer? Eur J Cancer. 2017;81:66–73. https://doi.org/10.1016/j.ejca.2017.04.022.
Bracht JWP, Mayo-de-Las-Casas C, Berenguer J, Karachaliou N, Rosell R. The present and future of liquid biopsies in non-small cell lung cancer: combining four biosources for diagnosis, prognosis, prediction, and disease monitoring. Curr Oncol Rep. 2018;20:70. https://doi.org/10.1007/s11912-018-0720-z.
Wong SQ, Raleigh JM, Callahan J, Vergara IA, Ftouni S, Hatzimihalis A, Colebatch AJ, Li J, Semple T, Doig K, Mintoff C, Sinha D, Yeh P, Silva MJ, Alsop K, Thorne H, Bowtell DD, Gyorki DE, Arnau GM, Cullinane C, Kee D, Brady B, Kelleher F, Dawson MA, Papenfuss AT, Shackleton M, Hicks RJ, McArthur GA, Sandhu S, Dawson S-J. Circulating tumor DNA analysis and functional imaging provide complementary approaches for comprehensive disease monitoring in metastatic melanoma. JCO Precis Oncol. 2017. https://doi.org/10.1200/po.16.00009.
De Mattos-Arruda L, Mayor R, Ng CKY, Weigelt B, Martínez-Ricarte F, Torrejon D, Oliveira M, Arias A, Raventos C, Tang J, Guerini-Rocco E, Martínez-Sáez E, Lois S, Marín O, de la Cruz X, Piscuoglio S, Towers R, Vivancos A, Peg V, Ramon y Cajal S, Carles J, Rodon J, González-Cao M, Tabernero J, Felip E, Sahuquillo J, Berger MF, Cortes J, Reis-Filho JS, Seoane J. Cerebrospinal fluid-derived circulating tumour DNA better represents the genomic alterations of brain tumours than plasma. Nat Commun. 2015;6:8839. https://doi.org/10.1038/ncomms9839.
De Mattos-Arruda L. Liquid biopsy for HER2-positive breast cancer brain metastasis: the role of the cerebrospinal fluid. ESMO Open. 2017;2:e000270. https://doi.org/10.1136/esmoopen-2017-000270.
Huang W-T, Lu N-M, Hsu W-Y, Chang S-E, Atkins A, Mei R, Javey M. CSF-ctDNA SMSEQ analysis to tailor the treatment of a patient with brain metastases: a case report. Case Rep Oncol. 2018;11:68–74. https://doi.org/10.1159/000486568.
Singh AP, Cheng H, Guo X, Levy B, Halmos B. Circulating tumor DNA in non-small-cell lung cancer: a primer for the clinician. JCO Precis Oncol. 2017. https://doi.org/10.1200/po.17.00054.
We are grateful to Katrien Van Den Bossche for her contribution to perform experiments and data management.
This scientific work was mainly funded by the non-profit organization “Kom Op Tegen Kanker” (http://www.komoptegenkanker.be/).
Ethics approval and consent to participate
Blood samples were prospectively collected after obtaining informed consent with an Ethical Committee (UZ Brussel) approved document.
Consent for publication
All parties consented for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Median and mean for BRAF/NRAS mutant copy number for responders (CR/PR) versus non-responders (SD/PD). Box plots showing the median (C) and mean (D) baseline ctDNA copy numbers compared in objective responders (CR/PR) and non-responders (SD/PD).
Comparison between the groups of patients with detectable versus undetectable ctDNA for OS/PFS and BRAF/NRAS mutant median copy number for responders (CR/PR) versus non-responders (SD/PD) during follow up at week 3, 6 and 9. Box plots detailing the ctDNA median copy numbers for responders (CR/PR) versus non-responders (SD/PD) at the second pembrolizumab cycle in week 3 (A), the third pembrolizumab cycle in week 6 (D) and the fourth pembrolizumab cycle in week 9 (G). Kaplan-Meyer curves comparing PFS (B, E, H) and OS (C, F, I) in patients with detectable or undetectable ctDNA levels at that time point, respectively at weeks 3, 6 and 9.
S1. Longitudinal monitoring of patients with progressive disease exclusively in the central nervous system. Table S2. Absolute values of longitudinal monitoring of patients with > 500 copies of mutated ctDNA/ml of plasma at baseline and/or at week 3.
S3. Kaplan–Meier survival curves for OS and PFS for patients with exclusively CNS PD and patients with BRAF/NRAS mutant copy number of > 500 copies/mL of plasma. (A) OS and PFS for the subgroup with PD exclusively in the CNS (n = 10), time is shown on the horizontal axis in weeks; (B) For Patient MEL36 the BRAFV600E/D mutation copy number is shown on the left vertical axis and the fractional abundance (% from the total cell-free DNA) on the right vertical axis. Time (in days) is shown on the horizontal axis. Illustrative CT or PET/CT images are shown on the upper part on the chart. The ctDNA was undetectable pretreatment with pembrolizumab when PD brain was observed on brain MRI and became detectable when the disease progressed as well in the visceral metastatic sites. (C) OS and PFS (in weeks) for the subgroup of patients with a baseline ctDNA copy number of > 500 copies/mL of plasma (n = 11) and (D) the subgroup of patients with a copy number of > 500 copies/mL of plasma at week 3 (n = 8).
About this article
Cite this article
Seremet, T., Jansen, Y., Planken, S. et al. Undetectable circulating tumor DNA (ctDNA) levels correlate with favorable outcome in metastatic melanoma patients treated with anti-PD1 therapy. J Transl Med 17, 303 (2019). https://doi.org/10.1186/s12967-019-2051-8