Main

Circulating cell-free DNA (cfDNA) is small DNA fragments circulating in the blood stream (Jahr et al, 2001). Elevated levels of cfDNA have been identified in cancer patients compared with healthy individuals (Leon et al, 1977). This elevation is believed to origin from both malignant tumour cells as well as from non-malignant cells constituting the tumour microenviroment (Mouliere and Thierry, 2012; Thierry et al, 2016). The precise mechanism for release of the DNA into the blood stream is not fully elucidated. It has been demonstrated that the amount of DNA in the circulation increases when cells are undergoing apoptosis or necrosis (Jahr et al, 2001). Other alternative mechanisms such as active secretion from tumour cells (Schwarzenbach et al, 2011) and release from circulating tumour cells (Bettegowda et al, 2014) have been suggested and may contribute to the total cfDNA level. In the recent years, it has been revealed that the contribution of cfDNA origin from tumour cells can be identified in the total cfDNA level by measuring cancer-specific mutations in the cfDNA also known to be present in the original tumour (Wan et al, 2017).

It has been proposed that the level of mutated cfDNA correlates with tumour burden. The level of tumour necrosis is thought to increase when tumour cells are rapidly proliferating as enlarging tumours are more likely to outgrow their blood supply (Leek et al, 1999). This causes tumour areas that are poorly vascularised and ischaemic, and if this ischaemia is sustained, ultimately necrosis will develop. As these necrotic areas are believed to release tumour DNA into the circulation, the amount of mutated cfDNA could be related to tumour burden (Diehl et al, 2005). Previous studies have supported this hypothesis by showing increased mutated cfDNA levels in advanced-stage patients compared with early-stage patients (Fleischhacker and Schmidt, 2007; Bettegowda et al, 2014; Newman et al, 2014). Moreover, correlations have been identified between changes in the amount of mutated cfDNA and changes in treatment response (Dawson et al, 2013; Newman et al, 2014; Sorensen et al, 2014) as well as relapse of disease (Tie et al, 2015). In addition to tumour volume, tumour aggressiveness may also be an important factor affecting the level of mutated cfDNA. An association between increasing necrosis grade and higher tumour grade has been identified, indicating that the amount of mutated cfDNA is associated with an aggressive tumour phenotype (Leek et al, 1999).

18F-fluoro-D-glucose (18F-FDG) positron emission tomography/computed tomography (PET/CT) is a well-established nuclear imaging modality for use in cancer patients. 18F-FDG is a radiolabelled glucose analogue accumulating in tissue with high glucose utilisation. FDG activity thereby reflects the metabolic activity and aggressiveness of a tumour. In addition, data on tumour volume can be assessed on an 18F-FDG-PET/CT scan. New volume-based PET parameters, as tumour lesion glycolysis (TLG), have been introduced capable of measuring the entire metabolic tumour burden (MTB) by incorporating both metabolic activity and volumetric data on all tumours in the patient (Larson et al, 1999). Thus TLG assessment on an 18F-FDG-PET/CT scan reflects both tumour burden and, to some extent, tumour aggressiveness of the disease and could be a more informative parameter for assessing the total active tumour burden in patients than tumour burden assessment on a CT scan.

Identification of oncogenic alterations in cfDNA is challenging due to the low concentration of target cfDNA in plasma. With the emergence of next-generation sequencing (NGS), it has become possible to assess a large amount of different mutations from a low input of DNA. NGS-based targeted panels directed against the most frequent alterations in cancer have been designed. In non-small cell lung cancer (NSCLC), a 22-gene NGS panel has been developed, which includes the most frequent mutations in NSCLC (Ion AmpliSeq Colon and Lung Cancer Panel v2). This panel has proven useful for mutation detection in plasma cfDNA (Frenel et al, 2015; Zonta et al, 2016).

In the present study, we evaluated the correlation between the amount of mutated cfDNA in plasma, and the MTB was measured by TLG in advanced-stage NSCLC patients. Moreover, we evaluated whether the presence of mutated cfDNA was correlated with patient survival.

Materials and methods

Patients

In this retrospective study, advanced-stage NSCLC patients prospectively enrolled in a single-centre study from April 2013 to August 2015 at the Department of Oncology, Aarhus University Hospital, Denmark were evaluated. Details on inclusion and exclusion criteria have been described previously (Winther-Larsen et al, 2016). At the time of inclusion, patients underwent an 18F-FDG-PET/CT scan and a blood sample was drawn. All patients were treated with erlotinib in a palliative setting. Testing for somatic epidermal growth factor receptor (EGFR) mutations had been performed on the diagnostic tumour biopsy in all adenocarcinoma patients as part of the routine diagnostic work-up (see Supplementary S1). For the purpose of this study, we included only patients from the total cohort who were EGFR mutation wild type (wt), had progression on first or second line of chemotherapy, and in whom both an 18F-FDG-PET/CT scan and a blood sample were available and evaluable. Selection of patients is illustrated in Supplementary S2. All patients gave informed written consent before inclusion and the study was approved by the Central Denmark Region Committees on Biomedical Research Ethics (no. 1-10-72-19-12). The study was reported to ClinicalTrial.gov (NCT02043002).

18F-fluoro-D-glucose positron emission tomography/computed tomography

All 18F-FDG-PET/CT scans were performed on a combined PET/CT scanner (Siemens Biograph TruePoint 40, Siemens Healthcare GMbH, Erlangen, Germany) at the Department of Nuclear Medicine and PET-Centre, Aarhus University Hospital, Denmark. The imaging protocol is described in Supplementary S1. Same scanner model, protocol for acquisition, and reconstruction software was used for all patients. An experienced nuclear medicine physician blinded to the patient outcome analysed all PET/CT scans using the Siemens Syngo.via software, Siemens Healthcare GMbH. In all scans, TLG was calculated for all evaluable lesions according to the Positron Emission Tomography Response Criteria in Solid Tumours 1.0 guideline (Wahl et al, 2009) (described in Supplementary S1 and illustrated in Supplementary S3). The sum of TLG for all evaluable lesions was defined as the MTB.

DNA extraction from blood

A peripheral blood sample of 10 ml was collected in an EDTA-containing tube for each patient. Within 2 h of withdrawal, samples were centrifuged (1400 g for 15 min) and plasma was isolated from the peripheral blood cells. Plasma was subsequently frozen at −80 °C until further analysis. cfDNA was purified from 2 ml of plasma by use of the QIAamp Circulating Nucleic Acid Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol and eluted in a volume of 100 μl TE buffer. The amount of cfDNA was quantified by measuring the beta-2-microglobulin gene by droplet digital PCR (Pallisgaard et al, 2015) as described in Supplementary S1.

Next-generation sequencing

NGS was carried out on all cfDNA samples. The Oncomine Solid Tumor DNA Kit was used to prepare libraries on 1.1–10 ng of cfDNA following the manufacturer’s instructions (Thermo Fisher Scientific, Waltham, MA, USA). NGS was conducted using the Ion Chef Instrument and the Ion Personal Genome Machine (PGM) System (Thermo Fisher Scientific). Sequencing was performed using Ion 316 v2 BC chips containing eight cfDNA samples per chip. Primary data processing was carried out using the Torrent Suite Software (version 5.0.4) on a Torrent Server and the Ion AmpliSeq Colon and Lung Cancer panel v2 template (all Thermo Fisher Scientific). Variant calling was performed using the Ion Reporter Software (version 5.0) and the AmpliSeq CHPv2 peripheral/CTC/CF DNA single sample workflow (Thermo Fisher Scientific). Default settings were used except from the reference and hotspot BED files, which were supplied in the kit. Sequencing was considered successful if the mean sequencing depth was ⩾2000. If this criterion could not be met, the sample was disqualified. Called variants were only accepted if the allele frequency (AF) was ⩾1%. The Integrative Genomics Viewer (Broad Institute, Cambridge, MA, USA) was used for visualisation of variants (Thorvaldsdóttir et al, 2013).

Statistical analysis

Correlations between detection of mutated cfDNA and clinical characteristics were calculated using χ2 or Fisher’s exact tests, where appropriate. Kruskal–Wallis test was used to calculate the correlation between the presence of mutated cfDNA and the levels of TLG and total cfDNA. Correlation between TLG and concentration of mutated cfDNA was evaluated using the Pearson’s correlation coefficient after log-transformation of both variables. Follow-up time was calculated using the reverse Kaplan–Meier method. Overall survival (OS) was measured from the start of enrolment until death of any cause or last follow-up date (19 April 2016). If patients were still alive on the last follow-up date, they were censored at that day. Estimates of median OS were calculated by Kaplan–Meier method and compared by the log-rank test. The Cox proportional hazards model was used to calculate crude and adjusted hazard ratios (HRs). Clinical variables were dichotomised except for age (continuous). All tests were two-sided, and P-values <0.05 were considered to be statistically significant. Statistical analyses were performed using SPSS statistics version 20.0 for windows (IBM SPSS Statistics, Chicago, IL, USA). STATA version 13 (Stata Corporation, College Station, TX, USA) was used for preparation of Kaplan–Meier curve and GraphPad Prism 7.0 (GraphPad Software Inc., San Diego, CA, USA) for preparation of scatter plot.

Results

Patients

A total of 46 patients were analysed. Patient characteristics are shown in Table 1. All patients had stage IV disease except for three patients with stage III disease. The majority of patients had progression on first-line chemotherapy. No patients were lost to follow-up. The median time interval between blood sampling and PET/CT scanning was 2 days (range: 0–9 days). After a median follow-up time of 20.9 months (95% CI: 11.2–30.6 month), five patients were still alive. The median OS in all patients was 5.9 months (95% CI: 2.9–8.9 months).

Table 1 Patient and tumour characteristics (N=46)

Correlation between the level of mutated cfDNA and MTB

cfDNA isolation was successful in all patients and a median value of 4587 copies ml−1 (range: 297–170 720 copies ml−1) was measured. Sequencing of cfDNA succeeded in 41 out of 46 patients (89%). A least one mutation was detected in 24 patients (59%). The number of identified mutations in each patient ranged from 1 to 4 with a median of 1 mutation. In total, 43 mutations were identified. The most frequently mutated genes were: TP53 (18 patients), KRAS (12 patients), and EGFR (2 patients; Table 2). The median AF of the most frequent mutation in each patient was 2.7% (range: 1.1–62.5%).

Table 2 Mutations identified in circulating cell-free DNA by use of next-generation sequencing (N=24)

TLG was assessable in all FDG-PET/CT. A significant correlation between the AF of the most frequent mutation (%) and the level of TLG was found as shown in Figure 1 (r=0.628, P=0.001). By contrast, no significant correlation could be identified between MTB and the number of mutations identified in each patient (data not shown).

Figure 1
figure 1

Scatter plot showing the correlation between the allele frequency of the most frequent mutation and the level of TLG in patients with circulating mutated cell-free DNA present in plasma ( N =24). The P-value was calculated using the Pearson’s correlation coefficient.

Evaluation of patients with and without detectable mutated cfDNA

No mutations could be identified in the cfDNA of 17 patients. A comparison of clinical and tumour characteristics between the mutated cfDNA-positive and the mutated cfDNA-negative group was performed, and no significant differences were identified (Supplementary S4). In addition, no significant difference in the level of total cfDNA was found between the mutated cfDNA-positive (5170 copies ml−1, range 1518–170 720) and the mutated cfDNA-negative (4510 copies ml−1, range 1232–39 160) patients (χ2 (2)=0.018, P=0.895). However, a significant higher median MTB was found in the mutated cfDNA-positive patients (689 (range 73–7034)) compared with the mutated cfDNA-negative patients (269 (range 67–3958)) (P=0.047).

When evaluating OS, patients with mutated cfDNA presence in the blood had a significant shorter median OS (3.7 months (95% CI: 0.0–7.6)) compared with mutated cfDNA-negative patients (10.6 months (95% CI: 7.1–14.1), P=0.019; Figure 2). In a univariate Cox regression analysis, the presence of mutated cfDNA together with high level of MTB was found to be associated with decreased OS (Table 3). To evaluate the independent impact of the presence of mutated cfDNA on survival, a multivariate Cox regression analysis was performed. Being mutated cfDNA positive remained an independent predictor of shorter OS with an adjusted HR of 2.11 (95% CI: 1.02–4.38) (Table 3).

Figure 2
figure 2

Kaplan–Meier curve for overall survival (OS) according to the presence or absence of circulating mutated cell-free DNA (cfDNA) in plasma. A median OS of 3.7 months (95% confidence interval (CI): 0–7.6 months) was found in patients with detectable mutated cfDNA, while a median OS of 12.6 months (95% CI: 8.4–16.9 months) was found in patients without detectable mutated cfDNA. The difference between the groups was calculated using the log-rank test.

Table 3 Univariate and multivariate Cox regression analyses of overall survival (N=46)

No difference in OS was found between mutated cfDNA-positive patients, with only one mutation identified compared with a higher number of mutations (P=0.355). Neither was any difference in survival found between mutated cfDNA-positive patients with an AF above or below the median 2.7% (P=0.949).

Discussion

Our study shows a significant correlation between the frequency of mutated cfDNA and the MTB measured by PET/CT in advanced-stage NSCLC patients. The study thereby supports the hypothesis of mutated cfDNA as a surrogate marker of tumour burden and tumour aggressiveness. In addition, we observed that mutated cfDNA-positive patients had a higher median MTB than patients without mutated cfDNA present in plasma and could confirm from previous reports showing that mutated cfDNA positivity is correlated to an inferior survival.

To our knowledge, only one previous study has, similar to ours, compared the exact level of mutated cfDNA with radiological tumour burden in NSCLC patients. Newman et al (2014) evaluated the correlation between amount of mutated cfDNA and tumour volume, defined as visible tumour on a CT scan, in nine newly diagnosed or recurrent NSCLC patients with different stages (I–IV). In line with ours, they found a significant correlation. In a study by Pécuchet et al (2016), the concentration of mutated cfDNA found in 75 advanced-stage NSCLC patients was divided in tertiles and correlated to the tumour burden defined as the sum of the RECIST target lesions. A significant higher median tumour burden was found in the groups with high and intermediate concentration of mutated cfDNA compared with the low group. In addition, there was also a significantly higher median tumour burden when comparing the group with high concentration of mutated cfDNA to the group with intermediate concentration. Both of these studies used a CT-defined tumour volume and no measurement of tumour aggressiveness was therefore incorporated. In contrast, we measured the MTB by use of the PET parameter TLG. This parameter not only measures the entire tumour burden by incorporating all measurable tumours in the patient (and not only RECIST target lesions) but also integrates the metabolic activity of each tumour. Thereby, the parameter becomes a measurement of both tumour burden and aggressiveness of the tumour. MTB has previously been used for assessment of tumour burden in a study evaluating the correlation between tumour burden and total amount of total cfDNA in 53 advanced NSCLC patients (Nygaard et al, 2014). No correlation could be found in this study indicating that total cfDNA is a more unspecific marker of tumour burden than mutated cfDNA. The majority of cfDNA originates from non-malignant tissue. Hence, the total amount of cfDNA not simply reflects tumour burden and tumour aggressiveness but rather reflect much more complicated biological mechanisms yet not fully understood.

Additionally, we found that patients with detectable mutated cfDNA in plasma had an inferior median survival compared with patients without mutated cfDNA present in plasma. A prognostic value of mutated cfDNA has previously been evaluated in several types of cancers, such as colorectal (Lecomte et al, 2002), breast (Dawson et al, 2013; Bettegowda et al, 2014), pancreatic (Hadano et al, 2016), and gastric (Gao et al, 2016). In NSCLC patients, a study of EGFR mutation-positive patients treated in the EURTAC trail found that the subgroup of L858R-positive patients with detectable mutations in plasma had a significant shorter OS as compared with L858R-positive patients with no plasma mutations, also supporting our finding (Karachaliou et al, 2015). However, the importance of mutated cfDNA positivity in EGFR mutation wt patients is less investigated. A recent study performed NGS of cfDNA and tumour samples in 105 unselected advanced NSCLC patients using the same NGS panel as used in the present study (Pécuchet et al, 2016). Mutations in cfDNA were identified in 71% of patients and a significantly reduced OS was seen in the mutated cfDNA-positive patients. Though, 43% of the patients harboured an EGFR mutation either in the tumour or in the cfDNA making the patient cohort heterogeneous. In our cohort, all patients had tumours that were EGFR wt and only 2 patients (2%) showed EGFR mutations in plasma. Owing to the homogenous nature of our cohort, we could show that mutated cfDNA positivity is a valuable tool for predicting survival in advanced EGFR mutation wt patients. The prognostic value of mutated cfDNA positivity could in theory be explained by the close correlation to tumour burden. However, we found that the impact of mutated cfDNA on survival was independent of the MTB in the multivariate analysis, indicating that mutated cfDNA positivity is an independent marker of survival.

Our study is the first of its kind and strengthened by the prospective nature of the original study allowing the PET scanning and the blood sampling to be performed in very close proximity to each other. The same PET scanner model, protocol for acquisition, and reconstruction software was used for all patients, thereby reducing the possible interindividual variability of the scans. Additionally, we had a homogeneous patient cohort receiving the same treatment and with complete clinical data in all patients. On the other hand, our study had some limitations to consider. Even though our study is the largest in the field, the number of patients is limited and our results are primarily hypothesis generating. Furthermore, we performed a targeted sequencing using a 22-gene panel with the risk that rare lung cancer mutations would not be detected. As we used the AF for quantification of mutated DNA, a possibility exist that a mutation could have been rejected (if the AF<1%) owing to a high amount of non-mutated cfDNA in the sample. Finally, we chose to define the most frequently mutation as the mutation defining the level of mutated cfDNA. This may underestimate the true level, but it would be more problematic to combine the frequency of different mutation types as it is not possible to determine whether different mutations exist in the same clone or arise from several different clones contributing their individual mutations.

In conclusion, we showed for the first time that a correlation exists between the level of mutated cfDNA and MTB in EGFR mutation wt NSCLC patients with advanced-stage disease. This indicates that mutated plasma cfDNA can be used as an indirect measure of tumour biology. Furthermore, we observed an inferior survival in patients with mutated cfDNA present in plasma compared with patients without. This correlation was independent of tumour burden and, if further validated, mutated cfDNA positivity could be a valuable tool for prediction of survival in patients with advanced NSCLC in future.