Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “single cell integration”

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 73 records · Page 4Linked to original sources
PubMed2026

Unveiling the power of TIIC: A prognostic tool for esophageal adenocarcinoma.

BACKGROUND: Esophageal adenocarcinoma (EAC) remains a lethal malignancy with limited prognostic tools for guiding immunotherapy. Tumor-infiltrating immune cells (TIICs) play a critical role in EAC prognosis and treatment response. METHODS: We integrated single-cell RNA sequencing and bulk transcriptome data from TCGA and GEO databases. TIIC-specific RNAs were identified via tissue specificity index calculation combined with machine learning feature selection. Twenty machine learning algorithms were benchmarked to construct an optimal TIIC signature score (TIIC-Score) based on the comprehensive C-index. Immunotherapy response, genomic mutation, and copy number variation were analyzed. Summary-data-based Mendelian randomization (SMR) and two-sample Mendelian randomization (MR) were performed to explore genetic associations. Core prognostic TIIC-related genes were functionally validated in esophageal cancer cell lines through loss-of-function assays. RESULTS: The TIIC-Score demonstrated robust prognostic value for 1-, 2-, and 3-year overall survival across multiple cohorts, outperforming 22 published models. High TIIC-Score was associated with poor survival and increased chromosomal instability. Mutation profiling revealed high frequencies of TP53 (78.2%), TTN (48.7%), and SYNE1 (30.8%). MR analysis identified a significant association between gastro-oesophageal reflux and EAC risk at SNP rs8130507. Functionally, CCNI was upregulated in esophageal cancer cells, and its knockdown suppressed malignant phenotypes while promoting apoptosis, supporting its pro-tumorigenic role. CONCLUSION: The TIIC-Score provides a novel prognostic framework for EAC that effectively stratifies patient risk and may help identify individuals most likely to benefit from immunotherapy.

Esophageal adenocarcinoma↗
PubMed2026

Proteomics-based approaches to neutrophil biology.

INTRODUCTION: Neutrophils are central effectors of innate immunity and key contributors to inflammation, host defense, and tissue injury across a wide range of physiological and pathological contexts. Due to their short lifespan, rapid activation, and extensive post-translational regulation, comprehensive molecular characterization of neutrophil function requires approaches that go beyond transcriptomics or marker-based analyses. AREAS COVERED: This review summarizes how proteomic technologies have advanced the understanding of neutrophil biology by enabling unbiased, system-wide profiling of protein abundance, subcellular organization, post-translational modifications, and functional heterogeneity. We discuss global and subcellular proteomics, PTM-centric analyses, and emerging low-input and single-cell proteomic strategies, highlighting recent studies of infection, cancer, metabolic disorders, aging, autoimmune disease, and inflammation. The literature covered includes current large-scale quantitative proteomics, targeted PTMs, and integrative multi-omics studies in both human samples and relevant experimental models. EXPERT OPINION: Proteomics has established neutrophils as highly plastic and context-dependent cells whose functions are governed by coordinated remodeling of signaling, metabolism, and effector pathways. Future progress will depend on expanding neutrophil-specific PTM maps, improving low-input workflows, and integrating single-cell and spatial proteomics. Together, these advances are expected to redefine neutrophil functional states and accelerate translation toward clinically meaningful biomarkers and therapeutic strategies.

Humans↗
PubMed2026

Single-cell and spatial transcriptomics define a progenitor subpopulation and fibroinflammatory niche at the leading edge of parathyroid carcinoma.

Parathyroid carcinoma (PC) is a rare but clinically aggressive endocrine malignancy with limited treatment options and a poorly defined tumor microenvironment (TME). To elucidate its cellular heterogeneity and spatial architecture, we integrated single-cell and spatial transcriptomic profiling with whole-exome sequencing and multiplex immunohistochemistry on eight parathyroid neoplasm specimens, including PC, parathyroid adenoma, and atypical parathyroid tumor. We identified a distinct progenitor-like endocrine subpopulation (Ca-1) enriched in CDC73-mutant PC, exhibiting stem-like properties, elevated cell cycle activity, and pronounced genomic instability. Spatial mapping revealed that Ca-1 cells preferentially localize at the leading edge, forming a fibroinflammatory niche characterized by the enrichment of inflammatory cancer-associated fibroblasts (iCAFs) and SPP1+ macrophages. Within this niche, the dipeptidyl peptidase 4 (DPP4) is selectively expressed in Ca-1 cells and iCAFs, implicating a potential paracrine axis driving stromal remodeling and immunosuppression. These findings suggest that a spatially organized ecosystem may promote PC progression through TME remodeling and highlight the DPP4-CXCL2 axis as a candidate pathway for future investigation in aggressive parathyroid neoplasms.

Humans↗
PubMed2005

Microfluidic selection and retention of a single cardiac myocyte, on-chip dye loading, cell contraction by chemical stimulation, and quantitative fluorescent analysis of intracellular calcium.

A microfluidic method to study the contraction of a single cardiac myocyte (heart muscle cell) has been developed. This method integrates various single-cell operations as well as on-chip dye loading, and quantitative analysis of intracellular calcium concentration, [Ca2+]i. After the channel enlargement by on-chip etching to accommodate large-sized cardiac myocytes, a single cell is selected and retained at a V-shaped cell retention structure within the microchip. Owing to the fragile property of the cardiac myocytes that could easily be damaged by centrifugation, the calcium-sensitive fluorescent dye was loaded in the cell by on-chip dye loading. This on-chip method minimized the damage to the cells from the use of a centrifuge in the conventional method and provided a way of cellular analysis of fragile cells. Subsequently, quantitative analysis of [Ca2+]i of a single cardiac myocyte by fluorescence measurement was achieved for the first time in a microfluidic chip, thanks to the intracellular calcium stimulant of ionomycin. The resting [Ca2+]i of the cardiomyocyte determined was consistent with the literature value. From the spontaneous contraction study, it was found that fluorescence intensity cannot represent the [Ca2+]i variation accurately, which implied the importance of the quantitative analysis of [Ca2+]i.

Animals↗
PubMed2026

Functional screening and single-cell cultivation of marine CO2-fixing bacteria via flow-mode Raman-activated cell sorting.

Most marine CO2-fixing microorganisms remain uncultivated due to strong culture bias and low throughput of conventional approaches, which fail to link in situ function with isolated strains and render slow-growing or low-abundance taxa virtually inaccessible. This study presents an integrated single-cell workflow that incorporates 13C-NaHCO3 labeling, high-throughput flow-mode Raman-activated cell sorting (RACS) and microwell cultivation for the isolation of active CO2-fixing bacteria from the Yellow Sea. Function-guided sorting was achieved by monitoring the 13C-induced Raman shifts of carotenoids (ν1 band: ∼1507 to ∼ 1503.78 cm-1 at 24 h). Genomic and physiological analyses identified Paraburkholderia aromaticivorans FR-4 as a novel facultative chemoautotrophic nitrite-oxidizing bacterium (NOB). Its genome encodes complete nitrite oxidation and Calvin cycle pathways, together with key carbon acquisition genes (carbonic anhydrase, bicarbonate transporter). FR-4 grows autotrophically using NO2- as the electron donor and CO2/HCO3- as the carbon source, confirming its ability to couple nitrite oxidation with carbon fixation, while retaining metabolic flexibility for heterotrophic growth. By directly linking in situ carbon-fixing activity, genotype, and phenotype, this workflow provides a targeted strategy for exploring elusive marine CO2-fixing bacteria and overcomes critical limitations of conventional cultivation.

Carbon-fixing↗
PubMed2026

Mapping the immune-genetic architecture of Epstein-Barr virus-related phenotypes and multiple sclerosis through a single-cell genetic framework for target prioritization and pharmacologic hypothesis generation.

BACKGROUND: Multiple sclerosis (MS) is a severe neuroinflammatory disease causing substantial long-term disability. Strong epidemiologic evidence links Epstein-Barr virus (EBV) exposure with MS risk, but genetic evidence for immune target prioritization in EBV-related phenotypes remains limited. METHODS: We integrated single-cell cis-eQTL data from 14 immune cell types with GWASs of an EBV-related clinical phenotype and MS using a single-cell Mendelian randomization framework with colocalization analyses. Candidate eGenes were evaluated in independent cohorts. For multi-SNP instruments, we performed heterogeneity, pleiotropy, MR-Egger, weighted median, mode-based, and MR-PRESSO sensitivity analyses. We also conducted phenome-wide association analyses and queried DrugBank to annotate candidate compounds targeting prioritized genes. RESULTS: We prioritized 43 immune-cell-specific candidate eGenes with convergent genetic support, including 6 for the EBV-related phenotype and 37 for MS. SERPINB1 in NK cells was associated with increased risk of the EBV-related phenotype, whereas HLA-G was associated with decreased risk. For MS, APOM and MSH5 showed protective associations, while AHI1 showed cell-type-dependent, bidirectional associations across immune lineages. Colocalization and independent cohort evaluation supported these findings. Among FDR-significant multi-SNP associations, MR-Egger intercept tests did not indicate directional pleiotropy, although a small subset showed heterogeneity or MR-PRESSO signals. Phenome-wide analyses identified no significant adverse phenotypic associations among evaluable genes at the prespecified threshold. DrugBank annotation nominated sodium nitroprusside, fasudil, artenimol, and choline as hypothesis-generating compounds for experimental follow-up. CONCLUSIONS: This study provides a single-cell genetic framework for prioritizing immune-cell-specific candidate targets for EBV-related phenotypes and MS, and nominates genetically supported targets and pharmacologic hypotheses for experimental investigation.

Humans↗
PubMed2026

Integration of Genome-Wide Association Studies With Single-Cell and Bulk Expression Quantitative Trait Locus to Identify Stroke Susceptibility Genes.

BACKGROUND: Previous studies have integrated genome-wide association studies with expression quantitative trait locus (eQTL) data from bulk tissues to identify stroke susceptibility genes. However, eQTL data exhibit high cell-type specificity, and genetic variants may have distinct effects across stroke subtypes. METHODS: We applied the summary-data-based Mendelian randomization (MR) method to integrate eQTL data from 7 brain cell types with genome-wide association studies data for 5 stroke phenotypes (stroke, ischemic stroke, cardioembolic stroke, large artery stroke, and small vessel stroke). Results were compared with summary-data-based MR using eQTL data from 49 tissues in the Genotype-Tissue Expression project. Robustness of significant single-cell summary-data-based MR associations was assessed via MR and colocalization analyses. Further evaluations included single-cell RNA-seq differential expression, protein-protein interaction, druggability, and phenome-wide association studies. RESULTS: Single-cell summary-data-based MR identified many novel significant genes not detected using bulk tissue eQTL data. Validated associations revealed 2 stroke risk genes (LRCH1, ICA1L), 3 stroke protective genes (AHI1, LYRM9, CENPQ), 2 large artery stroke risk genes (LIPA, ELL), and 1 ischemic stroke protective gene (CENPQ). Single-cell RNA-seq showed significantly increased LIPA expression in mouse stroke samples compared with controls. Protein-protein interaction and druggability analyses, along with phenome-wide association studies, prioritized LIPA and LRCH1 as potential therapeutic targets for stroke while indicating possible adverse effects. CONCLUSIONS: Integrating single-cell eQTL with stroke-subtype genome-wide association studies uncovers novel cell-type-specific causal genes and highlights promising therapeutic targets, advancing understanding of stroke pathogenesis.

Genome-Wide Association Study↗
PubMed2026

Cell-of-origin Discovery in Infant Leukemia through Integration of 3D Models and Patient Transcriptomic Data.

Pediatric hematological malignancies remain challenging to investigate and model due to the age group-specificity of certain genetic abnormalities. In utero origin has been demonstrated for a subset of pediatric leukemias, placing their respective cell of origin (CoO) during embryonic development. We recently reported a 3D hemogenic gastruloid (haemGx) model of embryonic blood formation derived from mouse embryonic stem cells, resolving the spatio-temporal complexity of developmental hematopoiesis. Importantly, it allows genetic engineering to introduce disease-relevant mutations. Using haemGx, we modeled the most common acute myeloid leukemia exclusive to infants (infAML), subtype t(7;12)(q36;p13), which arises in utero and is characterized by MNX1 overexpression. Here, we detail a method to define susceptibility to specific mutations that integrate phenotypic and transcriptional changes in the haemGx system and compares them with patient data. By proxy of our MNX1-overexpression haemGx, we show a pipeline from cell engineering to downstream analyses of leukemogenic potential. In particular, we focus on the clinical relevance of the model by integrating single-cell and/or bulk RNA sequencing from the haemGx platform with patient data to extract cellular composition and temporal placement of the putative CoO. This method is adaptable to the introduction of other oncogenic mutations, chromosomal rearrangements, or epigenetic modifications, as well as to chemical perturbations, including drug vulnerability and growth factor dependence. This flexibility allows for broad application across diverse disease contexts, enabling mechanistic dissection of how specific alterations disrupt early developmental trajectories with clinical relevance.

Humans↗
PubMed2026

Hox/Meis-dependent gene-regulatory transition underlies cardiopharyngeal neural crest diversification.

Neural crest cells (NCCs) are multipotent migratory cells essential for cardiac development, yet the lineage trajectories and gene regulatory networks underlying their differentiation in the cardiopharyngeal region remain unclear. Here, we integrate single-cell RNA-seq, spatial transcriptomics, and multiomic analyses to construct a comprehensive map of NCC lineages in developing mouse cardiopharyngeal tissues. We identify a transition from Hox-positive pharyngeal NCCs to Hox-negative intracardiac populations associated with the outflow tract cushion, accompanied by a shift in Meis transcription factor binding and gene-regulatory network architecture. By contrast, NCCs forming the aorticopulmonary septum and great vessel smooth muscle retain distinct Hox-codes. A Meis2-Sox9-Scx gene-regulatory network defines a skeletogenic progenitor-like intermediate state that gives rise to coronary artery smooth muscle and semilunar valves. Our findings suggest that the loss of Hox-dependent regional identity enables pharyngeal NCCs to acquire new fates upon entering the cardiac cushion, providing insight into the developmental origins of coronary and valvular calcification.

Journal Article↗
PubMed2026

Identification of ultrasound-associated gene candidates in myeloid cells and construction of a prognostic risk model for acute myeloid leukemia.

BACKGROUND: Incorporating ultrasound (US) treatment sensitivity analysis may improve the treatment of acute myeloid leukemia (AML). METHODS: This study integrated single-cell and bulk datasets for analysis. Differential expression analysis between US-treated and control samples was performed using limma package. The AUCell package was used to calculate US-associated scores in the single-cell dataset. Differentially expressed genes (DEGs) between the specific groups were identified, followed by intersection analysis with previously identified DEGs. Univariate regression, Least Absolute Shrinkage and Selection Operator (LASSO) analysis (using the glmnet package), and stepwise multivariate regression (using the MASS package) were used to refine the candidate genes and to construct a risk model. The model genes were validated using in vitro experiments. Enrichment analysis was conducted using gene set enrichment analysis (GSEA), and immune infiltration was evaluate by single-sample GSEA (ssGSEA) and ESTIMATE algorithms. The correlations between RiskScores and drug sensitivity were analyzed by oncoPredict package. Finally, tumor mutational burden (TMB) and genomic mutations were compared between the risk groups. RESULTS: Nine prognostic signatures (SPINK2, HNRNPAB, SH3BGRL3, CLEC11A, ITGA4, RPL39L, MX1, HEXIM1, and MAP4K4) were identified. Particularly, low expression of SPINK2 attenuated the activity and invasion of AML cells. High-risk group had higher immune cell infiltration. Eight drugs were predicted to be correlated with the RiskScore model. DNMT3A and RUNX1 showed higher mutation frequencies in the high-risk group, whereas KIT and MUC16 showed higher mutation frequencies in the low-risk group. CONCLUSION: The RiskScore model established in this study provides a theoretical basis for clinically screening responsive populations and optimizing treatment strategies.

Humans↗
PubMed2026

Pan-cancer single-cell atlas of immunotherapy response identifies ZNF385A as a regulator of immune evasion in small cell lung cancer.

Although immune checkpoint inhibitors (ICIs) have revolutionized the treatment landscape of solid tumors, response rates in patients with small cell lung cancer (SCLC) remain limited, and acquired resistance is highly prevalent. The underlying mechanisms of this immunotherapy resistance remain to be fully elucidated. Clinically, SCLC typically manifests as an "immune-cold" tumor, characterized by a low abundance of CD8+ T cell infiltration and the rare formation of tertiary lymphoid structures (TLS). While DNA damage repair (DDR) is closely linked to innate immune responses, how DDR networks orchestrate the SCLC immune microenvironment remains obscure. In this study, we integrated single-cell transcriptomic data (comprising 344,447 high-quality cells) from six cancer types (BCC, CRC, HCC, HNSCC, iCCA, and SCLC). Our comparative analysis revealed a fundamental depletion of TLS-associated cellular subpopulations (e.g., CXCL13+ CD8+ T cells, HLA-DRB5+ B cells, and CXCL9+ dendritic cells) in SCLC, which was significantly correlated with aberrant DDR activity. Through high-dimensional weighted gene co-expression network analysis (hdWGCNA), we identified ZNF385A as the core hub gene within the DDR-associated module. ZNF385A is highly expressed in SCLC and is associated with poorer prognosis. In vitro, ZNF385A depletion suppressed SCLC cell proliferation and induced apoptosis, accompanied by R-loop accumulation and activation of cGAS-STING signaling, indicating a potential link between ZNF385A, genomic stability and tumor-intrinsic innate immune signaling. Collectively, these findings identify ZNF385A as a potential regulator associated with TLS deficiency and immune evasion in SCLC.

Immunotherapy resistance↗
PubMed2026

Genome-wide association analyses highlight the neuronal contribution to multiple sclerosis susceptibility.

Multiple sclerosis (MS) is a chronic inflammatory and neurodegenerative disease. Previous genetic studies have identified susceptibility loci that primarily impact immune cells and microglia. Here we performed a multi-ancestry genome-wide association study of 20,831 MS cases and 729,220 controls and identified 236 susceptibility variants outside of the major histocompatibility complex, including four novel genomic loci. We also derived a polygenic score for MS; while optimized for European ancestry, it is informative for African American and Latino individuals. Integrating single-cell data from blood and brain tissue, we identified 76 candidate causal genes. Inhibitory neurons emerged as a key target cell type for MS-associated variants, with seven loci, including STAT3, displaying altered expression only in these cells. The STAT3 variant is also associated with cognition and white matter integrity in individuals with no MS and greater sNfL levels in individuals with MS, suggesting that MS susceptibility may reflect reduced central nervous system resilience to inflammatory challenges.

Humans↗
PubMed2026

A distinct effector B cell population drives autoantibody production in SARS-CoV-2 infection.

Autoantibodies (autoAbs) are linked to mortality and Long COVID, yet their cellular origins remain unclear. We analyzed the INCOV cohort and identified 12 age- and sex-matched participants with varying autoAb abundance and integrated single-cell RNA-seq and ATAC-seq data from B cells, plasma proteomics, proteome-wide autoAb profiling, clinical data, and in vitro assays. AutoAb abundance inversely correlated with neutralizing IgG and declined as infection resolved, paralleling the contraction of atypical memory B cells (AtMs). In vitro, AtMs preferentially differentiated into autoAb-producing antibody-secreting cells upon TLR7/8 stimulation. CD11c+ AtMs (double-negative 2, DN2s) in autoAb-high individuals exhibited increased TLR7 signaling, oxidative stress, and isotype switching, regulated by transcription factors T-bet and XBP1. Integrated genetic and genomic analyses showed that DN2s had the strongest enrichment for autoimmune trait heritability and inferred regulatory effects of autoimmune risk variants among B cell subsets. These findings identify DN2s as key precursors of autoAb-producing cells during SARS-CoV-2 infection.

B cell↗
PubMed2006

Asynchronous progression through the lytic cascade and variations in intracellular viral loads revealed by high-throughput single-cell analysis of Kaposi's sarcoma-associated herpesvirus infection.

Kaposi's sarcoma-associated herpesvirus (KSHV or human herpesvirus-8) is frequently tumorigenic in immunocompromised patients. The average intracellular viral copy number within infected cells, however, varies markedly by tumor type. Since the KSHV-encoded latency-associated nuclear antigen (LANA) tethers viral episomes to host heterochromatin and displays a punctate pattern by fluorescence microscopy, we investigated whether accurate quantification of individual LANA dots is predictive of intracellular viral genome load. Using a novel technology that integrates single-cell imaging with flow cytometry, we found that both the number and the summed immunofluorescence of individual LANA dots are directly proportional to the amount of intracellular viral DNA. Moreover, combining viral (immediate early lytic replication and transcription activator [RTA] and late lytic K8.1) and cellular (syndecan-1) staining with image-based flow cytometry, we were also able to rapidly and simultaneously distinguish among cells supporting latent, immediate early lytic, early lytic, late lytic, and a potential fourth "delayed late" category of lytic replication. Applying image-based flow cytometry to KSHV culture models, we found that de novo infection results in highly varied levels of intracellular viral load and that lytic induction of latently infected cells likewise leads to a heterogeneous population at various stages of reactivation. These findings additionally underscore the potential advantages of studying KSHV biology with high-throughput analysis of individual cells.

Antigens, Viral↗
PubMed2026

Nrsn1-Smarcc1 Coupling Regulates Neural Stem Cell Differentiation and Chronic-phase Recovery After Ischemic Stroke.

Stroke remains a leading cause of long-term neurological disability worldwide, largely due to irreversible neuronal loss and the limited regenerative capacity of the adult mammalian brain. Neural stem cells (NSCs) in the adult brain possess the potential to generate new neurons after injury, yet the molecular mechanisms regulating their neuronal differentiation following ischemic insult remain incompletely understood. Here, integrating single-cell multi-omics analyses with spatial transcriptomics, we systematically delineated cell type-specific spatiotemporal dynamics in the striatum of a mouse model of ischemia-reperfusion injury. We identified Neurensin 1 (Nrsn1) as a gene markedly upregulated during NSC-derived neuronal differentiation in the recovery phase. Mechanistically, Foxa2 directly activates Nrsn1 transcription, whereas Nrsn1 promotes neuronal differentiation by facilitating the nuclear translocation of the chromatin-remodeling factor Smarcc1 in vitro. In vivo, both endogenous NSCs and transplanted NSCs overexpressing Nrsn1 significantly enhanced neuronal regeneration and improved functional recovery in mice subjected to middle cerebral artery occlusion and reperfusion (MCAO/R). Collectively, these findings identify Nrsn1 as a key regulator of NSC neuronal differentiation and uncover a Nrsn1-Smarcc1 coupling mechanism that promotes neural regeneration after ischemic brain injury, highlighting a potential molecular target for strategies aimed at enhancing post-stroke recovery.

Foxa2↗
PubMed2026

scGPA: an LLM-assisted workflow for directional virtual gene perturbation analysis from single-cell transcriptomes.

BACKGROUND: Existing virtual perturbation methods can often infer directional changes by comparing predicted post-perturbation expression profiles with control cells. However, workflows that directly return direction-specific downstream candidate genes together with confidence scores, evidence support and interpretable summaries remain limited. We developed scGPA, an LLM-assisted workflow system for directional single-cell virtual gene perturbation analysis. METHODS: scGPA starts from raw single-cell RNA sequencing data and performs quality control, normalization, dimensionality reduction, clustering and cell-group selection. It then constructs cell-group-specific wild-type regulatory networks using repeated subsampling, principal component regression (PCR)/Ridge-based network inference and CP tensor denoising. Based on these networks, scGPA simulates dose-aware virtual knockdown of the target gene and applies signed perturbation propagation to estimate the magnitude and direction of downstream transcriptional responses. LLM assistance is used for marker-based cell-type annotation, evidence-guided candidate prioritization and user-facing biological summarization. RESULTS: We benchmarked scGPA across five public Perturb-seq datasets and compared its performance with GEARS, scGPT and a random baseline. The overall correct prediction rate of scGPA was 23.0%, exceeding those of GEARS (20.7%), scGPT (15.1%) and the random baseline (13.6%). These results indicate that scGPA achieved a higher correct prediction rate than the two comparator models and the random baseline. We subsequently evaluated scGPA using a public osteosarcoma single-cell dataset and performed qRT-PCR validation in 143B osteosarcoma cells. Among genes with significant experimental changes, scGPA achieved a directional concordance of 76.9%. When all tested downstream genes were counted, 37.0% were directionally correct, 51.9% showed no significant change and 11.1% changed in the opposite direction. CONCLUSIONS: scGPA provides a practical workflow system for predicting and prioritizing direction-specific downstream transcriptional responses after target-gene perturbation. By integrating single-cell regulatory network inference, signed virtual perturbation and LLM-assisted interpretation, scGPA supports target-gene function inference and downstream mechanistic investigation from single-cell transcriptomic data.

Single-Cell Gene Expression Analysis↗
PubMed2026

Single-cell transcriptional profiling identifies the swimming crab Portunus trituberculatus in response to bacterial infection.

Crustaceans rely entirely on innate immunity, yet the cellular composition, functional specialization, and pathogen-induced remodeling of their immune system remain poorly resolved. Here, we generated a high-resolution single-cell transcriptomic atlas of hemocytes from the swimming crab Portunus trituberculatus following Vibrio parahaemolyticus infection using 10× Genomics scRNA-seq. Seven putatively distinct hemocyte clusters were identified, including granulocytes, semigranular hemocytes, prohemocytes, unresolved hemocytes, hyalinocyte-like hemocytes, biosynthetically active secretory hemocytes, and regulatory hemocytes. Although the overall cellular composition remained relatively stable after infection, hemocytes exhibited pronounced cluster-specific transcriptional reprogramming involving Toll/NF-κB signaling, antimicrobial peptide synthesis and metabolic rewiring. By integrating single-cell and bulk transcriptomes, we identified multiple anti-lipopolysaccharide factors (ALFs) as key secretory effectors and experimentally validated their antibacterial activities. FITC-based bacterial engulfment assays and RNA-seq of sorted phagocytes demonstrated that phagocytic capability was shared across multiple hemocyte clusters. Notably, the immunoglobulin superfamily receptor DSCAM displayed extensive alternative splicing and strong infection-induced activation in unresolved hemocytes. Immune-training experiments showed that prior bacterial exposure was associated with altered DSCAM expression and reduced early cumulative mortality upon secondary challenge, suggesting a memory-like immune phenotype. These findings provide a foundational framework for understanding crustacean immunity and advancing disease-resistant breeding in aquaculture.

Antimicrobial peptides↗
PubMed2026

Moderated designs can balance between batch-effect mitigation and cell loss due to hashtag-assisted pooling in single-cell experiments.

Minimizing experimental noise is integral to robust data generation in single-cell omics. The current standard for avoiding batch effects during sample processing is barcode- or hashtag-assisted combining of different experimental treatments into one pool, allowing all samples to be subject to the technical protocols uniformly. The final data points for each treatment group are then computationally separated based on the original hashtag labels. Clearly, whereas hashtagging all groups and pooling them in a single well is expected to minimize batch effects, the procedure can also lead to a loss of cells that cannot be confidently decoded during the computational demultiplexing step. Here, we examine four alternate experimental designs, namely, compound, reference, chain, and confounded, that could be used instead of a single-pool approach and quantify the batch effects as well as cell loss in each case. We find a linear relationship-the percentage of cells lost is double the number of hashtags used in the experiment. We use these analyses to identify experimental designs that can successfully mitigate batch effects while minimizing multiplexing, hence the cell loss, in each well. Although a reference design offers the best overall performance, this study can help individual investigators choose particular approaches that are best suited for their biological questions.

Journal Article↗