ORIGINAL RESEARCH article

Front. Immunol., 01 December 2025

Sec. Cancer Immunity and Immunotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1653218

Overall survival according to time-of-day of combined immuno-chemotherapy for advanced gastric cancer: a propensity score-matched analysis

  • 1. Department of Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China

  • 2. Breast Cancer Center, Hubei Cancer Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

  • 3. National Key Clinical Specialty Discipline Construction Program, Hubei Provincial Clinical Research Center for Breast Cancer, Wuhan Clinical Research Center for Breast Cancer, Wuhan, Hubei, China

Abstract

Purpose:

The impact of time-of-day administration (ToDA) of immune checkpoint inhibitor (ICI) on patient outcomes across multiple cancer types is increasingly being elucidated. Following the results of the CheckMate-649 study, chemotherapy combined with immunotherapy has been established as the standard first-line treatment for advanced gastric cancer (GC). A thorough investigation of the relationship between ICI infusion timing and patient outcomes within this combined regimen holds considerable potential to enhance and optimize clinical treatment strategies.

Methods:

We conducted a retrospective analysis of patients with advanced GC from West China Hospital, Sichuan University, who received first-line chemotherapy combined with ICIs between January 2020 and September 2024. Patients who received fewer than two doses of ICIs were excluded. Follow-up continued until May 2025. The primary endpoint was overall survival (OS), defined as the time from initial ICI infusion to death from any cause. Secondary endpoints included progression-free survival (PFS), objective response rate (ORR), and adverse events (AEs). Propensity score matching (1:2 ratio, caliper width 0.1) mitigated confounding factors. The impact of ICI infusion timing (after 1630h) on OS and PFS was evaluated using Cox proportional hazards regression. ORR and AEs were assessed with either the Chi-square test or Fisher’s exact test. To further assess the robustness of our findings, sensitivity analyses were conducted using infusion proportion thresholds of 30%, 40%, and 50% at the fixed time point of 1630h, along with time-point sensitivity analyses at 30-minute intervals across the 1500h to 1700h period.

Results:

Among 214 patients, 50 patients received ≥20% of ICI after 1630h (Group A), while 164 patients received <20% (Group B). Before propensity score matching, patients in Group B exhibited significantly shorter OS compared to those in Group A (median 15.4 vs. 21.4 months, HR = 1.64, P = 0.014). After matching, the Group B (44 patients) continued to demonstrate significantly shorter OS than Group A (86 patients) (median 15.4 vs. 22.4 months, HR = 1.82, P = 0.013). Multivariable analysis confirmed these findings. No significant differences were found in PFS, ORR, or AEs (all P>0.050). Sensitivity analysis further validated the robustness of the results.

Conclusion:

Early ToDA of ICIs is associated with longer OS in patients with advanced GC receiving first-line chemotherapy combined with ICIs. Randomized clinical trials are required to validate these findings.

1 Introduction

Circadian rhythms constitute endogenous timing systems in organisms, characterized by 24-hour oscillations in physiological processes (). These evolutionarily conserved mechanisms maintain organism-environment homeostasis by coordinating essential functions including sleep-wake cycles, energy metabolism, and immune responses (–). Chronotherapy, also known as circadian medicine, leverages endogenous circadian rhythms through administering treatments during biologically optimal time windows, thereby maximize therapeutic outcomes ().

Previous studies have gradually elucidated the mechanisms through which circadian rhythms influence chemotherapy toxicity and the efficacy of therapeutic cytokines (, ). In recent years, the regulatory relationship between circadian rhythms and immune function has become increasingly well-defined. Research demonstrates that various immune cells, including monocytes, macrophages, neutrophils, and natural killer cells, possess intrinsic circadian clocks and through these regulatory mechanisms play vital roles in both innate and adaptive immune responses (–).Preclinical studies have demonstrated significant circadian fluctuations in the number of peripheral blood lymphocytes in healthy adults, with secretion patterns of various pro-inflammatory cytokines closely linked to cancer progression and treatment response (, 13). CD8+ T cells, which are key targets of immune checkpoint inhibitors (ICIs), have recently been shown to possess intrinsic circadian rhythms, with their immune response intensity demonstrating significant diurnal variation and peaking within specific time windows (). Collectively, these findings underscore that synchronizing immunotherapy with the body’s intrinsic circadian rhythms is critical for maximizing therapeutic efficacy.

A single-center cohort study published in 2021(MEMOIR) was the first to establish a significant correlation between the time-of-day administration (ToDA) of ICIs and patient outcomes (). In melanoma patients, those who received ≥20% of their ICI infusions after 1630h had significantly worse survival compared to those who received <20% of infusions after this time. This finding has since been corroborated in a variety of solid tumors, including lung cancer, head and neck cancer, and renal cell carcinoma (RCC) (–). Real-world data from eighteen retrospective studies indicated that early ToDA of ICIs could improve progression-free survival (PFS) and overall survival (OS) by up to fourfold compared to late ToDA (). However, research on advanced gastric cancer (GC) has been limited to nivolumab in later-line therapies (, ).

Given that combination chemotherapy and ICIs have been incorporated into first-line treatment guidelines for advanced GC, clarifying the impact of ICI infusion timing on survival outcomes with this specific regimen is of substantial clinical significance. The findings are expected to provide evidence-based insights for optimizing clinical dosing schedules.

2 Methods

2.1 Study design and participants

This single-center, retrospective cohort study consecutively enrolled patients with advanced GC who received first-line chemotherapy combined with ICIs at West China Hospital, Sichuan University between January 2020 and September 2024. Patients were included if they had completed at least two cycles of ICI treatment. Follow-up data were collected until May 1, 2025. Clinical and pathological characteristics, ICI infusion timing, and survival outcomes were retrieved from the hospital’s electronic medical records. Survival data were supplemented via telephone follow-up when necessary. The key inclusion criteria were: (1) histologically confirmed gastric cancer; (2) diagnosis of locally advanced or metastatic disease at the initiation of immunotherapy; and (3) completion of at least two cycles of first-line chemotherapy combined with immunotherapy. Patients with incomplete clinical data were excluded (Figure 1). The study protocol was approved by the Ethics Committee of West China Hospital, Sichuan University (Approval No. 20242530).

Figure 1

2.2 Procedures

Treatment plans were determined by a multidisciplinary team. Immune checkpoint inhibitors were administered at standard regimens: Sintilimab 200mg every 3 weeks (q3w), Nivolumab 240mg q2w/360mg q3w, Pembrolizumab 200mg q3w, and Tislelizumab 200mg q3w. Doses of concomitant chemotherapy or targeted therapy were adjusted based on patient tolerance. In accordance with our institutional medication protocols, ICIs were typically administered prior to chemotherapy, unless specific circumstances dictated otherwise. Tumor response was assessed by imaging every 2–3 cycles and evaluated per RECIST version 1.1.

All treatment-related adverse events (AEs) from treatment initiation until the last follow-up were retrospectively reviewed and graded according to the Common Terminology Criteria for Adverse Events (CTCAE) version 5.0. Immune-related adverse events (irAEs) were specifically identified and classified based on the criteria established in the ORIENT-16 trial ().

2.3 Outcomes

The primary endpoint was OS, defined as the time from the initiation of immunotherapy to death from any cause. Secondary endpoints included PFS, objective response rate (ORR), and AEs. PFS was defined as the time from treatment initiation to either disease progression or death from any cause, whichever occurred first. ORR was defined as the proportion of patients achieving a complete or partial response. AEs were assessed from the start of ICI treatment until the last follow-up, graded according to CTCAE version 5.0. The incidence of specific irAEs, such as thyroid dysfunction, myocarditis, and rash, was documented separately using the classification criteria from the ORIENT-16 study ().

2.4 Statistical analysis

Based on previous vaccination studies and retrospective research (, ), patients were stratified into two groups using a cut-off point of 1630h and an infusion threshold of 20%. Group A (Early) comprised patients receiving <20% of their ICI infusions after 1630h, while Group B (Late) comprised those receiving ≥20% of infusions after this time.

To improve group comparability, we performed 1:2 propensity score matching (PSM) using logistic regression (caliper=0.1), with covariates including concurrent targeted drugs, primary tumor location, degree of differentiation, MMR, pre-treatment CA19-9, and CA72-4. Because we found a trend towards association between these factors and OS in univariate analysis (P<0.10). Balance was assessed post-matching. Survival was analyzed using Kaplan-Meier curves and Cox proportional hazards regression. Unadjusted Hazard ratios (HRs) for OS were derived from Univariate Cox proportional hazards regression in all cohorts. Multivariate Cox proportional hazards regression, adjusted for factors with P < 0.10 from univariate analysis, provided adjusted HRs for OS. PFS, ORR, and AEs were analyzed in the unmatched cohort using Chi-square or Fisher’s exact tests.

Sensitivity analyses employed alternative infusion thresholds (30%, 40%, 50%) and time points (1530h, 1600h, 1700h). All statistical analyses were performed using R software (version 4.0.3). A two-sided significance level of P < 0.05 was applied for all statistical tests.

3 Results

3.1 Patient characteristics and study flow

A total of 246 patients were initially screened. After excluding 24 patients who received fewer than two cycles of ICI and 8 with incomplete data, 214 were included in the final analysis. Based on a 20% infusion threshold after 1630h, patients were stratified into Group B (Late, n=50) and Group A (Early, n=164).

The cohort consisted of 148 patients (69%) aged <65 years and 66 (31%) aged ≥65 years, with a male predominance (144 [67%] vs. 70 [33%] female). The most commonly used ICIs were sintilimab (75%), nivolumab (13%), an additional 12% received tislelizumab, pembrolizumab, toripalimab or camrelizumab. All patients received ICIs combined with chemotherapy, primarily FOLFOX, SOX, or CAPOX regimens. The majority (194, 91%) received doublet chemotherapy, while 17 (8%) received multi-drug regimens, and 3 (1%) received monotherapy due to poor performance status. Primary tumors were located in the stomach or gastroesophageal junction, and the predominant histology was adenocarcinoma (94%), with a few cases of mixed or other types (e.g., neuroendocrine carcinoma, poorly differentiated carcinoma). 93% of patients had poorly differentiated tumors, while the remaining 7% exhibited moderately differentiated histology.

Univariate Cox regression was performed to identify prognostic factors. Following 1:2 propensity score matching, a balanced cohort of 130 patients was generated, comprising 44 in Group B and 86 in Group A. The baseline characteristics of both the unmatched and matched cohorts are presented in Table 1.

Table 1

VariableUnmatched populationP valueMatched populationP value
Received ≥20% of infusions after 1630 h (n=50)Received <20% of infusions after 1630 h (n=164)Received ≥20% of infusions after 1630 h (n=44)Received <20% of infusions after 1630 h (n=86)
Age0.041*0.151
<6541 (82%)107 (65%)37 (84%)61 (71%)
≥659 (18%)57 (35%)7 (16%)25 (29%)
Sex0.6910.431
Male32 (64%)112 (68%)28 (64%)62 (72%)
Female18 (36%)52 (32%)16 (36%)24 (28%)
ECOG performance status0.3910.331
038 (76%)112 (68%)35 (80%)60 (70%)
≥112 (24%)52 (32%)9 (20%)26 (30%)
Disease stage0.5320.722
Locally advanced or recurrent2 (4%)12 (7%)2 (5%)7 (8%)
Metastatic48 (96%)152 (93%)42 (95%)79 (92%)
ICI0.3710.502
Nivolumab4 (8%)24 (15%)3 (7%)11 (13%)
Sintilimab41 (82%)119 (73%)36 (82%)62 (72%)
Others#5 (10%)21 (13%)5 (11%)13 (15%)
Concurrent targeted drugs0.7611.001
Yes13 (26%)37 (23%)12 (27%)23 (27%)
No37 (74%)127 (77%)32 (73%)63 (73%)
Chemotherapy regimen administered0.1320.222
Single agent0 (0%)3 (2%)0 (0%)2 (2%)
Dual therapy49 (98%)145 (88%)43 (98%)75 (87%)
Polychemotherapy1 (2%)16 (10%)1 (2%)9 (10%)
Primary tumor location1.0011.001
Gastric43 (86%)140 (85%)37 (85%)72 (84%)
Gastroesophageal junction7 (14%)24 (15%)7 (16%)14 (16%)
Histological type0.4021.002
Adenocarcinoma46 (92%)155 (95%)41 (93%)79 (92%)
Mixed2 (4%)7 (4%)2 (5%)5 (6%)
Others‡2 (4%)2 (1%)1 (2%)2 (2%)
Degree of differentiation0.2221.002
Poorly-differentiated44 (88%)154 (94%)43 (98%)85 (99%)
Moderately-differentiated6 (12%)10 (6%)1 (2%)1 (1%)
Prior surgery0.6810.771
Yes14 (28%)39 (24%)12 (27%)20 (23%)
No36 (72%)125 (76%)32 (73%)66 (77%)
Liver metastasis0.4110.881
Yes13 (26%)55 (34%)12 (27%)26 (30%)
No37 (74%)109 (66%)32 (73%)60 (70%)
Peritoneal metastasis0.021*0.161
Yes32 (64%)72 (44%)27 (61%)40 (47%)
No18 (36%)92 (56%)17 (39%)46 (53%)
Metastasis to other sites0.3310.151
Yes14 (28%)33 (20%)12 (27%)13 (15%)
No36 (72%)131 (80%)32 (73%)73 (85%)
Number of metastatic organs0.08510.0821
07 (14%)33 (20%)7 (16%)19 (22%)
127 (54%)102 (62%)23 (52%)54 (63%)
≥216 (32%)29 (18%)14 (32%)13 (15%)
MMR status0.3621.002
pMMR47 (94%)160 (98%)43 (98%)85 (99%)
dMMR3 (6%)4 (2%)1 (2%)1 (1%)
Her20.7720.732
Positive5 (10%)13 (8%)4 (9%)6 (7%)
Negative45 (90%)151 (92%)40 (91%)80 (93%)
CEA0.9311.001
≥5ng/mL19 (38%)59 (36%)16 (36%)31 (36%)
<5ng/mL31 (62%)105 (64%)28 (64%)55 (64%)
CA1250.2810.861
≥24U/mL30 (60%)82 (50%)25 (57%)46 (53%)
<24U/mL20 (40%)82 (50%)19 (43%)40 (47%)
CA19-90.5211.001
≥30 U/mL14 (28%)56 (34%)14 (32%)26 (30%)
< 30U/mL36 (72%)108 (66%)30 (68%)60 (70%)
CA72-41.0011.001
≥6.9U/mL31 (62%)100 (61%)26 (59%)51 (59%)
<6.9U/mL19 (38%)64 (39%)18 (41%)35 (41%)
AFP0.5210.631
≥7ng/mL6 (12%)28 (17%)5 (11%)14 (16%)
<7ng/mL44 (88%)136 (83%)39 (89%)72 (84%)
NLR0.1010.161
≥332 (64%)81 (49%)27 (61%)40 (47%)
<318 (36%)83 (51%)17 (39%)46 (53%)
PLR0.4210.151
≥20016 (32%)65 (40%)11 (25%)34 (40%)
<20034 (68%)99 (60%)33 (75%)52 (60%)
Charlson Comorbidity Index0.4110.461
041 (82%)123 (75%)37 (84%)66 (77%)
≥19 (18%)41 (25%)7 (16%)20 (23%)
number of concomitant medication types†0.0710.161
045 (90%)126 (77%)39 (89%)66 (77%)
≥15 (10%)38 (23%)5 (11%)20 (23%)
Antibiotic therapy1.0021.002
Yes1 (2%)6 (4%)1 (2%)2 (2%)
No49 (98%)158 (96%)43 (98%)84 (98%)

Patient characteristics in the unmatched and propensity score-matched population.

ICI, Immune checkpoint inhibitor; ECOG, Eastern Cooperative Oncology Group; AFP, Alpha Fetoprotein; CA19-9, Carbohydrate Antigen 19-9; CEA, Carcinoembryonic Antigen; CA125, Carbohydrate Antigen 125; CA72-4, Carbohydrate Antigen 72-4; NLR, Neutrophil-to-Lymphocyte Ratio; PLR, Platelet-to-Lymphocyte Ratio; MMR, Mismatch Repair.

The others# for ICI included tislelizumab, pembrolizumab, toripalimab and camrelizumab. The others‡ for histological type included neuroendocrine carcinoma and poorly differentiated carcinoma. Number of concomitant medication types† included antibiotics, anticoagulants, analgesics, and hypnotics.

1Chi-square test; 2Fisher’s Exact Test. *P ≤ 0.05

3.2 Clinical outcome

Before matching, the median follow-up time for patients in Group A was 14.6 months, while for those in Group B, it was 15.1 months. Kaplan-Meier analysis showed that Group B had significantly shorter OS than Group A (median 15.4 vs. 21.4 months, HR = 1.64, P = 0.014; Figure 2A). After adjusting for factors such as concurrent targeted drugs, primary tumor location, degree of differentiation, MMR status, pre-treatment CA19–9 and CA72-4, the adjusted HR was 1.58 (95% CI: 1.05-2.39, P = 0.028; Figure 2A). Multivariable cox proportional hazards regression indicated that concurrent targeted drugs was an independent prognostic factor (Table 2).

Figure 2

Table 2

VariableHR (95% CI) univariableP valueHR (95% CI) multivariableP value
Timing1630 (Reference= infusion<20%)
Received ≥20% of infusions after 1630h1.64 (1.11-2.44)0.0141.58 (1.05-2.39)0.028*
Age (Reference= <65)
≥650.99 (0.67-1.47)0.978
Sex (Reference= Male)
Female0.87 (0.60-1.28)0.483
ECOG (Reference= 0)
≥11.12 (0.76-1.65)0.554
Disease stages (Reference=Locally recurrent/advanced)
Metastatic1.39 (0.61-3.16)0.436
ICI (Reference= Sintilimab)
Nivolumab0.74(0.46-1.19)0.217
Others#0.67 (0.34-1.30)0.236
Concurrent targeted drugs (Reference=No)
Yes0.40 (0.23-0.70)0.0010.41 (0.24-0.72)0.002*
Chemotherapy regimen administered (Reference= Single agent)
Dual therapy0.65 (0.16-2.65)0.549
Polychemotherapy1.20 (0.27-5.30)0.810
Previous radical resection (Reference=No)
Yes0.96 (0.64-1.43)0.823
Primary tumor location (Reference= gastric)
Gastroesophageal junction0.62 (0.36-1.09)0.0950.67 (0.38-1.19)0.173
Histological type (Reference= Adenocarcinoma)
Mixed-type0.66 (0.24-1.79)0.416
Others‡0.89 (0.22-3.60)0.867
Degree of differentiation (Reference= Poorly-differentiate)
Moderately-differentiated1.96 (1.07-3.56)0.0281.74 (0.93-3.27)0.083
Liver metastasis (Reference=No)
Yes0.80 (0.55-1.18)0.261
Peritoneal metastasis (Reference=No)
Yes1.19 (0.84-1.70)0.331
Metastasis to other sites (Reference=No)
Yes1.38 (0.93-2.06)0.114
Number of metastatic organs (Reference=0)
10.79 (0.49-1.25)0.311
≥21.25 (0.74-2.11)0.403
MMR (Reference= pMMR)
dMMR2.23 (0.98-5.09)0.0571.49 (0.63-3.57)0.366
Her2 (Reference= Negative)
Positive0.73 (0.35-1.49)0.381
CEA (Reference=<5ng/mL)
≥5ng/mL1.33 (0.93-1.92)0.118
CA125 (Reference=<24U/mL)
≥24U/mL1.35 (0.94-1.93)0.100
CA19-9(Reference=<30U/mL)
≥30U/mL1.46 (1.02-2.09)0.0401.38 (0.95-2.00)0.095
CA72-4(Reference= <6.9U/mL)
≥6.9U/mL1.38 (0.95-1.99)0.0861.24 (0.84-1.82)0.284
AFP (Reference= <7ng/mL)
≥7ng/mL1.01 (0.62-1.64)0.979
NLR (Reference= NLR<3)
NLR≥31.12 (0.79-1.60)0.521
PLR (Reference= PLR<200)
PLR≥2000.80 (0.55-1.17)0.245
Charlson Comorbidity Index (Reference= 0)
≥11.04 (0.68-1.59)0.847
number of concomitant medication types† (Reference= 0)
≥11.23 (0.80-1.89)0.337
Antibiotic therapy (Reference= No)
Yes0.96 (0.35-2.61)0.935

Univariate and multivariable cox proportional hazards regression of overall survival in the unmatched population.

ICI, Immune checkpoint inhibitor; ECOG, Eastern Cooperative Oncology Group; AFP, Alpha Fetoprotein; CA19-9, Carbohydrate Antigen 19-9; CEA, Carcinoembryonic Antigen; CA125, Carbohydrate Antigen 125; CA72-4, Carbohydrate Antigen 72-4; NLR, Neutrophil-to-Lymphocyte Ratio; PLR, Platelet-to-Lymphocyte Ratio; MMR, Mismatch Repair.

The others# for ICI included tislelizumab, pembrolizumab, toripalimab and camrelizumab. The others‡ for histological type included neuroendocrine carcinoma and poorly differentiated carcinoma. Number of concomitant medication types† included antibiotics, anticoagulants, analgesics, and hypnotics.

Variables with P<0.1 were included in the multivariable cox proportional hazards regression. The adjustment factors include concurrent targeted drugs, primary tumor location, degree of differentiation, MMR, pre-treatment CA19-9, and pre-treatment CA72-4.

*P<0.05 in multivariable cox proportional hazards regression.

After matching, Group B again showed worse OS (median 15.4 vs. 22.5 months, HR = 1.82, P = 0.013; Figure 2B). Adjusting for the factors include concurrent targeted drugs, primary tumor location, liver metastasis, number of metastatic organs, pre-treatment CEA and CA19-9, subsequent multivariable analysis of the matched cohort confirmed the robustness of this finding (HR:1.65, 95% CI: 1.01-2.70, P = 0.044; Figure 2B) and revealed that concurrent targeted drugs was independent prognostic factors (Table 3).

Table 3

VariableHR (95% CI) univariableP valueHR (95% CI) multivariableP value
Timing1630 (Reference= infusions<20%)
Received ≥20% infusions after 1630 h1.82 (1.14-2.90)0.0131.65 (1.01-2.70)0.044*
Concurrent targeted drugs (Reference=No)
Yes0.42 (0.22-0.80)0.0080.44 (0.21-0.91)0.027*
Primary tumor location (Reference= gastric)
Gastroesophageal junction0.51 (0.26-1.03)0.0610.51 (0.25-1.03)0.060
Liver metastasis (Reference= No)
Yes0.55 (0.33-0.93)0.0260.68 (0.36-1.30)0.244
Member of metastatic organs (Reference=0)
10.55 (0.32-0.97)0.0400.74 (0.41-1.35)0.325
≥20.65 (0.33-1.25)0.193
CEA(Reference=<5ng/mL)
≥5ng/mL1.50 (0.95-2.39)0.0841.40 (0.80-2.46)0.240
CA19-9 (Reference= <30U/mL)
≥30U/mL1.83 (1.16-2.88)0.0101.57 (0.88-2.80)0.129

Univariate and multivariable cox proportional hazards regression of overall survival in the matched population.

CA19-9, Carbohydrate Antigen 19-9; CA72-4, Carbohydrate Antigen 72-4; MMR, Mismatch Repair.

Variables with P<0.1 were included in the multivariable cox proportional hazards regression. The adjustment factors include concurrent targeted drugs, primary tumor location, liver metastasis, member of metastatic organs, pre-treatment CEA, pre-treatment CA19-9.

*P<0.05 in multivariable cox proportional hazards regression.

No significant between-group differences were found in PFS (HR = 1.33, P = 0.190; Figure 3), ORR (P = 0.653; Table 4), nor the incidence of any category of AEs (all P > 0.050; Table 5). The most frequent irAEs were thyroid dysfunction (28.5%), elevated thyroid stimulating hormone (6.5%), and myocarditis (3.7%).

Figure 3

Table 4

GroupORR, n (%)χ2P value
YesNo
Received ≥20% infusions after 1630 h18 (36%)32 (64%)0.2020.6531
Received <20% infusions after 1630 h67 (41%)97 (59%)

ORR comparison in the unmatched population.

1Chi-square test; 2Fisher’s Exact Test.

Table 5

GroupUnmatched populationχ2P value
Received ≥20% infusions after 1630 h (n=50)Received <20% infusions after 1630 h (n=164)
AE45 (90%)149 (91%)0.0001.0001
Grade 3/4 AE14 (28%)49 (30%)0.0060.9381
irAE22 (44%)61 (37%)0.4880.4851
Grade 3/4 irAE0 (0)2 (1%)/1.0002
AE leading to discontinuation of medication0 (0)6 (4%)/0.3402

Adverse events comparison in the unmatched population.

AE, Adverse event; irAE, Immune-related adverse event.

1Chi-square test; 2Fisher’s Exact Test.

3.3 Sensitivity analysis

3.3.1 Sensitivity analysis of different proportion thresholds (1630h)

To validate the robustness of our findings, we further evaluated the impact of different infusion proportion thresholds (30%, 40%, and 50%) at 1630h on clinical outcomes.

In the pre-matched cohort, all three thresholds were significantly associated with OS:30% (HR:2.15, 95% CI: 1.39–3.33, P < 0.001), 40% (HR: 2.01, 95% CI: 1.24–3.26, P = 0.004), and 50% thresholds (HR: 2.37, 95% CI: 1.43–3.92, P < 0.001) (Supplementary Table S1). These associations remained significant via multivariate adjustment, with adjusted HRs of 2.08 (95% CI: 1.33–3.25, P = 0.001), 1.93 (95% CI: 1.19–3.15, P = 0.008), and 2.30 (95% CI: 1.38–3.82, P = 0.001), respectively (Supplementary Table S1).

In the propensity score-matched cohort, significant associations with OS were similarly observed at the 30% (HR: 1.99, 95% CI: 1.17–3.40, P = 0.011), 40% (HR: 1.94, 95% CI: 1.05–3.59, P = 0.035), and 50% thresholds (HR: 2.16, 95% CI: 1.15–4.07, P = 0.017) (Supplementary Table S2). After multivariable adjustment, the corresponding HRs were 2.17 (95% CI: 1.22–3.85, P = 0.009), 2.26 (95% CI: 1.18–4.31, P = 0.014), and 2.61 (95% CI: 1.34–5.08, P = 0.005), respectively (Supplementary Table S2).

In the unmatched cohort, significant correlations with PFS were observed at the 30% (HR: 1.66, 95% CI: 1.01–2.73, P = 0.045), 40% (HR: 1.88, 95% CI: 1.11–3.20, P = 0.020), and 50% thresholds (HR: 2.48, 95% CI: 1.41–4.36, P = 0.002) (Supplementary Table S3). These remained significant after multivariable adjustment (all P < 0.05; Supplementary Table S3). However, no significant differences in ORR or irAEs were observed at any threshold (all P > 0.05; Supplementary Tables S4, S5).

3.3.2 Sensitivity analysis of different time cut-off points

We further analyzed the impact of ICI infusion timing using 1530h, 1600h, and 1700h as alternative cut-off points.

In the unmatched sample, OS showed no significant association with infusion timing at 1530h. However, significant associations were observed at 1600h (HR;1.49, 95% CI: 1.03-2.15, P = 0.035) and 1700h (HR: 1.99, 95% CI: 1.24-3.20, P = 0.005) (Supplementary Table S6). OS was significantly associated infusion timing. These findings were confirmed in multivariable Cox regression analyses for the unmatched cohort, with adjusted HRs of 1.63 (95% CI: 1.11-2.37, P = 0.012) at 1600h and 1.79 (95% CI: 1.06-3.01, P = 0.029) at 1700h (Supplementary Table S6).

In the propensity score-matched cohort, infusion timing was significantly associated with OS at 1530h (HR: 1.74, 95% CI: 1.17-2.59, P = 0.007), 1600h (HR: 1.80, 95% CI: 1.16-2.79 P = 0.009), and 1700h (HR: 1.91, 95% CI: 1.01-3.62, P = 0.047) (Supplementary Table S7). Subsequent multivariable analysis of the matched cohort yielded adjusted hazard ratios of 1.72 (95% CI: 1.13-2.61; P = 0.012) for 1530h, 1.77(95% CI: 1.09-2.86, P = 0.021) for 1600h, and 1.88 (95% CI: 0.99-3.57, P = 0.053) for 1700h (Supplementary Table S7).

In the unmatched cohort, only the 1700h cut-off was significantly associated with PFS (HR: 1.82, 95% CI: 1.08-3.06, P = 0.024), and this persisted after multivariable adjustment (HR: 1.81, 95% CI: 1.05-3.09, P = 0.031; Supplementary Table S8). No significant differences were observed for 1530h and 1600h. No differences in ORR or irAEs were observed at any cut-off point (all P > 0.050; Supplementary Tables S9, S10).

4 Discussion

This retrospective study demonstrates a significant association between the timing of ICI administration and survival outcomes in patients with advanced GC receiving first-line immunotherapy combined with chemotherapy. The primary analysis revealed that in the unmatched cohort, patients who received ≥20% of their ICI infusions after 1630h had significantly shorter overall survival. This finding remained consistent after controlling for potential confounding factors through propensity score matching.

Sensitivity analyses further confirmed the robustness of our findings. When using 1630h as the reference time point with different thresholds (30%, 40%, and 50%), all thresholds demonstrated significant associations with both OS and PFS in the pre-matched cohort. These associations remained statistically significant in the matched cohort. Furthermore, analyses using alternative time cut-offs (1530h, 1600h, and 1700h) revealed significant associations with OS in the matched cohort. It should be noted that for the 1700h cut-off, the substantial sample size disparity in the pre-matched cohort and the subsequent reduction in effective sample size after 1:2 propensity score matching warrant cautious interpretation. The generalizability of these particular findings requires validation in prospective studies.

The biological plausibility of our findings is supported by a growing body of evidence on circadian immunology. Previous studies have demonstrated that the timing of vaccination could influence the magnitude of antibody responses; for instance, morning administration of influenza or COVID-19 vaccines has been associated with a more robust antiviral antibody response (, ). The circadian clock regulates innate immunity by controlling the rhythmic production and transport of cytokines and chemokines, as well as the maturation and tissue infiltration of immune cells. Within 24 hours following antibody infusion, tumor-infiltrating CD8+ T cells exhibit significant time-dependent phenotypic alterations, indicating their participation in the primary circadian response. It is noteworthy that the antitumor effect was significantly higher in the morning compared to the afternoon (). Studies in immunodeficient mice have shown that dendritic cells (DCs) rhythmically migrate to tumor-draining lymph nodes, thereby orchestrating the circadian rhythm of tumor antigen-specific CD8+ T cell responses (). This mechanism indicates that synchronizing immunotherapy with the peak functional phase of DCs could significantly enhance its antitumor efficacy. Conversely, disruption of the epithelial cell circadian clock alters cytokine secretion, exacerbates inflammatory responses, promotes neutrophil infiltration, and can ultimately lead to the expansion of myeloid-derived suppressor cells (MDSCs) (). Importantly, the therapeutic efficacy of PD-L1 inhibitors appears to be optimal when administered in alignment with the circadian fluctuations of MDSCs (). Consequently, strategically timing immune system activation based on circadian principles may not only maximize the intensity of immune responses but also optimize antitumor treatment outcomes.

Two previous retrospective studies investigating nivolumab as a later-line treatment for GC both showed superior survival outcomes for patients receiving morning infusions (, ). Following the publication of the CheckMate649 trial, the era of first-line chemotherapy combined with immunotherapy became the standard of care for advanced GC (). This regimen has improved mOS to approximately 15 months, compared to about 10 months with chemotherapy alone. Subsequent trials, including ORIENT-16, RATIONALE-305, and GEMSTONE-303, have further established the efficacy of Sintilimab, Tislelizumab, and Sugemalimab in first-line setting (, , ). Given the prevalent use of sintilimab in our clinical practice and study cohort, we performed a subgroup analysis focusing on patients receiving sintilimab plus chemotherapy. In the unmatched sintilimab subgroup, a significant association between infusion timing and OS was observed, which persisted after multivariable adjustment. However, in the propensity score-matched sintilimab subgroup, univariate analysis indicated only a borderline significant association (HR:1.53, 95% CI: 0.94-2.49, P = 0.088), which lost statistical significance after multivariable adjustment (HR:1.50, 95% CI: 0.92-2.45, P = 0.101).This discrepancy may be attributed to several factors: First, while propensity score matching balanced known prognostic variables, there may be unmeasured confounding variables, unmeasured confounders including tumor microenvironment heterogeneity, genomic profiles, or lifestyle differences may persist. The influence of such residual confounding can be magnified in subgroup analyses, potentially obscuring a true treatment effect. Second, the reduced sample size post-matching limited statistical power, increasing the risk of a Type II error. Notably, the hazard ratios from both univariate and multivariate analyses consistently ranged between 1.50 and 1.53, suggesting a stable, albeit not statistically significant, effect size. Although the P-values slightly exceeded the conventional threshold, this consistent trend may still hold clinical relevance, particularly in light of established circadian regulation of immune function.

Previous studies have reported circadian rhythm and gender differences in immune therapy responses, with female patients showing better prognoses than males (, , ). For instance, studies of hepatitis A and influenza vaccines have shown that female individuals often mount stronger peak antibody responses (). This heightened response may be partly attributed to sex-specific resilience to immune dysfunction induced by sleep disturbances and psychological stress (43). Additionally, women tend to display more pronounced circadian rhythmicity, characterized by higher plasma levels of melatonin and cortisol (44, 45). In contrast, a separate systematic review and meta-analysis demonstrated that while ICIs improve overall survival in patients with advanced cancers like melanoma and non-small cell lung cancer, the survival benefit was more pronounced in male patients (46). In our study, we found no significant association between patient sex and survival outcomes. Further research is warranted to clarify the potential role of sex as a biological variable and its interaction with circadian timing in shaping immunotherapy efficacy.

Our study has several limitations. First, although we observed a significant benefit in OS, no between-group differences were observed in PFS or ORR. This discrepancy may be explained by several factors. Patients in the control group might have received subsequent immunotherapy or other effective treatments after disease progression. Such crossover could mitigate between-group differences in PFS while concurrently improving long-term survival in the control group. Additionally, pseudoprogression or delayed clinical responses may not be accurately captured by conventional RECIST criteria, thereby influencing PFS and ORR assessments. The relatively limited sample size, coupled with closely matched median PFS between groups, may also have reduced the statistical power to detect a significant difference in PFS. Second, thresholds for defining early and late ICI infusion vary across previous studies, with cutoff times ranging from 1137h to 1630h and infusion proportion thresholds from 20% to 75% (, ). The 1630h cut-off adopted in our primary analysis is the most commonly applied, potentially due to its alignment with findings from earlier vaccination studies. To better identify optimal daily infusion timing, Catozzi et al. proposed a circadian mortality risk model; however, this model was derived from a heterogeneous cohort encompassing multiple cancer types and treatment lines (). While our sensitivity analyses investigated various infusion cut-offs, they were insufficient to pinpoint a single optimal administration time. This highlights the need for future research specifically designed to identify the most therapeutic time window. Previous studies often determined assessment timing based on initial imaging or overall patient condition and typically excluded patients receiving fewer than four ICI doses (47, 48). In our study, considering standard efficacy assessment practices in China and the generally good performance status (ECOG score ≤1) of enrolled patients, we included all patients who received at least two infusions. The persistence of therapeutic antibody effects over weeks suggests that circadian regulation during the initial ICI cycles may be clinically more impactful than in later cycles. Thus, our study design is particularly relevant for evaluating the effect of infusion timing during the early phase of immunotherapy. Including patients with fewer infusion cycles also improves the real-world applicability of our findings, as early treatment discontinuation due to toxicity or progression is common in clinical practice. We nevertheless acknowledge that this approach introduces limitations, including relatively limited treatment exposure and potential heterogeneity in immune responses, which may affect the robustness of the results.

Adjusting the administration timing of ICIs to the biologically favorable morning window may not only enhance antitumor efficacy but also improve healthcare resource allocation and patient throughput. In accordance with our institutional treatment pathway, an “ICI-first” sequence is generally adopted, wherein immunotherapy is administered prior to chemotherapy. This approach aims to fully activate the immune system before it is subjected to chemotherapy-induced suppression, thereby potentially improving the body’s tolerance to and recovery from the immunotoxic effects of subsequent chemotherapy. However, uniformly concentrating all ICI infusions in the morning hours presents two major challenges. First, it fails to account for individual differences in circadian rhythm phenotypes, which may diminish the potential benefits of chronotherapy for some patients. Second, such highly consolidated scheduling can lead to overutilization of infusion center capacity, nursing staff, and pharmacy resources during the morning, while resulting in underutilization in the afternoon, ultimately undermining sustained gains in operational efficiency. Therefore, future implementation of chronotherapy-based ICI administration requires careful consideration of both individual chronobiological profiles and institutional resource management to optimize therapeutic outcomes while maintaining healthcare system sustainability.

5 Conclusions

In conclusion, our study demonstrates a significant link between earlier ICI infusion timing and improved survival in advanced gastric cancer patients treated with first-line chemoimmunotherapy. Confirming this association in prospective trials and elucidating the specific circadian mechanisms involved are critical next steps.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by The Ethics Committee of West China Hospital, Sichuan University. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because due to the retrospective nature of this study and the fact that all analyses were conducted after de-identification of patient data, informed consent was waived.

Author contributions

CC: Data curation, Formal analysis, Investigation, Writing – original draft, Writing – review & editing. XQC: Formal analysis, Methodology, Supervision, Visualization, Writing – review & editing. TJW: Data curation, Funding acquisition, Methodology, Validation, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1653218/full#supplementary-material

Abbreviations

ToDA, Time-of-day administration; GC, Gastric cancer; OS, Overall survival; PFS, Progression-free survival; ORR, Objective response rate; AEs, Adverse events; ICI, Immune checkpoint inhibitor; RCC, Renal cell carcinoma; CTCAE, Common Terminology Criteria for Adverse Events; irAEs, Immune-related adverse events; HRs, Hazard ratios; DCs, Dendritic cells; MDSCs, Myeloid-derived suppressor cells; ECOG, Eastern Cooperative Oncology Group; AFP, Alpha Fetoprotein; CA19-9, Carbohydrate Antigen 19-9; CEA, Carcinoembryonic Antigen; CA125, Carbohydrate Antigen 125; CA72-4, Carbohydrate Antigen 72-4; NLR, Neutrophil-to-Lymphocyte Ratio; PLR, Platelet-to-Lymphocyte Ratio; MMR, Mismatch Repair.

References

Summary

Keywords

immune checkpoint inhibitors, gastric cancer, chronotherapy, circadian rhythm, immunotherapy

Citation

Cheng C, Cheng X and Wang T (2025) Overall survival according to time-of-day of combined immuno-chemotherapy for advanced gastric cancer: a propensity score-matched analysis. Front. Immunol. 16:1653218. doi: 10.3389/fimmu.2025.1653218

Received

24 June 2025

Revised

07 November 2025

Accepted

12 November 2025

Published

01 December 2025

Volume

16 - 2025

Edited by

Jeffrey J. Pu, Upstate Medical University, United States

Reviewed by

Danyi Lu, Guangzhou University of Chinese Medicine, China

Pasquale F. Innominato, University of Warwick, United Kingdom

Seline Ismail-Sutton, Betsi Cadwaladr University Health Board, United Kingdom

Updates

Copyright

*Correspondence: Tiejun Wang,

†

ORCID: Tiejun Wang, orcid.org/0000-0003-1244-0794

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics