SEARCH · Search PubMed

Results for “Transcriptomic profiling”

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 271 records · Page 15Linked to original sources

PubMed2025

An anti-androgen resistance-related gene signature acts as a prognostic marker and increases enzalutamide efficacy via PLK1 inhibition in prostate cancer.

BACKGROUND: Anti-androgen resistance remains a major clinical challenge in the treatment of prostate cancer (PCa), leading to disease progression and treatment failure. Despite extensive research on resistance mechanisms, a reliable prognostic model for predicting patient outcomes and guiding therapeutic strategies is still lacking. This study aimed to develop a novel gene signature related to anti-androgen resistance and evaluate its prognostic and therapeutic implications. METHODS: Anti-androgen resistance-related differentially expressed genes (ARRDEGs) were identified through transcriptomic analysis of enzalutamide- and dual enzalutamide abiraterone-resistant PCa cell lines from the GEO database. Functional enrichment analysis was performed to determine the biological roles of these genes. A prognostic gene signature was developed using univariate Cox regression, LASSO, and multivariate Cox regression models. The model was validated in independent PCa cohorts from The Cancer Genome Atlas (TCGA). Additionally, we assessed the correlation between the signature, immune infiltration, immune checkpoint expression, and drug sensitivity. The efficacy of PLK1 inhibition combined with enzalutamide was further explored using in vitro and in vivo experiments. RESULTS: We identified 304 ARRDEGs, from which three key genes (LMNB1, SSPO, and PLK1) were selected to construct a prognostic signature. This gene signature effectively stratified PCa patients into high- and low-risk groups, with the high-risk group exhibiting shorter recurrence-free survival and distinct immune characteristics. High-risk patients demonstrated elevated immune checkpoint expression (B7H3, CTLA-4, B7-1, and TIGIT), increased M2 macrophage infiltration, and enhanced sensitivity to chemotherapy and targeted therapy. Mechanistically, PLK1 inhibition potentiated the antitumor effect of enzalutamide by downregulating SLC7A11 and inducing ferroptosis, providing a potential therapeutic strategy to overcome anti-androgen resistance. CONCLUSION: We established a novel ARRDEGs-based prognostic signature that predicts PCa progression and response to chemotherapy and targeted therapy. The integration of this signature with immune profiling and drug sensitivity analysis provides a valuable tool for precision oncology in PCa. Our findings highlight the potential of PLK1 inhibition as a therapeutic strategy to enhance enzalutamide efficacy and overcome resistance.

Humans↗

PubMed2026

Liver transcriptome analysis revealed multiple immune processes and lipid metabolism pathways involved in the defense response of the turbot (Scophthalmus maximus) against Aeromonas salmonicida.

Aeromonas salmonicida is a significant pathogen causing notable economic losses in Scophthalmus maximus aquaculture. This study utilized Illumina sequencing technology to examine the transcriptional response characteristics of S. maximus liver at 24 h following A. salmonicida infection. A total of 2363 differentially expressed genes (DEGs) were identified when compared to the negative control group. The immunity-related Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, as well as metabolism-related PPAR signaling pathway and insulin signaling pathway, were notably enriched. Significant differences exist in the expression of key genes within the PPAR pathway, particularly cd36, acsl4a, pparαa, and plin2, all of which mediate the interaction between lipid metabolism and the immune response. These results offer valuable insights into the immunometabolic regulatory mechanism of S. maximus response to A. salmonicida infection.

Animals↗

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↗

PubMed2025

Higher Expression of HPV16 Derived E7_LI Transcript Observed in Men With HIV and Recurrent Anal Cancer.

Squamous cell carcinoma of the anus (SCCA) or anal cancer (AC) is an understudied cancer with a high occurrence rate in people with HIV (PWH), especially men having sex with men (MSM). Furthermore, AC recurs in approximately one-fourth of patients who undergo standard care with chemoradiation therapy (CRT). Using bulk RNA sequencing data of AC obtained from 12 patients with non-recurrent (NR, N = 9) or recurrent (R, N = 3) cancer, we previously showed upregulated expression of key immune genes in the NR compared to the R group. Although the main causative agent of AC is high-risk human papillomavirus (HPV), association of host and viral RNA transcript expression contributing to AC recurrence has not been extensively studied. The objective of the current study was to determine whether enrichment of specific HPV genotypes and/or HPV gene expression patterns differentiate the two groups and if any specific viral (HPV) and host (human) immune mediators correlate with each other. Using bulk RNA sequencing data and VIRTUS 2, we detected viral RNA reads mapping to seven high-risk and six low-risk HPV types, of which the high-risk HPV16 observed in 83% (10/12) AC tumors (7/9 NR and 3/3 R). Rate of all HPV genomes trended toward a decrease in NR AC isolates and correlation between HPV types was more commonly observed in low-risk ones. Analysis of HPV 16 gene expression profile showed a significantly lower positivity rate for a polycistronic transcript encoding for E7^L1 in the NR group (1/9, NR vs. 3/3, R, p < 0.05). An unbiased correlation analysis of HPV-human transcript expression showed a direct correlation between HPV transcripts and human genes involved in cell growth. The data also identified human transcripts showing an inverse correlation with HPV gene expression. These included genes involved in negative regulation of growth, proliferation, and immune response. Taken together, these data indicate that concurrent analyses of viral and host factors in the same tumor can identify potential new therapeutic targets to ameliorate cancer recurrence post-treatment.

Humans↗

PubMed2025

iSoMAs: Finding isoform expression and somatic mutation associations in human cancers.

Aberrant alternative splicing, prevalent in cancer, impacts various cancer hallmarks involving proliferation, angiogenesis, and invasion. Splicing disruption often results from somatic point mutations rewiring functional pathways to support cancer cell survival. We introduce iSoMAs (iSoform expression and somatic Mutation Association), an efficient computational pipeline leveraging principal component analysis technique, to explore how somatic mutations influence transcriptome-wide gene expression at the isoform level. Applying iSoMAs to 33 cancer types comprising 9,738 tumor samples in The Cancer Genome Atlas, we identified 908 somatically mutated genes significantly associated with altered isoform expression across three or more cancer types. Mutations linked to differential isoform expression occurred through both cis- and trans-acting mechanisms, involving well-known oncogenes/suppressor genes, RNA binding protein and splicing factor genes. With wet-lab experiments, we verified direct association between TP53 mutations and differential isoform expression in cell cycle genes. Additional iSoMAs genes have been validated in the literature with independent cohorts and/or methods. Despite the complexity of cancer, iSoMAs attains computational efficiency via dimension reduction strategy and reveals critical associations between regulatory factors and transcriptional landscapes.

Humans↗

PubMed2025

PoweREST: Statistical power estimation for spatial transcriptomics experiments to detect differentially expressed genes between two conditions.

Recent advancements in spatial transcriptomics (ST) have significantly enhanced biological research in various domains. However, the high cost for current ST data generation techniques restricts the large-scale application of ST. Consequently, maximization of the use of available resources to achieve robust statistical power for ST data is a pressing need. One fundamental question in ST analysis is detection of differentially expressed genes (DEGs) under different conditions using ST data. Such DEG analyses are performed frequently, but their power calculations are rarely discussed in the literature. To address this gap, we developed PoweREST, a power estimation tool designed to support the power calculation for DEG detection with 10X Genomics Visium data. PoweREST enables power estimation both before any ST experiments and after preliminary data are collected, making it suitable for a wide variety of power analyses in ST studies. We also provide a user-friendly, program-free web application that allows users to interactively calculate and visualize study power along with relevant parameters.

Gene Expression Profiling↗

PubMed2025

Partners in root nodule symbiosis respond uniquely to heavy metal stresses in a host genotype-dependent manner.

The mutualistic symbiosis between legume roots and soil rhizobia culminates in the formation of root nodules, where nitrogen is fixed. Root nodule symbiosis is inhibited by heavy metal stress. In this study, we investigated the relative responses of the symbiotic partners to a non-essential heavy metal cadmium (Cd) and an essential heavy metal zinc (Zn) stress and identified patterns in gene expression. We performed dual transcriptomics in nodules, using the Medicago truncatula-Sinorhizobium meliloti symbiotic system. Phenotypes were measured in the wild-type Medicago truncatula and a mutant in an ABC transporter gene (Mtabcg36), which showed compromised nodule formation in control conditions and further after heavy metal treatment. We observed that the rhizobia were particularly sensitive to Zn in mutant nodules. The greatest degree of differential gene expression in the host plant were observed under Cd and Zn treatments in wild-type nodules. Most Cd-regulated host genes were also differentially regulated by Zn, revealing little discernment between an essential and a non-essential ion under increased exposure. Furthermore, the host response to both the stresses affected auxin and iron homeostasis genes in a host genotype-dependent manner. Our results suggested impaired cadmium export from the mutant nodules. These results have potential implications in agricultural management systems and bioremediation strategies.

Symbiosis↗

PubMed2026

Differential expression of neuronal function genes follows a tissue-specific temporal dynamic during Deformed Wing Virus infection in honey bees.

Deformed Wing Virus type A (DWV-A) is one of the primary threats to honeybees (Apis mellifera), significantly impacting their nervous system physiology and behavior. While its neurotropic nature is well-recognized, the temporal dynamics of the neuronal transcriptomic response following oral infection, the natural transmission route, remains poorly understood. In this study, we analyzed gene expression in the heads of worker bees orally inoculated with DWV-A over a 16-day time course (1, 4, 7, 10, 13, and 16 days post-inoculation). RNA-seq analysis at day 10 identified 147 differentially expressed genes associated with different biological processes that are critical to the organism, including cellular metabolism and neuronal activity. RT-qPCR validation revealed a persistent downregulation of key genes related to glutamatergic system (eaat-2, neto, and kainate) and sensory perception-related genes in the antennae. Notably, the simultaneous co-expression of nurse-associated and forager-associated marker genes suggests that DWV-A infection induces an asynchrony in behavioral maturation. Our findings demonstrate that DWV-A disrupts neuronal homeostasis and peripheral sensory perception in a tissue-specific and time-dependent manner, providing a molecular framework to understand the behavioral impairment and the loss of coordination at the colony level.

Animals↗

PubMed2026

Multi-omics unveils seasonal remodeling and metabolic crosstalk between testis and abdominal fat body in a non-amplexus stream frog Nanorana taihangnica (Anura: Dicroglossidae).

BACKGROUND: Energy allocation between reproduction and survival represents a fundamental life-history challenge for animals in seasonal environments. Using integrated transcriptomics and metabolomics, we investigated Nanorana taihangnica (Anura: Dicroglossidae), a non-amplexus stream frog endemic to China, to elucidate the seasonal morphological and molecular coordination between the testis and abdominal fat body. RESULTS: Morphological analysis showed that fat body adipocyte cross-sectional area minimized at the end of the breeding season but rapidly recovered thereafter, while testicular volume continued declining post-breeding and only recovered during the non-breeding period. During breeding season, multi-omics analyses revealed that the fat body enhanced fatty acid oxidation, upregulated histidine-carnosine metabolism, activated NAD+ metabolism and FOXO3-mediated antioxidative responses to mitigate metabolic stress, and regulated adipocyte survival and apoptosis via sphingolipid signaling. Seasonal testicular development was centrally regulated by the mTOR signaling pathway, whose activity integrated autophagy levels, NAD+ availability, and aspartate metabolism to coordinate spermatogonial proliferation and spermatogenesis. CONCLUSIONS: This study demonstrates that N. taihangnica optimizes seasonal energy storage, allocation, and reproductive investment through molecular and metabolic crosstalk between the fat body and testis, providing empirical insights into the physiological integration of life-history strategies in animals inhabiting fluctuating environments.

Animals↗

PubMed2026

Nonadditive gene expression and reduced homoeolog expression bias in an intraspecific hexaploid wheat hybrid.

BACKGROUND: Intraspecific hybridization in allopolyploid plants can generate additive and nonadditive changes in gene expression through interactions between divergent parental genomes. However, how it simultaneously affects gene expression and the relative expression of homoeologs in higher-order polyploids is less well understood. To study this, we sequenced seedling leaf transcriptomes and profiled gene body methylation in two hexaploid wheat (Triticum aestivum L.) cultivars and their F₁ hybrids. RESULTS: Although only 4.3% of genes differed in expression between the parents, 22.3% deviated from mid-parent expression in the hybrids, with many showing transgressive expression. 32.1% of triads contained at least one homoeolog that deviated from mid-parent expression, and all three homoeologs deviated in 11% of triads, substantially more than expected by chance. Triads in which all three homoeologs were overexpressed also showed reduced differences in expression among homoeologs. Greater parental divergence in relative homoeolog expression was associated with nonadditive expression. Genes lacking gene body methylation were also more likely to show dominant or transgressive expression, whereas gene body methylation was associated with more balanced homoeolog expression and additive or conserved expression. CONCLUSIONS: Intraspecific hybridization in hexaploid wheat, even without a change in ploidy, was associated with widespread nonadditive gene expression and altered relative homoeolog expression within triads. These responses were associated with parental differences in homoeolog expression and the absence of gene body methylation. Although our findings are limited to seedling leaves from a single intraspecific cross, they provide a basis for testing the generality of these patterns across tissues, developmental stages, and genetic backgrounds.

Triticum↗

PubMed2026

Comparative transcriptomic analysis of the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress.

To investigate the differences in molecular responses between the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress, a 30-day chronic stress experiment was conducted with a control group and a stress group. Transcriptomic analysis of the gills and hepatopancreas was performed using Illumina sequencing; differentially expressed genes (DEGs) were identified, and GO, KEGG, GSEA, PPI, and RT-qPCR validation were carried out. The results showed that 196 DEGs (161 up-regulated and 35 down-regulated) were identified in the gills, and 287 DEGs (199 up-regulated and 88 down-regulated) in the hepatopancreas, with only 18 DEGs shared between the two tissues. DEGs in the gills were enriched in oxidoreductase activity, glycerophospholipid metabolism, and tyrosine metabolism; DEGs in the hepatopancreas were enriched in lipid transporter activity, phagosome, ECM-receptor interaction, and riboflavin metabolism. GSEA revealed significant suppression of the mTOR pathway in the gills and the Polycomb complex pathway in the hepatopancreas. PPI network analysis identified hub genes P5CS and eEF2 in the gills, and PER, TUBB1, SHMT, and TUBB4B in the hepatopancreas. RT-qPCR validation was consistent with the RNA-seq results (R2 = 0.764). This study indicates that, under chronic nitrite stress, the gill response is centered on redox regulation and inhibition of growth metabolism, whereas the hepatopancreas response primarily involves lipid transport, cytoskeletal remodeling, and phagosome activation. The two tissues synergistically adapt through fundamental biosynthetic and motor protein pathways. This research provides molecular evidence for deciphering the nitrite tolerance mechanisms in freshwater-cultured shrimp.

Animals↗

PubMed2026

Identification of Critical Genes for Recurrent Aphthous Ulcer by Transcriptome Data Analysis and Mendelian Randomization.

PURPOSE: Recurrent aphthous ulcer (RAU) is a common oral mucosal disorder with a poorly understood etiology, significantly affecting patients' quality of life. This study aims to investigate critical genes linked to RAU and explore their biological mechanisms using transcriptomic data and Mendelian randomization (MR) analysis. MATERIALS AND METHODS: RAU-related gene expression data from the GEO database (GSE37265) were analyzed to identify differentially expressed genes (DEGs). A two-sample MR approach was used to assess the causal impact of expression quantitative trait loci (eQTL) on RAU. Critical genes were identified by intersecting DEGs with significant MR findings. GO and KEGG pathway enrichment analyses were performed, along with GSEA and immune cell infiltration analysis, to investigate the functions and mechanisms of these genes in RAU. RESULTS: A total of 184 differentially expressed genes (DEGs) were identified, while 339 RAU-associated genes were screened through MR analysis. Cross-validation further identified 7 critical genes. Among these, CCR1, ERP27, HCK, MICB, and SLC2A3 showed protective associations with RAU risk, whereas CD177 and IFITM1 were positively associated with increased risk. Enrichment analysis revealed that these genes are involved in specific biological processes, including cell migration, immune response, and metabolic regulation, which are closely linked to RAU pathogenesis. CONCLUSION: This systematic study comprehensively investigates the critical causative genes underlying RAU, emphasizing the intricate relationships between immune regulation and metabolic disturbances in its pathology. These findings lay a solid foundation for the development of novel biomarkers and may inform future research on targeted therapeutic strategies for RAU.

Stomatitis, Aphthous↗

PubMed2025

Inflammatory pathways and immune dysregulation in pediatric postoperative septic shock: A study integrating transcriptomics, machine learning and molecular docking.

This study elucidates the molecular and immune regulatory mechanisms of pediatric postoperative septic shock. Transcriptomic data were obtained from the Gene Expression Omnibus database. Differentially expressed genes were identified using the limma package, and gene co-expression modules were constructed using Weighted Gene Co-expression Network Analysis. Functional enrichment was performed via gene set enrichment analysis, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses. Immune cell infiltration was assessed using ESTIMATE and CIBERSORT. Mendelian randomization was applied to explore causal relationships between gene expression and septic shock. Feature genes were selected using machine learning algorithms, and a diagnostic nomogram model was constructed. Finally, molecular docking analysis was performed to screen and evaluate the binding affinity of traditional Chinese medicine monomers to core target proteins. A total of 1331 differentially expressed genes were identified, and the turquoise module was strongly correlated with septic shock. Enrichment analysis revealed significant activation of IL-6/JAK/STAT3, TNF-α/NF-κB, and PI3K/Akt/mTOR pathways. Immune infiltration analysis indicated suppressed immune scores and imbalances in neutrophils, macrophages, T cells, and B cells. Mendelian randomization confirmed causal associations for 6 genes, including PIM3. The predictive model based on feature genes demonstrated high diagnostic performance. Molecular docking suggested that quercetin and astramembrannin I could stably bind PIM3. This study systematically identified core genes, dysregulated immune pathways, and candidate small-molecule interventions in pediatric septic shock, providing novel insights for early diagnosis and targeted therapy.

Humans↗

PubMed2026

Integrative machine learning and transcriptomic analysis reveals molecular mechanisms underlying low survival rate in larval Chinese Bahaba (Bahaba taipingensis).

Chinese Bahaba (Bahaba taipingensis) is a Class I protected marine fish endemic to China. Low larvae survival during artificial breeding severely hinder population recovery. To investigate the molecular mechanism of high mortality in larval fish, this study performed RNA-seq on liver from naturally deceased (ND) and mass-dead (MD) individuals, combined with least absolute shrinkage and selection operator (LASSO) regression and random forest (RF) algorithms to screen for core signature genes. A total of 873 differentially expressed genes (DEGs) were identified, including 112 upregulated and 761 downregulated genes. GO and KEGG enrichment analyses revealed significant enrichment in amino acid metabolism disorders, one‑carbon folate pool impairment, PPAR signaling abnormalities, ECM-receptor interaction, focal adhesion pathway, indicating widespread metabolic suppression accompanied by extracellular matrix remodeling and signaling disturbances in the livers of MD fish. MAD pre-filtering combined with dual machine learning algorithms yielded 18 robust core signature genes, among which SLC38A4, MMP1, FADD, FKBP5, and APOB were consistently identified as high-frequency core genes by both algorithms. SLC38A4 exhibited the highest importance score in the RF model and was significantly downregulated, making it the primary molecule distinguishing ND from MD phenotypes. ROC curve analysis showed that both models achieved an AUC of 1.000 (95% CI lower bound: 0.610), confirming the precise discriminatory ability of the core genes. GSEA further demonstrated significant enrichment of this core gene set in ND samples. This study provides the first systematic elucidation of the molecular mechanisms underlying liver dysfunction in low survival rate B. taipingensis, characterized by amino acid transport impairment, metabolic reprogramming, and structural remodeling, offering theoretical foundations for health assessment, early mortality risk warning, and artificial breeding conservation of this species.

Animals↗

PubMed2025

Transcriptome-wide analysis reveals potential roles of CFD and ANGPTL4 in fibroblasts regulating B cell lineage for extracellular matrix-driven clustering and novel avenues for immunotherapy in breast cancer.

BACKGROUND: The remodeling of the extracellular matrix (ECM) plays a pivotal role in tumor progression and drug resistance. However, the compositional patterns of ECM in breast cancer and their underlying biological functions remain elusive. METHODS: Transcriptome and genome data of breast cancer patients from TCGA database was downloaded. Patients were classified into different clusters by using non-negative matrix factorization (NMF) based on signatures of ECM components and regulators. Weighted Gene Co-expression Network Analysis (WGCNA) was used to identify core genes related to ECM clusters. Additional 10 independent public cohorts including Metabric, SCAN_B, GSE12276, GSE16446, GSE19615, GSE20685, GSE21653, GSE58644, GSE58812, and GSE88770 were collected to construct Training or Testing cohort, following machine learning calculating ECM correlated index (ECI) for survival analysis. Pathway enrichment and correlation analysis were used to explore the relationship among ECM clusters, ECI and TME. Single-cell transcriptome data from GSE161529 was processed for uncovering the differences among ECM clusters. RESULTS: Using NMF, we identified three ECM clusters in the TCGA database: C1 (Neuron), C2 (ECM), and C3 (Immune). Subsequently, WGCNA was employed to pinpoint cluster-specific genes and develop a prognostic model. This model demonstrated robust predictive power for breast cancer patient survival in both the Training cohort (n = 5,392, AUC = 0.861) and the Testing cohort (n = 1,344, AUC = 0.711). Upon analyzing the tumor microenvironment (TME), we discovered that fibroblasts and B cell lineage were the core cell types associated with the ECM cluster phenotypes. Single-cell RNA sequencing data further revealed that angiopoietin like 4 (ANGPTL4)+ fibroblasts were specifically linked to the C2 phenotype, while complement factor D (CFD)+ fibroblasts characterized the other ECM clusters. CellChat analysis indicated that ANGPTL4+ and CFD+ fibroblasts regulate B cell lineage via distinct signaling pathways. Additionally, analysis using the Kaplan-Meier Plotter website showed that CFD was favorable for immunotherapy response, whereas ANGPTL4 negatively impacted the outcomes of cancer patients receiving immunotherapy. CONCLUSION: We identified distinct ECM clusters in breast cancer patients, irrespective of molecular subtypes. Additionally, we constructed an effective prognostic model based on these ECM clusters and recognized ANGPTL4+ and CFD+ fibroblasts as potential biomarkers for immunotherapy in breast cancer.

Humans↗

PubMed2025

Colorectal Liver Metastasis Pathomics Model: Integrating Single-Cell and Spatial Transcriptome Analysis With Pathomics for Predicting Liver Metastasis in Colorectal Cancer.

The liver is the primary target organ for hematologic metastasis of colorectal cancer (CRC), and CRC liver metastasis (CRLM) often precludes radical resection, making it the leading cause of death in patients with CRC. To improve the identification and prediction of liver metastasis risk, we identified a cell type of liver metastasis--triggering malignant cells (LMTMCs) through integrating single-cell RNA sequencing and spatial transcriptome analysis. Multiomics cell communication analysis indicated that the interaction between fibroblasts and LMTMCs through the COL1A1-CD44/SDC4 and LAMA4-CD44 signaling axes could promote CRLM. By applying the one-class logistic regression algorithm, we developed a CRLM scoring system in the bulk RNA-sequencing data according to the abundance of LMTMCs in each individual. Using the grouping labels derived from the CRLM scoring system in the bulk data and the corresponding whole-slide images without any manual annotations at the region or pixel level, processed via slide-level weakly supervised learning, a deep-learning model based on the ResNet18 architecture, called Colorectal Liver Metastasis Pathomics Model, was developed to predict the risk of liver metastasis in patients with CRC. The Colorectal Liver Metastasis Pathomics Model achieved an area under the curve of 0.84 at the internal test set of The Cancer Genome Atlas-CRC histology images. In the external independent validation sets, namely the Affiliated Hospital of Southwest Medical University and the Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University cohorts, the areas under the curve were 0.89 and 0.72, respectively, indicating effective classification performances. This study provided new insights and tools for the early identification of CRLM and demonstrated the potential of combining multiomics with deep learning-based pathomics in cancer research.

Humans↗

PubMed2026

Spatial transcriptomics of primary and metastatic ALK-rearranged NSCLC reveals site-specific adaptations.

INTRODUCTION: Genetic alterations and the tumor microenvironment (TME) influence treatment response in anaplastic lymphoma kinase-rearranged non-small cell lung cancer (ALK+ NSCLC). This study maps site-specific TME adaptations and exploratory risk-associated signatures in lymph node metastases (LNT) to investigate metastatic evolution. METHOD: We applied spatial transcriptomics to profile tumor (PanCK+) and stromal (PanCK-) compartments in a pilot cohort of 16 cases: primary lung tumors (LT, n = 3), LNT (n = 10), and brain metastases (BT, n = 3), with three site-matched non-tumor controls. LNT-derived prognostic signatures were evaluated using The Cancer Genome Atlas-Lung Adenocarcinoma (TCGA LUAD) cohorts. RESULTS: Distinct, site-specific TME features were observed. LNT stroma was enriched in fibroblasts and macrophages, while tumor segments showed increased neutrophils. BT exhibited a macrophage-associated immunosuppressive TME. Tumor cells evolved divergently: LT retained pulmonary identity and showed trend towards translation-associated programs, LNT cells shifted toward senescence and epigenetic remodeling, and BT cells showed activation of Class A/1 (Rhodopsin-like) receptor, GPCR and drug metabolism pathways. In LNT, exploratory risk-associated differences were observed. Low-risk cases (n = 6) showed adaptive immune signatures, whereas high-risk cases (n = 4) showed enrichment for stromal MET signaling and stress-response pathways. Because treatment exposure differed markedly between the risk groups, these observations should be interpreted as hypothesis-generating. TCGA LUAD analysis suggested the broader biological relevance of immune-associated markers, but reflected general LUAD rather than ALK+ specific biology. Discordant associations for GCLC and TIMP1 underscored the importance of spatial context. CONCLUSION: Site-specific microenvironments may influence tumor adaptation across metastatic niches in ALK+ NSCLC. The exploratory risk-associated findings require validation in larger, uniformly treated cohorts.

Humans↗

PubMed2026

An oxidative stress - and immunotherapy-related six-gene signature defines immune subtypes and predicts prognosis and immunotherapy response in hepatocellular carcinoma.

BACKGROUND: Oxidative stress and the tumor immune microenvironment jointly shape hepatocellular carcinoma (HCC) progression and response to immunotherapy, yet integrated biomarkers linking these processes are lacking. METHODS: Transcriptomic and clinical data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were used to identify oxidative stress- and immunotherapyrelated differentially expressed genes (OSIRDEGs). Functional enrichment, weighted gene co-expression network analysis (WGCNA) and LASSO-Cox regression were used to construct a prognostic signature. Consensus clustering, TIDE, CIBERSORT and ssGSEA characterized immune phenotypes. Somatic mutation, copy-number and drug-response data were integrated to assess genomic alterations and drug sensitivity. Expression of model genes was validated by qRT-PCR and western blotting in HCC cell lines. RESULTS: We identified 24 OSIRDEGs enriched in cell-cycle and mitotic pathways. WGCNA intersection yielded 18 module genes, from which a six-gene signature (BUB1B, CDKN2A, CENPE, HMMR, PTTG1, SPP1) was derived. The signature robustly stratified patients into high- and low-risk groups with significantly different progression-free and disease-free survival in both TCGA-LIHC and GSE14520. Based on signature expression, two molecular subtypes were defined, exhibiting distinct survival, immune landscapes and predicted immunotherapy responsiveness. Model genes harbored recurrent alterations and showed significant correlations with anticancer agents. All six genes were upregulated at mRNA and protein levels in metastatic HCC cell lines versus normal hepatocytes. CONCLUSIONS: We systematically explored the landscape of OSIRDEGs in HCC, and proposed a validated six-gene signature that refines prognostic stratification, delineates immunerelevant HCC subtypes and highlights candidate biomarkers for therapeutic selection and mechanistic investigation.

Humans↗