SEARCH · PubMed Health

Results for “tumor microenvironment (TME)”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2026

Pan-cancer characterization of HMGA1 reveals its oncogenic role in tumor microenvironment and stemness: functional validation in pancreatic cancer migration and invasion.

BACKGROUND: HMGA1 is a chromatin-associated oncogenic factor implicated in tumor progression, epithelial-mesenchymal transition (EMT), stemness, and metastasis. However, its pan-cancer expression and prognostic patterns, epigenetic activation, and relationship with malignant-cell stemness/plasticity and tumor microenvironment (TME) remodeling in pancreatic adenocarcinoma (PAAD) remain incompletely defined. This study aimed to systematically characterize HMGA1 across cancers and clarify its clinical and biological relevance in PAAD. METHODS: Pan-cancer transcriptomic, clinical, genetic, methylation, immune, and stemness data were integrated from multiple public databases. PAAD single-cell RNA sequencing data were analyzed to localize HMGA1 expression, infer malignant-cell pseudotime, calculate a stemness module score, and assess ligand-receptor communication using CellChat. Public HMGA1-knockdown RNA sequencing data were reanalyzed to evaluate transcriptional remodeling. The Cancer Genome Atlas (TCGA)-PAAD expression and methylation data were used to assess TME-remodeling, immune-suppression, stemness/plasticity, cytokine/chemokine, checkpoint, and promoter-methylation features. HMGA1 expression and function were further examined using immunohistochemistry (IHC), quantitative real-time polymerase chain reaction, Western blotting, wound-healing assays, and Transwell migration and invasion assays. RESULTS: HMGA1 was upregulated in most tumor types, and high expression was associated with unfavorable survival in multiple cancers, including PAAD. In PAAD, HMGA1 was enriched in malignant epithelial cells and positively correlated with pseudotime (Spearman's rho =0.594), while the stemness module score increased along pseudotime (rho =0.748). HMGA1-high malignant cells showed markedly stronger CellChat-inferred outgoing communication, predominantly involving extracellular matrix (ECM)-receptor, adhesion-related, and selected immunomodulatory ligand-receptor axes. HMGA1 knockdown was associated with broad remodeling of EMT, TGF-β, Hedgehog, IL6/JAK/STAT3, KRAS, and cancer stem cell/stemness-related programs rather than uniform suppression of these programs. HMGA1 promoter methylation was inversely correlated with HMGA1 expression (rho =-0.633) and the TME-remodeling score (rho =-0.347). HMGA1 was associated with selected mediators, including PPIA, PLAU, ANXA1, LGALS9, TGFB1, CD276, and CD47, but not with a generalized checkpoint-high phenotype. Functionally, HMGA1 knockdown significantly reduced pancreatic cancer (PC) cell migration and invasion. CONCLUSIONS: These findings support an association-based model in which promoter hypomethylation-associated HMGA1 activation is linked to malignant epithelial stemness/plasticity, ECM/adhesion-dominant TME remodeling, selected immunomodulatory programs, and aggressive PAAD phenotypes. Further mechanistic and clinical validation is required before HMGA1 can be used for therapeutic stratification or immunotherapy-response prediction.

HMGA1↗

PubMed2025

Digital Immunophenotyping of Lung Atypical Carcinoids and Large Cell Neuroendocrine Carcinomas Identifies Three Subtypes With Specific Tumor-Immune Microenvironment Features.

Atypical carcinoids (ACs) and large cell neuroendocrine carcinomas (LCNECs) are defined by the WHO as intermediate- and high-grade lung neuroendocrine neoplasms, respectively, based on morphological criteria; however, treatment strategies remain debated. Given the emerging role of the tumor microenvironment (TME) and tumor-infiltrating lymphocytes (TILs) in cancer prognosis and therapy response, this study aimed to characterize the immune landscape of ACs and LCNECs comprehensively. Immunohistochemistry for T-cell markers (CD3, CD8), immune checkpoints (PD-1, PD-L1), HLA molecules (HLA-DR, HLA-I), and fibroblasts (α-SMA) was performed on a re-evaluated cohort of 56 ACs and 104 LCNECs. Digital image analysis quantified intra-tumor (iTILs) and stromal (sTILs) CD3 and CD8 TILs in the whole slide and in specific tumor regions (invasive margin [IM] and central tumor [CT]). LCNECs exhibited significantly higher stromal T-cell infiltration, immune checkpoint expression, and HLA compared to ACs (p < 0.001), while α-SMA was more prominent in ACs. No ACs showed PD-L1 tumor expression. Digital quantification confirmed greater iTILs and sTILs in LCNECs across all regions, with moderate concordance to manual counts. Interestingly, TIL parameters were higher at the IM than in the CT (p < 0.001). Using Boruta feature selection algorithm, Principal Component Analysis and Hierarchical Clustering, three patient clusters were identified: Cluster 1 (mainly ACs, low TILs, favorable prognosis), Cluster 2 (mixed histology, intermediate TILs, moderate prognosis), and Cluster 3 (mostly LCNECs, high TILs, poor prognosis), with distinct TME marker profiles. PD-L1 tumor expression was strongly linked to Cluster 3. These findings suggest that ACs and LCNECs may be stratified into three distinct immune clusters, highlighting the heterogeneity of their tumor microenvironment and providing a rationale for further translational studies.

Humans↗

PubMed2026

Oxidative stress and cancer: current insights and therapeutic implications.

OXIDATIVE STRESS: good or evil? Oxidative stress occurs when the balance between reactive oxygen species (ROS) and antioxidant defenses shifts toward an excess of ROS; while essential in physiological processes, it plays a context-dependent role in cancer, contributing to both the promotion and inhibition of tumorigenesis. Small to moderate amounts of ROS activate pathways supporting tumor progression and proliferation, while large amounts lead to genomic instability and cell death. ROS are generated endogenously and exogenously. In cancer, ROS activate pathways that prompt tumor development (KRAS, MYC, PI3K-Akt-mTOR) and block tumor suppressors (p53, BRCA1), allowing tumorigenesis and drug resistance. They also modulate the tumor microenvironment (TME) by altering tumor, stromal and immune cell interactions, which initiate angiogenesis, epithelial-mesenchymal transition (EMT), inflammation and metastasis. Myeloid-derived suppressor cells (MDSCs) and cancer-associated fibroblasts (CAFs) contribute to ROS-driven immunosuppression. Cancer cells mainly rely on glycolysis and oxidative phosphorylation (OXPHOS) to sustain their energetic and metabolic requirements. Generated ROS act as metabolic byproducts and signaling molecules supporting proliferation and tumorigenesis. Cancer stem cells (CSCs) produce low ROS levels by activating antioxidant pathways and mitochondria remodeling, ensuring recurrence and persistence. There is a redox duality that presents challenges and opportunities for therapies. Pro-oxidant approaches attempt to overwhelm the tumor's defenses, while antioxidants preserve healthy tissues. Advances in targeted redox modulation with immunotherapies improve therapy effectiveness. We propose a new "Adaptive Directed Redox Therapy" (ADRT), which involves a dynamic, feedback-controlled methodology that alternates pro- and antioxidant phases to selectively collapse tumor redox balance while preserving normal tissues.

Humans↗

PubMed2025

Integrative multi-omics analysis proposes a metabolic classification of gliomas: distinct metabolic states, immune infiltration, and prognosis.

BACKGROUND: The tumor microenvironment (TME) of glioma harbors diverse cell types; however, cell metabolic heterogeneity remains to be explored. This study aims to characterize the metabolic features of different cell types in the TME by integrating multiple datasets, including genomics, bulk and single-cell transcriptomics, and metabolomics. METHODS: Unsupervised machine learning was used to construct an energy metabolic classifier based on the metabolic pathways identified from bulk RNA-seq of gliomas in the TCGA dataset. The classifier was externally validated using multiple datasets, including genomics, bulk RNA-seq, snRNA-seq, and the metabolomics data. Furthermore, metabolic heterogeneity associated with the classifier was further characterized at single-cell resolution. RESULTS: The energy metabolism-based classifier stratified patients into two prognostic clusters: patients in cluster 1 were characterized by high pathway activity of glycolysis, the pentose phosphate pathway (PPP), and fatty acid oxidation (FAO), whereas patients in cluster 2 exhibited higher activity in glutaminolysis. This metabolic classifier revealed both intratumoral and intertumoral metabolic heterogeneity, and the complexity was further validated by the metabolomics profiling and snRNA-seq data from the CPTAC dataset. Notably, OSMR, highly expressed in cluster 1, showed significant co-expression with key glycolytic enzyme genes. The OSM/OSMR/JAK1/STAT3 axis potently drives malignant progression of glioma cells, specially enhancing their invasive and migratory capabilities. Single-cell resolution analyses demonstrated that tumor metabolic heterogeneity is primarily driven by malignant cells rather than non-malignant components, while tumor microenvironment (TME) factors were also found to modulate malignant cell metabolism. Significantly, glycolytic activity in glioma cells increased during the phenotypic transition from PN (proneural) to MES (mesenchymal), with cluster 1 metabolic phenotypes predominating in the tumor core. Compared to cluster 2, cluster 1 patients exhibited higher mRNA expression of immunosuppressive checkpoint genes, which correlated with pronounced immunosuppression in the TME. Furthermore, various immune cells demonstrated distinct metabolic preferences at single-cell resolution. CONCLUSIONS: This study developed an energy metabolic-based classifier for gliomas with prognostic and therapeutic potential. Metabolic reprogramming was linked with the PN-to-MES transition of glioma cells and immunosuppression in the tumor microenvironment. Multi-omics data, especially snRNA-seq, offered insights into metabolism heterogeneity at single-cell resolution, enabling personalized treatment strategies.

Humans↗

PubMed2026

Oncogenic Mutations and Tumor Microenvironment Alterations in Diffuse Large B-Cell Lymphoma With Bulky Disease.

BACKGROUND: Bulky disease represents a clinically aggressive subset of diffuse large B-cell lymphoma (DLBCL) associated with adverse clinical outcomes. The aim of this study was to investigate the influence of oncogenic mutations and tumor microenvironment alterations on bulky disease in DLBCL. METHODS: We analyzed a cohort of 939 patients with newly diagnosed DLBCL. Using DNA (n = 934) and RNA (n = 524) sequencing, we compared oncogenic mutations and tumor microenvironment (TME) alterations based on tumor diameter, with cutoff values at 5.0 cm and 10.0 cm. Further stratification by mutations in key genes (CD58, STAT6, EBF1) correlated with tumor diameter revealed distinct transcriptomic and immunologic profiles. Subsequent single-cell RNA sequencing, guided by these mutational signatures, resolved the cellular heterogeneity within the TME. RESULTS: Integrative analysis revealed that tumor diameter correlated with increased incidence of mutations in CD58, STAT6, and EBF1; adverse genetic subtypes such as EZB-like MYC+ and TP53Mut; activation of oncogenic pathways (JAK/STAT, BCR, PI3K, and MYC); and an immunosuppressive tumor microenvironment. Notably, immune checkpoint molecules varied across the bulky stages, with CTLA-4, TIGIT, ICOS, and CD28 expression inversely correlated with tumor diameter, while CD70 and 4-1BBL expression positively correlated. Single-cell RNA sequencing further revealed mutation-specific tumor microenvironment insights. CD58-mutated tumor exhibited a profoundly immune-deserted microenvironment dominated by malignant B cells with minimal immune infiltration, whereas STAT6-mutated tumor was associated with increased fibroblasts and CD4 + T cells, particularly regulatory T cells (Treg) and Th1-like cells; EBF1-mutated tumor was characterized by increased proportions of malignant B cells. CONCLUSIONS: Collectively, our findings highlight the biological complexity of bulky disease, identifying candidate molecular targets and providing a biological framework for future therapeutic hypothesis generation in this clinically aggressive subset of DLBCL.

Humans↗

PubMed2026

Context-Dependent cGAS-STING Activation Shapes Metastatic Progression and Dormancy.

Cancer cells survive, proliferate, and metastasize in part because the immune system fails to detect and eliminate them. Moreover, the tumor microenvironment (TME) that surrounds the tumor supports cancer cell survival and resistance to chemo- and immunotherapies by inhibiting antitumor immune responses and thereby reducing the efficacy of immunotherapeutic interventions. cGAS-STING signaling senses cytoplasmic DNA and coordinates innate immune responses that shape tumor-intrinsic outcomes and the TME. Emerging evidence reveals a context-dependent, dualistic role for cGAS-STING in metastatic progression and cancer dormancy. Acute, robust activation in antigen-presenting cells promotes type I interferon responses, leading to suppression of tumor growth. By contrast, chronic, low-level cancer-intrinsic STING signaling can engage inflammatory programs that foster immune suppression and therapy resistance. Dormant disseminated tumor cells exploit niche cues to downregulate STING signaling and evade immune detection, whereas reactivation of dormant cells often involves restoration of STING activity that can promote immune elimination. In this article, we review mechanisms linking genome instability and cytoplasmic DNA to STING activation, summarize evidence for tumor-suppressive versus tumor-promoting functions across metastatic niches, and discuss how STING agonists and combination strategies may be optimized to maximize antitumor immunity while avoiding protumorigenic effects.

Humans↗

PubMed2025

Development and validation of a novel risk stratification signature derived from migrasome and tumor microenvironment-related genes for molecular subtyping and improving clinical outcomes in head and neck squamous cell carcinoma.

BACKGROUND: The tumor microenvironment (TME) and migrasomes released by tumor cells significantly influence carcinogenesis and immune evasion. However, our understanding of the prognostic and therapeutic implications of migrasome and tumor microenvironment-related genes (mtmRGs) in head and neck squamous cell carcinoma (HNSCC) remains limited. METHODS: We explored the relationship between mtmRGs and HNSCC prognosis by utilizing The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. Subsequently, we developed an innovative prognostic signature, and assessed its prognostic significance using the Kaplan-Meier method, time-dependent receiver operating characteristic (ROC), and Cox regression analyses. To explore the underlying mechanisms, we conducted gene set variation analysis (GSVA), gene set enrichment analysis (GESA), and immune infiltration analysis. A nomogram was developed to estimate the overall survival (OS) rates for HNSCC patients. Lastly, we chose P4HA1, which was part of the signature, for additional experimental validation in vitro and in vivo. RESULTS: The mtmRGs signature effectively classifies HNSCC patients into two distinct risk subgroups, with the high-risk cohort demonstrating significantly poorer OS. The risk score serves as an independent prognostic factor for HNSCC patients; those with lower risk scores are more likely to exhibit favorable responses to immunotherapy, particularly with CTLA4 inhibitors. Furthermore, a lower risk score is significantly correlated with the sensitivity of HNSCC patients to cyclophosphamide, gemcitabine, and axitinib. CONCLUSION: This study presents an innovative gene signature associated with mtmRGs, which may be utilized both for predicting survival and directing personalized chemotherapy and immunotherapy regiments for patients with HNSCC.

Humans↗

PubMed2026

Tumor-associated macrophages display differential protein cargo sorting in extracellular vesicles associated with poor survival in ovarian cancer.

Ovarian cancer (OC) progression and metastasis are promoted by ascites, which constitutes a central part of the tumor microenvironment (TME). In this fluid, tumor-associated macrophages (TAMs) represent a prominent immune cell type. In addition to tumor and other host cells such as TAMs, ascites is highly enriched in soluble factors as well as extracellular vesicles (EVs). How TAMs contribute to the EV compartment of the OC TME remains, however, underexplored. In this work peripheral blood monocytes from healthy donors were differentiated into monocyte-derived macrophages (MDMs) and polarized into classically activated (M1-like), alternatively activated (M2-like) and TAM-like (by ascites incubation). For all subtypes, serum-free conditioned medium was collected for 24 h and EVs were isolated and characterized by nano-flow cytometry (nFC), label-free mass spectrometry-based proteomics and electron microscopy, among others. Our results demonstrated distinct traits for EV release and cargo across the different macrophage subtypes. Specifically, TAM-like macrophages exhibited impaired release of small EVs and reduced frequency of tetraspanin-positive particles. These EV subpopulations displayed sizing profiles closer to M1-like than to M2-like samples. Also, the low EV release in TAM-like MDMs was accompanied by altered expression of biogenesis-related markers like flotillin-1 (FLOT1) and a decreased N-glycosylation of CD63 protein, which was validated in patient-derived samples. Remarkably, the EV-associated proteome of TAMs displayed significant enrichment in both pro- and anti-inflammatory molecules with clinical value. Markers significantly enriched in the ascites TAM-EV signature were mostly associated with poor prognosis, whereas M1-like EV-related markers (pro-inflammatory) were mostly associated with longer survival. Our results confirmed previous data for proteins like CD163 and MRC1 to be associated to TAM-EVs, while also describing novel candidates with diagnostic (i.e., COLEC12) and/or prognostic (i.e., MSR1) value in plasma. Taken together, our data support a unique secretory profile of TAMs in OC and provide new EV-associated biomarkers with translational impact. Our results pave the way for a better understanding of the mechanisms behind TAM-EV cargo loading and function, and how these cells participate in the TME landscape.

Humans↗

PubMed2026

Construction and validation of a β-hydroxybutyrylation-related molecular model for predicting prognosis of papillary thyroid carcinoma.

BACKGROUND: Papillary thyroid carcinoma (PTC) usually has a favorable prognosis, yet a subset of patients develops persistent, recurrent, or biologically aggressive disease. The clinical relevance of lysine β-hydroxybutyrylation (Kbhb)-related transcriptional programs in PTC remains unclear. Accordingly, this study aimed to characterize Kbhb-related molecular heterogeneity in PTC, construct a prognostic signature, and explore its association with the tumor microenvironment (TME). METHODS: Transcriptomic and clinical data from PTC samples within The Cancer Genome Atlas Thyroid Carcinoma (TCGA-THCA) cohort were analyzed to identify Kbhb-related differentially expressed genes (DEGs), define molecular subtypes, construct a prognostic signature, and characterize tumor microenvironmental features. Single-cell RNA-sequencing data from PTC were further used to explore the cellular distribution of representative genes. RESULTS: We identified 51 Kbhb-related DEGs in PTC and defined two Kbhb molecular subtypes. The Kbhb_C2 subtype showed shorter progression-free interval (PFI) and a more immune- and stroma-enriched microenvironment. A six-gene prognostic signature comprising TARID, CDSN, PIMREG, KLRC1, SYT13, and NPR3 was then established. High-risk patients had significantly worse PFI in the full, training, and testing cohorts, with 1-, 3-, and 5-year areas under the curve (AUCs) of 0.715, 0.793, and 0.771, respectively, in the full cohort. High-risk tumors also exhibited higher stromal, immune, and ESTIMATE scores, altered immune infiltration, and increased expression of multiple immune checkpoint molecules. Single-cell analysis confirmed distinct cell-type-specific expression patterns of representative genes. CONCLUSIONS: Kbhb-related transcriptional programs define clinically relevant molecular heterogeneity in PTC and are closely associated with prognosis and TME remodeling. The identified six-gene signature provides a biologically interpretable framework for risk stratification in PTC.

Papillary thyroid carcinoma (PTC)↗

PubMed2027

Systematic Analysis of Tumor Microenvironment Using IOBR.

The Immuno-Oncology Biological Research (IOBR) package is an R-based analysis tool for exploring the tumor microenvironment (TME) and its influence on anti-tumor immunity. Built for high-throughput data-spanning both transcriptomic and genomic profiles-IOBR integrates six analytical modules, including transcriptomic data preprocessing, TME profiling, TME pattern identification, ligand-receptor interaction analysis, genome-TME interaction assessment, and visualization. In this chapter, we walk through a multi-omics workflow using example datasets, illustrating data preparation, distribution analyses, result interpretation, and graphical output. IOBR is open source and is available at https://github.com/IOBR/IOBR and a detailed GitBook ( https://iobr.github.io/book/ ) offers a complete manual and analysis guide for each function.

Tumor Microenvironment↗

PubMed2026

The role of KIAA1467 in breast cancer: insights from pan-cancer and single-cell sequencing analysis.

BACKGROUND: Improving the response rate of single-agent immune checkpoint blockade (ICB) urgently requires the discovery of new therapeutic targets for combinatorial regimens. Analyses of tumor microenvironment (TME)-associated biomarkers have verified that KIAA1467 drives the formation of an immune-excluded, non-inflamed TME in breast cancer (BRCA). This study systematically explores the expression pattern, prognostic value, immune regulatory function, biological effects, and drug resistance relevance of FAM234B (also known as KIAA1467) in BRCA. METHODS: We performed pan-cancer survival analysis using The Cancer Genome Atlas (TCGA) datasets. Multi-omics bioinformatics analyses were conducted to evaluate KIAA1467 expression across malignancies. Single-cell RNA sequencing (scRNA-seq) data from GSE176078 was utilized to localize KIAA1467 expression at the cellular level. Immunohistochemistry and western blot assays validated KIAA1467 expression in BRCA clinical specimens. Correlation analyses were implemented to assess relationships between KIAA1467 expression, clinicopathological features, immune modulators, tumor-infiltrating immune cells, and p53 mutation status. Functional enrichment analysis uncovered relevant signaling pathways. Bioinformatic half maximal inhibitory concentration (IC50) prediction and in vitro cellular experiments were applied to evaluate associations between KIAA1467 and chemotherapeutic drug sensitivity. RESULTS: TCGA pan-cancer survival analysis demonstrated that elevated KIAA1467 expression significantly predicted shortened overall survival in BRCA and multiple other tumor types. KIAA1467 displayed distinct expression patterns across cancers, with prominent upregulation in BRCA. scRNA-seq confirmed enriched KIAA1467 expression within BRCA cells, and its upregulation in BRCA tissues was further verified by immunohistochemistry and western blot. High KIAA1467 expression was positively correlated with advanced tumor grade and lymphatic metastasis. KIAA1467 showed negative correlations with most immune modulators and core immune checkpoint molecules, as well as tumor-infiltrating immune cells in the TME, implying its potential function in tumor immune evasion. Low KIAA1467 expression was tightly linked to p53 mutations. Enrichment analysis indicated participation of KIAA1467 in epithelial-mesenchymal transition, apoptosis and cell cycle arrest. Furthermore, high KIAA1467 expression corresponded to higher estimated IC50 values of cisplatin, gefitinib, paclitaxel and gemcitabine, consistent with reduced chemosensitivity observed in vitro. CONCLUSIONS: This study reveals the multifaceted oncogenic role of KIAA1467 in BRCA. KIAA1467 participates in remodeling an immunosuppressive TME, correlates with malignant progression and chemoresistance, and may serve as a promising candidate target to optimize ICB-based combination therapy for BRCA. These findings offer new perspectives for the clinical treatment and comprehensive management of BRCA.

KIAA1467↗

PubMed2026

Targeting USP22 reprograms the tumor microenvironment and sensitizes KRAS/p53-driven lung cancer to anti-PD-1 immunotherapy.

RATIONALE: Ubiquitin-specific peptidase 22 (USP22), a deubiquitinase and component of the "Death-from-Cancer" 11-gene signature, is overexpressed in multiple malignancies and linked to recurrence, therapy resistance, and poor prognosis. Its role in KRAS/p53-driven lung cancer and the response to immune checkpoint inhibitors (ICIs) remains poorly defined. Here, we investigated USP22 as a potential therapeutic target in KRAS/p53-driven lung cancer. METHODS: A conditional Usp22 knockout (Usp22-KO) was generated in the KRASG12D; p53-/- (KP) mouse model. Cancer progression was monitored by micro-computed tomography (micro-CT). Multiplex immunofluorescence (mIF), RNA sequencing, and spatial transcriptomics profiled cancer and tumor microenvironment (TME) changes. Responses to anti-PD-1/PD-L1 therapies were compared between KP and Usp22-KO KP (KPU-) lung cancers. RESULTS: USP22 was highly expressed in early-stage KRAS/p53-driven mouse lung cancers and strongly correlated with proliferation marker Ki67. Usp22 deletion suppressed cancer growth, prolonged survival, and promoted cancer differentiation. Spatial transcriptomics and mIF revealed reduced CD206+ M2 macrophages, myeloid-derived suppressor cells (MDSCs), TGF-β1, and angiogenesis, along with increased functional CD8+ T cells. Mechanistically, USP22 regulated gene expression and protein stability, reducing c-Myc, PD-L1, TGF-β1, and SPARC upon Usp22 loss. Compared with KP cancer, KPU- and SPARC-knockdown KP cancers showed reduced macrophage chemotaxis and impaired basal- and TGF-β1-induced M2 polarization of RAW264.7 cells, suggesting that TGF-β1 and SPARC downregulation partially contributes to decreased M2 macrophage infiltration in KPU- cancers. Notably, Usp22 loss enhanced the efficacy of anti-PD-L1 and anti-PD-1 therapies in orthotopic and subcutaneous KP lung cancer models, respectively. USP22 and SPARC expression were also strongly correlated in human lung cancers. CONCLUSIONS: USP22 promotes progression and immune evasion in KRAS/p53-driven lung cancer. Targeting USP22 reprograms the TME, suppresses oncogenic signaling, and sensitizes tumors to ICI, establishing USP22 as a promising therapeutic target.

Animals↗

PubMed2025

Artificial intelligence-powered spatial analysis of tumor microenvironment in patients with non-small cell lung cancer with acquired resistance to EGFR tyrosine kinase inhibitor.

PURPOSE: This study evaluated the dynamic changes in the tumor microenvironment (TME) in patients with non-small cell lung cancer (NSCLC) and acquired resistance to epidermal growth factor receptor (EGFR)-tyrosine kinase inhibitors (TKIs) using an artificial intelligence (AI)-powered spatial TME analyzer. We then assessed the predictive efficacy of immune-checkpoint inhibitors (ICIs)-based treatment. EXPERIMENTAL DESIGN: An AI-powered whole-slide image analyzer was used to segment cancer areas (CAs) and cancer stroma and to identify tumor-infiltrating lymphocytes (TILs), tertiary lymphoid structures, fibroblasts, and endothelial cells (ECs) in the tumor tissue. We analyzed 143 NSCLC samples after resistance to EGFR-TKIs from two cohorts: (1) 89 patients treated with ICI monotherapy and (2) 54 patients from the ATTLAS phase III trial comparing atezolizumab plus bevacizumab, paclitaxel, and carboplatin (ABCP) versus pemetrexed plus carboplatin. RESULTS: Post-TKI samples showed reduced TILs in the CA (p=0.045) and increased ECs in the CA (p=0.005) compared with pre-TKI samples. These changes differed according to EGFR mutation subtype. Higher TILs in CA were associated with a better overall response rate (ORR) and progression-free survival (PFS). Similarly, higher EC levels in CA correlated with improved ORR and PFS. In the ATTLAS cohort, these factors were associated with clinical benefits from ABCP, with a significant association with TILs and a marginal association with ECs. CONCLUSION: Our findings suggest that EGFR-TKIs affect the immune landscape of patients with EGFR-mutated NSCLC. Higher TILs or ECs in the CA were significantly associated with a favorable response to subsequent ICI-based treatment. TRIAL REGISTRATION NUMBER: NCT03991403.

Aged↗

PubMed2026

Deciphering CD8+ T cell exhaustion in human cancers through single-cell and spatial transcriptomics.

Exhausted CD8+ T cells (Tex) within the tumor microenvironment (TME) represents a critical barrier limiting anti-tumor immune responses. Tex cells are characterized by upregulated inhibitory immune checkpoint receptors, reduced cytotoxicity, and functional heterogeneity. Their genomic features and regulatory networks remain poorly defined, and only a minority of patients respond to immune checkpoint blockade (ICB) therapy. Single-cell RNA sequencing (scRNA-seq), through high-resolution transcriptomic profiling, has revealed diverse Tex subpopulations, identified subpopulation-specific marker genes and regulatory pathways. Spatial transcriptomics has further mapped the spatial distribution of Tex and their interaction networks with immune cells, tumor cells, and stromal cells, elucidating the impact of spatial heterogeneity on Tex functionality. Current studies indicate that the exhausted state of Tex is dynamic and modifiable, with functional differences among subpopulations closely associated with tumor progression and therapeutic response. However, the genomic characteristics, epigenetic regulation, and spatial interaction mechanisms of Tex require further exploration. This review summarizes recent advances in high-resolution omics technologies for precisely dissecting Tex heterogeneity, functional features, and interactions with other cells. It emphasizes the central value of optimizing Tex-targeted tumor immunotherapy strategies, providing theoretical foundations and directional guidance for developing more effective anti-tumor immunotherapies.

Humans↗

PubMed2025

Histology-Based Virtual RNA Inference Identifies Pathways Associated With Metastasis Risk in Colorectal Cancer.

Colorectal cancer (CRC) remains a major health concern, with >150,000 new diagnoses and >50,000 deaths annually in the United States, underscoring an urgent need for improved screening, prognostication, disease management, and therapeutic approaches. The tumor microenvironment (TME)-comprising cancerous and immune cells interacting within the tumor's spatial architecture-plays a critical role in disease progression and treatment outcomes, reinforcing its importance as a prognostic marker for metastasis and recurrence risk. However, traditional methods for TME characterization, such as bulk transcriptomics and multiplex protein assays, lack sufficient spatial resolution. Although spatial transcriptomics (ST) allows for the high-resolution mapping of whole transcriptomes at near-cellular resolution, current ST technologies (eg, Visium and Xenium) are limited by high costs, low throughput, and issues with reproducibility, preventing their widespread application in large-scale molecular epidemiology studies. In this study, we refined and implemented virtual RNA inference (VRI) to derive ST-level molecular information directly from hematoxylin and eosin (H&E)-stained tissue images. Our VRI models were trained on the largest matched CRC ST data set to date, comprising 45 patients and >300,000 Visium spots from primary tumors. Using state-of-the-art deep learning models (UNI, ResNet-50, Vision Transformer, and Vision Mamba), we achieved a median Spearman's correlation coefficient of 0.546 between predicted and measured spot-level expression. As validation, VRI-derived gene signatures linked to specific tissue regions (tumor, interface, submucosa, stroma, serosa, muscularis, and inflammation) showed strong concordance with signatures generated via direct ST, and VRI performed accurately in estimating cell-type proportions spatially from H&E slides. In an expanded CRC cohort controlling for tumor invasiveness and clinical factors, we further identified VRI-derived gene signatures significantly associated with key prognostic outcomes, including metastasis status. Although certain tumor-related pathways are not fully captured by histology alone, our findings highlight the ability of VRI to infer a wide range of "histology-associated" biological pathways at near-cellular resolution without requiring ST profiling. Future efforts will extend this framework to expand TME phenotyping from standard H&E tissue images, with the potential to accelerate translational CRC research at scale.

Humans↗

PubMed2026

A novel lactylation-related gene signature deciphers the immunosuppressive microenvironment and stratifies precision therapy in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) remains a leading cause of cancer mortality, largely due to the heterogeneity of the tumor microenvironment (TME) and the limited efficacy of immunotherapy in microsatellite stable (MSS) tumors. Histone lactylation, a post-translational modification derived from the Warburg effect, serves as a critical bridge linking metabolic reprogramming to gene regulation and immune evasion; however, its specific prognostic value and clinical implications in CRC remain to be fully elucidated. METHODS: In this study, we systematically analyzed transcriptome profiling data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts, supplemented by single-cell RNA sequencing (scRNA-seq) analysis and Human Protein Atlas (HPA) protein-level validation. By integrating univariate Cox regression, Least Absolute Shrinkage and Selection Operator (LASSO) analysis, and multivariate Cox regression, we constructed a novel lactylation-related gene (LRG) risk signature. We extensively evaluated the association between this risk signature and patient prognosis, immune infiltration patterns, somatic mutations, and therapeutic sensitivity. RESULTS: A robust 9-gene prognostic signature (DHRS7, SPR, MBD2, RBM17, CSRP2, S100A4, TMSB4X, TKT, COPS4) was identified and corroborated at the protein level. Patients with high risk scores exhibited significantly worse overall survival (OS) across the training and two independent validation cohorts. Immunogenomic and scRNA-seq analyses revealed that high-risk tumors were characterized by an immunosuppressive and stromal-dense microenvironment-with stromal cells exhibiting the highest lactylation risk scores-enriched with regulatory T cells (Tregs), and frequently harbored PIK3CA mutations. Differential expression analysis indicated that this immune exclusion is structurally maintained by enriched extracellular matrix (ECM) organization and TGF-β signaling. Conversely, low-risk tumors displayed an inflamed phenotype with active antitumor immunity. Pharmacogenomic prediction identified distinct therapeutic stratifications: low-risk patients exhibited significant sensitivity to standard chemotherapeutics (fluorouracil, oxaliplatin) and EGFR/HER2 inhibitors (e.g., lapatinib, erlotinib). In contrast, high-risk patients showed specific vulnerabilities to novel targeted agents, including PI3K pathway inhibitors (TG-100-115, XL765), microenvironment-modulating agents (sildenafil, GANT-61), and epigenetic inhibitors (UNC0638). CONCLUSION: We established a novel lactylation-related risk signature that effectively stratifies CRC patients by prognosis and TME characteristics. By elucidating the crosstalk between metabolic dysregulation, stromal barriers, and immune exclusion, this study provides potential biomarkers and stratified therapeutic strategies-ranging from standard chemotherapy to targeted metabolic and stromal interventions-to optimize precision medicine for CRC patients.

Colorectal cancer↗

PubMed2026

Malignant epithelial states drive immune dysfunction in ampulla of Vater carcinoma.

BACKGROUND: Ampulla of Vater (AoV) carcinoma is a rare malignancy arising at the junction of intestinal and pancreatobiliary epithelium. Its heterogeneous clinical behavior and histological diversity have hindered therapeutic advances, and the cellular basis of this heterogeneity remains unclear. We aimed to construct a single-cell transcriptomic atlas of AoV carcinoma, with a focus on identifying epithelial subtypes and their interactions with the tumor microenvironment (TME). METHODS: We performed single-cell RNA sequencing on eight primary AoV tumors and four matched normal tissues. Comprehensive clustering and transcriptomic analyses identified cell-type composition, epithelial heterogeneity, and tumor-immune interactions. Findings were validated using deconvolution of bulk RNA-seq data from 62 AoV carcinoma patients. Results Malignant epithelial cells were categorized into four distinct subtypes: Int-Wnt, PB-KRAS, Int-Hypoxia, and Cycling stage. PB-KRAS cells exhibited stem-like transcriptional programs and high genomic instability. Deconvolution analysis of bulk RNA-seq data from the independent AoV cohort revealed that enrichment of the PB-KRAS subtype correlated with tumor recurrence and poor survival. Our immune profiling analysis discovered a significant association between PB-KRAS subtype and GZMK+ CD8+ T cells, which are in a pre-dysfunctional state, alongside SPP1+ macrophages exhibiting immunosuppressive traits. Spatial transcriptome data further supports the immunosuppressive natures of TME around PB-KRAS subtype malignant epithelial cells in AoV carcinoma. CONCLUSIONS: Our study presents a single-cell atlas of AoV carcinoma, highlighting the molecular diversity of malignant epithelium and its association with the immune microenvironment. The PB-KRAS subtype emerges as a stem-like, immunosuppressive tumor state associated with poor prognosis, providing insights for future therapeutic targeting.

Ampulla of Vater carcinoma↗

PubMed2026

A pan-cancer single-cell atlas uncovers the role of sex hormones and chromosomes in sex-divergent reprogramming of the tumor microenvironment.

BACKGROUND: Sex bias is pervasive in tumors; however, how sex chromosomes and hormone-responsive signaling shape the tumor microenvironment (TME) remains insufficiently characterized. Considering the critical impact of the TME on tumor progression and response to immunotherapy, a pan-cancer investigation of sex-specific and cancer-context-dependent TME features is warranted. METHOD: Based on stringent inclusion criteria, we constructed a high-resolution pan-cancer single-cell sequencing atlas by integrating 31 publicly available single-cell RNA-seq datasets, comprising a total of 1,831,436 cells by integrating 468 samples from eight types of non-sex-specific solid tumors (282 males and 186 females). After correcting for batch effects, we identified major and minor cellular subsets. Multiple computational approaches were applied to investigate sex-associated differences in cellular composition, gene expression, pathway activity, malignant cell states and intercellular communication. RESULTS: We systematically compared sex-specific TME features across eight common solid malignancies. Male-biased CD8+ T cell exhaustion emerged as a recurrent but non-uniform feature, with its magnitude varying across cancer types and being modified by tissue-specific contexts. This pattern was associated with androgen-response signature scores and expression-based loss of the Y chromosome (LOY) scores. M2-like macrophage polarization showed a more cancer-type-dependent pattern; although female-biased enrichment was observed in selected malignancies, it did not represent a uniform pan-cancer feature. Expression-based X chromosome inactivation (XCI)/XCI escape-related programs, estrogen-response signature scores and stromal components, including fibroblasts and endothelial cells, were associated with macrophage and immune-regulatory states in specific tumor contexts. Tumor cells of male origin displayed higher genomic instability and more aggressive phenotypes, with androgen-response signatures and LOY contributing to the development of a male biased malignant state. Furthermore, expression-based LOY scores in malignant cells were associated with CD8+ T cell exhaustion based on transcriptomic proxies. CONCLUSION: Our study uncovers extensive but heterogeneous sex-specific differences in the TME across multiple cancer types. We propose a regulatory framework linking sex chromosomes, hormone-responsive signaling and TME interactions, which is consistent with recurrent male-biased CD8⁺ T cell exhaustion and context-dependent M2-like macrophage polarization. Importantly, the magnitude and, in some cancers, the direction of these sex-biased features are modified by tissue-specific contexts. These findings underscore the need to include sex chromosome and hormone status as essential biological variables in studies of the tumor microenvironment and the design of immunotherapies.

Tumor Microenvironment↗