Search PubMed⌕ Search

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.

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

At least 415 records · Page 23Linked to original sources
PubMed2026

Metab8D: a metabolic regulome network from multiomics and machine learning.

To explore multiomic regulation of the metabolome, we used machine learning to predict metabolomic variation across ~1000 different cancer cell lines with matched omics data from eight biomolecular classes: genomic copy number variation, mutations, DNA methylation, histone post-translational modifications (PTMs), transcriptomics and RNA splice variants, non-coding transcriptomics (miRNA and lncRNA), proteomics, and phosphoproteomics. Overall, the metabolome is tightly associated with the transcriptome, with coding and non-coding RNAs emerging as top predictors. Peripheral metabolites are predictable via levels of corresponding enzymes, while those in central metabolism require combinatorial predictors in signaling and redox pathways, and may not reflect corresponding pathway expression. We reconstruct multiomic interaction subnetworks for highly predictable metabolites, and YAP1 signaling emerged as a top global predictor across four omic layers. We prioritize predictive multiomic features for single-cell and spatial metabolomics assays. Top predictors were enriched for synthetic-lethal interactions and synergistic combination therapies that target compensatory metabolic modulators.

Machine Learning↗
PubMed2026

Transcriptome Analysis and Experimental Validation of Palmitoylation- Related Biomarkers in Atherosclerosis.

INTRODUCTION: Protein palmitoylation contributes to membrane localisation, signal transduction, and cell-fate regulation. It is closely associated with lipid metabolic dysfunction, immune inflammation, and vascular remodelling in atherosclerosis (AS). However, key palmitoylation-related transcriptomic markers and their potential causal associations with AS remain incompletely defined. METHODS: The Gene Expression Omnibus (GEO) dataset GSE100927 was used as the training cohort, and GSE43292 was used as an external validation cohort. Differentially expressed genes were identified using limma and intersected with palmitoylation-related genes to obtain palmitoylation-related differentially expressed genes (PRDEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were then performed using clusterProfiler. Two-sample Mendelian randomisation was used to evaluate potential causal relationships between characteristic genes and AS. Feature selection was conducted using random forest and support vector machine recursive feature elimination (SVM-RFE), and the overlapping genes selected by both methods were retained. Receiver operating characteristic (ROC) curves were used to assess diagnostic performance. A five-gene nomogram was constructed, and its clinical utility was evaluated using calibration curves and decision curve analysis (DCA). Gene set variation analysis (GSVA) was applied to compare pathway activity between high- and low-expression groups for each core gene. Single-cell analysis using Seurat and expression-based cell-cell communication analysis using CellChat were conducted with GSE159677, and upstream transcription factors were predicted using NetworkAnalyst. For in vivo validation, an AS model was established in ApoE⁸/⁸ mice fed a high-fat diet, and aortic gene and protein expression were assessed by RT-qPCR and western blotting. RESULTS: In GSE100927, 51 PRDEGs were identified. GO and KEGG enrichment analyses highlighted pathways associated with regulation of monoatomic ion transport, sarcomere and myofibril organisation, and immune inflammation. Mendelian randomisation suggested a potential protective causal association between SLC7A7 and AS. By integrating MR with random forest and SVM-RFE feature selection, we prioritised five core genes: PLCB2, GMIP, NEXN, PLN, and SLC7A7. These genes showed good diagnostic performance in GSE43292. The resulting nomogram was well calibrated and demonstrated stable net benefit in decision curve and clinical impact curve analyses. Single-gene GSVA identified consistently activated pathways across multiple genes, including innate and adaptive immune recognition, calcium signalling and myocardial contraction/cardiomyopathy, extracellular matrix-receptor interaction, cell junction pathways, autophagy-lysosome pathways, and several metabolic programmes. At the single-cell level, PLCB2 and GMIP were predominantly expressed in T cells and macrophages, NEXN and PLN were enriched in vascular smooth muscle cells, and SLC7A7 was mainly expressed in macrophages. CellChat analysis indicated increased signals for immune-related ligand-receptor interactions. In ApoE⁸/⁸ mice fed a high-fat diet, PLCB2, GMIP, and SLC7A7 were upregulated, whereas NEXN and PLN were downregulated; protein-level changes were concordant with the transcriptomic trends. DISCUSSION: These findings indicate that palmitoylation-related dysregulation in AS converges on immune inflammation, calcium signalling/contractile programmes, ECM remodelling, and autophagy-linked metabolism. The five-gene panel is supported by external validation, single-cell localisation to immune and vascular compartments, and concordant results in ApoE⁸/⁸ mice. CONCLUSION: This study identified and validated five palmitoylation-related genes associated with AS. SLC7A7 showed a potential protective causal signal in MR analysis. The enriched pathway patterns linked these genes to immune inflammation, calcium signalling-contraction coupling, ECM remodelling, cell adhesion, and autophagy- associated metabolic reprogramming. The five-gene nomogram showed potential utility for diagnostic classification and decision support, nominating candidate biomarkers and pathway targets for AS molecular subtyping, diagnosis, and mechanistic investigation.

Atherosclerosis (AS)↗
PubMed2025

Self-organization of mouse embryonic stem cells into reproducible pre-gastrulation embryo models via CRISPRa programming.

Embryonic stem cells (ESCs) can self-organize into structures with spatial and molecular similarities to natural embryos. During development, embryonic and extraembryonic cells differentiate through activation of endogenous regulatory elements while co-developing via cell-cell interactions. However, engineering regulatory elements to self-organize ESCs into embryo models remains underexplored. Here, we demonstrate that CRISPR activation (CRISPRa) of two regulatory elements near Gata6 and Cdx2 generates embryonic patterns resembling pre-gastrulation mouse embryos. Live single-cell imaging revealed that self-patterning occurs through orchestrated collective movement driven by cell-intrinsic fate induction. In 3D, CRISPRa-programmed embryo models (CPEMs) exhibit morphological and transcriptomic similarity to pre-gastrulation mouse embryos. CPEMs allow versatile perturbations, including dual Cdx2-Elf5 activation to enhance trophoblast differentiation and lineage-specific activation of laminin and matrix metalloproteinases, uncovering their roles in basement membrane remodeling and embryo model morphology. Our findings demonstrate that minimal intrinsic epigenome editing can self-organize ESCs into programmable pre-gastrulation embryo models with robust lineage-specific perturbation capabilities.

Animals↗
PubMed2026

Advancing insect research through cell line transcriptomics.

This review emphasizes the significance of insect cell lines in transcriptomic research, highlighting their role as vital tools for uncovering cellular and molecular mechanisms of insect physiology, immune responses, and adaptation to environmental stressors. Cell lines derived from tissues such as the midgut, fat body, nervous system, and reproductive organs enable researchers to examine gene expression changes in a controlled setting, making discoveries that are difficult to achieve through whole-organism studies. High-throughput sequencing and single-cell RNA sequencing (scRNA-seq) have identified genes linked to detoxification, stress response, development, and immune defense, offering valuable insights for future applications in agriculture, pest control, and biotechnology. To organize this information clearly, we have summarized key findings in a table, providing an accessible overview of each cell line's important roles in transcriptomic research. This method not only highlights the adaptability of insect cell lines in functional genomics but also underscores their usefulness as model systems in pest management, virology, and bioengineering. Through utilizing transcriptomics, insect cell lines continue to advance our understanding of insect biology and foster the development of innovative strategies for sustainable crop protection and biotechnological use.

Animals↗
PubMed2026

Placenta-derived Exosomes Mitigate Hypoxia-Induced Trophoblast Apoptosis and Inflammatory Progression via SASH1.

SASH1 is a signal adaptor protein involved in cell growth, apoptosis, and immune regulation, and has been increasingly studied in tumor and immune cells. Emerging evidence suggests that SASH1 plays an important role in inflammatory responses and cellular homeostasis, processes that are closely associated with the development of PE. This study aimed to determine whether SASH1 contributes to trophoblast apoptosis and inflammatory responses in PE and whether P-EXOS exerts protective effects through SASH1 regulation. In this study, three PE-related transcriptomic datasets (GSE75010, GSE10588, and GSE60438) were analyzed to identify shared differentially expressed genes (DEGs), followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Machine learning algorithms were further applied to screen key candidate genes, and single-cell RNA sequencing data were used to characterize cellular heterogeneity in placental tissue and to determine cell type-specific expression patterns. SASH1 was identified as a consensus candidate gene and was significantly upregulated in trophoblast cells from PE samples. In vitro, a hypoxia-treated HTR-8/SVneo trophoblast cell model was established, combined with SASH1 knockdown, SASH1 overexpression, and co-culture with P-EXOS. Functional experiments showed that knockdown of SASH1 significantly suppressed hypoxia-induced trophoblast apoptosis and reduced the secretion of pro-inflammatory cytokines, including IL-6, IL-1β, and TNF-α, whereas SASH1 overexpression promoted apoptosis and inflammatory responses. In addition, P-EXOS treatment markedly reduced SASH1 expression at both mRNA and protein levels and attenuated hypoxia-induced trophoblast injury, while SASH1 overexpression largely abolished these protective effects. Taken together, these findings indicate that SASH1 plays a critical role in trophoblast apoptosis and inflammatory responses in PE. P-EXOS may alleviate hypoxia-induced trophoblastic injury by suppressing SASH1 expression, providing new insights into the molecular mechanisms and potential therapeutic targets for PE.

Trophoblasts↗
PubMed2026

Multiancestry genome-wide association and multiomics analyses elucidate spatiocellular features of multiple sclerosis genetics.

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system characterized by demyelination disseminated in space and time. Here we performed a genome-wide association study (GWAS) using 688 MS cases and 205,199 controls from the Japanese population and identified significant associations in the major histocompatibility complex region and a population-specific risk variant in 11q24. Through cross-population GWAS meta-analyses using a total of 29,374 cases and 1,843,563 controls from 4 ancestral populations, we identified 22 novel susceptibility loci. Integration of GWAS and single-cell and single-nucleus RNA sequencing of peripheral blood mononuclear cells and subcortical lesions from patients with MS revealed enrichment of genetic risk factors for MS in CD4+ T helper cell lineage and regulatory T cells, as well as in endothelial cells. Furthermore, spatial transcriptomics of subcortical lesions demonstrated spatial and temporal heterogeneity in associations with MS genetic risk. Our study demonstrates the value of investigation of spatiocellular features of disease genetics across diverse populations and omics modalities.

Humans↗
PubMed2026

Unraveling lung cancer complexity: Spatial omics in tumor microenvironment characterization and precision medicine.

Heterogeneous tumor microenvironment (TME) in lung cancer plays a crucial role in disease progression and resistance to therapy. Despite advances in single-cell and bulk omics profiling, these methods often overlook spatial context, which is vital for understanding cell-cell interactions and regional heterogeneity. In recent years, spatial omics technologies-including spatial genomics, transcriptomics, proteomics, and metabolomics-have revolutionized the ability to map molecular landscapes while maintaining tissue architecture. These advancements have become essential components of next-generation lung cancer management. By providing unprecedented resolution in characterizing the lung cancer TME, spatial omics could reveal prognostic and predictive biomarkers and identify new therapeutic vulnerabilities. This review will provide the first critical evaluation of spatial multi-omics approaches for lung cancer prognosis. It will also assess various integration strategies for multi-omics data to explore the clinical translational potential of these tools for therapy selection and patient stratification. Therefore, a deeper understanding of spatial omics technologies and their application in lung cancer can significantly improve precision diagnostics and therapeutic decision-making.

Lung cancer↗
PubMed2026

Translation of scRNA-seq to a clinical blood test for infection diagnostics.

INTRODUCTION: Early and accurate triage of patients with febrile illness is crucial for appropriate treatment. While standard inflammatory biomarkers are often nonspecific, transcriptome analysis of peripheral blood has diagnostic potential. However, bulk gene expression data is often confounded by changes in cell count proportions, a more robust quantification of gene expression in specific single-cell types, such as monocytes, is required to serve as a reliable clinical biomarker. AREAS COVERED: Various methods to obtain single-cell-type gene expression results, including the gold standard of gene expression analysis after cell sorting and single-cell RNA sequencing, which are difficult to implement in the routine settings are discussed. Other method to interrogate gene expression of a single cell-type is needed. Finally, monocyte cell-type specific ratio-based biomarker (RBB, called Direct Leukocyte Single cell-type Transcript Abundance, or DIRECT LS-TA) which can estimate single cell-type (monocyte) specific gene expression without cell sorting is introduced. EXPERT OPINION: Traditional diagnostic test for differentiating infection has several limitations requiring breakthrough including turn-around time and cost. DIRECT LS-TA provides a reliable way to quantify monocyte-specific gene expression that strongly correlates with gold-standard methods. It is more affordable than single-cell RNA sequencing and can be readily implemented in clinical laboratories using widely available quantitative PCR or digital PCR machines.

Humans↗
PubMed2026

Comparative phylogenomics and transcriptional regulatory networks of AQPs, HSPs, and LEA proteins in salt-stressed Portulaca oleracea.

Soil salinization severely threatens global food security, necessitating systematic investigations of halophytes like Portulaca oleracea to decode the molecular mechanisms of environmental resilience. Utilizing an integrated framework of deep learning-based genome annotation (58,817 predicted genes; 96.5% BUSCO completeness), multi-tissue RNA-Seq, phylogenomics, and gene regulatory network (GRN) inference, the synergistic orchestration of 78 aquaporins (AQPs), 525 heat shock proteins (HSPs), and 119 late embryogenesis abundant (LEA) proteins was elucidated. The active transcriptome, encompassing 39,065 expressed loci, revealed a systemic growth-defense trade-off. Tissues displayed distinct adaptive mechanisms: leaves modulated intracellular water balance via specialized AQPs, whereas adult roots maintained proteostasis through robust HSP20/HSP70 induction. Phylogenomic clustering across 154 species demonstrated that salinity tolerance constitutes an evolutionary mosaic, identifying 81 halophyte-exclusive orthogroups and 1129 species-specific clusters. Comparative topology across six independent GRNs (4.2M-5.3 M edges) unmasked a highly modular transcriptional reprogramming strategy governed by a core apparatus of 22 stress-exclusive regulators, with functional enrichment heavily prioritizing protein dimerization and chromatin remodeling. Theoretically, the distinct convergence of Trihelix transcription factors with guard cell differentiation pathways offers a candidate transcriptomic framework to explain the plant's characteristic C4-CAM photosynthetic plasticity under severe osmotic pressure. Practically, these evolutionary blueprints and specific master switches transcend single-gene transgenic limitations. Utilizing these root-sustained and stress-inducible targets under localized promoters provides a naturally optimized, network-level precision engineering roadmap to transfer robust, compartmentalized halotolerance to sensitive glycophytic crops.

Gene Regulatory Networks↗
PubMed2026

Four-dimensional molecular mapping from a spatial snapshot reveals the dynamics of hair follicle organogenesis.

Understanding organ formation requires capturing molecular information simultaneously in three-dimensional (3D) space and across developmental time. To this end, we developed 3D DNase-Enhanced Expression Profiling (3DEEP), a tissue-clearing approach that removes genomic DNA to extend spatial transcriptomic profiling hundreds of microns into intact tissues. We applied 3DEEP to neonatal mouse skin, capturing hundreds of developing hair follicles across their organogenesis trajectory. Ordering follicles by molecularly inferred developmental age transformed this single spatial snapshot into a four-dimensional (3D + time) molecular map of organogenesis. This map revealed developmental dynamics spanning stem cell compartment stratification, emergence of new cell subtypes within the follicle, and cascading structural transformations leading to hair canal formation. Comparative analysis of Foxn1-deficient nude mice, a hairlessness model, revealed organ-wide changes in developmental dynamics, including delayed molecular progression, reduced coordination, and increased developmental instability, preceding overt structural defects. This work demonstrates how deep-tissue spatial transcriptomics can uncover hidden dynamics of organ formation.

Animals↗
PubMed2026

A transcription factor regulatory atlas for activity inference and perturbation prediction.

Inferring transcription factor (TF) activity from transcriptomes and predicting transcriptome-wide responses to TF perturbations remain challenging, in part because available TF-mRNA resources often face a trade-off between precision and coverage and typically lack signed regulatory information. Here, we present TFActProfiler, a TF-mRNA resource and computational framework that learns signed, quantitative TF-mRNA regulatory coefficients by integrating heterogeneous prior evidence (ChIP-based, motif-based, and curated TF-mRNA annotations) with large-scale bulk and single-cell RNA-seq atlases. TFActProfiler contains 2 606 176 signed TF-mRNA interactions and improves TF activity inference in TF knockdown benchmarks relative to widely used regulon resources while retaining broad TF and target coverage. In addition, because the same learned regulatory coefficients can be used to model downstream transcriptional effects, TFActProfiler enables prediction of transcriptome-wide gene expression responses to TF knockdown without training on task-matched perturbation data. When perturbation datasets are available, TFActProfiler can be further refined to achieve performance comparable to state-of-the-art machine-learning baselines. By providing a direction-aware representation of TF-mRNA regulation for both activity inference and perturbation-response modeling, TFActProfiler supports systematic dissection of gene regulatory programs across diverse cellular contexts.

Transcription Factors↗
PubMed2025

Charting Postnatal Heart Development Using In Vivo Single-Cell Functional Genomics.

The transition at birth, marked by increased circulatory demands and rapid growth, necessitates extensive remodeling of the heart's structure, function, and metabolism. This transformation requires precise spatial and temporal coordination among diverse cardiac cell types; central to this process is cardiomyocyte maturation, yet the regulatory mechanisms driving these changes remain poorly understood. Here, we present a temporal and spatial atlas of postnatal hearts by integrating single-nucleus transcriptomics with image-based spatial transcriptomics, which uncovers the dynamic regulatory networks of cardiomyocyte maturation. To functionally interrogate candidate regulators in vivo , we developed Probe-based Indel-detectable Perturb-seq (PIP-seq), a high-throughput platform that uses probe-based chemistry to directly capture sgRNA expression, perturbation status, and transcriptomic profiles at single-nucleus resolution. Applying PIP-seq to postnatal cardiac development identified 21 novel regulators of cardiomyocyte maturation, highlighting critical nodal points in this process. Our study establishes a high-resolution framework for dissecting postnatal heart development, underscoring the integrative and highly ordered roles of microenvironment and intercellular communication in cardiomyocyte maturation. Importantly, PIP-seq enables systematic, high-throughput exploration of gene function and networks underlying complex biological processes in their native in vivo context.

Journal Article↗
PubMed2026

Single-cell expression quantitative trait locus Mendelian randomization reveals immune cell-specific causal regulatory networks and actionable targets in polycystic ovary syndrome.

ObjectiveTo systematically investigate whether the pathogenesis of polycystic ovary syndrome (PCOS) is causally related to dysregulated gene expression in specific immune cell subsets, and to evaluate the potential of these causal genes as actionable drug targets.MethodsThis study employed a two-sample Mendelian randomization (MR) framework using publicly available genome-wide association study (GWAS) summary statistics. The participant data included 797 PCOS cases and 140,558 controls (no direct patient recruitment was involved). Instrumental variables were derived from high-resolution immune cell-specific single-cell expression quantitative trait locus (sc-eQTL) data (OneK1K project) across 14 immune cell types. Primary analyses utilized the inverse-variance weighted (IVW) method. Shared causal variants were validated using Bayesian colocalization. Phenome-wide association analysis (PheWAS), external transcriptomic dataset validation (GSE8157), and DrugBank database screening were conducted for pleiotropy assessment and drug repositioning.ResultsMR analysis revealed genome-wide significant causal associations for GLIPR1 in non-classical monocytes (Mono NC) and XBP1 in CD4+ effector memory T cells (CD4 ET) with PCOS risk. Higher GLIPR1 expression was associated with a decreased PCOS risk (OR = 0.669, P = 4.34×10-6), whereas higher XBP1 expression was associated with an increased risk (OR = 1.406, P = 9.53×10-8). Colocalization analysis confirmed that GLIPR1 shares a causal variant with PCOS (PP.H4 = 96.73%). PheWAS and external validation confirmed the safety profile and significant upregulation (P = 0.03) of GLIPR1. Drug repositioning identified SOT-107, a Phase III protein therapy drug, as a potential interacting agent for GLIPR1.ConclusionsThis sc-eQTL MR study reveals immune cell-specific causal regulatory networks in PCOS. GLIPR1 in non-classical monocytes represents a high-confidence protective target, while XBP1 provides suggestive evidence for immune-mediated pathogenesis. The candidate drug SOT-107 highlights theoretical repositioning opportunities, though rigorous preclinical validation remains required.

Female↗
PubMed2026

scRNA-seq and bulk RNA-seq reveal the characteristics of macrophage copper metabolism and establish a risk signature in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is a prevalent malignancy with an urgent need for improved prognostic stratification and treatment-response prediction. This study aimed to explore a macrophage copper metabolism-associated prognostic model and to investigate the relationship between this risk model and the tumor immune microenvironment. METHODS: The FindClusters function was used to analyze cell clusters, and CellChat and CellPhoneDB/LIANA were employed for cell-cell communication analysis. Copper metabolism-related genes were sourced from the MSigDB database. A prognostic risk model was established using least absolute shrinkage and selection operator (LASSO) analysis and multivariate Cox regression analysis, and a nomogram was constructed by integrating the prognostic model with clinicopathological factors. Additional analyses were performed to map the seven model genes in single-cell data, assess model uncertainty and robustness, evaluate macrophage/copper/cuproptosis-related transcriptional programs, and examine the correlations between risk score, immune infiltration and predicted drug sensitivity. RESULTS: Using single-cell RNA sequencing (scRNA-seq) data, we identified four macrophage subpopulations. Macrophages with high SPP1 expression showed close interaction with T cell populations and were associated with copper ion metabolism. By incorporating 141 copper metabolism-related genes and using The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort, we constructed a seven-gene risk prediction model. Additional single-cell mapping showed that the model genes were detectable in the HCC single-cell dataset and showed a macrophage-associated expression pattern. The model showed moderate prognostic discrimination in TCGA-LIHC, whereas its external performance was heterogeneous and remained evaluable across external cohorts, with performance varying among datasets. Immune and mechanism-related analyses suggested that the risk signature was associated with macrophage-related infiltration, copper metabolism and cuproptosis-related transcriptional programs. Drug sensitivity analysis nominated Daporinad as a computationally predicted candidate compound, supporting Daporinad as a pharmacogenomic candidate for follow-up investigation. CONCLUSIONS: By integrating scRNA-seq and bulk RNA sequencing (RNA-seq) data, we constructed a macrophage copper metabolism-associated prognostic signature for HCC. The risk score was associated with survival, immune microenvironment features and predicted drug response, providing a transcriptomic framework for risk stratification and therapeutic hypothesis generation.

Hepatocellular carcinoma (HCC)↗
PubMed2026

Moving Beyond Morphology to Multiplexed Molecular Imaging as the Next Frontier in Diagnostic Pathology.

Diagnostic pathology has long relied on the morphologic interpretation of hematoxylin and eosin-stained tissues to guide diagnosis and assess prognostic features. Although pathologists intuitively recognize spatial patterns and architectural organization, these assessments remain largely qualitative and difficult to quantify systematically. Immunohistochemistry and immunofluorescence have introduced molecular specificity but are limited in multiplexing capacity, whereas bulk genomic and transcriptomic assays provide high molecular depth but lose spatial context by averaging signals across heterogeneous cell populations. Recent advances in spatial proteomics-including mass spectrometry-based imaging and cyclic immunofluorescence-now enable multiplexed, single-cell protein analysis within intact tissue architecture. These technologies have revealed complex immune and stromal microenvironments, spatially organized biomarkers predictive of therapeutic response, and molecular gradients underlying disease progression. By integrating histologic and molecular information, spatial proteomics bridges traditional microscopy with high-dimensional omics, allowing quantitative, spatially resolved insights into tissue organization and disease mechanisms. This review summarizes recent developments in multiplexed spatial proteomics from both scientific and pathologic perspectives, highlighting how these technologies extend beyond morphology to quantify histologic patterns, refine biomarker discovery, and facilitate clinical translation. The review also examines translational challenges and barriers to clinical implementation, including costs, standardization requirements, and workflow integration.

Humans↗
PubMed2026

Multimodal computational framework resolves B cell maturation in autoimmunity and ageing.

Identification of the origin of pathogenic immune cells is crucial for therapeutic interventions and diagnosis but pseudotime methods struggle to trace immune cells accurately. Current trajectory inference methods for B cell development and response in health and disease either ignore or underutilize antigen receptor sequence information, limiting their ability to resolve developmental pathways, particularly for pathogenic populations. Widely used methods such as Monocle 3 reconstruct developmental paths from transcriptomic similarity alone, discarding the features from immune receptors. Dandelion has combined the immune receptor features with transcriptomics but it struggles to simulate the trajectory path of B cells. Here we present ClonoTrace, a computational framework that integrates BCR sequence features with transcriptomic trajectory inference through gated fusion of multimodal embeddings. In fetal B cell development and germinal centre development, ClonoTrace demonstrates closer concordance with the canonical reference ordering than Monocle 3 and Dandelion. Applied to systemic lupus erythematosus, ClonoTrace indicates a memory B cell extrafollicular maturation route alongside the naïve B cell route, accompanied by induction of ZEB2 with a concomitant decline of BACH2 along the trajectory, as a candidate alternative route to pathogenic double negative 2 B cells (DN2) in systemic lupus erythematosus (SLE) patients. In healthy ageing, ClonoTrace resolved three candidate age-related B cell maturation routes, from naïve, IgM+ memory and switched-memory B cells, each passing through a DN2-associated transcriptional state that is ordered before age-associated B cells along the inferred trajectory. ClonoTrace's fate probability algorithm indicated that IgM+ memory B cell to ABC transition as the leading candidate age-associated transition, which may be distinct from SLE DN2 maturation. ClonoTrace provides a generalizable framework for receptor-informed trajectory inference, describing candidate developmental routes of pathogenic B cell populations in autoimmunity and ageing.

Humans↗
PubMed2023

GBMdeconvoluteR accurately infers proportions of neoplastic and immune cell populations from bulk glioblastoma transcriptomics data.

BACKGROUND: Characterizing and quantifying cell types within glioblastoma (GBM) tumors at scale will facilitate a better understanding of the association between the cellular landscape and tumor phenotypes or clinical correlates. We aimed to develop a tool that deconvolutes immune and neoplastic cells within the GBM tumor microenvironment from bulk RNA sequencing data. METHODS: We developed an IDH wild-type (IDHwt) GBM-specific single immune cell reference consisting of B cells, T-cells, NK-cells, microglia, tumor associated macrophages, monocytes, mast and DC cells. We used this alongside an existing neoplastic single cell-type reference for astrocyte-like, oligodendrocyte- and neuronal progenitor-like and mesenchymal GBM cancer cells to create both marker and gene signature matrix-based deconvolution tools. We applied single-cell resolution imaging mass cytometry (IMC) to ten IDHwt GBM samples, five paired primary and recurrent tumors, to determine which deconvolution approach performed best. RESULTS: Marker-based deconvolution using GBM-tissue specific markers was most accurate for both immune cells and cancer cells, so we packaged this approach as GBMdeconvoluteR. We applied GBMdeconvoluteR to bulk GBM RNAseq data from The Cancer Genome Atlas and recapitulated recent findings from multi-omics single cell studies with regards associations between mesenchymal GBM cancer cells and both lymphoid and myeloid cells. Furthermore, we expanded upon this to show that these associations are stronger in patients with worse prognosis. CONCLUSIONS: GBMdeconvoluteR accurately quantifies immune and neoplastic cell proportions in IDHwt GBM bulk RNA sequencing data and is accessible here: https://gbmdeconvoluter.leeds.ac.uk.

Humans↗
PubMed2026

SpatialRNA: a Python package for easy application of Graph Neural Network models on single-molecule spatial transcriptomics dataset.

SUMMARY: Image-based spatial transcriptomics (iST) deliver gene expression measurements of RNA transcripts in tissue slices with single-molecule resolution and spatial context preserved. Modern Graph Neural Network (GNN) models are promising methods for capturing the complex molecular and cellular phenotypes in tissues at single-transcript and single-cell levels. A key application of GNNs is the detection of spatial domains or niches, that is, groups of molecules and/or cells that collaboratively work together to produce complex phenotypes. Due to the vast number of detected transcripts in (iST) dataset, applying GNNs on RNA molecule graphs is not trivial. We present a Python package, SpatialRNA, for easy (sub)graph generation from tissue samples and provide comprehensive tutorials for convenient and efficient application of Graph Neural Network models under the PyG framework. This highly scalable tool comprehensively segments tissue into spatial domains, aiding in biological interpretation of iST data and its underlying molecular microenvironments. AVAILABILITY AND IMPLEMENTATION: The SpatialRNA package is freely accessible from online repository https://github.com/ruqianl/spatialrna and can be installed via pip. Comprehensive tutorials, guidance on parameter selection, and complete workflows of case studies are available from the documentation website https://ruqianl.github.io/spatialrna_docs/, and uploaded on Zenodo with a DOI 10.5281/zenodo.17339575.

Neural Networks, Computer↗