Index of Cancer-Associated Fibroblasts Is Superior to the Epithelial-Mesenchymal Transition Score in Prognosis Prediction

The tumor microenvironment: shaping cancer progression and treatment response

The tumor microenvironment (TME) plays a crucial role in cancer progression and treatment response. It comprises a complex network of stromal cells, immune cells, extracellular matrix, and blood vessels, all of which interact with cancer cells and influence tumor behaviour. This review article provides an in-depth examination of the TME, focusing on stromal cells, blood vessels, signaling molecules, and ECM, along with commonly available therapeutic compounds that target these components. Moreover, we explore the TME as a novel strategy for discovering new anti-tumor drugs. The dynamic and adaptive nature of the TME offers opportunities for targeting specific cellular interactions and signaling pathways. We discuss emerging approaches, such as combination therapies that simultaneously target cancer cells and modulate the TME. Finally, we address the challenges and future prospects in targeting the TME. Overcoming drug resistance, improving drug delivery, and identifying new therapeutic targets within the TME are among the challenges discussed. We also highlight the potential of personalized medicine and the integration of emerging technologies, such as immunotherapy and nanotechnology, in TME-targeted therapies. This comprehensive review provides insights into the TME and its therapeutic implications. Understanding the TME's complexity and targeting its components offer promising avenues for the development of novel anti-tumor therapies and improved patient outcomes.

The Way to Malignant Transformation: Can Epigenetic Alterations Be Used to Diagnose Early-Stage Head and Neck Cancer?

Identifying and treating tumors early is the key to secondary prevention in cancer control. At present, prevention of oral cancer is still challenging because the molecular drivers responsible for malignant transformation of the 11 clinically defined oral potentially malignant disorders are still unknown. In this review, we focused on studies that elucidate the epigenetic alterations demarcating malignant and nonmalignant epigenomes and prioritized findings from clinical samples. Head and neck included, the genomes of many cancer types are largely hypomethylated and accompanied by focal hypermethylation on certain specific regions. We revisited prior studies that demonstrated that sufficient uptake of folate, the primary dietary methyl donor, is associated with oral cancer reduction. As epigenetically driven phenotypic plasticity, a newly recognized hallmark of cancer, has been linked to tumor initiation, cell fate determination, and drug resistance, we discussed prior findings that might be associated with this hallmark, including gene clusters (11q13.3, 19q13.43, 20q11.2, 22q11-13) with great potential for oral cancer biomarkers, and successful examples in screening early-stage nasopharyngeal carcinoma. Although one-size-fits-all approaches have been shown to be ineffective in most cancer therapies, the rapid development of epigenome sequencing methods raises the possibility that this nonmutagenic approach may be an exception. Only time will tell.

A Novel Gene List Identifies Tumors with a Stromal-Mesenchymal Phenotype and Worse Prognosis in Gastric Cancer

BACKGROUND

Molecular biomarkers that predict disease progression can help identify tumor subtypes and shape treatment plans. In this study, we aimed to identify robust biomarkers of prognosis in gastric cancer based on transcriptomic data obtained from primary gastric tumors.

METHODS

Microarray, RNA sequencing, and single-cell RNA sequencing-based gene expression data from gastric tumors were obtained from public databases. Freshly frozen gastric tumors (n = 42) and matched FFPE (formalin-fixed, paraffin-embedded) (n = 40) tissues from a Turkish gastric cancer cohort were used for quantitative real-time PCR and immunohistochemistry-based assessments of gene expression, respectively.

RESULTS

A novel list of 20 prognostic genes was identified and used for the classification of gastric tumors into two major tumor subgroups with differential stromal gene expression ("Stromal-UP" (SU) and "Stromal-DOWN" (SD)). The SU group had a more mesenchymal profile with an enrichment of extracellular matrix-related gene sets and a poor prognosis compared to the SD group. Expression of the genes within the signature correlated with the expression of mesenchymal markers ex vivo. A higher stromal content in FFPE tissues was associated with shorter overall survival.

CONCLUSIONS

A stroma-rich, mesenchymal subgroup among gastric tumors identifies an unfavorable clinical outcome in all cohorts tested.

Dissecting the functions of cancer-associated fibroblasts to therapeutically target head and neck cancer microenvironment

Head and neck cancers (HNC) are a diverse group of aggressive malignancies with high morbidity and mortality, leading to almost half-million deaths annually worldwide. A better understanding of the molecular processes governing tumor formation and progression is crucial to improve current diagnostic and prognostic tools as well as to develop more personalized treatment strategies. Tumors are highly complex and heterogeneous structures in which growth and dissemination is not only governed by the cancer cells intrinsic mechanisms, but also by the surrounding tumor microenvironment (TME). Cancer-associated fibroblasts (CAFs) emerge as predominant TME components and key players in the generation of permissive conditions that ultimately impact in tumor progression and metastatic dissemination. Although CAFs were initially considered a consequence of tumor development, it is now well established that they actively contribute to numerous cancer hallmarks i.e., tumor cell growth, migration and invasion, cancer cell stemness, angiogenesis, metabolic reprograming, inflammation, and immune system modulation. In this scenario, therapeutic strategies targeting CAF functions could potentially have a major impact in cancer therapeutics, providing avenues for new treatment options or for improving efficacy in established approaches. This review is focused on thoroughly dissecting existing evidences supporting the contribution of CAFs in HNC biology with an emphasis on current knowledge of the key molecules and pathways involved in CAF-tumor crosstalk, and their potential as novel biomarkers and/or therapeutic targets to effectively interfere the tumor-stroma crosstalk for HNC patients benefit. involved in CAF-tumor crosstalk, and their potential as novel biomarkers and/or therapeutic targets to effec- tively interfere the tumor-stroma crosstalk for HNC patients benefit.

Mimicking Tumor Cell Heterogeneity of Colorectal Cancer in a Patient-derived Organoid-Fibroblast Model

BACKGROUND & AIMS

Patient-derived organoid cancer models are generated from epithelial tumor cells and reflect tumor characteristics. However, they lack the complexity of the tumor microenvironment, which is a key driver of tumorigenesis and therapy response. Here, we developed a colorectal cancer organoid model that incorporates matched epithelial cells and stromal fibroblasts.

METHODS

Primary fibroblasts and tumor cells were isolated from colorectal cancer specimens. Fibroblasts were characterized for their proteome, secretome, and gene expression signatures. Fibroblast/organoid co-cultures were analyzed by immunohistochemistry and compared with their tissue of origin, as well as on gene expression levels compared with standard organoid models. Bioinformatics deconvolution was used to calculate cellular proportions of cell subsets in organoids based on single-cell RNA sequencing data.

RESULTS

Normal primary fibroblasts, isolated from tumor adjacent tissue, and cancer associated fibroblasts retained their molecular characteristics in vitro, including higher motility of cancer associated compared with normal fibroblasts. Importantly, both cancer-associated fibroblasts and normal fibroblasts supported cancer cell proliferation in 3D co-cultures, without the addition of classical niche factors. Organoids grown together with fibroblasts displayed a larger cellular heterogeneity of tumor cells compared with mono-cultures and closely resembled the in vivo tumor morphology. Additionally, we observed a mutual crosstalk between tumor cells and fibroblasts in the co-cultures. This was manifested by considerably deregulated pathways such as cell-cell communication and extracellular matrix remodeling in the organoids. Thrombospondin-1 was identified as a critical factor for fibroblast invasiveness.

CONCLUSION

We developed a physiological tumor/stroma model, which will be vital as a personalized tumor model to study disease mechanisms and therapy response in colorectal cancer.

Identification and verification of eight cancer-associated fibroblasts related genes as a prognostic signature for head and neck squamous cell carcinoma

Cancer-associated fibroblasts (CAFs) can exert their immunosuppressive effects by secreting various effectors that are involved in the regulation of tumor-infiltrating immune cells as well as other immune components in the tumor immune microenvironment (TIME), thereby promoting tumorigenesis, progression, metastasis, and drug resistance. Although a large number of studies suggest that CAFs play a key regulatory role in the development of head and neck squamous cell carcinoma (HNSCC), there are limited studies on the relevance of CAFs to the prognosis of HNSCC. In this study, we identified a prognostic signature containing eight CAF-related genes for HNSCC by univariate Cox analysis, lasso regression, stepwise regression, and multivariate Cox analysis. Our validation in primary cultures of CAFs from human HNSCC and four human HNSCC cell lines confirmed that these eight genes are indeed characteristic markers of CAFs. Immune cell infiltration differences analysis between high-risk and low-risk groups according to the eight CAF-related genes signature hinted at CAFs regulatory roles in the TIME, further revealing its potential role on prognosis. The signature of the eight CAF-related genes was validated in different independent validation cohorts and all showed that it was a valid marker for prognosis. The significantly higher overall survival (OS) in the low-risk group compared to the high-risk group was confirmed by Kaplan-Meier (K-M) analysis, suggesting that the signature of CAF-related genes can be used as a non-invasive predictive tool for HNSCC prognosis. The low-risk group had significantly higher levels of tumor-killing immune cell infiltration, as confirmed by CIBERSORT analysis, such as CD8⁺ T cells, follicular helper T cells, and Dendritic cells (DCs) in the low-risk group. In contrast, the level of infiltration of pro-tumor cells such as M0 macrophages and activated Mast cells (MCs) was lower. It is crucial to delve into the complex mechanisms between CAFs and immune cells to find potential regulatory targets and may provide new evidence for subsequently targeted immunotherapy. These results suggest that the signature of the eight CAF-related genes is a powerful indicator for the assessment of the TIME of HNSCC. It may provide a new and reliable potential indicator for clinicians to predict the prognosis of HNSCC, which may be used to guide treatment and clinical decision-making in HNSCC patients. Meanwhile, CAF-related genes are expected to become tumor biomarkers and effective targets for HNSCC.

Driving role of head and neck cancer cell secretome on the invasion of stromal fibroblasts: Mechanistic insights by phosphoproteomics

BACKGROUND

Cancer-associated fibroblasts (CAFs) are major players in tumor-stroma communication, and participate in several cancer hallmarks to drive tumor progression and metastatic dissemination. This study investigates the driving effects of tumor-secreted factors on CAF biology, with the ultimate goal of identifying effective therapeutic targets/strategies for head and neck squamous cell carcinomas (HNSCC).

METHODS

Functionally, conditioned media (CM) from different HNSCC-derived cell lines and normal keratinocytes (Kc) were tested on the growth and invasion of populations of primary CAFs and normal fibroblasts (NFs) using 3D invasion assays in collagen matrices. The changes in MMPs expression were evaluated by RT-qPCR and kinase enrichment was analyzed using mass spectrometry phosphoproteomics.

RESULTS

Our results consistently demonstrate that HNSCC-secreted factors (but not Kc CM) specifically and robustly promoted pro-invasive properties in both CAFs and NFs, thereby reflecting the plasticity of fibroblast subtypes. Concomitantly, HNSCC-secreted factors massively increased metalloproteinases levels in CAFs and NFs. By contrast, HNSCC CM and Kc CM exhibited comparable growth-promoting effects on stromal fibroblasts. Mechanistically, phosphoproteomic analysis predominantly revealed phosphorylation changes in fibroblasts upon treatment with HNSCC CM, and various promising kinases were identified: MKK7, MKK4, ASK1, RAF1, BRAF, ARAF, COT, PDK1, RSK2 and AKT1. Interestingly, pharmacologic inhibition of RAF1/BRAF using sorafenib emerged as the most effective drug to block tumor-promoted fibroblast invasion without affecting fibroblast viability CONCLUSIONS: Our findings demonstrate that HNSCC-secreted factors specifically fine tune the invasive potential of stromal fibroblasts, thereby generating tumor-driven pro-invasive niches, which in turn to ultimately facilitate cancer cell dissemination. Furthermore, the RAF/BRAF inhibitor sorafenib was identified as a promising candidate to effectively target the onset of pro-invasive clusters of stromal fibroblasts in the HNSCC microenvironment.

A cancer-associated fibroblast gene signature predicts prognosis and therapy response in patients with pancreatic cancer

Pancreatic cancer is a lethal malignancy with a 5-year survival rate of about 10% in the United States, and it is becoming an increasingly prominent cause of cancer death. Among pancreatic cancer patients, pancreatic ductal adenocarcinoma (PDAC) accounts for more than 90% of all cases and has a very poor prognosis with an average survival of only 1 year in about 18% of all tumor stages. In the past years, there has been an increasing interest in cancer-associated fibroblasts (CAFs) and their roles in PDAC. Recent data reveals that CAFs in PDAC are heterogeneous and various CAF subtypes have been demonstrated to promote tumor development while others hinder cancer proliferation. Furthermore, CAFs and other stromal populations can be potentially used as novel prognostic markers in cancer. In the present study, in order to evaluate the prognostic value of CAFs in PDAC, CAF infiltration rate was evaluated in 4 PDAC datasets of TCGA, GEO, and ArrayExpress databases and differentially expressed genes (DEGs) between CAF-high and CAF-low patients were identified. Subsequently, a CAF-based gene expression signature was developed and studied for its association with overall survival (OS). Additionally, functional enrichment analysis, somatic alteration analysis, and prognostic risk model construction was conducted on the identified DEGs. Finally, oncoPredict algorithm was implemented to assess drug sensitivity prediction between high- and low-risk cohorts. Our results revealed that CAF risk-high patients have a worse survival rate and increased CAF infiltration is a poor prognostic indicator in pancreatic cancer. Functional enrichment analysis also revealed that "extracellular matrix organization" and "vasculature development" were the top enriched pathways among the identified DEGs. We also developed a panel of 12 genes, which in additional to its prognostic value, could predict higher chemotherapy resistance rate. This CAF-based panel can be potentially utilized alone or in conjunction with other clinical parameters to make early predictions and prognosticate responsiveness to treatment in PDAC patients. Indeed, it is necessary to conduct extensive prospective investigations to confirm the clinical utility of these findings.

Interplay Between Immune and Cancer-Associated Fibroblasts: A Path to Target Metalloproteinases in Penile Cancer

Extracellular matrix (ECM) remodeling and inflammation have been reported in penile carcinomas (PeCa). However, the cell types and cellular crosstalk involved in PeCa are unexplored. We aimed to characterize the complexity of cells and pathways involved in the tumor microenvironment (TME) in PeCa and propose target molecules associated with the TME. We first investigated the prognostic impact of cell types with a secretory profile to identify drug targets that modulate TME-enriched cells. The secretome analysis using the PeCa transcriptome revealed the enrichment of inflammation and extracellular matrix pathways. Twenty-three secreted factors were upregulated, mainly collagens and matrix metalloproteinases (MMPs). The deregulation of collagens and MMPs was confirmed by Quantitative reverse transcription - polymerase chain reaction (RT-qPCR). Further, the deconvolution method (digital cytometry) of the bulk samples revealed a high proportion of macrophages and dendritic cells (DCs) and B cells. Increased DCs and B cells were associated with better survival. A high proportion of cancer-associated fibroblasts (CAFs) was observed in low-survival patients. Patients with increased CAFs had decreased immune cell proportions. The treatment with the MMP inhibitor GM6001 in CAF cells derived from PeCa resulted in altered cell viability. We reported a crosstalk between immune cells and CAFs, and the proportion of these cell populations was associated with prognosis. We demonstrate that a drug targeting MMPs modulates CAFs, expanding the therapeutic options of PeCa.

ANXA1: An Important Independent Prognostic Factor and Molecular Target in Glioma

Objective: The expression, prognosis, and related mechanisms of ANXA1 are investigated in glioma, with the objective to find potential therapeutic molecular targets for glioma. Methods: We analyzed the gene expression of ANXA1 using glioma-related databases, including the Chinese Glioma Genome Atlas (CGGA) database, The Cancer Genome Atlas (TCGA) database, and the Gene Expression Omnibus (GEO) database. Moreover, we collected the sample tissues and corresponding paracancerous tissues of 23 glioma patients and then conducted a Western blot experiment to verify the expression and correlate survival of ANXA1. Moreover, we generated survival ROC curves, performing univariate and multivariate Cox analyses and the construction of the nomogram. Differential expression analysis was conducted by high and low grouping based on the median of the ANXA1 gene expression values. We conducted Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and Gene Set Enrichment Analysis (GSEA) to explore possible mechanisms, and gene co-expression analysis was also performed. Results: The results showed that the ANXA1 expression level was higher in gliomas than in normal tissues, and a high expression level of ANXA1 in gliomas was associated with poorer prognosis. The independent prognosis analysis showed that the ANXA1 gene was an independent prognostic factor of glioma. In the analysis of KEGG and Gene Set Enrichment Analysis (GSEA), it is shown that ANXA1 may play an important role in glioma patients by affecting extracellular matrix (ECM)-receptor interaction and the focal adhesion signal pathway. The core genes, including COL1A1, COL1A2, FN1, ITGA1, and ITGB1, were screened for gene correlation and prognosis analysis. The expression level of the five genes was verified by qPCR in glioma. We concluded that these five core genes and ANXA1 could play a synergistic role in gliomas. Conclusion: The results indicated that a high expression level of ANXA1 leads to worse prognosis and ANXA1 is an independent prognostic factor and a potentially important target for the treatment of gliomas.

Identification of Prognostic Biomarkers Originating From the Tumor Stroma of Betel Quid-Associated Oral Cancer Tissues

Background

Partial epithelial-mesenchymal transition (p-EMT) is a distinct clinicopathological feature prevalent in oral cavity tumors of The Cancer Genome Atlas. Located at the invasion front, p-EMT cells require additional support from the tumor stroma for collective cell migration, including track clearing, extracellular matrix remodeling and immune evasion. The pathological roles of otherwise nonmalignant cancer-associated fibroblasts (CAFs) in cancer progression are emerging.

Methods

Gene set enrichment analysis was used to reveal differentially enriched genes and molecular pathways in OC3 and TW2.6 xenograft tissues, representing mesenchymal and p-EMT tumors, respectively. R packages of genomic data science were executed for statistical evaluations and data visualization. Immunohistochemistry and Alcian blue staining were conducted to validate the bioinformatic results. Univariate and multivariate Cox proportional hazards models were performed to identify covariates significantly associated with overall survival in clinical datasets. Kaplan–Meier curves of estimated overall survival were compared for statistical difference using the log-rank test.

Results

Compared to mesenchymal OC3 cells, tumor stroma derived from p-EMT TW2.6 cells was significantly enriched in microvessel density, tumor-excluded macrophages, inflammatory CAFs, and extracellular hyaluronan deposition. By translating these results to clinical transcriptomic datasets of oral cancer specimens, including the Puram single-cell RNA-seq cohort comprising ~6000 cells, we identified the expression of stromal TGFBI and HYAL1 as independent poor and protective biomarkers, respectively, for 40 Taiwanese oral cancer tissues that were all derived from betel quid users. In The Cancer Genome Atlas, TGFBI was a poor marker not only for head and neck cancer but also for additional six cancer types and HYAL1 was a good indicator for four tumor cohorts, suggesting common stromal effects existing in different cancer types.

Conclusions

As the tumor stroma coevolves with cancer progression, the cellular origins of molecular markers identified from conventional whole tissue mRNA-based analyses should be cautiously interpreted. By incorporating disease-matched xenograft tissue and single-cell RNA-seq results, we suggested that TGFBI and HYAL1, primarily expressed by stromal CAFs and endothelial cells, respectively, could serve as robust prognostic biomarkers for oral cancer control.

Pan-Cancer Analysis of the Associations of TGFBI Expression With Prognosis and Immune Characteristics

Transforming growth factor-beta-induced (TGFBI) protein has important roles in tumor growth, metastasis, and immunity. However, there is currently no pan-cancer evidence regarding TGFBI. In this study, we conducted a pan-cancer analysis of TGFBI mRNA and protein expression and prognoses of various cancer types using public databases. We also investigated the associations of TGFBI expression with tumor microenvironment (TME) components, immune cell infiltration, tumor mutational burden (TMB), and microsatellite instability (MSI), along with the TGFBI genetic alteration types. The results showed that TGFBI expression varied among different cancer types, and it was positively or negatively related to prognosis in various cancers. TGFBI expression was also significantly correlated with TME components, TMB, MSI, immune cell infiltration, and immunoinhibitory and immunostimulatory gene subsets. These findings indicate that TGFBI participates in various immune responses and it may function as a prognostic marker in various cancers. The findings may be useful for developing immunotherapies that target TGFBI.

A Prognostic Model for Brain Glioma Patients Based on 9 Signature Glycolytic Genes

Objective

To screen glycolytic genes linked to the glioma prognosis and construct the prognostic model.

Methods

The relevant data of glioma were downloaded from TCGA and GTEx databases. GSEA of glycolysis-related pathways was carried out, and enriched differential genes were extracted. Screening out prognostic-related genes with conspicuous significance and construction of the prognostic model were conducted by multivariate Cox regression analysis and Lasso regression analysis. The model was evaluated, and cBioPortal was used to analyze the mutation of the model gene. The expression of the model gene in tumor and normal colon tissue was analyzed. The model was used to evaluate the prognosis of patients in different groups to verify the applicability of the model.

Results

339 differentially glycolytic-related genes were enriched in REACTOME_GLYCOLYSIS, GLYCOLYTIC_PROCESS, HALLMARK_GLYCOLYSIS, and other pathways. We obtained 9 key prognostic genes and constructed the prognostic evaluation model. The 3-year AUC values of the ROC curve display model are greater than 0.75, which indicates that the accuracy of the model is good. The relation of age and risk score to prognosis is shown by univariate and multivariate Cox analysis. The expression of SRD5A3, MDH2, and B3GAT3 genes was significantly upregulated in the tumor tissues, while the HDAC4 and G6PC2 genes were downregulated. The mutation rate of MDH2 and HDAC4 genes was the highest. This model could effectively distinguish the risk of poor prognosis of patients in any age stage.

Conclusion

The prognostic assessment models based on glycolysis-related nine-gene signature could accurately predict the prognosis of patients with GBM.