Choline metabolism reprogramming mediates an immunosuppressive microenvironment in non-small cell lung cancer (NSCLC) by promoting tumor-associated macrophage functional polarization and endothelial cell proliferation.

Lung cancer is a prevalent malignancy globally, and immunotherapy has revolutionized its treatment. However, resistance to immunotherapy remains a challenge. Abnormal cholinesterase (ChE) activity and choline metabolism are associated with tumor oncogenesis, progression, and poor prognosis in multiple cancers. Yet, the precise mechanism underlying the relationship between ChE, choline metabolism and tumor immune microenvironment in lung cancer, and the response and resistance of immunotherapy still unclear. Firstly, 277 advanced non-small cell lung cancer (NSCLC) patients receiving first-line immunotherapy in Sun Yat-sen University Cancer Center were enrolled in the study. Pretreatment and the alteration of ChE after 2 courses of immunotherapy and survival outcomes were collected. Kaplan-Meier survival and cox regression analysis were performed, and nomogram was conducted to identify the prognostic and predicted values. Secondly, choline metabolism-related genes were screened using Cox regression, and a prognostic model was constructed. Functional enrichment analysis and immune microenvironment analysis were also conducted. Lastly, to gain further insights into potential mechanisms, single-cell analysis was performed. Firstly, baseline high level ChE and the elevation of ChE after immunotherapy were significantly associated with better survival outcomes for advanced NSCLC. Constructed nomogram based on the significant variables from the multivariate Cox analysis performed well in discrimination and calibration. Secondly, 4 choline metabolism-related genes (MTHFD1, PDGFB, PIK3R3, CHKB) were screened and developed a risk signature that was found to be related to a poorer prognosis. Further analysis revealed that the choline metabolism-related genes signature was associated with immunosuppressive tumor microenvironment, immune escape and metabolic reprogramming. scRNA-seq showed that MTHFD1 was specifically distributed in tumor-associated macrophages (TAMs), mediating the differentiation and immunosuppressive functions of macrophages, which may potentially impact endothelial cell proliferation and tumor angiogenesis. Our study highlights the discovery of ChE as a prognostic marker in advanced NSCLC, suggesting its potential for identifying patients who may benefit from immunotherapy. Additionally, we developed a prognostic signature based on choline metabolism-related genes, revealing the correlation with the immunosuppressive microenvironment and uncovering the role of MTHFD1 in macrophage differentiation and endothelial cell proliferation, providing insights into the intricate workings of choline metabolism in NSCLC pathogenesis.

Xiao B ，Li G ，Gulizeba H ，Liu H ，Sima X ，Zhou T ，Huang Y ... -  《Journal of Translational Medicine》

Identification and Validation of a Tumor Microenvironment-Related Gene Signature for Prognostic Prediction in Advanced-Stage Non-Small-Cell Lung Cancer.

The development of immunotherapy has greatly changed the advanced-stage non-small-cell lung cancer (NSCLC) treatment landscape. The complexity and heterogeneity of tumor microenvironment (TME) lead to discrepant immunotherapy effects among patients at the same pathologic stages. This study is aimed at exploring potential biomarkers of immunotherapy and accurately predicting the prognosis for advanced NSCLC patients. RNA-seq data and clinical information on stage III/IV NSCLC were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). In TCGA-NSCLC with stage III/IV (n = 192), immune scores and stromal scores were calculated by using the ESTIMATE algorithms. Univariate, LASSO, and multivariate Cox regression analyses were performed to screen prognostic TME-related genes (TMERGs) and constructed a gene signature risk score model. It was validated in external dataset including GSE41271 (n = 91) and GSE81089 (n = 36). Additionally, a nomogram incorporating TMERG signature risk score and clinical characteristics was established. Further, we accessed the proportion of 22 types of tumor-infiltrating immune cells (TIIC) from the CIBERSORT website and analyzed the difference between two risk groups. OS of patients with high immune/stromal scores were higher (log-rank P = 0.044/log-rank P = 0.048). Multivariate Cox regression identified six prognostic TMERGs, including CD200, CHI3L2, CNTN1, CTSL, FYB1, and SLC52A1. We developed a six-gene risk score model, which was validated as an independent prognostic factor for OS (HR: 3.32, 95% CI: 2.16-5.09). Time-ROC curves showed useful discrimination for TCGA-NSCLC cohort (1-, 2-, and 3-year AUCs were 0.718, 0.761, and 0.750). The predictive robustness was validated in the external dataset. The C-index and 1-, 2-, and 3-year AUCs of nomogram were the largest, which demonstrated the nomogram had the greatest predictive accuracy and effectiveness and could be used for clinical guidance. Besides, the increased infiltration of T cells regulatory (Tregs) and macrophages M2 in the high-risk group suggested that chronic inflammation may reduce survival probability in patients with advanced NSCLC. We conducted a comprehensive analysis of the tumor microenvironment and identified the TMERG signature, which could predict prognosis accurately and provide a reference for the personalized immunotherapy for advanced NSCLC patients.

Zhang X ，Shi X ，Zhao H ，Jia X ，Yang Y ... -  《-》

Construction of a novel choline metabolism-related signature to predict prognosis, immune landscape, and chemotherapy response in colon adenocarcinoma.

Colon adenocarcinoma (COAD) is a common digestive system malignancy with high mortality and poor prognosis. Accumulating evidence indicates that choline metabolism is closely related to tumorigenesis and development. However, the efficacy of choline metabolism-related signature in predicting patient prognosis, immune microenvironment and chemotherapy response has not been fully clarified. Choline metabolism-related differentially expressed genes (DEGs) between normal and COAD tissues were screened using datasets from The Cancer Genome Atlas (TCGA), Kyoto Encyclopedia of Genes and Genomes (KEGG), AmiGO2 and Reactome Pathway databases. Two choline metabolism-related genes (CHKB and PEMT) were identified by univariate and multivariate Cox regression analyses. TCGA-COAD was the training cohort, and GSE17536 was the validation cohort. Patients in the high- and low-risk groups were distinguished according to the optimal cutoff value of the risk score. A nomogram was used to assess the prognostic accuracy of the choline metabolism-related signature. Calibration curves, decision curve analysis (DCA), and clinical impact curve (CIC) were used to improve the clinical applicability of the prognostic signature. Gene Ontology (GO) and KEGG pathway enrichment analyses of DEGs in the high- and low-risk groups were performed. KEGG cluster analysis was conducted by the KOBAS-i database. The distribution and expression of CHKB and PEMT in various types of immune cells were analyzed based on single-cell RNA sequencing (scRNA-seq). The CIBERSORT and ESTIMATE algorithms evaluated tumor immune cell infiltration in the high- and low-risk groups. Evaluation of the half maximal inhibitory concentration (IC50) of common chemotherapeutic drugs based on the choline metabolism-related signature was performed. Small molecule compounds were predicted using the Connectivity Map (CMap) database. Molecular docking is used to simulate the binding conformation of small molecule compounds and key targets. By immunohistochemistry (IHC), Western blot, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) experiments, the expression levels of CHKB and PEMT in human, mouse, and cell lines were detected. We constructed and validated a choline metabolism-related signature containing two genes (CHKB and PEMT). The overall survival (OS) of patients in the high-risk group was significantly worse than that of patients in the low-risk group. The nomogram could effectively and accurately predict the OS of COAD patients at 1, 3, and 5 years. The DCA curve and CIC demonstrate the clinical utility of the nomogram. scRNA-seq showed that CHKB was mainly distributed in endothelial cells, while PEMT was mainly distributed in CD4+ T cells and CD8+ T cells. In addition, multiple types of immune cells expressing CHKB and PEMT differed significantly. There were significant differences in the immune microenvironment, immune checkpoint expression and chemotherapy response between the two risk groups. In addition, we screened five potential small molecule drugs that targeted treatment for COAD. Finally, the results of IHC, Western blot, and qRT-PCR consistently showed that the expression of CHKB in human, mouse, and cell lines was elevated in normal samples, while PMET showed the opposite trend. In conclusion, we constructed a choline metabolism-related signature in COAD and revealed its potential application value in predicting the prognosis, immune microenvironment, and chemotherapy response of patients, which may lay an important theoretical basis for future personalized precision therapy.

Liu C ，Liu D ，Wang F ，Liu Y ，Xie J ，Xie J ，Xie Y ... -  《Frontiers in Immunology》

Single-Cell RNA Sequencing Reveals the Heterogeneity of Tumor-Associated Macrophage in Non-Small Cell Lung Cancer and Differences Between Sexes.

Non-Small Cell Lung Cancer (NSCLC) is a disease with high morbidity and mortality, which has sex-related differences in prognosis and immunotherapy efficacy. However, the difference in the mechanisms remains unclear. Macrophages, characterized by high plasticity and heterogeneity, act as one of the key cells that exert anti-tumor effects in the tumor microenvironment (TME) and play a complicated role in the process of tumor progression. To elucidate the subtype composition and functional heterogeneity of tumor-associated macrophages (TAMs) in NSCLC and further compare the sex-mediated differences, we conducted a single-cell level analysis in early-stage smoking NSCLC patients, combined with ssGSEA analysis, pseudotime ordering, and SCENIC analysis. We found two universally presented immune-suppressive TAMs with different functional and metabolic characteristics in the TME of NSCLC. Specifically, CCL18+ macrophages exerted immune-suppressive effects by inhibiting the production of inflammatory factors and manifested high levels of fatty acid oxidative phosphorylation metabolism. Conversely, the main metabolism pathway for SPP1+ macrophage was glycolysis which contributed to tumor metastasis by promoting angiogenesis and matrix remodeling. In terms of the differentially expressed genes, the complement gene C1QC and the matrix remodeling relevant genes FN1 and SPP1 were differentially expressed in the TAMs between sexes, of which the male upregulated SPP1 showed the potential as an ideal target for adjuvant immunotherapy and improving the efficacy of immunotherapy. According to the early-stage TCGA-NSCLC cohort, high expression of the above three genes in immune cells were associated with poor prognosis and acted as independent prognostic factors. Moreover, through verification at the transcription factor, transcriptome, and protein levels, we found that TAMs from women showed stronger immunogenicity with higher interferon-producing and antigen-presenting ability, while men-derived TAMs upregulated the PPARs and matrix remodeling related pathways, thus were more inclined to be immunosuppressive. Deconstruction of the TAMs at the single-cell level deepens our understanding of the mechanism for tumor occurrence and progress, which could be helpful to achieve the precise sex-specific tumor treatment sooner.

Yang Q ，Zhang H ，Wei T ，Lin A ，Sun Y ，Luo P ，Zhang J ... -  《Frontiers in Immunology》

Spatial interaction and functional status of CD68(+)SHP2(+) macrophages in tumor microenvironment correlate with overall survival of NSCLC.

Tumor-associated macrophages (TAMs) constitute a plastic and heterogeneous cell population of the tumor microenvironment (TME) that can regulate tumor proliferation and support resistance to therapy, constituting promising targets for the development of novel anticancer agents. Our previous results suggest that SHP2 plays a crucial role in reprogramming the phenotype of TAMs. Thus, we hypothesized that SHP2+ TAM may predict the treatment efficacy of non-small cell lung cancer NSCLC patients as a biomarker. We analyzed cancer tissue samples from 79 NSCLC patients using multiplex fluorescence (mIF) staining to visualize various SHP-2+ TAM subpopulations (CD68+SHP2+, CD68+CD86+, CD68 + 206+, CD68+ CD86+SHP2+, CD68+ CD206+SHP2+) and T cells (CD8+ Granzyme B +) of immune cells. The immune cells proportions were quantified in the tumor regions (Tumor) and stromal regions (Stroma), as well as in the overall tumor microenvironment (Tumor and Stroma, TME). The analysis endpoint was overall survival (OS), correlating them with levels of cell infiltration or effective density. Cox regression was used to evaluate the associations between immune cell subsets infiltration and OS. Correlations between different immune cell subsets were examined by Spearman's tests. In NSCLC, the distribution of different macrophage subsets within the TME, tumor regions, and stroma regions exhibited inconsistency. The proportions of CD68+ SHP2+ TAMs (P < 0.05) were higher in tumor than in stroma. And the high infiltration of CD68+SHP2+ TAMs in tumor areas correlated with poor OS (P < 0.05). We found that the expression level of SHP2 was higher in M2-like macrophages than in M1-like macrophages. The CD68+SHP2+ subset proportion was positively correlated with the CD68+CD206+ subset within TME (P < 0.0001), tumor (P < 0.0001) and stroma (P < 0.0001). The high infiltration of CD68+SHP2+ TAMs predict poor OS in NSCLC. Targeting SHP2 is a potentially effective strategy to inhibit M2-phenotype polarization. And it provides a new thought for SHP2 targeted cancer immunotherapy.

Liu X ，Zhang Z ，Yuan J ，Yu J ，Chen D ... -  《Frontiers in Immunology》