An 8-Gene Bevacizumab Resistance Signature Predicts Prognosis and Reveals Immunosuppressive Microenvironment in Colorectal Cancer
Discuss this preprint
Start a discussion What are Sciety discussions?Listed in
This article is not in any list yet, why not save it to one of your lists.Abstract
Background
Bevacizumab resistance severely limits long-term efficacy in metastatic colorectal cancer (CRC). This study aimed to develop and validate a bevacizumab resistance-associated gene signature for prognosis prediction and immune microenvironment characterization in CRC.
Methods
Two GEO datasets (GSE19862, GSE86582) with bevacizumab response data and TCGA-COAD/READ RNA-seq data were analyzed. Overlapping differentially expressed genes (DEGs) linked to both CRC progression and bevacizumab resistance were identified. An 8-gene signature (AXIN2, PSORS1C1, KRT74, SLC2A3, STIL, IL33, GALNT6, HSD11B2) was constructed via univariate Cox and LASSO-Cox regression.
Results
In the TCGA cohort, high-risk patients had shorter overall survival (OS; log-rank P < 0.0001). Time-dependent ROC yielded 1-year AUC = 0.638, 3-year AUC = 0.657, and 5-year AUC = 0.757. Multivariate Cox regression confirmed the risk score as an independent prognostic factor. External validation in GSE39582 (optimal cutoff = −1.49) replicated these findings: high-risk patients had inferior OS ( P = 0.0016) with acceptable 1/3/5-year AUCs and retained independent prognostic value (HR = 1.634, P = 0.00415). CIBERSORT and ESTIMATE analyses showed that the high-risk group was characterized by increased M2 macrophages and neutrophils, higher immune and stromal scores, and reduced activated memory CD4 + T cells, monocytes, and activated dendritic cells (all P < 0.05). GSEA highlighted enrichment of TNF-α/NF-κB, IL-6/JAK/STAT3, and immune checkpoint pathways in the high-risk group. AXIN2 (HR = 0.829, P = 0.032) was an independent protective factor, while PSORS1C1 (HR = 1.356, P = 0.048) was an independent risk factor.
Conclusion
The 8-gene bevacizumab resistance signature robustly predicts prognosis and reflects an immunosuppressive microenvironment closely linked to bevacizumab failure in CRC. These findings provide novel insights into immune-mediated resistance and support clinical risk stratification.