Search PubMed⌕ Search

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.

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

At least 541 records · Page 30Linked to original sources
PubMed2026

Multi-Omics Biomarker Signatures for Precision Diagnosis and Prognosis in Primary Liver Cancer: A Literature Review.

Primary liver cancer (PLC) is a biologically heterogeneous group of malignancies dominated by hepatocellular carcinoma (HCC), intrahepatic cholangiocarcinoma (iCCA), and a smaller subset of combined hepatocellular-cholangiocarcinoma (cHCC-CCA), and its clinical burden remains high because current diagnostic and prognostic tools do not adequately capture molecular diversity. Conventional imaging, serum markers, and histopathological assessment remain insufficient for precise early diagnosis, subtype-resolved classification, and outcome stratification, while tissue and liquid biopsy approaches have expanded the range of analytes available for clinical assessment. Recent studies have identified candidate biomarker signatures across genomic, epigenomic, transcriptomic, proteomic, metabolomic, and circulating layers, suggesting that integrated multi-omics profiling may better represent tumor lineage, clonal evolution, immune context, and therapeutic vulnerability than isolated molecular readouts. However, these layers are not equally mature for clinical use: genomic testing is closest to routine therapeutic application in iCCA, plasma methylation assays are advancing for HCC surveillance augmentation, and many proteomic or metabolomic panels remain validation-stage tools. Their clinical value remains constrained by sampling bias, biospecimen-dependent signal loss, assay standardization, cost, and the need for prospective validation across clinically diverse populations. This narrative review critically synthesizes current evidence on multi-omics biomarker signatures for precision diagnosis and prognosis in primary liver cancer and argues that clinically useful signatures should be question-specific, stage-aware, and specimen-aware rather than universal multi-analyte panels.

Humans↗
PubMed2026

Integrative multi-omics profiling of insomnia-related molecular features reveals microbiome, immune, and therapy-relevant heterogeneity in colorectal cancer.

Emerging evidence implicates insomnia as a potential risk factor in carcinogenesis, potentially involving systemic inflammation, circadian disruption, and microbiome alterations. However, the molecular associations linking insomnia-related features to colorectal cancer (CRC), particularly with respect to tumor biology, immune microenvironmental states, and therapy-relevant phenotypes, remain largely unexplored. Multi-omics integration of genomic, transcriptomic, and microbiome data from 3,026 CRC patients across seven independent cohorts, including a large, well-annotated Clinical Omics study of Colorectal Cancer in China (COCC) cohort, enabled insomnia-based molecular classification through unsupervised non-negative matrix factorization (NMF) clustering. The insomnia subtype (IS) was biologically characterized via pathway enrichment, immune deconvolution, microbial profiling, and single-cell transcriptomics. Furthermore, an insomnia score (ISscore) was developed and validated in multiple cohorts for risk stratification and assessment of treatment-response-related indicators in CRC. Unsupervised clustering revealed two distinct molecular subtypes (IS1/IS2), with IS2 demonstrating significantly poorer survival. IS2 exhibited marked activation of EMT/angiogenesis pathways versus cell cycle activation in IS1. The IS2 microenvironment showed increased immunosuppression-related infiltration and exhausted T cell signatures, together with intratumoral microbiome variation characterized by depletion of Ruminococcaceae UCG-002 and enrichment of Hungatella/Selenomonas. The ISscore system stratified survival risk and was associated with computational indicators of immunotherapy response. Single-cell analysis nominated PPIA-BSG as a potential cell-cell communication signal involving high-ISscore tumor cells, CXCL12+ endothelial cells, and CLEC9A+ dendritic cell subsets. This multi-omics characterization of insomnia-CRC interplay suggests that insomnia-related molecular features are associated with an immunologically distinct and microbiome-altered tumor ecosystem. The ISscore provides a reproducible framework for capturing insomnia-related molecular heterogeneity, supporting risk stratification and future evaluation of therapy-relevant phenotypes.IMPORTANCEChronic insomnia affects millions, but it is not typically considered a cancer risk factor. Our study, analyzing vast biological data from over 3,000 colorectal cancer patients, uncovers a potential link between a person's predisposition to insomnia and their risk of developing this disease. This suggests that the biological pathways related to sleep may play a role in cancer development. Understanding this connection opens up new avenues for identifying individuals at higher risk and developing novel prevention strategies for colorectal cancer.

colorectal cancer↗
PubMed2026

Next-Generation Disease Profiling by Integrating Histopathology with Spatial Multi-Omics Data.

The field of pathology has experienced several transformative changes in recent years with the advent of digital pathology and spatial multi-omics. These technologies have enhanced every aspect of pathology practice, from streamlining daily workflows to generating high-fidelity multi-omics data that provide pathologists with novel tools to refine disease profiling and clinical diagnosis. Each layer of multimodal data (genomic, metabolomic, proteomic, or transcriptomic) has uncovered a distinct facet of disease pathologies, and combined with machine learning/artificial intelligence-based data analysis and pattern recognition models, has provided holistic understanding of regulatory mechanisms underpinning them. However, high-dimensional data have far exceeded the volume, scale, and complexity of immunostaining methods implemented by pathologists and, thus, have generated significant challenges related to deconvolution, interpretation, and clinical translation. Furthermore, these multimodal studies have predominantly relied on computational methods to process data and extract disease-relevant insights, thus raising questions around relevance or role of a pathologist in this new era of multi-omics. This review will provide a perspective on the evolving fields of molecular histopathology and spatial -omics, leveraging them to approach disease profiling, and redefining the role of a pathologist during this process.

Humans↗
PubMed2026

Progressive salinity drives flavonoid branch reprogramming in Anoectochilus roxburghii.

Flavonoids play critical roles in plant adaptation to abiotic stress; however, how salt stress modulates metabolic flux distribution within flavonoid branches remains poorly understood, particularly in non-model medicinal plants. Here, we integrated targeted metabolomics, transcriptomics, and proteomics to examine flavonoid regulation in Anoectochilus roxburghii under 0, 50, 100, and 200 mmol·L- 1 NaCl. Metabolite profiling showed that salinity reshaped flavonoid composition rather than uniformly increasing flavonoid abundance. A metabolite-derived branch bias index (MI), representing the balance between reductive branch metabolites and flavonol products, increased under salt treatment, peaked at 100 mmol·L- 1 NaCl, and declined at 200 mmol·L- 1, indicating maximal branch bias under moderate stress followed by partial rebalancing under severe stress. Transcriptomic analysis showed induction of upstream phenylpropanoid and flavonoid entry genes, including PAL, 4CL, and CHS, whereas F3H was suppressed and FLS showed no induction. Furthermore, several short-chain dehydrogenase/reductase homologs (IFR-like SDR homologs) were upregulated, and the transcript-derived reductive branch index (EI) increased progressively across the salt gradient. EI was positively associated with MI, although the relationship was not strictly proportional under severe stress (200 mmol·L- 1 NaCl). Proteomic profiling further provided supportive evidence for sustained activation of upstream flavonoid biosynthesis, such as salt-induced accumulation of chalcone synthase (CHS) protein, complementing the transcriptomic and metabolomic datasets. Together, these results indicate that salt stress reorganizes flavonoid metabolism in A. roxburghii through persistent upstream activation and branch-specific regulation, favoring the reductive branch under moderate salinity.

Orchidaceae↗
PubMed2026

Transcriptomic responses to developmental temperature in two field-collected Spodoptera exigua populations from Korea.

The beet armyworm, Spodoptera exigua, is a polyphagous insect whose development and seasonal occurrence are strongly influenced by temperature. However, transcriptomic responses to developmental thermal regimes remain insufficiently characterized in field-collected populations. In this study, we compared two Korean field-collected populations of S. exigua: a Haenam population collected in May and initially maintained at 15 ± 1 °C (HN), and a Jeju population collected in July and initially maintained at 27 ± 1 °C (JJ). F1 larvae from each population were reared under three fluctuating developmental temperature regimes: low (15-21 °C), middle (21-27 °C), and high (27-33 °C), followed by RNA-seq analysis. Differential expression analysis revealed population-associated variation in transcriptomic responses across developmental temperatures. HN exhibited a larger number of differentially expressed genes under the high-temperature regime, suggesting stronger transcriptomic sensitivity to elevated developmental temperature. Functional enrichment analyses identified population-associated differences in pathways related to heat response, oxidative metabolism, cytoskeletal organization, cuticle-associated processes, lipid metabolism, and immune-related functions. In JJ, heat-response and cuticle-related expression patterns were more prominent under warmer developmental conditions, whereas HN showed broader changes in stress- and metabolism-associated pathways under high temperature. Overall, this study provides a comparative transcriptomic analysis of two field-collected S. exigua populations under different developmental temperature regimes and identifies RNA-seq-based molecular response patterns associated with population-specific thermal response profiles.

Animals↗
PubMed2026

Esketamine multi-omic biomarker evaluation in major depressive disorder (EMBER-MDD): concept, objectives and methodologies of a non-clinical investigator-initiated study.

Treatment resistance (TR) in major depressive disorder (MDD) affects a substantial minority of patients and is hard to recognize early, delaying intensified care. The Esketamine multi-omic biomarker evaluation in MDD (EMBER-MDD) is a non-interventional, investigator-initiated, in-vitro study within the EU Psych-STRATA programme, analyzing biospecimens collected in the randomized INTENSIFY study and the mirror OBS-TR cohort after participants complete treatment. EMBER-MDD aims to discover individual-omic and integrated multi-omic (hypothesis-free) biomarkers and signatures associated with TR risk, and molecular correlates of clinical response to esketamine nasal spray versus treatment as usual (TAU). Biomaterials will derive from approximately 420 adults with MDD (estimated n = 210 esketamine; n = 210 TAU) and include whole blood, RNA-stabilized whole blood, plasma and serum, sampled at baseline and, when feasible, during and after treatment (up to ~ 5,040 aliquots stored at - 80 °C). Genomics will use baseline DNA genotyping on Illumina Infinium GSA v3.0+MD arrays; epigenomics will profile genome-wide DNA methylation across time points using MethylationEPIC v2.0; transcriptomics will employ mRNA-seq (NovaSeq X/ X Plus); and proteomics/ metabolomics will be generated using high-throughput Olink and/ or Biocrates platforms. Each layer will undergo state-of-the-art preprocessing and analyses (e.g., GWAS/ PRS, EWAS, differential expression, WGCNA, pathway and network analyses), followed by integrative strategies including QTL mapping (meQTL/ eQTL/ pQTL/ mQTL) and intermediate-fusion machine learning with nested cross-validation, explainable AI (SHAP/ LIME) and treatment-effect modelling. All outputs are research-only and will not support individual efficacy, tolerability, or clinical decision-making. The study will deliver robust biosignatures and mechanistic hypotheses to guide future validation and inform stratified, molecularly guided intervention strategies in subsequent prospective trials. Trial registration number: 2023-506617-21-00 and 2025-178-f-S.

Humans↗
PubMed2026

Polyethylene transformation by a psychrotolerant Rhodococcus strain assessed by transcriptomics and 13C-isotope tracing.

Polyethylene is increasingly accumulating in nature, including remote places like the Arctic. While abiotic processes fragment polyethylene in situ, biotic transformation by microorganisms is assumed to occur. However, the enzymes and pathways involved remain poorly characterized. In this study, we used an in-house biobank from cold environments to screen for potential bacteria capable of degrading polyethylene by screening the strains in silico using the database PlasticDB and in vivo using a fluorescence-based assay. Using transcriptomic and proteomic analyses to identify genes in promising candidate strains that encode extracellular enzymes potentially capable of degrading PE, we selected a Rhodococcus erythropolis strain and two of its enzymes: a hypothetical protein (Hypr1) and a lipase family protein (Lip2). Expressing the candidate genes heterologously in Escherichia coli resulted in positive results in the fluorescence-based assay for polyethylene transformation. Applying 13C-labelled polyethylene for assessing and estimating polyethylene transformation and carbon assimilation, we found that R. erythropolis and both untransformed and recombinant E. coli extracellularly transformed the initially added polyethylene after 70 days. In addition, untransformed E. coli and R. erythropolis converted small, but significant amounts of polyethylene-derived carbon to carbon dioxide. The 13C-label was also traced into the bacterial biomass of R. erythropolis. Overall, our results provide evidence for biotic transformation of untreated polyethylene and suggests a hypothetical protein and a lipase family protein as two novel enzyme candidates associated with PE transformation.

Rhodococcus↗
PubMed2026

Dynamic effects of short-term storage temperature and duration on the transcriptome and functional pathways of umbilical cord blood-derived NK cells.

OBJECTIVE: This study aims to elucidate the dynamic impact of different storage conditions (temperature and time) on the function of natural killer (NK) cells derived from umbilical cord blood (UCB) at the transcriptome level, providing a theoretical basis for optimizing the standardized post-collection processing protocol of UCB in clinical settings. METHODS: Four healthy full-term UCB samples were collected and assigned to a fresh control group (0H) and experimental groups stored at 4 °C (4C) or 25 °C (RT) for 24 h (24H) and 72 h (72H). Umbilical cord blood mononuclear cells (CBMCs) were isolated and expanded in vitro to derive NK cells. Using RNA sequencing (RNA-seq) technology, combined with principal component analysis (PCA), screening of differentially expressed genes (DEGs), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, the effects of storage conditions on the transcriptome of NK cells were comprehensively analyzed. RESULTS: The results indicated that storage time was the primary factor causing transcriptomic differences. Compared to the 0H group, storage for 24 h induced a limited number of Differentially Expressed Genes (DEGs) (only 11 in the 24H-RT group), indicating overall transcriptome stability. In contrast, storage for 72 h triggered profound transcriptomic reprogramming, with DEGs significantly enriched in immune activation-related pathways. Specifically, at the 72-h time point, storage at 4 °C compared to 25 °C prevented the downregulation of key metabolic pathways such as 'response to decreased oxygen levels' and 'regulation of leukocyte cell-cell adhesion'. CONCLUSION: CBMC-derived NK cells can tolerate up to 24 h of ex vivo storage while maintaining transcriptome and core functional stability. However, extending the ex vivo storage period to 72 h leads to significant reprogramming in the expression of immune and metabolism-related genes. Notably, the degree of metabolic suppression under refrigerated conditions at 4 °C was less pronounced than that at 25 °C room temperature storage. This study recommends that the ex vivo storage time of UCB should ideally be limited to within 24 h, with room temperature storage being a viable option. If prolonged storage is necessary, refrigeration at 4 °C is advised.

Humans↗
PubMed2025

Ossifying Spindled and Epithelioid Tumor: A Novel Soft Tissue Tumor.

This investigation describes the clinicoradiologic, pathologic, and molecular features of a unique soft tissue tumor characterized by a peripheral shell of bone and composed of bland myoid spindle and epithelioid cells that are keratin-positive. Our study cohort consists of 6 men and 6 women, with a mean age of 32 years. The tumors arose in the extremities (n = 9) and proximal limb girdle (n = 3) and were equally distributed between deep and superficial soft tissues. Patients reported dull painless masses of several months to >10 years duration (mean: 2.9 years). Imaging demonstrated a complete or partial peripheral shell of bone that could extend centrally, and the tumor's mean size was 5.7 cm. Histologically, the tumors were composed of uniform, eosinophilic myoid spindled cells growing in sheets and intersecting fascicles, surrounded by mature lamellar and/or woven bone. Also present was an admixed component of intermediate-sized epithelioid cells with eosinophilic cytoplasm. Mitotic activity was consistently low. Immunohistochemistry showed strong multifocal staining for keratins, and 50% (5/10) showed focal staining for S100; however, all were negative for SMA, desmin, SOX10, ERG, and CD34. Genetic analysis by multiple targeted RNA sequencing panels was negative (n = 10); however, whole transcriptome sequencing (n = 8) revealed a recurrent and novel in-frame SRSF7::NFATC3 fusion in 4 tumors. Dual fluorescence in situ hybridization probes for SRSF7::NFATC3 successfully confirmed this fusion and identified a fifth case, which had not undergone whole transcriptome sequencing but was negative by a targeted RNA fusion panel. Methylation profiling (n = 8) demonstrated a shared epigenetic profile distinct from other entities. Clinical follow-up (n = 11) showed no evidence of recurrence after primary excision with a mean of 41.6 months. In summary, we describe a novel soft tissue tumor designated "ossifying spindled and epithelioid tumor" as a descriptive histologic term that also emphasizes its close radiologic mimic, ossifying fibromyxoid tumor. All cases have behaved in a benign fashion without recurrence following simple excision. Awareness of this entity is important, so that it can be distinguished from other neoplasms that have more aggressive biological potential.

Humans↗
PubMed2026

A pancreatic cancer organoid-macrophage co-culture using starPEG-heparin hydrogel deciphers tumor-immune cell interactions.

Macrophages are among the most abundant immune cells in the pancreatic ductal adenocarcinoma (PDAC) tumor microenvironment (TME) and play a key role in regulating the immunosuppressive niche that facilitates tumor growth. Although recent three-dimensional (3D) culture systems using patient-derived materials have advanced our understanding of tumor biology, most models lack key cellular TME components and thus fail to capture tumor-immune cell interactions. To address this gap, we developed an in-vitro 3D co-culture model incorporating PDAC patient-derived organoids (PDOs) and macrophages within a synthetic hydrogel matrix. We optimized culture conditions by tuning medium and matrix conditions to support both cell lineages. Flow cytometry and transcriptomic analyses revealed that initially undifferentiated macrophages adopt an M2-like profile upon exposure to PDAC PDOs in starPEG-heparin hydrogels, mirroring the macrophage phenotypes observed by multiplex immunohistochemistry in the matched primary PDAC tissues. Cytokine secretome profiling revealed PDO-specific differences, indicating distinct underlying macrophage polarization subtypes. Collectively, our starPEG-heparin hydrogel-based 3D co-culture enables hypothesis-driven and physiologically relevant studies of tumor-macrophage interactions and may advance immune-modulatory treatment strategies in patients with PDAC.

Journal Article↗
PubMed2024

Construction of a molecular diagnostic system for neurogenic rosacea by combining transcriptome sequencing and machine learning.

Patients with neurogenic rosacea (NR) frequently demonstrate pronounced neurological manifestations, often unresponsive to conventional therapeutic approaches. A molecular-level understanding and diagnosis of this patient cohort could significantly guide clinical interventions. In this study, we amalgamated our sequencing data (n = 46) with a publicly accessible database (n = 38) to perform an unsupervised cluster analysis of the integrated dataset. The eighty-four rosacea patients were partitioned into two distinct clusters. Neurovascular biomarkers were found to be elevated in cluster 1 compared to cluster 2. Pathways in cluster 1 were predominantly involved in neurotransmitter synthesis, transmission, and functionality, whereas cluster 2 pathways were centered on inflammation-related processes. Differential gene expression analysis and WGCNA were employed to delineate the characteristic gene sets of the two clusters. Subsequently, a diagnostic model was constructed from the identified gene sets using linear regression methodologies. The model's C index, comprising genes PNPLA3, CUX2, PLIN2, and HMGCR, achieved a remarkable value of 0.9683, with an area under the curve (AUC) for the training cohort's nomogram of 0.9376. Clinical characteristics from our dataset (n = 46) were assessed by three seasoned dermatologists, forming the NR validation cohort (NR, n = 18; non-neurogenic rosacea, n = 28). Upon application of our model to NR diagnosis, the model's AUC value reached 0.9023. Finally, potential therapeutic candidates for both patient groups were predicted via the Connectivity Map. In summation, this study unveiled two clusters with unique molecular phenotypes within rosacea, leading to the development of a precise diagnostic model instrumental in NR diagnosis.

Humans↗
PubMed2025

HighDimMixedModels.jl: Robust high-dimensional mixed-effects models across omics data.

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm's behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl.

Algorithms↗
PubMed2025

Transcriptome-wide root causal inference.

Root causal genes correspond to the first gene expression levels perturbed during pathogenesis by genetic or non-genetic factors. Targeting root causal genes has the potential to alleviate disease entirely by eliminating pathology near its onset. No existing algorithm has been designed to discover root causal genes from observational data alone. We therefore propose the Transcriptome-Wide Root Causal Inference (TWRCI) algorithm that identifies root causal genes and their causal graph using a combination of genetic variant and unperturbed bulk RNA sequencing data. TWRCI uses a novel competitive regression procedure to annotate cis and trans-genetic variants to the gene expression levels they directly cause. The algorithm simultaneously determines the sequence in which gene expression changes propagate through the system to pinpoint the underlying causal graph and estimate root causal effects. TWRCI outperforms alternative approaches across a diverse group of metrics by directly targeting root causal genes while accounting for distal relations, linkage disequilibrium, patient heterogeneity and widespread pleiotropy. We demonstrate the algorithm by uncovering the root causal mechanisms of two complex diseases, which we confirm by replication using independent genome-wide summary statistics.

Algorithms↗
PubMed2026

Systematic characterization of neurotransmitter receptor dysregulation identifies a neural-related prognostic signature associated with biochemical recurrence in prostate cancer.

BACKGROUND: The nervous system is increasingly recognized to play a critical role in tumor initiation and progression. Central to this complex relationship are the interactions between neurotransmitters secreted by neurons and their receptors (neurotransmitter receptors, NTRs) expressed on cancer cells, which activate multiple intracellular signaling pathways. However, the spectrum of NTR dysregulation and its association with biochemical recurrence (BCR) in prostate cancer (PCa) has not been explored. Therefore, the aim of this study was to fill this gap. METHODS: We systematically characterized the expression profiles of 130 NTR genes by integrating bulk and single-cell transcriptomic data. Consistently dysregulated NTR (cdNTR) genes were identified and used to construct a PCa signature (PCaSig) using elastic-net regression. The robustness of PCaSig was evaluated across three independent cohorts. In addition, the associations of PCaSig with clinicopathological characteristics, genomic alterations, tumor immune-related characteristics, and biological pathways were comprehensively investigated. RESULTS: Thirteen cdNTR genes with strong cell-type specificity, particularly in luminal epithelial cells, were identified. PCaSig robustly stratified patients into distinct BCR risk groups across multiple independent cohorts and remained an independent predictor after adjustment for clinicopathological factors. High PCaSig scores were associated with aggressive clinicopathological features, elevated tumor mutation burden (TMB), suppression of neurotransmitter-related signaling, and activation of cell-cycle and immune-related pathways. Notably, PCaSig refined prognostic stratification regardless of TMB status and was associated with distinct immune-related characteristics, including immune checkpoint expression and immune cell infiltration. Incorporation of PCaSig into a clinical nomogram significantly improved prognostic accuracy and clinical net benefit. CONCLUSIONS: These findings establish NTR dysregulation as a previously underappreciated dimension of PCa and support PCaSig as a clinically relevant tool for personalized management.

Neurotransmitter receptor (NTR)↗
PubMed2026

Cyclin-dependent kinase 4 and 6 inhibitors and the breast cancer immune ecosystem: immune remodeling, resistance, and therapeutic reprogramming.

Cyclin-dependent kinase 4 and 6 inhibitors (CDK4/6 inhibitors) combined with endocrine therapy have become a therapeutic backbone for hormone receptor-positive, human epidermal growth factor receptor 2-negative breast cancer, yet durable disease control is frequently limited by intrinsic and acquired resistance. Canonical tumor-cell mechanisms, including retinoblastoma-pathway escape, cyclin E-cyclin-dependent kinase 2 (CDK2) activation, endocrine adaptation, and phosphoinositide 3-kinase (PI3K)-AKT-mechanistic target of rapamycin (mTOR) signaling, explain only part of this failure because they do not fully capture dynamic immune and stromal remodeling. Preclinical and translational studies indicate that early CDK4/6 inhibition can enhance antigen presentation, activate interferon-related programs, restrain regulatory T cells, and promote a T-cell-inflamed state. These effects are conditional and may not persist during prolonged treatment. Sustained therapy can instead drive heterogeneous resistant niches characterized by stromal remodeling, myeloid recruitment, checkpoint adaptation, and T-cell dysfunction. This immune-state dependence provides a rationale for immune checkpoint blockade, although clinical combinations have shown mixed efficacy and clinically relevant hepatic, pulmonary, and hematologic toxicities. Sequential or lead-in strategies therefore warrant prospective evaluation. Oxidative phosphorylation (OXPHOS) and redox adaptation may sustain selected resistant states and expose context-dependent ferroptotic vulnerabilities. Ferroptosis may connect tumor-cell killing with immune regulation, whereas nanomedicine may improve tumor-selective delivery. Both strategies remain largely preclinical and require further evaluation of pharmacokinetics, biodistribution, toxicity, manufacturability, and immune-cell safety. This Review distinguishes intrinsic from acquired resistance across interpatient, intratumoral, spatial, and temporal dimensions. It integrates tumor-cell escape with cytokine, immune, stromal, vascular, and metabolic remodeling and summarizes emerging therapeutic strategies. We further propose a candidate biomarker-informed framework that integrates genomic profiling, spatial immune architecture, circulating biomarkers, T-cell receptor (TCR) dynamics, transcriptomic and single-cell analyses, artificial intelligence (AI)-assisted multimodal integration, and longitudinal sampling. This framework is intended to support biomarker development and prospective trial design rather than current clinical decision-making, providing a translational basis for testing state-informed and sequence-aware therapeutic strategies.

Humans↗
PubMed2026

Multistage Genetic, Transcriptomic, and Single-Cell Evidence Prioritizes MAP1LC3A among Ferroptosis-Related Genes in Glioblastoma.

Glioblastoma (GBM) remains a highly aggressive malignancy, and the contribution of ferroptosis-related genes to disease susceptibility remains incompletely understood. A genetically anchored, multistage framework was applied to prioritize ferroptosis-related genes associated with GBM. Among 483 genes curated from FerrDb V2, 315 had candidate cis-expression quantitative trait loci (cis-eQTLs) in eQTLGen, 250 retained at least three independent instruments after linkage disequilibrium clumping, and 226 yielded valid inverse-variance weighted (IVW) Mendelian randomization estimates using a GBM genome-wide association study comprising 6,183 cases and 18,169 controls. Thirty-four genes met the exploratory discovery criteria of P < 0.05 and a Benjamini-Hochberg false discovery rate (BH-FDR) < 0.20, with directionally concordant Bayesian weighted Mendelian randomization (BWMR) estimates. Replication-stage Mendelian randomization using GTEx V10 whole-blood cis-eQTLs supported four genes: ATG7, RPTOR, MAP1LC3A, and CHMP6. Evaluation across three independent tumor-control transcriptomic cohorts demonstrated that MAP1LC3A was consistently downregulated in tumor tissue and showed a significant random-effects pooled estimate (log&#x2082; fold change, -1.273; 95% confidence interval, -1.625 to -0.920; false discovery rate = 0.016), whereas the other three genes lacked comparable cross-cohort statistical support. Single-cell virtual knockout analysis was subsequently performed in a patient-balanced subset of 2,400 malignant cells selected from 4,916 eligible cells across 20 adult IDH-wild-type GBM tumors. Across five independently seeded runs, 3, 15, 4, and 7 robust downstream genes were identified for ATG7, RPTOR, MAP1LC3A, and CHMP6, respectively. The resulting consensus sets comprised 17 unique genes, with RND3 shared across all four targets. Gene Ontology analysis indicated enrichment of cell-adhesion and cell-surface processes, whereas no KEGG or Reactome pathways remained significant after multiple-testing correction. Collectively, these findings prioritize MAP1LC3A for future experimental investigation while distinguishing genetic association, tumor-expression concordance, and computational perturbation from definitive evidence of causality or mechanism.

Humans↗
PubMed2026

Deep FLASH-seq profiling of purified canine sensory neurons uncovers species-specific signatures relevant to pain and itch.

Naturally occurring pain and itch disorders in the domestic dog represent an important and underexploited opportunity for translational sensory neuroscience. These conditions largely mirror human disease, highlighting the need for detailed comparative understanding of canine somatosensory neurobiology. Here, we present a single-cell transcriptomic characterisation of the canine dorsal root ganglion (DRG), providing molecular insights into sensory neuron diversity in a species of direct veterinary and biomedical relevance. We develop a novel mechanical dissociation and fluorescence-activated cell sorting strategy enabling purification of intact whole neurons from adult canine DRG, followed by deep, full-length RNA sequencing using FLASH-seq. This approach yields high-quality transcriptional profiles with molecular depth analogous to deep neuronal profiling in human DRG, enabling resolution of neuronal identities and subtype-specific gene programs. Using these data, we identify canine sensory neuron clusters conforming to conserved principles of DRG molecular organization observed across species, including peptidergic and noncanonical peptidergic nociceptors, low-threshold mechanoreceptors, proprioceptors, and thermosensory populations. Cross-species comparisons with human and mouse DRG datasets reveal broad conservation of pain- and itch-relevant pathways and therapeutic targets, alongside biologically meaningful divergence. We further identify species-specific differences in subtype-restricted expression of the pharmacologically relevant receptors IL31RA and SSTR2 , which we validate using in situ hybridization and contextualize with human spatial transcriptomic data. Finally, we provide evidence that domestication-associated genes are nonrandomly enriched in specific sensory neurons, suggesting that evolutionary history may have shaped somatosensory function. These data represent a resource for comparative sensory neuroscience and inform translational interpretation of pain and itch therapeutics across species.

Animals↗
PubMed2026

DPAS-Graph: adaptive spatial-feature relation learning for spatial RNA-to-protein prediction and virtual protein profiling.

Paired spatial multi-omics provides a supervised basis for learning RNA-protein correspondence in situ, but predicting protein abundance from spatial transcriptomic data alone remains challenging across tissue contexts and protein panels. Here, we present DPAS-Graph, an adaptive relation-learning framework for spatial RNA-to-protein prediction. Rather than directly merging spatial proximity and transcriptomic similarity as fixed graph priors, DPAS-Graph represents them as two relation channels on a shared edge support and updates their contributions during representation learning for protein prediction. Its Niche-Coupled Field Encoder combines layer-wise edge-relation modeling, intra-branch relation refinement, and cross-branch residual correction to learn spot representations for protein abundance prediction. In a leave-one-dataset-out benchmark across seven paired spatial multi-omics datasets, DPAS-Graph achieved lower aggregate prediction errors and improved spot-level agreement of protein expression profiles, with gains mainly reflected in error-based metrics and PCC-Spot. Spatial autocorrelation and protein-derived domain agreement analyses were further used to characterize the spatial behavior of the predicted protein maps. When applied to external RNA-only spatial sections, DPAS-Graph generated qualitatively interpretable marker-level virtual protein maps, illustrating its use as a complementary tool for protein-level interpretation of transcriptomics-only spatial data.

RNA↗