Next Article in Journal
The Relationship between Exercise Self-Efficacy, Intention, and Structural Barriers for Physical Activity after a Cancer Diagnosis
Next Article in Special Issue
Genomic Aberrations in Circulating Tumor DNAs from Palbociclib-Treated Metastatic Breast Cancer Patients Reveal a Novel Resistance Mechanism
Previous Article in Journal
How Phosphofructokinase-1 Promotes PI3K and YAP/TAZ in Cancer: Therapeutic Perspectives
Previous Article in Special Issue
Baseline Plasma Tumor DNA (ctDNA) Correlates with PSA Kinetics in Metastatic Castration-Resistant Prostate Cancer (mCRPC) Treated with Abiraterone or Enzalutamide
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Early Assessment of Chemotherapy Response in Advanced Non-Small Cell Lung Cancer with Circulating Tumor DNA

by
Stephanie J. Yaung
1,*,
Corinna Woestmann
2,
Christine Ju
3,
Xiaoju Max Ma
1,
Sandeep Gattam
3,
Yiyong Zhou
3,
Liu Xi
1,
Subrata Pal
3,
Aarthi Balasubramanyam
3,
Nalin Tikoo
4,
Claus Peter Heussel
5,6,7,
Michael Thomas
7,8,
Mark Kriegsmann
9,
Michael Meister
7,8,
Marc A. Schneider
7,8,
Felix J. Herth
7,8,
Birgit Wehnl
10,
Maximilian Diehn
11,
Ash A. Alizadeh
11,
John F. Palma
1 and
Thomas Muley
7,8
add Show full author list remove Hide full author list
1
Roche Sequencing Solutions, Inc., Pleasanton, CA 94588, USA
2
Roche Sequencing Solutions, Inc., 14473 Potsdam, Germany
3
Roche Molecular Systems, Inc., Pleasanton, CA 94588, USA
4
Alector, Inc., South San Francisco, CA 94080, USA
5
Diagnostic and Interventional Radiology with Nuclear Medicine, Thoraxklinik, University Hospital, 69126 Heidelberg, Germany
6
Diagnostic and Interventional Radiology, University Hospital, 69120 Heidelberg, Germany
7
Translational Lung Research Centre (TLRC) Heidelberg, Member of the German Centre for Lung Research (DZL), 69120 Heidelberg, Germany
8
Translational Research Unit, Thoraxklinik at Heidelberg University Hospital, 69126 Heidelberg, Germany
9
Institute of Pathology, University Hospital Heidelberg, 69120 Heidelberg, Germany
10
Roche Diagnostics GmbH, 82377 Penzberg, Germany
11
Stanford Cancer Institute, Stanford University, Stanford, CA 94305, USA
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(10), 2479; https://doi.org/10.3390/cancers14102479
Submission received: 22 April 2022 / Revised: 12 May 2022 / Accepted: 15 May 2022 / Published: 18 May 2022
(This article belongs to the Special Issue Cell-Free DNA as Prognostic and Predictive Biomarker in Solid Cancers)

Abstract

:

Simple Summary

An early assessment of response to treatment is crucial to informing appropriate therapeutic management. Using a plasma-only strategy, we measured changes in circulating tumor DNA (ctDNA) levels after one or two cycles of chemotherapy in 92 patients with advanced non-small-cell lung cancer (NSCLC) treated with first-line chemo- or chemoradiation therapies. A ≤50% decrease in ctDNA level after one cycle of chemotherapy was associated with shorter progression-free survival and overall survival. A ≤50% decrease in ctDNA level after two cycles of chemotherapy also had shorter survival. Our findings demonstrate that using liquid biopsies to measure early changes in ctDNA levels in response to chemotherapy may help identify non-responders before standard-of-care imaging in advanced NSCLC. Monitoring treatment efficacy earlier and accurately can enable more personalized regimens to improve patient outcomes.

Abstract

Monitoring treatment efficacy early during therapy could enable a change in treatment to improve patient outcomes. We report an early assessment of response to treatment in advanced NSCLC using a plasma-only strategy to measure changes in ctDNA levels after one cycle of chemotherapy. Plasma samples were collected from 92 patients with Stage IIIB-IV NSCLC treated with first-line chemo- or chemoradiation therapies in an observational, prospective study. Retrospective ctDNA analysis was performed using next-generation sequencing with a targeted 198-kb panel designed for lung cancer surveillance and monitoring. We assessed whether changes in ctDNA levels after one or two cycles of treatment were associated with clinical outcomes. Subjects with ≤50% decrease in ctDNA level after one cycle of chemotherapy had a lower 6-month progression-free survival rate (33% vs. 58%, HR 2.3, 95% CI 1.2 to 4.2, log-rank p = 0.009) and a lower 12-month overall survival rate (25% vs. 70%, HR 4.3, 95% CI 2.2 to 9.7, log-rank p < 0.001). Subjects with ≤50% decrease in ctDNA level after two cycles of chemotherapy also had shorter survival. Using non-invasive liquid biopsies to measure early changes in ctDNA levels in response to chemotherapy may help identify non-responders before standard-of-care imaging in advanced NSCLC.

1. Introduction

Lung cancer remains the leading cause of cancer-related mortality worldwide [1], with only 10–30% of patients surviving five years after diagnosis [2]. The latest analysis of epidemiological data on non-small cell lung cancer (NSCLC) in the United States suggests the introduction of more effective treatments and earlier treatment associated with earlier diagnosis may have increased the five-year survival estimate from 23.3% in 2014 to 26.4% in 2017 [3]. With the growing number of available treatment options, such as molecular targeted therapies and immunotherapies, monitoring response to therapy during initial cycles of treatment may enable an earlier change in treatment to further improve patient outcomes. However, treatment response in solid tumors is primarily based on radiological assessments, which are recommended from every 6 to 8 weeks or from 6 to 12 weeks for advanced NSCLC [4]. Performing serial computed tomography (CT) scans at short intervals from 3 to 4 weeks can detect changes in tumor size after one cycle of chemotherapy and potentially change clinical management [5], though standard methods such as Response Evaluation Criteria in Solid Tumors (RECIST) [6] are limited by practicability, radiation exposure, and challenges in morphological evaluation [7].
Liquid biopsies have emerged as a promising tool for precision oncology [8]. In 2018, the International Association for the Study of Lung Cancer (IASLC) issued a statement paper concluding that clinical implementation of liquid biopsy approaches is justified in several therapeutic settings relevant to NSCLC, particularly for initial molecular diagnosis when tumor tissue is unavailable and for sparing an invasive re-biopsy to determine subsequent therapy after progression on a targeted therapy [9]. Furthermore, quantifying changes in the levels of tumor-specific mutations in circulating tumor DNA (ctDNA) enables early response prediction and monitoring in various advanced solid tumors [10], including colorectal cancer [11,12,13] and breast cancer [14,15]. In advanced lung cancer, changes in ctDNA levels can be used as a biomarker to monitor early response to immunotherapy [16,17,18,19,20] or tyrosine kinase inhibitors (TKIs) [21,22,23,24].
Despite novel treatment advances, chemo- and chemoradiation therapy remain standard-of-care options in certain settings for treating advanced NSCLC [4]. A limited number of studies have explored the use of changes in ctDNA levels to monitor early response to chemotherapy. Previous work tracked changes in mutant levels of single genes in ctDNA, such as KRAS and TP53 [25,26,27], or examined ctDNA positivity rather than a change from baseline [28,29]. Moreover, analyses of ctDNA typically leverage prior knowledge of variants from matched tissue biopsies [29,30]. Since a tissue biopsy is often inaccessible in many late-stage lung cancer patients [31], we sought to determine whether variants found in baseline plasma ctDNA can also enable longitudinal monitoring without any prior knowledge of variants from matched tissue biopsies. Using a targeted next-generation sequencing (NGS) panel optimized for surveillance and monitoring of lung cancer, we hypothesized that a plasma-only approach measuring changes in ctDNA levels can assess early molecular response to chemo- or chemoradiation therapies in advanced NSCLC.

2. Materials and Methods

2.1. Study Design and Population

The Lung Cancer Multi-Marker Study was an observational, prospective clinical research study at Thoraxklinik, Heidelberg University Hospital, Germany. Eligible treatment-naïve subjects with histologically or cytologically proven late-stage adenocarcinoma or squamous cell carcinoma of the lung were enrolled in this study between December 2015 and March 2018 (Table 1). Blood samples from subjects were collected during regular clinical assessments at baseline, prior to the start of therapy, and throughout treatment. Response Evaluation Criteria in Solid Tumors guidelines (RECIST 1.1) [6] were followed when assessing the efficacy of treatment. The study was approved by the local ethics committee (S-611/2017) and health authority and carried out in accordance with Good Clinical Practice and the Declaration of Helsinki: Ethical Principles for Medical Research Involving Human Patients. Informed consent was obtained from all donors.
For our assessment of plasma-based measurements within the first one or two cycles of chemotherapy, we analyzed up to three plasma samples per subject (Figure 1A). The baseline (B0) plasma sample was the last sample collected prior to the treatment start date. The treatment start date was defined as the date of initiation of cycle 1 of chemotherapy in first line. The first post-treatment (P1) plasma sample was the first sample collected between 21 and 49 days after treatment started. The P1 sample had to be collected on or before the initiation of cycle 2 of chemotherapy in first-line, or on or before initiation of second-line therapy. The second post-treatment (P2) plasma sample was the first sample collected between 21 and 49 days after the start of chemotherapy cycle 2. The P2 sample had to be collected on or before the initiation of cycle 3 of chemotherapy in first-line or on or before initiation of second-line therapy. For all time points, samples were not included if the patient received any radiotherapy within 14 days prior to plasma collection.

2.2. ctDNA Sequencing and Bioinformatics Analysis

Cell-free DNA (cfDNA) was isolated from plasma samples and prepared for sequencing using the AVENIO ctDNA Surveillance Kit (Roche; for research use only; not for use in diagnostic procedures). The analytical quality of the extracted DNA was evaluated using a quantitative polymerase chain reaction (qPCR) method as described by Saelee et al. [32] to maximize cfDNA input into the hybrid capture-based NGS workflow, which was based on cancer personalized profiling by deep sequencing (CAPP-Seq) technology [33]. The workflow utilizes a 198-kb panel designed to cover recurrently mutated genomic loci in lung and colorectal cancers as well as genes in the U.S. National Comprehensive Cancer Network (NCCN) Guidelines. Analysis was carried out as described in the instructions of the kit. For the 223 samples from 92 patients used in this study, the median input mass was 25.7 ng (range 5–290 ng), median paired-end reads per sample were 54.9 million (range 31.9–164.9 million), median percentage of on-target reads was 69.0% (range 31.5–78.9%), and median unique molecular depth was 4747 (range 781–11,883).
NGS analysis was completed with the AVENIO Oncology Analysis Software [34] (Roche; for research use only; not for use in diagnostic procedures), which incorporates bioinformatics methods from CAPP-Seq technology and integrated digital error suppression [35] to remove random PCR errors and stereotypical errors from technical artifacts. Somatic calling was performed on the unfiltered list of single nucleotide variants (SNVs) from the software. The somatic calling approach used a previously described machine learning model [36] that takes into account the variation in variant allele fraction (AF) and information on common germline variants in public databases (Exome Aggregation Consortium [37] release 0.3.1, 1000 Genomes [38] phase3v5b, and dbSNP [39] build 144). The model was trained to remove germline variants based on the lower variation in AFs of germline variants than somatic variants in longitudinal plasma samples from the same subject. In addition to the tumor-specific SNVs, filtered insertion/deletions (indels) from the AVENIO Oncology Analysis Software were included in downstream analyses as tumor-specific variants for each subject.
We performed longitudinal ctDNA analysis and measured changes in mean mutant molecules per milliliter (MMPM) in post-treatment plasma compared to baseline by tracking the AF of the tumor-specific variants identified in the baseline plasma (Figure 1A). The MMPM for each variant was quantified as the variant AF multiplied by the extracted mass (ng) and 330, then divided by the plasma volume (mL). The factor of 330 is the estimated number of haploid genome equivalents per ng of DNA based on the molecular weight of the human genome. In ctDNA analyses not informed by baseline plasma information, we calculated a mean MMPM from all somatic variants identified in post-treatment plasma.

2.3. ctDNA Thresholds for Molecular Response

A pre-specified threshold of >50% vs. ≤50% decrease in ctDNA level from baseline was selected to classify molecular responders vs. molecular non-responders. Selection of the threshold was based on review of prior studies, which included ctDNA decreases of either >50% or >30% for immunotherapy response [16,18,19] or >98% or >90% for EGFR TKI response [22,23]. Zou et al. [20] used a threshold of >50% ctDNA decrease at 6 weeks after second-line chemotherapy or immunotherapy.
We also explored the possibility of establishing a threshold based on a single post-treatment sample, not informed by a subject-specific baseline sample. This approach to molecular monitoring of response is well-established in the clinical management of chronic myeloid leukemia treated with TKIs, where molecular response is measured by the log decrease in BCR-ABL1 transcript level from a standardized baseline [40,41]. Instead of requiring analysis of the baseline sample, analysis of ctDNA levels in a single post-treatment plasma could provide an assessment of treatment response. For this method, the P1 samples in our study were randomly split into 60% training and 40% test sets with stratification by stage, smoking history, performance status, and radiotherapy status. The optimal threshold identified by maximally selected LogRank statistics [42] in the training set for classifying responders vs. non-responders was then applied to the test set.

2.4. Statistical Analysis

Subject baseline characteristics were summarized using proportions for categorical variables and median and interquartile ranges (IQR) for continuous variables. Tumor response assessment on the first CT scan after the post-treatment plasma collection date was used to evaluate the association with mean MMPM change after one cycle of chemotherapy using Fisher’s exact test and calculation of the odds ratio and associated 95% confidence interval.
Kaplan–Meier estimates and log-rank tests were used to evaluate classification of subjects as molecular responders vs. non-responders based on the developed algorithms using pre-defined or optimal thresholds. The differences in median progression-free survival (PFS) and overall survival (OS) between classification groups were estimated based on the Kaplan–Meier method. Six-month PFS was defined as time from treatment initiation to the date of first documented progressive disease (PD) in first line or death within 182 days after treatment initiation. Subjects who did not have a documented PD or remained alive at 182 days were censored at either 182 days or the date last known alive at time of data cut, whichever was earlier. Twelve-month OS was defined as time from treatment initiation to the date of death from any cause within 365 days after treatment initiation. Subjects who remained alive at 365 days were censored at either 365 days or the date last known alive at time of data cut, whichever was earlier. PFS rate at the six-month time point and OS rate at the twelve-month time point were estimated using Kaplan–Meier methodology.
Association between the classification groups and survival (PFS and OS) was assessed using both unadjusted and adjusted Cox proportional hazards models. Univariable Cox models were performed to assess the association between each potential adjustment variable (histology, age, sex, smoking history [smoker vs. ex-smoker vs. never smoked], UICC [III vs. IV], ECOG at baseline [0 vs. 1 or 2], and radiotherapy status) and survival. Variables with clinical relevance (p-value < 0.1) were included as multivariable adjustments.
Data were analyzed using R 3.6.1 (The R Foundation for Statistical Computing Platform, Vienna, Austria). All tests were two-sided and considered statistically significant at p-value < 0.05 unless specified otherwise.

3. Results

3.1. Cohort Baseline Characteristics

A total of 223 plasma samples from 92 stage IIIB-IV NSCLC subjects were analyzed. In total, 80 subjects had a baseline (B0) and post-chemotherapy cycle 1 (P1) plasma available for analysis. In total, 51 subjects had a B0 and post-chemotherapy cycle 2 (P2) plasma available for analysis. Approximately 60% of subjects had lung adenocarcinoma, and 40% had lung squamous cell carcinoma (Table 1). Most subjects received first-line doublet chemotherapy with a platinum-based agent (carboplatin or cisplatin) in combination with vinorelbine or pemetrexed (Table S1).

3.2. Quantification of Circulating Tumor DNA in Plasma

P1 plasma was collected on treatment days 21–41 (median: 23) after chemotherapy initialization, and P2 plasma was collected on days 42–91 (median: 45) of chemotherapy. All subjects had at least one tumor-specific SNV, insertion, or deletion detected in baseline plasma, with a median of 10 somatic mutations that were most frequently found in lung cancer genes, including TP53, KRAS, and NPAP1 (Table S2). Each variant was quantified with the MMPM metric, and the ctDNA level in a plasma sample was summarized as the mean MMPM across the tumor-specific variants in that sample. The median mean MMPM in plasma samples was 37.4 at baseline, 7.6 at P1, and 4.9 at P2 (Figure 1B). The distribution of changes in mean MMPM at P1 and P2 compared to baseline shows that most subjects experienced a decrease in ctDNA levels after one cycle of chemotherapy, and this decrease persisted after two cycles of chemotherapy (Figure 1C).

3.3. Association of Change in ctDNA Level with Clinical Outcome

The tumor response assessment from the first available CT scan after treatment initiation (median of 20 days after P1 plasma collection) was compared with the change in ctDNA level between B0 and P1 plasma (Figure 2A). Most subjects who had a partial response had a decrease in ctDNA level. A >50% decrease (equivalent to <−50% change) in ctDNA level from baseline was significantly associated with partial response (OR 13.15, 95% CI 2.75–127.66, p < 0.001; Figure 2B).
Subjects with ≤50% decrease in ctDNA level after one cycle of chemotherapy had a lower 6-month progression-free survival rate (Figure 3A; 33% vs. 58%, HR 2.3, 95% CI 1.2 to 4.2, log-rank p = 0.009) and a lower 12-month overall survival rate (Figure 3B; 25% vs. 70%, HR 4.3, 95% CI 2.2 to 8.7, log-rank p < 0.001). Unadjusted Cox proportional hazards models were performed for each potential adjustment variable to determine which ones to include in multivariate models determining associations between treatment response classification and PFS and OS (Table S3). A significant association between a ≤50% decrease in ctDNA level and PFS was seen after adjusting for stage (HR 1.9, 95% CI 1.0 to 3.7, log-rank p = 0.041; Figure 4A). Similarly, after adjusting for radiotherapy (received at any time during the study), ECOG performance status at baseline, and smoking history, a significant association was observed between a ≤50% decrease in ctDNA level and OS (HR 4.43, 95% CI 2.14 to 9.2, log-rank p < 0.001; Figure 4B).
In the 51 subjects with both a baseline plasma and post-chemotherapy cycle 2 plasma, those with ≤50% decrease in ctDNA level in the P2 plasma had a lower 6-month progression-free survival rate (29% vs. 73%, HR 4.3, 95% CI 1.8 to 10, log-rank p < 0.001) and a lower 12-month overall survival rate (37% vs. 75%, HR 3.2, 95% CI 1.3 to 8.2, log-rank p = 0.009) (Figure 5).

3.4. Association of Single Time Point ctDNA Level with Clinical Outcome

The 80 subjects with a P1 plasma sample were randomly split into 60% training and 40% test sets with stratification by stage, smoking history, performance status, and radiotherapy status. For the stratification factors as well as other variables such as histology and age, there were no significant differences between the training and test sets (Figure S1A). By maximally selected LogRank statistics, an optimal threshold of 10 mean MMPM was identified for classifying OS in the 49 subjects in the training cohort. Subjects with a ≥10 mean MMPM had a worse 12-month OS (HR 4.7, 95% CI 1.7 to 13, log-rank p < 0.001) (Figure S1B). Using the same threshold for the test set of 31 subjects, those with ≥10 mean MMPM had worse 12-month OS (HR 4.6, 95% CI 0.99 to 21, log-rank p = 0.033) (Figure S1C). For PFS, an optimal threshold of 227 mean MMPM was identified. Subjects with ≥227 mean MMPM had worse 6-month PFS in the training cohort (HR 2.8, 95% CI 1.1 to 7.7, log-rank p = 0.031) and the test cohort (HR 5.5, 95% CI 1.6 to 19, log-rank p = 0.003) (Figure S2A). Given the distribution of mean MMPM values at P1 (Figure 1B), 227 is a very high cutoff that would classify only 14% of subjects as non-responders. Therefore, we also tried the threshold of 10 mean MMPM from the OS analysis to classify PFS. Although not statistically significant, subjects with ≥10 mean MMPM had worse 6-month PFS (Figure S2B).
There were not enough P2 samples (n = 51) for a train/test split, and, therefore, we did not analyze post-chemotherapy cycle 2 in single time point analyses.

4. Discussion

Our results suggest that ctDNA analysis using a targeted 197-gene NGS panel is a feasible approach for monitoring early treatment effect in subjects with advanced NSCLC treated with chemotherapy. We found that a ≤50% decrease in ctDNA levels after one or two cycles of chemotherapy was significantly associated with worse 6-month PFS and 12-month OS (Figure 3 and Figure 5). This finding is consistent with the results in Zou et al. [20], which used the same targeted NGS panel and reported that a ≤50% decrease from baseline in maximum AF at 6 weeks after second-line treatment initiation in the docetaxel cohort had worse OS, though the separation was not significant (p = 0.061). We observed a larger separation in the Kaplan–Meier curves in our study for a ≤50% decrease in ctDNA from baseline at the same time point (P2), possibly because the treatment effect could be more pronounced in a first-line setting.
Another study utilizing a 22-gene NGS panel [28] found that the absence of ctDNA, not the decrease in ctDNA level, from 4 to 8 weeks after treatment initiation was the best prognostic marker; however, the analysis was based on a heterogeneous cohort of patients treated with either EGFR TKI or chemotherapy, precluding a direct comparison of chemotherapy response assessment. The smaller gene panel may limit the ability to detect tumor burden in plasma. In fact, this may explain the smaller decreases in ctDNA levels observed in single-gene studies, such as Jiang et al. [27], where patients with partial responses averaged a 31% drop in TP53 ctDNA levels after two cycles of chemotherapy.
Patient coverage is also an important consideration; studies tracking KRAS in patients with advanced NSCLC could detect KRAS mutations in only 10% or 48% of baseline plasma [25,26]. In our study, 100% of the 92 subjects had at least one trackable baseline mutation with the 197-gene NGS panel. For subjects treated with targeted therapies, one can take advantage of the mutations in genes that are relevant to the corresponding targeted therapies for response monitoring [43]. For tumors with no known driver genes, there are no readily available tools to monitor the disease burden of these subjects due largely in part to the molecular heterogeneity of the disease. Earlier work described a way to monitor these subjects, but customized assays are necessary [44]. The current work demonstrated the feasibility of an NGS panel approach to achieve high mutation detection rates for lung cancers, which is consistent with other work that utilized CAPP-Seq technology [33,45].
Importantly, we relied on baseline plasma instead of tissue biopsy at diagnosis to identify tumor-specific variants to enable post-treatment monitoring. Previous studies on longitudinal monitoring have mainly used baseline tissue biopsies as a starting point to identify tumor-derived variants [29,33,44]. Since tissue is not always accessible in late-stage lung cancer subjects, a complementary approach is needed in the absence of a tissue biopsy. In this study, we demonstrate the feasibility of identifying tumor-specific variants in baseline plasma ctDNA that could enable longitudinal monitoring. Consistent with prior studies [33,45], most of the somatic variants we identified in the plasma of lung cancer subjects were in non-driver genes. However, the presence and level of these non-driver mutations correlated with disease burden and were useful for post-treatment response assessment.
If baseline plasma is unavailable, a single-time point analysis of ctDNA level in post-treatment plasma could monitor therapy response. We identified a potential threshold of 10 mean MMPM in plasma after one cycle of chemotherapy to classify responders vs. non-responders. The same threshold was not significantly associated with 6-month PFS, suggesting that predicting PFS but not OS after treatment initiation requires quantification relative to a subject-specific baseline. This is consistent with the fact that progressive disease per RECIST criteria is measured by comparing tumor size from a baseline CT scan. Longitudinal measurement relative to baseline has been explored with other methods as well, such as a >20% reduction in standardized uptake volume by positron emission tomography after one cycle of chemotherapy [46], a >27% drop in serum levels of cytokeratin 19 fragment (CYFRA 21-1) after one cycle of chemotherapy [47], or a >20% decrease in serum levels of both carcinoembryonic antigen (CEA) and CYFRA 21-1 after two cycles of chemotherapy [48]. Besides measuring changes in ctDNA levels in plasma, other approaches could be considered for a multimodal approach for predicting progression.

5. Conclusions

We found that quantification of changes in ctDNA level using a targeted NGS panel is a promising biomarker reflecting tumor dynamics and could identify subjects who had a higher risk of failing therapy one treatment cycle earlier than routine CT imaging. Moreover, an early assessment of the treatment effect can potentially be measured in plasma within the first treatment cycle. Our results demonstrate that ctDNA quantitation is potentially feasible for early assessment of the treatment response of solid tumors, especially when tissue is not available. We present a framework whereby longitudinal ctDNA analysis could enable real-time, non-invasive monitoring of cancer, providing novel insights into whether an early response assessed by changes in the mutant molecule level has the potential to predict treatment effect in a tissue-independent fashion. The cutoffs and parameters described in this study may be suitable for different clinical research settings depending on sample availability relative to treatment.

6. Patents

A.B., C.J., X.M.M., T.M., F.J.H., N.T., and B.W. have filed a patent (WO2019211418) on utilizing mutant molecule count for quantifying ctDNA.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14102479/s1, Figure S1: Association between overall survival and ctDNA level after one chemotherapy cycle, Figure S2: Association between progression-free survival and ctDNA level after one chemotherapy cycle; Table S1: Summary of First Line Treatment Description, Table S2: Most Frequently Mutated Genes in Baseline Plasma (n = 92), Table S3: Hazard Ratios and Confidence Intervals for Potential Adjustment Variables.

Author Contributions

Experiment conception and design: X.M.M., T.M. and J.F.P. Data collection: C.P.H., M.T., M.K., M.M., M.A.S., F.J.H. and B.W. Study selection and data analysis: S.J.Y., C.W., C.J., S.G., Y.Z., L.X., S.P., A.B., N.T., M.D. and A.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Roche Sequencing Solutions, Inc.

Institutional Review Board Statement

This study was approved by the Ethics Committee of the Heidelberg Medical Faculty (S-611/2017) and health authority and carried out in accordance with Good Clinical Practice and the Declaration of Helsinki: Ethical Principles for Medical Research Involving Human patients.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data analyzed are available from the corresponding author on reasonable request.

Acknowledgments

The authors would like to acknowledge the valuable contribution of Andrea Bopp, Christa Stolp, Ingrid Heinzmann-Groth, and others in the expert study nurse team at the Lungbiobank Heidelberg (member of the Biomaterialbank Heidelberg (BMBH) and the biobank platform of the German Center for Lung Research (DZL)). The authors would also like to thank Katrina Mayol and Joshua Lefkowitz from Roche Sequencing Solutions, Pleasanton, CA, US and Sandra Siemann, Christina Reinsch, Maria Lange, Mirjam Feldkamp, and Anja Blueher from Roche Sequencing Solutions, Potsdam, Germany for their valuable contributions to this publication.

Conflicts of Interest

C.P.H.: received consultation or other fees from Schering-Plough, Pfizer, Basilea, Boehringer Ingelheim, Novartis, Roche, Astellas, Gilead, MSD, Lilly, Intermune, and Fresenius; received lecture fees from Gilead, Essex, Schering-Plough, Astrazeneca, Lilly, Roche, MSD, Pfizer, Bracco, MEDA Pharma, Intermune, Chiesi, Siemens, Covidien, Pierre Fabre, Boehringer Ingelheim, Grifolis, Novartis, Basilea, and Bayer. M.T.: received grants, personal fees, and non-financial support from AstraZeneca; personal fees and non-financial support from Beigene; personal fees and non-financial support from Boehringer Ingelheim; grants, personal fees, and non-financial support from Bristol-Myers Squibb; personal fees and non-financial support from Celgene; personal fees and non-financial support from Chugai; personal fees and non-financial support from Daiichi Sankyo; personal fees and non-financial support from GlaxoSmithKline; personal fees and non-financial support from Janssen Oncology; personal fees and non-financial support from Lilly; grants, personal fees and non-financial support from Merck; personal fees and non-financial support from MSD; personal fees and non-financial support from Novartis; personal fees and non-financial support from Pfizer; grants, personal fees, and non-financial support from Roche; personal fees and non-financial support from Sanofi; grants, personal fees, and non-financial support from Takeda outside the submitted work. M.K.: received honoraria/fees from GSK for consultation. M.M.: received research grants for his institution from Roche Diagnostics. M.A.S.: received research grants for his institution from Roche Diagnostics. F.J.H.: received grants/research support from Roche Diagnostics and personal fees for advisory boards and lectures from Pulmonx, BTG, Uptake Medical, Olympus Medical, and Holaira. M.D.: served as a consultant for Roche, Genentech, AstraZeneca, Novartis, BioNTech, Boehringer Ingelheim, and Gritstone Oncology, is a shareholder in CiberMed and Forsight Diagnostics, and has received research support from Varian Medical Systems, AstraZeneca, and Genentech, outside the submitted work. A.A.A.: served as a consultant for Roche, Genentech, Janssen, Pharmacyclics, Gilead, Celgene, and Chugai, is a shareholder in CiberMed and FortySeven, and has received research support from Pfizer, outside the submitted work. T.M.: received grants, personal fees, and non-financial support from Roche Diagnostics GmbH, Penzberg, Germany, during the conduct of the study; grants from Cellzome, Heidelberg, Germany, grants from Oncohost, Binyamina, Israel, outside the submitted work. S.J.Y., C.W., C.J., X.M.M., S.G., Y.Z., L.X., S.P., A.B., and J.F.P. are employees of F. Hoffmann-La Roche. N.T. is an employee of Alector, Inc. B.W. is an employee of Roche Diagnostics GmbH, Germany.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLO-BOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA. Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Allemani, C.; Matsuda, T.; Di Carlo, V.; Harewood, R.; Matz, M.; Nikšić, M.; Bonaventure, A.; Valkov, M.; Johnson, C.J.; Estève, J.; et al. Global Surveillance of Trends in Cancer Survival 2000–14 (CONCORD-3): Analysis of Individual Records for 37 513 025 Patients Diagnosed with One of 18 Cancers from 322 Population-Based Registries in 71 Countries. Lancet 2018, 391, 1023–1075. [Google Scholar] [CrossRef] [Green Version]
  3. Ganti, A.K.; Klein, A.B.; Cotarla, I.; Seal, B.; Chou, E. Update of Incidence, Prevalence, Survival, and Initial Treatment in Patients with Non-Small Cell Lung Cancer in the US. JAMA Oncol. 2021, 7, 1824–1832. [Google Scholar] [CrossRef]
  4. National Comprehensive Cancer Network Non-Small Cell Lung Cancer (Version 3.2022). Available online: https://www.nccn.org/professionals/physician_gls/pdf/nscl.pdf (accessed on 31 March 2022).
  5. Bruzzi, J.F.; Truong, M.; Zinner, R.; Erasmus, J.J.; Sabloff, B.; Munden, R. Short-Term Restaging of Patients with Non-Small Cell Lung Cancer Receiving Chemotherapy. J. Thorac. Oncol. 2006, 1, 425–429. [Google Scholar] [CrossRef] [PubMed]
  6. Eisenhauer, E.A.; Therasse, P.; Bogaerts, J.; Schwartz, L.H.; Sargent, D.; Ford, R.; Dancey, J.; Arbuck, S.; Gwyther, S.; Mooney, M.; et al. New Response Evaluation Criteria in Solid Tumours: Revised RECIST Guideline (Version 1.1). Eur. J. Cancer 2009, 45, 228–247. [Google Scholar] [CrossRef] [PubMed]
  7. Ko, C.C.; Yeh, L.R.; Kuo, Y.T.; Chen, J.H. Imaging Biomarkers for Evaluating Tumor Response: RECIST and Beyond. Biomark. Res. 2021, 9, 52. [Google Scholar] [CrossRef]
  8. Siravegna, G.; Marsoni, S.; Siena, S.; Bardelli, A. Integrating Liquid Biopsies into the Management of Cancer. Nat. Rev. Clin. Oncol. 2017, 14, 531–548. [Google Scholar] [CrossRef]
  9. Rolfo, C.; Mack, P.C.; Scagliotti, G.V.; Baas, P.; Barlesi, F.; Bivona, T.G.; Herbst, R.S.; Mok, T.S.; Peled, N.; Pirker, R.; et al. Liquid Biopsy for Advanced Non-Small Cell Lung Cancer (NSCLC): A Statement Paper from the IASLC. J. Thorac. Oncol. 2018, 13, 1248–1268. [Google Scholar] [CrossRef] [Green Version]
  10. Boonstra, P.A.; Wind, T.T.; van Kruchten, M.; Schuuring, E.; Hospers, G.A.P.; van der Wekken, A.J.; de Groot, D.J.; Schröder, C.P.; Fehrmann, R.S.N.; Reyners, A.K.L. Clinical Utility of Circulating Tumor DNA as a Response and Follow-up Marker in Cancer Therapy. Cancer Metastasis Rev. 2020, 39, 999–1013. [Google Scholar] [CrossRef]
  11. Tie, J.; Kinde, I.; Wang, Y.; Wong, H.L.; Roebert, J.; Christie, M.; Tacey, M.; Wong, R.; Singh, M.; Karapetis, C.S.; et al. Circu-lating Tumor DNA as an Early Marker of Therapeutic Response in Patients with Metastatic Colorectal Cancer. Ann. Oncol. 2015, 26, 1715–1722. [Google Scholar] [CrossRef]
  12. Garlan, F.; Laurent-Puig, P.; Sefrioui, D.; Siauve, N.; Didelot, A.; Sarafan-Vasseur, N.; Michel, P.; Perkins, G.; Mulot, C.; Blons, H.; et al. Early Evaluation of Circulating Tumor DNA as Marker of Therapeutic Efficacy in Metastatic Colorectal Cancer Pa-tients (PLACOL Study). Clin. Cancer Res. 2017, 23, 5416–5425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ma, X.M.; Bendell, J.C.; Hurwitz, H.I.; Ju, C.; Lee, J.J.; Lovejoy, A.; Mancao, C.; Nicholas, A.; Price, R.; Sommer, N.; et al. Disease Monitoring Using Post-Induction Circulating Tumor DNA Analysis Following First-Line Therapy in Patients with Metastatic Colorectal Cancer. Clin. Cancer Res. 2020, 26, 4010–4017. [Google Scholar] [CrossRef] [Green Version]
  14. Dawson, S.-J.; Tsui, D.W.Y.; Murtaza, M.; Biggs, H.; Rueda, O.M.; Chin, S.-F.; Dunning, M.J.; Gale, D.; Forshew, T.; Mah-ler-Araujo, B.; et al. Analysis of Circulating Tumor DNA to Monitor Metastatic Breast Cancer. N. Engl. J. Med. 2013, 368, 1199–1209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Hrebien, S.; Citi, V.; Garcia-Murillas, I.; Cutts, R.; Fenwick, K.; Kozarewa, I.; McEwen, R.; Ratnayake, J.; Maudsley, R.; Carr, T.H.; et al. Early CtDNA Dynamics as a Surrogate for Progression-Free Survival in Advanced Breast Cancer in the BEECH Trial. Ann. Oncol. 2019, 30, 945–952. [Google Scholar] [CrossRef]
  16. Goldberg, S.B.; Narayan, A.; Kole, A.J.; Decker, R.H.; Teysir, J.; Carriero, N.J.; Lee, A.; Nemati, R.; Nath, S.K.; Mane, S.M.; et al. Early Assessment of Lung Cancer Immunotherapy Response via Circulating Tumor DNA. Clin. Cancer Res. 2018, 24, 1872–1880. [Google Scholar] [CrossRef] [Green Version]
  17. Guibert, N.; Jones, G.; Beeler, J.F.; Plagnol, V.; Morris, C.; Mourlanette, J.; Delaunay, M.; Keller, L.; Rouquette, I.; Favre, G.; et al. Targeted Sequencing of Plasma Cell-Free DNA to Predict Response to PD1 Inhibitors in Advanced Non-Small Cell Lung Cancer. Lung Cancer 2019, 137, 1–6. [Google Scholar] [CrossRef]
  18. Nabet, B.Y.; Esfahani, M.S.; Moding, E.J.; Hamilton, E.G.; Chabon, J.J.; Rizvi, H.; Steen, C.B.; Chaudhuri, A.A.; Liu, C.L.; Hui, A.B.; et al. Noninvasive Early Identification of Therapeutic Benefit from Immune Checkpoint Inhibition. Cell 2020, 183, 363–376. [Google Scholar] [CrossRef]
  19. Van der Leest, P.; Hiddinga, B.; Miedema, A.; Azpurua, M.L.A.; Rifaela, N.; ter Elst, A.; Timens, W.; Groen, H.J.M.; van Kempen, L.C.; Hiltermann, T.J.N.; et al. Circulating Tumor DNA as a Biomarker for Monitoring Early Treatment Responses of Patients with Advanced Lung Adenocarcinoma Receiving Immune Checkpoint Inhibitors. Mol. Oncol. 2021, 15, 2910–2922. [Google Scholar] [CrossRef]
  20. Zou, W.; Yaung, S.J.; Fuhlbrück, F.; Ballinger, M.; Peters, E.; Palma, J.F.; Shames, D.S.; Gandara, D.; Jiang, Y.; Patil, N.S. CtDNA Predicts Overall Survival in Patients With NSCLC Treated With PD-L1 Blockade or With Chemotherapy. JCO Precis. Oncol. 2021, 5, 827–838. [Google Scholar] [CrossRef]
  21. Kato, K.; Uchida, J.; Kukita, Y.; Kumagai, T.; Nishino, K.; Inoue, T.; Kimura, M.; Oba, S.; Imamura, F. Numerical Indices Based on Circulating Tumor DNA for the Evaluation of Therapeutic Response and Disease Progression in Lung Cancer Patients. Sci. Rep. 2016, 6, 29093. [Google Scholar] [CrossRef] [Green Version]
  22. Phallen, J.; Leal, A.; Woodward, B.D.; Forde, P.M.; Naidoo, J.; Marrone, K.A.; Brahmer, J.R.; Fiksel, J.; Medina, J.E.; Cristiano, S.; et al. Early Noninvasive Detection of Response to Targeted Therapy in Non–Small Cell Lung Cancer. Cancer Res. 2019, 79, 1204–1213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Cheng, M.L.; Lau, C.J.; Milan, M.S.D.; Supplee, J.G.; Riess, J.W.; Bradbury, P.A.; Jänne, P.A.; Oxnard, G.R.; Paweletz, C.P. Plasma CtDNA Response Is an Early Marker of Treatment Effect in Advanced NSCLC. JCO Precis. Oncol. 2021, 5, 393–402. [Google Scholar] [CrossRef] [PubMed]
  24. Angeles, A.K.; Christopoulos, P.; Yuan, Z.; Bauer, S.; Janke, F.; Ogrodnik, S.J.; Reck, M.; Schlesner, M.; Meister, M.; Schneider, M.A.; et al. Early Identification of Disease Progression in ALK-Rearranged Lung Cancer Using Circulating Tumor DNA Analysis. NPJ Precis. Oncol. 2021, 5, 100. [Google Scholar] [CrossRef] [PubMed]
  25. Nygaard, A.D.; Spindler, K.L.G.; Pallisgaard, N.; Andersen, R.F.; Jakobsen, A. Levels of Cell-Free DNA and Plasma KRAS during Treatment of Advanced NSCLC. Oncol. Rep. 2014, 31, 969–974. [Google Scholar] [CrossRef] [Green Version]
  26. Zulato, E.; Attili, I.; Pavan, A.; Nardo, G.; Del Bianco, P.; Boscolo Bragadin, A.; Verza, M.; Pasqualini, L.; Pasello, G.; Fassan, M.; et al. Early Assessment of KRAS Mutation in CfDNA Correlates with Risk of Progression and Death in Advanced Non-Small-Cell Lung Cancer. Br. J. Cancer 2020, 123, 81–91. [Google Scholar] [CrossRef]
  27. Jiang, T.; Li, X.; Wang, J.; Su, C.; Han, W.; Zhao, C.; Wu, F.; Gao, G.; Li, W.; Chen, X.; et al. Mutational, Andscape of CfDNA Identifies Distinct Molecular Features Associated with Therapeutic Response to First-Line Platinum-Based Doublet Chemo-therapy in Patients with Advanced NSCLC. Theranostics 2017, 7, 4753–4762. [Google Scholar] [CrossRef]
  28. Pécuchet, N.; Zonta, E.; Didelot, A.; Combe, P.; Thibault, C.; Gibault, L.; Lours, C.; Rozenholc, Y.; Taly, V.; Laurent-Puig, P.; et al. Base-Position Error Rate Analysis of Next-Generation Sequencing Applied to Circulating Tumor DNA in Non-Small Cell Lung Cancer: A Prospective Study. PLoS Med. 2016, 13, 1–19. [Google Scholar] [CrossRef]
  29. Jiang, J.; Adams, H.-P.; Lange, M.; Siemann, S.; Feldkamp, M.; McNamara, S.; Froehler, S.; Yaung, S.J.; Yao, L.; Bal-asubramanyam, A.; et al. Plasma-Based Longitudinal Mutation Monitoring as a Potential Predictor of Disease Progression in Subjects with Adenocarcinoma in Advanced Non-Small Cell Lung Cancer. BMC Cancer 2020, 20, 885. [Google Scholar] [CrossRef]
  30. Sholl, L.M.; Aisner, D.L.; Allen, T.C.; Beasley, M.B.; Cagle, P.T.; Capelozzi, V.L.; Dacic, S.; Hariri, L.P.; Kerr, K.M.; Lantuejoul, S.; et al. Liquid Biopsy in Lung Cancer: A Perspective from Members of the Pulmonary Pathology Society. Arch. Pathol. Lab. Med. 2016, 140, 825–829. [Google Scholar] [CrossRef] [Green Version]
  31. Gutierrez, M.E.; Choi, K.; Lanman, R.B.; Licitra, E.J.; Skrzypczak, S.M.; Pe Benito, R.; Wu, T.; Arunajadai, S.; Kaur, S.; Harper, H.; et al. Genomic Profiling of Advanced Non–Small Cell Lung Cancer in Community Settings: Gaps and Opportunities. Clin. Lung Cancer 2017, 18, 651–659. [Google Scholar] [CrossRef] [Green Version]
  32. Saelee, S.L.; Lovejoy, A.F.; Hinzmann, B.; Mayol, K.; Huynh, S.; Harrell, A.; Lefkowitz, J.; Deodhar, N.; Garcia-Montoya, G.; Yaung, S.J.; et al. Quantitative PCR Based Method to Assess CfDNA Quality, Adjust Input Mass, and Improve next-Generation Sequencing Assay Performance. J. Mol. Diagnostics 2022. [Google Scholar] [CrossRef] [PubMed]
  33. Newman, A.M.; Bratman, S.V.; To, J.; Wynne, J.F.; Eclov, N.C.W.; Modlin, L.A.; Liu, C.L.; Neal, J.W.; Wakelee, H.A.; Merritt, R.E.; et al. An Ultrasensitive Method for Quantitating Circulating Tumor DNA with Broad Patient Coverage. Nat. Med. 2014, 20, 548–554. [Google Scholar] [CrossRef] [PubMed]
  34. Deveson, I.W.; Gong, B.; Lai, K.; LoCoco, J.S.; Richmond, T.A.; Schageman, J.; Zhang, Z.; Novoradovskaya, N.; Willey, J.C.; Jones, W.; et al. Evaluating the Analytical Validity of Circulating Tumor DNA Sequencing Assays for Precision Oncology. Nat. Biotechnol. 2021, 39, 1115–1128. [Google Scholar] [CrossRef] [PubMed]
  35. Newman, A.M.; Lovejoy, A.F.; Klass, D.M.; Kurtz, D.M.; Chabon, J.J.; Scherer, F.; Stehr, H.; Liu, C.L.; Bratman, S.V.; Say, C.; et al. Integrated Digital Error Suppression for Improved Detection of Circulating Tumor DNA. Nat. Biotechnol. 2016, 34, 547–555. [Google Scholar] [CrossRef]
  36. Yao, L.; Lal, P.; Fang, L.-T.; Lee, J.; Palma, J.; Rosenthal, A.; Hinzmann, B.; Lovejoy, A.; Lam, H.Y.K. Abstract 5293: A Method to Identify Somatic Mutations from Tumor Samples in the Absence of Matched Normal Tissue. Cancer Res. 2018, 78, 5293. [Google Scholar] [CrossRef]
  37. Karczewski, K.J.; Weisburd, B.; Thomas, B.; Solomonson, M.; Ruderfer, D.M.; Kavanagh, D.; Hamamsy, T.; Lek, M.; Samocha, K.E.; Cummings, B.B.; et al. The ExAC Browser: Displaying Reference Data Information from over 60,000 Exomes. Nucleic Acids Res. 2017, 45, D840–D845. [Google Scholar] [CrossRef] [Green Version]
  38. Auton, A.; Abecasis, G.R.; Altshuler, D.M.; Durbin, R.M.; Bentley, D.R.; Chakravarti, A.; Clark, A.G.; Donnelly, P.; Eichler, E.E.; Flicek, P.; et al. A Global Reference for Human Genetic Variation. Nature 2015, 526, 68–74. [Google Scholar] [CrossRef] [Green Version]
  39. Smigielski, E.M.; Sirotkin, K.; Ward, M.; Sherry, S.T. DbSNP: A Database of Single Nucleotide Polymorphisms. Nucleic Acids Res. 2000, 28, 352–355. [Google Scholar] [CrossRef] [Green Version]
  40. Izzo, B.; Gottardi, E.M.; Errichiello, S.; Daraio, F.; Baratè, C.; Galimberti, S. Monitoring Chronic Myeloid Leukemia: How Mo-lecular Tools May Drive Therapeutic Approaches. Front. Oncol. 2019, 9, 833. [Google Scholar] [CrossRef] [Green Version]
  41. Hochhaus, A.; Baccarani, M.; Silver, R.T.; Schiffer, C.; Apperley, J.F.; Cervantes, F.; Clark, R.E.; Cortes, J.E.; Deininger, M.W.; Guilhot, F.; et al. European LeukemiaNet 2020 Recommendations for Treating Chronic Myeloid Leukemia. Leukemia 2020, 34, 966–984. [Google Scholar] [CrossRef] [Green Version]
  42. Hothorn, T.; Lausen, B. On the Exact Distribution of Maximally Selected Rank Statistics. Comput. Stat. Data Anal. 2003, 43, 121–137. [Google Scholar] [CrossRef]
  43. Chabon, J.J.; Simmons, A.D.; Lovejoy, A.F.; Esfahani, M.S.; Newman, A.M.; Haringsma, H.J.; Kurtz, D.M.; Stehr, H.; Scherer, F.; Karlovich, C.A.; et al. Circulating Tumour DNA Profiling Reveals Heterogeneity of EGFR Inhibitor Resistance Mechanisms in Lung Cancer Patients. Nat. Commun. 2016, 7, 11815. [Google Scholar] [CrossRef] [PubMed]
  44. Diehl, F.; Schmidt, K.; Choti, M.A.; Romans, K.; Goodman, S.; Li, M.; Thornton, K.; Agrawal, N.; Sokoll, L.; Szabo, S.A.; et al. Circulating Mutant DNA to Assess Tumor Dynamics. Nat. Med. 2008, 14, 985–990. [Google Scholar] [CrossRef] [PubMed]
  45. Chaudhuri, A.A.; Chabon, J.J.; Lovejoy, A.F.; Newman, A.M.; Stehr, H.; Azad, T.D.; Khodadoust, M.S.; Esfahani, M.S.; Liu, C.L.; Zhou, L.; et al. Early Detection of Molecular Residual Disease in Localized Lung Cancer by Circulating Tumor DNA Profiling. Cancer Discov. 2017, 7, 1394–1403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Novello, S.; Vavalà, T.; Levra, M.G.; Solitro, F.; Pelosi, E.; Veltri, A.; Scagliotti, G.V. Early Response to Chemotherapy in Patients with Non-Small-Cell Lung Cancer Assessed by [18F]-Fluoro-Deoxy-d-Glucose Positron Emission Tomography and Computed Tomography. Clin. Lung Cancer 2013, 14, 230–237. [Google Scholar] [CrossRef]
  47. Vollmer, R.T.; Govindan, R.; Graziano, S.L.; Gamble, G.; Garst, J.; Kelley, M.J.; Christenson, R.H. Serum CYFRA 21-1 in Ad-vanced Stage Non-Small Cell Lung Cancer: An Early Measure of Response. Clin. Cancer Res. 2003, 9, 1728–1733. [Google Scholar]
  48. Ardizzoni, A.; Cafferata, M.A.; Tiseo, M.; Filiberti, R.; Marroni, P.; Grossi, F.; Paganuzzi, M. Decline in Serum Carcinoem-bryonic Antigen and Cytokeratin 19 Fragment during Chemotherapy Predicts Objective Response and Survival in Patients with Advanced Nonsmall Cell Lung Cancer. Cancer 2006, 107, 2842–2849. [Google Scholar] [CrossRef]
Figure 1. Study design and ctDNA level distributions. (A) ctDNA sequencing was performed on baseline plasma (B0) and first available plasma after one cycle of chemotherapy (P1) and first available plasma after two cycles of chemotherapy (P2). ctDNA levels were measured as mean mutant molecules per milliliter (MMPM). (B) Distribution of ctDNA levels in B0, P1, and P2 samples. (C) Distribution of changes in ctDNA level at P1 from B0 and P2 from B0 in matched samples. Center line indicates the median value. **** Wilcoxon p ≤ 0.0001. ns: not significant, Wilcoxon p > 0.05.
Figure 1. Study design and ctDNA level distributions. (A) ctDNA sequencing was performed on baseline plasma (B0) and first available plasma after one cycle of chemotherapy (P1) and first available plasma after two cycles of chemotherapy (P2). ctDNA levels were measured as mean mutant molecules per milliliter (MMPM). (B) Distribution of ctDNA levels in B0, P1, and P2 samples. (C) Distribution of changes in ctDNA level at P1 from B0 and P2 from B0 in matched samples. Center line indicates the median value. **** Wilcoxon p ≤ 0.0001. ns: not significant, Wilcoxon p > 0.05.
Cancers 14 02479 g001
Figure 2. Association between first response assessment and ctDNA decrease after one chemotherapy cycle. (A) Waterfall plot showing change in ctDNA level at post-chemotherapy cycle 1 (P1) from baseline (B0). ctDNA levels were measured as mean mutant molecules per milliliter (MMPM). Each bar is colored by the first available response assessment after initiation of therapy: partial response (PR), stable disease (SD), progressive disease (PD), or CT scan not available (NA). Dotted line indicates the threshold of 50% decrease in mean MMPM. (B) Contingency table shows classification of response by ctDNA using a >50% decrease threshold vs. response by CT scan. p-value and odds ratio were calculated with Fisher’s Exact Test.
Figure 2. Association between first response assessment and ctDNA decrease after one chemotherapy cycle. (A) Waterfall plot showing change in ctDNA level at post-chemotherapy cycle 1 (P1) from baseline (B0). ctDNA levels were measured as mean mutant molecules per milliliter (MMPM). Each bar is colored by the first available response assessment after initiation of therapy: partial response (PR), stable disease (SD), progressive disease (PD), or CT scan not available (NA). Dotted line indicates the threshold of 50% decrease in mean MMPM. (B) Contingency table shows classification of response by ctDNA using a >50% decrease threshold vs. response by CT scan. p-value and odds ratio were calculated with Fisher’s Exact Test.
Cancers 14 02479 g002
Figure 3. Association between survival and ctDNA decrease after one chemotherapy cycle. Progression-free survival (A) and overall survival (B) for >50% vs. ≤50% decrease in mean mutant molecules per milliliter (MMPM) from baseline in the first available plasma after one cycle of chemotherapy.
Figure 3. Association between survival and ctDNA decrease after one chemotherapy cycle. Progression-free survival (A) and overall survival (B) for >50% vs. ≤50% decrease in mean mutant molecules per milliliter (MMPM) from baseline in the first available plasma after one cycle of chemotherapy.
Cancers 14 02479 g003
Figure 4. Multivariate Cox regression analysis of survival. Forest plots displaying hazard ratios from multivariate analysis of progression-free survival (A) and overall survival (B) for >50% vs. ≤50% decrease in mean mutant molecules per milliliter (MMPM) ctDNA. * p ≤ 0.05 *** p ≤ 0.001.
Figure 4. Multivariate Cox regression analysis of survival. Forest plots displaying hazard ratios from multivariate analysis of progression-free survival (A) and overall survival (B) for >50% vs. ≤50% decrease in mean mutant molecules per milliliter (MMPM) ctDNA. * p ≤ 0.05 *** p ≤ 0.001.
Cancers 14 02479 g004
Figure 5. Association between survival and ctDNA decrease after two chemotherapy cycles. Progression-free survival (A) and overall survival (B) for >50% vs. ≤50% decrease in mean mutant molecules per milliliter (MMPM) from baseline in the first available plasma after two cycles of chemotherapy.
Figure 5. Association between survival and ctDNA decrease after two chemotherapy cycles. Progression-free survival (A) and overall survival (B) for >50% vs. ≤50% decrease in mean mutant molecules per milliliter (MMPM) from baseline in the first available plasma after two cycles of chemotherapy.
Cancers 14 02479 g005
Table 1. Baseline characteristics of study subjects.
Table 1. Baseline characteristics of study subjects.
CharacteristicSubjects (n = 92)
Histology Subtype, n (%)
Adenocarcinoma56 (61)
Squamous cell carcinoma36 (39)
Age, median (IQR)65 (59, 71)
Sex, n (%)
Female23 (25)
Male69 (75)
Smoking History, n (%)
Ex-smoker48 (52)
Never smoked5 (5)
Smoker39 (42)
UICC, n (%)
IIIB29 (32)
IIIC9 (10)
IV34 (37)
IVA14 (15)
IVB6 (7)
ECOG, n (%)
055 (60)
136 (39)
21 (1)
Tumor Stage, n (%)
T13 (3)
T215 (16)
T331 (34)
T442 (46)
TX1 (1)
N Stage, n (%)
N04 (4)
N16 (7)
N235 (38)
N346 (50)
NX1 (1)
M Stage, n (%)
M034 (37)
M154 (59)
MX4 (4)
Note: IQR, interquartile range; UICC, Union for International Cancer Control; ECOG—0 = Fully active, able to carry on all pre-disease performance without restriction, 1 = Restricted in physically strenuous activity but ambulatory and able to carry out work of a light or sedentary nature, 2 = Ambulatory and capable of all self-care but unable to carry out any work activities; up and about more than 50% of waking hours.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yaung, S.J.; Woestmann, C.; Ju, C.; Ma, X.M.; Gattam, S.; Zhou, Y.; Xi, L.; Pal, S.; Balasubramanyam, A.; Tikoo, N.; et al. Early Assessment of Chemotherapy Response in Advanced Non-Small Cell Lung Cancer with Circulating Tumor DNA. Cancers 2022, 14, 2479. https://doi.org/10.3390/cancers14102479

AMA Style

Yaung SJ, Woestmann C, Ju C, Ma XM, Gattam S, Zhou Y, Xi L, Pal S, Balasubramanyam A, Tikoo N, et al. Early Assessment of Chemotherapy Response in Advanced Non-Small Cell Lung Cancer with Circulating Tumor DNA. Cancers. 2022; 14(10):2479. https://doi.org/10.3390/cancers14102479

Chicago/Turabian Style

Yaung, Stephanie J., Corinna Woestmann, Christine Ju, Xiaoju Max Ma, Sandeep Gattam, Yiyong Zhou, Liu Xi, Subrata Pal, Aarthi Balasubramanyam, Nalin Tikoo, and et al. 2022. "Early Assessment of Chemotherapy Response in Advanced Non-Small Cell Lung Cancer with Circulating Tumor DNA" Cancers 14, no. 10: 2479. https://doi.org/10.3390/cancers14102479

APA Style

Yaung, S. J., Woestmann, C., Ju, C., Ma, X. M., Gattam, S., Zhou, Y., Xi, L., Pal, S., Balasubramanyam, A., Tikoo, N., Heussel, C. P., Thomas, M., Kriegsmann, M., Meister, M., Schneider, M. A., Herth, F. J., Wehnl, B., Diehn, M., Alizadeh, A. A., ... Muley, T. (2022). Early Assessment of Chemotherapy Response in Advanced Non-Small Cell Lung Cancer with Circulating Tumor DNA. Cancers, 14(10), 2479. https://doi.org/10.3390/cancers14102479

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop