Targeting inhibition of prognosis-related lipid metabolism genes including CYP19A1 enhances immunotherapeutic response in colon cancer
Journal of Experimental & Clinical Cancer Research volume 42, Article number: 85 (2023)
Lipid metabolic reprogramming in colon cancer shows a potential impact on tumor immune microenvironment and is associated with response to immunotherapy. Therefore, this study aimed to develop a lipid metabolism-related prognostic risk score (LMrisk) to provide new biomarkers and combination therapy strategies for colon cancer immunotherapy.
Differentially expressed lipid metabolism-related genes (LMGs) including cytochrome P450 (CYP) 19A1 were screened to construct LMrisk in TCGA colon cancer cohort. The LMrisk was then validated in three GEO datasets. The differences of immune cell infiltration and immunotherapy response between LMrisk subgroups were investigated via bioinformatic analysis. These results were comfirmed by in vitro coculture of colon cancer cells with peripheral blood mononuclear cells, human colon cancer tissue microarray analysis, multiplex immunofluorescence staining and mouse xenograft models of colon cancer.
Six LMGs including CYP19A1, ALOXE3, FABP4, LRP2, SLCO1A2 and PPARGC1A were selected to establish the LMrisk. The LMrisk was positively correlated with the abundance of macrophages, carcinoma-associated fibroblasts (CAFs), endothelial cells and the levels of biomarkers for immunotherapeutic response including programmed cell death ligand 1 (PD-L1) expression, tumor mutation burden and microsatellite instability, but negatively correlated with CD8+ T cell infiltration levels. CYP19A1 protein expression was an independent prognostic factor, and positively correlated with PD-L1 expression in human colon cancer tissues. Multiplex immunofluorescence analyses revealed that CYP19A1 protein expression was negatively correlated with CD8+ T cell infiltration, but positively correlated with the levels of tumor-associated macrophages, CAFs and endothelial cells. Importantly, CYP19A1 inhibition downregulated PD-L1, IL-6 and TGF-β levels through GPR30-AKT signaling, thereby enhancing CD8+ T cell-mediated antitumor immune response in vitro co-culture studies. CYP19A1 inhibition by letrozole or siRNA strengthened the anti-tumor immune response of CD8+ T cells, induced normalization of tumor blood vessels, and enhanced the efficacy of anti-PD-1 therapy in orthotopic and subcutaneous mouse colon cancer models.
A risk model based on lipid metabolism-related genes may predict prognosis and immunotherapeutic response in colon cancer. CYP19A1-catalyzed estrogen biosynthesis promotes vascular abnormality and inhibits CD8+ T cell function through the upregulation of PD-L1, IL-6 and TGF-β via GPR30-AKT signaling. CYP19A1 inhibition combined with PD-1 blockade represents a promising therapeutic strategy for colon cancer immunotherapy.
Colon cancer is one of the most common human malignancies in the world . Despite the fact that clinical treatment for colon cancer has been improved with the development of immune checkpoint inhibitors (ICIs), the majority of colon cancer patients have limited response to ICI therapies . Emerging biomarkers such as tumor mutational burden (TMB), inflammatory tumor microenvironment (TME) and microsatellite instability (MSI) have been identified to predict therapeutic benefit in colon cancer . Unfortunately, drawbacks still remained in current biomarkers . Therefore, it is important to excavate the predictive biomarkers for immunotherapy response and find a novel strategy for sensitizing ICIs in colon cancer.
Lipid metabolic reprogramming promotes tumor growth, angiogenesis and metastasis . A prognostic signature of nine lipid metabolism-related genes (LMGs), including CDIPT, MTMR7, PIK3CB, PIK3C2G, ARSE, ARSJ, GLA, GLB1 and UGCG, has been established in diffuse gliomas . The model including the four LMGs (ABCA1, ACSL1, AGPAT1 and SCD) is proposed as a prognostic marker of colon cancer with stage II , but not effective in all stages of colon cancer . A very recent study has shown an 8-gene prognostic signature based on LMGs in colon adenocarcinoma . However, the prognostic prediction model of multiple LMGs in colon cancer has only begun to be appreciated.
Metabolic reprogramming of tumor cells induces metabolic stress in tumor-infiltrating immune cells and stromal cells, and thereby impairs antitumor immune responses . Targeted reprogramming of lipid metabolism inhibits tumor cell growth, alleviates the immunosuppressive TME and improves response to ICI therapy . A study has shown that cyclooxygenase (COX) enzyme inhibitor aspirin in combination with anti-programmed cell death protein 1 (PD-1) treatment exhibits a synergistic effect in reducing tumor growth . Our recent study has also shown that inhibition of COX-2 catalyzed metabolism of arachidonic acid (AA) by melafolone promotes anti-PD-1 therapy in lung cancer through PD-L1 downregulation . Avasimibe, an inhibitor of cholesterol esterification enzyme, significantly empowers the anti-tumor response of CD8+ T cells, and its combination with anti-PD-1 antibody has a better efficacy in melanoma . Some clinical studies have further confirmed that targeted metabolic reprogramming enhances the anti-tumor efficacy of ICIs . Therefore, it is essential to investigate the new strategy of targeting tumor lipid metabolism to alleviate the immunosuppressive TME and to enhance anti-tumor immunotherapy.
Cytochrome P450 (CYP) 19A1 encodes aromatase, an isoenzyme of estrogen biosynthesis, and is overexpressed in colon cancer tissues . Aromatase inhibitors, including letrozole, anastrozole and exemestane, inhibit the synthesis of estrogen from androgen by binding to the aromatase, and thereby block tumor cell proliferation . Aromatase inhibitors are widely used to treat post-menopausal hormone-sensitive breast cancer . Exemestane exhibits great growth inhibitory potentials in gastric cancer when administered in combination therapy with 5-fluorouracil . Letrozole plus mTOR inhibitor everolimus results in a better progression free survival in patients with breast cancer, accompanied by decreases in Ki-67 index and tumor-infiltrating regulatory T cells, and an increase in tumor-specific CD8+ T cells [18,19,20]. Recently, the combination of exemestane and CTLA-4 monoclonal antibody tremelimumab in breast cancer has entered a phase II clinical trial . However, the effect of aromatase inhibitors on the immunotherapeutic response of colon cancer is not well understood.
In this study, we utilized The Cancer Genome Atlas (TCGA) cohort as a training set to develop a lipid metabolism-related gene prognostic risk score (LMrisk) in colon cancer, and the significant prognostic values of the model were then validated in three Gene Expression Omnibus (GEO) testing sets. The relationships of LMrisk with 36 immune cell signatures and predictive biomarkers for immunotherapeutic response were evaluated. In particular, the bioinformatic findings were confirmed using human colon cancer tissue samples, in vitro coculture of colon cancer cells with peripheral blood mononuclear cells (PBMCs), and mouse xenograft models of colon cancer. Our study develops a novel lipid metabolism-related risk model to predict immunotherapeutic response, and elucidates the molecular mechanism by which CYP19A1-catalyzed estrogen synthesis mediates immune escape, providing new targets and candidates for sensitizing colon cancer immunotherapy.
Data and resources
The RNA-sequencing data and clinical information of colon cancer patients were downloaded from the TCGA database (https://portal.gdc.cancer.gov). The RNA expression profiles contained 494 samples. We obtained 453 colon cancer samples and 41 adjacent normal tissue samples. 435 colon cancer patients with survival data were included in the following study. Gene expression profiles of GSE41258 dataset based on the Affymetrix human genome U133A array platform, GSE38832 and GSE39582 datasets based on the Affymetrix human genome U133 Plus 2.0 array platform and clinical data were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). A list of LMGs was collected from the “metabolism of lipids” in Reactome (https://reactome.org/download-data/). Genes not included in TCGA or GEO databases were excluded, and 731 genes related to lipid metabolism were obtained.
Differential expression analysis of LMGs and functional annotation
The “edgeR”  package was used to identify differentially expressed LMGs between the tumor and adjacent normal tissue samples. Adjusted P-value < 0.05 and |log2 (fold change) |> 1 were chosen as the cut-off threshold. The protein–protein interaction network of differentially expressed LMGs was analyzed by STRING database (https://string-db.org/). Principal component analysis (PCA) was utilized to analyze the expression pattern in the colon cancer and normal tissues. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed with the differentially expressed LMGs by using the “clusterProfiler” R package . False discovery rate < 0.05 was considered statistically significant.
Construction and verification of the LMrisk
Univariate Cox regression analysis was performed to identify the prognosis-related LMGs. The P-value < 0.05 in univariate Cox regression analysis was considered statistically significant. The “glmnet” R package was used to perform a least absolute shrinkage and selection operator (LASSO)-Cox regression model analysis. The weighted LASSO-Cox coefficients based on individual gene expression levels were used to calculate the lipid metabolism-related risk score (LMrisk) as follows: LMrisk = ∑ (expression of genei × Coefficient of genei). The patients in the TCGA database were stratified into the low- and high-LMrisk groups according to median LMrisk value, and their survival were analyzed using the Kaplan–Meier method. Log-rank test was used to compare the survival curves of two or more groups. The specificity and sensitivity of the LMrisk in predicting 3-, 5- and 10-year survival were determined by receiver operating characteristic (ROC) analysis using the “survivalROC” R package, and the areas under curves (AUC) were calculated. We used a similar approach in the GSE41258, GSE38832 and GSE39582 datasets to verify the applicability of the LMrisk. To study whether the LMrisk is an independent predictor for overall survival of colon cancer patients, univariate and multivariate Cox regression analyses were conducted. The LMrisk, age, gender and TNM stage were used as covariates.
Construction and verification of prognostic nomogram
After testing for collinearity, all independent prognostic parameters were included in the construction of a nomogram to predict 3- and 5-year overall survival of colon cancer patients. Age, TNM stage and LMrisk were used to construct the nomogram using the “rms” and “survival” packages in R. The discrimination power of the predictive model was evaluated by Harrell’s concordance index (C-index). We calculated the C-index with 95% confidence interval using the bootstrap approach with 1000 resamples. Then, calibration curves were drawn to assess the consistency between actual and predicted survival. The nomogram performance was validated by ROC curves at 3 and 5 year using the TNM stage as control.
Analysis of tumor immune signatures and function enrichment for LMrisk
The Estimation of Stromal and Immune cells in Malignant Tumors using Expression data (ESTIMATE) was used to evaluate the immune score, stromal score and ESTIMATE score of each sample . Based on the gene expression data in the cancer tissues, the xCell and TIMER algorithms were applied to estimate the infiltration of immune cells in each sample . The R package “maftools” was used to evaluate and sum the mutation data.
Chemicals and reagents
Letrozole (T1590) was purchased from TargetMol (Bellingham, WA, USA). CYP19A1 siRNA (sc-41498) and CYP19A1 antibody (sc-374176) were purchased from Santa Cruz Biotechnology (Santa Cruz, CA, USA). PD-L1 (13684) antibody was purchased from Cell Signaling Technology (Danvers, MA, USA). PE-conjugated anti-CD3 antibody (555340), APC-conjugated anti-CD8 antibody (566852), PE-Cy7-conjugated anti-IFN-γ antibody (557643) were purchased from BD Biosciences (San Jose, CA, USA). FITC-conjugated anti-CD107a antibody (328606) and Zombie NIR (423105) were purchased from BioLegend (San Diego, CA, USA). Opal Polaris 7 Color IHC Detection Kit (NEL861001KT) was purchased from Akoya Bioscience (Menlo Park, CA, USA). Chitosan (CS, deacetylation 98%, Mw = 50 KDa) and sodium tripolyphosphate (TPP) were purchased from Sigma Chemical Co (St. Louis, MO, USA). Hyaluronic acid (HA) was purchased from Meryer Chemical Technology Co Ltd (Shanghai, China).
Human tissue microarray analyses for CYP19A1 and PD-L1 expression
Human colon carcinoma tissue microarray (HCol-Ade180Sur-05) containing 90 patients with colon cancer and adjacent normal tissues were obtained from Shanghai Outdo Biotech Co., Ltd. Each sample dot with a diameter of 1.5 mm and a thickness of 4 μm was prepared according to a standard method. Immunohistochemistry (IHC) of human colon cancer tissue microarrays were performed by using the CYP19A1 (1:100) and PD-L1 (1:300) antibodies. Immunostaining was graded using a two-score system based on intensity score and proportion score as previously described . The two scores were then multiplied to yield a total immunoreactivity score regarding the protein expression in a sample. IHC score was assessed independently by two pathologists, and a consensus of grading was reached. Patients were divided into high- and low-expression groups according to the median IHC score.
Multiplex immunofluorescence analysis
Tumor tissues from 20 patients with colon cancer were obtained from Zhongnan Hospital, Wuhan University (Hubei, China) between 2020 and 2021. The histological diagnosis for each sample was reconfirmed using microscopic examination of hematoxylin/eosin-stained sections. All samples were collected from patients with informed consent, and all procedures were conducted with the approval of the Ethical Committee of the Medical School of Wuhan University and performed in accordance with relevant regulations and guidelines.
Quantitative multiplex immunofluorescence (mIF) was performed to characterize the immune landscape in human colon cancer tissues using Opal Polaris 7 Color IHC Detection Kits (NEL861001KT, Akoya Bioscience, CA, USA) as previously described . Briefly, formalin-fixed paraffin-embedded sections were deparaffinized, followed by antigen retrieval with citrate acid buffer (pH 6.0)/Tris–EDTA buffer (pH 9.0) and then blocking with blocking/Ab diluent. Next, slides were incubated with primary antibodies against CYP19A1 (1:200), CD68 (Abcam, ab955, 1:800), CD31 (Abcam, ab24590, 1:300), α-SMA (Abcam, ab7817, 1:300), CD8 (DAKO, IR623, 1:5) and CK (AE1/AE3) (DAKO, IS05330-2, 1:3). Primary antibody was visualized using tyramide signal amplification linked to a specific fluorochrome from the mIF kit for each primary antibody. A stripping procedure based on microwaves was performed for each consecutive antibody staining. The stained slides were scanned using a Vectra® 3 multispectral microscope (Akoya Bioscience). From each slide, Vectra automatically captured the fluorescent spectra from 420 to 720 nm at 20-nm intervals with the same exposure time and then combined the captured images to create a single stack image that retained the particulate spectral signature of all markers. These data were analyzed using InForm 2.6 software.
Cell culture, cell proliferation assays, aromatase activity assays, CYP19A1 siRNA transfection, isolation of PBMCs and their coculture with colon cancer cells, flow cytometry analyses for tumor cell death and CD8+ T cell function, western blot analyses for CYP19A1 and PD-L1 expression, tumor models and therapeutic efficacy study, and statistical analysis were described in Supplementary Materials and Methods.
Differential expression analysis of LMGs and functional annotation
Expression data and clinical information of four cohorts were downloaded from the TCGA and GEO databases. The TCGA cohort consisted of 453 colon cancer and 41 normal samples, and was used for differential expression analysis. Seven hundred thirty-one LMGs were screened, and 178 LMGs were differentially expressed with 86 upregulation and 92 downregulation between colon cancer and normal samples (Fig. 1A). PCA revealed an obvious difference in a distinct cluster of these differentially expressed LMGs between the colon cancer tissues and adjacent normal tissues (Fig. 1B). The hub genes including CYP19A1 were identified by the protein–protein interaction network of differentially expressed LMGs in the colon cancer (Fig. 1C). Next, GO and KEGG pathway annotation analyses were performed to explore the biological functions of the differentially expressed LMGs in the colon cancer. GO analysis showed that these LMGs were enriched in the basic biological processes, including steroid metabolism, fatty acid metabolism, lipid catabolism and glycerolipid metabolism, and the most highly enriched terms for the cellular component and molecular function were lipid droplet and oxidoreductase activity, respectively (Fig. 1D). KEGG analysis showed that 178 LMGs were primarily related to AA metabolism, PPAR signaling pathway and bile secretion (Fig. 1E). These results suggest significant alterations of LMGs in human colon cancer.
Construction and verification of the LMrisk
To develop prognostic models in the TCGA colon cancer cohort, machine learning methods such as random survival forest, CoxBoost, support vector machine, gradient boosting machine, elastic net and LASSO-Cox were employed. A univariate combined LASSO-Cox regression model with the maximum C-index (Fig. 2A) was used for subsequent analysis. Univariate Cox regression analysis revealed that the LRP2, CYP19A1, SLCO1A2, FABP4, ALOXE3, PLAAT5 and SPHK1 were risk genes with hazard ratio (HR) > 1, while PPARGC1A was protective gene with HR < 1 among 178 differentially expressed LMGs (Fig. 2B). By using LASSO-Cox regression analysis to further reduce overfitting, the six LMGs, including LRP2, CYP19A1, SLCO1A2, FABP4, ALOXE3 and PPARGC1A, were identified to construct the optimal lipid metabolism-related prognostic signature (Fig. 2C and D). The expression levels and LASSO coefficients of six LMGs were extracted to calculate the LMrisk for each patient with the following formula: LMrisk = (0.124 × LRP2 expression) + (0.145 × CYP19A1 expression) + (0.0551 × SLCO1A2 expression) + (0.00237 × FABP4 expression) + (0.117 × ALOXE3 expression) + (-0.182 × PPARGC1A expression). Based on the median LMrisk, the patients were assigned to low LMrisk (n = 217) and high LMrisk groups (n = 218). We observed more deaths and shorter survival time in the colon cancer patients with the high LMrisk as compared with the low LMrisk (Fig. 2E and F). Next, time-dependent ROC curves were used to evaluate the predictive value of the LMrisk. As shown in Fig. 2G, the AUC of LMrisk at 3, 5 and 10 years were 0.68, 0.73 and 0.81, respectively. Consistently, Kaplan–Meier survival curves indicated that the survival time was shorter in the high LMrisk group than in the low LMrisk group in the GSE41258 (P = 0.0072), GSE38832 (P = 0.0099) and GSE39582 (P = 0.022) datasets (Fig. 2H–J). To validate the predictive ability of the LMrisk on colon cancer prognosis, we conducted a stratified analysis based on clinicopathological features including age, gender, TNM stage, N stage, M stage and T stage in the high and low LMrisk groups. Kaplan–Meier survival analyses revealed that the high-risk group had worse overall survival compared to the low LMrisk group in different strata of clinical characteristics including younger (< 66 years) or older (≥ 66 years), male or female, stage I–II or stage III–IV, N0 or N1–N2, M0 or M1 and T3–T4 patients (Fig. S1). These results suggest that the LMrisk has good robustness for predicting prognosis of colon cancer.
To identify the factors affecting the survival of patients with colon cancer, we analyzed the prognostic value of the LMrisk and clinicopathological parameters in the TCGA dataset. As shown in Fig. S2A, the LMrisk, in addition to clinical factors (age, T stage, N stage, M stage and TNM stage), was closely associated with survival outcomes. Multivariate analysis manifested that the LMrisk, age and TNM stage were independent prognostic indicators (Fig. S2B). The AUC of LMrisk was larger than other indicators including TNM stage, suggesting that the LMrisk has a better ability to predict prognosis in the colon cancer (Fig. S2C). The LMrisk was significantly correlated with the age, T stage, N stage, M stage and TNM stage. We also observed relatively low expression of PPARGC1A and high expression of CYP19A1, FABP4, ALOXE3, LRP2 and SLCO1A2 in the patients with high LMrisk (Fig. S2D). The patients had higher LMrisk in stage III–IV than in stage I–II (P = 0.023), and in T3–T4 than in T1–T2 (P = 0.016). The patients with lymph node metastasis (N1–N2) or distant metastasis (M1) had higher LMrisk than those without lymph node metastasis (N0) or distant metastasis (M0) (Fig. S2E). These data suggest that the LMrisk is an independent prognostic risk factor in the colon cancer.
Construction and validation of a prognostic nomogram based on LMrisk
To predict overall survival probabilities of colon cancer patients, a prognostic nomogram was established by integrating the LMrisk with age and TNM stage in the TCGA dataset (Fig. S3A). We used calibration curve to verify the accuracy of the prediction model, and found that the predicted probability of 3- and 5-year survival fitted well with the observed survival (Fig. S3B and C). The AUC values of nomogram were higher than those of TNM stage in 3 and 5 year (Fig. S3D and E). The nomogram also showed satisfactory discrimination and calibration with a C-index of 0.79 (95% confidence interval: 0.74–0.84). These data suggest that our nomogram presents superior prognostic value than that of TNM stage in colon cancer.
The association of LMrisk with immune cell infiltration
The tumor immune microenvironment innately modulates tumor progression . Thus, the ESTIMATE algorithm was used to evaluate the immune score and stromal score of colon cancer patients. Figure 3A showed that the LMrisk was positively correlated with the immune score, stroma score and ESTIMATE score in the TCGA dataset. Next, immune cell and stromal cell infiltration in colon cancer tissues were inferred by xCell and EPIC algorithms. We observed less infiltration of CD8+ T cells, but more proportions of immunosuppressive cells including macrophages, monocytes, CAFs and endothelial cells in the high LMrisk group compared to the low LMrisk group (Fig. 3B). These data confirm that the infiltration of CD8+ T cells is decreased but the infiltration of immunosuppressive cells is increased in the colon cancer with high LMrisk.
Association of the LMrisk with immunotherapy response in colon cancer
The TME includes immune checkpoint regulators as well as inflammatory mediators besides immune cells . To understand the molecular mechanisms altering immunotherapy responsiveness, we evaluated the association of the LMrisk with immunotherapy response in colon cancer. As shown in Fig. 4A, PD-1, PD-L1, CTLA-4, LAG3 and HAVCR2 (TIM3) were significantly overexpressed in the high LMrisk group than in the low LMrisk group. The patients with low LMrisk and high PD-L1 had better survival than those with high LMrisk and high PD-L1, and patients with low LMrisk and low PD-L1 had better survival than those with high LMrisk and low PD-L1 (Fig. 4B). Similar results were obtained in the PD-1 and CTLA-4 stratifying groups (Fig. 4C and D). MSI status, consensus molecular subtypes (CMS) heterogeneity and TMB display different immune landscapes in the TME, and therefore are used as predictive biomarkers for response of ICIs [29,30,31,32]. Accordingly, we compared the LMrisk between different MSI status, and found that the LMrisk was higher in the colon cancer with MSI-H and MSI-L as compared with MSS (Fig. 4E). The LMrisk was higher in CMS1 (immune activated phenotype) or CMS4 (immune inflamed phenotype) than in CMS2 (immune desert phenotype) and CMS3 (immune excluded phenotype) (Fig. 4F). A previous study shows that, in six immune subtypes, IFN-γ dominant subtype has better ICI outcomes [33, 34]. We observed the significantly increased LMrisk in IFN-γ dominant subtype than in wound healing subtype, though there was no significant difference of LMrisk between inflammatory and lymphocyte depleted subtypes (Fig. 4G). The TMB was significantly increased in the patients with high LMrisk as compared with low LMrisk (Fig. 4H). Maftools analysis results showed that TTN gene mutated more frequently in the high LMrisk than in the low LMrisk group (Fig. 4I). Because of unavailable datasets for colon cancer patients treated with immune checkpoint blockade, we analyzed publicly available datasets from urothelial cancer patients treated with anti-PD-1 therapy. As shown in Fig. 4J, the high LMrisk was correlated with a complete or partial response (CR/PR), while low LMrisk correlated with stable disease or progressive disease (SD/PD) when the patients received anti-PD-1 therapies. The high LMrisk also predicted significant improvement in overall survival for the patients received anti-PD-1 therapy compared with the low LMrisk. These results indicate that the LMrisk is a potential predictive biomarker of immunotherapeutic response.
CYP19A1, a risk factor of the LMrisk model, is associated with immune signature
Among the differentially expressed lipid metabolism genes in colon cancer and adjacent normal tissues, CYP19A1 is a key molecular node and has the largest prognostic hazard ratio (Figs. 1C and 2B). Therefore, CYP19A1, a risk factor of the LMrisk model, was used for subsequently experimental verification. We found that CYP19A1 gene was highly expressed in the colon cancer tissues compared to the adjacent normal tissues in the Gene Expression Profiling Interactive Analysis (GEPIA) webserver (Fig. S4A). CYP19A1 gene expression was a strong prognostic risk factor in the colon cancer (Fig. S4B), and positively correlated with PD-L1 expression in the GEPIA webserver (Fig. S4C). Subsequently, we conducted functional verification by using a tissue microarray of human colon cancer. Similarly, CYP19A1 protein expression was significantly increased in human colon cancer tissues (Fig. 5A and B), and positively correlated with PD-L1 expression (Fig. 5D). CYP19A1 protein expression in the colon cancer tissues was negatively correlated with overall survival (Fig. 5C), and positively associated with N stage, and an independent prognostic parameter of overall survival (Tables S1 and S2). Next, mIF staining was performed to explore the relationship between CYP19A1 expression and immune landscape in the colon cancer. We observed an apparent colocalization of CYP19A1 with CK (tumor cell marker), indicating that CYP19A1 is mainly expressed in tumor cells. The infiltration of CD8+ T cells was more but the proportions of CD68+ cells (macrophages), α-SMA+ cells (CAFs) and CD31+ cells (endothelial cells) were fewer in human colon cancer tissues with low CYP19A1 expression as compared with high CYP19A1 expression (Fig. 5E). We analyzed PD-L1 expression in the same biopsies, and found that PD-L1+ CAFs and macrophages were significantly higher in CYP19A1 high expression group as compared with low CYP19A1 expression group (Fig. 5F). These data indicate that high CYP19A1 expression leads to poor prognosis, accompanied by the low infiltration level of CD8+ T cells and high infiltration levels of immunosuppressive cells in colon cancer.
CYP19A1 inhibition potentiates CD8+ T cell-mediated anti-tumor immune response in vitro
To gain insight into the role of CYP19A1 in tumor immunity, the inhibitor (letrozole) and siRNAs of CYP19A1 were used in the in vitro co-culture of human HT29 and HCT116 cells with PBMCs. We found that the letrozole (5 μM) significantly decreased aromatase activity in human HT29 and HCT116 cells (Fig. 6A) without affecting their proliferation (Fig. 6B). The CYP19A1 siRNAs significantly inhibited CYP19A1 expression (Fig. 6C). Unsurprisingly, CYP19A1 inhibition by the letrozole or siRNA in tumors increased the proportion of 7-AAD+ tumor cells after co-culture with PBMCs (1.5- and 2-fold, respectively) (Fig. 6D). The letrozole or CYP19A1 siRNA-treated HT29 and HCT116 cells significantly increased CD8+ T cell proliferation and the proportion of CD8+ CD107a+ and CD8+ IFN-γ+ T cells (Fig. 6E–G), and promoted pericyte cell migration and reduced endothelial cell migration (Fig. 6H). These data indicate that CYP19A1 inhibition in colon cancer cells enhances CD8+ T cell-mediated antitumor immunity.
CYP19A1 inhibits CD8+ T cells through upregulation of PD-L1, IL-6 and TGF-β
Tumor cells upregulate the expression of immune checkpoint molecules or secrete cytokines to induce abnormal angiogenesis, thereby promoting tumor immune escape [35, 36]. We found that CYP19A1 expression level was significantly positively correlated with PD-L1, IL-6 and TGF-β levels in the TCGA colon cancer dataset (Fig. 7A). Next, we verified the factors derived from tumor cells in HT29 and HCT116 colon cancer cells, and comfirmed that CYP19A1 inhibition significantly reduced PD-L1, IL-6, and TGF-β expression in tumor cells (Fig. 7B–D). Neutralizing antibodies against PD-L1, IL-6 or TGF-β partially ameliorated the effects of CYP19A1 overexpression on the proliferation and effector function of CD8+ T cells (Fig. 7E and F). Neutralizing IL-6 and TGF-β also promoted pericyte cell migration and inhibited endothelial cell migration (Fig. 7G). Importantly, combined blockade of PD-L1, IL-6 and TGF-β exerted better efficacy than single blockade (Fig. 7E–G). These data highlight compelling evidence that CYP19A1 inhibits the proliferation and effector function of CD8+ T cells by upregulating PD-L1, IL-6 and TGF-β expression in colon cancer cells.
CYP19A1/estradiol mediates immunosuppression via GPR30-Akt signaling
Activation of G-protein coupled estrogen receptor 30 (GPR30) promotes colon cancer growth through mediating estrogenic activity in colon cancer cells , and Akt signaling is one of the key downstream pathways of GPR30 in cancer progression . We discovered that CYP19A1 inhibition by letrozole or si-CYP19A1 downregulated GPR30, p-Akt and PD-L1 in HCT116 and HT29 colon cancer cells (Fig. 8A), which were greatly enhanced by exogenous addition of estradiol (Fig. 8B). GPR30 knockdown decreased p-Akt and PD-L1 levels (Fig. 8C and D), and Akt knockdown also inhibited PD-L1 protein expression in HCT116 and HT29 colon cancer cells (Fig. 8E and F). Exogenous supplementation of estradiol increased IL-6 and TGF-β expression in HCT116 and HT29 colon cancer cells (Fig. 8G). Importantly, GPR30 or Akt siRNA reversed the increased expression of IL-6 and TGF-β in HCT116 and HT29 colon cancer cells in response to estradiol stimulation (Fig. 8G). These results imply that GPR30-Akt signaling is crucial for CYP19A1/estradiol-mediated immunosuppression in colon cancer.
Nanoparticle-encapsulated CYP19A1 siRNA synergizes with anti-PD-1 therapy
The established hyaluronic acid-modified chitosan nanoparticles are recently shown to be effective for the systemic administration of siRNA to tumors [39, 40]. To explore the possibility of using CYP19A1 as a target for sensitizing anti-PD-1 therapy, the effect of nanoparticle-encapsulated CYP19A1 siRNA on PD-1 blockade therapy was investigated in the orthotopic and subcutaneous colon cancer models. We detected obvious decreases in CYP19A1 protein expression and estradiol production in siCYP19A1-treated MC38 tumor as compared with the si-Scram (Fig. 9A-C). Knockdown of CYP19A1 by specific siRNA significantly impeded growth of orthotopic and subcutaneous MC38 colon cancer, but its combination with anti-PD-1 treatment was far more effective than monotherapy (Fig. 9D-H). CYP19A1 knockdown boosted tumor-infiltrating CD8+ T cells and enhanced GzmB and IFNγ production, indicative of CD8+ T cells-elicited adaptive immunity against colon cancer, which were amplyfied in combination therapy (Fig. 9I). Additionally, the CYP19A1 siRNA remarkably inhibited PD-L1 expression on tumor cells (Fig. 9J), induced vascular normalization, as manifested by an obvious decrease in the hypoxia-inducible factor (HIF) -1α level and an increase in CD31+ αSMA+ cells number (Fig. 9K and L). Together, our results demonstrate that CYP19A1 inhibition drastically enhances anti-PD-1 therapy for colon cancer.
Letrozole facilitates anti-PD-1 therapy by promoting CD8+ T cell-mediated anti-tumor immune response in murine models
Next, we further explored the possibility of using CYP19A1 as a target for sensitizing anti-PD-1 therapy, and found that the CYP19A1 inhibitor letrozole sensitized anti-PD-1 therapy in subcutaneous tumors derived from the MC38 and CT26 murine colon cancer (Fig. 10A and B). The combination of letrozole with PD-1 blockade promoted expression of IFNγ and CD107a and intratumoral infiltration of CD8+ T cells compared to the monotherapy (Fig. 10C and D). The combination treatment reduced GPR30 and p-Akt levels in MC38 and CT26 colon cancer tissues with an obvious decrease in estradiol production compared to the anti-PD-1 monotherapy (Fig. 10E and F). The letrozole inhibited PD-L1 expression in tumor cells (Fig. 10G), boosted tumor perfusion, promoted vascular normalization, and downregulated HIF-1α in CT26 and MC38 tumor (Fig. 10H–J). Similar results were obtained in the mice model of orthotopic MC38 colon cancer (Supplementary Fig. S5). The combination of letrozole with anti-PD-1 treatment did not result in any extra toxicity compared with the anti-PD-1 monotherapy (Supplementary Table S3). Overall, these findings strongly indicate that CYP19A1 inhibition greatly facilitates anti-PD-1 therapy for colon cancer.
Previous studies have shown that the alterations in lipid metabolism including cholesterol metabolism and AA metabolism affect immunotherapeutic response, and are promising biomarkers to predict the efficacy of immunotherapy [41, 42]. In the present study, two novel observations have been made. First, we used differentially expressed LMGs, including CYP19A1, FABP4, LRP2, SLCO1A2, PPARGC1A and ALOXE3, to construct the LMrisk for predicting prognosis of colon cancer, and found that the LMrisk was associated with the immunosuppressive TME and predictive biomarkers of immunotherapeutic response in colon cancer. To our knowledge, this is the first study that directly demonstrated that a prognostic risk model based on LMGs may predict immunotherapeutic response in colon cancer. Second, we uncovered for the first time that CYP19A1 inhibition downregulated PD-L1, IL-6 and TGF-β expression in colon cancer cells, and thereby enhanced the tumor-killing ability of CD8+ T cells. These results implicate the LMrisk based on the prognostic model as a predictive biomarker of immunotherapeutic response in colon cancer, highlighting its therapeutic potential for optimizing anti-PD-1 immunotherapy.
ALOXE3, PPARGC1A or FABP4 are prognostic biomarkers in colon cancer [43,44,45]. However, a single gene is difficult to provide powerful predictive performance for colon cancer patients . Therefore, it is a trend to use a multiple gene model to predict prognosis of cancer patients . Indeed, several studies have used multiple LMGs to construct prognostic models in patients with breast cancer, gastric cancer or osteosarcoma [48,49,50]. Here we integrated CYP19A1, FABP4, LRP2, SLCO1A2, PPARGC1A and ALOXE3 with clinicopathological information of patients to construct a prognostic nomogram in colon cancer. A recent study shows an 8-gene signature based on LMGs including RTN2, FYN, HEYL, FAM69A, FBXL5, HMGN2, LGALS4, STOX1 as a novel marker to predict colon cancer patients’ survival . We reason that this discrepancy in prognosis-related LMGs could be due to different algorithms (differentially expressed LMGs vs. all LMGs). Our lipid metabolism-related nomogram has larger C-index than the model of Jiang et al. , suggesting its ability to better predict colon cancer prognosis as compared with gene signature-derived risk score. A previous study shows that nomogram-derived prognosis, together with user-friendly digital interfaces, elevates the accuracy of cancer survival prediction, and thereby allows for seamless incorporation to aid clinical decision making .
High level of MSI is closely related to the improved prognosis of colon cancer patients receiving immunotherapy , and higher TMB results in more neo-antigens, increasing chances for T cell recognition, accompanied by better ICI outcomes . Here we demonstrated that the LMrisk was positively correlated with TMB and MSI, suggesting that the patients with high LMrisk are more likely to benefit from immunotherapy. Further analysis revealed that TTN mutation was the major reason for the high TMB in the colon cancer patients with high LMrisk. Previous studies show that TTN mutation is related to high immunogenicity and inflammatory tumor immune microenvironment in lung adenocarcinoma, accompanied by favorable objective response and survival with ICI administration [51, 52]. Along with immunogenicity, inflammatory TME makes colon cancer patients amenable to respond to ICIs . We also observed more infiltration of macrophages, monocytes, NK cells and CAFs, and higher PD-L1 expression in colon cancer tissues from the patients with high LMrisk. Our extensive functional studies demonstrated that CYP19A1 protein expression was positively correlated with PD-L1 expression and infiltration of macrophages, CAFs and endothelial cells in human colon cancer tissues. Given that strong correlations of the LMrisk with PD-L1 expression, TMB and MSI, we conclude that the LMrisk including CYP19A1 is a promising biomarker for predicting immunotherapeutic response to colon cancer.
Several studies showed that men have a higher incidence of colon cancer than women, suggesting that estrogen may play a protective role in the development of colon cancer . However, recent epidemiological studies have shown that hormone replacement therapy in postmenopausal women does not inhibit the development of colon cancer [55, 56]. Moreover, anti-PD-1 or anti-PD-L1 treatment resulted in a trend to greater overall survival and better response rates in individual males with colon cancer as compared with females . Aromatase encoded by CYP19A1 and GPR30 are highly expressed in colon cancer tissues, and high expression of GPR30 predicts poor prognosis in colon cancer patients [58, 59]. Here, we demonstrated that CYP19A1 inhibition by letrozole or siRNA reduced production of IL-6 and TGF-β and downregulated PD-L1 expression by inactivating GPR30-Akt signaling, and thereby promoted the proliferation and cytotoxic activity of CD8+ T cells (Fig. S6). These data delineate that CYP19A1 inhibition combined with PD-1 antibody represents a promising therapeutic strategy for colon cancer.
According to our analyses, the LMrisk may predict response and outcome of immunotherapy in colon cancer. However, there are still some limitations that should be acknowledged. Firstly, because the mRNA expression data from colon cancer patients with immunotherapy is not available, the predictive ability of LMrisk for immunotherapeutic response is estimated indirectly by urothelial cancer cohorts and biomarkers. Therefore, further well-powered prospective studies are still needed. Secondly, we studied the expression and pathophysiological significance of CYP19A1 in vitro and in vivo, and further studies are needed on other LMGs including FABP4, LRP2, SLCO1A2, PPARGC1A and ALOXE3 in colon cancer.
Collectively, we screened the six LMGs including CYP19A1, FABP4, LRP2, SLCO1A2, PPARGC1A and ALOXE3, and constructed a signature to predict prognosis and immunotherapeutic response in colon cancer, which was extendedly and externally validated. Importantly, CYP19A1 inhibition improves anti-PD-1 immunotherapy for colon cancer and blunts PD-L1-induced anergy of CD8+ T cells. Our findings facilitate prediction for prognosis and immunotherapeutic response in colon cancer, and targeting lipid metabolism reprogramming in the TME may be promising strategies for synergy with anti-PD-1 treatment.
Availability of data and materials
The datasets generated or analysed during the current study are available on reasonable request.
Areas under curve
Consensus molecular subtypes
Gene Expression Omnibus
Immune checkpoint inhibitor
Kyoto Encyclopedia of Genes and Genomes
Least absolute shrinkage and selection operator
Lipid metabolism-related gene
Lipid metabolism-related gene prognostic risk score
Peripheral blood mononuclear cells
Principal component analysis
Programmed cell death protein 1
Receiver operating characteristic
The Cancer Genome Atlas
Tumor mutational burden
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.
Tolba MF. Revolutionizing the landscape of colorectal cancer treatment: the potential role of immune checkpoint inhibitors. Int J Cancer. 2020;147:2996–3006.
Gibney GT, Weiner LM, Atkins MB. Predictive biomarkers for checkpoint inhibitor-based immunotherapy. Lancet Oncol. 2016;17:e542–51.
Bian X, Liu R, Meng Y, Xing D, Xu D, Lu Z. Lipid metabolism and cancer. J Exp Med. 2021;218:e20201606.
Wu F, Zhao Z, Chai RC, Liu YQ, Li GZ, Jiang HY, et al. Prognostic power of a lipid metabolism gene panel for diffuse gliomas. J Cell Mol Med. 2019;23:7741–8.
Vargas T, Moreno-Rubio J, Herranz J, Cejas P, Molina S, González-Vallinas M, et al. ColoLipidGene: signature of lipid metabolism-related genes to predict prognosis in stage-II colon cancer patients. Oncotarget. 2015;6:7348–63.
Fernández LP, Ramos-Ruiz R, Herranz J, Martín-Hernández R, Vargas T, Mendiola M, et al. The transcriptional and mutational landscapes of lipid metabolism-related genes in colon cancer. Oncotarget. 2018;9:5919–30.
Jiang C, Liu Y, Wen S, Xu C, Gu L. In silico development and clinical validation of novel 8 gene signature based on lipid metabolism related genes in colon adenocarcinoma. Pharmacol Res. 2021;169:105644.
Hao Y, Li D, Xu Y, Ouyang J, Wang Y, Zhang Y, et al. Investigation of lipid metabolism dysregulation and the effects on immune microenvironments in pan-cancer using multiple omics data. BMC Bioinformatics. 2019;20:195.
Gottfried E, Kreutz M, Mackensen A. Tumor metabolism as modulator of immune response and tumor progression. Semin Cancer Biol. 2012;22:335–41.
Zelenay S, van der Veen AG, Bottcher JP, Snelgrove KJ, Rogers N, Acton SE, et al. Cyclooxygenase-dependent tumor growth through evasion of immunity. Cell. 2015;162:1257–70.
Tang H, Liu Y, Wang C, Zheng H, Chen Y, Liu W, et al. Inhibition of COX-2 and EGFR by melafolone improves anti-PD-1 therapy through vascular normalization and PD-L1 downregulation in lung cancer. J Pharmacol Exp Ther. 2019;368:401–13.
Yang W, Bai Y, Xiong Y, Zhang J, Chen S, Zheng X, et al. Potentiating the antitumour response of CD8(+) T cells by modulating cholesterol metabolism. Nature. 2016;531:651–5.
Weng CY, Kao CX, Chang TS, Huang YH. Immuno-metabolism: the role of cancer niche in immune checkpoint inhibitor resistance. Int J Mol Sci. 2021;22:1258.
Simpson ER, Clyne C, Rubin G, Boon WC, Robertson K, Britt K, et al. Aromatase–a brief overview. Annu Rev Physiol. 2002;64:93–127.
Molehin D, Rasha F, Rahman RL, Pruitt K. Regulation of aromatase in cancer. Mol Cell Biochem. 2021;476:2449–64.
Yang J-C, Chang N, Wu D-C, Cheng W-C, Chung W-M, Chang W-C, et al. Preclinical evaluation of exemestane as a novel chemotherapy for gastric cancer. J Cell Mol Med. 2019;23:7417–26.
Fan Y, Sun T, Shao Z, Zhang Q, Ouyang Q, Tong Z, et al. Effectiveness of adding everolimus to the first-line treatment of advanced breast cancer in premenopausal women who experienced disease progression while receiving selective estrogen receptor modulators: a phase 2 randomized clinical trial. JAMA Oncol. 2021;7:e213428.
Wu W, Chen J, Deng H, Jin L, He Z, Rao N, et al. Neoadjuvant everolimus plus letrozole versus fluorouracil, epirubicin and cyclophosphamide for ER-positive, HER2-negative breast cancer: a randomized pilot trial. BMC Cancer. 2021;21:862.
Baselga J, Semiglazov V, van Dam P, Manikhas A, Bellet M, Mayordomo J, et al. Phase II randomized study of neoadjuvant everolimus plus letrozole compared with placebo plus letrozole in patients with estrogen receptor-positive breast cancer. J Clin Oncol. 2009;27:2630–7.
Vonderheide RH, LoRusso PM, Khalil M, Gartner EM, Khaira D, Soulieres D, et al. Tremelimumab in combination with exemestane in patients with advanced breast cancer and treatment-associated modulation of inducible costimulator expression on patient T cells. Clin Cancer Res. 2010;16:3485–94.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (N Y). 2021;2:100141.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Sturm G, Finotello F, Petitprez F, Zhang JD, Baumbach J, Fridman WH, et al. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics. 2019;35:i436–45.
Wu M, Wang Y, Yang D, Gong Y, Rao F, Liu R, et al. A PLK1 kinase inhibitor enhances the chemosensitivity of cisplatin by inducing pyroptosis in oesophageal squamous cell carcinoma. EBioMedicine. 2019;41:244–55.
Chen X, Liu Y, Wang Y, Wang C, Chen X, Xiong Y, et al. CYP4F2-catalyzed metabolism of arachidonic acid promotes stromal cell-mediated immunosuppression in non-small cell lung cancer. Cancer Res. 2022;82:4016–30.
Hinshaw DC, Shevde LA. The tumor microenvironment innately modulates cancer progression. Cancer Res. 2019;79:4557–66.
Bao X, Zhang H, Wu W, Cheng S, Dai X, Zhu X, et al. Analysis of the molecular nature associated with microsatellite status in colon cancer identifies clinical implications for immunotherapy. J Immunother Cancer. 2020;8:e001437.
Samstein RM, Lee CH, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, et al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. 2019;51:202–6.
Marisa L, Blum Y, Taieb J, Ayadi M, Pilati C, Le Malicot K, et al. Intratumor CMS heterogeneity impacts patient prognosis in localized colon cancer. Clin Cancer Res. 2021;27:4768–80.
Becht E, de Reynies A, Giraldo NA, Pilati C, Buttard B, Lacroix L, et al. Immune and stromal classification of colorectal cancer is associated with molecular subtypes and relevant for precision immunotherapy. Clin Cancer Res. 2016;22:4057–66.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The immune landscape of cancer. Immunity. 2018;48:812-830.e14.
Soldevilla B, Carretero-Puche C, Gomez-Lopez G, Al-Shahrour F, Riesco MC, Gil-Calderon B, et al. The correlation between immune subtypes and consensus molecular subtypes in colorectal cancer identifies novel tumour microenvironment profiles, with prognostic and therapeutic implications. Eur J Cancer. 2019;123:118–29.
Kornepati AVR, Vadlamudi RK, Curiel TJ. Programmed death ligand 1 signals in cancer cells. Nat Rev Cancer. 2022;22:174–89.
Anderson TS, Wooster AL, Piersall SL, Okpalanwaka IF, Lowe DB. Disrupting cancer angiogenesis and immune checkpoint networks for improved tumor immunity. Semin Cancer Biol. 2022;86:981–96.
Gilligan LC, Rahman HP, Hewitt A-M, Sitch AJ, Gondal A, Arvaniti A, et al. Estrogen activation by steroid sulfatase increases colorectal cancer proliferation via gper. J Clin Endocrinol Metab. 2017;102:4435–47.
Yang W, Tan W, Zheng J, Zhang B, Li H, Li X. MEHP promotes the proliferation of cervical cancer via GPER mediated activation of Akt. Eur J Pharmacol. 2018;824:11–6.
Cao Y, Tan YF, Wong YS, Liew MWJ, Venkatraman S. Recent advances in chitosan-based carriers for gene delivery. Mar Drugs. 2019;17:381.
Zhang W, Xu W, Lan Y, He X, Liu K, Liang Y. Antitumor effect of hyaluronic-acid-modified chitosan nanoparticles loaded with siRNA for targeted therapy for non-small cell lung cancer. Int J Nanomedicine. 2019;14:5287–301.
Li X, Wenes M, Romero P, Huang SC-C, Fendt S-M, Ho P-C. Navigating metabolic pathways to enhance antitumour immunity and immunotherapy. Nat Rev Clin Oncol. 2019;16:425–41.
Zhang H, Zhao W, Li X, He Y. Cholesterol metabolism as a potential therapeutic target and a prognostic biomarker for cancer immunotherapy. Onco Targets Ther. 2021;14:3803–12.
Ruan GT, Gong YZ, Zhu LC, Gao F, Liao XW, Wang XK, et al. The perspective of diagnostic and prognostic values of lipoxygenases mRNA expression in colon adenocarcinoma. Onco Targets Ther. 2020;13:9389–405.
Uddin MN, Li M, Wang X. Identification of transcriptional signatures of colon tumor stroma by a meta-analysis. J Oncol. 2019;2019:8752862.
Lin K, Huang J, Luo H, Luo C, Zhu X, Bu F, et al. Development of a prognostic index and screening of potential biomarkers based on immunogenomic landscape analysis of colorectal cancer. Aging. 2020;12:5832–57.
Collins GS, Reitsma JB, Altman DG, Moons KGM. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ. 2015;350:g7594.
Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16:e173-180.
Ye Z, Zou S, Niu Z, Xu Z, Hu Y. A novel risk model based on lipid metabolism-associated genes predicts prognosis and indicates immune microenvironment in breast cancer. Front Cell Dev Biol. 2021;9:691676.
Wei X-L, Luo T-Q, Li J-N, Xue Z-C, Wang Y, Zhang Y, et al. Development and validation of a prognostic classifier based on lipid metabolism-related genes in gastric cancer. Front Mol Biosci. 2021;8:691143.
Qian H, Lei T, Hu Y, Lei P. Expression of lipid-metabolism genes is correlated with immune microenvironment and predicts prognosis in osteosarcoma. Front Cell Dev Biol. 2021;9:673827.
Wang Z, Wang C, Lin S, Yu X. Effect of TTN mutations on immune microenvironment and efficacy of immunotherapy in lung adenocarcinoma patients. Front Oncol. 2021;11:725292.
Su C, Wang X, Zhou J, Zhao J, Zhou F, Zhao G, et al. Titin mutation in circulatory tumor DNA is associated with efficacy to immune checkpoint blockade in advanced non-small cell lung cancer. Transl Lung Cancer Res. 2021;10:1256–65.
Sadeghi Rad H, Monkman J, Warkiani ME, Ladwa R, O’Byrne K, Rezaei N, et al. Understanding the tumor microenvironment for effective immunotherapy. Med Res Rev. 2021;41:1474–98.
Barzi A, Lenz AM, Labonte MJ, Lenz H-J. Molecular pathways: Estrogen pathway in colorectal cancer. Clin Cancer Res. 2013;19:5842–8.
Manson JE, Chlebowski RT, Stefanick ML, Aragaki AK, Rossouw JE, Prentice RL, et al. Menopausal hormone therapy and health outcomes during the intervention and extended poststopping phases of the Women’s Health Initiative randomized trials. JAMA. 2013;310:1353–68.
Lavasani S, Chlebowski RT, Prentice RL, Kato I, Wactawski-Wende J, Johnson KC, et al. Estrogen and colorectal cancer incidence and mortality. Cancer. 2015;121:3261–71.
Haupt S, Caramia F, Klein SL, Rubin JB, Haupt Y. Sex disparities matter in cancer development and therapy. Nat Rev Cancer. 2021;21:393–407.
Bustos V, Nolan ÁM, Nijhuis A, Harvey H, Parker A, Poulsom R, et al. GPER mediates differential effects of estrogen on colon cancer cell proliferation and migration under normoxic and hypoxic conditions. Oncotarget. 2017;8:84258–75.
Sato R, Suzuki T, Katayose Y, Miura K, Shiiba K, Miki Y, et al. Aromatase in colon carcinoma. Anticancer Res. 2012;32:3069–75.
The authors thank the Experimental Teaching Center at the Wuhan University School of Basic Medical Sciences for technical assistance.
This work was partly supported by the National Natural Science Foundation of China under Grant [82073399 and 81872443 (to Jing Yang), 82102934 (to Yang Wang), 82203831 (to Xuewei Chen)] and Medical Science Advancement Program (Basic Medical Science) of Wuhan University under Grant [TFJC 2018003].
Ethics approval and consent to participate
The animal studies were approved by the Animal Research Committee of Wuhan University and maintained in accordance with the guidelines by the Association for Assessment and Accreditation of Laboratory Animal Care International. All human tissue and blood samples were collected with informed consent from patients, and all related procedures were performed with the approval of ethics boards of the Medical School of Wuhan University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Stratified analysis of the LMrisk based on clinicopathological features including age, gender, T stage, N stage, M stage and TNM stage in TCGA dataset. Fig. S2. The LMrisk is an independent prognostic indicator in colon cancer. Fig. S3. Establishment and validation of the prognostic nomogram for colon cancer patients. Fig. S4. High CYP19A1 expression predicts poor prognosis and positively correlated with PD-L1 expression in the GEPIA webserver. Fig. S5. CYP19A1 inhibitor letrozole facilitates anti-PD-1 therapy in mice bearing orthotopic MC38 colon tumor. Fig. S6. A proposed mechanism to explain the role of CYP19A1 in tumor immune microenvironment in colon cancer. Supplementary Materials and Methods. Supplementary Table S1. The relationships between CYP19A1 expression and clinicopathological features including age, gender, T stage, N stage M stage and TNM stage in the tissue microarray. Supplementary Table S2. Univariate and multivariate Cox regression analyses of CYP19A1 expression in the human colon cancer tissue microarray. Supplementary Table S3. Effects of letrozole on body weight, biochemical profile and complete blood counts in the orthotopic MC38 tumor model.
About this article
Cite this article
Liu, L., Mo, M., Chen, X. et al. Targeting inhibition of prognosis-related lipid metabolism genes including CYP19A1 enhances immunotherapeutic response in colon cancer. J Exp Clin Cancer Res 42, 85 (2023). https://doi.org/10.1186/s13046-023-02647-8
- Lipid metabolism
- Tumor immune microenvironment
- Immunotherapeutic response
- Colon cancer