SEARCH · PubMed Health

Results for “Single cell transcriptomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 343 records · Page 19Linked to original sources

PubMed2026

Cellular transcriptomic signatures underpinning the heterogeneity of depression in Alzheimer's disease.

INTRODUCTION: Late-onset Alzheimer's disease (LOAD) and major depressive disorder (MDD) share genetic etiologies. Here, we investigated brain transcriptomic landscapes to gain insights into shared and divergent molecular and biological etiologies across LOAD and MDD. METHODS: Brain single-nucleus RNA sequencing (snRNA-seq) datasets from cognitively normal older and young individuals and LOAD patients stratified by comorbid MDD were analyzed to identify differential expressed genes (DEGs). Using cell type-specific DEGs we performed biological pathway and intercellular-communication networks analyses. We investigated shared DEGs across MDD and LOAD cohorts and sex-specific DEGs. Results were validated by comparison with four transcriptomic and proteomic studies of MDD and depression. RESULTS: MDD-associated dysregulated genes and pathways were shared between LOAD and cognitive-normal individuals, including JUNB and DUSP1 in glutamatergic neurons, and PRAM1 and SNX9 in microglia. DEGs shared between the MDD and LOAD cohorts included HSPA1A and NDUFB7 in glutamatergic neurons. Sex interaction analysis identified numerous new DEGs in the MDD cohorts, whereas there were ≈5 to 10 times more DEGs in female than in male individuals. LOAD and MDD common microglial pathways included neuronal injury, stress, peroxisome proliferator-activated receptor (PPAR) signaling and interferon alpha/beta signaling. DISCUSSION: LOAD and MDD exhibited common molecular profiles, dysregulated pathways, and cellular communication changes. MDD develops earlier in life, thus, our findings provide a window into early molecular and biological processes preceding LOAD-onset.

Humans↗

PubMed2026

The pseudouridine epitranscriptomic landscape of advanced prostate cancer therapeutic resistance identifies TIMM17A as a key player.

BACKGROUND: Resistance to androgen receptor signaling inhibitors (ARSIs) remains a major barrier of advanced prostate cancer (PCa) treatment. While RNA epitranscriptomic modifications are increasingly recognized as key regulators of tumor biology, the role of pseudouridine (Ψ) in therapeutic resistance is largely unexplored. METHODS: A darolutamide-resistant PCa cell model was established and subjected to integrated multi-omics profiling using bulk RNA sequencing and photo-crosslinking-assisted Ψ sequencing (PA-Ψ-seq). Differential expression and pseudouridylation analyses were combined to identify Ψ-associated genes. Public datasets validated expression and prognosis. Functional assays including RNA knockdown, cell proliferation, colony formation, and xenograft models were conducted. Single-cell RNA sequencing investigated tumor microenvironment (TME) interactions. RESULTS: We identified extensive transcriptomic and pseudouridylation alterations associated with ARSI resistance, with a significant positive correlation between Ψ modification and mRNA expression. Integrated analysis highlighted a subset of "hyper-up" genes enriched in resistance-related pathways. Thus, TIMM17A was identified as a novel candidate. TIMM17A expression was significantly elevated in PCa and correlated with disease progression and poor prognosis. Experimental validations demonstrated that TIMM17A promoted tumor growth and resistance, while its knockdown restored sensitivity to darolutamide both in vitro and in vivo. Mechanistically, TIMM17A expression may be regulated by PUS1‑mediated pseudouridylation. Single-cell analysis further revealed that TIMM17A is enriched in malignant epithelial cells and associated with enhanced cell-cell communication within the TME. CONCLUSIONS: This study delineates the pseudouridine epitranscriptomic landscape in advanced PCa and identifies TIMM17A as a key mediator of therapeutic resistance. Targeting the Ψ-TIMM17A axis may offer a novel strategy to overcome ARSI resistance.

Advanced prostate cancer↗

PubMed2026

Genomic Structural Equation Modeling Identifies a Shared Inflammatory Genetic Dimension Across Inflammatory Arthritis Phenotypes and Biomarkers.

BACKGROUND: Inflammatory arthritis (IA), including rheumatoid arthritis (RA), psoriatic arthritis (PsA) and gout, shares systemic inflammatory features indexed by C-reactive protein (CRP) and interleukin-6 (IL-6), yet the extent of their common genetic basis remains unclear. AIMS: We aimed to delineate the shared genetic architecture across IA phenotypes and inflammatory biomarkers. MATERIALS AND METHODS: We applied genomic structural equation modelling (Genomic SEM) to GWAS summary statistics for RA, PsA, gout, CRP and IL-6, fitted a single common factor, and performed multivariate GWAS followed by fine-mapping, transcriptome-wide association, gene-based analysis, pathway enrichment, and cell-type and spatial mapping. RESULTS: A single common factor was fitted (CFI = 0.990, SRMR = 0.045). The multivariate GWAS identified 56 genome-wide significant SNPs across 10 independent lead loci, including one novel signal. Fine-mapping prioritized high-confidence variants near PTPN22, the CRP gene cluster and a urate-associated locus. Gene-level analyses converged on DCLRE1B, PTPN22, IL6R, NLRP3 and HNF1A, with pathway enrichment implicating inflammasome assembly and metabolic-inflammatory overlap. Cell-type enrichment highlighted myeloid populations, and spatial mapping localized signals to lung, kidney, mucosal epithelium and gastrointestinal tissues. DISCUSSION: These results delineate a shared inflammatory genetic dimension across IA phenotypes and biomarkers, anchored in immune, inflammasome, cytokine-receptor and metabolic pathways. CONCLUSION: Together, these findings provide a valuable framework for prioritizing candidate genes and cellular contexts for future investigation.

TWAS↗

PubMed2026

FAP+ pericyte-like cells promote monocyte differentiation into tumor-associated macrophages in glioblastoma.

Glioblastoma (GBM) is a highly aggressive primary brain tumor characterized by profound immunosuppression that facilitates tumor progression and promotes therapeutic resistance. Fibroblast activation protein (FAP), a recognized theranostic target in multiple cancers, is upregulated in GBM and predominantly expressed by pericyte-like stromal cells. Here we identify a role for FAP⁺ pericyte-like cells in shaping the GBM immune microenvironment through monocyte recruitment and differentiation. Analysis of The Cancer Genome Atlas (TCGA) datasets, supported by reverse-transcription quantitative PCR and immunohistochemistry, revealed that elevated FAP expression-serving as a proxy for the abundance of FAP⁺ pericyte-like cells-is associated with an immune-enriched tumor microenvironment characterized by higher macrophage abundance and elevated expression of M2 polarization markers. Spatial analyses, including immunofluorescence and spatial transcriptomics, demonstrated that immunosuppressive macrophages preferentially localize in proximity to FAP⁺ pericytes. Single-cell RNA sequencing identified these FAP⁺ cells as a distinct perivascular stromal subset with a unique expression pattern of extracellular matrix components and cytokines, including CCL2 and CSF1, with corresponding receptors expressed on myeloid cells. Functional assays using patient-derived FAP⁺ pericyte-like cells confirmed their ability to attract monocytes via soluble mediators and to promote their differentiation and polarization into tumor-associated macrophages with immunoregulatory features, partly mediated by the CSF1-CSF1R axis. Orthotopic co-implantation experiments in mice further supported their capacity to enhance myeloid infiltration in vivo. Consistent with these biological effects, a transcriptional signature characteristic of FAP⁺ pericytes correlated with worse overall survival in patients with GBM. Together, these findings position FAP⁺ pericyte-like cells as modulators of the GBM immune landscape, fostering a tumor-permissive niche by promoting the differentiation of circulating monocytes into immunoregulatory macrophages. Targeting this stromal population may offer new therapeutic avenues to reprogram tumor-associated immune responses in GBM.

Journal Article↗

PubMed2026

SCMO: a deep learning model integrating the single-cell resolution TME ecosystem and multi-omics for survival prediction in CRC patients.

BACKGROUND: Colorectal cancer (CRC) remains a leading cause of global cancer mortality, highlighting the need for precise survival prediction to guide clinical decisions. Although tissue-level multi-omics is widely utilized for survival prediction, its limited resolution cannot capture tumor heterogeneity. Single-cell RNA sequencing (scRNA-seq) enables dissection of the tumor microenvironment (TME) at cellular resolution, supporting personalized prognostic assessment. METHODS: We collected 213 CRC scRNA-seq samples and established a CRC-specific TME atlas comprising 339,060 cells. Using this atlas as a reference, we deconvolved bulk RNA-seq data from TCGA-CRC cohort with the EcoTyper algorithm to reconstruct TME features. Clinical, genomic, and transcriptomic data were obtained from the Xena platform; microbial data were sourced from the BIC database. We integrated TME and multi-omics features through a self-normalizing neural network to construct a deep learning model (single-cell resolution TME ecosystem with multi-omics data [SCMO]) for survival prediction. To enhance interpretability, we utilized the Integrated Gradients algorithm and spatial transcriptomic data to analyze multi-omics and TME features. We performed anticancer drug screening with tumor necrosis factor receptor-associated protein 1 (TRAP1), a critical feature according to the Integrated Gradients algorithm, as a potential target. RESULTS: We identified 13 survival-related TME features from the CRC-specific atlas: 12 cell states and one multi-cellular ecosystem. SCMO, which combined TME and multi-omics features, improved survival prediction and outperformed existing methods, achieving a concordance index of 0.762. The SCMO demonstrated robust performance for long-term predictions, achieving areas under the curve (AUCs) of 0.752, 0.772, and 0.869 for 1-, 3-, and 5-year predictions in the training set, with corresponding test set AUCs of 0.639, 0.756, and 0.772. TME features from the SCMO model revealed that ecosystem density increased with CRC malignancy. Multi-omics features included TRAP1 as a potential drug target. Drug screening identified saikosaponin A as a novel TRAP1 inhibitor, and its anticancer activity was validated in vitro. We developed SCMO-Lite, a simplified model incorporating 12 high-attribution-weight multi-omics features, which demonstrated robust risk stratification. CONCLUSIONS: SCMO combines analytical precision with biological interpretability, offering novel insights for oncology survival prediction.

Humans↗

PubMed2025

The offonome reveals on and off states of gene expression near the detection limit of RNA-seq.

RNA-seq, widely used for gene expression profiling, provides nucleotide level genome coverage and summary gene expression values. Generally, low-expressed genes are ignored due to their unfavorable signal-to-noise ratio, however, these genes may offer crucial information, such as detecting rare cells in bulk tissues. In this study, we applied an approach that transforms the expression levels of low-expressed genes into a robust dichotomized on/off state by leveraging similarities in transcript coverage shape. Applied to three human cancer cohorts from the Cancer Genome Atlas (TCGA), chosen based on tissue morphology and anatomic site, we identified genes, the "offonome" near the detection limit, consistently or occasionally off across samples. Genes in the offonome spectrum proved useful for supervised and unsupervised applications, including characterizing oncogenic pathways, and identifying rare populations of cells in bulk tissue. Interrogating the offonome is relevant to bulk tumor analyses like TCGA, potentially expediting gene investigation in low-input situations like single cell RNA-seq.

Humans↗

PubMed2026

Multi-omics characterization of a GPRC5A+ epithelial subpopulation associated with malignant features in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) exhibits marked cellular heterogeneity, and the cellular context of malignancy-associated epithelial programs remains incompletely defined. METHODS: We integrated 2,993 CRC samples spanning bulk RNA-seq (n = 2,568; two OS/RFS cohorts), scRNA-seq (281,961 cells/152 specimens), spatial transcriptomics (n = 6), and proteomics (n = 267). Analyses included single-cell integration/annotation, GSVA/HALLMARK, interactome, pseudotime, and ligand-receptor mapping; functional CRISPR assays, EMT immunoblotting, and xenografts; TF profiling (SCENIC/JASPAR/ChIP-qPCR); and exploratory drug-response prediction (OncoPredict), cell-sensitivity assays, and docking/MD modeling. RESULTS: We constructed a stage-stratified single-cell atlas and resolved eleven malignant epithelial subsets, characterizing Epi_4 as late-stage-enriched with EMT, hypoxia, and inflammatory programs and adverse OS/RFS. GPRC5A marked this subset, which we define as GPRC5A+Epi; its expression rose from stage I→IV and was associated with poor outcomes across cohorts, with concordant spatial/proteomic observations. GPRC5A perturbation affected CRC proliferation, migration/invasion, EMT, and xenograft tumorigenicity, supporting a functionally important role in the tested models. SCENIC and ChIP-qPCR supported FOSL1 as an upstream regulator that occupies the GPRC5A promoter. Spatial and ligand-receptor analyses predicted close association and potentially reciprocal signaling between GPRC5A+Epi and POSTN+fibroblasts (COL1A1-SDC4, COL1A1/1A2-ITGA2/ITGB1, PPIA-BSG); concurrent high GPRC5A+Epi/POSTN+Fib signatures were associated with inferior OS/RFS. Drug-response analyses identified an association between GPRC5A status and trametinib sensitivity. Docking/MD produced a computational model of a possible trametinib-GPRC5A interaction, which remains experimentally unvalidated. CONCLUSIONS: GPRC5A⁺Epi is a malignancy-associated epithelial state in CRC, and GPRC5A is functionally important for malignant phenotypes in the tested models. Its inferred relationships with POSTN⁺ fibroblasts and the trametinib findings should be regarded as hypothesis-generating pending functional crosstalk, direct-binding, and therapeutic validation.

Humans↗

PubMed2025

Droplet-Based Single-Cell 3' mRNA Sequencing of Marburg Virus-Infected Samples.

Single-cell technologies are continually evolving with emerging methods that are gradually uncovering the central DNA-RNA-protein dogma. Single-cell RNA sequencing is one arm of a multi-omic approach that achieves an astounding level of granularity to reveal the complexity of virus-host interactions at the transcriptomic level. Cell tropism, virus replication, pathogenesis, and gene expression changes mediated by the virus and the host's immune response to infection are just some areas of study that are gaining better clarity due to the high-resolution analysis afforded by the technology.We describe a single-cell sequencing protocol for Marburg virus infection in vivo using nonhuman primate blood and the 10× Chromium Next GEM single-cell genomics methodology. Working with pathogens of high consequence is logistically complicated, requiring containment in biosafety level (BSL)-4 laboratories and harsh inactivation procedures before samples can safely be removed to lower biosafety conditions. We provide procedural insight into sample isolation and processing conducted in BSL-4 and describe the requirements for safe sample removal without jeopardizing quality for down-stream sequencing and analysis in BSL-2 conditions. Characterization of complicated biological processes mediated by high-containment pathogens, typically restricted to analogous model systems, e.g., minigenome, can be achieved using live virus.

Animals↗

PubMed2025

scnanoseq: an nf-core pipeline for Oxford Nanopore single-cell RNA-sequencing.

MOTIVATION: Recent advancements in long-read single-cell RNA sequencing (scRNA-seq) have facilitated the quantification of full-length transcripts and isoforms at the single-cell level. Historically, long-read data would need to be complemented with short-read single-cell data in order to overcome the higher sequencing errors to correctly identify cellular barcodes and unique molecular identifiers. Improvements in Oxford Nanopore sequencing, and development of novel computational methods have removed this requirement. Though these methods now exist, the limited availability of modular and portable workflows remains a challenge. RESULTS: Here, we present, nf-core/scnanoseq, a secondary analysis pipeline for long-read single-cell and single-nuclei RNA that delivers gene and transcript-level quantification. The scnanoseq pipeline is implemented using Nextflow and is built upon the nf-core framework, enabling portability across computational environments, scalability and reproducibility of results across pipeline runs. The nf-core/scnanoseq workflow follows best practices for analyzing single-cell and single-nuclei data, performing barcode detection and correction, genome and transcriptome read alignment, unique molecular identifier deduplication, gene and transcript quantification, and extensive quality control reporting. AVAILABILITY AND IMPLEMENTATION: The source code, and detailed documentation are freely available at https://github.com/nf-core/scnanoseq and https://nf-co.re/scnanoseq under the MIT License. Documentation for the version of nf-core/scnanoseq used for this paper, including default parameters and descriptions of output files are available at https://nf-co.re/scnanoseq/1.1.0.

Single-Cell Analysis↗

PubMed2026

Multiomics approaches to cardiovascular disease: technological innovations and clinical translation.

Cardiovascular diseases (CVDs) remain the leading cause of global morbidity and mortality, reflecting a persistent gap between clinical phenotyping and the molecular mechanisms that govern disease initiation, progression, and interindividual variability. Recent advances in emerging technologies have fundamentally reshaped cardiovascular physiology by enabling high-resolution, cross-layer profiling of the heart and vasculature across genomic, epigenomic, transcriptomic, proteomic, metabolomic, lipidomic, glycomic, and fluxomic layers, increasingly at single-cell and spatial resolution. These approaches reveal CVD as a coordinated, multilayered process driven by dynamic interactions among cell types, regulatory programs, and metabolic states, rather than isolated gene-level defects. In this review, we synthesize how emerging multiomic, computational, and functional genomic technologies are redefining the study of cardiovascular disease across molecular, cellular, and tissue levels. We highlight recent innovations in single-cell and spatial atlases, long-read sequencing, proteomics and metabolomics, integrative data modeling, and functional omics approaches, including genome-scale perturbation screens and single-cell perturbation frameworks. These platforms enable mechanistic dissection of regulatory circuits, distinguish primary disease drivers from secondary adaptations, and directly assess therapeutic reversibility, advancing the field beyond associative biomarker discovery toward mechanism-guided target prioritization. We further discuss key methodological and translational challenges accompanying high-dimensional cardiovascular data, including preanalytical variability, control selection, temporal misalignment across molecular layers, population diversity, and reference bias. By integrating technological innovation with computational rigor and functional validation, this review frames emerging omics-enabled strategies as a unified, physiologically grounded framework for translating molecular insight into clinically meaningful cardiovascular phenotypes and advancing precision cardiovascular medicine.

Humans↗

PubMed2005

The maize root transcriptome by serial analysis of gene expression.

Serial Analysis of Gene Expression was used to define number and relative abundance of transcripts in the root tip of well-watered maize seedlings (Zea mays cv FR697). In total, 161,320 tags represented a minimum of 14,850 genes, based on at least two tags detected per transcript. The root transcriptome has been sampled to an estimated copy number of approximately five transcripts per cell. An extrapolation from the data and testing of single-tag identifiers by reverse transcription-PCR indicated that the maize root transcriptome should amount to at least 22,000 expressed genes. Frequency ranged from low copy number (2-5, 68.8%) to highly abundant transcripts (100-->1,200; 1%). Quantitative reverse transcription-PCR for selected transcripts indicated high correlation with tag frequency. Computational analysis compared this set with known maize transcripts and other root transcriptome models. Among the 14,850 tags, 7,010 (47%) were found for which no maize cDNA or gene model existed. Comparing the maize root transcriptome with that in other plants indicated that highly expressed transcripts differed substantially; less than 5% of the most abundant transcripts were shared between maize and Arabidopsis (Arabidopsis thaliana). Transcript categories highlight functions of the maize root tip. Significant variation in abundance characterizes transcripts derived from isoforms of individual enzymes in biochemical pathways.

Base Sequence↗

PubMed2026

scAmp enables focal gene amplification analysis from single-cell data.

Oncogene amplification on extrachromosomal DNA is a common driver of tumor progression and is associated with acquired drug resistance and poor patient survival. While bulk whole genome sequencing studies have revealed the landscape of genes amplified on extrachromosomal DNA in tumors, it remains challenging to study the subclonal heterogeneity and functional (e.g., transcriptomic) consequences of extrachromosomal DNA on tumors. To address this, we introduce scAmp: a probabilistic algorithm for detecting and analyzing extrachromosomal DNA from single-cell datasets. Using well-characterized cell lines, we demonstrate that scAmp has improved specificity over bulk genome sequencing in predicting extrachromosomal DNA status and can resolve the status of chromosomal amplifications that were historically extrachromosomal. We further showcase scAmp by analyzing 73 patient tumors profiled with single-cell assay for transposase-accessible chromatin by sequencing, where we characterize the subclonal evolution of subclones with extrachromosomal DNA and identify the effect of these amplifications on the chromatin accessibility landscape of cancer cells. Finally, we provide proof-of-concept analyses that scAmp aids in the detection of extrachromosomal DNA from clinical histopathology assays. Together, we anticipate that scAmp will broadly enable further studies - both retrospective and prospective - that dissect critical questions of how extrachromosomal DNAs affect cancer cells and the tumors in which they reside.

Humans↗

PubMed2025

CAFs activated by YAP1 upregulate cancer matrix stiffness to mediate hepatocellular carcinoma progression.

BACKGROUND: The stiffness of the matrix is closely related to the progression of hepatocellular carcinoma (HCC). Although direct targeting of stromal rigidity in HCC remains a clinical challenge, cancer-associated fibroblasts (CAFs) are considered key contributors to this process. Given the heterogeneity of CAFs, this study explored the relationship between specific CAF subsets and liver cancer matrix stiffness, aiming to identify novel therapeutic targets for HCC patients. METHODS: Single-cell sequencing datasets were leveraged to identify cell types within liver cancer and characterize the transcriptomic profiles of CAFs. Prognostic analysis, utilizing the Gene Expression Profiling Interactive Analysis (GEPIA) and The Cancer Genome Atlas (TCGA) liver cancer datasets, assessed the correlation between matrix stiffness-related genes and HCC patient outcomes. Pseudo-time analysis was applied to trace the developmental trajectories of CAFs. By calculating intercellular communication probabilities and analyzing transcription factor activity, the functions and interactions of different CAF subsets were elucidated. Gene Ontology (GO) analysis was used to explore the functional roles of CAFs in distinct Yes-associated protein 1 (YAP1) groups. Finally, cellular experiments and animal experiments were further conducted to validate the hypotheses of this study. RESULTS: This study identified CAF subpopulations based on single-cell sequencing data and analyzed transcriptional changes within these subpopulations. Key findings include the identification of collagen type I alpha 1 (COL1A1), collagen type III alpha 1 (COL3A1), and lysyloxidase (LOX) as pivotal node genes during CAF development. Moreover, the expression of matrix stiffness-related genes was inversely correlated with the prognosis of HCC patients. Notably, the YAP1-positive CAF subpopulation emerged as the primary contributor to matrix stiffness in liver cancer. This subpopulation upregulates the expression of matrix stiffness-related genes and promotes tumor progression by activating signaling pathways such as autophagy and GTPase activity regulation. Cellular experiments and animal studies further validated this conclusion. CONCLUSION: This single-cell analysis uncovered the functional roles of CAFs in liver cancer. The YAP1-positive CAF subpopulation, in particular, was shown to contribute to matrix stiffness by upregulating the expression of relevant genes and promoting tumor progression through the activation of specific signaling pathways.

Carcinoma, Hepatocellular↗

PubMed2024

Nonlethal deleterious mutation-induced stress accelerates bacterial aging.

Random mutagenesis, including when it leads to loss of gene function, is a key mechanism enabling microorganisms' long-term adaptation to new environments. However, loss-of-function mutations are often deleterious, triggering, in turn, cellular stress and complex homeostatic stress responses, called "allostasis," to promote cell survival. Here, we characterize the differential impacts of 65 nonlethal, deleterious single-gene deletions on Escherichia coli growth in three different growth environments. Further assessments of select mutants, namely, those bearing single adenosine triphosphate (ATP) synthase subunit deletions, reveal that mutants display reorganized transcriptome profiles that reflect both the environment and the specific gene deletion. We also find that ATP synthase α-subunit deleted (ΔatpA) cells exhibit elevated metabolic rates while having slower growth compared to wild-type (wt) E. coli cells. At the single-cell level, compared to wt cells, individual ΔatpA cells display near normal proliferation profiles but enter a postreplicative state earlier and exhibit a distinct senescence phenotype. These results highlight the complex interplay between genomic diversity, adaptation, and stress response and uncover an "aging cost" to individual bacterial cells for maintaining population-level resilience to environmental and genetic stress; they also suggest potential bacteriostatic antibiotic targets and -as select human genetic diseases display highly similar phenotypes, - a bacterial origin of some human diseases.

Escherichia coli↗

PubMed2026

Integrative multi-omics and single-cell analysis identifies EGFR pathway activation and metabolic reprogramming as potential synthetic lethal vulnerabilities in resistance to the FGFR inhibitor AZD4547.

BACKGROUND: Although fibroblast growth factor receptor (FGFR) inhibitors (FGFRi) have demonstrated clinical promise, the inevitable emergence of acquired resistance remains a critical bottleneck, severely compromising their long-term clinical efficacy. The pan-cancer molecular landscape and heterogeneous mechanisms driving this resistance, ranging from genetic alterations to dynamic network rewiring, remain poorly understood. METHODS: We integrated large-scale pharmacogenomic profiling of the FGFR inhibitor AZD4547 from the GDSC2 and PRISM databases with single-cell RNA sequencing to dissect the multi-omics landscape of FGFRi resistance across 312 cell lines from 8 cancer types. This multi-omics framework was further extended by machine learning modeling and systematic synthetic lethality screening to uncover actionable therapeutic targets. In vitro viability assays and western blot analysis were subsequently conducted to experimentally evaluate the predicted FGFR-EGFR synthetic lethality. RESULTS: Our dual-database analysis unveiled a multi-dimensional atlas of FGFRi resistance. We identified cancer-specific genomic drivers, such as ELF4 amplification in glioblastoma, alongside key transcriptomic markers including UCP2 and FSCN1, highlighting a shift towards metabolic reprogramming and epithelial-mesenchymal transition (EMT). Single-cell analysis unveiled that resistance is linked to the heterogeneous enrichment of baseline subpopulations characterized by distinct metaprograms, including cell-cycle dysregulation. Furthermore, a random forest model built on a LASSO-derived transcriptomic signature was constructed, demonstrating promising predictive capability for AZD4547 sensitivity (mean test-set AUC = 0.73, 95% CI [0.63, 0.80]); the signature generalized well to erdafitinib but showed limited transferability to some other FGFR inhibitors (e.g. pemigatinib, BGJ398). Most notably, our synthetic lethal screening revealed a convergent reliance on compensatory RTK signaling (specifically EGFR pathway enrichment) and downstream MAPK/PI3K cascades in resistant phenotypes, providing converging computational evidence for EGFR pathway activation as an adaptive bypass mechanism. This predicted synthetic lethality was experimentally supported in two FGFR-dependent cell line models (RT112 and CCLP1), in which combined FGFR-EGFR inhibition produced marked synergistic antiproliferative effects. CONCLUSIONS: This study establishes a comprehensive multi-omics atlas of resistance to the FGFR inhibitor AZD4547, delineating convergent mechanisms of metabolic reprogramming and EGFR-mediated bypass signaling. Our findings characterize the resistance as a dynamic network rewiring and nominate rational combination strategies to overcome this therapeutic bottleneck. While FGFR-EGFR co-inhibition is experimentally supported, metabolic co-targeting remains a computationally derived, hypothesis-generating strategy.

Benzamides↗

PubMed2025

Single-cell RNA sequencing defines developmental progression and reproductive transitions of Pneumocystis carinii.

UNLABELLED: Pneumocystis species are host-obligate fungal pathogens that cause severe pneumonia in immunocompromised individuals. Despite their clinical importance, their life cycle remains poorly understood, in part because Pneumocystis depends on the host environment for most nutrients and requires sexual reproduction for survival, which occurs exclusively in vivo. This study presents the first single-cell RNA sequencing (scRNA-seq) atlas of Pneumocystis carinii, generated from isolated organisms recovered from the bronchoalveolar lavage fluid of infected rats to map the life cycle of P. carinii. Transcriptomes from 87,716 cells were analyzed using the 10× Genomics platform, revealing 13 transcriptionally distinct clusters representing key developmental stages, including biosynthetically active trophic forms, mating-competent intermediates, and asci undergoing sporulation. These states were characterized by expression of MAPK signaling components, β-glucan-modifying enzymes, and spore-associated genes, respectively. The scRNA-seq data support previous evidence that these host-obligate fungi undergo sexual reproduction and provide new insights into the gene expression patterns associated with different life cycle phases. Biomarkers associated with ascus formation identified by scRNA-seq were validated by RT-qPCR, showing decreased expression levels in ascus-depleted populations treated with anidulafungin, a drug that halts ascus formation. More broadly, this approach provides a strategy for studying the full life cycles of fungal pathogens that cannot be continuously cultured. IMPORTANCE: Pneumocystis species (spp.) are clinically significant fungal pathogens that cannot be sustainably cultured in vitro due to their host-obligate nature. This longstanding limitation has impeded progress in understanding their life cycle and identifying therapeutic vulnerabilities. Here, we apply scRNA-seq to P. carinii isolated directly from infected rat lungs, generating the first transcriptional map of its developmental progression. Our results define discrete gene expression states associated with trophic growth, mating activation, and ascus formation and provide transcriptional evidence for a structured life cycle, clarifying key developmental transitions and identifying potential regulatory targets for therapeutic intervention. Importantly, this study demonstrates that scRNA-seq can resolve the developmental biology of host-restricted fungal pathogens that cannot be cultured in vitro. This approach offers a generalizable framework for investigating other unculturable or obligate microbial pathogens directly within their native host environments, where traditional experimental tools are limited.

Pneumocystis carinii↗

PubMed2025

Alevin-fry-atac enables rapid and memory frugal mapping of single-cell ATAC-seq data using virtual colors for accurate genomic pseudoalignment.

Ultrafast mapping of short reads via lightweight mapping techniques such as pseudoalignment has significantly accelerated transcriptomic and metagenomic analyses, often with minimal accuracy loss compared to alignment-based methods. However, applying pseudoalignment to large genomic references, like chromosomes, is challenging due to their size and repetitive sequences. We introduce a new and modified pseudoalignment scheme that partitions each reference into "virtual colors…. These are essentially overlapping bins of fixed maximal extent on the reference sequences that are treated as distinct "colors" from the perspective of the pseudoalignment algorithm. We apply this modified pseudoalignment procedure to process and map single-cell ATAC-seq data in our new tool alevin-fry-atac . We compare alevin-fry-atac to both Chromap and Cell Ranger ATAC . Alevin-fry-atac is highly scalable and, when using 32 threads, is approximately 2.8 times faster than Chromap (the second fastest approach) while using approximately one third of the memory and mapping slightly more reads. The resulting peaks and clusters generated from alevin-fry-atac show high concordance with those obtained from both Chromap and the Cell Ranger ATAC pipeline, demonstrating that virtual colorenhanced pseudoalignment directly to the genome provides a fast, memory-frugal, and accurate alternative to existing approaches for single-cell ATAC-seq processing. The development of alevin-fry-atac brings single-cell ATAC-seq processing into a unified ecosystem with single-cell RNA-seq processing (via alevin-fry ) to work toward providing a truly open alternative to many of the varied capabilities of CellRanger . Furthermore, our modified pseudoalignment approach should be easily applicable and extendable to other genome-centric mapping-based tasks and modalities such as standard DNA-seq, DNase-seq, Chip-seq and Hi-C.

Journal Article↗

PubMed2006

Exploring molecular networks directly in the cell.

The hierarchy of cell function comprises at least four distinct functional levels: genome, transcriptome, proteome, and toponome. The toponome is the entirety of all protein networks traced out directly as patterns on the single cell level in the natural environment of cells in situ (e.g. tissues). In this work a photonic microscopic robot technology (MELK) capable of tagging and imaging hundreds (and possibly thousands) of different molecular components (e.g. proteins) of morphologically-intact fixed cells and tissue have been developed. MELK data sets represent multidimensional vectors of the topologically determined arrangements of proteins within the cell. The data, assembled in a toponome dictionary of the cell, give rise to a new concept for target and drug lead discovery.

Eukaryotic Cells↗