Analyses of lncRNA-GC1 expression profiles in patients with GC undergoing ICI treatments
Initially, differential centrifugation was used to purify EVs from patient plasma samples. TEM and NTA analyses then confirmed that the particles were 30–200 nm in size (eFigures 1A–B). Western immunoblotting also confirmed the particles to be positive for the EV surface biomarkers CD9, CD63, and CD81, confirming the high quality of the isolated EVs (eFigures 1C). To exclude potential interference caused by variable EV content concentrations on our findings, absolute lncRNA-GC1 expression was normalized to λ polyA+ RNA for establishing relative expression levels. For patients in cohort 1, no significant changes in the expression of the EV-related biomarkers CD9/CD63/CD81 were observed when comparing baseline samples with those obtained during treatment (eFigure 2), suggesting that EV content remained stable throughout ICI treatments. These findings thus confirmed that EV-derived lncRNA-GC1 can be reliably quantified in patients undergoing ICI treatments.
The concordance between lncRNA-GC1 expression in EV and tissue samples was then examined. EV-derived lncRNA-GC1 and PD-L1 both showed no significant differences between CPS negative (CPSneg, combined positive score, CPS = 0) and CPS positive (CPSpos, CPS ≥ 1) tissue samples (eFigures 3A–B). Moreover, no significant correlations were observed between EV-derived expression of the T cell marker CD3 or the B cell marker CD19 and circulating T or B cell populations (eFigures 3C–D). These results suggest that circulating EVs exhibit a transcriptional spectrum that is distinct from those of tissues or circulating cells in patients with GC.
Baseline lncRNA-GC1 levels predict the survival outcomes of patients with GC receiving ICI treatments independent of PD-L1 expression status and CD8 + T cell infiltration
In a previous study, we identified a 5.6-fold EV-derived lncRNA-GC1 cutoff value as a reliable threshold for the detection of GC in a relatively large patient cohort.29 This cutoff value was also used for all subsequent analyses in the current study, as it was derived from a GC patient population of almost 1,000—which was considered to be a sufficiently large population to allow for the minimization of batch effects between cohorts, thus ensuring that we could perform simple and reliable predictive analyses.
Based on this cutoff value, the patients in cohort 1 were assigned to either the lncRNA-GC1high group (> 5.6-fold, n = 41, 30.1%) or lncRNA-GC1low group (≤ 5.6-fold, n = 95, 69.9%). Correlations between EV-derived lncRNA-GC1 and ICI treatment responses were then assessed. Using the RECIST criteria, patients were classified as either responders (CR + PR) or non-responders (SD + PD). EV-derived lncRNA-GC1 levels were found to be significantly lower in the responders vs the non-responders (P = 0.0002, Fig. 2A). In addition, the proportion of responders in the lncRNA-GC1high group was significantly lower than that among the lncRNA-GC1low patients (16.0% vs. 52.5%, P < 0.0001, Fig. 2B). The patients who had high baseline levels of EV-derived lncRNA-GC1 had poor PFS (hazard ratio [HR], 2.628; 95% confidence interval [CI], 1.523–4.534; P < 0.0001; Fig. 2C) and OS (HR, 3.432; 95% CI; 1.894–6.217, P < 0.0001; Fig. 2D).
We then proceeded to evaluate the correlations between EV-derived lncRNA-GC1 and patient survival outcomes as a function of the PD-L1 expression status and CD8+ T cell infiltration density. In total, 57.4% of the samples were negative for PD-L1 (CPS = 0; 78/136; PD-L1neg), while 42.6% were PD-L1-positive (CPS ≥ 1; 58/136; PD-L1pos). The levels of EV-derived lncRNA-GC1 showed no significant differences between PD-L1neg and PD-L1pos GC (P = 0.6347, Fig. 3A). The correlations between EV-derived lncRNA-GC1 and survival outcomes as a function of CD8+ T cell infiltration (CD8high vs. CD8low) were then assessed. In total, 36.0% had high CD8+ T cell infiltration density (49/136; CD8high), while 64.0% had low CD8+ T cell infiltration density (89/136; CD8low). No significant differences in EV-derived lncRNA-GC1 expression were observed when comparing patients with low vs high levels of CD8+ T cell infiltration (P = 0.6311; Fig. 3B).
Among the patients who were PD-L1pos, the median PFSs in the lncRNA-GC1high and lncRNA-GC1low groups were 14.30 and 24.30 months, respectively (P = 0.0007), while the median OSs were 17.60 and 27.30 months, respectively (P = 0.0102, Fig. 3C). Among the PD-L1neg patients, the median PFSs in the lncRNA-GC1low and lncRNA-GC1high groups were 11.40 and 18.70 months, respectively (P = 0.0076, Fig. 3D), while the median OSs were 16.70 and 29.80 months, respectively (P = 0.0002, Fig. 3D). These results suggest that the patients with low EV-derived lncRNA-GC1 levels exhibited better survival outcomes than those with high levels, independent of PD-L1 expression status. The patients with low EV-derived lncRNA-GC1 levels were associated with superior PFS and OS outcomes than those with high levels, independently of the density of CD8+ T cell infiltration (all P < 0.05, Figs. 3E–F). These data suggest that EV-derived lncRNA-GC1 can be used to predict the survival outcomes of patients with GC receiving ICI treatments, independently of PD-L1 expression status or CD8+ T cell infiltration density.
Dynamic changes of EV-derived lncRNA-GC1 enable the monitoring of ICI treatment responses
Based on the observed predictive value of baseline EV-derived lncRNA levels, we then assessed the monitoring role of this biomarker throughout ICI treatments. By comparing EV-derived lncRNA-GC1 levels at baseline and one month, the patients in cohort 1 were classified into four groups: group low (from baseline-high to one month-low), group high (from baseline-low to one month-high), group remain-low (from baseline-low to one month-low), and group remain-high (from baseline-high to one month-high). During the first month of ICI treatments, EV-derived lncRNA-GC1 was able to effectively predict the survival outcomes of these patients (PFS: HR, 3.104; 95% CI, 1.612–5.976; P < 0.0001; OS: HR, 2.818; 95% CI, 1.473–5.391; P < 0.0001; Fig. 4A).
Considering all available time points, EV-derived lncRNA-GC1 levels were found to be lower in responders (CR + PR) vs non-responders (PD + SD; P < 0.05, Fig. 4B), supporting the prognostic relevance of this biomarker. When assessing EV-derived lncRNA-GC1 levels at baseline, between 1–3, 4–6, and ≥ 6 months following the start of treatment, these levels were lower in the patients who achieved overall response (ORR), disease control (DCR), non-progression (PFS), or remained alive (OS) relative to those that did not meet these endpoints (Figs. 4C–F). Notably, timeline-scaled EV-derived lncRNA-GC1 levels remained high in the high and remain-high groups while remaining low in the low and remain-low groups, which highlighted the ability of EV-derived lncRNA-GC1 to represent a stable indicator that can be used to stratify patients throughout treatment (Fig. 4G). These results confirm the value of EV-derived lncRNA-GC1 as a powerful biomarker capable of dynamically monitoring the response to treatment of patients with GC who are treated via ICI therapy.
Retrospective and prospective validation confirms the value of lncRNA-GC1 as a biomarker for predicting and monitoring ICI treatment response
To confirm the reliability of the predictive value of EV-derived lncRNA-GC1 in patients with GC who undergo ICI treatment, two retrospective cohorts and one prospective cohort were formed, which were then used for external validation using the same criteria and cutoff threshold used for cohort 1. For both retrospective validation cohorts (cohorts 2 and 3), baseline EV-derived lncRNA-GC1 levels were negatively correlated with patient PFS (eFigures 4A and 5A) and OS (eFigures 4B and 5B). In addition, patients who achieved ORR (eFigures 4C and 5C), DCR (eFigures 4D and 5D), non-progression (eFigures 4E and 5E), or remained alive (eFigures 4F and 5F) exhibited lower levels of EV-derived lncRNA-GC1 relative to those that did not. Notably, timeline-scaled EV-derived lncRNA-GC1 levels remained high for patients in the high and remain-high groups, and low for those in the low and remain-low groups (eFigures 4G and 5G). Furthermore, the proportion of patients with DCB in the lncRNA-GC1low group was significantly higher than that in the lncRNA-GC1high group (eFigures 4H and 5H). The area under the Kaplan-Meier curve (AUC) values for EV-derived lncRNA-GC1 in these two retrospective cohorts were 0.8821 and 0.8748 (eFigures 4I and 5I).
In our prospective validation analyses (cohort 5), compared to the lncRNA-GC1high group, PFS and OS improved in patients with low levels of EV-derived lncRNA-GC1 (Figs. 5A–B). In addition, patients who achieved ORR (Fig. 5C), DCR (Figs. 5D and 5H), non-progression (Fig. 5E), or remained alive (Fig. 5F) exhibited lower levels of EV-derived lncRNA-GC1 for treatment. The timeline-scaled EV-derived lncRNA-GC1 remained high in patients in the high and remain-high groups, and low in those in the low and remain-low groups (Fig. 5G). Critically, EV-derived lncRNA-GC1 remained independently associated with patient PFS and OS (eTable 2), with an AUC of 0.8778, confirming its predictive value (Fig. 5I).
EV-derived lncRNA-GC1 is associated with unique tumor microenvironmental features
Given that EV-derived lncRNA-GC1 was found to be related to the response of our patients with GC to ICI treatments, correlations between EV-derived lncRNA-GC1 levels and the intratumoral immune landscape were further evaluated. Tertiary lymphoid structures (TLS), T follicular helper cells, and B cells were found to be significantly correlated with EV-derived lncRNA-GC1 (eFigures 6A–B). Moreover, low levels of EV-derived lncRNA-GC1 were significantly correlated with high levels of activated CD8+ T/NK cells and an increased TH1/TH2 ratio, thus indicating robust antitumor immunity (eFigures 6A–B). In addition, significantly more TLS-positive patients were observed in the lncRNA-GC1low group vs the lncRNA-GC1high one (50% vs. 24%, P = 0.003, eFigure 6C). The levels of EV-derived lncRNA-GC1 may correspond to the unique features of the tumor microenvironment that are associated with ICI treatment outcomes.
EV-derived lncRNA-GC1 remains stable under routine clinical testing conditions
To further establish the potential real-world clinical utility of monitoring dynamic changes of EV-derived lncRNA-GC1 as a means of guiding clinical ICI treatment planning, the stability of this biomarker was then assessed—under the presumption that poor stability would severely restrict any potential for the clinical implementation of these results. A retrospective expanding cohort (cohort 4) comprised of 30 patients was enrolled for the purpose. No significant changes in EV-derived lncRNA-GC1 levels were observed when these samples were treated with RNase (P = 0.524, eFigure 7A), incubated at room temperature for prolonged periods (P = 0.353, eFigure 7B), or subjected to repeated freeze-thaw cycles (P = 0.413, eFigure 7C). Notably, circulating lncRNA-GC1 levels were primarily derived from EVs, as was confirmed by the absence of EV-depleted sera (P < 0.0001, eFigure 7D). A significant correlation between EV-derived lncRNA-GC1 and total circulating lncRNA-GC1 levels was also observed (r = 0.7615, P = 0.0011, eFigure 7E). Together, these data indicate that EV encapsulation can shield lncRNA-GC1 and protect it from degradation, providing the levels of stability necessary for effective clinical use.