3.1 Search results and study characteristics
The literature search originally yielded 96 identified articles from three following electronic databases (PubMed, EMBASE, and Cochrane central). There were 60 records left after the duplicate was discarded. After the screening of titles and abstracts, 18 articles were further scrutinized by a full-text screen. We ultimately conceal 13 articles for different reasons (i.e. Conferences or meeting posters, reviews, letters, comments, and subgroups of included RCTs) based on the inclusion and exclusion criteria. 5 trials were finally included in the quantitative meta-analysis. The flow diagram of the search and selection process is shown in Fig. 1.
The 5 included articles (1 phase II and 4 phase III) were all designed as RCTs with low bias. The sample size varied from 190 to 512, with a total of 1,970 subjects. Among them, 987 received PD-1 inhibitors including Pembrolizumab, Nivolumab, Camrelizumab, Tirelizumab, and Sintilimab, and the other 983 received standard chemotherapy including paclitaxel, irinotecan, and docetaxel. The total number of males in the study subjects was 1,705 (86.5%). Table 1 presents the study characteristics of these five studies.
Table 1
characteristics of the studies included in the meta-analysis
Trial ID
|
Design(phase), histology
|
Period of enrollment
|
Treatment arms
|
No. of patients for analysis (male%)
|
OS
|
PFS
|
ORR
|
DCR
|
Incidence of TRAEs
|
Any Grade
|
3–5 Grade
|
ATTRACTION-3 (NCT02569242), 2019
|
RCT(III), Squamous
|
Jan 7, 2016-May 25, 2017
|
Nivolumab
|
210(85.2%)
|
10.9 (9.2–13.3)
|
1.7 (1.5–2.7)
|
19.2%
|
37.4%
|
65.6%
|
18.2%
|
(1) Paclitaxel
(2) Docetaxel
|
209 (88.5%)
|
8.4 (7.2–9.9)
|
3.4 (3.0-4.2)
|
21.5%
|
62.7%
|
95.2%
|
63.9%
|
ESCORT (NCT03099382), 2020
|
RCT(III), Squamous
|
May 10, 2017-July 24, 2018
|
Camrelizumab
|
228 (91.2%)
|
8.3 (6.8–9.7)
|
1.9 (1.9–2.4)
|
20.2%
|
44.7%
|
94.3%
|
19.3%
|
(1) Docetaxel
(2) Irinotecan
|
220 (87.2%)
|
6.2 (5.7–6.9)
|
1.9 (1.9–2.1)
|
6.4%
|
34.5%
|
88.5%
|
39.5%
|
KEYNOTE-181 (NCT02564263), 2020
|
RCT(III), Squamous and adenomatous
|
Dec 8, 2015-June 16, 2017
|
Pembrolizumab
|
198 (337, 84.0%)
|
8.2 (6.7–10.3)
|
2.2 (2.1–3.2)
|
16.7%
|
46.0%
|
64.3%
|
18.2%
|
(1) Paclitaxel
(2) Docetaxel
(3) Irinotecan
|
203
|
7.1 (6.1–8.2)
|
3.1 (2.2–3.9)
|
7.4%
|
49.8%
|
88.5%
|
40.9%
|
RATIONALE 302 (NCT03430843), 2021
|
RCT(III), Squamous
|
January 2018-March 2020
|
Tirelizumab
|
256 (84.8%)
|
8.6 (7.5–10.4)
|
1.6 (1.4–2.7)
|
20.3%
|
46.9%
|
73.3%
|
18.8%
|
(1) Paclitaxel
(2) Docetaxel
(3) Irinotecan
|
256 (84.0%)
|
6.3 (5.3-7.0)
|
2.1 (1.5–2.7)
|
9.8%
|
41.8%
|
93.8%
|
55.8%
|
ORIENT-2 (NCT03116152), 2022
|
RCT(II), squamous
|
May 16, 2017-August 30, 2018
|
Sintilimab
|
95 (92.6%)
|
7.2 (5.8–9.7)
|
1.6 (1.5–2.8)
|
12.8%
|
44.7%
|
54.3%
|
20.2%
|
(1) Paclitaxel
(2) irinotecan
|
95 (88.4%)
|
6.2 (5.4–7.9)
|
2.9 (2.6–3.6)
|
6.3%
|
43.2%
|
90.8%
|
39.1%
|
3.2 Clinical efficacy and Safety
As depicted in Fig. 2 and Table 2, we conducted a meta-analysis and created a forest plot to uncover the association between the clinical efficacy and safety and PD-1 inhibitors. OS, PFS, ORR, and DCR were employed as indicators to assess clinical efficacy. Results revealed that the PD-1 inhibitors group obtained greater OS benefit compared with the standard chemotherapy group (HR = 0.73, 95%CI 0.66–0.81, P < 0.001; I2 = 0%, P = 0.94). PFS was nearly favorable in the PD-1 inhibitor group (HR = 0.89, 0.76–1.04, P = 0.13; I2 = 58%, P = 0.05). We also observed a significant increase in ORR among patients with PD-1 inhibitors (RR = 1.88, 95%CI 1.19–2.98, P < 0.001; I2 = 73%, P < 0.01). However, DCR was not significantly higher in the group treated with PD-1 inhibitors than in the standard chemotherapy group (RR = 0.96, 95%CI 0.74–1.25, P = 0.77; I2 = 84%, P < 0.01).
Table 2
pooled hazard ratio for clinical efficacy and safety.
|
Heterogeneity
|
Effect size
|
P-value
|
model
|
I2
|
P
|
Efficacy
|
|
|
|
|
|
OS
|
0%
|
0.94
|
HR: 0.73 (0.66–0.81)
|
< 0.001
|
common effect
|
PFS
|
58%
|
0.05
|
HR: 0.89 (0.76–1.04)
|
0.129
|
random effect
|
ORR
|
73%
|
< 0.01
|
RR: 1.88 (1.19–2.98)
|
< 0.001
|
random effect
|
DCR
|
84%
|
< 0.01
|
RR: 0.96 (0.74–1.25)
|
0.773
|
random effect
|
Safety
|
|
|
|
|
|
TRAEs
|
96%
|
< 0.01
|
RR: 0.76 (0.64–0.91)
|
< 0.001
|
random effect
|
TRAEs Grade3-5
|
57%
|
0.05
|
RR: 0.40 (0.32–0.49)
|
< 0.001
|
random effect
|
As for safety, TRAEs and the Grade 3–5 TRAEs were selected as indicators to evaluate safety. Of note, Data on TRAEs for separate subgroups of esophageal squamous carcinoma were not available in the KEYNOTE 181 study. Here, we pooled TRAEs data for esophageal squamous and adrenal carcinoma in KEYNOTE 181. Analyses of TRAEs suggested that the PD-1 inhibitors might decrease the rates of overall TRAEs (RR = 0.76, 95%CI 0.64–0.91, P = 0.004; I2 = 96%, P < 0.01) and especially the Grade 3–5 TRAEs (RR = 0.40, 95%CI 0.32–0.49, P < 0.001; I2 = 57.3%, P = 0.05).
Afterward, we performed a sensitivity analysis to ascertain whether the individual studies significantly affect the overall results. The results indicated that individual studies had little influence on the final results and demonstrated that the analysis was relatively robust and reliable (Supplement Fig. 1). Due to the comparatively small number of the included studies, we did not draw a funnel plot to demonstrate publication bias.
3.3 Subgroup analysis to explore factors modifying the clinical efficacy
In the subgroup analysis, no significant difference was observed in subgroups including age, sex, smoking status, race, region,ECOG status, organ metastasis, and lymph node metastasis. However, the results demonstrated the lower level of PD-L1 expression was associated with a more diminished response to PD-1 inhibitors.
The results showed that the overall HR for PD-1 inhibitors versus. Chemotherapy was 0.61 (95%CI 0.48–0.76, P < 0.001; I2 = 0%, P = 0.61) for CPS of PD-L1 ≥ 10 subgroups, and 0.84 (95%CI 0.70-1.00, P < 0.001; I2 = 0%, P = 0.77) for CPS of PD-L1 < 10 subgroups. The HR for interaction effect between CPS of PD-L1 ≥ 10 and PD-L1 < 10 subgroups was 0.72 (95%CI 0.53–0.96, P = 0.02; I2 = 0%, P = 0.56). As regards to the TPS of PD-L1 status, we separately compared the three groups (HR interaction: ≥10% vs. <10%, 0.91; ≥5% vs. <5%, 0.83; ≥1% vs. <1%, 0.80) and ultimately identified no statistically significant differences. Regardless of the degree of PD-L1 expression, as evaluated by CPS or TPS, the PD-1 inhibitor groups all obtained a greater survival advantage compared with chemotherapy groups.
Table 3 preserved the subgroup analysis of overall survival according to different modifying factors. Apart from the PD-L1 status, sex was probably an underlying significant factor modifying the PD-1 inhibitors clinical efficiency. The overall HR for PD-1 inhibitors versus. Chemotherapy was 0.76 (95%CI 0.68–0.84, P < 0.001; I2 = 0%, P = 0.97) for males and 0.62 (95%CI 0.46–0.82, P < 0.001; I2 = 0%, P = 0.57) for females. There was a nearly significant interaction effect between sexes on the relationship between OS and PD-1 inhibitors (Male versus. Female, HR = 1.26, 95%CI 0.92–1.71, P = 0.14; I2 = 0%, P = 0.66).
Table 3
Subgroup analysis of meta-analysis concerning overall survival.
Subgroup
|
Included studies
|
No. of patients
|
HR (95% CI)
|
HR (interaction action, 95% CI)
|
P (subgroups)
|
Sex13–17
|
|
|
|
|
|
Male
|
5
|
1705
|
0.76 (0.68–0.84)
|
1.26 (0.92–1.71)
|
0.14
|
Female
|
5
|
265
|
0.62 (0.46–0.82)
|
Age13–17
|
|
|
|
|
|
<65 years
|
5
|
1215
|
0.73 (0.64–0.82)
|
1.00 (0.81–1.24)
|
0.97
|
≥65 years
|
5
|
755
|
0.73 (0.62–0.87)
|
Race rationale16–17
|
|
|
|
|
|
Asian
|
2
|
|
0.75 (0.64–0.87)
|
1.38 (0.84–2.25)
|
0.20
|
Non-Asian
|
2
|
|
0.53 (0.34–0.84)
|
Region15,17
|
|
|
|
|
|
Asia
|
2
|
635
|
0.70 (0.59–0.83)
|
0.94 (0.48–1.81)
|
0.99
|
Ex-Asia
|
2
|
278
|
0.74 (0.43–1.28)
|
ECOG13–17
|
|
|
|
|
|
0
|
5
|
624
|
0.79 (0.66–0.95)
|
1.12 (0.89–1.40)
|
0.84
|
1
|
5
|
1346
|
0.70 (0.62–0.79)
|
Smoking status13,16,17
|
|
|
|
|
|
Former/Current
|
3
|
863
|
0.75(0.64–0.88)
|
1.03 (0.74–1.44)
|
0.85
|
Never
|
3
|
257
|
0.72(0.54–0.96)
|
Number of organs with metastases14,16
|
|
|
|
|
|
1
|
2
|
359
|
0.81 (0.63–1.03)
|
1.19 (0.87–1.63)
|
0.24
|
>=2
|
2
|
508
|
0.68 (0.56–0.82)
|
Lymph node metastasis14,16
|
|
|
|
|
|
No
|
2
|
405
|
0.68 (0.55–0.86)
|
0.85 (0.63–1.15)
|
0.36
|
Yes
|
2
|
462
|
0.80 (0.65–0.98)
|
PD-L1 status (CPS)13, 15,17
|
|
|
|
|
|
≥ 10
|
3
|
384
|
0.61(0.48–0.76)
|
0.72(0.53–0.96)
|
0.02
|
< 10
|
3
|
565
|
0.84(0.70-1.00)
|
PD-L1 status (TPS)13,14,16
|
|
|
|
|
|
≥ 10%
|
3
|
339
|
0.70 (0.51–0.96)
|
0.91 (0.63–1.31)
|
0.61
|
< 10%
|
3
|
931
|
0.75 (0.64–0.88)
|
≥ 5%
|
2
|
229
|
0.65 (0.48–0.87)
|
0.83 (0.59–1.18)
|
0.30
|
< 5%
|
2
|
628
|
0.76 (0.63–0.91)
|
≥ 1%
|
3
|
448
|
0.66 (0.53–0.81)
|
0.80 (0.56–1.15)
|
0.24
|
< 1%
|
3
|
550
|
0.82(0.68-1.00)
|