SEARCH · Search PubMed

Results for “Immune cell subtypes”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 19 recordsLinked to original sources

PubMed2026

Identification of molecular subtypes in clear cell renal cell carcinoma based on chromatin regulators and tumor immune microenvironment profiling.

In the histological classification of renal cell carcinoma, clear cell renal cell carcinoma (ccRCC) accounts for the highest proportion and is the most common subtype. Despite advances in management, it continues to be associated with considerable incidence and mortality. Although surgery and systemic therapies are available, their efficacy is constrained by pronounced intratumoral heterogeneity and treatment resistance. Identifying robust biomarkers and clarifying the underlying biological mechanisms are therefore essential to improving diagnosis, risk stratification and therapeutic decision-making. In this work, we identified two ccRCC molecular subtypes displaying divergent chromatin regulator (CR) profiles and different clinical prognoses. Using the genes differentially expressed between these subgroups, we constructed a CR-related score (CRS) that effectively stratified patients according to survival. More analysis concluded that the low expression of CR was more linked with the immune-activated tumors, which encompassed the immune pathway enrichment, as well as the elevation of numerous immune cell subtypes. Moreover, elevated CRS was associated with improved immunotherapy responsiveness. Drug-sensitivity analyses nominated several candidate agents, and SMARCD3 knockdown in 786-O cells inhibited proliferation and migration and reduced sensitivity to masitinib. Collectively, these findings support the prognostic and therapeutic relevance of CR-related states in ccRCC and provide a framework for future experimental validation of chromatin-regulated tumor-immune interactions.

Humans↗

PubMed2026

Rethinking immune studies: population-level immune variations and the path forward.

Shifts in immune cell proportions underlie disease progression and immunotherapy response, positioning them as promising diagnostic biomarkers and therapeutic targets. However, these shifts also occur naturally across the lifespan and vary with demographic factors such as age and sex, which may complicate their interpretation and clinical utility. While demographic associations have been explored previously, there have been mixed results likely due to small sample sizes and cross-cohort population-specific differences. To address these limitations, we conducted a meta-analysis across 3 large, diverse cohort studies to evaluate associations between 20 immune cell subtypes, 3 informative cell ratios, and a range of sociodemographic variables such as age, sex, self-identified race and ethnicity (SIRE), and socioeconomic status. We find consistent and significant associations across all sociodemographic dimensions. Cytomegalovirus (CMV)-a key driver of immune senescence-emerged as a major contributor to variation in immune composition and CMV antibody levels were higher among women, individuals of lower socioeconomic status, and marginalized racial and ethnic groups. In addition, male sex showed similar patterns of association with immune profiles as aging, whereas race did not. These findings underscore the need to account for diverse sociodemographic factors in immunology study design and participant recruitment to avoid population-specific biases and ensure broadly generalizable results.

Humans↗

PubMed2026

Mapping the Immune cell-specific gene regulatory network in bipolar disorder: A framework from scTWMR to exploratory drug-target annotation.

BACKGROUND: Although the involvement of the immune system in the genetic susceptibility of bipolar disorder (BD) is widely acknowledged, the causal relationship between gene expression in specific immune cell subtypes and BD requires systematic elucidation. METHODS: We implemented an analytical framework integrating single-cell transcriptome-wide Mendelian randomization (scTWMR) with colocalization analysis. This approach utilized cis-expression quantitative trait loci (cis-eQTLs) derived from 14 distinct immune cell types as instrumental variables to interrogate BD genome-wide association study (GWAS) summary statistics (comprising 41,917 cases and 371,549 controls). Subsequent investigations encompassed functional enrichment analysis, protein-protein interaction (PPI) network construction, phenome-wide association study (PheWAS), and performed an exploratory drug-target annotation. RESULTS: Our analysis identified 33 gene-immune cell associations. Colocalization analysis provided robust evidence (PPH4 > 90%) for shared causal variants implicating the MAD1L1, APOM, and NFKBIL1 loci. Significantly enriched biological pathways included cell cycle regulation, circadian rhythm entrainment, and neuroinflammation. The PPI network revealed a core regulatory module centered on histone-encoding and immune-related genes. Exploratory drug-target annotation nominated compounds for further investigation for compounds targeting APOM, TMEM258, and NFKBIL1. CONCLUSION: This study systematically delineates a genetically supported regulatory network of immune cell-specific gene expression in BD, predominantly implicating CD8⁺ effector T cells, plasma cells, and B cells. The findings corroborate established pathological pathways while uncovering novel cell type-specific therapeutic targets, thereby providing a genetic framework for prioritizing candidate targets for future investigation.

Bipolar disorder↗

PubMed2026

Transcriptomic profiling across stages of non-muscle-invasive bladder cancer identifies fibroblast activation protein-alpha as a stromal biomarker associated with progression.

BACKGROUND: T1 non-muscle-invasive bladder cancer (NMIBC) represents a biologically aggressive subgroup with substantial heterogeneity in recurrence and progression risk. Current clinicopathological risk stratification tools lack sufficient precision to identify patients at the highest risk of progression to muscle-invasive bladder cancer (MIBC). OBJECTIVE: To characterize transcriptomic differences between T1 and < T1 (Ta/Tis) NMIBC and to explore the association of fibroblast activation protein-α (FAP) gene expression with disease progression. METHODS: Transcriptomic profiling was performed on formalin-fixed paraffin-embedded (FFPE) tumor tissue from 66 patients with primary, treatment-naïve NMIBC and 5 patients with T2 disease (included for exploratory comparisons). Analyses included differential gene expression, gene set enrichment analysis (GSEA), molecular subtyping, immune cell deconvolution, and evaluation of FAP expression in relation to recurrence and progression. External validation of FAP was conducted in three independent NMIBC cohorts. RESULTS: T1 tumors demonstrated a distinct transcriptomic profile compared with < T1 tumors, characterized by enrichment of cell cycle-related and metabolic pathways and a higher prevalence of aggressive molecular subtypes. Despite these molecular differences, no statistically significant differences in recurrence-free, progression-free, cancer-specific, and overall survival were observed, likely reflecting limited event numbers. Among recurrent tumors, early recurrences (≤ 24 months) were associated with epithelial-mesenchymal transition signatures. FAP expression increased with tumor stage (p = 0.0005) and was associated with progression (p = 0.002) and mortality (p = 0.01). Patients with tumors in the highest quartile of FAP expression had worse progression-free survival. This association was consistently observed in three external NMIBC cohorts. CONCLUSIONS: T1 NMIBC exhibits distinct transcriptomic features suggestive of increased biological aggressiveness. Elevated FAP expression is reproducibly associated with progression risk across multiple cohorts, supporting its potential role as a biomarker of aggressive disease. Given the limited number of progression events, these findings should be considered hypothesis-generating and warrant prospective validation before clinical implementation.

Humans↗

PubMed2026

Cell Type-Resolved Causal Inference and Spatial Transcriptomic Integration Reveal Immune-Specific Genetic Drivers of Autoimmune and Malignant Thyroid Disease.

BACKGROUND: Thyroid diseases, including autoimmune thyroid disease (AITD) and thyroid cancer, are characterized by immune dysregulation, yet the cell type-specific genetic mechanisms underlying these conditions remain poorly understood. Most genome-wide association studies (GWAS) have relied on bulk tissue expression quantitative trait loci (eQTL), which cannot resolve the heterogeneity of immune cell populations. METHODS: We performed two-sample Mendelian randomization (MR) analyses using single-cell cis-eQTLs from 14 immune cell subtypes (OneK1K cohort) as instrumental variables against GWAS summary statistics for four thyroid outcomes: autoimmune hyperthyroidism, autoimmune hypothyroidism, thyroid cancer and autoimmune thyroiditis. Causal associations were validated through Bayesian colocalization, phenome-wide association analysis (PheWAS) and multi-layered transcriptomic validation encompassing spatial transcriptomics of AITD tissue (GSE248205), bulk RNA-seq of thyroid cancer (GSE3678) and single-cell RNA-seq of thyroid tumours (GSE250521). gsMap spatial LD score regression was applied to map disease heritability onto spatial tissue architecture. RESULTS: We identified six Bonferroni-significant causal gene-cell type pairs for autoimmune hyperthyroidism, including protective effects of ABHD16A in naïve/immature B cells (OR = 0.440), HIST1H3H in CD8 NC T cells (OR = 0.324), HMGN4 in NK recruiting cells (OR = 0.556) and ZKSCAN4 in CD8 S100B T cells (OR = 0.427), with five pairs showing strong colocalization (PP.H4 ≥ 86%). Three pairs reached significance for autoimmune hypothyroidism, including a risk association of HLA-F in CD4 NC T cells (OR = 1.139). For autoimmune thyroiditis, FAM134B/RETREG1 showed consistent suggestive protective associations across both CD4 and CD8 NC T cells (PP.H4 ≥ 90% for both), suggesting a possible involvement of ER phagy regulation in thyroiditis susceptibility. Thyroid cancer showed a suggestive association with HLA-G in classical monocytes (OR = 1.899, PP.H4 = 53%). Spatial transcriptomic validation demonstrated progressive immune infiltration from control tissue to Graves' disease to Hashimoto's thyroiditis (7.7%-15.7%, 46.1%-54.1%, respectively) and strong spatial correlation between target gene expression and corresponding cell type enrichment (e.g., plasma cell-HLA-DQB1: r = 0.491, p < 10-300). HLA-G was independently validated in thyroid cancer bulk (log2fc = 0.542, p = 9.51 × 10-3, AUC = 0.857) and single-cell datasets. PheWAS revealed no significant associations detected for the core candidates. gsMap identified significant enrichment of autoimmune hypothyroidism heritability in gastrointestinal tract, adrenal gland and adipose tissue (all Bonferroni p < 0.002). CONCLUSIONS: This study establishes a multi-scale analytical framework integrating cell type-resolved genetic inference with spatial tissue validation, revealing distinct immunogenetic architectures underlying autoimmune versus malignant thyroid disease. Protective genetic programs in autoimmune hyperthyroidism converge on chromatin remodelling (HIST1H3H, HMGN4, ZKSCAN4) and lipid metabolism (ABHD16A) across lymphocyte subsets, whereas thyroid cancer risk involves immune escape mediated by HLA-G in myeloid cells. The ER-phagy receptor RETREG1 represents a candidate pathway warranting further investigation in autoimmune thyroiditis. These findings provide genetically supported, cell type-specific therapeutic targets and demonstrate a generalizable strategy for dissecting the immune-mediated mechanisms of complex thyroid diseases.

Mendelian randomization↗

PubMed2026

Nested co-expression network analysis identifies compact gene clusters in a black box.

MOTIVATION: Digital analysis of biological systems requires methods capable of identifying both broad and nested gene modules reflecting complex biological processes. Existing transcriptomic methods often miss compact gene sets corresponding to subprocesses in specialized cell types, limiting insights into functional heterogeneity. RESULTS: We present Nested-WGCNA, a two-stage unsupervised network analysis algorithm designed to identify coarse-grained and fine-grained gene modules. Applied to bulk RNA-Seq data, Nested-WGCNA reveals stable modules reproducible across datasets. When validated against scRNA-Seq data, these modules correspond to both major and minor immune cell subtypes. Application to immunotherapy response datasets uncovers predictive and prognostic biomarkers, highlighting its utility in treatment stratification and biomarker discovery. AVAILABILITY: The NestedWGCNA source code and analysis pipeline are available on GitHub (https://github.com/ilyada/NestedWGCNA) and archived on Zenodo (https://doi.org/10.5281/zenodo.18959244).

Algorithms↗

PubMed2026

Plasma cell-CD8+ T cell co-enrichment distinguishes immunotherapy-responsive hepatocellular carcinoma subtypes.

BACKGROUND: Hepatocellular carcinoma (HCC) is characterised by significant racial disparities in incidence and outcomes, yet whether these reflect distinct tumour biology or differential distribution of molecular subtypes among immunotherapy patients remains unclear. METHODS: We characterised molecular heterogeneity among 46 patients with HCC of differing background population from the NCI-CLARITY cohort receiving immune checkpoint inhibitor therapy, using transcriptomic and genomic profiling, with validation across multiple independent cohorts. RESULTS: Differential expression analysis comparing African American versus non-African American patients identified 126 genes, of which 55 demonstrated tumour-specific expression across independent validation cohorts with paired tumour-normal samples. Consensus clustering revealed two molecular subtypes with no significant race association, indicating these clusters capture tumour-intrinsic biology rather than ancestry. The genomic landscape showed minimal differences between subtypes. A prognostic signature derived from these expression profiles demonstrated significant risk stratification in the NCI-CLARITY cohort and TCGA-LIHC, but not in Asian cohorts, suggesting population-specific applicability. Immune deconvolution revealed that the two subtypes represent distinct immune microenvironments: one subtype exhibited markedly elevated plasma cell infiltration with strong plasma cell-CD8+T cell correlation suggesting coordinated adaptive immunity, along with elevated tertiary lymphoid structure signatures. The other subtype showed regulatory T cell-macrophage correlation and enrichment for immune-excluded phenotypes. The immune-enriched subtype trended towards higher immunotherapy response rates. CONCLUSIONS: Molecular heterogeneity in HCC reveals distinct tumour-immune ecosystems that transcend racial classification. Tumour immune heterogeneity in HCC reflects distinct molecular patterns, with immune hot tumours characterised by elevated tertiary lymphoid structure signatures and enriched plasma cell and CD8+T cells. These patterns may serve as prognostic biomarkers for immunotherapy patient stratification and demonstrate the value of diverse cohort representation in identifying clinically relevant therapeutic targets.

Gastrointestinal Cancer↗

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↗

PubMed2025

Integration of multiple omics reveals key targets and cellular mechanisms for intervention in sarcopenia.

BACKGROUND: Sarcopenia, an age-related syndrome characterized by progressive loss of muscle mass, strength, and function, presents a significant global health burden with limited therapeutic interventions. This study integrates genomic causality, multi-tissue omics, and cellular mediation analyses to identify and prioritize mechanistically grounded therapeutic targets. METHODS: A multi-tiered analytical framework was applied, beginning with two-sample Mendelian randomization (MR) to infer causal relationships between 4907 plasma proteins (cis-pQTLs from 35,559 individuals) and sarcopenia traits in Pan-UK Biobank participants. Bayesian colocalization and transcriptomic validation in human sarcopenia muscle biopsies were employed to prioritize targets. Cellular mediation analysis quantified contributions of immune and stromal cell subtypes to protein-trait pathways using transcriptomic deconvolution. RESULTS: MR identified 1237 plasma proteins causally associated with sarcopenia traits, with six targets (HGFAC, GATM, HMOX2, F2, LMAN2L, HPGDS) validated through colocalization, transcriptomic expression, and sarcopenia-related dysregulation. Cellular mediation revealed immune mechanisms underlying HGFAC's effects, with CD4+ regulatory T cells mediating 3.49 % of its impact on sarcopenia traits. Prothrombin exhibited muscle-protective effects independent of coagulation. CONCLUSION: This study establishes a causal map linking plasma proteins to sarcopenia through immune-stromal interactions. The integration of MR, multi-omics validation, and cellular mediation prioritizes six proteins as actionable targets, supporting repurposing of thrombin inhibitors and development of immunometabolic therapies. The framework bridges genomic causality with cellular pathophysiology, advancing precision strategies for age-related muscle decline.

Humans↗

PubMed2025

Tumor-immune partitioning and clustering algorithm for identifying tumor-immune cell spatial interaction signatures within the tumor microenvironment.

BACKGROUND: Growing evidence supports the importance of characterizing the organizational patterns of various cellular constituents in the tumor microenvironment in precision oncology. Most existing data on immune cell infiltrates in tumors, which are based on immune cell counts or nearest neighbor-type analyses, have failed to fully capture the cellular organization and heterogeneity. METHODS: We introduce a computational algorithm, termed Tumor-Immune Partitioning and Clustering (TIPC), that jointly measures immune cell partitioning between tumor epithelial and stromal areas and immune cell clustering versus dispersion. As proof-of-principle, we applied TIPC to a prospective cohort incident tumor biobank containing 931 colorectal carcinoma cases. TIPC identified tumor subtypes with unique spatial patterns between tumor cells and T lymphocytes linked to certain molecular pathologic and prognostic features. T lymphocyte identification and phenotyping were achieved using multiplexed (multispectral) immunofluorescence. In a separate hepatocellular carcinoma cohort, we replaced the stromal component with specific immune cell types-CXCR3+CD68+ or CD8+-to profile their spatial relationships with CXCL9+CD68+ cells. RESULTS: Six unsupervised TIPC subtypes based on T lymphocyte distribution patterns were identified, comprising two cold and four hot subtypes. Three of the four hot subtypes were associated with significantly longer colorectal cancer (CRC)-specific survival compared to a reference cold subtype. Our analysis showed that variations in T-cell densities among the TIPC subtypes did not strictly correlate with prognostic benefits, underscoring the prognostic significance of immune cell spatial patterns. Additionally, TIPC revealed two spatially distinct and cell density-specific subtypes among microsatellite instability-high colorectal cancers, indicating its potential to upgrade tumor subtyping. TIPC was also applied to additional immune cell types, eosinophils and neutrophils, identified using morphology and supervised machine learning; here two tumor subtypes with similarly low densities, namely 'cold, tumor-rich' and 'cold, stroma-rich', exhibited differential prognostic associations. Lastly, we validated our methods and results using The Cancer Genome Atlas colon and rectal adenocarcinoma data (n = 570). Moreover, applying TIPC to hepatocellular carcinoma cases (n = 27) highlighted critical cell interactions like CXCL9-CXCR3 and CXCL9-CD8. CONCLUSIONS: Unsupervised discoveries of microgeometric tissue organizational patterns and novel tumor subtypes using the TIPC algorithm can deepen our understanding of the tumor immune microenvironment and likely inform precision cancer immunotherapy.

Humans↗

PubMed2026

Single-cell profiling reveals a novel CAF subpopulation linking stromal heterogeneity to immune suppression in breast cancer subtypes.

BACKGROUND: The tumor microenvironment critically influences breast cancer (BC) progression, immune surveillance, and therapeutic response. Cancer-associated fibroblasts (CAFs), a heterogeneous stromal population, are key regulators of these processes, yet their subtype-specific contributions in BC remain insufficiently defined. METHODS: We integrated three single-cell RNA sequencing datasets from 29 BC patients to characterize stromal populations. Bulk RNA-seq data from The Cancer Genome Atlas (TCGA) were analyzed to assess correlations between CAF subsets and immune infiltration. Gene signatures were derived to identify subtype-specific CAF-immune interactions, prognostic markers, and potential predictors of chemotherapy response. RESULTS: Three conserved stromal populations (iCAFs, myCAFs, and pericytes) were identified, along with a previously unrecognized subset, the cluster 3 (CL3) CAF-like cells, referred as metabolic stressed CAF (msCAF). msCAF cells displayed transcriptional programs associated with antigen presentation, stress response, glycolysis, and extracellular matrix remodeling. Their abundance was inversely correlated with T-cell infiltration and function, in a subtype-specific manner: triple negative breast cancer (TNBC) was enriched for msCAFs in immune-infiltrated but functionally constrained microenvironments, whereas Luminal A tumors exhibited weaker immune infiltration with heterogeneous CAF-immune associations. msCAFs were characterized by a conserved gene signature (HLA-A, HLA-C, IL32, EMP3) and subtype-specific genes related to T-cell exhaustion. Several genes demonstrated prognostic relevance with distinct patterns in Luminal A (IER3, TIMP1, TBX3, SEC61G) and TNBC (ADM, C4orf3, LDHA) tumors, as well as shared biomarkers (FN1, LOXL2, P4HA1). Multiple msCAF genes also predicted chemotherapy response, suggesting utility as treatment stratification biomarkers. CONCLUSION: msCAFs represent a clinically relevant CAF subset that drives immune suppression, impacts subtype-specific prognosis, and influences therapy response in BC. These findings highlight msCAFs as promising targets for enhancing immunotherapy and personalizing treatment strategies.

Humans↗

PubMed2026

Immune subtyping of colorectal adenoma identifies a subtype with activated adaptive immunity ahead of progressing to cancer.

BACKGROUND: Colorectal adenomas (CRA) represent precursor lesions with varying risks of malignant transformation. However, molecular subtyping, particularly immune-related classification, remains underexplored in adenomas. This study aims to characterize the immune landscape of CRA through immune subtyping and evaluate its association with cancer progression, gene expression signatures, and functional pathways. METHODS: We conducted a retrospective analysis of transcriptomic data from multiple cohorts of CRA samples. Immune subtypes were identified using non-negative matrix factorization (NMF) based on immune-related genes. Diverse deconvolution algorithms were used to estimate immune cell infiltration. The immune status alteration in premalignant lesion was further consolidated by single-cell transcriptome data. Differential gene expression analysis was performed between subtypes, followed by functional enrichment analyses (Gene Ontology [GO] and Kyoto Encyclopedia of Genes and Genomes [KEGG]). RESULTS: Two distinct immune subtypes were identified: an immune-enriched subtype characterized by high lymphocyte infiltration and elevated expression of immune-related genes, and an immune-deficient subtype with suppressed immune activity. Differential expression analysis revealed significant upregulation of immune response genes (e.g., CD4, CD86, HLA-DRA) in the immune-enriched subtype. GO and KEGG analyses highlighted enrichments in leukocyte transendothelial migration, chemokine signaling, and antigen processing and presentation pathways. Single-cell result revealed an early occurrence of TIGIT activation and exhausted CD8 T cell features in adenoma when compared to normal tissue. CONCLUSION: This study delineates distinct immune subtypes within CRAs. The immune-enriched subtype demonstrates activated adaptive immunity and may reflect a higher potential for immune surveillance, while the immune-deficient subtype exhibits stromal features suggestive of progressive transformation. These findings provide insights into early immune microenvironment alterations and may inform strategies for risk stratification and immunoprevention in colorectal carcinogenesis.

Colorectal adenoma↗

PubMed2026

Profiling tumor immune microenvironment of epithelial ovarian carcinoma.

BACKGROUND: Epithelial ovarian carcinoma (EOC) comprises five main histological subtypes: high-grade serous (HGSOC), low-grade serous (LGSOC), clear cell (CCOC), mucinous (MOC), and endometrioid (ENOC). Each histotype harbors specific genomic alterations and clinical outcome. Few studies systematically compared the tumor immune microenvironment across the five subtypes. METHODS: We performed 7-plex (CD45, CD8, CD68, CD163, FoxP3, CD20, and cytokeratin) sequential immunohistochemistry on a clinically annotated tissue microarray including 139 EOC representing the five subtypes and 26 borderline tumors (serous and mucinous). Digital pathology was used to quantify immune cell abundance, their spatial distribution (stroma vs tumor core), and correlation with survival. RESULTS: Immune cells were dominated by macrophages and more abundant in the stroma than tumor core across the five subtypes, consistent with immune excluded pattern. Compared to HGSOC, CCOC displayed the highest infiltration by CD45+ leukocytes and CD68+ macrophages, particularly M2-like CD163+ cells, suggesting a macrophage-rich, immunosuppressive phenotype. LGSOC exhibited the highest infiltration by intraepithelial FoxP3+ regulatory T cells. Comparison of borderline tumors with invasive carcinoma (LGOSC and MOC) revealed that malignant progression is accompanied by loss of CD8+ T cells, enrichment in regulatory T cells and increase of CD163+/CD68+ ratio, consistent with immune evasion during tumorigenesis. There was a trend toward better survival in HGSOC highly infiltrated by lymphocytes, either intraepithelial (CD8+ and FoxP3+) or stromal (FoxP3+ and CD20+). CONCLUSIONS: EOC is characterized by histotype-specific immune milieux defined by macrophage dominance, epithelial immune exclusion and dynamic immune remodeling during progression from borderline tumors to invasive carcinomas.

Humans↗

PubMed2025

m6A regulator-based molecular classification and hub genes associated with immune infiltration characteristics and clinical outcomes in diffuse gliomas.

BACKGROUND: m6A methylation modification is a new regulatory mechanism involved in tumorigenesis and tumor-immunity interaction. However, its impact on glioma immune microenvironment and clinical outcomes remains unclear. METHODS: Comprehensive expression profiles of 18 m6A regulators were used to identify molecular subtypes exhibiting distinct m6A modification patterns in 1673 glioma samples sourced from public datasets. A multi-genes signature was constructed for predicting clinical outcomes and response to immunotherapy in glioma patients. Immunohistochemistry and cellular experiments were performed for validation. RESULTS: Two m6A subtypes of gliomas were identified. The m6A-low-risk subtype was characterized by paucity of immune infiltrates; While the m6A-high-risk subtype had higher abundances of multiple immune cells including lymphocyte and macrophage as well as increased expression of PD-L1, corresponding to an immunosuppressive phenotype. The m6A-high-risk subtype had poorer survival than the m6A-low-risk subtype in both the glioblastoma and lower grade gliomas cohorts. Eight m6A-related hub genes of high prognostic significances were identified and selected for developing a scoring signature termed as m6Ascore. Elevated m6Ascore indicated worse survival for glioma patients under standard care, but showed enhanced response to immunotherapy. Moreover, we demonstrated that overexpression of FTO, a m6A demethylase, inhibited the expressions of m6A-related hub genes (PTX3, SPAG4), impaired glioma cell viability and reduced macrophage chemotaxis. CONCLUSION: This work develops an immune- and clinical-relevant m6A subtyping and a scoring model, which enhances our understanding of the role of m6A modification in regulating immune infiltration microenvironment in gliomas and helps to identify patients who are more likely to benefit from immunotherapy.

Humans↗

PubMed2026

Molecular subtyping of adrenocortical carcinoma reveals distinct subtypes with prognostic and therapeutic implications.

Adrenocortical carcinoma (ACC) is a rare but aggressive malignancy with poor survival and limited treatment options. To comprehensively characterize its molecular landscape and identify clinically relevant subtypes, we performed an integrated genomic analysis - including whole-exome sequencing, RNA sequencing, and copy number variation profiling - on 61 Chinese patients with ACC. We identified recurrent mutations in TP53 (25%), CTNNB1 (15%), ZNRF3 (10%), and MEN1 (8%). Unsupervised clustering of transcriptomic data revealed four distinct molecular subtypes: cortisol-driven (CD, 14%), immune-suppressed (IS, 40%), cell cycle-altered (CCA, 22%), and immunomodulatory (IM, 24%). The CD subtype exhibited steroidogenic pathway activation; the IS subtype showed T cell receptor downregulation and the worst disease-free survival; the CCA subtype was marked by chromosomal instability and cell cycle gene overexpression; and the IM subtype displayed enriched immune signaling and favorable outcomes. Copy number analysis further uncovered focal amplifications (e.g. TERT, CDK4) and HLA-II deletions. This study establishes a novel molecular classification of ACC, providing a framework for subtype-specific therapeutic strategies, such as CDK4/6 inhibition for CCA and immunotherapy for IM tumors, while highlighting the clinical challenges of immune-cold IS tumors.

Humans↗

PubMed2026

Integrative evidence-knowledge marker selection enhances LLM-based cell type annotation in single-cell RNA-seq analysis.

BACKGROUND: Cell type annotation is essential for gaining biological insight from single-cell RNA sequencing data, yet manual labeling remains time-consuming and difficult to reproduce. Various computational approaches have been developed to automate this process, and recent studies suggest that large language models can infer cell types with promising accuracy in single-cell analysis. However, most workflows still rely on cluster-specific markers derived from gene expression alone or manual curation. As a result, marker selection can be sensitive to statistical criteria and dataset-dependent bias, which may lead to the selection of less informative genes or missing important markers, while providing limited biological context. RESULTS: To address this limitation, we introduce CELLIA, an LLM-based workflow for automated and robust cell type annotation. CELLIA employs an integrative evidence-knowledge marker selection strategy that combines statistical differential expression criteria with curated tissue-specific marker resources to identify informative marker genes. In benchmarking analyses of 102 cell types, this approach improved agreement with manual annotations. In addition, CELLIA achieved higher agreement in subtype-level analyses of closely related immune populations and was further evaluated in a non-immune stromal subtype setting, covering 25 cell types in total. CONCLUSION: By integrating evidence-knowledge from gene expression with curated biological prior knowledge, CELLIA provides a more stable marker selection and improves the reliability of LLM-cell type annotation.

Cell type annotation↗

PubMed2026

Integrative subtyping by bile acid metabolism identifies CLCA1/UGT2A3/ZG16 as markers of immune dysfunction and poor prognosis in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is the primary driver of cancer-related death and illness across the world. Despite the full-scale shift of the treatment approach for some colorectal cancer patients due to the use of immune checkpoint inhibitors (ICIs), primary resistance still poses a huge challenge to clinicians. Bile acid metabolism is involved in the pathogenesis of CRC. However, its particular function in shaping the tumor immune microenvironment (TIME) and its effect on prognosis and immune treatment response remain unclear. METHODS: Based on the transcriptome and clinical data from The Cancer Genome Atlas-Colon Adenocarcinoma (TCGA-COAD) cohort, we performed unsupervised consensus clustering and classified patients into different molecular subtypes according to bile acid metabolism. We subsequently compared overall survival (OS), immune cell infiltration levels, and differentially expressed genes among the subtypes. In addition, protein-protein interaction (PPI) network and Cox proportional hazards regression were used to identify key hub genes. Finally, the expression of these crucial hub genes was validated in the Gene Expression Omnibus (GEO) cohort and independent clinical patients. RESULTS: The bile-low group showed a significant reduction in OS time (p = 0.0049). The infiltration levels of CD8+ T cells (p < 0.05) and M1 macrophages (p < 0.01) were significantly higher in the bile-low group than in the bile-high group. We identified three key genes-CLCA1, UGT2A3, and ZG16-and found that they all were downregulated in tumor tissues across the TCGA-COAD and GEO datasets, as well as in independent clinical samples. Survival analysis showed that high CLCA1 expression was significantly associated with favorable overall survival (p < 0.001), whereas UGT2A3 (p = 0.23) and ZG16 (p = 0.17) did not reach statistical significance. The three hub genes were negatively correlated with the (TIDE) score (CLCA1: R = - 0.24, p < 0.001; UGT2A3: R = - 0.15, p = 0.0022; ZG16: R = - 0.14, p = 0.0039). CONCLUSION: Our findings suggest that bile acid metabolism could shape the TIME via key genes CLCA1, UGT2A3, and ZG16, and subsequently modify CRC prognosis and immunotherapy responses. These genes may serve as potential prognostic indicators and mechanistic mediators linking bile acid metabolism to T-cell dysfunction, offering insights for future combination strategies targeting the metabolism-barrier-immunity axis.

CLCA1↗

PubMed2026

Distinct spatial immune microenvironment features of different EGFR mutation subtypes in early-stage lung adenocarcinoma.

Epidermal growth factor receptor (EGFR) mutations are common in lung adenocarcinoma (LUAD), yet their influence on the spatial tumor immune microenvironment (TIME) in early-stage disease remains unclear. We characterized the spatial TIME in 144 treatment-naïve, early-stage LUADs using integrated genomic sequencing and multiplex immunohistochemistry (mIHC). Although EGFR-mutant tumors overall displayed reduced CD8 + T-cell infiltration compared with EGFR-wild-type tumors, substantial heterogeneity was observed among EGFR subtypes. Specifically, L858R and rare-variant subtypes exhibited higher tumor mutational burden, greater CD8 + T-cell density, and enrichment of T-cell-dominant cellular neighborhoods relative to 19del subtype, consistent with a comparatively immune-infiltrated phenotype. In contrast, 19del tumors showed lower T-cell infiltration. TP53 co-mutation was also associated with enhanced CD8 + T-cell infiltration. These cross-sectional findings identify hypothesis-generating spatial immune phenotypes across EGFR-mutant LUAD subtypes; their potential relevance to perioperative treatment selection requires prospective validation in outcome-annotated treatment cohorts.

Humans↗