Impact of the immuno-oncology era on cerebrovascular disease and cancer-specific mortality in older adults with advanced non-small cell lung cancer: a competing-risks cohort study
Highlight box
Key findings
• Post-immuno-oncology era (Post-I-O) vs. pre-immuno-oncology era (Pre-I-O): non-small cell lung cancer (NSCLC)-specific death sub-distribution hazard ratio (sHR) 0.624 (↓37.6%), cerebrovascular disease (CVD) death sHR 1.424 (↑42.4%).
What is known and what is new?
• Immune checkpoint inhibitors (ICIs) improve NSCLC survival; older adults have high vascular risk.
• First population-level competing-risks evidence quantifying a significant rise in CVD mortality alongside cancer survival gains in patients ≥70 years.
What is the implication, and what should change now?
• A "success paradox"—better cancer outcomes unmask growing cerebrovascular competing risk.
• Integrate routine vascular risk assessment and interdisciplinary survivorship care (oncology + neurology/cardiology) for older NSCLC patients starting immunotherapy.
Introduction
Non-small cell lung cancer (NSCLC) remains a leading cause of cancer-related mortality worldwide, and outcomes are particularly poor among older adults with advanced-stage disease (1,2). Since 2015, the clinical introduction of immune checkpoint inhibitors (ICIs) has fundamentally altered the treatment landscape for advanced NSCLC (3). In the past decade, ICIs have fundamentally reshaped the treatment landscape of advanced NSCLC, improving survival in multiple clinical settings (4). Moreover, clinical trials have demonstrated the efficacy of these agents in older adults; however, real-world outcomes in the “oldest-old” population remain complex due to the presence of significant comorbidities (5). The survival gains observed in routine practice may be accompanied by changes in non-cancer mortality patterns, especially among older populations who carry a higher baseline burden of vascular disease and frailty.
Cerebrovascular disease (CVD), including ischemic stroke and intracranial hemorrhage, is a leading cause of death and disability in the elderly (6,7). In addition to traditional vascular risk factors, contemporary cancer therapies may influence thrombo-inflammation, endothelial function, and immune-mediated vascular injury (8,9). In patients with advanced cancer, CVD serves as a formidable “competing risk”—a phenomenon where a patient may die from a non-cancer cause before the cancer-specific outcome occurs. The prolongation of life facilitated by ICIs may inadvertently increase the cumulative incidence of cerebrovascular events by extending the “at-risk” window or through potential immune-related inflammatory effects on the vasculature (10,11). As a result, an “improved cancer survival but increased competing mortality” scenario is plausible in the immuno-oncology (I-O) era: patients may live longer and thus accumulate more time at risk for non-cancer events, while ICI-related inflammatory or prothrombotic mechanisms may further increase vascular risk.
Despite increasing attention to cardiotoxicity in oncology, population-level evidence describing how the introduction of ICIs has reshaped cause-specific mortality—particularly cerebrovascular mortality—among older adults with advanced NSCLC remains limited. Importantly, standard survival models may be inadequate when competing events (e.g., NSCLC death) preclude observation of CVD death. Therefore, competing-risks methods are essential to quantify the real-world cumulative burden of CVD death since the I-O era.
In this study, we used the Surveillance, Epidemiology, and End Results (SEER) database to evaluate whether the introduction of ICI therapy corresponds to changes in CVD mortality and NSCLC-specific mortality among adults aged ≥70 years diagnosed with advanced NSCLC. We compared outcomes between the pre-immuno-oncology era (Pre-I-O; 2007–2014) and a post-immuno-oncology era (Post-I-O; 2015–2022) using propensity score matching (PSM), Fine-Gray competing-risks models, cause-specific Cox regression, and multistate modeling, with subgroup analyses to assess consistency across clinically relevant strata. We present this article in accordance with the STROBE reporting checklist (available at https://tcr.amegroups.com/article/view/10.21037/tcr-2026-1-0169/rc).
Methods
Data source and study design
We conducted a retrospective cohort study using the SEER-Medicare linked database. SEER provides population-based cancer registry data, including tumor characteristics and cause of death, and Medicare claims provide information on demographics and treatments (12). The study was designed to compare mortality outcomes before and after the introduction and widespread adoption of I-O therapies for advanced NSCLC in routine care. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments.
Study population
We identified patients who: (I) were aged ≥70 years at diagnosis; (II) had advanced NSCLC at diagnosis (as defined in your analytic datasets using SEER staging variables); (III) were diagnosed between January 1, 2007 and December 31, 2022.
Patients were categorized by diagnosis period into Pre-I-O (January 2007 to December 2014) and Post-I-O (January 2015 to December 2022).
The index date was the date of NSCLC diagnosis. Follow-up started at diagnosis and continued until death or end of follow-up.
Exposure definition
The primary exposure was the I-O era, operationalized as diagnosis period: Pre-I-O vs. Post-I-O. This approach captures system-level changes in treatment availability and clinical adoption associated with ICIs, rather than requiring confirmed ICI receipt in claims (which can be incomplete or influenced by coding changes over time).
Outcomes
The primary mortality outcomes were CVD death, NSCLC-specific death and other-cause death (as a competing event category). Cause of death was determined using SEER cause-of-death fields and categorized into CVD, NSCLC, or other causes according to your predefined coding rules.
Covariates
We adjusted for the following covariates: (I) demographics and social factors: age, sex, race, marital status, rural urban distribution, median household income; (II) tumor characteristics: tumor grade, tumor (T) stage, node (N) stage, metastasis (M) stage (all strictly analyzed as categorical variables); and (III) treatment variables: surgery recode, radiotherapy recode, chemotherapy recode.
PSM
To reduce confounding between eras, we estimated each patient’s propensity to be diagnosed in the Post-I-O using the covariates above. We then performed 1:1 PSM to generate balanced cohorts for the primary analyses. Covariate balance was assessed using standardized differences and/or distributional checks (as implemented in your analysis).
Statistical analysis
Cumulative incidence (absolute risk)
We estimated and visualized the cumulative incidence functions (CIFs) for CVD death and NSCLC death, treating other causes of death as competing events. Differences in CIF patterns were interpreted as changes in the population-level absolute burden of each cause of death across eras.
Primary model: Fine-Gray competing-risks regression
We used Fine-Gray sub-distribution hazard models to estimate the association between I-O era and CVD death [sub-distribution hazard ratio (sHR)] as well as NSCLC-specific death (sHR). This approach directly models the CIF and is appropriate when the clinical question concerns absolute risk accumulation over time under competing risks. To aid clinical interpretation, the sHR estimates the relative probability of a specific event occurring, while appropriately accounting for the fact that some patients will experience a competing event (e.g., dying from cancer) and are thus removed from the “at-risk” pool for the primary event
Sensitivity analysis: cause-specific Cox regression
As sensitivity analysis, we fitted cause-specific Cox proportional hazards models for: CVD death (treating non-CVD deaths as competing censoring) and NSCLC death (treating non-NSCLC deaths as competing censoring). Cause-specific hazard ratios (HRs) provide complementary inference regarding instantaneous hazard among those still event-free.
Additional sensitivity/subgroup framework: multistate model
To further evaluate robustness and to support subgroup inference within a coherent event-process framework, we implemented a multistate model representing transitions from “Alive after diagnosis” to each absorbing death state (NSCLC death, CVD death, other death). Era effects were estimated as hazard ratios for transitions under the multistate formulation.
Subgroup analyses
We performed subgroup analyses (pre-specified in your results) across clinically relevant strata, including age groups, sex, race categories, and treatment-related variables, reporting hazard ratios and confidence intervals (CIs). Consistency of effect direction and statistical interaction patterns were assessed descriptively and/or by interaction terms as implemented. All tests were two-sided, with significance set at P<0.05.
Results
Cohort characteristics and matching
The screening process for NSCLC patients is shown in Figure 1. Before matching, the study included 40,585 patients in the Post-I-O and 41,772 patients in the Pre-I-O. As shown in Table 1, baseline characteristics differed between eras across several domains, including sociodemographic variables (e.g., income distribution), tumor features (grade and stage distributions), and treatment patterns (surgery, radiotherapy, chemotherapy), supporting the need for confounding control.
Table 1
| Variable | Post-I-O (n=40,585) | Pre-I-O (n=41,772) | P value |
|---|---|---|---|
| Age, years | 77.19 [5.45] | 77.28 [5.35] | 0.02 |
| Sex | |||
| Male | 20,598 (50.8) | 20,599 (49.3) | <0.001 |
| Female | 19,987 (49.2) | 21,172 (50.7) | |
| Race | |||
| American Indian/Alaska Native | 197 (0.5) | 146 (0.3) | <0.001 |
| Asian or Pacific Islander | 4,625 (11.4) | 3,608 (8.6) | |
| Black | 3,929 (9.7) | 3,619 (8.7) | |
| White | 31,834 (78.4) | 34,398 (82.3) | |
| Marital status | |||
| Married | 20,603 (50.8) | 21,312 (51.0) | 0.47 |
| Others | 19,982 (49.2) | 20,459 (49.0) | |
| Rural urban | |||
| Metropolitan areas | 35,243 (86.8) | 36,056 (86.3) | 0.03 |
| Nonmetropolitan counties | 5,342 (13.2) | 5,715 (13.7) | |
| Median household income | |||
| $50,000–$100,000 | 27,255 (67.2) | 32,204 (77.1) | <0.001 |
| Less than $50,000 | 2,224 (5.5) | 2,900 (6.9) | |
| Over $100,000 | 11,106 (27.4) | 6,667 (16.0) | |
| Tumor grade | |||
| G1 | 1,263 (3.1) | 1,481 (3.5) | <0.001 |
| G2 | 5,942 (14.6) | 6,198 (14.8) | |
| G3 | 8,937 (22.0) | 11,065 (26.5) | |
| G4 | 232 (0.6) | 583 (1.4) | |
| Unknown | 24,211 (59.7) | 22,444 (53.7) | |
| T stage | |||
| T1 | 6,024 (14.8) | 5,111 (12.2) | <0.001 |
| T2 | 9,580 (23.6) | 10,978 (26.3) | |
| T3 | 8,198 (20.2) | 7,139 (17.1) | |
| T4 | 12,552 (30.9) | 14,421 (34.5) | |
| Unknown | 4,231 (10.4) | 4,122 (9.9) | |
| N stage | |||
| N0 | 8,447 (20.8) | 9,094 (21.8) | <0.001 |
| N1 | 3,749 (9.2) | 3,377 (8.1) | |
| N2 | 18,635 (45.9) | 20,923 (50.1) | |
| N3 | 7,521 (18.5) | 6,107 (14.6) | |
| Unknown | 2,233 (5.5) | 2,270 (5.4) | |
| M stage | |||
| M0 | 14,226 (35.1) | 16,306 (39.0) | <0.001 |
| M1 | 26,359 (64.9) | 25,465 (61.0) | |
| Surgery record | |||
| Surgery of primary tumor | 3,359 (8.3) | 3,687 (8.8) | 0.005 |
| None | 37,226 (91.7) | 38,084 (91.2) | |
| Radiotherapy record | |||
| Radiotherapy | 18,766 (46.2) | 18,938 (45.3) | 0.01 |
| No/unknown | 21,819 (53.8) | 22,833 (54.7) | |
| Chemotherapy record | |||
| Chemotherapy | 24,840 (61.2) | 24,308 (58.2) | <0.001 |
| No/unknown | 15,745 (38.8) | 17,463 (41.8) |
Data are presented as n (%) or mean [SD]. G, grade; M, metastasis; N, node; NSCLC, non-small cell lung cancer; Post-I-O, post-immuno-oncology era; Pre-I-O, pre-immuno-oncology era; PSM, propensity score matching; SD, standard deviation; T, tumor.
After 1:1 PSM, 37,169 patients were retained in each era cohort (Table 2). Matching improved balance across key demographic, tumor, and treatment variables. Some variables (notably income distribution and certain staging or race categories) showed residual differences, which is expected in very large samples and was further addressed through multivariable adjustment in regression models.
Table 2
| Variable | Pre-I-O (n=37,169) | Post-I-O (n=37,169) | P value |
|---|---|---|---|
| Age, years | 77.17 [5.32] | 77.21 [5.45] | 0.37 |
| Sex | |||
| Male | 18,595 (50.0) | 18,416 (49.5) | 0.19 |
| Female | 18,574 (50.0) | 18,753 (50.5) | |
| Race | |||
| American Indian/Alaska Native | 145 (0.4) | 165 (0.4) | 0.001 |
| Asian or Pacific Islander | 3,457 (9.3) | 3,684 (9.9) | |
| Black | 3,468 (9.3) | 3,651 (9.8) | |
| White | 30,099 (81.0) | 29,669 (79.8) | |
| Marital status | |||
| Married | 18,838 (50.7) | 18,744 (50.4) | 0.495 |
| Others | 18,331 (49.3) | 18,425 (49.6) | |
| Rural urban | |||
| Metropolitan areas | 32,016 (86.1) | 31,962 (86.0) | 0.58 |
| Nonmetropolitan counties | 5,153 (13.9) | 5,207 (14.0) | |
| Median household income | |||
| $50,000–$100,000 | 28,223 (75.9) | 26,913 (72.4) | <0.001 |
| Less than $50,000 | 2,305 (6.2) | 2,215 (6.0) | |
| Over $100,000 | 6,641 (17.9) | 8,041 (21.6) | |
| Tumor grade | |||
| G1 | 1,201 (3.2) | 1,210 (3.3) | 0.17 |
| G2 | 5,566 (15.0) | 5,557 (15.0) | |
| G3 | 8,856 (23.8) | 8,599 (23.1) | |
| G4 | 252 (0.7) | 232 (0.6) | |
| Unknown | 21,294 (57.3) | 21,571 (58.0) | |
| T stage | |||
| T1 | 4,994 (13.4) | 5,281 (14.2) | <0.001 |
| T2 | 9,332 (25.1) | 9,045 (24.3) | |
| T3 | 6,942 (18.7) | 7,269 (19.6) | |
| T4 | 12,036 (32.4) | 11,778 (31.7) | |
| Unknown | 3,865 (10.4) | 3,796 (10.2) | |
| N stage | |||
| N0 | 7,935 (21.3) | 7,906 (21.3) | <0.001 |
| N1 | 3,222 (8.7) | 3,371 (9.1) | |
| N2 | 17,954 (48.3) | 17,428 (46.9) | |
| N3 | 5,984 (16.1) | 6,480 (17.4) | |
| Unknown | 2,074 (5.6) | 1,984 (5.3) | |
| M stage | |||
| M0 | 13,460 (36.2) | 13,432 (36.1) | 0.84 |
| M1 | 23,709 (63.8) | 23,737 (63.9) | |
| Surgery record | |||
| Surgery of primary tumor | 3,149 (8.5) | 3,155 (8.5) | 0.948 |
| None | 34,020 (91.5) | 34,014 (91.5) | |
| Radiotherapy record | |||
| Radiotherapy | 17,136 (46.1) | 17,223 (46.3) | 0.53 |
| No/unknown | 20,033 (53.9) | 19,946 (53.7) | |
| Chemotherapy record | |||
| Chemotherapy | 22,285 (60.0) | 22,447 (60.4) | 0.23 |
| No/unknown | 14,884 (40.0) | 14,722 (39.6) |
Data are presented as n (%) or mean [SD]. G, grade; M, metastasis; N, node; NSCLC, non-small cell lung cancer; Post-I-O, post-immuno-oncology era; Pre-I-O, pre-immuno-oncology era; PSM, propensity score matching; SD, standard deviation; T, tumor.
Cumulative incidence patterns
Cumulative incidence analyses demonstrated a divergent pattern after the introduction of I-O: the cumulative incidence of NSCLC-specific death was lower in the Post-I-O than in the Pre-I-O, consistent with improved cancer outcomes in routine practice (Figure 2A). In contrast, the cumulative incidence of CVD death was higher in the Post-I-O, indicating an increased absolute burden of cerebrovascular mortality despite improved cancer survival (Figure 2B). These CIF patterns motivated formal competing-risks regression to quantify era-associated differences while appropriately accounting for competing causes of death.
Primary competing-risks analysis
The result of multivariable Fine-Gray competing-risks regression is shown in Table 3. Diagnosis in the Post-I-O was associated with a substantially lower sub-distribution hazard of NSCLC-specific death with a sHR of 0.624 (95% CI: 0.613–0.635; P<0.001). This corresponds to an approximate 38% relative reduction in the sub-distribution hazard of NSCLC death in the Post-I-O, consistent with lower CIF observed in Figure 3. Conversely, diagnosis in the Post-I-O was associated with a higher sub-distribution hazard of CVD death with a sHR of 1.424 (95% CI: 1.293–1.657, P<0.001). This indicates an approximate 42% relative increase in the sub-distribution hazard of CVD death, aligning with the higher CVD CIF in Figure 2.
Table 3
| Variable | NSCLC-specific death | CVD death | |||
|---|---|---|---|---|---|
| sHR (95% CI) | P value | sHR (95% CI) | P value | ||
| Period_post-I-O | 0.624 (0.613–0.635) | <0.001 | 1.424 (1.293–1.657) | <0.001 | |
| Age | 1.008 (1.006–1.01) | <0.001 | 1.01 (1.006–1.014) | <0.001 | |
| Sex_female | 0.897 (0.882–0.913) | <0.001 | 0.872 (0.837–0.909) | <0.001 | |
| Race_Asian or Pacific Islander | 0.857 (0.747–0.982) | 0.03 | 0.997 (0.735–1.352) | 0.98 | |
| Race_Black | 0.965 (0.842–1.106) | 0.61 | 1.142 (0.843–1.546) | 0.39 | |
| Race_White | 1.008 (0.882–1.153) | 0.90 | 1.025 (0.761–1.38) | 0.87 | |
| Marital_status_others | 1.033 (1.014–1.051) | <0.001 | 1.097 (1.052–1.143) | <0.001 | |
| Rural_urban_Nonmetropolitan counties | 1.005 (0.977–1.034) | 0.74 | 1.039 (0.976–1.106) | 0.23 | |
| Median_household_income_less than $50,000 | 1.063 (1.021–1.107) | 0.003 | 1.094 (1.003–1.193) | 0.04 | |
| Median_household_income_over $100,000 | 0.957 (0.937–0.978) | <0.001 | 0.943 (0.896–0.992) | 0.02 | |
| Tumor_grade_G2 | 1.253 (1.192–1.317) | <0.001 | 0.86 (0.771–0.959) | 0.006 | |
| Tumor_grade_G3 | 1.328 (1.265–1.394) | <0.001 | 0.896 (0.806–0.996) | 0.04 | |
| Tumor_grade_G4 | 1.437 (1.284–1.608) | <0.001 | 1.04 (0.81–1.336) | 0.76 | |
| Tumor_grade_unknown | 1.228 (1.172–1.287) | <0.001 | 0.953 (0.861–1.055) | 0.36 | |
| T_stage_T2 | 1.259 (1.224–1.295) | <0.001 | 0.773 (0.728–0.821) | <0.001 | |
| T_stage_T3 | 1.279 (1.24–1.318) | <0.001 | 0.765 (0.717–0.816) | <0.001 | |
| T_stage_T4 | 1.295 (1.26–1.331) | <0.001 | 0.741 (0.699–0.785) | <0.001 | |
| T_stage_unknown | 1.214 (1.169–1.26) | <0.001 | 0.855 (0.789–0.927) | <0.001 | |
| N_stage_N1 | 1.185 (1.144–1.226) | <0.001 | 0.823 (0.761–0.89) | <0.001 | |
| N_stage_N2 | 1.254 (1.225–1.284) | <0.001 | 0.802 (0.761–0.846) | <0.001 | |
| N_stage_N3 | 1.255 (1.22–1.29) | <0.001 | 0.747 (0.698–0.799) | <0.001 | |
| N_stage_unknown | 1.111 (1.064–1.161) | <0.001 | 0.937 (0.85–1.033) | 0.19 | |
| M_stage_M1 | 1.528 (1.498–1.558) | <0.001 | 0.673 (0.644–0.704) | <0.001 | |
| Surgery_record_none | 1.69 (1.634–1.749) | <0.001 | 0.845 (0.79–0.904) | <0.001 | |
| Radiotherapy_record_no/unknown | 0.978 (0.961–0.996) | 0.02 | 0.971 (0.933–1.011) | 0.16 | |
| Chemotherapy_record_no/unknown | 1.218 (1.195–1.241) | <0.001 | 1.277 (1.226–1.33) | <0.001 | |
CI, confidence interval; CVD, cerebrovascular disease; G, grade; M, metastasis; N, node; NSCLC, non-small cell lung cancer; I-O, immuno-oncology; sHR, sub-distribution hazard ratio; T, tumor.
Sensitivity analysis
Cause-specific Cox models yielded consistent directions of association. Diagnosis in the Post-I-O was still significantly associated with a substantially lower risk of NSCLC-specific death (HR =0.739; 95% CI: 0.727–0.752, P<0.001) after adjusting for other confounding factors (Table 4). After adjusting for other confounding factors, patients diagnosed in the Post-I-O correlated with a substantially higher risk of CVD death (HR =1.425; 95% CI: 1.217–1.699; P<0.001). Moreover, female, married cases, metropolitan area, higher median household income, lower tumor grade, without distant disease, surgery of primary tumor and chemotherapy were independent protective factors for both NSCLC-specific death and CVD death. Thus, both the CIF-focused Fine-Gray approach and the instantaneous-risk cause-specific approach supported the conclusion that the Post-I-O is associated with lower NSCLC mortality but higher cerebrovascular mortality.
Table 4
| Variable | NSCLC-specific death | CVD-related death | |||
|---|---|---|---|---|---|
| HR (95% CI) | P value | HR (95% CI) | P value | ||
| Pre-I-O | Ref. | Ref. | |||
| Post-I-O | 0.739 (0.727–0.752) | <0.001 | 1.425 (1.217–1.699) | <0.001 | |
| Age | 1.015 (1.013–1.016) | <0.001 | 1.027 (1.023–1.031) | <0.001 | |
| Male | Ref. | Ref. | |||
| Female | 0.832 (0.817–0.847) | <0.001 | 0.726 (0.697–0.757) | <0.001 | |
| American Indian/Alaska Native | Ref. | Ref. | |||
| Asian or Pacific Islander | 0.888 (0.777–1.015) | 0.08 | 0.98 (0.719–1.336) | 0.90 | |
| Black | 1.087 (0.951–1.243) | 0.22 | 1.346 (0.989–1.833) | 0.058 | |
| White | 1.089 (0.955–1.241) | 0.20 | 1.192 (0.88–1.615) | 0.26 | |
| Married | Ref. | Ref. | |||
| Others | 1.073 (1.053–1.092) | <0.001 | 1.177 (1.129–1.227) | <0.001 | |
| Metropolitan areas | Ref. | Ref. | |||
| Nonmetropolitan counties | 1.042 (1.013–1.072) | 0.004 | 1.105 (1.037–1.177) | 0.002 | |
| $50,000–$100,000 | Ref. | Ref. | |||
| Less$50,000 | 1.057 (1.016–1.1) | 0.006 | 1.103 (1.011–1.205) | 0.03 | |
| Over$100,000 | 0.949 (0.928–0.97) | <0.001 | 0.911 (0.865–0.96) | <0.001 | |
| G1 | Ref. | Ref. | |||
| G2 | 1.392 (1.318–1.471) | <0.001 | 1.192 (1.065–1.334) | 0.002 | |
| G3 | 1.514 (1.435–1.596) | <0.001 | 1.32 (1.184–1.471) | <0.001 | |
| G4 | 1.635 (1.459–1.833) | <0.001 | 1.51 (1.174–1.944) | 0.001 | |
| Unknown | 1.418 (1.347–1.493) | <0.001 | 1.362 (1.227–1.512) | <0.001 | |
| T1 | Ref. | Ref. | |||
| T2 | 1.287 (1.249–1.325) | <0.001 | 0.944 (0.888–1.004) | 0.07 | |
| T3 | 1.353 (1.311–1.396) | <0.001 | 1.004 (0.94–1.073) | 0.90 | |
| T4 | 1.38 (1.341–1.42) | <0.001 | 0.994 (0.937–1.055) | 0.85 | |
| Unknown | 1.313 (1.265–1.363) | <0.001 | 1.081 (0.997–1.171) | 0.059 | |
| N0 | Ref. | Ref. | |||
| N1 | 1.238 (1.196–1.282) | <0.001 | 1.033 (0.954–1.119) | 0.42 | |
| N2 | 1.33 (1.299–1.362) | <0.001 | 1.051 (0.996–1.109) | 0.07 | |
| N3 | 1.328 (1.29–1.366) | <0.001 | 0.973 (0.909–1.042) | 0.44 | |
| Unknown | 1.155 (1.108–1.205) | <0.001 | 1.068 (0.97–1.177) | 0.18 | |
| M0 | Ref. | Ref. | |||
| M1 | 1.7 (1.666–1.735) | <0.001 | 1.054 (1.007–1.103) | 0.02 | |
| Surgery of primary tumor | Ref. | Ref. | |||
| None | 2.027 (1.951–2.106) | <0.001 | 1.562 (1.454–1.677) | <0.001 | |
| Radiotherapy | Ref. | Ref. | |||
| No/unknown | 0.992 (0.974–1.01) | 0.38 | 1.01 (0.968–1.053) | 0.65 | |
| Chemotherapy | Ref. | Ref. | |||
| No/unknown | 1.402 (1.376–1.428) | <0.001 | 1.691 (1.621–1.765) | <0.001 | |
CI, confidence interval; CVD, cerebrovascular disease; G, grade; HR, hazard ratio; M, metastasis; N, node; NSCLC, non-small cell lung cancer; I-O, immuno-oncology; Ref., reference.
Multistate model analysis
In multistate modeling of the death process after diagnosis, Post-I-O diagnosis was associated with a lower overall hazard of all caused death in the modeled transition framework (HR =0.746; 95% CI: 0.735–0.758; P<0.001) (Table 5). This result is consistent with improved overall survival dynamics in the I-O era while still allowing for cause-specific increases in CVD mortality observed in competing-risks models.
Table 5
| Variable | HR (95% CI) | P value |
|---|---|---|
| Period_post-I-O | 0.746 (0.735–0.758) | <0.001 |
| Age | 1.017 (1.015–1.019) | <0.001 |
| Sex_female | 0.813 (0.799–0.826) | <0.001 |
| Race_Asian or Pacific Islander | 0.889 (0.786–1.006) | 0.06 |
| Race_Black | 1.117 (0.988–1.262) | 0.08 |
| Race_White | 1.088 (0.965–1.227) | 0.1701 |
| Marital_status_others | 1.092 (1.074–1.11) | <0.001 |
| Rural_urban_non-metropolitan counties | 1.056 (1.029–1.084) | <0.001 |
| Median_household_income_less than $50,000 | 1.074 (1.036–1.114) | <0.001 |
| Median_household_income_over $100,000 | 0.94 (0.921–0.96) | <0.001 |
| Tumor_grade_G2 | 1.323 (1.26–1.39) | <0.001 |
| Tumor_grade_G3 | 1.443 (1.376–1.514) | <0.001 |
| Tumor_grade_G4 | 1.582 (1.426–1.755) | <0.001 |
| Tumor_grade_unknown | 1.38 (1.318–1.445) | <0.001 |
| T_stage_T2 | 1.179 (1.148–1.21) | <0.001 |
| T_stage_T3 | 1.241 (1.206–1.276) | <0.001 |
| T_stage_T4 | 1.258 (1.226–1.291) | <0.001 |
| T_stage_unknown | 1.23 (1.19–1.273) | <0.001 |
| N_stage_N1 | 1.173 (1.136–1.211) | <0.001 |
| N_stage_N2 | 1.248 (1.221–1.275) | <0.001 |
| N_stage_N3 | 1.226 (1.194–1.259) | <0.001 |
| N_stage_unknown | 1.131 (1.088–1.175) | <0.001 |
| M_stage_M1 | 1.496 (1.469–1.524) | <0.001 |
| Surgery_record_none | 1.84 (1.779–1.902) | <0.001 |
| Radiotherapy_record_no/unknown | 0.991 (0.975–1.008) | 0.28 |
| Chemotherapy_record_no/unknown | 1.451 (1.426–1.476) | <0.001 |
This model accounts for competing risks as transitions from “alive” to all caused death with stratified baseline hazards. CI, confidence interval; G, grade; HR, hazard ratio; I-O, immuno-oncology; M, metastasis; N, node; T, tumor.
Subgroup analyses
Across most subgroups, the era effect on mortality demonstrated consistent directionality, with hazard ratios generally favoring the Post-I-O for the all caused mortality pattern in the subgroup framework (HRs approximately ranging from ~0.71 to ~0.87). Effects were observed across age strata (70–74 years, 75–79 years, 80–84 years, and ≥85 years), across sex, and across most race categories. For smaller race subgroups such as American Indian/Alaska Native, CIs were wider and statistical significance was not consistently achieved, likely reflecting limited sample size and reduced precision. Collectively, these subgroup results suggest that the observed era-associated changes are broadly consistent across clinically relevant patient groups.
Discussion
In this large SEER cohort of adults aged ≥70 years with advanced NSCLC, the introduction of I-O (post-2015) was associated with a dual mortality shift. (I) Lower NSCLC-specific mortality in the Post-I-O, supported by both competing-risks regression (Fine-Gray sHR =0.624) and cause-specific Cox models (HR =0.739), and reflected in lower cumulative incidence curves. Importantly, while standard P values indicate statistical significance, evaluating the minimal clinically important difference (MCID) provides necessary context for clinical relevance. Based on the recently proposed Horita’s MCID framework for evaluating effect sizes other than mean differences, an sHR of 0.624 represents a robust and clinically meaningful reduction in cancer-specific mortality risk, confirming the profound real-world impact of the I-O era (13). (II) Higher CVD mortality in the Post-I-O, again consistent across Fine-Gray (sHR =1.424) and cause-specific Cox (HR =1.425), with higher cumulative incidence of CVD death. These findings indicate that while cancer outcomes improved in the I-O era for older adults with advanced NSCLC, the population-level burden of cerebrovascular mortality increased, highlighting the importance of integrated oncologic and vascular risk management in this vulnerable group.
Several non-mutually exclusive mechanisms may explain the observed rise in CVD death. Improved cancer control reduces early NSCLC deaths, leading to longer survival time and increased exposure window for non-cancer events such as CVD. In competing-risk terms, when the dominant competing event (NSCLC death) becomes less frequent or occurs later, the CIF for other causes (including CVD) can rise because more patients remain alive long enough to experience those events. The concordance between lower NSCLC CIF and higher CVD CIF in Figure 2 supports this explanation.
Immunotherapy can trigger systemic immune activation, endothelial inflammation, and prothrombotic pathways (14-16). Even if ICI-related vascular events are uncommon on an individual basis, small increases in risk may translate to meaningful population-level changes in older adults—particularly those with baseline vascular vulnerability. The observed increase in CVD death in both Fine-Gray and cause-specific Cox frameworks is compatible with either true risk elevation or increased detection/attribution of CVD as cause of death in longer-term survivors.
The post-2015 era includes broader changes: improved imaging, evolving coding, modifications in chemotherapy use, radiotherapy techniques, and supportive care practices. Some of these changes could influence vascular outcomes indirectly. While matching and multivariable adjustment mitigate measured confounding, unmeasured secular factors may contribute to the observed pattern.
Our results emphasize a critical “success paradox” in modern oncology: as cancer mortality declines, non-cancer causes of death become increasingly visible and clinically relevant. For older adults with advanced NSCLC in the I-O era, cerebrovascular mortality appears to represent a growing competing risk. These findings support routine vascular risk assessment at diagnosis and during systemic therapy planning for older NSCLC patients and coordination between oncology and primary care/cardiology/neurology to optimize blood pressure control, diabetes management, lipid therapy, smoking cessation, and antithrombotic strategies when indicated.
This study has several strengths. This study was a large, population-based cohort focused on older adults (≥70 years), who are underrepresented in clinical trials. Moreover, competing-risks methodology (Fine-Gray) aligned with the clinical question of cumulative cause-specific mortality. We also performed robustness checks via cause-specific Cox and multistate models, plus subgroup evaluation of effect consistency.
Several limitations merit consideration: era-based exposure definition does not confirm individual ICI receipt. The Post-I-O reflects availability and adoption but includes treatment heterogeneity. Residual confounding is possible due to unmeasured factors (e.g., smoking history, performance status, frailty, comorbidities not captured or incompletely modeled, detailed regimen lines, and supportive care). Cause-of-death misclassification may occur, particularly between CVD and other causes in older adults. Changes in clinical practice over time (diagnostics, coding, supportive care) may contribute to observed differences independent of immunotherapy. Despite these limitations, the consistent direction of associations across analytic frameworks and alignment between CIF curves and model estimates strengthen confidence in the overall conclusion.
Conclusions
Among Medicare-linked older adults with advanced NSCLC, the Post-I-O [2015–2022] was associated with substantially lower NSCLC-specific mortality but higher CVD mortality compared with 2007–2014. These results underscore the need for integrated cancer and vascular risk management to maximize long-term benefits of I-O in older patients. Specifically, these findings highlight the increasingly vital role of the cardio-oncology discipline in providing comprehensive, interdisciplinary survivorship care for the aging cancer population.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://tcr.amegroups.com/article/view/10.21037/tcr-2026-1-0169/rc
Peer Review File: Available at https://tcr.amegroups.com/article/view/10.21037/tcr-2026-1-0169/prf
Funding: This study was funded by
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tcr.amegroups.com/article/view/10.21037/tcr-2026-1-0169/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Liu L, Soler J, Reckamp KL, et al. Emerging Targets in Non-Small Cell Lung Cancer. Int J Mol Sci 2024;25:10046. [Crossref] [PubMed]
- Xue M, Ma L, Zhang P, et al. New insights into non-small cell lung cancer bone metastasis: mechanisms and therapies. Int J Biol Sci 2024;20:5747-63. [Crossref] [PubMed]
- Desai A, Peters S. Immunotherapy-based combinations in metastatic NSCLC. Cancer Treat Rev 2023;116:102545. [Crossref] [PubMed]
- Hendriks LEL, Remon J, Faivre-Finn C, et al. Non-small-cell lung cancer. Nat Rev Dis Primers 2024;10:71. [Crossref] [PubMed]
- Tsukita Y, Tozuka T, Kushiro K, et al. Immunotherapy or Chemoimmunotherapy in Older Adults With Advanced Non-Small Cell Lung Cancer. JAMA Oncol 2024;10:439-47. [Crossref] [PubMed]
- Adams HP Jr. Cancer and Cerebrovascular Disease. Curr Neurol Neurosci Rep 2019;19:73. [Crossref] [PubMed]
- Goldstein LB. Introduction for Focused Updates in Cerebrovascular Disease. Stroke 2020;51:708-10. [Crossref] [PubMed]
- Kulkarni P, Do D, Shrestha S, et al. Cancer Immunotherapy-An Overview. Cancer Treat Res 2025;129:1-16. [Crossref] [PubMed]
- Michot JM, Bigenwald C, Champiat S, et al. Immune-related adverse events with immune checkpoint blockade: a comprehensive review. Eur J Cancer 2016;54:139-48. [Crossref] [PubMed]
- Taylor X, Noristani HN, Fitzgerald GJ, et al. Amyloid-β (Aβ) immunotherapy induced microhemorrhages are linked to vascular inflammation and cerebrovascular damage in a mouse model of Alzheimer's disease. Mol Neurodegener 2024;19:77. [Crossref] [PubMed]
- Xiong M, Jiang H, Serrano JR, et al. APOE immunotherapy reduces cerebral amyloid angiopathy and amyloid plaques while improving cerebrovascular function. Sci Transl Med 2021;13:eabd7522. [Crossref] [PubMed]
- Che WQ, Li YJ, Tsang CK, et al. How to use the Surveillance, Epidemiology, and End Results (SEER) data: research design and methodology. Mil Med Res 2023;10:50. [Crossref] [PubMed]
- Horita N, Yamamoto S, Mizuki Y, et al. Minimal Clinically Important Difference (MCID) of Effect Sizes other than Mean Difference. Journal of Clinical Question 2024;1:116-27.
- Martins F, Sofiya L, Sykiotis GP, et al. Adverse effects of immune-checkpoint inhibitors: epidemiology, management and surveillance. Nat Rev Clin Oncol 2019;16:563-80. [Crossref] [PubMed]
- Mahalingam P, Newsom-Davis T. Cancer immunotherapy and the management of side effects. Clin Med (Lond) 2023;23:56-60. [Crossref] [PubMed]
- Brown TJ, Mamtani R, Bange EM. Immunotherapy Adverse Effects. JAMA Oncol 2021;7:1908. [Crossref] [PubMed]

