A Patient-Derived, Pan-Cancer EMT Signature Identifies Global Molecular Alterations and Immune Target Enrichment Following Epithelial-to-Mesenchymal Transition

Differential role of glucocorticoid receptor based on its cell type specific expression on tumor cells and infiltrating lymphocytes

BACKGROUND

The glucocorticoid receptor (GR) is frequently expressed in breast cancer (BC), and its prognostic implications are contingent on estrogen receptor (ER) status. To address conflicting reports and explore therapeutic potential, a GR signature (GRsig) independent of ER status was developed. We also investigated cell type-specific GR protein expression in BC tumor epithelial cells and infiltrating lymphocytes.

METHODS

GRsig was derived from Dexamethasone treated cell lines through a bioinformatic pipeline. Immunohistochemistry assessed GR protein expression. Associations between GRsig and tumor phenotypes (proliferation, cytolytic activity (CYT), immune cell distribution, and epithelial-to-mesenchymal transition (EMT) were explored in public datasets. Single-cell RNA sequencing data evaluated context-dependent GR roles, and a cell type-specific prognostic role was assessed in an independent BC cohort.

RESULTS

High GRsig levels were associated with a favorable prognosis across BC subtypes. Tumor-specific high GRsig correlated with lower proliferation, increased CYT, and anti-tumorigenic immune cells. Single-cell data analysis revealed higher GRsig expression in immune cells, negatively correlating with EMT while a positive correlation was observed with EMT primarily in tumor and stromal cells. Univariate and multivariate analyses demonstrated the robust and independent predictive capability of GRsig for favorable prognosis. GR protein expression on immune cells in triple-negative tumors indicated a favorable prognosis.

CONCLUSION

This study underscores the cell type-specific role of GR, where its expression on tumor cells is associated with aggressive features like EMT, while in infiltrating lymphocytes, it predicts a better prognosis, particularly within TNBC tumors. The GRsig emerges as a promising independent prognostic indicator across diverse BC subtypes.

Metabolism pathway-based subtyping in endometrial cancer: An integrated study by multi-omics analysis and machine learning algorithms

Endometrial cancer (EC), the second most common malignancy in the female reproductive system, has garnered increasing attention for its genomic heterogeneity, but understanding of its metabolic characteristics is still poor. We explored metabolic dysfunctions in EC through a comprehensive multi-omics analysis (RNA-seq datasets from The Cancer Genome Atlas [TCGA], Cancer Cell Line Encyclopedia [CCLE], and GEO datasets; the Clinical Proteomic Tumor Analysis Consortium [CPTAC] proteomics; CCLE metabolomics) to develop useful molecular targets for precision therapy. Unsupervised consensus clustering was performed to categorize EC patients into three metabolism-pathway-based subgroups (MPSs). These MPS subgroups had distinct clinical prognoses, transcriptomic and genomic alterations, immune microenvironment landscape, and unique patterns of chemotherapy sensitivity. Moreover, the MPS2 subgroup had a better response to immunotherapy. Finally, three machine learning algorithms (LASSO, random forest, and stepwise multivariate Cox regression) were used for developing a prognostic metagene signature based on metabolic molecules. Thus, a 13-hub gene-based classifier was constructed to predict patients' MPS subtypes, offering a more accessible and practical approach. This metabolism-based classification system can enhance prognostic predictions and guide clinical strategies for immunotherapy and metabolism-targeted therapy in EC.

Integrative molecular and spatial analysis reveals evolutionary dynamics and tumor-immune interplay of in situ and invasive acral melanoma

In acral melanoma (AM), progression from in situ (AMis) to invasive AM (iAM) leads to significantly reduced survival. However, evolutionary dynamics during this process remain elusive. Here, we report integrative molecular and spatial characterization of 147 AMs using genomics, bulk and single-cell transcriptomics, and spatial transcriptomics and proteomics. Vertical invasion from AMis to iAM displays an early and monoclonal seeding pattern. The subsequent regional expansion of iAM exhibits two distinct patterns, clonal expansion and subclonal diversification. Notably, molecular subtyping reveals an aggressive iAM subset featured with subclonal diversification, increased epithelial-mesenchymal transition (EMT), and spatial enrichment of APOE+/CD163+ macrophages. In vitro and ex vivo experiments further demonstrate that APOE+CD163+ macrophages promote tumor EMT via IGF1-IGF1R interaction. Adnexal involvement can predict AMis with higher invasive potential whereas APOE and CD163 serve as prognostic biomarkers for iAM. Altogether, our results provide implications for the early detection and treatment of AM.

A Histone Methylation-MAPK Signaling Axis Drives Durable Epithelial-Mesenchymal Transition in Hypoxic Pancreatic Cancer

The tumor microenvironment in pancreatic ductal adenocarcinoma (PDAC) plays a key role in tumor progression and response to therapy. The dense PDAC stroma causes hypovascularity, which leads to hypoxia. Here, we showed that hypoxia drives long-lasting epithelial-mesenchymal transition (EMT) in PDAC primarily through a positive-feedback histone methylation-MAPK signaling axis. Transformed cells preferentially underwent EMT in hypoxic tumor regions in multiple model systems. Hypoxia drove a cell autonomous EMT in PDAC cells, which, unlike EMT in response to growth factors, could last for weeks. Furthermore, hypoxia reduced histone demethylase KDM2A activity, suppressed PP2 family phosphatase expression, and activated MAPKs to post-translationally stabilize histone methyltransferase NSD2, leading to an H3K36me2-dependent EMT in which hypoxia-inducible factors played only a supporting role. Hypoxia-driven EMT could be antagonized in vivo by combinations of MAPK inhibitors. Collectively, these results suggest that hypoxia promotes durable EMT in PDAC by inducing a histone methylation-MAPK axis that can be effectively targeted with multidrug therapies, providing a potential strategy for overcoming chemoresistance.

SIGNIFICANCE

Integrated regulation of histone methylation and MAPK signaling by the low-oxygen environment of pancreatic cancer drives long-lasting EMT that promotes chemoresistance and shortens patient survival and that can be pharmacologically inhibited. See related commentary by Wirth and Schneider, p. 1739.

Silibinin is a suppressor of the metastasis-promoting transcription factor ID3

BACKGROUND

ID3 (inhibitor of DNA binding/differentiation-3) is a transcription factor that enables metastasis by promoting stem cell-like properties in endothelial and tumor cells. The milk thistle flavonolignan silibinin is a phytochemical with anti-metastatic potential through largely unknown mechanisms.

HYPOTHESIS/PURPOSE

We have mechanistically investigated the ability of silibinin to inhibit the aberrant activation of ID3 in brain endothelium and non-small cell lung cancer (NSCLC) models.

METHODS

Bioinformatic analyses were performed to investigate the co-expression correlation between ID3 and bone morphogenic protein (BMP) ligands/BMP receptors (BMPRs) genes in NSCLC patient datasets. ID3 expression was assessed by immunoblotting and qRT-PCR. Luciferase reporter assays were used to evaluate the gene sequences targeted by silibinin to regulate ID3 transcription. In silico computational modeling and LanthaScreen TR-FRET kinase assays were used to characterize and validate the BMPR inhibitory activity of silibinin. Tumor tissues from NSCLC xenograft models treated with oral silibinin were used to evaluate the in vivo anti-ID3 effects of silibinin.

RESULTS

Analysis of lung cancer patient datasets revealed a top-ranked positive association of ID3 with the BMP9 endothelial receptor ACVRL1/ALK1 and the BMP ligand BMP6. Silibinin treatment blocked the BMP9-induced activation of the ALK1-phospho-SMAD1/5-ID3 axis in brain endothelial cells. Constitutive, acquired, and adaptive expression of ID3 in NSCLC cells were all significantly downregulated in response to silibinin. Silibinin blocked ID3 transcription via BMP-responsive elements in ID3 gene enhancers. Silibinin inhibited the kinase activities of BMPRs in the micromolar range, with the lower IC50 values occurring against ACVRL1/ALK1 and BMPR2. In an in vivo NSCLC xenograft model, tumoral overexpression of ID3 was completely suppressed by systematically achievable oral doses of silibinin.

CONCLUSIONS

ID3 is a largely undruggable metastasis-promoting transcription factor. Silibinin is a novel suppressor of ID3 that may be explored as a novel therapeutic approach to interfere with the metastatic dissemination capacity of NSCLC.

Integrated bioinformatics investigation of adenylyl cyclase family co-expression network in bladder cancer followed by preliminary validation of member 2 (ADCY2) in tumorigenesis and prognosis

Background

The adenylyl cyclase (ADCY) gene family encodes enzymes responsible for the synthesis of cyclic adenosine monophosphate (cAMP) from adenosine triphosphate (ATP), which comprises nine transmembrane isoforms (ADCYs 1-9). Although ADCYs correlate with intracellular signalling and tumorigenesis in different malignancies, their roles in bladder cancer remain unclear.

Methods

Utilizing the bladder urothelial carcinoma (BLCA) dataset from The Cancer Genome Atlas (TCGA), we employed the R package 'limma' to identify differential genes. Subsequent correlation analysis with corresponding clinical data was conducted. Prognostic significance of ADCY family genes was assessed through survival analysis. Univariate and multivariate Cox regression determined ADCY2 as a potential independent risk factor for BLCA. Validation was performed using immunohistochemistry results from independent cohorts. Additionally, we delved into the mechanism of genetic variations, methylation modifications, and signalling pathways of ADCY family genes. Evaluation of their role in the immune microenvironment was achieved through R packages single-sample gene set enrichment analysis (ssGSEA), CIBERPORT, and ESTIMATE.

Results

Cases of bladder cancer were retrieved from TCGA, and the transcriptionally differentially expressed members of ADCY were identified (members 2, 4, and 5). Genomic alteration, epigenomic modification, clinicopathological characteristics and clinical survival were systematically investigated. A co-expression network was established based on the intersection of correlated genes, which was centred around ADCY2, ADCY4, and ADCY5. Enrichment analysis revealed that correlated genes were involved in epithelial-mesenchymal transition (EMT). The ADCY2 was selected as the most representative biomarker for prognosis in bladder cancer. Bladder tumour with higher ADCY2 expression had higher prognostic risk and worse survival outcomes. Moreover, ADCY2 was correlated with classic immune checkpoints, and a better responsiveness to immunotherapy was exhibited in high-expression subsets. To ameliorate universality of the conclusion, our study also included several real-world cohorts into the preliminary validation, using datasets from the Gene Expression Omnibus (GEO; GSE13507), tissue microarray (TMA) with 80 bladder cancer inclusion and clinical trial IMvigor210, which were associated with immunotherapy sensitivity, prognosis, and common biomarker presentation.

Conclusions

Our study reveals that ADCY family has prognostic value in patients with bladder cancer; the ADCY2 is a prominent prognostic biomarker. The bioinformatics analyses and validation provide direction for further functional and mechanistic studies on the screened members of ADCY family.

A novel subtype based on driver methylation-transcription in lung adenocarcinoma

AIMS

To identify driver methylation genes and a novel subtype of lung adenocarcinoma (LUAD) by multi-omics and elucidate its molecular features and clinical significance.

METHODS

We collected LUAD patients from public databases, and identified driver methylation genes (DMGs) by MethSig and MethylMix algrothms. And novel driver methylation multi-omics subtypes were identified by similarity network fusion (SNF). Furthermore, the prognosis, tumor microenvironment (TME), molecular features and therapy efficiency among subtypes were comprehensively evaluated.

RESULTS

147 overlapped driver methylation were identified and validated. By integrating the mRNA expression and methylation of DMGs using SNF, four distinct patterns, termed as S1-S4, were characterized by differences in prognosis, biological features, and TME. The S2 subtype showed unfavorable prognosis. By comparing the characteristics of the DMGs subtypes with the traditional subtypes, S3 was concentrated in proximal-inflammatory (PI) subtype, and S4 was consisted of terminal respiratory unit (TRU) subtype and PI subtype. By analyzing TME and epithelial mesenchymal transition (EMT) features, increased immune infiltration and higher expression of immune checkpoint genes were found in S3 and S4. While S4 showed higher EMT score and expression of EMT associated genes, indicating S4 may not be as immunosensitive as the S3. Additionally, S3 had lower TIDE and higher IPS score, indicating its increased sensitivity to immunotherapy.

CONCLUSION

The driver methylation-related subtypes of LUAD demonstrate prognostic predictive ability that could help inform treatment response and provide complementary information to the existing subtypes.

EVA1A reverses lenvatinib resistance in hepatocellular carcinoma through regulating PI3K/AKT/p53 signaling axis

Lenvatinib is a commonly used first-line drug for the treatment of advanced hepatocellular carcinoma (HCC). However, its clinical efficacy is limited due to the drug resistance. EVA1A was a newly identified tumor suppressor, nevertheless, the impact of EVA1A on resistance to lenvatinib treatment in HCC and the potential molecular mechanisms remain unknown. In this study, the expression of EVA1A in HCC lenvatinib-resistant cells is decreased and its low expression was associated with a poor prognosis of HCC. Overexpression of EVA1A reversed lenvatinib resistance in vitro and in vivo, as demonstrated by its ability to promote cell apoptosis and inhibit cell proliferation, invasion, migration, EMT, and tumor growth. Silencing EVA1A in lenvatinib-sensitive parental HCC cells exerted the opposite effect and induced resistance to lenvatinib. Mechanistically, upregulated EVA1A inhibited the PI3K/AKT/MDM2 signaling pathway, resulting in a reduced interaction between MDM2 and p53, thereby stabilizing p53 and enhancing its antitumor activity. In addition, upregulated EVA1A suppressed the PI3K/AKT/mTOR signaling pathway and promoted autophagy, leading to the degradation of mutant p53 and attenuating its oncogenic impact. On the contrary, loss of EVA1A activated the PI3K/AKT/MDM2 signaling pathway and inhibited autophagy, promoting p53 proteasomal degradation and mutant p53 accumulation respectively. These findings establish a crucial role of EVA1A loss in driving lenvatinib resistance involving a mechanism of modulating PI3K/AKT/p53 signaling axis and suggest that upregulating EVA1A is a promising therapeutic strategy for alleviating resistance to lenvatinib, thereby improving the efficacy of HCC treatment.

Epithelial-mesenchymal transition induced by tumor cell-intrinsic PD-L1 signaling predicts a poor response to immune checkpoint inhibitors in PD-L1-high lung cancer

BACKGROUND

We investigated the role of tumor cell-intrinsic PD-L1 signaling in the epithelial-mesenchymal transition (EMT) in non-small-cell lung cancer (NSCLC) and the role of EMT as a predictive biomarker for immune checkpoint inhibitor (ICI) therapy.

METHODS

PD-L1-overexpressing or PD-L1-knockdown NSCLC cells underwent RNA-seq and EMT phenotype assessment. Mouse lung cancer LLC cells were injected into nude mice. Two cohorts of patients with NSCLC undergoing ICI therapy were analyzed.

RESULTS

RNA-seq showed that EMT pathways were enriched in PD-L1-high NSCLC cells. EMT was enhanced by PD-L1 in NSCLC cells, which was mediated by transforming growth factor-β (TGFβ). PD-L1 promoted the activation of p38-MAPK by binding to and inhibiting the protein phosphatase PPM1B, thereby increasing the TGFβ production. Tumor growth and metastasis increased in nude mice injected with PD-L1-overexpressing LLC cells. In the ICI cohort, EMT signature was higher in patients with progressive disease than in those with responses, and EMT was significantly associated with poor survival in PD-L1-high NSCLC. In PD-L1-high NSCLC, EMT was associated with increased M2-macrophage and regulatory T-cell infiltrations and decreased cytotoxic T-cell infiltration.

CONCLUSIONS

Tumor cell-intrinsic PD-L1 function contributes to NSCLC progression by promoting EMT. EMT may predict an unfavorable outcome after ICI therapy in PD-L1-high NSCLC.

P. S, C. E. A, Rajarajan S, Sahoo S, Mohan G, Jaikumar VS, Ramesh RS, B. S. S, Jolly MK, Maliekal TT, Prabhu JS. Elucidating the Role of MicroRNA-18a in Propelling a Hybrid Epithelial-Mesenchymal Phenotype and Driving Malignant Progression in ER-Negative Breast Cancer

Epigenetic alterations that lead to differential expression of microRNAs (miRNAs/miR) are known to regulate tumour cell states, epithelial-mesenchymal transition (EMT) and the progression to metastasis in breast cancer. This study explores the key contribution of miRNA-18a in mediating a hybrid E/M cell state that is pivotal to the malignant transformation and tumour progression in the aggressive ER-negative subtype of breast cancer. The expression status and associated effects of miR-18a were evaluated in patient-derived breast tumour samples in combination with gene expression data from public datasets, and further validated in in vitro and in vivo breast cancer model systems. The clinical relevance of the study findings was corroborated against human breast tumour specimens (n = 446 patients). The down-regulated expression of miR-18a observed in ER-negative tumours was found to drive the enrichment of hybrid epithelial/mesenchymal (E/M) cells with luminal attributes, enhanced traits of migration, stemness, drug-resistance and immunosuppression. Further analysis of the miR-18a targets highlighted possible hypoxia-inducible factor 1-alpha (HIF-1α)-mediated signalling in these tumours. This is a foremost report that validates the dual role of miR-18a in breast cancer that is subtype-specific based on hormone receptor expression. The study also features a novel association of low miR-18a levels and subsequent enrichment of hybrid E/M cells, increased migration and stemness in a subgroup of ER-negative tumours that may be attributed to HIF-1α mediated signalling. The results highlight the possibility of stratifying the ER-negative disease into clinically relevant groups by analysing miRNA signatures.

KRAS Allelic Variants in Biliary Tract Cancers

Importance

Biliary tract cancers (BTCs) contain several actionable molecular alterations, including FGFR2, IDH1, ERBB2 (formerly HER2), and KRAS. KRAS allelic variants are found in 20% to 30% of BTCs, and multiple KRAS inhibitors are currently under clinical investigation.

Objectives

To describe the genomic landscape, co-sequence variations, immunophenotype, genomic ancestry, and survival outcomes of KRAS-mutated BTCs and to calculate the median overall survival (mOS) for the most common allelic variants.

Design, Setting, and Participants

This retrospective, multicenter, pooled cohort study obtained clinical and next-generation sequencing data from multiple databases between January 1, 2017, and December 31, 2022. These databases included Princess Margaret Cancer Centre, MD Anderson Cancer Center, Foundation Medicine, American Association for Cancer Research Project GENIE, and cBioPortal for Cancer Genomics. The cohort comprised patients with BTCs who underwent genomic testing.

Main Outcome and Measure

The main outcome was mOS, defined as date of diagnosis to date of death, which was measured in months.

Results

A total of 7457 patients (n = 3773 males [50.6%]; mean [SD] age, 63 [5] years) with BTCs and genomic testing were included. Of these patients, 5813 had clinical outcome data available, in whom 1000 KRAS-mutated BTCs were identified. KRAS allelic variants were highly prevalent in perihilar cholangiocarcinoma (28.6%) and extrahepatic cholangiocarcinoma (36.1%). Thirty-six KRAS allelic variants were identified, and the prevalence rates in descending order were G12D (41%), G12V (23%), and Q61H (8%). The variant G12D had the highest mOS of 25.1 (95% CI, 22.0-33.0) months compared with 22.8 (95% CI, 19.6-31.4) months for Q61H and 17.8 (95% CI, 16.3-23.1) months for G12V variants. The majority of KRAS-mutated BTCs (98.9%) were not microsatellite instability-high and had low tumor mutational burden (ranging from a median [IQR] of 1.2 (1.2-2.5) to a mean [SD] of 3.3 [1.3]). Immune profiling through RNA sequencing of KRAS and NRAS-mutated samples showed a pattern toward a more immune-inflamed microenvironment with higher M1 macrophage activation (0.16 vs 0.12; P = .047) and interferon-γ expression compared with wild-type tumors. The G12D variant remained the most common KRAS allelic variant in all patient ancestries. Patients with admixed American ancestry had the highest proportion of G12D variant (45.0%).

Conclusions and Relevance

This cohort study found that KRAS allelic variants were relatively common and may be potentially actionable genomic alterations in patients with BTCs, especially perihilar cholangiocarcinoma and extrahepatic cholangiocarcinoma. The findings add to the growing data on genomic and immune landscapes of KRAS allelic variants in BTCs and are potentially of value to the planning of specific therapies for this heterogeneous patient group.

Immune checkpoints between epithelial-mesenchymal transition and autophagy: A conflicting triangle

Inhibitory immune checkpoint (ICP) molecules are pivotal in inhibiting innate and acquired antitumor immune responses, a mechanism frequently exploited by cancer cells to evade host immunity. These evasion strategies contribute to the complexity of cancer progression and therapeutic resistance. For this reason, ICP molecules have become targets for antitumor drugs, particularly monoclonal antibodies, collectively referred to as immune checkpoint inhibitors (ICI), that counteract such cancer-associated immune suppression and restore antitumor immune responses. Over the last decade, however, it has become clear that tumor cell-associated ICPs can also induce tumor cell-intrinsic effects, in particular epithelial-mesenchymal transition (EMT) and macroautophagy (hereafter autophagy). Both of these processes have profound implications for cancer metastasis and drug responsiveness. This article reviews the positive or negative cross-talk that tumor cell-associated ICPs undergo with autophagy and EMT. We discuss that tumor cell-associated ICPs are upregulated in response to the same stimuli that induce EMT. Moreover, ICPs themselves, when overexpressed, become an EMT-inducing stimulus. As regards the cross-talk with autophagy, ICPs have been shown to either stimulate or inhibit autophagy, while autophagy itself can either up- or downregulate the expression of ICPs. This dynamic equilibrium also extends to the autophagy-apoptosis axis, further emphasizing the complexities of cellular responses. Eventually, we delve into the intricate balance between autophagy and apoptosis, elucidating its role in the broader interplay of cellular dynamics influenced by ICPs. In the final part of this article, we speculate about the driving forces underlying the contradictory outcomes of the reciprocal, inhibitory, or stimulatory effects between ICPs, EMT, and autophagy. A conclusive identification of these driving forces may allow to achieve improved antitumor effects when using combinations of ICIs and compounds acting on EMT and/or autophagy. Prospectively, this may translate into increased and/or broadened therapeutic efficacy compared to what is currently achieved with ICI-based clinical protocols.

Bioinformatic Analysis of m6A Regulator-Mediated RNA Methylation Modification Patterns and Immune Microenvironment Characterization in Endometriosis

Epigenetic regulation plays an essential role in immunity and inflammation in endometriosis. In this study, we aimed to explore differences in m6A regulators between endometriosis patients and normal women and analyze the effect of m6A modification on immune and inflammatory microenvironment. The samples for analysis were downloaded from the Gene Expression Omnibus database, including ectopic endometrium (EC), eutopic endometrium (EU), and normal eutopic endometrium (NM) samples from non-endometriosis women. The validation process involved utilizing our previous RNA-sequencing data. Subsequently, a correlation analysis was performed to ascertain the relationship between m6A and the inflammatory microenvironment profile, encompassing infiltrating immunocytes, immune-inflammation reaction gene sets, and human leukocyte antigen genes. LASSO analyses were used to develop risk signature. The findings of this study indicate that the m6A regulators FTO were observed to be significantly up-regulated, while YTHDF2, CBLL1, and METTL3 were down-regulated in endometriosis tissues. The CIBERSORT analysis revealed that the local inflammatory microenvironment of ectopic lesions plays a crucial role in the development of endometriosis. Notably, M2 macrophages exhibited a significant difference between the EC and NM groups. Moreover, M2 macrophages demonstrated a positive correlation with FTO (0.39) and a negative correlation with CBLL1 (- 0.35). Furthermore, consistent clustering of EC and EU samples resulted in the identification of three distinct cell subtypes. Among different cell subtypes, significant differences were in immunoinfiltrating cells, plasma cells, naive CD4 T cells, memory activated CD4 T cells, gamma delta T cells, resting NK cells and activated NK cells but not in macrophages. Furthermore, the identification of various compounds capable of targeting these m6A genes was achieved. In conclusions, our integrated bioinformatics analysis results demonstrated that m6A-related genes METTL3, CBLL1 and YTHDF2 may be useful biomarkers for endometriosis in ectopic endometrium. The potential therapeutic approach of targeting m6A regulators holds promise for the treatment of endometriosis.

Epithelial-to-mesenchymal transition in cancer progression: unraveling the immunosuppressive module driving therapy resistance

Cancer cells undergo phenotypic switching (cancer cell plasticity) in response to microenvironmental cues, including exposure to therapy/treatment. Phenotypic plasticity enables the cancer cells to acquire more mesenchymal traits promoting cancer cells' growth, survival, therapy resistance, and disease recurrence. A significant program in cancer cell plasticity is epithelial-to-mesenchymal transition (EMT), wherein a comprehensive reprogramming of gene expression occurs to facilitate the translational shift from epithelial-to-mesenchymal phenotypes resulting in increased invasiveness and metastasis. In addition, EMT plays a pivotal role in facilitating cancer cells' escape from the body's immune system using several mechanisms, such as the downregulation of major histocompatibility complex-mediated antigen presentation, upregulation of immune checkpoint molecules, and recruitment of immune-suppressive cells. Cancer cells' ability to undergo phenotypic switching and EMT-driven immune escape presents a formidable obstacle in cancer management, highlighting the need to unravel the intricate mechanisms underlying these processes and develop novel therapeutic strategies. This article discusses the role of EMT in promoting immune evasion and therapy resistance. We also discuss the ongoing research on developing therapeutic approaches targeting intrinsic and induced cell plasticity within the immune suppressive microenvironment. We believe this review article will update the current research status and equip researchers, clinicians, and other healthcare professionals with valuable insights enhancing their existing knowledge and shedding light on promising directions for future cancer research. This will facilitate the development of innovative strategies for managing therapy-resistant cancers and improving patient outcomes.

Comprehensive Immunogenomic Profiling of IDH1-/2-Altered Cholangiocarcinoma

PURPOSE

Isocitrate dehydrogenase (IDH)1/2 genomic alterations (GA) occur in 20% of intrahepatic cholangiocarcinoma (iCCA); however, the immunogenomic landscape of IDH1-/2-mutated iCCA is largely unknown.

METHODS

Comprehensive genomic profiling (CGP) was performed on 3,067 cases of advanced iCCA. Tumor mutational burden (TMB), PD-L1 expression (Dako 22C3), microsatellite instability (MSI), and genomic loss of heterozygosity (gLOH) as a surrogate marker for homologous recombination deficiency were examined. RNA sequencing of 73 patient samples was analyzed for differences in stromal/immune cell infiltration, immune marker expression, and T-cell inflammation. Tissue microarray arrays were subjected to multiplex immunohistochemistry and colocalization analysis in 100 surgical samples. Retrospective clinical data were collected for 501 patients with cholangiocarcinoma to examine median overall survival (mOS) in IDH1/2+ versus IDHwt.

RESULTS

Of 3,067 iCCA cases subjected to CGP, 426 (14%) were IDH1+ and 125 (4%) were IDH2+. IDH1 GA included R132C (69%) and R132L/G/S/H/F (16%/7%/4%/3%/<1%). IDH2 GA occurred at R172 (94.4%) and R140 (6.6%). No significant difference was seen in median gLOH between IDH1+ versus IDHwt iCCA (P = .37), although patterns of comutations differed. MSI-High (P = .009), TMB ≥10 mut/Mb (P < .0001), and PD-L1 positivity were lower in IDH1/2+ versus IDHwt iCCA. Resting natural killer cell population, CD70, and programmed cell death 1 expression were significantly higher in non-IDH1-mutated cases, whereas V-set domain containing T-cell activation inhibitor 1 (B7-H4) expression was significantly higher in IDH1+. No significant difference in mOS was observed between IDH1/2+ versus IDHwt patients.

CONCLUSION

Significant differences in GA and immune biomarkers are noted between IDH1/2+ and IDHwt iCCA. IDH1-/2-mutated tumors appear immunologically cold without gLOH. These immunogenomic data provide insight for precision targeting of iCCA with IDH alterations.

Revealing neuropilin expression patterns in pancreatic cancer: From single‑cell to therapeutic opportunities (Review)

Pancreatic cancer, one of the most fatal types of human cancers, includes several non-epithelial and stromal components, such as activated fibroblasts, vascular cells, neural cells and immune cells, that are involved in different cancers. Vascular endothelial cell growth factor 165 receptors 1 [neuropilin-1 (NRP-1)] and 2 (NRP-2) play a role in the biological behaviors of pancreatic cancer and may appear as potential therapeutic targets. The NRP family of proteins serve as co-receptors for vascular endothelial growth factor, transforming growth factor β, hepatocyte growth factor, fibroblast growth factor, semaphorin 3, epidermal growth factor, insulin-like growth factor and platelet-derived growth factor. Investigations of mechanisms that involve the NRP family of proteins may help develop novel approaches for overcoming therapy resistance in pancreatic cancer. The present review aimed to provide an in-depth exploration of the multifaceted roles of the NRP family of proteins in pancreatic cancer, including recent findings from single-cell analysis conducted within the context of pancreatic adenocarcinoma, which revealed the intricate involvement of NRP proteins at the cellular level. Through these efforts, the present study endeavored to further reveal their relationships with different biological processes and their potential as therapeutic targets in various treatment modalities, offering novel perspectives and directions for the treatment of pancreatic cancer.

Stromal Signals Dominate Gene Expression Signature Scores That Aim to Describe Cancer Cell-intrinsic Stemness or Mesenchymality Characteristics

Epithelial-to-mesenchymal transition (EMT) in cancer cells confers migratory abilities, a crucial aspect in the metastasis of tumors that frequently leads to death. In multiple studies, authors proposed gene expression signatures for EMT, stemness, or mesenchymality of tumors based on bulk tumor expression profiling. However, recent studies suggested that noncancerous cells from the microenvironment or macroenvironment heavily influence such signature profiles. Here, we strengthen these findings by investigating 11 published and frequently referenced gene expression signatures that were proposed to describe EMT-related (EMT, mesenchymal, or stemness) characteristics in various cancer types. By analyses of bulk, single-cell, and pseudobulk expression data, we show that the cell type composition of a tumor sample frequently dominates scores of these EMT-related signatures. A comprehensive, integrated analysis of bulk RNA sequencing (RNA-seq) and single-cell RNA-seq data shows that stromal cells, most often fibroblasts, are the main drivers of EMT-related signature scores. We call attention to the risk of false conclusions about tumor properties when interpreting EMT-related signatures, especially in a clinical setting: high patient scores of EMT-related signatures or calls of "stemness subtypes" often result from low cancer cell content in tumor biopsies rather than cancer cell-specific stemness or mesenchymal/EMT characteristics.

SIGNIFICANCE

Cancer self-renewal and migratory abilities are often characterized via gene module expression profiles, also called EMT or stemness gene expression signatures. Using published clinical tumor samples, cancer cell lines, and single cancer cells, we highlight the dominating influence of noncancer cells in low cancer cell content biopsies on their scores. We caution on their application for low cancer cell content clinical cancer samples with the intent to assign such characteristics or subtypes.

A comprehensive clinically informed map of dependencies in cancer cells and framework for target prioritization

Genetic screens in cancer cell lines inform gene function and drug discovery. More comprehensive screen datasets with multi-omics data are needed to enhance opportunities to functionally map genetic vulnerabilities. Here, we construct a second-generation map of cancer dependencies by annotating 930 cancer cell lines with multi-omic data and analyze relationships between molecular markers and cancer dependencies derived from CRISPR-Cas9 screens. We identify dependency-associated gene expression markers beyond driver genes, and observe many gene addiction relationships driven by gain of function rather than synthetic lethal effects. By combining clinically informed dependency-marker associations with protein-protein interaction networks, we identify 370 anti-cancer priority targets for 27 cancer types, many of which have network-based evidence of a functional link with a marker in a cancer type. Mapping these targets to sequenced tumor cohorts identifies tractable targets in different cancer types. This target prioritization map enhances understanding of gene dependencies and identifies candidate anti-cancer targets for drug development.

Molecular profile of bladder cancer progression to clinically aggressive subtypes

Bladder cancer is a histologically and clinically heterogenous disease. Most bladder cancers are urothelial carcinomas, which frequently develop distinct histological subtypes. Several urothelial carcinoma histological subtypes, such as micropapillary, plasmacytoid, small-cell carcinoma and sarcomatoid, show highly aggressive behaviour and pose unique challenges in diagnosis and treatment. Comprehensive genomic characterizations of the urothelial carcinoma subtypes have revealed that they probably arise from a precursor subset of conventional urothelial carcinomas that belong to different molecular subtypes - micropapillary and plasmacytoid subtypes develop along the luminal pathway, whereas small-cell and sarcomatoid subtypes evolve along the basal pathway. The subtypes exhibit distinct genomic alterations, but in most cases their biological properties seem to be primarily determined by specific gene expression profiles, including epithelial-mesenchymal transition, urothelial-to-neural lineage plasticity, and immune infiltration with distinct upregulation of immune regulatory genes. These breakthrough studies have transformed our view of bladder cancer histological subtype biology, generated new hypotheses for therapy and chemoresistance, and facilitated the discovery of new therapeutic targets.

How important is EMT for cancer metastasis?

Epithelial-to-mesenchymal transition (EMT), a biological phenomenon of cellular plasticity initially reported in embryonic development, has been increasingly recognized for its importance in cancer progression and metastasis. Despite tremendous progress being made in the past 2 decades in our understanding of the molecular mechanism and functional importance of EMT in cancer, there are several mysteries around EMT that remain unresolved. In this Unsolved Mystery, we focus on the variety of EMT types in metastasis, cooperative and collective EMT behaviors, spatiotemporal characterization of EMT, and strategies of therapeutically targeting EMT. We also highlight new technical advances that will facilitate the efforts to elucidate the unsolved mysteries of EMT in metastasis.

Immunosurveillance encounters cancer metabolism

Tumor cells reprogram nutrient acquisition and metabolic pathways to meet their energetic, biosynthetic, and redox demands. Similarly, metabolic processes in immune cells support host immunity against cancer and determine differentiation and fate of leukocytes. Thus, metabolic deregulation and imbalance in immune cells within the tumor microenvironment have been reported to drive immune evasion and to compromise therapeutic outcomes. Interestingly, emerging evidence indicates that anti-tumor immunity could modulate tumor heterogeneity, aggressiveness, and metabolic reprogramming, suggesting that immunosurveillance can instruct cancer progression in multiple dimensions. This review summarizes our current understanding of how metabolic crosstalk within tumors affects immunogenicity of tumor cells and promotes cancer progression. Furthermore, we explain how defects in the metabolic cascade can contribute to developing dysfunctional immune responses against cancers and discuss the contribution of immunosurveillance to these defects as a feedback mechanism. Finally, we highlight ongoing clinical trials and new therapeutic strategies targeting cellular metabolism in cancer.

Integration of Pan-Cancer Single-Cell and Spatial Transcriptomics Reveals Stromal Cell Features and Therapeutic Targets in Tumor Microenvironment

Stromal cells are physiologically essential components of the tumor microenvironment (TME) that mediates tumor development and therapeutic resistance. Development of a logical and unified system for stromal cell type identification and characterization of corresponding functional properties could help design antitumor strategies that target stromal cells. Here, we performed a pan-cancer analysis of 214,972 nonimmune stromal cells using single-cell RNA sequencing from 258 patients across 16 cancer types and analyzed spatial transcriptomics from 16 patients across seven cancer types, including six patients receiving anti-PD-1 treatment. This analysis uncovered distinct features of 39 stromal subsets across cancer types, including various functional modules, spatial locations, and clinical and therapeutic relevance. Tumor-associated PGF+ endothelial tip cells with elevated epithelial-mesenchymal transition features were enriched in immune-depleted TME and associated with poor prognosis. Fibrogenic and vascular pericytes (PC) derived from FABP4+ progenitors were two distinct tumor-associated PC subpopulations that strongly interacted with PGF+ tips, resulting in excess extracellular matrix (ECM) abundance and dysfunctional vasculature. Importantly, ECM-related cancer-associated fibroblasts enriched at the tumor boundary acted as a barrier to exclude immune cells, interacted with malignant cells to promote tumor progression, and regulated exhausted CD8+ T cells via immune checkpoint ligand-receptors (e.g., LGALS9/TIM-3) to promote immune escape. In addition, an interactive web-based tool (http://www.scpanstroma.yelab.site/) was developed for accessing, visualizing, and analyzing stromal data. Taken together, this study provides a systematic view of the highly heterogeneous stromal populations across cancer types and suggests future avenues for designing therapies to overcome the tumor-promoting functions of stromal cells.

SIGNIFICANCE

Comprehensive characterization of tumor-associated nonimmune stromal cells provides a robust resource for dissecting tumor microenvironment complexity and guiding stroma-targeted therapy development across multiple human cancer types.

TERT Promoter Mutations Frequency Across Race, Sex, and Cancer Type

BACKGROUND

Telomerase reverse transcriptase (TERT) gene promoter mutations have been explored, as biomarkers of improved survival for patients with cancer receiving immune checkpoint inhibitors. We sought to investigate their prevalence by race and sex across different cancer types to inform patient selection in clinical trials.

RESULTS

In this observational study, 31 925 patients with cancer underwent next-generation sequencing of their tumors with 88% (27 970) patients self-reported being Whites, 7.1% (2273) Asians, and 5.3% (1682) Blacks. Examining the distribution of TERT promoter mutations by race, White patients with melanoma harbored more TERT promoter mutations than Asian and Black patients (OR = 25.83; 95%CI, 6.84-217.42; P < .001). In contrast, Asian patients with head and neck cancer (HNC) harbored more TERT promoter mutations compared to White patients (OR = 2.47; 95%CI, 1.39-4.37; P = .004). In addition, the distribution of TERT promoter mutations differed by sex. Males were enriched for TERT gene promoter mutations compared to females with melanoma (OR = 1.82; 95%CI, 1.53-2.16; P < .001), cancer of unknown primary (OR = 1.96; 95%CI, 1.43-2.69; P < .001), hepatobiliary (OR = 3.89; 95%CI, 2.65-5.69; P < .001), and thyroid cancers (OR = 1.42; 95%CI, 1.10-1.84; P = .0087), while females were more enriched for TERT promoter mutations compared to males for HNC (OR = 0.56; 95%CI, 0.39-0.81; P = .0021).

CONCLUSIONS

The prevalence of TERT gene promoter mutations varies among patients with cancer based on race and sex. These findings inform our understanding of cancer biology and can assist in the design of future clinical trials that leverage drugs targeting TERT promoter dependencies.

Mutant p53 gains oncogenic functions through a chromosomal instability-induced cytosolic DNA response

Inactivating TP53 mutations leads to a loss of function of p53, but can also often result in oncogenic gain-of-function (GOF) of mutant p53 (mutp53) proteins which promotes tumor development and progression. The GOF activities of TP53 mutations are well documented, but the mechanisms involved remain poorly understood. Here, we study the mutp53 interactome and find that by targeting minichromosome maintenance complex components (MCMs), GOF mutp53 predisposes cells to replication stress and chromosomal instability (CIN), leading to a tumor cell-autonomous and cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING)-dependent cytosolic DNA response that activates downstream non-canonical nuclear factor kappa light chain enhancer of activated B cell (NC-NF-κB) signaling. Consequently, GOF mutp53-MCMs-CIN-cytosolic DNA-cGAS-STING-NC-NF-κB signaling promotes tumor cell metastasis and an immunosuppressive tumor microenvironment through antagonizing interferon signaling and regulating genes associated with pro-tumorigenic inflammation. Our findings have important implications for understanding not only the GOF activities of TP53 mutations but also the genome-guardian role of p53 and its inactivation during tumor development and progression.

Dual role of CASP8AP2/FLASH in regulating epithelial-to-mesenchymal transition plasticity (EMP)

BACKGROUND

Epithelial-to-mesenchymal transition (EMT) is a developmental program that consists of the loss of epithelial features concomitant with the acquisition of mesenchymal features. Activation of EMT in cancer facilitates the acquisition of aggressive traits and cancer invasion. EMT plasticity (EMP), the dynamic transition between multiple hybrid states in which cancer cells display both epithelial and mesenchymal markers, confers survival advantages for cancer cells in constantly changing environments during metastasis.

METHODS

RNAseq analysis was performed to assess genome-wide transcriptional changes in cancer cells depleted for histone regulators FLASH, NPAT, and SLBP. Quantitative PCR and Western blot were used for the detection of mRNA and protein levels. Computational analysis was performed on distinct sets of genes to determine the epithelial and mesenchymal score in cancer cells and to correlate FLASH expression with EMT markers in the CCLE collection.

RESULTS

We demonstrate that loss of FLASH in cancer cells gives rise to a hybrid E/M phenotype with high epithelial scores even in the presence of TGFβ, as determined by computational methods using expression of predetermined sets of epithelial and mesenchymal genes. Multiple genes involved in cell-cell junction formation are similarly specifically upregulated in FLASH-depleted cells, suggesting that FLASH acts as a repressor of the epithelial phenotype. Further, FLASH expression in cancer lines is inversely correlated with the epithelial score. Nonetheless, subsets of mesenchymal markers were distinctly up-regulated in FLASH, NPAT, or SLBP-depleted cells.

CONCLUSIONS

The ZEB1low/SNAILhigh/E-cadherinhigh phenotype described in FLASH-depleted cancer cells is driving a hybrid E/M phenotype in which epithelial and mesenchymal markers coexist.

Simultaneously targeting extracellular vesicle trafficking and TGF-β receptor kinase activity blocks signaling hyperactivation and metastasis

Metastasis is the leading cause of cancer-related deaths. Transforming growth factor beta (TGF-β) signaling drives metastasis and is strongly enhanced during cancer progression. Yet, the use of on-target TGF-β signaling inhibitors in the treatment of cancer patients remains unsuccessful, highlighting a gap in the understanding of TGF-β biology that limits the establishment of efficient anti-metastatic therapies. Here, we show that TGF-β signaling hyperactivation in breast cancer cells is required for metastasis and relies on increased small extracellular vesicle (sEV) secretion. Demonstrating sEV's unique role, TGF-β signaling levels induced by sEVs exceed the activity of matching concentrations of soluble ligand TGF-β. Further, genetic disruption of sEV secretion in highly-metastatic breast cancer cells impairs cancer cell aggressiveness by reducing TGF-β signaling to nearly-normal levels. Otherwise, TGF-β signaling activity in non-invasive breast cancer cells is inherently low, but can be amplified by sEVs, enabling invasion and metastasis of poorly-metastatic breast cancer cells. Underscoring the translational potential of inhibiting sEV trafficking in advanced breast cancers, treatment with dimethyl amiloride (DMA) decreases sEV secretion, TGF-β signaling activity, and breast cancer progression in vivo. Targeting both the sEV trafficking and TGF-β signaling by combining DMA and SB431542 at suboptimal doses potentiated this effect, normalizing the TGF-β signaling in primary tumors to potently reduce circulating tumor cells, metastasis, and tumor self-seeding. Collectively, this study establishes sEVs as critical elements in TGF-β biology, demonstrating the feasibility of inhibiting sEV trafficking as a new therapeutic approach to impair metastasis by normalizing TGF-β signaling levels in breast cancer cells.

Tumor microenvironment remodeling plus immunotherapy could be used in mesenchymal-like tumor with high tumor residual and drug resistant rate

Epithelial-mesenchymal transition (EMT) is a common process during tumor progression and is always related to residual tumor, drug resistance and immune suppression. However, considering the heterogeneity in EMT process, there is still a need to establish robust EMT classification system with reasonable molecular, biological and clinical implications to investigate whether these unfavorable survival factors are common or unique in different individuals. In our work, we classify tumors with four EMT status, that is, EMTlow, EMTmid, EMThigh-NOS (Not Otherwise Specified), and EMThigh-AKT (AKT pathway overactivation) subtypes. We find that EMThigh-NOS subtype is driven by intrinsic somatic alterations. While, EMThigh-AKT subtype is maintained by extrinsic cellular interplay between tumor cells and macrophages in an AKT-dependent manner. EMThigh-AKT subtype is both unresectable and drug resistant while EMThigh-NOS subtype can be treated with cell cycle related drugs. Importantly, AKT activation in EMThigh-AKT not only enhances EMT process, but also contributes to the immunosuppressive microenvironment. By remodeling tumor immune-microenvironment by AKT inhibition, EMThigh-AKT can be treated by immune checkpoint blockade therapies. Meanwhile, we develop TumorMT website ( http://tumormt.neuroscience.org.cn/ ) to apply this EMT classification and provide reasonable therapeutic guidance.

Propionate reinforces epithelial identity and reduces aggressiveness of lung carcinoma

The epithelial-to-mesenchymal transition (EMT) plays a central role in the development of cancer metastasis and resistance to chemotherapy. However, its pharmacological treatment remains challenging. Here, we used an EMT-focused integrative functional genomic approach and identified an inverse association between short-chain fatty acids (propionate and butanoate) and EMT in non-small cell lung cancer (NSCLC) patients. Remarkably, treatment with propionate in vitro reinforced the epithelial transcriptional program promoting cell-to-cell contact and cell adhesion, while reducing the aggressive and chemo-resistant EMT phenotype in lung cancer cell lines. Propionate treatment also decreased the metastatic potential and limited lymph node spread in both nude mice and a genetic NSCLC mouse model. Further analysis revealed that chromatin remodeling through H3K27 acetylation (mediated by p300) is the mechanism underlying the shift toward an epithelial state upon propionate treatment. The results suggest that propionate administration has therapeutic potential in reducing NSCLC aggressiveness and warrants further clinical testing.

RSK3 switches cell fate: from stress-induced senescence to malignant progression

BACKGROUND

TGFβ induces several cell phenotypes including senescence, a stable cell cycle arrest accompanied by a secretory program, and epithelial-mesenchymal transition (EMT) in normal epithelial cells. During carcinogenesis cells lose the ability to undergo senescence in response to TGFβ but they maintain an EMT, which can contribute to tumor progression. Our aim was to identify mechanisms promoting TGFβ-induced senescence escape.

METHODS

In vitro experiments were performed with primary human mammary epithelial cells (HMEC) immortalized by hTert. For kinase library screen and modulation of gene expression retroviral transduction was used. To characterize gene expression, RNA microarray with GSEA analysis and RT-qPCR were used. For protein level and localization, Western blot and immunofluorescence were performed. For senescence characterization crystal violet assay, Senescence Associated-β-Galactosidase activity, EdU staining were conducted. To determine RSK3 partners FLAG-baited immunoprecipitation and mass spectrometry-based proteomic analyses were performed. Proteosome activity and proteasome enrichment assays were performed. To validate the role of RSK3 in human breast cancer, analysis of METABRIC database was performed. Murine intraductal xenografts using MCF10DCIS.com cells were carried out, with histological and immunofluorescence analysis of mouse tissue sections.

RESULTS

A screen with active kinases in HMECs upon TGFβ treatment identified that the serine threonine kinase RSK3, or RPS6KA2, a kinase mainly known to regulate cancer cell death including in breast cancer, reverted TGFβ-induced senescence. Interestingly, RSK3 expression decreased in response to TGFβ in a SMAD3-dependent manner, and its constitutive expression rescued SMAD3-induced senescence, indicating that a decrease in RSK3 itself contributes to TGFβ-induced senescence. Using transcriptomic analyses and affinity purification coupled to mass spectrometry-based proteomics, we unveiled that RSK3 regulates senescence by inhibiting the NF-κΒ pathway through the decrease in proteasome-mediated IκBα degradation. Strikingly, senescent TGFβ-treated HMECs display features of epithelial to mesenchymal transition (EMT) and during RSK3-induced senescence escaped HMECs conserve EMT features. Importantly, RSK3 expression is correlated with EMT and invasion, and inversely correlated with senescence and NF-κΒ in human claudin-low breast tumors and its expression enhances the formation of breast invasive tumors in the mouse mammary gland.

CONCLUSIONS

We conclude that RSK3 switches cell fate from senescence to malignancy in response to TGFβ signaling.

Stem Cell Origin of Cancer: Clinical Implications for Cancer Immunity and Immunotherapy

A simple way to understand the immune system is to separate the self from non-self. If it is self, the immune system tolerates and spares. If it is non-self, the immune system attacks and destroys. Consequently, if cancer has a stem cell origin and is a stem cell disease, we have a serious problem and a major dilemma with immunotherapy. Because many refractory cancers are more self than non-self, immunotherapy may become an uphill battle and pyrrhic victory in cancer care. In this article, we elucidate cancer immunity. We demonstrate for whom, with what, as well as when and how to apply immunotherapy in cancer care. We illustrate that a stem cell theory of cancer affects our perspectives and narratives of cancer. Without a pertinent theory about cancer's origin and nature, we may unwittingly perform misdirected cancer research and prescribe misguided cancer treatments. In the ongoing saga of immunotherapy, we are at a critical juncture. Because of the allure and promises of immunotherapy, we will be treating more patients not immediately threatened by their cancer. They may have more to lose than to gain, if we have a misconception and if we are on a wrong mission with immunotherapy. According to the stem cell theory of cancer, we should be careful with immunotherapy. When we do not know or realize that cancer originates from a stem cell and has stem-ness capabilities, we may cause more harm than good in some patients and fail to separate the truth from the myth about immunotherapy in cancer care.

Stochastic epithelial-mesenchymal transitions diversify non-cancerous lung cell behaviours

Epithelial-mesenchymal plasticity (EMP) is a hallmark of cancer. By enabling cells to shift between different morphological and functional states, EMP promotes invasion, metastasis and therapy resistance. We report that near-diploid non-cancerous human epithelial lung cells spontaneously shift along the EMP spectrum without genetic changes. Strikingly, more than half of single cell-derived clones adopt a mesenchymal morphology. We independently characterise epithelial-like and mesenchymal-like clones. Mesenchymal clones lose epithelial markers, display larger cell aspect ratios and lower motility, with mostly unaltered proliferation rates. Stemness marker expression and metabolic rewiring diverge independently of phenotypes. In 3D culture, more epithelial clones become mesenchymal-like. Thus, non-cancerous epithelial cells may acquire cancer metastasis-associated features prior to genetic alterations and cancerous transformation.

An androgen receptor regulated gene score is associated with epithelial to mesenchymal transition features in triple negative breast cancers

BACKGROUND

Androgen receptor (AR) is considered a marker of better prognosis in hormone receptor positive breast cancers (BC), however, its role in triple negative breast cancer (TNBC) is controversial. This may be attributed to intrinsic molecular differences or scoring methods for AR positivity. We derived AR regulated gene score and examined its utility in BC subtypes.

METHODS

AR regulated genes were derived by applying a bioinformatic pipeline on publicly available microarray data sets of AR+ BC cell lines and gene score was calculated as average expression of six AR regulated genes. Tumors were divided into AR high and low based on gene score and associations with clinical parameters, circulating androgens, survival and epithelial to mesenchymal transition (EMT) markers were examined, further evaluated in invitro models and public datasets.

RESULTS

53% (133/249) tumors were classified as AR gene score high and were associated with significantly better clinical parameters, disease-free survival (86.13 vs 72.69 months, log rank p = 0.032) when compared to AR low tumors. 36% of TNBC (N = 66) were AR gene score high with higher expression of EMT markers (p = 0.024) and had high intratumoral levels of 5α-reductase, enzyme involved in intracrine androgen metabolism. In MDA-MB-453 treated with dihydrotestosterone, SLUG expression increased, E-cadherin decreased with increase in migration and these changes were reversed with bicalutamide. Similar results were obtained in public datasets.

CONCLUSION

Deciphering the role of AR in BC is difficult based on AR protein levels alone. Our results support the context dependent function of AR in driving better prognosis in ER positive tumors and EMT features in TNBC tumors.

Intrinsic resistance to ROS1 inhibition in a patient with CD74-ROS1 mediated by AXL overexpression

BACKGROUND

The vast majority of patients with ROS1 positive non-small cell lung cancer (NSCLC) derive clinical benefit from currently approved ROS1 therapies, including crizotinib and entrectinib. However, a small proportion of patients treated with ROS1 inhibitors fail to derive any clinical benefit and demonstrate rapid disease progression. The biological mechanisms underpinning intrinsic resistance remain poorly understood for oncogene-driven cancers.

METHODS

We generated a patient-derived cell line, CUTO33, from a ROS1 therapy naive patient with CD74-ROS1+ NSCLC, who ultimately did not respond to a ROS1 inhibitor. We evaluated a panel of ROS1+ patient-derived NSCLC cell lines and used cell-based assays to determine the mechanism of intrinsic resistance to ROS1 therapy.

RESULTS

The CUTO33 cell line expressed the CD74-ROS1 gene fusion at the RNA and protein level. The ROS1 fusion protein was phosphorylated at baseline consistent with the known intrinsic activity of this oncogene. ROS1 phosphorylation could be inhibited using a wide array of ROS1 inhibitors, however these inhibitors did not block cell proliferation, confirming intrinsic resistance in this model and consistent with the patient's lack of response to a ROS1 inhibitor. CUTO33 expressed high levels of AXL, which has been associated with drug resistance. Combination of an AXL inhibitor or AXL knockdown with a ROS1 inhibitor partially reversed resistance.

CONCLUSIONS

In summary, we demonstrate that AXL overexpression is a mechanism of intrinsic resistance to ROS1 inhibitors.

Microdissection of cancer-associated fibroblast infiltration subtypes unveils the secreted SERPINE2 contributing to immunosuppressive microenvironment and immuotherapeutic resistance in gastric cancer: A large-scale study integrating bulk and single-cell transcriptome profiling

In the era of immunotherapy, the suboptimal response rate and the development of acquired resistance among the initial beneficiaries continue to present significant challenges across multiple malignancies, including gastric cancer (GC). Considering that the interactions of tumor stroma, especially the cancer-associated fibroblasts (CAFs), with immune and tumor cells, play indispensable roles in tumor progression, tumor microenvironment remodeling and therapeutic responsiveness, in-depth exploration on the roles of CAFs and pivotal mediators of their functions may provide novel clues to increase the effectiveness of current immunotherapeutic drugs and further achieve synergistic antitumor response. Herein, through the consensus clustering of canonical biomarkers, three GC subclasses with different abundance of CAFs were virtually microdissected in four integrated bulk cohorts encompassing 2148 GC patients from 11 independent datasets. An extensive immunogenomic analysis revealed that tumors with high CAFs infiltration were characterized with unfavorable outcomes, aggressive phenotypes, decreased tumor immunogenicity, high risk of immune evasion and thus immunotherapeutic resistance. By leveraging large-scale single-cell transcriptomic profiling, a series of CAF-secreted proteins were identified, among which the SERPINE2 was confirmed to be restrictively enriched in stromal fibroblasts of GC tissues and contribute to promoting a protumor milieu and fostering an immunosuppressive microenvironment via bioinformatics computations and tissue microarray analysis. Moreover, pan-cancer investigations generalized the immunological roles of SERPINE2, especially in pan-gastrointestinal malignancies, with multiple real-world immunotherapy cohorts further confirming its implications on predicting immunotherapeutic efficacy. In conclusion, these findings suggest that the CAF-derived SERPINE2 is a promising immune-oncology target with therapeutic implications to further synergize the immunotherapeutic combinations.

Comprehensive profiling identifies tumour and immune microenvironmental differences in clinical subsets of cutaneous squamous cell carcinoma

BACKGROUND

Cutaneous squamous cell carcinomas (cSCCs) are the second most diagnosed skin cancer worldwide; however, little is known about the pathobiological factors that contribute to the diverse clinical outcomes seen.

OBJECTIVES

To profile cSCCs comprehensively and identify the pathological processes that contribute to the disparities seen in their clinical behaviour.

METHODS

We characterized the genomic, transcriptomic and immunohistochemical profiles of 211 cSCC tumours, including 37 cSCCs from immunocompromised patients.

RESULTS

cSCCs from immunocompromised patients were characterized by a lack of B cells in the peritumoral stroma compared with immunocompetent patients. Further, an abundance of a memory B-cell-like population in the peritumoral stroma was associated with a better prognosis in all patients (immunocompetent and immunocompromised), as well as only immunocompetent patients. No differences in genetic -variants, tumour mutational burden or mutational signatures were observed between cSCCs from immunocompetent and immunocompromised patients. Thus, differences in survival between cSCCs from immunocompromised patients and immunocompetent patients are not likely to be driven by tumour genomic factors, but may be associated with differential host immune response. cSCC not from a primary head and neck site had lower tumour mutational burden and exhibited upregulation of the epithelial-mesenchymal transition programme compared with head and neck cSCC. Both factors were implicated with poorer responses to immune checkpoint inhibition, and the latter with poorer survival.

CONCLUSIONS

We identified tumour and host immune factors that contribute to the disparate clinical behaviour of cSCC, with broad translational application, including prognostication, treatment prediction to current therapies and the identification of novel anticancer therapy approaches in cSCC.

Tumor Microenvironment Modifications Induced by Afatinib in Squamous Cell Carcinoma of the Head and Neck: A Window-of-Opportunity Study (EORTC-90111-24111)

PURPOSE

The EORTC-90111-24111 phase II window study evaluated afatinib versus no preoperative treatment in patients with primary squamous cell carcinoma of the head and neck (HNSCC). We investigated afatinib-induced tumor and microenvironment modifications by comparing pre- and posttreatment tumor biopsies.

PATIENTS AND METHODS

Thirty treatment-naïve patients with primary HNSCC were randomized. Twenty-five patients received afatinib for 14 days before surgery (40 mg 1×/day) and 5 patients were attributed to the control arm. Biopsies were taken at work-up and during surgery. Good quality RNA samples were used for omics analyses. The control arm was enlarged by samples coming from our previous similar window study.

RESULTS

IHC analyses of afatinib-treated tumor biopsies showed a decrease in pEGFR (P ≤ 0.05) and pERK (P ≤ 0.05); and an increase in CD3+ (P ≤ 0.01) and CD8+ (P ≤ 0.01) T-cell infiltration, and in CD3+ (P ≤ 0.05) T-cell density. RNA sequencing analyses of afatinib-treated tumor samples showed upregulation of inflammatory genes and increased expression scores of signatures predictive of response to programmed cell death protein 1 blockade (P ≤ 0.05). In posttreatment biopsies of afatinib-treated patients, two clusters were observed. Cluster 1 showed a higher expression of markers and gene sets implicated in epithelial-to-mesenchymal transition (EMT) and activation of cancer-associated fibroblasts (CAF) compared with cluster 2 and controls.

CONCLUSIONS

Short-term treatment with afatinib in primary HNSCC induces CD3+ and CD8+ tumor infiltration and, in some patients, EMT and CAF activation. These results open perspectives to overcome resistance mechanisms to anti-HER therapy and to potentiate the activity of immune checkpoint inhibitors.

TIM-3 Expression and M2 Polarization of Macrophages in the TGFβ-Activated Tumor Microenvironment in Colorectal Cancer

TGFβ signaling in the tumor microenvironment (TME) drives immune evasion and is a negative predictor of immune checkpoint inhibitor (ICI) efficacy in colorectal cancer (CRC). TIM-3, an inhibitory receptor implicated in anti-tumor immune responses and ICI resistance, has emerged as an immunotherapeutic target. This study investigated TIM-3, M2 macrophages and the TGFβ-activated TME, in association with microsatellite instability (MSI) status and consensus molecular subtypes (CMSs). Transcriptomic cohorts of CRC tissues, organoids and xenografts were examined (n = 2240). TIM-3 and a TGFβ-inducible stromal protein, VCAN, were evaluated in CRC specimens using immunohistochemistry (n = 45). TIM-3 expression on monocytes and generated M2 macrophages was examined by flow cytometry. We found that the expression of HAVCR2 (TIM-3) significantly correlated with the transcriptional signatures of TGFβ, TGFβ-dependent stromal activation and M2 macrophage, each of which were co-upregulated in CMS4, CMS1 and MSI CRCs across all datasets. Tumor-infiltrating TIM-3+ immune cells accumulated in TGFβ-responsive cancer stroma. TIM-3 was increased on M2-polarized macrophages, and on monocytes in response to TGFβ treatment. In conclusion, we identified a close association between TIM-3 and M2-like polarization of macrophages in the TGFβ-rich TME. Our findings provide new insights into personalized immunotherapeutic strategies based on the TME for CRCs.

A prognostic matrix code defines functional glioblastoma phenotypes and niches

Interactions among tumor, immune and vascular niches play major roles in driving glioblastoma (GBM) malignancy and treatment responses. The composition, heterogeneity, and localization of extracellular core matrix proteins (CMPs) that mediate such interactions, however, are not well understood. Here, we characterize functional and clinical relevance of genes encoding CMPs in GBM at bulk, single cell, and spatial anatomical resolution. We identify a "matrix code" for genes encoding CMPs whose expression levels categorize GBM tumors into matrisome-high and matrisome-low groups that correlate with worse and better patient survival, respectively. The matrisome enrichment is associated with specific driver oncogenic alterations, mesenchymal state, infiltration of pro-tumor immune cells and immune checkpoint gene expression. Anatomical and single cell transcriptome analyses indicate that matrisome gene expression is enriched in vascular and leading edge/infiltrative anatomic structures that are known to harbor glioma stem cells driving GBM progression. Finally, we identified a 17-gene matrisome signature that retains and further refines the prognostic value of genes encoding CMPs and, importantly, potentially predicts responses to PD1 blockade in clinical trials for GBM. The matrisome gene expression profiles provide potential biomarkers of functionally relevant GBM niches that contribute to mesenchymal-immune cross talk and patient stratification which could be applied to optimize treatment responses.

Construction of PANoptosis signature: Novel target discovery for prostate cancer immunotherapy

PANoptosis pathway gene sets encompassing pyroptosis, apoptosis, and necroptosis were identified from the MSigDB database. We analyzed the perturbations and crosstalk in the PANoptosis pathway in prostate adenocarcinoma (PRAD), including gene mutation, transcription, methylation, and clinical features. By constructing a PANoptosis signature, we accurately predicted the prognosis and immunotherapeutic response of PRAD patients. We further explored the molecular features and immunological roles of the signature, dividing patients into high- and low-score groups. Notably, the high-score group correlated with better survival outcomes and immunotherapeutic responses, as well as a higher mutation frequency and enrichment score in the PANoptosis and HALLMARK pathways. The PANoptosis signature also enhanced overall antitumor immunity, promoted immune cell infiltration, upregulated immune checkpoint regulators, and revealed the cold tumor characteristics of PRAD. We also identified potential drug targets based on the PANoptosis signature. These findings lead the way in identifying novel prognostic markers and therapeutic targets for patients with PRAD.

Three-Course Neoadjuvant Chemotherapy Associated with Unfavorable Survival of Non-responders to the First Two Courses for Locally Advanced Esophageal Cancer

BACKGROUND

Three-course neoadjuvant chemotherapy (NAC) followed by surgery has become a standard of care for locally advanced esophageal cancer (EC). However, some patients occasionally experience a poor tumor response to the third course and have a poor clinical outcome.

METHODS

An exploratory analysis of data from the authors' recent multicenter randomized phase 2 trial compared patients with locally advanced EC who received two courses (n = 78) and those who received three courses (n = 68) of NAC. The association between tumor response and clinico-pathologic factors, including survival, was evaluated to identify risk factors in the three-course group.

RESULTS

Of 68 patients who received three courses of NAC, 28 (41.2%) had a tumor reduction rate lower than 10% during the third course. This rate was associated with unfavorable overall survival (OS) and progression-free survival (PFS) compared with a tumor reduction rate of 10% or higher (2-year OS rate: 63.5% vs. 89.3%, P = 0.007; 2-year PFS rate: 52.6% vs. 79.7%, P = 0.020). The independent prognostic factors for OS were tumor reduction rate lower than 10% during the third course (hazard ratio [HR], 2.735; 95% confidence interval [CI] 1.041-7.188; P = 0.041) and age of 65 years or older (HR, 9.557, 95% CI 1.240-73.63; P = 0.030). Receiver operating characteristic curve and multivariable logistic regression analyses identified a tumor reduction rate lower than 50% after the first two courses as an independent predictor of a tumor reduction rate lower than 10% during the third course of NAC (HR, 4.315; 95% CI 1.329-14.02; P = 0.015).

CONCLUSION

Continuing NAC through a third course may worsen survival for patients who do not experience a response to the first two courses in locally advanced EC.

What role does PDL1 play in EMT changes in tumors and fibrosis?

Epithelial-mesenchymal transformation (EMT) plays a pivotal role in embryonic development, tissue fibrosis, repair, and tumor invasiveness. Emerging studies have highlighted the close association between EMT and immune checkpoint molecules, particularly programmed cell death ligand 1 (PDL1). PDL1 exerts its influence on EMT through bidirectional regulation. EMT-associated factors, such as YB1, enhance PDL1 expression by directly binding to its promoter. Conversely, PDL1 signaling triggers downstream pathways like PI3K/AKT and MAPK, promoting EMT and facilitating cancer cell migration and invasion. Targeting PDL1 holds promise as a therapeutic strategy for EMT-related diseases, including cancer and fibrosis. Indeed, PDL1 inhibitors, such as pembrolizumab and nivolumab, have shown promising results in clinical trials for various cancers. Recent research has also indicated their potential benefit in fibrosis treatment in reducing fibroblast activation and extracellular matrix deposition, thereby addressing fibrosis. In this review, we examine the multifaceted role of PDL1 in immunomodulation, growth, and fibrosis promotion. We discuss the challenges, mechanisms, and clinical observations related to PDL1, including the limitations of the PD1/PDL1 axis in treatment and PD1-independent intrinsic PDL1 signaling. Our study highlights the dynamic changes in PDL1 expression during the EMT process across various tumor types. Through interplay between PDL1 and EMT, we uncover co-directional alterations, regulatory pathways, and diverse changes resulting from PDL1 intervention in oncology. Additionally, our findings emphasize the dual role of PDL1 in promoting fibrosis and modulating immune responses across multiple diseases, with potential implications for therapeutic approaches. We particularly investigate the therapeutic potential of targeting PDL1 in type II EMT fibrosis: strike balance between fibrosis modulation and immune response regulation. This analysis provides valuable insights into the multifaceted functions of PDL1 and contributes to our understanding of its complex mechanisms and therapeutic implications.

The Transcriptional and Epigenetic Landscape of Cancer Cell Lineage Plasticity

Lineage plasticity, a process whereby cells change their phenotype to take on a different molecular and/or histologic identity, is a key driver of cancer progression and therapy resistance. Although underlying genetic changes within the tumor can enhance lineage plasticity, it is predominantly a dynamic process controlled by transcriptional and epigenetic dysregulation. This review explores the transcriptional and epigenetic regulators of lineage plasticity and their interplay with other features of malignancy, such as dysregulated metabolism, the tumor microenvironment, and immune evasion. We also discuss strategies for the detection and treatment of highly plastic tumors.

SIGNIFICANCE

Lineage plasticity is a hallmark of cancer and a critical facilitator of other oncogenic features such as metastasis, therapy resistance, dysregulated metabolism, and immune evasion. It is essential that the molecular mechanisms of lineage plasticity are elucidated to enable the development of strategies to effectively target this phenomenon. In this review, we describe key transcriptional and epigenetic regulators of cancer cell plasticity, in the process highlighting therapeutic approaches that may be harnessed for patient benefit.

Pharmacological targeting of netrin-1 inhibits EMT in cancer

Epithelial-to-mesenchymal transition (EMT) regulates tumour initiation, progression, metastasis and resistance to anti-cancer therapy1-7. Although great progress has been made in understanding the role of EMT and its regulatory mechanisms in cancer, no therapeutic strategy to pharmacologically target EMT has been identified. Here we found that netrin-1 is upregulated in a primary mouse model of skin squamous cell carcinoma (SCC) exhibiting spontaneous EMT. Pharmacological inhibition of netrin-1 by administration of NP137, a netrin-1-blocking monoclonal antibody currently used in clinical trials in human cancer (ClinicalTrials.gov identifier NCT02977195 ), decreased the proportion of EMT tumour cells in skin SCC, decreased the number of metastases and increased the sensitivity of tumour cells to chemotherapy. Single-cell RNA sequencing revealed the presence of different EMT states, including epithelial, early and late hybrid EMT, and full EMT states, in control SCC. By contrast, administration of NP137 prevented the progression of cancer cells towards a late EMT state and sustained tumour epithelial states. Short hairpin RNA knockdown of netrin-1 and its receptor UNC5B in EPCAM+ tumour cells inhibited EMT in vitro in the absence of stromal cells and regulated a common gene signature that promotes tumour epithelial state and restricts EMT. To assess the relevance of these findings to human cancers, we treated mice transplanted with the A549 human cancer cell line-which undergoes EMT following TGFβ1 administration8,9-with NP137. Netrin-1 inhibition decreased EMT in these transplanted A549 cells. Together, our results identify a pharmacological strategy for targeting EMT in cancer, opening up novel therapeutic interventions for anti-cancer therapy.

Netrin-1 blockade inhibits tumour growth and EMT features in endometrial cancer

Netrin-1 is upregulated in cancers as a protumoural mechanism1. Here we describe netrin-1 upregulation in a majority of human endometrial carcinomas (ECs) and demonstrate that netrin-1 blockade, using an anti-netrin-1 antibody (NP137), is effective in reduction of tumour progression in an EC mouse model. We next examined the efficacy of NP137, as a first-in-class single agent, in a Phase I trial comprising 14 patients with advanced EC. As best response we observed 8 stable disease (8 out of 14, 57.1%) and 1 objective response as RECIST v.1.1 (partial response, 1 out of 14 (7.1%), 51.16% reduction in target lesions at 6 weeks and up to 54.65% reduction during the following 6 months). To evaluate the NP137 mechanism of action, mouse tumour gene profiling was performed, and we observed, in addition to cell death induction, that NP137 inhibited epithelial-to-mesenchymal transition (EMT). By performing bulk RNA sequencing (RNA-seq), spatial transcriptomics and single-cell RNA-seq on paired pre- and on-treatment biopsies from patients with EC from the NP137 trial, we noted a net reduction in tumour EMT. This was associated with changes in immune infiltrate and increased interactions between cancer cells and the tumour microenvironment. Given the importance of EMT in resistance to current standards of care2, we show in the EC mouse model that a combination of NP137 with carboplatin-paclitaxel outperformed carboplatin-paclitaxel alone. Our results identify netrin-1 blockade as a clinical strategy triggering both tumour debulking and EMT inhibition, thus potentially alleviating resistance to standard treatments.

Predicting Outcomes of Atezolizumab and Bevacizumab Treatment in Patients with Hepatocellular Carcinoma

A combination of atezolizumab with bevacizumab (AB) is the first regimen that has shown superiority compared to sorafenib and is now being used as the systemic treatment of choice for hepatocellular carcinoma (HCC) patients with Barcelona Liver Cancer Clinic stage C. However, a considerable number of patients do not achieve survival or significant responses, indicating the need to identify predictive biomarkers for initial and on-treatment decisions in HCC patients receiving AB. In this manuscript, we summarized the current data from both experimental and clinical studies. This review will be beneficial for both clinicians and researchers in clinical practice as well as those designing experimental, translational, or clinical studies.

Immune Checkpoint Blockade Therapy for Breast Cancer: Lessons from Epithelial-Mesenchymal Transition

Immune checkpoint blockade therapies have generated efficacious responses in certain tumor types; however, the responses of breast carcinomas have been largely limited. Moreover, the identity of various parameters that can predict responses to immunotherapies, and at the same time, serve as putative biomarkers that can be therapeutically targeted to enhance the effectiveness of immunotherapies for breast cancers, remains to be comprehensively delineated. Activation of epithelial-mesenchymal plasticity in cancer cells, including those of the breast, increases their tumor-initiating potential and promotes their aggressiveness and resistance to multiple treatment regimens. Moreover, the residence of cancer cells in alternating epithelial or mesenchymal plastic phenotypic states can also influence their immuno-modulatory properties and susceptibilities to immune checkpoint blockade therapies. In this current opinion, we discuss the lessons that can be learnt from epithelial-mesenchymal transition to potentiate the efficacy of immunotherapy for breast cancers. We also discuss strategies to sensitize more-mesenchymal cancer cells to anti-tumor immunity and immune checkpoint blockade therapies, with the hope that these can serve as new translational avenues for the treatment of human breast tumors.

The Role of EMT-Related lncRNAs in Ovarian Cancer

Ovarian cancer (OC) is one of the deadliest cancers worldwide; late diagnosis and drug resistance are two major factors often responsible for high morbidity and treatment failure. Epithelial-to-mesenchymal transition (EMT) is a dynamic process that has been closely linked with cancer. Long non-coding RNAs (lncRNAs) have been also associated with several cancer-related mechanisms, including EMT. We conducted a literature search in the PubMed database in order to sum up and discuss the role of lncRNAs in regulating OC-related EMT and their underlying mechanisms. Seventy (70) original research articles were identified, as of 23 April 2023. Our review concluded that the dysregulation of lncRNAs is highly associated with EMT-mediated OC progression. A comprehensive understanding of lncRNAs' mechanisms in OC will help in identifying novel and sensitive biomarkers and therapeutic targets for this malignancy.

Unpacking the Complexity of Epithelial Plasticity: From Master Regulator Transcription Factors to Non-Coding RNAs

Cellular plasticity in cancer enables adaptation to selective pressures and stress imposed by the tumor microenvironment. This plasticity facilitates the remodeling of cancer cell phenotype and function (such as tumor stemness, metastasis, chemo/radio resistance), and the reprogramming of the surrounding tumor microenvironment to enable immune evasion. Epithelial plasticity is one form of cellular plasticity, which is intrinsically linked with epithelial-mesenchymal transition (EMT). Traditionally, EMT has been regarded as a binary state. Yet, increasing evidence suggests that EMT involves a spectrum of quasi-epithelial and quasi-mesenchymal phenotypes governed by complex interactions between cellular metabolism, transcriptome regulation, and epigenetic mechanisms. Herein, we review the complex cross-talk between the different layers of epithelial plasticity in cancer, encompassing the core layer of transcription factors, their interacting epigenetic modifiers and non-coding RNAs, and the manipulation of cancer immunogenicity in transitioning between epithelial and mesenchymal states. In examining these factors, we provide insights into promising therapeutic avenues and potential anti-cancer targets.

A prognostic matrix code defines functional glioblastoma phenotypes and niches

Interactions among tumor, immune and vascular niches play major roles in driving glioblastoma (GBM) malignancy and treatment responses. The composition, heterogeneity, and localization of extracellular core matrix proteins (CMPs) that mediate such interactions, however, are not well understood. Here, we characterize functional and clinical relevance of genes encoding CMPs in GBM at bulk, single cell, and spatial anatomical resolution. We identify a "matrix code" for genes encoding CMPs whose expression levels categorize GBM tumors into matrisome-high and matrisome-low groups that correlate with worse and better survival, respectively, of patients. The matrisome enrichment is associated with specific driver oncogenic alterations, mesenchymal state, infiltration of pro-tumor immune cells and immune checkpoint gene expression. Anatomical and single cell transcriptome analyses indicate that matrisome gene expression is enriched in vascular and leading edge/infiltrative anatomic structures that are known to harbor glioma stem cells driving GBM progression. Finally, we identified a 17-gene matrisome signature that retains and further refines the prognostic value of genes encoding CMPs and, importantly, potentially predicts responses to PD1 blockade in clinical trials for GBM. The matrisome gene expression profiles may provide biomarkers of functionally relevant GBM niches that contribute to mesenchymal-immune cross talk and patient stratification to optimize treatment responses.

Single-cell RNA sequencing reveals microenvironment context-specific routes for epithelial-mesenchymal transition in pancreas cancer cells

In the PDAC tumor microenvironment, multiple factors initiate the epithelial-mesenchymal transition (EMT) that occurs heterogeneously among transformed ductal cells, but it is unclear if different drivers promote EMT through common or distinct signaling pathways. Here, we use single-cell RNA sequencing (scRNA-seq) to identify the transcriptional basis for EMT in pancreas cancer cells in response to hypoxia or EMT-inducing growth factors. Using clustering and gene set enrichment analysis, we find EMT gene expression patterns that are unique to the hypoxia or growth factor conditions or that are common between them. Among the inferences from the analysis, we find that the FAT1 cell adhesion protein is enriched in epithelial cells and suppresses EMT. Further, the receptor tyrosine kinase AXL is preferentially expressed in hypoxic mesenchymal cells in a manner correlating with YAP nuclear localization, which is suppressed by FAT1 expression. AXL inhibition prevents EMT in response to hypoxia but not growth factors. Relationships between FAT1 or AXL expression with EMT were confirmed through analysis of patient tumor scRNA-seq data. Further exploration of inferences from this unique dataset will reveal additional microenvironment context-specific signaling pathways for EMT that may represent novel drug targets for PDAC combination therapies.