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 343 records · Page 19Linked to original sources
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↗
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↗
PubMed2026

Deciphering novel targets in salivary gland pleomorphic adenoma by integrating plasma proteomics and parotid transcriptomics analyses.

BACKGROUND/PURPOSE: Pleomorphic adenoma (PA) is the most common salivary gland benign tumor, with its molecular drivers elusive due to a lack of experimental models. This study aimed to decipher novel targets in PA by systematically integrating plasma protein quantitative trait loci (pQTL)-based Mendelian randomization (MR) with multi-omics profiling of parotid gland tissues. MATERIALS AND METHODS: We performed two-sample MR using 5450 plasma pQTLs and genome-wide association study summary for benign or broader salivary gland diseases from FinnGen consortium. Bulk RNA-sequencing (RNA-seq) and single-cell RNA-seq (scRNA-seq) comparing PA to normal tissue were used for transcriptomic validation. Immunohistochemistry (IHC) was applied for protein-level validation in human PA, adenoid cystic carcinoma (ACC), and murine inflammatory lesions. RESULTS: MR identified 12 plasma proteins associated with benign salivary gland tumor risk. Transmembrane serine protease 6 (TMPRSS6) was the only protein significantly risk-increasing for both benign and broader salivary gland diseases. Strikingly, mitogen-activated protein kinase kinase 4 (MAP2K4) showed opposite MR effects between benign and all-lesion outcomes. Bulk RNA-seq showed limited concordance with MR findings, while scRNA-seq revealed a unique plastic epithelium and partially validated candidates at cellular resolution. Critically, IHC confirmed MAP2K4 protein overexpression specifically in human PA, but not in ACC or inflammatory lesions, while TMPRSS6 was downregulated in established pathologies despite its genetic risk association. CONCLUSION: By integrating plasma proteome-based causal inference with parotid tissue multi-omics, this study unveils MAP2K4 as a potential PA-specific driver. This integrative framework provides novel, context-specific targets for further functional investigation in salivary gland tumorigenesis.

Gene expression profiling↗
PubMed2025

Using Callus as an Ex Vivo System for Chromatin Analysis.

Next-generation sequencing has revolutionized epigenetics research, enabling a comprehensive analysis of DNA methylation and histone modification profiles to explore complex biological systems at unprecedented depth. Deciphering the intricate epigenetic mechanisms that regulate gene activity presents significant challenges, including the issue of analyzing heterogeneous cell populations in bulk. Bulk analysis introduces bias and can obscure crucial information by averaging readouts from distinct cells. Various approaches have been developed to address this issue, such as cell-type-specific enrichment or single-cell sequencing techniques. However, the need for transgenic lines with fluorescent markers, along with technical challenges such as efficient protoplast isolation and low yield, limits their widespread adoption and use in multi-omic studies. This review discusses the pros and cons of these approaches, providing a valuable basis for selecting the most suitable strategy to minimize heterogeneity. We will also highlight the use of cotyledon-derived callus as an ex vivo system as a simple, accessible, and robust platform for enabling high-throughput multi-omic analyses.

Chromatin↗
PubMed2025

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

SUMMARY: Ultrafast mapping of short reads via lightweight mapping techniques such as pseudoalignment has significantly accelerated transcriptomic and metagenomic analyses 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 2.8 times faster than Chromap (the second fastest approach) while using only 33% of the memory required by Chromap. 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 color-enhanced 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. AVAILABILITY AND IMPLEMENTATION: Alevin-fry-atac is written in Rust and C++17, and is freely-available under a BSD 3-clause license. It is integrated into piscem (https://github.com/COMBINE-lab/piscem) and alevin-fry (https://github.com/COMBINE-lab/alevin-fry), and is also supported directly as part of simpleaf (https://github.com/COMBINE-lab/simpleaf).

Single-Cell Analysis↗
PubMed2026

Pan-cancer characterization of HMGA1 reveals its oncogenic role in tumor microenvironment and stemness: functional validation in pancreatic cancer migration and invasion.

BACKGROUND: HMGA1 is a chromatin-associated oncogenic factor implicated in tumor progression, epithelial-mesenchymal transition (EMT), stemness, and metastasis. However, its pan-cancer expression and prognostic patterns, epigenetic activation, and relationship with malignant-cell stemness/plasticity and tumor microenvironment (TME) remodeling in pancreatic adenocarcinoma (PAAD) remain incompletely defined. This study aimed to systematically characterize HMGA1 across cancers and clarify its clinical and biological relevance in PAAD. METHODS: Pan-cancer transcriptomic, clinical, genetic, methylation, immune, and stemness data were integrated from multiple public databases. PAAD single-cell RNA sequencing data were analyzed to localize HMGA1 expression, infer malignant-cell pseudotime, calculate a stemness module score, and assess ligand-receptor communication using CellChat. Public HMGA1-knockdown RNA sequencing data were reanalyzed to evaluate transcriptional remodeling. The Cancer Genome Atlas (TCGA)-PAAD expression and methylation data were used to assess TME-remodeling, immune-suppression, stemness/plasticity, cytokine/chemokine, checkpoint, and promoter-methylation features. HMGA1 expression and function were further examined using immunohistochemistry (IHC), quantitative real-time polymerase chain reaction, Western blotting, wound-healing assays, and Transwell migration and invasion assays. RESULTS: HMGA1 was upregulated in most tumor types, and high expression was associated with unfavorable survival in multiple cancers, including PAAD. In PAAD, HMGA1 was enriched in malignant epithelial cells and positively correlated with pseudotime (Spearman's rho =0.594), while the stemness module score increased along pseudotime (rho =0.748). HMGA1-high malignant cells showed markedly stronger CellChat-inferred outgoing communication, predominantly involving extracellular matrix (ECM)-receptor, adhesion-related, and selected immunomodulatory ligand-receptor axes. HMGA1 knockdown was associated with broad remodeling of EMT, TGF-β, Hedgehog, IL6/JAK/STAT3, KRAS, and cancer stem cell/stemness-related programs rather than uniform suppression of these programs. HMGA1 promoter methylation was inversely correlated with HMGA1 expression (rho =-0.633) and the TME-remodeling score (rho =-0.347). HMGA1 was associated with selected mediators, including PPIA, PLAU, ANXA1, LGALS9, TGFB1, CD276, and CD47, but not with a generalized checkpoint-high phenotype. Functionally, HMGA1 knockdown significantly reduced pancreatic cancer (PC) cell migration and invasion. CONCLUSIONS: These findings support an association-based model in which promoter hypomethylation-associated HMGA1 activation is linked to malignant epithelial stemness/plasticity, ECM/adhesion-dominant TME remodeling, selected immunomodulatory programs, and aggressive PAAD phenotypes. Further mechanistic and clinical validation is required before HMGA1 can be used for therapeutic stratification or immunotherapy-response prediction.

HMGA1↗
PubMed2025

A single-nucleus transcriptomic atlas of the adult Aedes aegypti mosquito.

The female Aedes aegypti mosquito's remarkable ability to hunt humans and transmit pathogens relies on her unique biology. Here, we present the Aedes aegypti Mosquito Cell Atlas, a comprehensive single-nucleus RNA sequencing dataset of more than 367,000 nuclei from 19 dissected tissues of adult female and male Aedes aegypti, providing cellular-level resolution of mosquito biology. We identify novel cell types and expand our understanding of sensory neuron organization of chemoreceptors to all sensory tissues. Our analysis uncovers male-specific cells and sexually dimorphic gene expression in the antenna and brain. In female mosquitoes, we find that glial cells in the brain, rather than neurons, undergo the most extensive transcriptional changes following blood feeding. Our findings provide insights into the cellular basis of mosquito behavior and sexual dimorphism. The Aedes aegypti Mosquito Cell Atlas resource enables systematic investigation of cell type-specific expression across all mosquito tissues.

Aedes aegypti↗
PubMed2026

Necroptosis in alveolar epithelium orchestrates lung ischemia-reperfusion injury: a multi-omics study.

BACKGROUND: Lung ischemia-reperfusion injury (LIRI) is a leading cause of early morbidity and mortality following lung transplantation and other cardiopulmonary procedures. It is characterized by acute sterile inflammation driven by regulated cell death (RCD). While various RCD modalities, including apoptosis, necroptosis, pyroptosis, and ferroptosis, have been implicated in lung injury, their relative contributions and distinct activation patterns in LIRI remain poorly defined. METHODS: We employed an integrated multi-omics approach combining transcriptomics and proteomics with histological and functional validations in a murine hilar clamping model of LIRI. Key findings were further corroborated using single-cell RNA sequencing (scRNA-seq) data from human lung transplant recipients. The functional role of necroptosis was validated using pharmacological inhibitors (Nec-1, GSK'872) and Mlkl-deficient (Mlkl-/-) mice. RESULTS: LIRI triggered acute, time-dependent lung injury peaking within 24 h of reperfusion. Although transcriptomic profiling suggested broad activation of multiple RCD pathways, proteomic and biochemical analyses revealed a distinct landscape in our experimental setting: markers of apoptosis, pyroptosis, and ferroptosis were either downregulated or showed no significant positive correlation with injury severity and inflammatory peaks. In contrast, the necroptotic pathway emerged as a highly activated modality. Specifically, necroptosis, marked by phosphorylated RIPK1, RIPK3, and MLKL, was localized primarily in alveolar epithelial cells, correlated strongly with cytokine release and histological lung injury, and preceded the inflammatory response. Pharmacological inhibition or genetic ablation of necroptosis significantly attenuated tissue damage and inflammation. This pronounced necroptotic signature appeared distinct from the broad multi-pathway activation observed in lipopolysaccharide (LPS)-induced lung injury. Translational analysis of human scRNA-seq data further confirmed the selective upregulation of necroptosis signatures in alveolar type 2 (AT2) cells following lung transplantation. CONCLUSION: Our multi-omics analysis identifies necroptosis, particularly in alveolar epithelial cells, as a critical driver of sterile inflammation and tissue injury in the early phase of LIRI. Targeting alveolar epithelial necroptosis may represent a precise and promising therapeutic strategy for lung transplantation and ischemia-reperfusion-associated pulmonary disorders.

Animals↗
PubMed2025

Multi-step genomics on single cells and live cultures in sub-nanoliter capsules.

Single-cell sequencing methods uncover natural and induced variation between cells. Many functional genomic methods, however, require multiple steps that cannot yet be scaled to high throughput, including assays on living cells. Here we develop capsules with amphiphilic gel envelopes (CAGEs), which selectively retain cells and large analytes while being freely accessible to media, enzymes and reagents. Capsules enable high-throughput multi-step assays combining live-cell culture with genome-wide readouts. We establish methods for barcoding CAGE DNA libraries, and apply them to measure persistence of gene expression programs in cells by capturing the transcriptomes of tens of thousands of expanding clones in CAGEs. The compatibility of CAGEs with diverse enzymatic reactions will facilitate the expansion of the current repertoire of single-cell, high-throughput measurements and extension to live-cell assays.

Journal Article↗