Evolution of tumor microenvironment in colorectal liver metastases under treatment stress
Authors: Na Huang, Dongqiang Zeng, Xiaoxiang Rong, Chunlin Wang, Zhenzhen Wu, Jian Guo, Yuqi Wang, Jin Li, Jing Li, Jiao Wang, Siting Zheng, Genjie Huang, Jianping Bin, Yulin Liao, Qian Li, Xin Yi, Wangjun Liao, Min Shi
Journal: Cancer Communications (2022)

The tumor microenvironment (TME) heavily impacts disease biology and may influence responses to systemic treatments, and thereby, affects patients' prognosis. In our previous study, we found that immune features could predict prognosis and guide the therapy choices for stage I-III colon cancer [1, 2]. Increasing evidence shows that therapy-induced TME changes can promote tumor progression, metastasis, and the development of resistance [3, 4]. However, the TME dynamics in colorectal liver metastases (CRLM) under treatment are still incompletely clear. In this study, we surveyed the evolution of the TME under treatment stress from 15 CRLM patients via RNA sequencing (RNA-seq), T-cell receptor (TCR) sequencing (TCR-seq), immunohistochemical (IHC) staining and immunofluorescence (Supplementary Materials and Methods) (Supplementary Figure S1A, and Supplementary Table S1). During treatment, the fraction of epithelial tumor cells of regressive tumors (partial response [PR]) decreased significantly, whereas it increased in stable and progressive tumors (SD and PD) (Supplementary Figure S1B-D). Consensus molecular subtypes, which is one of the most robust colorectal cancer classifications, changed dramatically during treatment (Supplementary Figure S1E). Comparison of treatment-naïve patients with response (PR) and non-response (SD or PD) showed that along with the upregulation of interferon response, the expression of immune checkpoint-related genes also increased (Supplementary Figure S2A-C). Immune-related genes were significantly enriched in sensitive tumors after treatment whereas progressive tumors contained highly expressed genes involved in fatty acid, amino acid metabolism and cholesterol homeostasis (Supplementary Figure S2D-F). Analysis of matched pre- and post-therapy samples revealed that upregulated genes in responders during therapy were immune-related, whereas these genes were downregulated in non-responders (Supplementary Figure S2G-I). Gene Set Enrichment Analysis (GSEA) revealed that the response of interferon and inflammation was upregulated during treatment in responders, but was downregulated in non-responders (Supplementary Figure S2J). Single sample gene set enrichment analysis (ssGSEA) [5] indicated that immune infiltration in regressive tumors was increased, whereas tumor cell cycle and DNA metabolic processes were decreased during therapy (Figure 1A, and Supplementary Figure S3A-B). Conversely, immune-activation relevant signature scores decreased, whereas tumor cell-cycle processes fluctuated in non-responders (Supplementary Figure S3A, C). As shown in Figure 1B, with the application of gene sets [6], immune infiltration in responsive tumors after treatment was higher than before treatment. Dendritic cells, antigen-presenting cells and T cells were upregulated in shrinking tumors but not changed significantly in progressive metastatic tumors. The good-prognosis angiogenesis genes (GPAGs) score [7], which may reflect vessel normalization, was increased after treatment in responders, whereas it was decreased in non-responders. Overall, the RNA-based immune signatures from tumor tissue suggested that therapy remodeled TME, which correlated with treatment benefits. IHC analysis showed that the densities of CD8+ T cells in tumors were significantly increased after treatment in responders (P = 0.011; Figure 1C). However, in non-responders, CD8+ T-cell densities were slightly but not significantly decreased. TCR-seq indicated that four TCR diversity indexes of tumors were significantly increased in responders during treatment but were unchanged in non-responders (Supplementary Figure S4A-C). TCR diversity scores which were defined according to the sum of four diversity indexes, increased as the tumors shrank, and decreased as the tumors progressed (Figure 1D-E, and Supplementary Figure S4D). These findings confirmed responsive tumors had immune remodeling, including elevated T-cell infiltration and activation, increased TCR diversity, as well as enhanced antigen presentation. Deconvoluted by the CIBERSORT algorithm [8], the fraction of neutrophils and monocytes were significantly more numerous in progressive tumors (Figure 1F, and Supplementary Figure S5A). IHC staining confirmed that neutrophil density was significantly increased after treatment in non-responders (P = 0.032; Figure 1G), whereas there were no changes in tumor-associated macrophage (Supplementary Figure S5B). Meanwhile, we deconvoluted bulk RNA-seq profiles of the TCGA Pan-Cancer dataset and validated the poor outcome of patients with high neutrophil infiltration (Supplementary Figure S5C). Remarkably, correlation analysis of the TCGA Pan-Cancer data revealed a significant positive correlation between neutrophils and fatty acid metabolism, and a negative correlation between neutrophils and T-cell immunity (Supplementary Figure S5D-E). These results suggested a major role of fatty acid metabolism–neutrophil crosstalk in treatment resistance. Transcriptome analyses showed that VEGF-A inhibition increased intra-tumor T-cell abundance and upregulated the expression of immune checkpoint-relevant genes after treatment (Supplementary Figure S6A). GSEA and ssGSEA analysis indicated that tumors treated with bevacizumab plus chemotherapy significantly upregulated the interferon response, TNF-α signaling pathways and T-cell receptor signaling pathway (Supplementary Figure S6B-C), suggesting that chemotherapy with bevacizumab promoted tumor T-cell infiltration and activation. Analysis of responders who had received chemotherapy plus bevacizumab revealed that VEGF-A inhibition activated T cells, facilitated leukocyte migration, and upregulated inflammatory responses (Supplementary Figure S6D-E). Responders had decreased vessel density (P = 0.043) and increased proportion of normalized vessels (α-SMA+ vessels; P = 0.023) after treatment (Figure 1H). The HIF1α expression in responders after treatment was lower than before treatment (P = 0.023) and in non-responders after treatment (P = 0.016; Supplementary Figure S6F). These results suggested that the infiltration of T-cell was associated with tumor vasculature normalization. In conclusion, under treatment stress, the TME was remodeled to promote T-cell immunity and vasculature normalization in responders, whereas in non-responders, both lipid metabolism in the tumor and neutrophil infiltration in the TME increased (Figure 1I). Based on the dynamic analysis of clinical samples, our findings are closer to the clinical settings, contributing to the improvement of therapy, and providing valuable clues for further exploration into the relationship between resistance, lipid metabolism and the TME. Not applicable. This work was supported by grants from the National Natural Science Foundation of China (81772580), Natural Science Foundation of Guangdong Province, China (2021A1515011705), and the Clinical Research Startup Program of Southern Medical University by High-level University Construction Funding of Guangdong Provincial Department of Education (LC2016ZD014). The authors declare that they have no competing interests. Patients' samples were collected and analyzed after informed consents were obtained and approved by the ethics committee (NFEC-2017-206) of Nanfang Hospital, Southern Medical University (Guangzhou, China). Study design: MS and WJL; Clinical data acquisition and sample annotation: NH, XXR, CLW, ZZW, JG and JW; Bioinformatics analysis and data visualization: DQZ, XY, YQW, Jin Li, Jing Li, and QL; Performance and analysis of immunohistochemistry and immunofluorescence: GJH and STZ; Statistical analysis: NH, DQZ, and XY; Drafting and critical revision of the manuscript: MS, WJL, NH, DQZ, JPB, and YLL. All authors have read and approved the final manuscript. The datasets generated and analyzed during the current study are available in the Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE136114; accession number: GES136114). Not applicable. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.