SEARCH · Search PubMed

Results for “single cell transcriptomics”

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 91 records · Page 5Linked to original sources

PubMed2026

Single-cell profiling reveals epithelial and immune responses in BK polyomavirus-infected human kidney biopsies.

INTRODUCTIONBK polyomavirus (BKV) infection is associated with injury and subsequent graft loss due to the extent of injury or rejection. However, the molecular mechanisms driving injury and subsequent adverse outcomes remain poorly understood.METHODSIn a cross-sectional study, single-cell RNA-seq from kidney allograft biopsies was used to assess cell type-specific responses between uninfected controls and 2 distinct phases of BKV infection: peaking (increasing viral blood titers) and resolving (decreasing viral titers following immunosuppression reduction).RESULTSGenes upregulated in BK viral nephropathy (BKVN) were enriched for polyomavirus infection hallmarks, including ribosome biogenesis, translation, and energy restructuring. Additionally, enriched pathways included wound healing, cellular stress, antigen presentation and immune signaling. Even without BKVN (peaking BK viremia alone), epithelial cells expressed signatures for wound healing, cellular stress, and extracellular matrix remodeling. In vivo tubular cell responses at single-cell resolution were validated against single cell transcriptomic data of BKV-infected cells in a cell culture model. Despite similarities, in vivo tubular cells underwent metabolic adaptation favoring fatty acid oxidation and proinflammatory responses not observed in culture models, likely due to an absent innate and adaptive immune system. Despite lymphopenia and immunosuppressive therapies, the proportion of recipient-derived intrarenal adaptive immune cells was increased in biopsies associated with peaking viremia alongside activation of innate immune responses. Adaptive immune cells exhibited persistent inflammatory signaling and remodeling of energy metabolism during the resolving phase of infection.CONCLUSIONThese not previously reported insights into BKV-associated injury may have implications for clinical management and improved allograft outcomes.

Humans↗

PubMed2026

Integrative genetic and transcriptomic analyses prioritize CDC16 as a candidate marker for gastric cancer.

BackgroundGastric cancer (GC) remains a major cause of cancer-related mortality, and biomarkers for early detection are needed.MethodsStomach and blood expression quantitative trait loci were integrated with two GC genome-wide association studies using Mendelian randomization (MR), Bayesian colocalization, and summary-data-based MR/heterogeneity in dependent instruments (SMR/HEIDI) testing. Bulk and single-cell transcriptomic analyses characterized candidate expression and lesion-associated patterns. CDC16 protein expression was evaluated by immunohistochemistry in 53 paired GC and non-neoplastic tissues, followed by paired and exploratory receiver operating characteristic analyses.ResultsMR prioritized PILRB, CDC16, and GABPB1-AS1; SMR/HEIDI provided complementary support, while colocalization for CDC16 and GABPB1-AS1 was suggestive and model-dependent. Bulk-tissue CDC16 abundance was higher in GC, but the modest TCGA-STAD tumor-normal difference (log2FC = 0.210, FDR = 0.019) was attenuated after proliferation adjustment (log2FC = -0.002, FDR = 0.987), indicating close coupling with proliferative activity. Single-cell analysis localized CDC16 predominantly to epithelial populations, and the proportion of CDC16-detectable epithelial cells increased across lesion categories (ρ = 0.735; permutation P = 0.031). CDC16 H-scores were higher in GC than in paired non-neoplastic tissues (161.15 ± 45.11 vs 102.15 ± 54.50; P < 0.001), with higher cancer-tissue scores in 41 of 53 cases. Exploratory AUC was 0.794 (95% CI, 0.704-0.874; sensitivity, 66.0%; specificity, 79.2%).ConclusionsConvergent genetic, transcriptomic, and protein-level evidence prioritizes CDC16 as a GC-associated candidate tissue marker whose expression is closely linked to proliferative activity. Prospective validation in independent cohorts, including appropriate disease controls and blood-based evaluation, is warranted.

Stomach Neoplasms↗

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

Integrated Bulk and Single-Cell RNA-Seq Analysis Reveals Transcriptional Activation of PTGS2 by FOS in Progression From T2DM to T2DM-Associated NAFLD.

Type 2 diabetes mellitus (T2DM) and nonalcoholic fatty liver disease (NAFLD) frequently coexist, exacerbating disease burden. However, the molecular mechanisms underlying the progression from T2DM to T2DM-associated NAFLD remain unclear. This study investigated the regulatory function of FOS-mediated PTGS2 activation in this transition. We integrated bulk RNA-seq data from GEO, single-cell transcriptomic data and transcriptomes from patients with T2DM-associated NAFLD. Differentially expressed genes were identified using the limma package, and T2DM-related gene modules were defined by weighted gene co-expression network analysis. LASSO regression and random forest identified 14 candidate genes, with PTGS2 and FOS prioritised. Single-cell analysis showed increased FOS and PTGS2 expression in monocytes, CD8+ T cells and Kupffer cells. Transcription factor prediction and dual-luciferase assays confirmed that FOS directly binds the PTGS2 promoter and drives its transcription. In vitro, FOS silencing decreased PTGS2 expression, cytokine secretion and apoptosis under high-glucose and free fatty acid conditions, whereas PTGS2 overexpression exacerbated inflammation and apoptosis independently of FOS expression. These findings demonstrate that FOS transcriptionally activates PTGS2, contributing to hepatic inflammation and apoptosis during the progression from T2DM to NAFLD. PTGS2 may serve as a promising biomarker and therapeutic target for T2DM-associated NAFLD.

Single-Cell Gene Expression Analysis↗

PubMed2023

Emerging trends in the study of spiralian larvae.

Many animals undergo indirect development, where their embryogenesis produces an intermediate life stage, or larva, that is often free-living and later metamorphoses into an adult. As their adult counterparts, larvae can have unique and diverse morphologies and occupy various ecological niches. Given their broad phylogenetic distribution, larvae have been central to hypotheses about animal evolution. However, the evolution of these intermediate forms and the developmental mechanisms diversifying animal life cycles are still debated. This review focuses on Spiralia, a large and diverse clade of bilaterally symmetrical animals with a fascinating array of larval forms, most notably the archetypical trochophore larva. We explore how classic research and modern advances have improved our understanding of spiralian larvae, their development, and evolution. Specifically, we examine three morphological features of spiralian larvae: the anterior neural system, the ciliary bands, and the posterior hyposphere. The combination of molecular and developmental evidence with modern high-throughput techniques, such as comparative genomics, single-cell transcriptomics, and epigenomics, is a promising strategy that will lead to new testable hypotheses about the mechanisms behind the evolution of larvae and life cycles in Spiralia and animals in general. We predict that the increasing number of available genomes for Spiralia and the optimization of genome-wide and single-cell approaches will unlock the study of many emerging spiralian taxa, transforming our views of the evolution of this animal group and their larvae.

Animals↗

PubMed2026

Genomic characterization of a hypervirulent Aeromonas veronii NN0115 from Nile tilapia and head kidney transcriptome of infected fish reveals B-cell-dominated immune response with specific immunoglobulin downregulation.

Aeromonas veronii is a pathogen of multiple fish species, yet systematic understanding of its infection in Nile tilapia (Oreochromis niloticus) remains limited. A dominant strain, NN0115, was isolated from a natural outbreak and identified as A. veronii by 16S rRNA and whole-genome average nucleotide identity (ANI, 96.33%). Experimental infection revealed high virulence (LD50 = 3.41 × 106 CFU/mL, equivalent to 8.53 × 104 CFU/fish). The genome is 4.58 Mb (58.57% GC) and encodes 4216 proteins. Virulence factor analysis identified 1253 genes, dominated by motility-related (264) and immune modulation (208) factors. Genomic island GI2 harbors 7 virulence genes and two dual-function resistance-virulence genes. The strain is resistant to 9 of 25 agents tested but carries three RND efflux pump genes whose predicted resistance was not phenotypically observed. The head kidney transcriptome of tilapia at 24 h post-bacterial infection identified 773 differentially expressed genes; among them, 57 were immunoglobulin (Ig) genes, and 56 were down-regulated. Integration of published single-cell transcriptomic data showed that non-Ig B-cell marker genes were down-regulated by 32%, whereas Ig genes were reduced by 63%, indicating selective transcriptional suppression of Ig genes rather than a general decrease in B-cell transcriptional activity. Together, this study provides a comprehensive characterization of a highly virulent A. veronii from Nile tilapia and reveals that selective downregulation of B-cell Ig genes is the dominant transcriptional feature of the host head kidney response.

Animals↗

PubMed2025

Trajectory inference from single-cell genomics data with a process time model.

Single-cell transcriptomics experiments provide gene expression snapshots of heterogeneous cell populations across cell states. These snapshots have been used to infer trajectories and dynamic information even without intensive, time-series data by ordering cells according to gene expression similarity. However, while single-cell snapshots sometimes offer valuable insights into dynamic processes, current methods for ordering cells are limited by descriptive notions of "pseudotime" that lack intrinsic physical meaning. Instead of pseudotime, we propose inference of "process time" via a principled modeling approach to formulating trajectories and inferring latent variables corresponding to timing of cells subject to a biophysical process. Our implementation of this approach, called Chronocell, provides a biophysical formulation of trajectories built on cell state transitions. The Chronocell model is identifiable, making parameter inference meaningful. Furthermore, Chronocell can interpolate between trajectory inference, when cell states lie on a continuum, and clustering, when cells cluster into discrete states. By using a variety of datasets ranging from cluster-like to continuous, we show that Chronocell enables us to assess the suitability of datasets and reveals distinct cellular distributions along process time that are consistent with biological process times. We also compare our parameter estimates of degradation rates to those derived from metabolic labeling datasets, thereby showcasing the biophysical utility of Chronocell. Nevertheless, based on performance characterization on simulations, we find that process time inference can be challenging, highlighting the importance of dataset quality and careful model assessment.

Single-Cell Analysis↗

PubMed2026

Machine learning-integrated multi-omics risk prediction for pulmonary fungal infection in COPD and lung cancer: a transcriptomic and immune profiling study.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) and lung cancer are major risk factors for invasive pulmonary fungal infection (IPFI), carrying an attributable mortality of 30%-80%. Their coexistence further amplifies immunosuppression, while current diagnostic criteria remain inadequate for early risk identification. METHODS: Transcriptomic data from the GEO dataset GSE296912 (scRNA-seq; 12,078 cells from normal and COPD lung tissue) and The Cancer Genome Atlas (TCGA)-lung adenocarcinoma (LUAD) bulk RNA-seq cohort (539 tumor and 59 normal samples) underwent differential expression and cross-omics integration analysis. Five machine learning models were constructed: logistic regression, SVM, random forest, XGBoost, and LASSO. Candidate genes were validated by qRT-PCR in A549 cells and THP-1-derived macrophages stimulated with heat-inactivated Aspergillus fumigatus conidia, a protocol selected to ensure BSL-2 biosafety compliance and isolate PAMP-mediated innate immune signaling. Model performance was evaluated using 5-fold stratified cross-validation with AUC, calibration curves, and decision curve analysis. RESULTS: Single-cell transcriptomic analysis of 12,078 cells identified 14 distinct cell populations, with marked myeloid expansion and immune dysregulation in COPD lung tissue. Cross-omics integration with TCGA-LUAD data identified 1,145 shared genes (79 immune-related), converging on NF-κB, TLR4, and cytokine receptor signaling. The random forest model achieved excellent discriminative performance (5-fold CV AUC = 0.988), with Treg infiltration, TLR4, and MMP9 as the top predictors. qRT-PCR confirmed significant upregulation of all five candidate genes (DEFB4A, S100A8, IL-8, MMP9, and TLR4) in both A549 and THP-1 cells following fungal stimulation. CONCLUSION: This multi-omics machine learning model integrating scRNA-seq and TCGA transcriptomic data demonstrates excellent discriminative performance (AUC = 0.988), with mechanistic convergence of NF-κB, TLR4, and oncogenic signaling pathways identified across shared immune gene signatures. In vitro qRT-PCR validation confirms the biological relevance of five key antifungal immune genes, providing a transcriptomic foundation for future prospective IPFI risk stratification in patients with COPD and lung cancer.

TLR4↗

PubMed2026

Expression patterns of potential targets for antibody-directed therapy in metastatic castration-resistant prostate cancer patients.

INTRODUCTION: Survival in metastatic castration-resistant prostate cancer (mCRPC) patients remains limited and treatment is complicated by tumor heterogeneity. As antibody-based therapeutics emerge, identifying actionable antigen targets and patient subgroups most likely to benefit is essential. MATERIALS & METHODS: Gene expression of 62 antibody-targetable proteins was analyzed in 296 mCRPC biopsies. These genes encode proteins targeted by approved or investigational antibody-based cancer therapeutics. Associations between target expression with genomic classifications and transcriptomic subtypes were evaluated. Target expression was also assessed in tumors with low expression of established mCRPC targets. Subgroup-specific targets were validated in an independent cohort and single-cell transcriptomics. RESULTS: Established targets KLK2, FOLH1 (PSMA) and STEAP1 showed the highest median expression across the cohort. Target expression did not correlate with genomic classifications, including homologous recombination deficiency, microsatellite instability, CDK12, TP53, PTEN or AR alterations Target expression did associate with transcriptomic subtypes: CRPC-AR (driven by androgen receptor-signaling) and CRPC-SCL (stem cell-like features, AP-1/YAP/TAZ-driven), displayed the highest expression of multiple targets, including KLK2, FOLH1, and SLC44A4. CRPC-NE (neuroendocrine phenotype) showed heterogeneous expression, with high CD46 expression, whereas CRPC-WNT (Wnt-signaling driven) generally showed low target expression. Notably, CD46 was highly expressed in tumors with low KLK2, FOLH1, and STEAP1 expression, a subgroup associated with poor prognosis. CONCLUSIONS: Although several antibody targets showed broad expression in mCRPC-tumors, expression varied by transcriptomic subtype. Subgroups such as CRPC-WNT expressed fewer targets, suggesting the need for alternative therapeutic strategies. CD46 emerged as a promising target, with wide expression across multiple subtypes, including clinically challenging CRPC-NE and mCRPC tumors lacking expression of established targets.

Humans↗

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

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

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

Humans↗

PubMed2026

Thyroid hormone deprivation creates an immunological signature in the mouse liver, involving Kupffer cell presentation as the mouse ages.

PURPOSE: Aging is associated with an increased prevalence of chronic liver diseases suggesting impaired immune and metabolic function. In addition, thyroid hormone (TH) impacts liver physiology and TH deprivation or excess negatively affect organ maintenance. However, whether age-dependent consequences of TH alterations are reflected in a liver-specific adaptation is unknown so far. The present study aimed to characterize the impact of TH deprivation or excess on the liver transcriptome during aging. METHODS: Five- and 21-month-old male C57BL/6 mice were exposed either to chronic TH deprivation or to chronic TH excess and compared to control treatment by microarray-based liver transcriptome analysis. RESULTS: Significant roles of both TH state and age became obvious: Bioinformatic analysis of the liver transcriptome data revealed an age-dependent immune signature by chronic TH deprivation, an age-dependent immune and metabolic signature independent of exogenous TH modulation, as well as an age-dependent metabolic signature by chronic TH excess. Published data of single cell transcriptomic atlas characterizing aging tissues in the mouse were compared with our data and revealed Kupffer cell presentation in the immunological signature by TH deprivation during aging. Literature data for four prominent differentially expressed genes, namely C1qb, C3ar1, Ctss, and Msr1, revealed that the complement system, extracellular matrix remodelling, as well as the proinflammatory phenotype of Kupffer cells are altered by TH deprivation during aging. CONCLUSION: In conclusion, our study illuminates the interplay between TH deprivation, aging, and liver transcriptome signatures, highlighting potential implications for immune function and tissue maintenance, particularly through the modulation of Kupffer cell presentation.

Animals↗

PubMed2024

Single cell RNA sequencing provides novel cellular transcriptional profiles and underlying pathogenesis of presbycusis.

Age-related hearing loss (ARHL) or presbycusis is associated with irreversible progressive damage in the inner ear, where the sound is transduced into electrical signal; but the detailed mechanism remains unclear. Here, we sought to determine the potential molecular mechanism involved in the pathogeneses of ARHL with bioinformatics methods. A single-cell transcriptome sequencing study was performed on the cochlear samples from young and aged mice. Detection of identified cell type marker allowed us to screen 18 transcriptional clusters, including myeloid cells, epithelial cells, B cells, endothelial cells, fibroblasts, T cells, inner pillar cells, neurons, inner phalangeal cells, and red blood cells. Cell-cell communications were analyzed between young and aged cochlear tissue samples by using the latest integration algorithms Cellchat. A total of 56 differentially expressed genes were screened between the two groups. Functional enrichment analysis showed these genes were mainly involved in immune, oxidative stress, apoptosis, and metabolic processes. The expression levels of crucial genes in cochlear tissues were further verified by immunohistochemistry. Overall, this study provides new theoretical support for the development of clinical therapeutic drugs.

Animals↗

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

Downregulated lysyl oxidase in plasma extracellular vesicles: a biomarker linked to brain metastasis risk in lung adenocarcinoma.

BACKGROUND: Brain metastasis (BrM) is a leading cause of mortality in patients with lung adenocarcinoma (LUAD). Extracellular vesicles (EVs), which carry bioactive molecules, play a critical role in tumor microenvironment remodeling and exhibit metastatic organotropism, holding promise as liquid biopsy biomarkers. This study aims to identify plasma EV-derived proteins associated with LUAD-BrM. METHODS: A multi-omics framework was applied. Plasma EVs from 59 stage IV LUAD patients (30 BrM vs 29 non-BrM) were profiled using data-independent acquisition mass spectrometry proteomics. Candidate proteins were screened via bioinformatics and machine learning (LASSO/RF/SVM). Initial validation included tissue proteomics (n = 13), single-cell transcriptomics (TISCH2), and Western blot analysis of a subset of the discovery samples. Functional experiments were conducted in vitro. The lead candidate was ultimately validated in an independent plasma cohort (n = 158) through ELISA. RESULTS: Proteomic analysis implicated collagen-containing extracellular matrix (ECM) pathways. Lysyl oxidase (LOX), a key ECM cross-linking enzyme, was identified as a lead candidate. LOX and its family member LOXL1 were consistently downregulated in BrM tissues and plasma EVs. Single-cell analysis revealed decreased LOX expression specifically in BrM-associated fibroblasts, which showed suppressed ECM-related pathways. In vitro experiments supported a PI3K/AKT-LOX-ECM regulatory axis. Plasma EV-derived LOX demonstrated strong diagnostic performance in the independent cohort, with an AUC of 0.786 (95% CI 0.713iated fi. CONCLUSIONS: Our study establishes plasma EV-derived LOX as a promising non-invasive biomarker for LUAD-BrM through a comprehensive multi-omics validation strategy. We propose a model wherein downregulation of LOX, potentially driven by PI3K/AKT signaling in tumor-associated fibroblasts, contributes to ECM degradation and may promote brain-tropic metastasis. This finding offers new insights for risk stratification and timely intervention in LUAD patients.

Humans↗

PubMed2026

From transcriptomic profiling to precision oncology: a bibliometric analysis of RNA sequencing in acute myeloid leukemia.

BACKGROUND: RNA sequencing (RNA-seq) has become an important tool for investigating the molecular heterogeneity of acute myeloid leukemia (AML); however, the global development and thematic evolution of this field remain inadequately characterized. OBJECTIVE: To map the global landscape of AML RNA-seq research and identify major knowledge domains, emerging themes, and temporal changes in research priorities. METHODS: Publications indexed in the Web of Science Core Collection and Scopus between January 1, 2007, and August 18, 2025, were retrieved. After database filtering, merging, and deduplication, 3,460 articles and reviews were included. CiteSpace, VOSviewer, the bibliometrix R package, and Microsoft Excel were used to analyze publication trends, collaboration networks, co-citation structures, keyword evolution, and citation bursts. RESULTS: Publication output increased steadily, accelerating after 2014. China contributed the largest number of publications (n = 547, 15.8%), whereas the United States had the highest total citation count. Major publication outlets spanned hematology, oncology, genomics, and molecular biology. Co-citation analysis identified prominent themes involving next-generation sequencing, gene mutations, KMT2A rearrangements, epigenetic dysregulation, leukemia-initiating cells, drug resistance, biomarkers, T-cell biology, and single-cell sequencing. Earlier literature emphasized sequencing technologies, gene expression profiling, and molecular alterations, whereas recent publications show increasing representation of cellular heterogeneity, single-cell transcriptomics, drug resistance, biomarker applications, immune-related research, and computational interpretation. CONCLUSION: While molecular characterization remains foundational, AML RNA-seq research has broadened to encompass increasingly prominent cellular, functional, computational, and translational dimensions. This study provides a structured overview of the field; nevertheless, bibliometric prominence should not be interpreted as direct evidence of clinical utility.

RNA sequencing↗

PubMed2026

A CFH- and SPINT2-based prognostic signature for cholangiocarcinoma.

BACKGROUND: Cholangiocarcinoma (CCA) is a highly malignant tumor with a poor prognosis, and reliable biomarkers for postoperative risk stratification remain limited. This study aimed to develop and validate a CFH- and SPINT2-based prognostic signature to support postoperative risk stratification and inform adjuvant therapy selection in CCA through integrative machine learning and single-cell transcriptomics. METHODS: Differentially expressed genes were screened from GSE26566. Integrative machine learning (least absolute shrinkage and selection operator-Cox, random forest, and univariate Cox regression) was performed in the training cohort (GSE89749; n=115) to construct a risk model, which was externally validated in two independent cohorts: cohort 1 (E-MTAB-6389; n=75) and cohort 2 [The Cancer Genome Atlas Cholangiocarcinoma (TCGA-CHOL) data set; n=36]. Systematic analysis was conducted and included examinations of immune infiltration [via single-sample gene set enrichment analysis (ssGSEA)], pathway enrichment (via hallmark GSEA), cellular localization (via single-cell RNA sequencing), and drug sensitivity (via the Genomics of Drug Sensitivity in Cancer 2 database). RESULTS: Two genes, CFH and SPINT2, were identified and incorporated into a prognostic risk score. High-risk patients in the training cohort had a significantly worse overall survival (log-rank P=0.02). External validation was performed in two independent cohorts. In validation cohort 1, the risk group was an independent prognostic factor [hazard ratio =2.27, 95% confidence interval (CI): 1.18-4.37; P=0.01]. In validation cohort 2, the model demonstrated acceptable discriminative ability (concordance index =0.721; 3-year area under the curve =0.692). The high-risk group exhibited an immunosuppressive microenvironment characterized by increased infiltration of macrophages and myeloid-derived suppressor cells, along with the activation of epithelial-mesenchymal transition, inflammatory response, and NF-κB signaling pathways. Single-cell analysis revealed a cell-type-specific expression pattern: CFH was predominantly expressed in fibroblasts, while SPINT2 was mainly expressed in malignant cells. Drug sensitivity analysis demonstrated that the high-risk group was more sensitive to gemcitabine, cisplatin, poly(ADP-ribose) polymerase (PARP) inhibitors, and mammalian target of rapamycin (mTOR) inhibitors, whereas the low-risk group was more sensitive to lapatinib. CONCLUSIONS: The CFH- and SPINT2-based prognostic signature may serve as an independent biomarker for postoperative risk stratification in CCA. High-risk patients, characterized by fibroblast-derived CFH enrichment and malignant-cell SPINT2 loss, exhibit an immunosuppressive microenvironment and may be more suitable for gemcitabine-based chemotherapy or PARP/mTOR inhibitors, whereas low-risk patients may benefit from less intensive adjuvant strategies or HER2/EGFR-targeted lapatinib. Prospective validation is warranted before clinical implementation.

Cholangiocarcinoma (CCA)↗

PubMed2026

Systematic discovery of retina-enriched Rik genes identifies 1190005I06Rik as a novel modulator of visual signalling.

BACKGROUND: High‑throughput transcriptome projects have revealed thousands of mammalian genes with little or no functional annotation. Among these are hundreds of loci assigned provisional “Rik” identifiers following discovery in the RIKEN cDNA annotation effort. Although often dismissed as genomic dark matter, such genes may encode tissue‑restricted proteins that modulate physiologic functions and influence disease. The retina is a highly specialised neural tissue and a common site of inherited disorders; understanding its molecular repertoire could illuminate novel therapeutic avenues. METHODS: We integrated bulk RNA‑seq from ten adult mouse tissues, evolutionary and domain analysis, single‑cell RNA‑seq, and CRISPR/Cas9 gene disruption to systematically catalogue protein‑coding Rik genes enriched in the retina and test the function of a representative gene. RESULTS: A rigorous differential expression analysis identified 44 Rik genes with robust retina‑specific expression compared with nine non‑retinal tissues. Many of these genes lack orthologues beyond rodents, while others show broad conservation, illustrating a continuum from lineage‑restricted to conserved retinopathy candidates. Single‑cell transcriptomics revealed that these genes are expressed across retinal cell types, with the highest aggregate expression in cone photoreceptors and inner interneurons. To evaluate physiological significance, we generated a 1190005I06Rik knockout mouse. Although retinal architecture appeared normal, loss of 1190005I06Rik enhanced electroretinogram b‑wave amplitudes and altered light‑avoidance behaviour, indicating that this previously uncharacterised gene acts as a negative modulator of visual signalling. CONCLUSIONS: We present a curated atlas of retina‑enriched Rik genes and demonstrate that 1190005I06RIK modulates retinal circuit function. This resource expands the molecular landscape of the retina and provides new candidates for the genetic basis of inherited retinal disease. Our findings underscore that unannotated genes may exert measurable effects on sensory processing and warrant systematic exploration in the context of human ocular disorders.

Animals↗