Glycolysis and interferon transcription in melanoma immunotherapy

Abstract

Glycolysis and interferon signalling have both been associated with responses to PD-1 blockade, but their joint transcriptional relationship with treatment outcome remains uncertain. We examined whether a state with high glycolysis and low interferon- response was associated with reduced disease control in melanoma. Fixed Hallmark gene sets were scored by single-sample gene set enrichment analysis in GSE91061, GSE78220 and an additional Liu cohort. Disease control analyses included 47, 26 and 100 patients. Conditional odds ratios for the joint state were 6.06 (exact 95% confidence interval, 1.04 to 65.61), 0.35 (0.026 to 3.11) and 1.03 (0.39 to 2.73), respectively. Holm-adjusted P values for the two initial comparisons were 0.073 and 0.365. The additional exploratory P value was 1.000. Continuous models and analyses that omitted each patient in turn did not provide consistent evidence of an adverse association. IFN scores were negatively correlated with DNA-derived tumour purity in 102 patients (Spearman {rho} = -0.671). Among 40 patients with paired biopsies, disease control was associated with higher on-treatment IFN scores after baseline adjustment, although the paired subset had lower baseline glycolysis than patients without paired profiles. In patient-level single-cell analyses, the baseline CD8 T cell correlation weakened when an alternative rank score used a common gene background. The defined joint state was not consistently associated with poorer disease control across cohorts. Bulk IFN scores were related to tumour purity, and cellular correlations were sensitive to scoring choices. These transcriptomic results do not measure pyruvate concentration, glycolytic flux or STAT1 pyruvylation.