Despite the initial successes of immunotherapy, there is an urgent clinical need for molecular assays that identify patients more likely to respond. Here, we report that ultrasensitive measures of circulating tumor DNA (ctDNA) and T-cell expansion can be used to assess responses to immune checkpoint blockade in metastatic lung cancer patients (N = 24). Patients with clinical response to therapy had a complete reduction in ctDNA levels after initiation of therapy, whereas nonresponders had no significant changes or an increase in ctDNA levels. Patients with initial response followed by acquired resistance to therapy had an initial drop followed by recrudescence in ctDNA levels. Patients without a molecular response had shorter progression-free and overall survival compared with molecular responders [5.2 vs. 14.5 and 8.4 vs. 18.7 months; HR 5.36; 95% confidence interval (CI), 1.57–18.35; P = 0.007 and HR 6.91; 95% CI, 1.37–34.97; P = 0.02, respectively], which was detected on average 8.7 weeks earlier and was more predictive of clinical benefit than CT imaging. Expansion of T cells, measured through increases of T-cell receptor productive frequencies, mirrored ctDNA reduction in response to therapy. We validated this approach in an independent cohort of patients with early-stage non–small cell lung cancer (N = 14), where the therapeutic effect was measured by pathologic assessment of residual tumor after anti-PD1 therapy. Consistent with our initial findings, early ctDNA dynamics predicted pathologic response to immune checkpoint blockade. These analyses provide an approach for rapid determination of therapeutic outcomes for patients treated with immune checkpoint inhibitors and have important implications for the development of personalized immune targeted strategies.

Significance: Rapid and sensitive detection of circulating tumor DNA dynamic changes and T-cell expansion can be used to guide immune targeted therapy for patients with lung cancer.

See related commentary by Zou and Meyerson, p. 1038

Despite the durable clinical benefit observed with immune checkpoint inhibitors for patients with non–small cell lung cancer (NSCLC), the majority of patients are either refractory or eventually develop acquired resistance after an initial response (1). Similar to the targeted therapy paradigm, success of immuno-oncology seems to depend on choosing patient populations most likely to benefit. The plasticity of the immune system under immunotherapy has weakened single biomarker-driven approaches (2) and currently used predictive biomarkers have been unable to accurately identify the subset of patients who benefit from these therapies.

We hypothesized that noninvasive molecular analyses that evaluate tumor-derived cell free circulating tumor DNA (ctDNA) and tumor-extrinsic [T-cell receptor (TCR) repertoire] parameters may be useful for rapidly determining which patients would ultimately benefit from immune checkpoint blockade. Such approaches may be of particular importance for immune-targeted agents as the therapeutic responses have been challenging to evaluate using radiographic imaging due to tumor immune infiltration (3). Conventional response criteria such as the RECIST do not consistently capture the unique patterns and timing of antitumor immune responses (4, 5).

The temporal relationship between detection of ctDNA and emergence of recurrent or progressive disease has been shown in patients with early-stage NSCLC (6, 7) and as we show in our companion study in patients with advanced NSCLC receiving targeted therapies (8). During treatment with immunotherapy, our group has shown that ctDNA may be predictive of outcome in melanoma patients treated with CTLA-4 blockade (9). ctDNA changes have been associated with therapeutic outcome during immune checkpoint blockade in NSCLC (10–12), however these analyses have been limited by the low sensitivity of the approaches, permitting analyses in approximately half of the cases analyzed. Even less is known about the dynamics of the peripheral T-cell repertoire during immune checkpoint blockade in NSCLC (13) and how these changes relate to ctDNA levels and tumor response.

To overcome these issues and to allow ultrasensitive evaluation of ctDNA during therapy, we have developed targeted error-correction sequencing (TEC-Seq), a custom capture and sequencing approach that permits sensitive and specific detection of low abundance sequence alterations using next-generation sequencing (14). We have also developed new methods of evaluating TCR clonal expansion in the tumor microenvironment during immune checkpoint blockade (13). Here, we use these approaches to investigate whether ctDNA and TCR dynamics are reflective of therapeutic outcome for patients with NSCLC treated with immune checkpoint blockade.

Patient characteristics

Our study group consisted of 24 metastatic patients with NSCLC treated with immune checkpoint blockade as a standard of care (n = 19) or in the setting of a clinical trial (n = 5) between October 2014 and August 2016 at Johns Hopkins Sidney Kimmel Cancer Center. In parallel, we evaluated 14 patients with stage I to IIIA surgically resectable NSCLC that received anti-PD1 therapy in the setting of a neoadjuvant nivolumab clinical trial (15). The studies were conducted in accordance with the Declaration of Helsinki, were approved by the Institutional Review Board (IRB), and patients provided written informed consent for sample acquisition for research purposes. Clinical characteristics for all patients are summarized in Supplementary Table S1.

Treatment and assessment of therapeutic response

Therapeutic responses were evaluated by the RECIST version 1.1 (16). Baseline disease burden was determined by the sum of the longest diameters of target lesions as determined by RECIST 1.1 criteria. After baseline imaging, radiographic evaluation was performed at 5- to 10-week intervals or as clinically indicated for the metastatic cohort and 7 days prior to surgery for the early-stage cohort. The timing of radiologic assessments typically followed the early timepoints of blood sample collection. Although this approach may be subject to lead-time bias, we sought to mirror the imaging schedule used in clinical practice. Furthermore, in contrast to chemotherapy or targeted therapy, where therapeutic response may be accurately evaluated by imaging early after treatment initiation, the unique nature and timing of response to immune checkpoint blockade mandates response assessments at later timepoints including confirmation of the radiologic response. Of the 19 metastatic patients with NSCLC with detectable ctDNA, one achieved complete response (CGLU111), three patients achieved partial response (CGLU135, CGLU337, and CGLU347), and 12 achieved SD (CGLU115, CGLU117, CGLU159, CGLU160, CGLU162, CGLU168, CGLU203, CGLU211, CGLU212, CGLU340, CGLU351, and CGLU357) as best overall response. Three patients (CGLU121, CGLU243, and CGLU348) experienced disease progression. Of the three patients with partial response, two eventually developed molecular resistance. In the neoadjuvant cohort, a repeat chest CT ≤7 days prior to surgery revealed stable disease for all patients with detectable ctDNA at baseline.

Progression-free survival (PFS) and overall survival (OS) were defined as the time elapsed between the date of treatment initiation and the date of disease progression or death from disease, or the date of death, respectively (Supplementary Table S1). For the early stage cases with detectable ctDNA, two patients demonstrated a major pathologic response (pMPR defined as ≥90% decrease in tumor burden; CGLU206 and CGLU249), three patients had a partial pathologic response (at least 30% decrease in the tumor burden; CGLU205, CCGLU219, and CGLU221) and two patients had a pathologic nonresponse (CGLU222 and CGLU225).

Blood sample collection

For all patients, at least two serial blood samples (range 2–8) were collected over the course of treatment for isolation of plasma and extraction of cell-free DNA for genomic analyses. We analyzed a total of 105 serial plasma samples that were obtained prior to anti-PD1, at 4 to 8 weeks and additional time points during therapy for all metastatic patients with NSCLC except for CGLU135 and CGLU161. For these two patients, blood from early timepoints was not available and blood samples from the time of radiographic response and the time of acquired resistance were analyzed. A detailed description of the time points analyzed is shown in Supplementary Tables S2 and S3. Baseline tumors were analyzed by whole exome sequencing or targeted next-generation sequencing for patients in the metastatic cohort, with the exception of CGLU168, for which a tumor specimen from the time of resistance to immune checkpoint blockade was used. For the early-stage cohort, tumor samples prior to therapy initiation or at the time of resection in the cases where baseline tumor was not available, were analyzed by whole exome sequencing (15).

Sample preparation and next-generation sequencing of cfDNA

Whole blood was collected in K2 EDTA tubes; plasma and cellular components were separated by centrifugation at 800 × g for 10 minutes at 4°C. Plasma was centrifuged a second time at 18,000 × g at room temperature to remove any remaining cellular debris and stored at −80°C until the time of DNA extraction. DNA was isolated from plasma using the Qiagen Circulating Nucleic Acids Kit (Qiagen GmbH). TEC-Seq next-generation sequencing cell-free DNA libraries were prepared from 12 to 125 ng of cfDNA. Genomic libraries were prepared as previously described and targeted capture was performed using the Agilent SureSelect reagents and a custom set of hybridization probes targeting 58 genes, described in Supplementary Table S4 (14). TEC-Seq libraries were sequenced using 100 bp paired end runs on the Illumina HiSeq 2500 (Illumina). The analytical performance and validation including sensitivity and specificity and limits of detection of our ctDNA platform have been recently reported (14).

Primary processing of cfDNA next-generation sequencing data and identification of putative somatic mutations

Primary processing of next-generation sequence data for cfDNA samples was performed as previously described (14) using Illumina CASAVA software (v1.8), including demultiplexing and masking of dual index adapter sequences. Sequence reads were aligned against the human reference genome (hg19) using Novoalign with additional realignment of select regions using the Needleman–Wunsch method (17). Next, candidate somatic mutations, consisting of point mutations, small insertions, and deletions were identified using VariantDx (17) across the targeted regions of interest. VariantDx examined sequence alignments of cfDNA plasma samples while applying filters to exclude alignment and sequencing artifacts as previously described (14). Specifically, an alignment filter was applied to exclude quality failed reads, unpaired reads, and poorly mapped reads in the plasma. A base quality filter was applied to limit inclusion of bases with reported Phred quality score >30. Criteria for calling alterations in cfDNA have been previously described (14). TEC-Seq characteristics are shown in Supplementary Table S5.

Identification of tumor-derived cfDNA

Genomic alterations in ctDNA were cross-referenced against each patient's tumor-specific genomic alterations to identify bona fide tumor-specific ctDNA variants. Variants identified in ctDNA as previously described (14) and in the matching tumor with a MAF of ≥2% were considered tumor specific. We focused on somatic variants that were identified both in the tumor sample and in ctDNA for each patient to exclude variants related to clonal hematopoiesis (Supplementary Table S6). Our dataset is deposited in the database of Genotypes and Phenotypes (dbGaP; study ID 32485).

Clonality estimates of cfDNA variants

To assess the cellular prevalence of plasma mutations in their corresponding tumors, tumor samples of each case were analyzed as follows. The density of reads mapping to target and off-target regions in tumor whole exome sequence data was corrected for GC content, target size, and sequence complexity and compared with a reference panel of normal samples to establish log copy ratio values as a measure of relative copy number across the genome (18). Bin-level copy ratio values were segmented using circular binary segmentation (19). Segment copy ratio values and minor allele frequency of germline heterozygous SNPs overlapping the segments were analyzed to determine the purity and ploidy of the sample, and allele-specific copy number for segments using an in-house pipeline. Next, we used SCHISM (20) to determine the cellular prevalence of mutations based on the observed variant allele frequency, estimated copy number, and sample purity by following an approach similar to that previously described (13). This approach to clonality assessment was not feasible for mutations in four cases (CGLU168, CGLU206, CGLU219, and CGLU249) where purity and ploidy could not be determined due to low tumor content. Mutation cellularity analysis is summarized in Supplementary Table S7.

TCR sequencing and differential expansion analyses

TCR clones were evaluated in pretreatment tumor tissue (with the exception of CGLU117, where tumor tissue from the time of resistance was also analyzed), and 40 serial peripheral blood lymphocytes (PBL) by next-generation sequencing in the metastatic NSCLC cohort (Supplementary Table S8). DNA from pretreatment tumor samples and PBLs was isolated by using the Qiagen DNA FFPE and Qiagen DNA Blood Mini Kit, respectively (Qiagen). TCR-β CDR3 regions were amplified using the survey (tumor) or deep (PBLs) ImmunoSeq assay in a multiplex PCR method using 45 forward primers specific to TCR Vβ gene segments and 13 reverse primers specific to TCR Jβ gene segments (Adaptive Biotechnologies; refs. 21, 22). Productive TCR sequences were further analyzed. TCR sequencing data from TILs were used to identify tumor-specific TCR clonotypes in the peripheral blood. Peripheral TCR clones achieving a frequency of at least 0.005% were evaluated for differential abundance between baseline and the time of radiographic response using Fisher exact test with FDR P-value correction (corrected P ≤ 0.05). Those differentially abundant clones also found in the tumor were further selected to determine their frequencies in peripheral blood prior to treatment, at the time of response and upon emergence of resistance (Supplementary Tables S9 and S18). We calculated the average productive frequency of differentially abundant clones and used it as a metric of TCR dynamics during therapy. To cluster significantly expanded intratumoral TCR-β CDR3s based on potential recognition specificity, we used the GLIPH method (Grouping of Lymphocyte Interactions by Paratope Hotspots; ref. 23).

CDR3 and VJ gene usage analyses

Subsequent to initial filtering, we further reduced noise by eliminating clones that did not have frequencies beyond a mean rate of five counts. Thus, when two points were examined (baseline and time of radiologic response), the total sum of counts were greater or equal to 10. Using these data, we examined the usage of CDR3b Variable (V) and Joining (J) regions, and their overall clonal composition by known significant clones at the two time points.

Multiplex cytokine immunoassay

We used a multiplex bead-based immunoassay on the Luminex platform that examines cytokines involved in T-cell activation, expansion, differentiation, and long-term proliferation (IFNγ, IL1, IL2), Th1 immune response (IL12), acquisition of the Th2 phenotype (IL4) as well as immunosuppressive cytokines important for regulatory T cells (IL10) in four patients of the early stage cohort where additional serum was available. Differences in concentration of cytokines were evaluated between baseline and on treatment (week 2–6) samples.

Statistical analyses

ctDNA values were dichotomized as detectable and undetectable. Characteristics for each group were compared using chi-square or Fischer exact test for categorical variables. Pearson correlation coefficient (R) was used to assess correlations between continuous variables. Differences between molecular responders and nonresponders were assessed by the Mann–Whitney test. Tumors were classified based on their nonsynonymous sequence alteration load in high and low mutators as previously described (24). The median point estimate and 95% confidence interval for PFS and OS were estimated by the Kaplan–Meier method. Survival curves were compared by using the log-rank test. Univariate Cox proportional hazards regression analysis was used to determine the impact of ctDNA molecular response on PFS and OS. All P values were based on two-sided testing and differences were considered significant at P < 0.05. Statistical analyses were done using the SPSS software program (version 25.0.0 for Windows; IBM).

Overall approach and patient characteristics

We analyzed 105 serial blood samples from 38 patients with NSCLC, including 24 patients with metastatic NSCLC during immune checkpoint blockade and 14 patients with stage I to IIIA surgically resectable NSCLC that received anti-PD1 therapy as part of a clinical trial of neoadjuvant nivolumab (Table 1; Supplementary Table S1; ref. 15). The median duration of follow-up was 12.7 months (range 3.0–37.8 months) and 16 months (range 2–30 months) for the metastatic and early-stage patients respectively, and median duration of treatment was 7 months (range 1–20 months) for the metastatic cohort. We evaluated response to immune checkpoint blockade using standard computed tomographic (CT) imaging and changes in tumor burden were assessed by RECIST 1.1. Blood samples for the metastatic patients with NSCLC were prospectively collected prior to therapy, at an early time point between 4 and 8 weeks from treatment initiation and at additional serial time points during therapy until the time of disease progression (Supplementary Tables S2 and S3). For the patients with early-stage NSCLC treated with anti-PD1 therapy in the neoadjuvant setting, blood samples were collected prior to immunotherapy, at 2 weeks, immediately prior to resection and post-resection (Supplementary Table S2). ctDNA was measured using the TEC-Seq approach (14) and the TCR repertoire was studied longitudinally by means of TCR sequencing (Fig. 1). Given the possibility of hematopoietic alterations that may be detected in the plasma (14), especially in heavily treated patients, we focused only on tumor-specific sequence alterations in cell-free DNA. Clinical characteristics, outcome, and liquid biopsy analyses are summarized in Table 1.

Table 1.

Summary of clinical and molecular characteristics

PatientImmune-targeted therapyMetastatic vs. early stageProgression-free survival (0: progression-free; 1: progression)PFS (months)Overall survival (0: censored; 1: dead of disease)OS (months)Number of tumor-specific variants at baselinectDNA molecular response (0 = no, 1 = yes)Clonal TCR expansion (0 = no, 1 = yes)Tumor mutation burden
Metastatic NSCLC (N = 24) 
CGLU111 Nivolumab IV N/A 31.5 174 
 CGLU115 Nivolumab IV 3.2 3.8 285 
 CGLU121 Nivolumab IV 1.3 9.4 68 
 CGLU159 Nivolumab IV 3.9 5.6 65 
 CGLU160 Nivolumab IV 13.6 13.6 N/A 50 
 CGLU161 Nivolumab–ipilimumab IV 8.6 13.2 ND 127 
 CGLU162 Nivolumab IV 7.0 7.0 169 
 CGLU168 Nivolumab IV 7.3 12.6 N/A 411 
 CGLU203 Nivolumab IV 3.9 3.9 90 
 CGLU211 Nivolumab IV 10.7 21.3 N/A 161 
 CGLU212 Nivolumab IV 12.3 12.8 368 
 CGLU135 Nivolumab IV 23.8 38.7 358 
 CGLU127 Nivolumab IV 9.9 25.4 N/A N/A 335 
 CGLU117 Nivolumab IV 7.8 13.8 296 
 CGLU243 Nivolumab IV 2.4 11.4 42 
 CGLU329 Pembrolizumab IV 14.0 14.0 ND N/A 91 
 CGLU337 Pembrolizumab-chemotherapy IV 14.0 14.0 N/A 846 
 CGLU340 nivo-anti-LAG3 IV 6.6 13.4 1* N/A N/Aa 
 CGLU341 Pembrolizumab-chemotherapy IV 13.0 13.0 ND N/A 624 
 CGLU347 Pembrolizumab IV 5.0 12.0 1* N/A N/Aa 
 CGLU348 Pembrolizumab IV 3.0 3.0 N/A 191 
 CGLU357 Pembrolizumab IV 5.2 13.0 1* N/A N/Aa 
 CGLU351 Pembrolizumab IV 12.0 12.0 1* N/A N/Aa 
 CGU368 Pembrolizumab-chemotherapy IV 12.0 12.0 ND N/A N/Aa 
Early-stage NSCLC (N = 14) 
CGLU204 Nivolumab IIA 20 20.0 ND N/A N/Ab 
 CGLU205 Nivolumab IIIA 30 30 N/A 99 
 CGLU206 Nivolumab IB 23 23 N/A N/Ab 
 CGLU215 Nivolumab IA ND N/A 310 
 CGLU217 Nivolumab IIIA 14 14 ND N/A 68 
 CGLU218 Nivolumab IB 17 17 ND N/A 
 CGLU219 Nivolumab IIIA N/A N/Ab 
 CGLU220 Nivolumab IIIA 17 17 ND N/A 26 
 CGLU221 Nivolumab IIA N/A 28 N/A 190 
 CGLU222 Nivolumab IIA 28 N/A 75 
 CGLU224 Nivolumab IB 11 11 ND N/A 105 
 CGLU225 Nivolumab IIIA 15 15 N/A N/Ab 
 CGLU249 Nivolumab IIB N/A N/Ab 
 CGLU279 Nivolumab IIB 12 12 ND N/A N/Ab 
PatientImmune-targeted therapyMetastatic vs. early stageProgression-free survival (0: progression-free; 1: progression)PFS (months)Overall survival (0: censored; 1: dead of disease)OS (months)Number of tumor-specific variants at baselinectDNA molecular response (0 = no, 1 = yes)Clonal TCR expansion (0 = no, 1 = yes)Tumor mutation burden
Metastatic NSCLC (N = 24) 
CGLU111 Nivolumab IV N/A 31.5 174 
 CGLU115 Nivolumab IV 3.2 3.8 285 
 CGLU121 Nivolumab IV 1.3 9.4 68 
 CGLU159 Nivolumab IV 3.9 5.6 65 
 CGLU160 Nivolumab IV 13.6 13.6 N/A 50 
 CGLU161 Nivolumab–ipilimumab IV 8.6 13.2 ND 127 
 CGLU162 Nivolumab IV 7.0 7.0 169 
 CGLU168 Nivolumab IV 7.3 12.6 N/A 411 
 CGLU203 Nivolumab IV 3.9 3.9 90 
 CGLU211 Nivolumab IV 10.7 21.3 N/A 161 
 CGLU212 Nivolumab IV 12.3 12.8 368 
 CGLU135 Nivolumab IV 23.8 38.7 358 
 CGLU127 Nivolumab IV 9.9 25.4 N/A N/A 335 
 CGLU117 Nivolumab IV 7.8 13.8 296 
 CGLU243 Nivolumab IV 2.4 11.4 42 
 CGLU329 Pembrolizumab IV 14.0 14.0 ND N/A 91 
 CGLU337 Pembrolizumab-chemotherapy IV 14.0 14.0 N/A 846 
 CGLU340 nivo-anti-LAG3 IV 6.6 13.4 1* N/A N/Aa 
 CGLU341 Pembrolizumab-chemotherapy IV 13.0 13.0 ND N/A 624 
 CGLU347 Pembrolizumab IV 5.0 12.0 1* N/A N/Aa 
 CGLU348 Pembrolizumab IV 3.0 3.0 N/A 191 
 CGLU357 Pembrolizumab IV 5.2 13.0 1* N/A N/Aa 
 CGLU351 Pembrolizumab IV 12.0 12.0 1* N/A N/Aa 
 CGU368 Pembrolizumab-chemotherapy IV 12.0 12.0 ND N/A N/Aa 
Early-stage NSCLC (N = 14) 
CGLU204 Nivolumab IIA 20 20.0 ND N/A N/Ab 
 CGLU205 Nivolumab IIIA 30 30 N/A 99 
 CGLU206 Nivolumab IB 23 23 N/A N/Ab 
 CGLU215 Nivolumab IA ND N/A 310 
 CGLU217 Nivolumab IIIA 14 14 ND N/A 68 
 CGLU218 Nivolumab IB 17 17 ND N/A 
 CGLU219 Nivolumab IIIA N/A N/Ab 
 CGLU220 Nivolumab IIIA 17 17 ND N/A 26 
 CGLU221 Nivolumab IIA N/A 28 N/A 190 
 CGLU222 Nivolumab IIA 28 N/A 75 
 CGLU224 Nivolumab IB 11 11 ND N/A 105 
 CGLU225 Nivolumab IIIA 15 15 N/A N/Ab 
 CGLU249 Nivolumab IIB N/A N/Ab 
 CGLU279 Nivolumab IIB 12 12 ND N/A N/Ab 

NOTE: Major pathologic response (MPR) was defined as ≤10% viable tumor cells at the time of surgical resection (15). TCR expansion was assessed at the time of radiographic response.

Abbreviations: PFS, progression-free survival; OS, overall survival; ND, not detected; N/A, not evaluable.

aExome sequencing was not performed; for these cases, CLIA-targeted NGS was performed for clinical purposes.

bThe baseline tumor was not available for whole exome sequencing; the resection sample was analyzed and used to identify tumor-specific variants in ctDNA.

Figure 1.

Overview of next-generation sequencing and T-cell analyses. We used serial blood samples collected at baseline, early after treatment initiation and at additional timepoints during immune checkpoint blockade to determine ctDNA and TCR repertoire dynamics. ctDNA trends were evaluated by TEC-Seq and the evolving TCR repertoire was assessed by TCR next-generation sequencing. Dynamic changes in ctDNA and TCR clonotypic expansions were used to identify molecular response patterns and compared with RECIST 1.1 tumor burden evaluations. T0–T4 denote serial timepoints from the time of treatment initiation (T0) to the time of molecular response (T1), radiologic response (T2), molecular resistance (T3), and radiologic progression (T4).

Figure 1.

Overview of next-generation sequencing and T-cell analyses. We used serial blood samples collected at baseline, early after treatment initiation and at additional timepoints during immune checkpoint blockade to determine ctDNA and TCR repertoire dynamics. ctDNA trends were evaluated by TEC-Seq and the evolving TCR repertoire was assessed by TCR next-generation sequencing. Dynamic changes in ctDNA and TCR clonotypic expansions were used to identify molecular response patterns and compared with RECIST 1.1 tumor burden evaluations. T0–T4 denote serial timepoints from the time of treatment initiation (T0) to the time of molecular response (T1), radiologic response (T2), molecular resistance (T3), and radiologic progression (T4).

Close modal

ctDNA dynamics and tumor response

In the metastatic NSCLC cohort, ctDNA was detected in 19 of 24 patients either at baseline (n = 14) or at other time points when baseline samples were not available (n = 5), with a median mutant allele fraction of 1.87% (range 0.09%–34.7%). In the early-stage cohort, ctDNA was detected at baseline in 7 of 14 patients, with a median allele fraction of 0.34% (range 0.15%–2.19%). For patients with detectable ctDNA, an average of one tumor-specific alterations were detected (median 1, range 1–4) affecting one or more of 12 driver genes, including those commonly altered in lung cancer (Supplementary Tables S4–S6). The vast majority of tumor-specific variants were clonal in the corresponding tumor samples (Supplementary Table S7).

We observed three patterns of molecular response in ctDNA for patients treated with immune checkpoint inhibitors. Among the patients with a molecular response (n = 9), individuals had a dramatic reduction in ctDNA to undetectable levels on average at 9 weeks from treatment initiation (Fig. 2A–E; Supplementary Fig. S1). As an example, for patient CGLU111 with a sustained clinical response, ctDNA-based molecular analyses showed a complete molecular response at week 4, more than 5 weeks prior to a radiologic partial response and 26 weeks earlier than complete radiologic response determined by RECIST 1.1 (Fig. 2). In contrast, for patients with a pattern of molecular resistance (n = 10), ctDNA levels had limited fluctuations or displayed a rise 3 to 16 weeks after therapeutic initiation. As a representative patient, ctDNA levels in CGLU121 continued to rise from the time of initiation of immune checkpoint blockade, consistent with radiographic disease progression (Fig. 3A–E). All patients with ctDNA features of primary molecular resistance had radiologic disease progression that followed molecular resistance by 5.5 weeks (Supplementary Fig. S2).

Figure 2.

ctDNA and TCR clonal dynamics for a patient with sustained response to anti-PD1. A, ctDNA (TP53 993+1G>T mutation, blue) decreased to undetectable levels, signifying a complete molecular response at week 4. In contrast, CT imaging did not accurately capture the rate (B) or timing (C) of tumor regression (RECIST tumor burden dynamics, green). A complete response by RECIST 1.1 was achieved 26 weeks later than the molecular response (C). In parallel, TCR repertoire dynamics revealed clonotypic amplifications of intratumoral TCR clones in peripheral blood at the time of radiographic response. TCR clones with statistically significant differential abundance were evaluated as individual clones (D) and as a composite of productive frequencies (E). Patient was off anti-PD1 therapy and on immunosuppressive therapy at week 30 (arrow) due to emergence of immune-related toxicity.

Figure 2.

ctDNA and TCR clonal dynamics for a patient with sustained response to anti-PD1. A, ctDNA (TP53 993+1G>T mutation, blue) decreased to undetectable levels, signifying a complete molecular response at week 4. In contrast, CT imaging did not accurately capture the rate (B) or timing (C) of tumor regression (RECIST tumor burden dynamics, green). A complete response by RECIST 1.1 was achieved 26 weeks later than the molecular response (C). In parallel, TCR repertoire dynamics revealed clonotypic amplifications of intratumoral TCR clones in peripheral blood at the time of radiographic response. TCR clones with statistically significant differential abundance were evaluated as individual clones (D) and as a composite of productive frequencies (E). Patient was off anti-PD1 therapy and on immunosuppressive therapy at week 30 (arrow) due to emergence of immune-related toxicity.

Close modal
Figure 3.

ctDNA and TCR clonal dynamics for a patient with primary resistance to anti-PD1. A, ctDNA levels [EGFR 745KELREA>T(blue) and TP53173V>L mutations (red)] continued to rise from the time of initiation of anti-PD1 therapy. For this patient, the change in the RECIST tumor burden was similar to the increase in ctDNA levels (RECIST tumor burden dynamics, green; B); however, molecular resistance was detected earlier than conventional CT imaging (C). There were no clones with statistically significant expansion at week 4 compared with baseline; top 10 intratumoral clones found in peripheral blood are shown as individual clones (D) and by their average productive frequency (E).

Figure 3.

ctDNA and TCR clonal dynamics for a patient with primary resistance to anti-PD1. A, ctDNA levels [EGFR 745KELREA>T(blue) and TP53173V>L mutations (red)] continued to rise from the time of initiation of anti-PD1 therapy. For this patient, the change in the RECIST tumor burden was similar to the increase in ctDNA levels (RECIST tumor burden dynamics, green; B); however, molecular resistance was detected earlier than conventional CT imaging (C). There were no clones with statistically significant expansion at week 4 compared with baseline; top 10 intratumoral clones found in peripheral blood are shown as individual clones (D) and by their average productive frequency (E).

Close modal

The third observed pattern, seen in five of the molecular responders, was one consistent with molecular acquired resistance, where ctDNA dynamics reflected clonal evolution under selective pressure of anti-PD1 therapy and emergence of immune escape. In such cases, tumor-specific variants were undetectable at the time of response followed by increase in mutant allele fraction at the time of acquired resistance (Supplementary Fig. S1). Emergence of molecular resistance preceded disease progression on imaging by an average of 10.8 weeks. Overall, ctDNA-based molecular responses were detected on average 8.7 weeks earlier than conventional RECIST1.1 response assessment (6.7 weeks vs. 15.4 weeks, P = 0.004; Supplementary Fig. S3).

Early ctDNA clearance was a significant prognostic factor for PFS and OS. Patients with a reduction of ctDNA to undetectable levels demonstrated a significantly longer PFS and OS compared with patients with no evidence of ctDNA elimination (log rank P = 0.001 and 0.008, respectively; Fig. 4A and B; Supplementary Fig. S4). The duration of the molecular responses tightly correlated with PFS and OS (Supplementary Fig. S5).

Figure 4.

Early ctDNA clearance predicts PFS and OS. Patients with reduction of ctDNA to undetectable levels demonstrated a significantly longer PFS and OS compared with patients with no evidence of ctDNA elimination (log-rank P = 0.001 and 0.008, respectively; A and B). Patients with undetectable ctDNA (molecular responders) clustered together independent of their tumor mutation burden (C) and the same pattern was observed for patients with detectable ctDNA (molecular nonresponders; D). Patients with ctDNA molecular response and either high or low tumor mutation burden had a significantly longer PFS and OS (log-rank P = 0.015 and 0.027, respectively).

Figure 4.

Early ctDNA clearance predicts PFS and OS. Patients with reduction of ctDNA to undetectable levels demonstrated a significantly longer PFS and OS compared with patients with no evidence of ctDNA elimination (log-rank P = 0.001 and 0.008, respectively; A and B). Patients with undetectable ctDNA (molecular responders) clustered together independent of their tumor mutation burden (C) and the same pattern was observed for patients with detectable ctDNA (molecular nonresponders; D). Patients with ctDNA molecular response and either high or low tumor mutation burden had a significantly longer PFS and OS (log-rank P = 0.015 and 0.027, respectively).

Close modal

Radiographic imaging at the time of first assessment was a worse predictor of outcome to anti-PD1 therapy compared with ctDNA molecular response for these patients (Fig. 4A and B; Supplementary Fig. S6). Patients with radiographically stable disease (n = 12) had differential responses to immune checkpoint blockade that were consistent with their molecular response pattern (Supplementary Fig. S7). More specifically, five patients with stable disease by imaging showed a clear molecular response pattern, with ctDNA elimination between week 4 and 13 from immune checkpoint blockade initiation (Supplementary Fig. S1). All five patients derived clinical benefit from PD-1 blockade (PFS and OS ranging from 7.3 to 13.6 and 12 to 21.3 months, respectively; Supplementary Fig. S7), suggesting that imaging failed to detect the magnitude of therapeutic response.

Interestingly, ctDNA molecular responses more accurately predicted PFS and OS compared with tumor mutation burden in our cohort (TMB; Supplementary Fig. S8). When TMB and ctDNA were combined the ctDNA-based molecular responders clustered together independent of the TMB for both PFS and OS (Fig. 4C and D). Given that clonal mutation burden may be a more accurate predictor of response to immune checkpoint blockade, we performed survival analyses incorporating clonal TMB and we again found that ctDNA dynamics predict survival independent of clonal TMB status (Supplementary Fig. S9).

Molecular responses predict pathologic response to immune checkpoint blockade

Given the challenges with radiologic response assessments to immune checkpoint blockade, we sought to validate our observations in a NSCLC cohort where the therapeutic effect was measured at a pathologic level instead of using conventional imaging. We analyzed serial plasma samples from a recently reported clinical trial of neoadjuvant nivolumab for early-stage operable NSCLC (15). For these patients, the therapeutic effect was rigorously measured by pathologic assessment of residual tumor after two doses of anti-PD1 therapy (15, 25). Similar to our initial analyses, we observed that all tumors with a major or partial pathologic response to anti-PD1 therapy demonstrated a molecular response pattern of elimination of tumor-specific mutations in the circulation (Fig. 5A; Supplementary Fig. S10). In contrast, tumors without a pathologic response demonstrated a molecular resistance pattern at the time of resection of the primary tumor (Fig. 5B; Supplementary Fig. S10).

Figure 5.

Early ctDNA clearance is associated with pathologic response to anti-PD1 therapy. Molecular responses were consistent with pathologic responses to anti-PD1 therapy in early-stage NSCLC. A, For a patient with a major pathologic response, ctDNA elimination (TP53 K132N mutation, blue) accurately captured the therapeutic effect compared with RECIST tumor burden dynamics (green) that showed stable disease. In contrast, ctDNA levels [KRAS G12C (blue) and ALK G875R mutations (purple)] increased from baseline for a patient that did not achieve a pathologic response to anti-PD1 therapy (B). Changes in RECIST tumor burden, shown on the secondary axis of each plot, did not accurately predict outcome as both patients were classified as stable disease. The timeline of anti-PD1 therapy dosing, radiographic assessments, and tumor resection is shown below each graph.

Figure 5.

Early ctDNA clearance is associated with pathologic response to anti-PD1 therapy. Molecular responses were consistent with pathologic responses to anti-PD1 therapy in early-stage NSCLC. A, For a patient with a major pathologic response, ctDNA elimination (TP53 K132N mutation, blue) accurately captured the therapeutic effect compared with RECIST tumor burden dynamics (green) that showed stable disease. In contrast, ctDNA levels [KRAS G12C (blue) and ALK G875R mutations (purple)] increased from baseline for a patient that did not achieve a pathologic response to anti-PD1 therapy (B). Changes in RECIST tumor burden, shown on the secondary axis of each plot, did not accurately predict outcome as both patients were classified as stable disease. The timeline of anti-PD1 therapy dosing, radiographic assessments, and tumor resection is shown below each graph.

Close modal

Peripheral TCR landscape and therapeutic outcome

We investigated how immune checkpoint blockade affects the peripheral TCR repertoire and whether there are TCR clonotype dynamic changes reflective of a systemic antitumor immune response. We focused our analyses on TCR clones found in the tumor microenvironment using TCR sequencing and investigated their dynamics in the peripheral blood, identifying those with a statistically significant differential abundance from baseline. Twelve of the 24 metastatic patients with NSCLC had available samples from both tumor infiltrating lymphocytes and peripheral blood lymphocytes for analysis (Supplementary Tables S3 and S8), including five that had previously undergone TCR sequencing (13) but had not been analyzed using this approach.

Similar to ctDNA analyses, we observed distinct patterns in TCR clonotype dynamics among the analyzed patients. For patients with clinical responses to immune checkpoint blockade, a statistically significant oligoclonal expansion of preexisting intratumoral T-cell clones was observed in peripheral blood at the time of radiologic response to PD1 blockade (CGLU111, CGLU117, CGLU127, and CGLU212; Fig. 2; Supplementary Fig. S11; Supplementary Tables S9–S12). For patients who developed acquired resistance, productive frequencies of intratumoral clones significantly decreased in peripheral blood at the time of acquired resistance (CGLU117, CGLU127, CGLU135, and CGLU161; Supplementary Fig. S11; Supplementary Tables S10, S11, S13, S14), with a timing that was similar to ctDNA analyses for most cases.

In contrast, for patients CGLU121 and CGLU115 that had primary resistance to immunotherapy, we did not identify any differentially abundant TCR clones among serial peripheral blood samples (Fig. 3; Supplementary Fig. S12). These patients progressed radiographically within 5 to 13 weeks from initiation of therapy and, in line with the clinical course, there was no evidence of TCR clonal expansion among the intratumoral TCR repertoire. A transient oligoclonal TCR expansion was observed for nonresponding patient CGLU159 at week 11; however, productive frequencies of differentially abundant clones quickly decreased to baseline levels at week 16, which coincided with disease progression (Supplementary Fig. S12; Supplementary Table S15). Patients CGLU162, CGLU203, and CGLU243 had 1 to 26 intratumoral TCR clones with differential abundance at the time of best radiographic response compared with baseline but were classified as ctDNA molecular nonresponders (Supplementary Fig. S12; Supplementary Tables S16–S18). Patients CGLU203 and CGLU243 had unfavorable outcome to anti-PD1 therapy, suggesting that for these patients, ctDNA kinetics may more accurately predict therapeutic outcome.

We did not identify any shared TCR clones among the differentially expanded ones for all patients analyzed, consistent with the notion that the mutation-associated neoantigen repertoires are largely private. We evaluated putative shared CDR3 motifs among significantly expanded TCR clones using the grouping of lymphocyte interactions by paratope hotspots algorithm (23). Interestingly, TCR clones CSARVGVGNTIYF and CSARSGVGNTIYF, which were differentially abundant at the time of response to immune checkpoint blockade for patient CGLU127 and CGLU135, respectively, clustered together, suggesting a common specificity to a tumor- or mutation-associated antigen. We subsequently investigated potential differential sequence features focusing on Variable (V) and Joining (J) gene usage and CDR3 lengths among different timepoints for each patient. Usage of specific V and J gene segments increased at the time of response compared with baseline for a patient with sustained response (CGLU111) in contrast to a representative patient with primary resistance (CGLU121; Supplementary Fig. S13). Our findings on differential V gene usage may suggest clonotypic amplifications of specific immune subsets (CD8+ vs. CD4+) during immune checkpoint blockade (26).

The unique nature of responses to immune checkpoint blockade (27, 28) and known limitations of conventional radiologic response assessments (3) highlight the need for development of biomarker-driven approaches to interpret therapeutic responses. Success of immunotherapy approaches depends on choosing patient populations most likely to benefit. There is therefore an urgent clinical need for molecular assays of response and resistance to immune-targeted agents. To this end, we analyzed ctDNA and TCR clonal dynamics during immune checkpoint blockade in NSCLC and assessed the value of longitudinal monitoring of liquid biopsies as a surrogate for response to therapy. Our findings indicate that ctDNA dynamics after treatment initiation may allow patients with primary resistance to immune checkpoint blockade to be rapidly identified and redirected to receive alterative options.

Noninvasive detection and monitoring of acquired resistance to EGFR-targeted therapy has been evaluated by serial sampling of ctDNA (29, 30) and as we have shown in a complementary study, changes in ctDNA levels may predict response to targeted therapy in NSCLC (8). Longitudinal assessment of ctDNA in patients with metastatic melanoma receiving anti-PD1 therapy has been demonstrated to be an accurate predictor of tumor response and therapeutic outcome (31) and early ctDNA clearance may correlate with durable clinical benefit to PD-1 blockade (10, 12). ctDNA dynamics may be also informative in differentiating pseudoprogression from disease progression during immunotherapy (32). However, these approaches have been limited by low sensitivity and specificity of ctDNA methods, with many patients lacking detectable alterations and potential admixture between tumor alterations and those involved in clonal hematopoiesis (33, 34).

Moreover, interpretation of ctDNA analyses without knowledge of tumor-specific somatic alterations may be difficult in the setting of heavily pretreated patient populations such as patients with late-stage lung cancer, given the mutagenic effects of systemic chemotherapy and ionizing radiation on cells of the myeloid lineage (35). To address the possible presence of alterations in cfDNA from clonal nonmalignant hematopoietic cells, we have focused our analyses of variants in ctDNA that were also identified through next-generation sequencing of the matched tumor, allowing distinction of tumor-specific from blood cell proliferation variants.

Clonal expansion of intratumoral T cells may predict therapeutic outcome for immune checkpoint blockade (36); however, little is known about the significance of peripheral expansion of TCR clones found in the tumor microenvironment during therapy. Expansion of peripheral CD8+ T-cell populations has been shown to precede immune-related adverse events in patients treated with ipilimumab (37). We investigated whether the patients with NSCLC in our cohort developed immune-related adverse events at the time of TCR clonal expansion and did not identify a definitive pattern with the exception of patient CGLU243, where pneumonitis emerged shortly after treatment initiation. Although there were cases for which TCR expansion preceded the development of a grade 2 to 4 immune-related adverse event (CGLU161, CGLU117), such events were also noted significantly later from the time of TCR expansion (CGLU111, CGLU135). These observations highlight the challenges with interpretation of the evolving peripheral TCR repertoire. Assessing the quality of the immune response in conjunction with clonotypic amplifications may provide additional information on the evolving TCR repertoire; to this end, we looked at differences in cytokine levels in selected early-stage patients with available serum at baseline and 2 to 6 weeks during anti-PD1 therapy. We did not identify any significant changes in cytokine levels in peripheral blood between baseline and week 2 to 6 on anti-PD1 therapy, however these analyses were limited by small number of cases tested (Supplementary Fig. S14).

In summary, we have developed dynamic assays that capture the tumor-immune system equilibrium and assess immune editing of neoantigens during immunotherapy. We have shown that these approaches have advantages compared with conventional radiologic response assessment and static molecular analyses such as baseline TMB. We believe that these methods are especially suited for the interpretation of unique responses seen with immune-targeted agents that are not adequately captured by traditional response criteria. In addition to more accurately predicting long-term response to immunotherapy, we were able to predict therapeutic outcome on average 8.7 weeks earlier than radiographic imaging. However, our work is limited by the small sample size, cohort heterogeneity, and retrospective nature of the analyses. Validation of these findings may lead to early therapeutic decisions to ensure that an ineffective treatment is discontinued as well as allow response adaptive combination and sequencing of subsequent therapies. Additional work will be needed to address the frequency of serial monitoring and feasibility of interpreting ctDNA dynamics without prior knowledge of tumor mutations. Prospective studies will be needed to assess whether switching therapy based on ctDNA dynamics prior to radiologic progression will improve outcome and ultimately whether a liquid biopsy approach can replace conventional imaging as a gold standard for early response assessment to immune checkpoint blockade.

P.M. Forde reports receiving a commercial research grant from AstraZeneca, BMS, Corvus, Kyowa, Merck, and Novartis, and is a consultant/advisory board member of Abbvie, AstraZeneca, BMS, Boehringer, EMD Serono, Lilly, Merck, and Novartis. J.R. White is a consultant at Personal Genome Diagnostics and has ownership interest (including stock, patents, etc.) in Resphera Biosciences. J. Naidoo reports receiving a commercial research grant from AstraZeneca/MedImmune and Merck, is a consultant/advisory board member of AstraZeneca/MedImmune, Bristol-Myers Squibb, Takeda, Genentech/Roche, and has provided expert testimony for Bristol-Myers Squibb and AstraZeneca/MedImmune. V. Adleff is a consultant/advisory board member of Personal Genome Diagnostics. C.L. Hann reports receiving other commercial research support from Bristol-Myers Squibb and is a consultant/advisory board member of Bristol-Myers Squibb and Genentech. B. Levy is a consultant/advisory board member of Genentech, AstraZeneca, Eli Lilly, Celgene, and Takeda. J. Feliciano is a consultant/advisory board member of AstraZeneca, Merck, Genentech, and Eli Lilly, and has provided expert testimony for legal review. P. Illei reports receiving other commercial research support from Bristol Myers-Squib and is a consultant/advisory board member of AstraZeneca, Bayer, Roche, and Abbvie. J.M. Isbell has ownership interest (including stock, patents, etc.) in LumaCyte, Inc. J.L. Sauter has ownership interest (including stock, patents, etc.) in Merck & Company, Inc. (New), Chemed Corporation, Pfizer, Thermo Fisher Scientific, Celgene, and Allergan PLC SHS. J. Taube reports receiving a commercial research grant from Bristol Myers Squibb and is a consultant/advisory board member of Bristol Myers Squibb, Astra Zeneca, Merck, and Amgen. D.M. Pardoll reports receiving a commercial research grant from Bristol Myers Squibb, Astra Zeneca, Compugen, Merck, has ownership interest (including stock, patents, etc.) in WindMil, Trieza Therapeutics, Potenza, Ervaxx, Dracen Pharmaceuticals, DNAtrix, Aduro Biotech, Five Prime, and Tizona, and is a consultant/advisory board member of Amgen, Bayer, Camden Partners, Jansesen, FLX Bio, Immunomics, and Rock Springs Capitol. J.E. Chaft is a consultant/advisory board member of AstraZeneca, Merck, BMS, and Genentech. M.D. Hellmann reports receiving a commercial research grant from BMS, is a consultant/advisory board member of Merck, BMS, AstraZeneca/MedImmune, Genentech/Roche, Janssen, Nektar, Syndax, Mirati, and Shattuck Labs, and has provided expert testimony for a patent that has been filed by MSK related to the use of tumor mutation burden to predict response to immunotherapy (PCT/US2015/062208), which has received licensing fees from PGDx. J.R. Brahmer reports receiving a commercial research grant from BMS and is a consultant/advisory board member of BMS, Merck, Genentech, Amgen, Janssen, Syndax, Celgene, and AstraZeneca. V.E. Velculescu has ownership interest (including stock, patents, etc.) in Ignyta and Personal Genome Diagnostics and is a consultant/advisory board member of Ignyta and Personal Genome Diagnostics. No potential conflicts of interest were disclosed by the other authors.

Conception and design: V. Anagnostou, P.M. Forde, J.E. Chaft, J.R. Brahmer, V.E. Velculescu

Development of methodology: V. Anagnostou, P.M. Forde, N. Niknafs, J. Phallen, V. Adleff, D.M. Pardoll, V.E. Velculescu

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): V. Anagnostou, P.M. Forde, C. Hruban, J. Naidoo, S. Rosner, J. Phallen, A. Leal, K.N. Smith, L. Rhymee, C.L. Hann, J. Feliciano, C. Georgiades, P. Illei, E. Gabrielson, M.V. Brock, J.M. Isbell, J.L. Sauter, J. Taube, M.D. Hellmann, J.R. Brahmer

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): V. Anagnostou, P.M. Forde, J.R. White, N. Niknafs, C. Hruban, K. Marrone, I.K.A. Sivakumar, D.C. Bruhm, J. Phallen, V. Adleff, K.N. Smith, T.R. Cottrell, D.N. Palsgrove, B. Levy, J. Feliciano, F. Verde, J.M. Isbell, J.L. Sauter, J. Taube, R.B. Scharpf, R. Karchin, D.M. Pardoll, M.D. Hellmann, J.R. Brahmer, V.E. Velculescu

Writing, review, and/or revision of the manuscript: V. Anagnostou, P.M. Forde, J.R. White, N. Niknafs, J. Naidoo, K. Marrone, I.K.A. Sivakumar, D.C. Bruhm, K.N. Smith, T.R. Cottrell, C.L. Hann, B. Levy, J. Feliciano, P. Illei, M.V. Brock, J.M. Isbell, J.L. Sauter, J. Taube, D.M. Pardoll, M.D. Hellmann, J.R. Brahmer, V.E. Velculescu

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J.R. White, S. Rosner, L. Rhymee, Q.K. Li, V.E. Velculescu

Study supervision: V. Anagnostou, P.M. Forde, R.B. Scharpf, R. Karchin, V.E. Velculescu

Other (biopsy performance, biopsy consent): J. Feliciano

This work was supported in part by U.S. NIH grants CA121113 (to V. Velculescu, V. Anagnostou), CA006973 (to D. Pardoll, V. Velculescu), CA180950 (to V. Velculescu), the Commonwealth Foundation (to V. Velculescu), the Bloomberg-Kimmel Institute for Cancer Immunotherapy (to V. Anagnostou, P. Forde, J. Brahmer, D. Pardoll, V. Velculescu), the Dr. Miriam and Sheldon G. Adelson Medical Research Foundation (to V. Velculescu), the Eastern Cooperative Oncology Group- American College of Radiology Imaging Network (to V. Anagnostou), MacMillan Foundation (to V. Anagnostou), the V Foundation (to V. Anagnostou, V. Velculescu), the ICTR-ATIP UL1TR001079 (to V. Anagnostou), the Pardee Foundation (to V. Anagnostou), Swim Across America (to V. Anagnostou), the William R. Brody Faculty Scholarship (to R. Karchin), the SU2C-ACS Lung Cancer Dream Team (to P. Forde and E. Gabrielson), PRIME Oncology (to J. Naidoo), the MSK Cancer Center Support Grant/Core Grant (P30 CA008747), the SU2C DCS International Translational Cancer Research Dream Team Grant (SU2C-AACR-DT1415; to V. Velculescu), the SU2C-LUNGevity-American Lung Association Lung Cancer Interception Dream Team, Translational Cancer Research Grant (SU2C-AACR-DT23-17 to J. Brahmer, V. Velculescu), the Allegheny Health Network – Johns Hopkins Research Fund (to V. Anagnostou, V. Velculescu), the LUNGevity Foundation (to V. Anagnostou and P. Forde), the Mark Foundation (to A. Leal, V. Velculescu), and Bristol Meyers Squibb (to P. Forde). Stand Up To Cancer is a program of the Entertainment Industry Foundation administered by the American Association for Cancer Research. This publication was made possible in part by the Johns Hopkins Institute for Clinical and Translational Research (ICTR), which is funded in part by Grant Number UL1TR001079 from the National Center for Advancing Translational Sciences (NCATS), a component of the National Institutes of Health (NIH), and NIH Roadmap for Medical Research. Its contents are solely the responsibility of the authors and do not necessarily represent the official view of the Johns Hopkins ICTR, NCATS, or NIH. We thank Dr. Suzanne Topalian and members of our laboratories for helpful discussions and critical review of the manuscript.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Horn
L
,
Spigel
DR
,
Vokes
EE
,
Holgado
E
,
Ready
N
,
Steins
M
, et al
Nivolumab versus docetaxel in previously treated patients with advanced non-small-cell lung cancer: two-year outcomes from two randomized, open-label, phase III trials (CheckMate 017 and CheckMate 057)
.
J Clin Oncol
2017
;
35
:
3924
33
.
2.
Sharma
P
,
Allison
JP
. 
The future of immune checkpoint therapy
.
Science
2015
;
348
:
56
61
.
3.
Anagnostou
V
,
Yarchoan
M
,
Hansen
AR
,
Wang
H
,
Verde
F
,
Sharon
E
, et al
Immuno-oncology trial endpoints: capturing clinically meaningful activity
.
Clin Cancer Res
2017
;
23
:
4959
69
.
4.
Hodi
FS
,
Hwu
WJ
,
Kefford
R
,
Weber
JS
,
Daud
A
,
Hamid
O
, et al
Evaluation of immune-related response criteria and RECIST v1.1 in patients with advanced melanoma treated with pembrolizumab
.
J Clin Oncol
2016
;
34
:
1510
7
.
5.
Bohnsack
O
,
Hoos
A
,
Ludajic
K
. 
Adaptation of the immune related response criteria: irRECIST
.
Ann Oncol
2014
;
Supplement 4
:
iv361
iv72
.
6.
Abbosh
C
,
Birkbak
NJ
,
Wilson
GA
,
Jamal-Hanjani
M
,
Constantin
T
,
Salari
R
, et al
Phylogenetic ctDNA analysis depicts early-stage lung cancer evolution
.
Nature
2017
;
545
:
446
51
.
7.
Jamal-Hanjani
M
,
Wilson
GA
,
McGranahan
N
,
Birkbak
NJ
,
Watkins
TBK
,
Veeriah
S
, et al
Tracking the evolution of non-small-cell lung cancer
.
N Engl J Med
2017
;
376
:
2109
21
.
8.
Phallen
J
,
Leal
A
,
Woodward
B
,
Forde
PM
,
Naidoo
J
,
Marrone
K
, et al
Early noninvasive detection of response to targeted therapy in non-small cell lung cancer
.
Cancer Res
2019
;
79
:
1204
13
.
9.
Lipson
EJ
,
Velculescu
VE
,
Pritchard
TS
,
Sausen
M
,
Pardoll
DM
,
Topalian
SL
, et al
Circulating tumor DNA analysis as a real-time method for monitoring tumor burden in melanoma patients undergoing treatment with immune checkpoint blockade
.
J Immunother Cancer
2014
;
2
:
42
.
doi: 10.1186/s40425-014-0042-0
.
10.
Cabel
L
,
Riva
F
,
Servois
V
,
Livartowski
A
,
Daniel
C
,
Rampanou
A
, et al
Circulating tumor DNA changes for early monitoring of anti-PD1 immunotherapy: a proof-of-concept study
.
Ann Oncol
2017
;
28
:
1996
2001
.
11.
Iijima
Y
,
Hirotsu
Y
,
Amemiya
K
,
Ooka
Y
,
Mochizuki
H
,
Oyama
T
, et al
Very early response of circulating tumour-derived DNA in plasma predicts efficacy of nivolumab treatment in patients with non-small cell lung cancer
.
Eur J Cancer
2017
;
86
:
349
57
.
12.
Goldberg
SB
,
Narayan
A
,
Kole
AJ
,
Decker
RH
,
Teysir
J
,
Carriero
NJ
, et al
Early assessment of lung cancer immunotherapy response via circulating tumor DNA
.
Clin Cancer Res
2018
;
24
:
1872
80
.
13.
Anagnostou
V
,
Smith
KN
,
Forde
PM
,
Niknafs
N
,
Bhattacharya
R
,
White
J
, et al
Evolution of neoantigen landscape during immune checkpoint blockade in non-small cell lung cancer
.
Cancer Discov
2017
;
7
:
264
76
.
14.
Phallen
J
,
Sausen
M
,
Adleff
V
,
Leal
A
,
Hruban
C
,
White
J
, et al
Direct detection of early-stage cancers using circulating tumor DNA
.
Sci Transl Med
2017
;
9
.
doi: 10.1126/scitranslmed.aan2415.
15.
Forde
PM
,
Chaft
JE
,
Smith
KN
,
Anagnostou
V
,
Cottrell
TR
,
Hellmann
MD
, et al
Neoadjuvant PD-1 blockade in resectable lung cancer
.
N Engl J Med
2018
;
378
:
1976
86
.
16.
Eisenhauer
EA
,
Therasse
P
,
Bogaerts
J
,
Schwartz
LH
,
Sargent
D
,
Ford
R
, et al
New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1)
.
Eur J Cancer
2009
;
45
:
228
47
.
17.
Jones
S
,
Anagnostou
V
,
Lytle
K
,
Parpart-Li
S
,
Nesselbush
M
,
Riley
DR
, et al
Personalized genomic analyses for cancer mutation discovery and interpretation
.
Sci Transl Med
2015
;
7
:
283ra53
.
18.
Talevich
E
,
Shain
AH
,
Botton
T
,
Bastian
BC
. 
CNVkit: genome-wide copy number detection and visualization from targeted DNA sequencing
.
PLoS Comput Biol
2016
;
12
:
e1004873
.
19.
Olshen
AB
,
Venkatraman
ES
,
Lucito
R
,
Wigler
M
. 
Circular binary segmentation for the analysis of array-based DNA copy number data
.
Biostatistics
2004
;
5
:
557
72
.
20.
Niknafs
N
,
Beleva-Guthrie
V
,
Naiman
DQ
,
Karchin
R
. 
SubClonal hierarchy inference from somatic mutations: automatic reconstruction of cancer evolutionary trees from multi-region next generation sequencing
.
PLoS Comput Biol
2015
;
11
:
e1004416
.
21.
Carlson
CS
,
Emerson
RO
,
Sherwood
AM
,
Desmarais
C
,
Chung
MW
,
Parsons
JM
, et al
Using synthetic templates to design an unbiased multiplex PCR assay
.
Nat Commun
2013
;
4
:
2680
.
doi: 10.1038/ncomms3680.
22.
Robins
HS
,
Campregher
PV
,
Srivastava
SK
,
Wacher
A
,
Turtle
CJ
,
Kahsai
O
, et al
Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells
.
Blood
2009
;
114
:
4099
107
.
23.
Glanville
J
,
Huang
H
,
Nau
A
,
Hatton
O
,
Wagar
LE
,
Rubelt
F
, et al
Identifying specificity groups in the T cell receptor repertoire
.
Nature
2017
;
547
:
94
8
.
24.
Wood
DE
,
White
JR
,
Georgiadis
A
,
Van Emburgh
B
,
Parpart-Li
S
,
Mitchell
J
, et al
A machine learning approach for somatic mutation discovery
.
Sci Transl Med
2018
;
10
.
doi: 10.1126/scitranslmed.aar7939.
25.
Cottrell
TR
,
Thompson
ED
,
Forde
PM
,
Stein
JE
,
Duffield
AS
,
Anagnostou
V
, et al
Pathologic features of response to neoadjuvant anti-PD-1 in resected non-small-cell lung carcinoma: a proposal for quantitative immune-related pathologic response criteria (irPRC)
.
Ann Oncol
2018
;
29
:
1853
60
.
26.
Riaz
N
,
Havel
JJ
,
Makarov
V
,
Desrichard
A
,
Urba
WJ
,
Sims
JS
, et al
Tumor and microenvironment evolution during immunotherapy with nivolumab
.
Cell
2017
;
171
:
934
49
e15
.
27.
Topalian
SL
,
Drake
CG
,
Pardoll
DM
. 
Targeting the PD-1/B7-H1(PD-L1) pathway to activate anti-tumor immunity
.
Curr Opin Immunol
2012
;
24
:
207
12
.
28.
Pardoll
DM
. 
The blockade of immune checkpoints in cancer immunotherapy
.
Nat Rev
2012
;
12
:
252
64
.
29.
Murtaza
M
,
Dawson
SJ
,
Tsui
DW
,
Gale
D
,
Forshew
T
,
Piskorz
AM
, et al
Non-invasive analysis of acquired resistance to cancer therapy by sequencing of plasma DNA
.
Nature
2013
;
497
:
108
12
.
30.
Oxnard
GR
,
Paweletz
CP
,
Kuang
Y
,
Mach
SL
,
O'Connell
A
,
Messineo
MM
, et al
Noninvasive detection of response and resistance in EGFR-mutant lung cancer using quantitative next-generation genotyping of cell-free plasma DNA
.
Clin Cancer Res
2014
;
20
:
1698
705
.
31.
Lee
JH
,
Long
GV
,
Boyd
S
,
Lo
S
,
Menzies
AM
,
Tembe
V
, et al
Circulating tumour DNA predicts response to anti-PD1 antibodies in metastatic melanoma
.
Ann Oncol
2017
;
28
:
1130
6
.
32.
Guibert
N
,
Mazieres
J
,
Delaunay
M
,
Casanova
A
,
Farella
M
,
Keller
L
, et al
Monitoring of KRAS-mutated ctDNA to discriminate pseudo-progression from true progression during anti-PD-1 treatment of lung adenocarcinoma
.
Oncotarget
2017
;
8
:
38056
60
.
33.
Genovese
G
,
Kahler
AK
,
Handsaker
RE
,
Lindberg
J
,
Rose
SA
,
Bakhoum
SF
, et al
Clonal hematopoiesis and blood-cancer risk inferred from blood DNA sequence
.
N Engl J Med
2014
;
371
:
2477
87
.
34.
Steensma
DP
,
Bejar
R
,
Jaiswal
S
,
Lindsley
RC
,
Sekeres
MA
,
Hasserjian
RP
, et al
Clonal hematopoiesis of indeterminate potential and its distinction from myelodysplastic syndromes
.
Blood
2015
;
126
:
9
16
.
35.
McNerney
ME
,
Godley
LA
,
Le Beau
MM
. 
Therapy-related myeloid neoplasms: when genetics and environment collide
.
Nat Rev
2017
;
17
:
513
27
.
36.
Tumeh
PC
,
Harview
CL
,
Yearley
JH
,
Shintaku
IP
,
Taylor
EJ
,
Robert
L
, et al
PD-1 blockade induces responses by inhibiting adaptive immune resistance
.
Nature
2014
;
515
:
568
71
.
37.
Subudhi
SK
,
Aparicio
A
,
Gao
J
,
Zurita
AJ
,
Araujo
JC
,
Logothetis
CJ
, et al
Clonal expansion of CD8 T cells in the systemic circulation precedes development of ipilimumab-induced toxicities
.
Proc Natl Acad Sci U S A
2016
;
113
:
11919
24
.

Supplementary data