The predicting role of circulating tumor DNA landscape in gastric cancer patients treated with immune checkpoint inhibitors
Molecular Cancer volume 19, Article number: 154 (2020)
A more common and noninvasive predicting biomarker for programmed cell death 1 (PD-1) antibody remains to be explored. We assessed 46 patients with advanced gastric cancer who received PD-1 antibody immunotherapy and 425-genes next-generation sequencing (NGS) testing. Patients who had a > 25% decline in maximal somatic variant allelic frequency (maxVAF) had a longer progression free survival (PFS) and higher response rate than those who did not (7.3 months vs 3.6 months, p = 0.0011; 53.3% vs 13.3%, p = 0.06). The median PFS of patients with undetectable and detectable post-treatment circulating tumor DNA (ctDNA) was 7.4 months vs. 4.9 months (p = 0.025). Mutation status of TGFBR2, RHOA, and PREX2 in baseline ctDNA influenced the PFS of immunotherapy (p < 0.05). Patients with alterations in CEBPA, FGFR4, MET or KMT2B (p = 0.09) gene had greater likelihood of immune-related adverse events (irAEs). ctDNA can serve as a potential biomarker of the response to immunotherapy in advanced gastric cancers, and its potential role in predicting irAEs worth further exploration.
Gastric cancer is the fifth most common malignant tumor and the third leading cause of cancer-related death worldwide . The immune Checkpoint Inhibitors (ICIs), mainly the antibodies against PD-1, have been recommended as a palliative therapy option for selective patients with metastatic gastric adenocarcinoma, with predictive biomarkers including the microsatellite instability (MSI) status, Combined Positive Score (CPS) of PD-L1 expression, and Epstein-Barr virus (EBV) status of tumor [2, 3]. Our team have identified tumor mutation burden (TMB) as a biomarker for OS benefit in chemo-refractory gastric cancer treated with PD-1 antibody . However, in real-world practice, patients with advanced gastric cancer may lack fresh tissue and have to obtained tissue samples from invasive re-biopsy.
ctDNA has been reported to have utility in identifying genetic alterations and predicting the prognosis, identifying resistance of target therapy, and monitoring relapse or progression of gastric cancer . What’s more, ctDNA could provide longitudinal and dynamic surveillance of the tumor-specific genetic characteristics without repeatedly performing invasive tumor biopsy that costs more time and money. Recently, the use of dynamic ctDNA to predict the response of melanoma and non-small cell lung cancer (NSCLC) to checkpoint inhibitors has been reported [6, 7], while in gastric cancer little has been explored. In this study, we aim to explore the predicting role of ctDNA in gastric cancer patients treated with immunotherapy.
Results and discussion
Clinicopathologic characteristics and genomic landscape of ctDNA-NGS
Totally, forty-six eligible patients with metastatic gastric cancer were enrolled between October of 2018 and December of 2019(Additional file 1: Supplementary materials and methods). Baseline clinicopathologic characteristics are summarized in Additional file 2: Table S1 and the procedure of data analysis was described in Additional file 3: Fig. S1; It comprised of 30 males and 16 females. The majority (78.3%) of patients had metastatic sites of 1 or 2, and common metastatic site was peritoneum (71.7%). In 26 patients with available PD-L1 status, 13 patients had tumors with PD-L1 CPS ≥ 1, and 7 with CPS ≥ 10; 25 patients had an EBV in situ hybridization test, all of which were negative. One patient was MSI-H, and 4 patients were HER-2 positive. More than half (58.7%) of the patients received chemotherapy and immunotherapy as first-line palliative therapy, while others received the treatment as second-line or late line therapy. The chemotherapy included fluorouracil, capecitabine, S-1, albumin-bound paclitaxel, irinotecan, and trastuzumab. The PD-1 antibody included Nivolumab, Pembrolizumab, Toripalimab, and Sintilimab.
Plasma circulating tumor DNA sequencing (ctDNA-NGS) and tissue tumor DNA sequencing (tissue-NGS) results were obtained using the commercially available 425-gene NGS panel (Fig. 1, Additional file 4: Fig. S2). Of all 46 patients, 38 had tissue samples and 43 had baseline blood samples (Additional file 5: Fig. S3). Among these 43 patients evaluated for baseline ctDNA, 88.4% had at least one genomic alteration (Additional file 5: Fig. S3, Additional file 6: Fig. S4) and 11.6% had undetectable baseline ctDNA. The consistency of ctDNA-NGS with tissue-NGS were also explored (Additional file 7: Fig. S5).
Decreasing ctDNA was correlated to higher response to immunotherapy
We further analyzed the ctDNA data of patients who had detectable baseline ctDNA and underwent serial ctDNA-NGS (N = 32) and found that patients who had a > 25% decline in maxVAF (N = 17) had a significantly longer PFS than those who did not (7.3 months; 95%CI, 2.4–4.8 months vs 3.6 months; 95%CI 4.6–10.0 months; p = 0.0011; time between serial blood samples collection: median 68 days; rang 19–252 days; Fig. 1). The ORR of patients who had a > 25% decline in maxVAF was 53.3%, higher than those who did not (13.3%; p = 0.06). These findings were in consistency with the results of a post hoc exploratory analysis of a prospective phase 2 clinical trial, the change in ctDNA post-treatment were reported to predict response to pembrolizumab, and 14 gastric cancer patients with decreasing ctDNA demonstrated significant improvements in DCR, ORR, and PFS . We identified 10 patients with undetectable and 30 patients with detectable ctDNA of post-treatment plasma samples. The median PFS of patients with undetectable and detectable post-treatment ctDNA was 7.4 months (95%CI, 4.7–10.1 months) vs. 4.9 months (95%CI, 3.1–6.6 months; p = 0.025; Fig. 1), respectively.
Molecular landscape of ctDNA predicts resistance to immunotherapy
Response evaluations were available for all patients with a median follow-up of 10.4 months. One (2.1%) patient achieved complete response (CR) with the pathological result of the resected tumor indicating a pathological CR (pCR), ten (21.7%) achieved partial response (PR), twenty-twelve (47.8%) had stable disease (SD), six (13.0%) were nonCR/nonPD, and seven (15.2%) patients were progressive disease (PD). In 26 patients with available PD-L1 status, the median PFS of those with PDL1 CPS ≥ 10 and CPS < 10 were 3.4 and 4.9 months (p > 0.05). In the patients received first-line treatment, the median PFS of those with PDL1 CPS ≥ 10 and CPS < 10 were 7.8 and 7.4 months(p > 0.05). In the second-line or late line setting, there’s also no significant difference between the median PFS in patients with different CPS, which may be due to the small number of patients in each subgroup. Of 43 baseline ctDNA-NGS, we excluded 5 cases with undetectable ctDNA to further explore the correlation of baseline ctDNA-NGS with the duration of immunotherapy, and found the mutation status of TGFBR2, RHOA, and PREX2 influenced the PFS of immunotherapy (Fig. 2). TGFBR2 p.P525L, p.P129Afs*3, p.E269* mutations had been identified and the median PFS of patients with TGFBR2wt and TGFBR2mt were 5.0 months (95%CI, 2.5–7.4 months) and 1.6 months (95%CI, 1.4–1.8 months; p < 0.001), respectively. The genomic alteration of the TGFBR2 gene that discovered in our study as a candidate biomarker of immunotherapy had been reported to be related to the responses of anti-pd-L1 treatment (atezolizumab) in patients with metastatic urothelial cancer . RHOA p.L57V, p.L69P, p.D45N mutations had been identified and the median PFS of patients with RHOAwt and RHOAmt were 4.9 months (95%CI, 2.5–7.4 months) and 2.4 months (95%CI, 0.8–4.0 months; p = 0.0056), respectively. PREX2 p.M627V, p.W519G, p.R1567Q mutations had been identified and the median PFS of patients with PREX2wt and PREX2mt were 5.0 months (95%CI,2.0–8.0 months) and 2.4 months (95%CI, 0.5–4.3 months; p = 0.037), respectively. The PREX2 mutations were reported to enhance cell proliferation in melanoma, hepatocellular carcinoma and pancreatic ductal adenocarcinoma, while its role in immunotherapy had not been explored. Patients with PREX2 mutation tended to have lower PD-L1 CPS and fewer CD8+ T cells (Although not statistically significant, it might due to the limited PREX2 mutant samples), which indicated PREX2 mutation might have potential influence on tumor immune microenvironment and this is worthy of further exploring (Additional file 8: Fig. S6).
To explore the acquired resistance of immunotherapy, we collected plasma samples at PD of the patients who had benefit from the treatment (achieved CR, PR, or SD) and the plasma samples of 16 patients were available. Most of the ctDNA-NGS at PD were largely increased compared with baseline or ctDNA from plasma collected last time (Additional file 9: Table S2), and new genomic alterations were observed in 13 patients. We further compared the newly emerged genomic alterations in ctDNA-NGS with their baseline tissue DNA-NGS. Of the 13 patients with new alterations in plasma, excluding 3 patients that did not have tissue samples, 7 patients had new genomic alterations not found in baseline tissue, and 3 patients had new alterations both could not be found and could be found in baseline tissue. Mutation of FOXL2 gene and CNV of FGFR2 were identified as new alterations in two patients, and RHOA mutations that were discussed above was also found in one patient with acquired resistance, suggesting that these alterations were candidate genes involving in the acquired resistance to immunotherapy.
Immune-related adverse events of immunotherapy
In this study, most common irAEs were endocrine (10.9%) and hepatic toxicities (8.7%). Severe irAEs that led to termination of immunotherapy were occurred in 2 (4.3%) patients. Some literatures suggested onset of irAE was a biomarker for immunotherapy response, and patients experienced irAEs showed improvements in PFS, OS and ORR . However, in our studies, no marked difference of PFS or ORR were observed between patients with and without irAEs. We investigated the correlations between ctDNA and irAEs, and found that patients with alterations in CEBPA (100% vs 27.3%, p = 0.09), FGFR4(100% vs 27.3%, p = 0.09), MET (100% vs 27.3%, p = 0.09), or KMT2B(100% vs 27.3%, p = 0.09) in baseline plasma had greater likelihood of experiencing irAEs, which requires further independent validation in a clinical trial setting with more patients enrolled. Previous studies showed that plentiful biomarkers, including gender, serum cytokines, immune cells, and autoantibodies, were associated with the occurrence of high-grade irAEs . However, there is not any biomarker with high grade of evidence level, more predictive biomarkers of irAEs deserved to be investigated.
In conclusion, the result of our study indicated that dynamic ctDNA can serve as a potential biomarker of the response to immunotherapy in advanced gastric cancers, and its potential role in detecting the resistance mechanisms and predicting irAEs worth further exploration.
Availability of data and materials
All the data that support the findings of this study are available from the corresponding authors upon reasonable request.
Programmed cell death 1
Maximal somatic variant allelic frequency
Circulating tumor DNA
Progression free survival
Immune-related adverse events
Immune Checkpoint Inhibitors
Combined Positive Score
Tumor mutation burden
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Kim ST, Cristescu R, Bass AJ, Kim KM, Odegaard JI, Kim K, et al. Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancer. Nat Med. 2018;24(9):1449–58.
Wang FH, Shen L, Li J, Zhou ZW, Liang H, Zhang XT, et al. The Chinese Society of Clinical Oncology (CSCO): clinical guidelines for the diagnosis and treatment of gastric cancer. Cancer Commun. 2019;39(1):10.
Wang F, Wei XL, Wang FH, Xu N, Shen L, Dai GH, et al. Safety, efficacy and tumor mutational burden as a biomarker of overall survival benefit in chemo-refractory gastric cancer treated with toripalimab, a PD-1 antibody in phase Ib/II clinical trial NCT02915432. Ann Oncol. 2019;30(9):1479–86..
Parikh AR, Leshchiner I, Elagina L, Goyal L, Levovitz C, Siravegna G, et al. Liquid versus tissue biopsy for detecting acquired resistance and tumor heterogeneity in gastrointestinal cancers. Nat Med. 2019;25(9):1415–21.
Forschner A, Battke F, Hadaschik D, Schulze M, Weissgraeber S, Han CT, et al. Tumor mutation burden and circulating tumor DNA in combined CTLA-4 and PD-1 antibody therapy in metastatic melanoma - results of a prospective biomarker study. J Immunother Cancer. 2019;7(1):180.
Li L, Wang Y, Shi W, Zhu M, Liu Z, Luo N, et al. Serial ultra-deep sequencing of circulating tumor DNA reveals the clonal evolution in non-small cell lung cancer patients treated with anti-PD1 immunotherapy. Cancer Med. 2019;8(18):7669–78.
Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, et al. TGFbeta attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. 2018;554(7693):544–8.
Das S, Johnson DB. Immune-related adverse events and anti-tumor efficacy of immune checkpoint inhibitors. J Immunother Cancer. 2019;7(1):306.
Yao L, Jia G, Lu L, Bao Y, Ma W. Factors affecting tumor responders and predictive biomarkers of toxicities in cancer patients treated with immune checkpoint inhibitors. Int Immunopharmacol. 2020;85:106628.
The authors thank Qi-nian Wu from Department of Pathology, and colleagues from State Key Laboratory of Oncology in South China, Sun Yat-sen University Cancer Center for study support, suggestion of study design and oversight.
Funding for this study was provided by National Key R&D Program of China (2018YFC1313300); Natural Science Foundation of Guangdong (2019A1515011786); Science and Technology Program of Guangdong (2019B020227002); CAMS Innovation Fund for Medical Sciences (CIFMS) (2019-I2M-5-036); Major Science and Technology Project for Significant New Drugs Creation (2018ZX09734003).
Ethics approval and consent to participate
This study was reviewed and approved by the Ethnics Committees of Sun Yat-sen University Cancer Center. The study was conducted in accordance with the International Ethical Guidelines for Biomedical Research Involving Human Subjects.
Consent for publication
All authors have agreed to publish this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary materials and methods.
Baseline clinicopathologic characteristics.
CONSORT diagram. CONSORT diagram of 46 patients enrolled and samples analyzed.
425-gene NGS panel that used for plasma circulating tumor DNA sequencing and tissue tumor DNA sequencing.
The landscape of high-frequency genomic alterations detected in 38 tissue samples and 38 baseline plasma samples. High-frequency means genomic alterations detected in tissue or baseline plasma > 5%.
The characteristics of 38 patients with available baseline plasma samples.
Scatter plot between tTMB and bTMB. The Spearman’s rank test showed significant correlations. (p < 0.05) a Scatter plot between tTMB and bTMB in 24 patients who received immunotherapy as first-line treatment. b Scatter plot between tTMB and bTMB in all 38 patients whose tissue and baseline plasma samples were available.
Pathological examination of the patient’s gastric carcinoma in PREX2 wild type and mutant patients. a,d Representative section of H.E. staining showing gastric carcinoma in wild type and mutant patients respectively. Scale bar: 100 μm. b,e The PREX2 expression in tumor and non-tumor regions in wild type and mutant patients respectively. Scale bar: 100 μm. c,f CD8+ T cell densities in wild type and mutant patients respectively. Scale bar: 100 μm. g Statistical analysis of PD-L1 expression in 25 MSS patients by combined positive score (CPS). h Comparisons of the number of CD8+ T cells in PREX2 wild type and mutant patients. Student’s t-test were used. p > 0.05 was not considered statistically significant.
The dynamic changes of ctDNA of patients had benefit from the treatment.
About this article
Cite this article
Jin, Y., Chen, DL., Wang, F. et al. The predicting role of circulating tumor DNA landscape in gastric cancer patients treated with immune checkpoint inhibitors. Mol Cancer 19, 154 (2020). https://doi.org/10.1186/s12943-020-01274-7