Clinical and biomarker analyses of sintilimab versus chemotherapy as second-line therapy for advanced or metastatic esophageal squamous cell carcinoma: a randomized, open-label phase 2 study (ORIENT-2)

This randomized, open-label, multi-center phase 2 study (NCT03116152) assessed sintilimab, a PD-1 inhibitor, versus chemotherapy in patients with esophageal squamous cell carcinoma after first-line chemotherapy. The primary endpoint was overall survival (OS), while exploratory endpoint was the association of biomarkers with efficacy. The median OS in the sintilimab group was significantly improved compared with the chemotherapy group (median OS 7.2 vs.6.2 months; P = 0.032; HR = 0.70; 95% CI, 0.50–0.97). Incidence of treatment-related adverse events of grade 3–5 was lower with sintilimab than with chemotherapy (20.2 vs. 39.1%). Patients with high T-cell receptor (TCR) clonality and low molecular tumor burden index (mTBI) showed the longest median OS (15.0 months). Patients with NLR < 3 at 6 weeks post-treatment had a significantly prolonged median OS (16.6 months) compared with NLR ≥ 3. The results demonstrate a significant improvement in OS of sintilimab compared to chemotherapy as second-line treatment for advanced or metastatic ESCC.


Results
Patient characteristics. Between May 16, 2017 and August 30, 2018, 253 patients were screened and 190 eligible patients were randomly assigned to receive either sintilimab or chemotherapy (n = 95 per group, Fig. 1). In total, 94 patients in the sintilimab group and 87 patients (15 with paclitaxel and 72 with irinotecan) in the chemo group received at least one dose of the assigned treatment ( Fig. 1). Both the demographics and disease characteristics at baseline were generally balanced between the two groups (Table 1).
After initial disease progression per RECIST v1.1, 37 patients in the sintilimab group continued to receive sintilimab (continuous subgroup), whereas treatment was discontinued for 28 patients (noncontinuous subgroup). A post hoc analysis of OS in the above subgroups indicated that the median OS of patients with continuous treatment was dramatically improved over those with discontinued treatment (median OS, 12 Table 2).
In the subgroup analysis for OS, the HRs favored sintilimab over chemo for the following five subgroups (sintilimab vs. chemo): age ≤65 years, male, previous treatment with paclitaxel, Eastern Cooperative Oncology Group performance-status (ECOG PS) score of 1, and current smokers (Fig. 3). However, no obvious differences were noted in either OS or PFS between the groups across PD-L1-expression subgroups (P > 0.05, Fig. 3 and Fig. S1).
Health-related quality of life was assessed in 87 patients in the sintilimab group and 84 patients in the chemo group at baseline. Based on the change from baseline, patients in the sintilimab group presented with better health-related quality of life from patients in the chemo group (Table S4).
Treatment-emergent adverse events (TEAEs) of any grade occurred in 88 (93.6%) patients with sintilimab, and 81 (93.1%) patients with chemo. Treatment-related adverse events (TRAEs) were reported in 54.3% of patients in the sintilimab group and 90.8% of patients in the chemo group (Table 3). The most common TRAEs were hypothyroidism (12.8%), pulmonary inflammation (10.6%), anemia (8.5%), and decreased whiteblood cell (WBC) counts (8.5%) in the sintilimab group. The incidence of serious TRAEs was similar between the sintilimab (18.1%) and chemo (20.7%) groups. The most common serious TRAEs were pulmonary inflammation (7.4%) and abnormal hepatic function (4.3%) in the sintilimab group. However, the incidence of TRAEs of grade 3 or worse was lower in the sintilimab group (20.2%) than in the chemo group (39.1%). The most common grade 3 or worse TRAEs were pulmonary inflammation (5.3%) and elevated lipase levels (4.3%) in sintilimab group. The number of patients discontinuing treatment because of TRAEs was 13 (13.8%) and 7 (8.0%) in the sintilimab and chemo groups, respectively. In addition, three (3.2%) deaths that occurred in the sintilimab group were attributed to treatment-related adverse events, including upper gastrointestinal bleeding, pneumonitis, and lung infection.
Biomarker analyses. In the sintilimab group, a post hoc analysis of efficacy in association with NLR was performed. Patients with NLR < 3 were assigned into the low NLR group and those with NLR ≥ 3 were assigned into the high NLR group. In total, 81 patients had evaluable NLR at baseline. Patients with a low NLR at baseline showed a significant improvement on OS (HR 0.54, P = 0.019; median 14.0 vs. 6.2 months) and PFS (HR 0.47, P = 0.002; median 2.9 vs. 1.5 months) compared with those with a high NLR at baseline. In total, 65 patients had evaluable NLR at 6 weeks post treatment. Notably, the survival benefit was also significant in patients with low versus high NLR, with a median OS of 16.6 vs. 6.2 months (HR 0.19, P < 0.001) and a median PFS of 4.3 vs. 2.3 months (HR 0.47, P = 0.006; Fig. 4).  In total, 118 patients (n = 64 in the sintilimab group and n = 54 in the chemo group) had qualifiable RNA-seq data for downstream analysis from tumor tissue. We analyzed an outcome association of 28 immune-cell populations and overall, only few potential positive correlations of immune-cell infiltrations with PFS but not OS could be observed. Only a significant correlation with PFS could be established for high infiltration of T-follicular helper cells (raw P = 0.007; HR, 0.46) and activated B cells (raw P = 0.030; HR, 0.54) in the sintilimab group (Table S1, Fig. 5).
The two signatures of T-follicular helper cells and activated B-cells were moderately correlated (raw P < 0.0001; Spearman's correlation, 0.48).
The outcome association of 45 signaling pathways was analyzed. Based on raw p-values, three pathways were significantly associated with improved OS only in the sintilimab group (Wnt, Met, and PI3K/Akt), one pathway was significantly associated with improved OS only in the chemo group (glucose metabolism) and one pathway was significantly associated with  To further explore indicators of clinical benefit in the periphery blood, T-cell receptor (TCR, n = 94) and circulating cell-free DNA (cfDNA, n = 83) sequencing was performed at baseline for patients in the sintilimab group in order to derive TCR clonality and molecular tumor-burden index (mTBI). Our analyses showed that patients with low mTBI had a significantly longer PFS (P = 0.018; HR, 0.55) compared with those with high mTBI. Similar trend was also observed for OS (P = 0.200; HR, 0.71). In contrast, TCR clonality alone could not efficiently predict the clinical outcome (Fig. S2). However, the prediction power for OS benefit could be enhanced when TCR clonality was combined with mTBI. Despite the limited sample size, the patients with high TCR clonality and low mTBI had the longest median OS (15.0 months) and PFS (4.1 months), respectively (Fig. 6).

Discussion
This study revealed that sintilimab monotherapy resulted in an evidently prolonged survival outcome with a favorable safety profile vs. chemotherapy as second-line treatment in Chinese patients with advanced or metastatic ESCC.
Up to now, four similar studies of PD-1 inhibitors have reported encouraging clinical efficacy in patients with ESCC after first-line therapy, including KEYNOTE-181 of pembrolizumab, ATTRACTION-3 of nivolumab, ESCORT of camrelizumab, and RETIONALE-302 of tislelizumab [12][13][14]19 . In these studies the median OS of PD-1 inhibitor groups are all greater than 7.0 months and the HRs for OS are all lesser than 0.8 [12][13][14]19 . The significant survival benefit of sintilimab over chemotherapy for patients with ESCC was consistent with the other studies. The Kaplan-Meier survival curves for OS were shown to diverge beyond 5 months in favor of sintilimab. This was similar to the results obtained in the other four studies, which was the Kaplan-Meier survival curves for OS diverged beyond 3-6 months [12][13][14]19 .  Data are presented in %. Grade 1-2 TRAE listed with an incidence of ≥5% of patients in either treatment group, and grade 3-5 TRAEs with an incidence of ≥2% in either group. IrAEs were occurred in ≥2% of patients in either treatment group. TRAE treatment-related adverse event, irAE immune-related adverse event, WBC white blood cell, ALT alanine aminotransferase, AST aspartate aminotransferase.
RMST is a statistical technique recently adopted widely as an alternative to median to provide an intuitive interpretation of the data. It is especially more appealing than median survival time when the data violate proportional hazard assumption when only evaluating median difference could be very misleading. The superior RMST of sintilimab over chemo was increased from 0.3 months at 9 months after randomization to 1.6 months at 18 months after randomization, further suggesting a delayed survival benefit of sintilimab. In addition, P-values of different weight coefficients (0, 0.2), (0, 0.5), and (0, 1) from Fleming-Harrington test are less than 0.01, indicating the delayed effect of the significant OS benefit with sintilimab over chemo. Consequently, it is suggested that there are a delayed response and more durable benefit of immunotherapies based on similar studies and different statistical methods. The survival advantage of sintilimab as second treatment for patients with ESCC may be demonstrated via a prolonged follow-up duration.
Moreover, the ORR of patients with sintilimab was twice higher than those with chemo, similar to the findings in the other four studies [12][13][14]19 . Notably, for patients who were initially evaluated as PD per RECIST v1.1 in the sintilimab group, those continuing to receive sintilimab had a remarkable increase on median OS, compared with those who discontinued treatment. This finding suggested that continuous treatment with sintilimab could have a long-term survival benefit and should be further explored for patients after initial PD by RECIST v1.1.
Sintilimab showed a favorable safety profile over chemo, with a substantially reduced incidence of TRAEs (54.3 vs. 90.8%) and TRAEs of grade 3 or worse (20.2 vs. 39.1%). Most toxic effects were comparable to the historical data on treatment with sintilimab in other indications and to that of treatments with other PD-1 inhibitors in esophagus cancer [12][13][14][15][16][17][18][19] . In particular, three treatment-related deaths, which were due to upper gastrointestinal bleeding, pneumonitis, and lung infection, occurred in the sintilimab group. Pneumonitis and lung infection were previously reported in PD-1 inhibitor-treated patients with esophagus cancer [12][13][14]19 . Nevertheless, upper gastrointestinal bleeding is one of the most common fatal causes for ESCC. Upper gastrointestinal bleeding was considered to be treatment-related, perhaps because this study was an open-label study.
In this study, comprehensive biomarker analyses were performed, including tumor PD-L1 expression, NLR, tumor transcriptome analysis and peripheral blood TCR, and cfDNA analysis.
Even if only the HR for OS in the combined-proportion score (CPS) ≥ 10 subgroups was less than the CPS < 10 subgroups, neither PD-L1-expression level assessed by tumor-proportion score (TPS) nor CPS could predict definitely benefit from the treatment with sintilimab unlike other studies [12][13][14]19 . That is potentially attributed to limited sample availability.
NLR is an inflammatory-response biomarker for the prognosis of ESCC in patients treated with chemoradiotherapy and immune-checkpoint inhibitors [20][21][22][23] . As is known, NLR can reflect differences in the immune status during cancer development and progression 22,23 . A meta-analysis showed that high NLR is associated with worse survival in esophageal cancer, as it is also shown in several solid tumors 24 . Our results revealed that low NLR (<3) at baseline or at 6 weeks post treatment with sintilimab was significantly correlated with a longer OS and PFS. The increased number of neutrophils in the peritumoral inflammatory-cell infiltrate can reduce the antitumor activity of natural-killer (NK) and activated T cells 25,26 . Moreover, many cytokines (e.g., TNF, IL-1, and IL-6) and VEGF produced as a result of neutrophil activation may increase tumor growth 27 . But   the mechanism is still unclear that NLR is related to survival in esophageal cancer and other solid tumors. The test of NLR is very simple and practicable, but the use of NLR in clinical practice is still underestimated.
Although general T-cell tumor infiltration did not show a predictive effect in the analysis of tumor transcriptomes, infiltration of T-follicular helper cells and activated B cells was significantly related to PFS in the sintilimab group and not in the chemo group. Especially, the predictive value of activated B lymphocytes is supported by an emerging role of B cells and tertiary lymphoid structures (TLS) in tumors, which have been recently reported in association with survival, even in tumor with low T-cell infiltration [28][29][30] . Furthermore, T-follicular helper cells have been shown to play important roles in the generation of TLS in the tumor microenvironment 31 , which might be relatively enriched for B-cell signatures and TLS in ESCC 32 . Thus, the potential relevance of TLS in the mode of action of cancer immunotherapy to mount an efficient antitumor response is supported by our study and warrants further investigation. However, due to tumor-specimen limitations, additional analyses to characterize TLS and their spatial distribution within the tumor microenvironment could not be performed. The correlation of several signaling pathways with clinical efficacy is more challenging to interpret. Different reports indicate a potentially context-dependent role for Wnt/β-catenin signaling in CD8 + Tcell differentiation and memory formation 33,34 . Wnt pathway also enhances the proliferation and cytolytic activity of human γδT cells 35 . Autocrine Wnt pathway activation induces vascular normalization and endows endothelial-cell basement membrane the capacity to boost T-cell infiltration 36 . Further evidence shows that, in contrast to early-stage tumors, advanced and metastatic tumors have downregulated Wnt expression. In some cases, endogenous Wnt signaling even generally counteracts tumor progression 37 . Also, the small sample size could contribute to these surprising findings. Such conflicting data clearly underline the need for additional analyses and validation of the data.
In contrast to the analysis of tumor tissue, the investigation of peripheral blood has practical advantages without the limitation of tumor tissue accessibility and size. We found that peripheral TCR clonality at baseline did not predict the efficacy of sintilimab. In addition, we also analyzed mTBI as an indicator of tumor burden at the molecular level 38 . Whereas mTBI is derived from peripheral blood, tumor mutation burden (TMB) is generally derived from mutation assessment in tumor tissue. In contrast to TCR clonality, mTBI index was significantly associated with clinical efficacy in the treatment with sintilimab. However, the predictive value of mTBI in combination with TCR clonality might be further enhanced as patients with low mTBI and high TCR clonality exhibited the longest survival in the sintilimab group. TCR-clonality analyses revealed conflicting results in previous studies, potentially due to the dilution effect of nontumor-specific TCRs since this analysis alone cannot distinguish between tumor vs. nontumor specificity [39][40][41] . However, a high clonality might be indicative of a preexisting and expanded tumor-specific T-cell response that has been blocked or exhausted. It is intriguing to speculate that such a block is mediated by the PD-1 axis and sintilimab treatment might release this block. However, reinvigoration and subsequent tumor eradication seems to mainly occur in patients with low mTBI. Detailed mechanisms for this correlation need to be examined, but apart from physical barriers in larger tumors it was previously reported that the ratio of T-cell reinvigoration to baseline tumor burden instead of the magnitude of the reinvigoration alone might be a better indicator of response 42,43 . Additionally, it has been reported that skewed ratios of infiltrating suppressive immune cells like myeloidderived suppressor cells and T regulatory cells vs. T-effector cells, dendritic cells, and natural-killer cells correlate with tumor progression and increased tumor burden. Such cellular changes are accompanied with increasingly immune-suppressive cytokine milieus in the growing tumors [44][45][46][47][48] . Our data further support the notion that an effective mounting of an antitumor response is multifactorial and does not only rely on successful reinvigoration of exhausted T cells, but also on tumor burden and preexisting TCR clonality in addition to previously described predictive biomarkers like PD-L1 expression, TMB, and T-cell infiltration. It remains to be investigated if these factors are independent contributors or surrogates of each other.
This study had several limitations. First, the sample size of patients was relatively small. Second, the open-label design of the study might have influenced the assessment of the incidence of adverse effects. Nevertheless, this design was deemed acceptable because both dose regimens and toxicities in the two groups were disparate. Third, it was hard to evaluate the correlation between the expression of PD-L1 and the survival benefit of sintilimab due to the limited tumor samples available for the assessment of PD-L1 expression.
In conclusion, the results of ORIENT-2 study favored the use of sintilimab over chemo in Chinese patients with advanced ESCC refractory to previous chemotherapies, as it suggested a prolonged survival benefit and a favorable safety profile. Both NLR <3 at 6 weeks post treatment and the combination of high TCR clonality with low mTBI may be potent biomarkers for the prediction of improved OS and PFS in patients with ESCC treated with sintilimab.

Methods
Study design and patients. This was a multicenter, randomized, parallel, openlabel phase-2 study conducted at 30 sites in China (Table S5). Patients with histopathologically or cytologically confirmed locally advanced or metastatic ESCC aged between 18 and 75 years old were enrolled, including male and female. Major eligibility criteria were at least one measurable lesion per Response Evaluation Criteria in Solid Tumors (RECIST), version 1.1, Eastern Cooperative Oncology Group (ECOG) performance status (PS) of 0 or 1, eligible to provide a fresh or archived tumor sample, and radiological or clinical evidence of disease progression during or after first-line chemotherapy. Patients were excluded who received prior PD-1 or PD-L1 inhibitors, or received any radiotherapy, immunosuppressive drugs, or the study drug within 4 weeks prior to the first dose. Full inclusion criteria are provided in the Supplementary. The study was performed in accordance with the Declaration of Helsinki and Good Clinical Practice guidelines and was approved by the ethics committee at each site. All patients provided written informed consent prior to enrollment. The trial is registered with Clinical-Trials.gov, number NCT03116152.
Procedures. Patients were randomly (1:1) assigned to receive either sintilimab or the investigator's choice of chemo (paclitaxel or irinotecan), using an interactive web-response system with a block size of a mixture of 2, 4, and 6, and with the stratification factor of the ECOG PS score (0 vs. 1). Neither investigators nor patients were blinded to treatment allocation.
After randomization, patients were treated with 200 mg of intravenous sintilimab, once every 3 weeks, or with the investigator's choice of chemo (175 mg/m 2 intravenous paclitaxel once every 3 weeks, or 180 mg/m 2 intravenous irinotecan once every 2 weeks), until disease progression, death, unacceptable toxicity, or withdrawal of informed consent.
Tumor assessment was performed by investigators according to RECIST version 1.1 every 6 weeks till 24 weeks after randomization and every 9 weeks thereafter until initiation of new antitumor therapy, disease progression, withdrawal of informed consent, or death. In the sintilimab group, patients with stable status and initial progressive disease per RECIST version 1.1 assessed by the investigator could continue to receive sintilimab. Then tumor assessment was performed according to modified RECIST v1.1 for immune-based therapeutics (iRECIST). Adverse events (AEs) and serious adverse events (SAEs) were monitored throughout the treatment period and for 90 days after the last dose, based on the National Cancer Institute Common Terminology Criteria for Adverse Events (NCI CTCAE v4.03). Measurements of NLR. The NLR inflammatory biomarker, which was calculated as neutrophil count/lymphocyte count, was obtained from hematology panels at baseline and 6 weeks post treatment. High and low NLR groups were classified using the cutoff of NLR = 3 49 .

Measurements
Biomarker analysis using next-generation sequencing and bioinformatics methods. Quality control of RNAseq raw data was performed using FastQC (version 0.11.8) 50 . Clean reads were obtained through trimming clean reads using Trimmomatic (version 0.36) 51 . STAR aligner (version 2.7.0a) was used to map clean reads to human reference genome (hg38) 52 . HTSeq (version 0.11.4) was used for counting reads mapped to genes 53 . The gene sets defining 28 immune-cell populations were extracted from the previous research 38 . The gene sets of 45 signaling pathways were downloaded from the tumor signaling 360 panel of NanoString (Seattle, Washington, US). The tumor samples of biomarker-evaluable population (BEP) with qualified RNA-seq data were used to conduct downstream analysis. To screen for survival-related genes, based on transcript-per-million (TPM)-normalized gene expression level of each protein-coding gene, the BEP was iteratively split into high or low expression groups (S high_low ) by median value of the whole cohort of the respective gene. Then, Cox regression model by the formula Surv (time, event)~S high_low in each of the treatment arm was run using the R survival package (version 2.44-1.1), to derive the respective hazard ratio and raw P-value 54 . All genes were ranked by raw P-value, whereas genes with P < 0.05 were considered as potentially survival-related genes. Gene-set signature scores were calculated using the GSVA algorithm (version 1.32.0) based on the gene expression levels of TPM 55 . High and low signature scores were split according to the median value of the whole cohort. Then the same screening procedure as described above was conducted.
Peripheral baseline T-cell receptor (TCR) CDR3 regions and cell-free DNA were sequenced to calculate TCR clonality and molecular tumor-burden index (mTBI) in each sample. The latter is a reflection of the percentage of circulating-tumor DNA (ctDNA) detected in whole cell-free DNA (cfDNA) and is considered to reflect the tumor burden at the molecular level, which is defined as the mean allele fraction of mutations in the trunk cluster as published earlier 38 . Sequencing data of ctDNA were mapped to human reference genome (hg19) by using BWA (version 0.7.17) 56 . GATK (version 4.1.4.1) was used to call variants of each sample 57 . PyClone (version 0.13.1) was used to analyze the clonal-population structures based on baseline-mutation profile 58 . Sequencing reads of TCRβ CDR3 region were merged to obtain contigs using Pear (version 0.9.11) 59 , and then contig sequences were aligned to reference TRB V/(D)/J gene sequences (http://www.imgt.org) using MiXCR (version 3.0.4) to determine the TRB V/(D)/J gene segment in each contig 60 . High and low TCR clonality or mTBI of baseline samples were split according to the median value of the whole cohort. Then, Cox regression model by the formula Surv(time, event)~S high_low in each of treatment arm was run using the R survival package (version 2.44-1.1) 54 , to derive the respective hazard ratio (HR) and raw P-value (P). The 'time' represents PFS or OS, 'event' for progression or death event, and 'S high_low ' for split groups by median value. Survival curves were plotted using the R survminer package (version 0.4.6) 61 .
Outcomes. The primary endpoint was OS, defined as the time from randomization until death. Surviving patients were censored using the data at the last follow-up. Secondary endpoints were PFS, which was the time from randomization to initial disease progression evaluated by the investigator per RECIST v1.1 or death; ORR, defined as the percentage of patients achieving complete response (CR) or partial response (PR); DCR, which was the percentage of patients with the best overall response of CR, PR, or stable disease (SD); time to response (TTR); and duration of response (DOR). As safety endpoint was considered the incidence of all AEs across the study. Exploratory endpoints were the association between efficacy (OS and PFS) and biomarkers, including PD-L1 expression, mRNA expression, and ctDNA results. The association between efficacy (OS, PFS) and NLR was post hoc analysis. Health-related quality of life (HRQoL) was assessed by the European Organisation for Research and Treatment of Cancer Quality-of-Life Questionnaire-Core 30 (EORTC QLQ-C30), and the five-level version of European Quality of Life-5 Dimensions (EQ-5D-5L) visual analogue scale (VAS), at the first dose, as well as at every radiological evaluation time and at the time of the first safety follow-up.
Statistical analysis. The required death-event number was 142, with an expected HR for OS of 0.7 (α = 0.2, two-sided) under a statistical power of 80%. Assuming a 21% dropout rate, 180 patients were required.
All efficacy endpoints were assessed in the intention-to-treat (ITT) population that included all randomly assigned patients. The safety profile was assessed in the safety set (SS), which included patients receiving at least one dose of the study drug.
The Kaplan-Meier method was utilized to evaluate the time-to-event endpoints (median OS, PFS, and DOR), and the differences in the survival curves between the two groups were analyzed using the stratified log-rank test by ECOG PS score. A stratified Cox proportional hazard regression model was used to estimate the HR and corresponding 95% confidence interval (CI) of OS and PFS between treatment groups, using the ECOG PS score (0 vs. 1) as a stratification factor. The RMST was evaluated to compare the differences in OS between groups at 9, 12, 15, and 18 months post treatment. The binomial distribution method was used to evaluate the 95% CI of ORR and DCR, and differences were compared using the Cochran Mantel-Haenszel test. Patient-reported health-related quality-of-life outcomes were analyzed for patients who had a baseline and at least one post-baseline assessment of EORTC QLQ-C30 (general health) and EQ-5D-5L (VAS and index).
To evaluate any factors correlated with the crossing of the survival curves in OS, different treatments were assessed using a weighted log-rank test from the FH G (ρ-γ) , which accounts for nonproportional hazards (ρ, γ = 0,1; 0, 0.5; 0, 0.2; 1, 0). Considering the potential pseudoprogression for PD-1 inhibitors, a post hoc analysis was performed to compare the difference in OS between patients with and without continuous treatment of sintilimab after initial PD per RECIST v1.1 in the sintilimab group.
Statistical analyses for the clinical part were performed using the SAS software (version 9.4) (Cary, North Carolina, US). The significance level for all endpoints was α = 0.05 (two-sided) and the significance threshold was P < 0.05. For biomarker identification using bioinformatics, survival analysis was performed using the survival package 43 in R, by which both the P-value and hazard ratio between two groups were calculated using the coxph function. Survival curves were plotted using the survminer package 44 .
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
All study data are presented in the paper and supplementary files. Source data are provided with this paper, including the baseline patient information, the results of efficacy in all patients, the NLR, and the results of RNA sequencing, TCR clonality, and molecular tumor-burden index (mTBI) in the sintilimab group. The raw RNA-seq, TCRseq, and ctDNA-seq reads generated in this paper are deposited in Genome Sequence Archive (GSA) under the accession code HRA001773. The raw sequence data are available under restricted access because of data-privacy laws. Source data are provided with this paper.