SEARCH · Search PubMed

Results for “single-cell transcriptome sequencing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3Linked to original sources

PubMed2026

Generation and validation of a Myh11Dre-Spp1Cre intersectional mouse model for lineage tracing of disease-associated smooth muscle cell states.

BACKGROUND: Phenotypic modulation of vascular smooth muscle cells (VSMCs) is a hallmark of vascular remodeling and cardiovascular disease. Recent lineage-tracing and single-cell transcriptomic studies have identified secreted phosphoprotein 1 (SPP1) as a prominent marker associated with disease-associated VSMC states, particularly those linked to fibrotic remodeling and vascular calcification. However, the cellular origins and fate of SPP1-associated VSMC populations remain incompletely understood. METHODS AND RESULTS: We generated a novel Spp1-rSTOPr-Cre (Spp1Cre) knock-in mouse line in which Cre recombinase is expressed from the endogenous Spp1 locus following Dre-mediated excision of a rox-flanked transcriptional STOP cassette. Correct targeting of the knock-in allele was validated by internal, 5' junction, 3' junction, and long-range PCR analyses, as well as Sanger sequencing. To establish an intersectional lineage-tracing strategy, Spp1Cre mice were crossed with Myh11DreERT2 and Rosa26-RSR-LSL-tdTomato-LSL-eGFP reporter mice, enabling permanent labeling of VSMC-derived populations following activation of the endogenous Spp1 locus. Under physiological conditions, eGFP-positive cells were detected at low frequency within the vascular wall and were predominantly negative for the contractile markers ACTA2 and MYH11. As a proof-of-principle application, eGFP-positive cells markedly expanded within atherosclerotic lesions induced by AAV-PCSK9D377Y and high-fat diet feeding. These lineage-traced cells remained largely ACTA2- and MYH11-negative, consistent with a modulated phenotype. Notably, only a minority of eGFP-positive cells expressed SPP1 or fibronectin at the time of analysis, demonstrating the utility of permanent lineage tracing for tracking cells with a history of endogenous Spp1 activation during vascular remodeling. CONCLUSION: We report the generation and validation of a novel Myh11Dre-Spp1Cre intersectional mouse model for lineage tracing of VSMC-derived populations that have activated the endogenous Spp1 locus. This genetic resource provides a valuable platform for investigating the origin, fate, and phenotypic evolution of Spp1-associated VSMC populations during vascular remodeling and cardiovascular disease.

Animals↗

PubMed2026

ScRNA-seq analysis reveals the effects of nitrite stress on the endocrine system of the eyestalk in Litopenaeus vannamei.

Nitrite is a harmful substance generated in Litopenaeus vannamei farming systems, largely originating from the inadequate breakdown of surplus feed and shrimp feces. Its accumulation in the water can affect the growth and physiological functions of shrimp, damage the immune system, and even cause mass mortality, thus becoming a key environmental factor restricting the green development of the industry. Under nitrite stress, the eyestalk, as an important neuroendocrine regulatory center in crustaceans, participates in the stress adaptation of the organism and exerts a protective effect by regulating energy metabolism and immune function. However, the molecular regulatory mechanism of the eyestalk in response to nitrite stress remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) technology was used to analyze the heterogeneity of eyestalk cells in L. vannamei under nitrite stress. A total of 18, 394 high-quality cells were obtained, and six major cell subpopulations, including Neurosecretory cell, Motor neuron, Sensory neuron, Interneuron, Neurogliocyte, and Support cell, were identified. Differential expression analysis identified 839 differentially expressed genes, and different cell types showed distinct specific responses to nitrite stress. Functional enrichment analysis indicated that pathways such as glycolysis, oxidative phosphorylation, ribosome function, and endoplasmic reticulum protein processing were significantly activated, while signal transduction and DNA repair-related pathways were inhibited. Further analysis revealed that nitrite stress could induce mitochondrial function changes and trigger oxidative stress, thereby affecting the neuroendocrine system function of the eyestalk. This study provided insights into transcriptomic responses of the eyestalk to nitrite stress at the single-cell level, laying a theoretical foundation for the management of aquaculture environments.

Animals↗

PubMed2026

Single-cell sequencing reveals synovial fluid γδ T-cell expansion in equine experimental osteoarthritis.

OBJECTIVE: Define temporal cellular changes following joint injury using single-cell RNA sequencing in experimental equine posttraumatic osteoarthritis (PTOA). METHODS: PTOA was induced in 4 Quarter Horses (3 to 5 years) via carpal osteochondral fragmentation and high-speed treadmill exercise. Synovial fluid (SF) cells and synovium were sampled over 18 weeks (November 2023 to April 2024). Single-cell suspensions were processed (10x Genomics Chromium iX), then aligned to the equine genome (Cell Ranger). Downstream analysis was completed in the R Seurat package. Differential gene expression (log2[fold change] > 1; P < .05) and differential abundance analyses were performed (P < .1). RESULTS: Cartilage injury had a modest impact on gene expression changes and cell abundance shifts in SF. Integrated analysis of 90,323 SF cells across 4 time points revealed 9 distinct cell types, primarily T cells (73 ± 19%) followed by myeloid cells (20 ± 13%). Subcluster analysis of T cells revealed 9 transcriptomically distinct subtypes (3 CD8, 2 CD4, 3 γδ, and 1 cycling). Differential abundance analyses of temporal changes identified increased γδ T and decreased CD4+ T-cell subsets in joints over time. Expanded populations of IL-23 receptor-positive γδ T cells exhibited increased T-helper 17 signatures. CONCLUSIONS: IL-23 receptor-positive γδ T-cell expansion, associated with joint inflammation, occurred in PTOA. Limitations include small sample size and individual heterogeneity; further investigation over extended timeframe is necessary to confirm whether later stages of the experimental model reflect natural chronic OA. CLINICAL RELEVANCE: Cellular immunotherapy targeting γδ T cells and IL-23/IL-17 blockade may warrant investigation to mitigate equine OA progression.

equine↗

PubMed2026

The current and future perspective of ChickenGTEx project and its applications in precision breeding.

The Chicken Genotype-Tissue Expression (ChickenGTEx) project was established to systematically characterize the regulatory landscape of the chicken genome and to accelerate the translation of functional genomics into precision breeding. By integrating whole-genome sequencing with multi-tissue transcriptomic profiling, ChickenGTEx provides a comprehensive atlas of gene expression regulation across diverse tissues and physiological systems. Current findings demonstrate that complex production traits are governed by coordinated regulatory networks rather than isolated loci, with substantial contributions from tissue-specific gene expression, structural variation, and genotype-by-sex interactions. Sex-dependent regulatory effects further refine the genetic architecture of metabolic, immune, and reproductive traits, highlighting the importance of incorporating sex as a biological variable in genomic analyses. Application of integrative omics frameworks within elite layer populations has revealed multilayer regulatory mechanisms underlying extended laying performance, feed efficiency, metabolic health, and eggshell quality. By partitioning phenotypic variance into genetic, regulatory, and host-microbiome components, these approaches move beyond association-based mapping toward causal inference and biological interpretation. Importantly, validated regulatory loci identified through ChickenGTEx and related analyses provide actionable markers for genomic selection and rational targets for precision genome modification. Looking forward, continued expansion of regulatory atlases, incorporation of single-cell and longitudinal data in diverse environmental conditions, and integration of functional annotation into breeding pipelines will further enhance prediction accuracy and sustainable genetic improvement. The ChickenGTEx project thus represents a foundational platform bridging functional genomics and practical poultry breeding.

Animals↗

PubMed2026

Identification of immune cell type-specific susceptibility genes in multiple cancers using transcriptome-wide association studies.

BACKGROUND: Transcriptome-wide association studies (TWAS) integrate gene expression and genome-wide association studies (GWAS) to identify disease susceptibility genes. Because gene expression varies substantially across cell types within tissues, cell type-specific prediction models may enhance the power of TWAS. METHODS: We conducted cell type-specific TWAS leveraging single-cell RNA sequencing data from the OneK1K cohort (14 immune cell types, 1.27 million cells) and GWAS summary statistics for 7 cancers (>290 000 cases in total). To improve prediction accuracy, we developed a modeling framework that incorporates shared gene expression effects across cell types. RESULTS: At a false discovery rate of 5%, we identified 106 (Bonferroni 5%: 13) previously unreported loci for breast cancer, 51 (4) loci for prostate cancer, 11 (4) loci for lung cancer, 39 (5) loci for melanoma, 9 (1) loci for ovarian cancer, and 2 (1) loci for diffuse large B-cell lymphoma, with most genes exhibiting cell type specificity. Gene set analyses confirmed joint associations of unreported genes with breast and prostate cancer risk in UK Biobank data. Additional lung tissue single-cell RNA sequencing data with 113 individuals validated 18 of 32 (56.3%) statistically significant genes for lung cancer. Across cancers, 139 statistically significant genes were shared by at least 2 cancer types and were primarily enriched in specific immune cell types. CONCLUSION: Cell type-specific TWAS improve the identification of novel cancer susceptibility loci and provide insights into the immune landscape of cancer etiology.

Humans↗

PubMed2026

From wild to domestic: Single-cell transcriptomic perspectives on hippocampal regulation and evolution.

How domestication shapes brain evolution remains an open question. In this study, we integrated single-nucleus RNA sequencing (snRNA-seq), population genomics, and machine learning to investigate the hippocampal evolution under domestication. Across-species comparisons revealed that hippocampal cell type profiles are largely conserved across vertebrate species, while supporting the presence of adult hippocampal neurogenesis in birds. We further found that domestication and selective breeding likely influence the cellular composition and molecular regulation of the hippocampus. Our findings provide cellular evidence supporting the hypothesis that domestication affects adult hippocampal neurogenesis. Additionally, we showed that genes associated with neural progenitor cells (NPC) states and cell-marker programs are enriched for signatures of selection. Many of these genes function as regulators of neurogenesis and pathways mediating stress and fear reduction. Specifically, we identified selection at the FKBP5 promoter that may influence its expression in the NPC lineage, potentially contributing to stress-response regulation during domestication. Collectively, these results suggest that domestication is associated with hippocampal remodeling as part of an adaptive response to human-managed environments. This study provides a cellular and genetic perspective on how domestication reshapes the brain and offers a basis for further investigation into the mechanisms of neural evolution within the context of microevolution.

Animals↗

PubMed2026

Lactylation-related immune-metabolic dysregulation defines prognostic and therapeutic stratification in lung adenocarcinoma.

BACKGROUND: Lactylation links lactate metabolism with inflammatory signaling and immune regulation in tumors. However, its cellular distribution and translational value in lung adenocarcinoma (LUAD) remain unclear. METHODS: Single-cell RNA-sequencing datasets GSE189357 and GSE171145 were integrated to characterize lactylation-related activity, intercellular communication, and malignant epithelial cell states in LUAD. Single-cell-derived lactylation-related differentially expressed genes were mapped to TCGA-LUAD and multiple GEO cohorts. Univariate Cox regression and machine learning algorithms were used to construct a lactylation-related prognostic signature (LRPS). The associations of LRPS with prognosis, immunotherapy response, drug sensitivity, genomic alterations, immune infiltration, and inflammation- and metabolism-related pathways were evaluated. KRT7 was further validated using virtual knockout analysis, spatial transcriptomics, and in vitro and in vivo experiments. RESULTS: lactylation-related transcriptional activity showed heterogeneous distribution across LUAD cell populations and was associated with altered cell-cell communication. In malignant epithelial cells, LRTS-high and LRTS-low states exhibited distinct metabolic, inflammatory, and tumor-related pathway activities. LRPS showed stable prognostic performance in TCGA-LUAD and multiple GEO cohorts and remained an independent prognostic factor. Low LRPS was associated with greater potential benefit from immunotherapy, whereas different LRPS groups displayed distinct drug sensitivity, genomic alteration, and immune microenvironment patterns. KRT7 was highly expressed in LUAD and associated with poor prognosis. KRT7 knockdown suppressed LUAD cell proliferation, migration, invasion, colony formation, and tumor growth in vivo. CONCLUSIONS: This study identifies lactylation-related immune-metabolic dysregulation as a clinically relevant feature of LUAD and develops a single-cell-guided LRPS for prognosis and therapeutic stratification. KRT7 emerged as an LRPS-related functional candidate with experimentally supported roles in malignant LUAD phenotypes.

Immunotherapy↗

PubMed2026

Immature Neutrophil Programs Associate With Burn Mortality and Extend Across Critical Illnesses.

Severe burns provoke a systemic "genomic storm," yet cell states associated with divergent outcomes remain unclear. We profiled blood cells by single-cell RNA-Sequencing (73 014 cells) from adult patients with burn injuries within postburn day 17 (n = 4) and healthy donors (n = 5), integrated data with bulk signatures of burn size, inhalation injury, and mortality, and evaluated clinical associations in the American Burn Association National Burn Repository. Burn was associated with emergency hematopoiesis marked by expansion of hematopoietic stem/progenitor-like cells, immature neutrophils, and plasmablast/plasma cell states, alongside depletion of naïve CD4+/CD8+ T cells and dendritic cells. Larger burns (>20% TBSA) showed enrichment of humoral transcriptional programs, including plasmablast/plasma cell activation and suppression of cytotoxic CD8+ T-cell states. In multivariable models, inhalation injury was a stronger predictor of death (adjusted odds ratio [OR] 1.9) than burn size (adjusted OR 1.1) and shared greater overlap with the most perturbed single cells in non-survivors; 55% of co-perturbed cells were neutrophils, implicating granulocyte dysregulation as a common lethal axis. We identified a neutrophil-specific 5-gene panel (OLFM4, RETN, LCN2, ARG1, and BTNL3) that discriminated survivors vs non-survivors after burns (area under the curve [AUC] > 0.9) and generalized to trauma (n = 158; AUC 0.81) and intensive care unit COVID-19 (n = 103; AUC 0.75), providing information orthogonal to conventional biomarkers and severity scores. Cytomorphology corroborated transcriptomic immaturity, with ~2-fold higher band neutrophils and larger neutrophil size in a fatal case. Computational drug-reversal analysis highlighted galectin-1 inhibition as a candidate modulator of mortality-associated neutrophil programs. Together, our findings suggest that immature neutrophils represent a shared immune feature across severe burns and other forms of critical illness.

Humans↗

PubMed2026

Atherosclerotic plaque fibroblasts derive from adventitial and medial Pdgfra-lineage-positive cells and predominantly maintain fibroblast identity.

AIMS: Fibroblasts are mesenchymal cells in the healthy vascular adventitia. In atherosclerosis, single-cell sequencing datasets suggest fibroblasts are abundant in plaques. However, their identity, origin, and fate during plaque progression remain unclear, which we aim to unravel here. APPROACH AND RESULTS: To robustly define fibroblast identity, origin, and fate, we employed meta-analyses of 54 single-cell RNA sequencing libraries, including murine smooth muscle cell (Myh11) and endothelial cell (EC) (Cdh5) lineage reporter mice with and without atherosclerosis; human control and atherosclerotic arteries; and murine adventitia and atherosclerotic plaques processed separately from low-density lipoprotein (LDL) receptor knockout (Ldlr-/-) mice. These meta-analyses showed that murine and human plaque fibroblast identity was robustly defined by Pdgfra, Pi16, Cygb, and Serpinf1 mRNA. Ninety-five percent of plaque fibroblasts do not derive from the Myh11 lineage, while no Cdh5-lineage-positive cells were present in the fibroblast cluster. We identified five murine arterial fibroblast subsets in atherosclerotic murine aorta: progenitor fibroblasts, matrix fibroblasts, inflammatory fibroblasts, an EC-like fibroblast subset, detected in both adventitia and plaques, and Col5a3+ fibroblasts, unique to the adventitia. We next studied fibroblast identity, origin, and fate using pseudotime analysis and Pdgfra-CreERT2/tdTomato lineage reporter mice (Pdgfra Lin+). Healthy Pdgfra Lin+ reporter mice showed predominant adventitial tdTomato expression, and infrequent medial and intimal Pdgfra Lin+ cells co-expressing MYH11 and PECAM1, respectively. The Pdgfra Lin+ plaque area increased with diet duration. Pdgfra Lin+ cells largely maintain fibroblast identity in the plaque, while <10% co-express SMC markers (MYH11, SM22α), or contribute to ACTA2+ cap cells. ECs gaining mesenchymal markers are transcriptionally distinct from Cdh5-lineage-negative fibroblasts gaining EC markers. Plaque-resident EC-like fibroblasts displayed a mesenchymal-to-endothelial transition transcriptome, which was induced in human primary fibroblasts in vitro by starvation, and dampened or reversed by IL1B, TGFB1, TGFB3, and oxidized LDL. Cross-species integration showed that all murine plaque fibroblasts were conserved in human atherosclerosis, with one additional subset partially resembling murine subsets, and three human-specific subsets. Importantly, human fibroblast subsets differentially correlated to human plaque traits, with EC-like fibroblasts correlating to plaque instability. CONCLUSION: Our results indicate that 95% of plaque-residing fibroblasts are Myh11 Lin- Plaque fibroblasts have a dual origin, predominantly adventitial Pdgfra Lin+ progenitor fibroblasts, with a minor contribution from medial Pdgfra Lin+  Myh11+ SMCs. Most plaque fibroblasts maintain fibroblast identity. Murine plaque fibroblast subsets were conserved in human atherosclerosis. EC-like fibroblasts are linked to human plaque instability. Intervening in progenitor-to-specific fibroblast transitions could present a new avenue to promote plaque stability in atherosclerosis.

Atherosclerosis↗

PubMed2026

Human Monocytic Models Reveal Genotype-Dependent Inflammatory Programs in VEXAS Syndrome.

OBJECTIVES: VEXAS syndrome is a severe X-linked autoinflammatory disorder caused by somatic mutations in ubiquitin-like modifier activating enzyme 1 (UBA1), with clinical outcomes that vary by UBA1 genotype. We aimed to elucidate genotype-specific inflammatory programs and identify potential therapeutic targets. METHODS: We conducted longitudinal deep phenotyping, including whole-blood RNA sequencing (RNA-seq) and clinical activity assessment. Peripheral blood samples were analyzed by single-cell RNA-seq. Human monocytic cell lines harboring each major UBA1 mutation (p.Met41Val, p.Met41Thr, or p.Met41Leu) were generated and subjected to transcriptomic and functional analyses. RESULTS: Thirteen patients with VEXAS syndrome contributed a total of 79 RNA-seq samples. Among genes upregulated in VEXAS syndrome, RNASE1 showed the strongest correlation with longitudinal disease activity (r = 0.70, FDR < 0.05) and was upregulated in patients' monocytes. In UBA1-mutant monocytic cell lines, genotype-dependent ubiquitination defects were observed in a graded manner (p.Met41Val > p.Met41Thr > p.Met41Leu), even in the absence of exogenous stimuli. These defects were accompanied by unfolded protein response activation, increased pro-inflammatory cytokine production, progressive cell death, and RNASE1 upregulation, all following the same graded pattern, recapitulating patient genotype-phenotype associations. Transcriptomic analyses demonstrated enrichment of pro-inflammatory, interferon, and necroptosis signatures in more severe genotypes. Notably, inhibition of receptor-interacting protein kinase 3 (RIPK3) markedly attenuated all pathological features, including RNASE1 upregulation. CONCLUSIONS: Our UBA1-mutant monocytic cell-line models, representing three distinct genotypes, recapitulate genotype-dependent inflammatory phenotypes that can be modulated by RIPK3 inhibition, providing a translational platform for mechanistic investigation and precision therapy development in VEXAS syndrome.

Journal Article↗

PubMed2026

Integrated Genomic and Proteomic Analysis Reveals T-B Lymphocyte Signatures in the MYCN Driven "Immune Desert" of Specific Neuroblastoma Subtypes.

AIMS: This study aims to systematically dissect how MYCN amplification shapes the immunosuppressive tumor microenvironment (TME) in high-risk neuroblastoma, elucidating key mechanisms underlying immune evasion. METHODS: We performed an integrated multi-omics analysis of bulk RNA-seq (n = 721), single-cell RNA-seq (n = 9), proteomic data (n = 49) and spatial transcriptomics (Visium, with external validation in melanoma). Analyses included unsupervised clustering, cell-cell communication inference, transcriptional regulatory network reconstruction, and spatial proximity assessment to map the immune landscape. RESULTS: A distinct molecular subtype (Class C), defined by MYCN amplification and poor prognosis, exhibited a comprehensive "immune desert" phenotype characterized by low immune scores and minimal leukocyte infiltration. Single-cell analysis confirmed significant depletion of T and B lymphocytes within the Class C TME. Dysregulated transcriptional networks were identified, including upregulation of REL and EOMES in T cells-with EOMES potentially driving exhaustion via regulation of Transient Receptor Potential (TRP) genes, and REL inhibition enhancing cytotoxic function in vitro. A unique immunosuppressive B-cell subset (B7) engaged in enhanced crosstalk with exhausted T cells and harbored a MYC-centered network linked to cell cycle dysregulation and poor survival. Spatial transcriptomics revealed significant proximity between B7-active regions and Treg/exhaustion-enriched areas, externally validated in melanoma. Proteomic data validated elevated REL expression in MYCN-amplified tumors. CONCLUSION: This work delineates the immunosuppressive architecture of MYCN-driven neuroblastoma, revealing novel regulatory nodes within specific lymphocyte compartments. Integrating single-cell, spatial, and proteomic evidence, we propose REL inhibition as a therapeutic candidate, the EOMES/TRP axis as a bioinformatically supported hypothesis, and the B7/MYC hub as a hypothesis supported by transcriptomic and spatial evidence.

Humans↗

PubMed2026

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

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

Animals↗

PubMed2026

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

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

Humans↗

PubMed2026

A single-nucleus transcriptome atlas of soybean anthers.

Anther development is crucial for plant sexual reproduction. However, a high-resolution, cell-type-specific transcriptomic atlas of this process is lacking for the legume crop soybean (Glycine max). Here, we construct a comprehensive transcriptional atlas of developing soybean anthers using single-nucleus RNA sequencing (snRNA-seq). We identify and characterize nine distinct cell types spanning both somatic and reproductive lineages. Our analysis reveals robust transcriptional continuity across anther developmental stages and dynamic reprogramming during key transitions. Notably, the shift from diploid meiocytes to haploid unicellular microspores is marked by the induction of previously inactive genes, despite an overall reduction in transcript abundance. Subsequently, within bicellular microspores, generative and vegetative cell lineages exhibit sharply divergent transcriptional programs: generative cells specialize in mRNA export and turnover, whereas vegetative cells up-regulate translational machinery. Evolutionary analysis further indicates that generative-cell-specific genes are subject to more relaxed purifying selection compared to those specific to vegetative cells. Functional validation using mutants generated by CRISPR/Cas9-mediated genome editing and EMS mutagenesis reveals the essential roles of OSD1A and PKSA in pollen development and fertility. This high-resolution atlas provides fundamental insights into the transcriptional regulation of soybean anther development and serves as a valuable resource for manipulating male fertility to advance hybrid breeding programs. The data are available at https://databases.genedenovo.com/pollen.

Glycine max↗

PubMed2026

Decoding the landscape of cell-type-specific co-expressed transcription factors in soybean.

Soybean (Glycine max) is an essential source of protein and oil with high nutritional value for human and animal consumption. To enhance our understanding of soybean biology, it is essential to have accurate information regarding the expression of each of its protein-coding genes. Here, we present Tabula Glycine max, a soybean single-cell resolution transcriptome atlas. This atlas comprises single-nucleus RNA-sequencing data from ten different G. max organs and morphological structures constituting the entire soybean plant. These nuclei are grouped into 156 different clusters based on their transcriptomic profiles. The breadth of various organs, tissues and cell types represented in Tabula Glycine max reveals that the pattern of co-expressed transcription factor genes is sufficient to define most cell types based on their function and organ of origin. Defining cell-type-specific co-expressed transcription factor genes offers a new perspective to engineer cell-type-specific programmes and enhance the biology of unique soybean cell types. This cellular resolution and breadth make the Tabula Glycine max an exceptional resource for the plant and soybean communities.

Journal Article↗

PubMed2026

Mapping ovarian cellular and molecular landscape across the lifespan of women: a scoping review.

BACKGROUND: With growing interest in ART, fertility preservation, and postmenopausal health of women, reproductive medicine is increasingly focused on characterizing oocytes and ovarian tissue composition, as well as understanding the molecular mechanisms that guide ovarian function throughout its lifecycle. High-throughput omics technologies have enabled the characterization of different molecular layers, leading to substantial advances in our understanding of their complex dynamics. However, not all molecular aspects are studied equally, and studies examining the same modalities often show inconsistencies, underscoring the need for data standardization and highlighting the potential for using transformative artificial intelligence and machine-learning (AI/ML) methods for ovary studies. OBJECTIVE AND RATIONALE: This study aims to evaluate how multi-omic studies have advanced our understanding of the ovarian lifecycle from fetal development to postmenopause. We systematically reviewed published studies that have investigated molecular/omic layers, including the genome, methylome, transcriptome, and proteome throughout ovarian development and aging. Our analysis identified key molecular and cellular patterns, highlighted inconsistencies across studies and addressed gaps in data analysis, interpretation, and reproducibility to guide future research. SEARCH METHODS: We conducted a systematic literature search of Medline (PubMed), Embase (Ovid), and Web of Science Core Collection (Clarivate) using a combination of controlled and free text terms for human ovary, oogenesis, folliculogenesis, ovary development and (epi)genome, transcriptome, proteome, and multi-omic mechanisms to find relevant articles published before August 2025. To focus the scope of the current review, studies of domesticated and farm animals, rodents and other model organisms, non-human primates, as well as those examining various human ovarian pathologies were excluded. OUTCOMES: The search identified 23 546 studies for screening, of which 637 full-text studies were assessed for eligibility. Subsequently, we extracted data from 121 studies. Most studies analyzed the transcriptome of oocytes, granulosa cells, and ovarian tissue from reproductive-age individuals (n = 91), with fewer studies examining samples from individuals of advanced reproductive age (n = 45) and fetal (n = 16) samples. Transcriptome analyses were most common (n = 103, 85%), followed by proteome (n = 19, 16%) and epigenome (n = 14, 12%) studies. We found substantial variation in how studies defined and reported participants' groups as well as in their sequencing technologies and data analysis methods, with a lack of standardized reporting of background clinical information, data analysis methods, and pipeline details. The key findings underscore the prevailing consensus on genes defining major ovarian cell types and their roles throughout the ovarian lifespan, from prenatal development to postmenopausal transformation. This review highlighted the underrepresentation of certain patient groups, particularly prepubertal and peri-/postmenopausal individuals, among researched populations, due to obvious clinical and ethical reasons. WIDER IMPLICATIONS: This scoping review offers a comprehensive overview and benchmark of the current state of high-throughput omics-based research on ovarian cellular composition and molecular dynamics. To address these shortcomings, we propose general recommendations for multi-omics ovary studies and emphasize the necessity for more thorough multi-omic data integration by effectively applying novel AI/ML approaches. They can potentially improve the quality of multi-omics analyses at both single-cell and tissue levels despite limited sample sizes and enable integration of molecular profiling data with clinical and radiology datasets, enabling a more comprehensive understanding of ovarian biology. Such advancements can enhance reproducibility of research findings and guide future research to deepen our understanding of ovarian biology and ultimately support the development of medical technologies for better preserving fertility and alleviating infertility. REGISTRATION NUMBER: A protocol was published a priori on the Open Science Framework (https://osf.io/z38gb/).

Female↗

PubMed2026

Amino acid reprogramming and biofilm-specific tricarboxylate transporters in PET-degrading Piscinibacter sakaiensis.

Plastic-degrading bacteria predominantly colonize polymer surfaces as biofilms, yet it remains unclear whether the biofilm phenotype contributes to metabolism beyond retaining extracellular enzymes. Here, we combine population-level RNA-sequencing across three conditions-biofilm cells on polyethylene terephthalate (PET), planktonic cells incubated with PET, and planktonic cells on maltose-with single-cell Raman spectroscopy to characterize the PET response of Piscinibacter sakaiensis (formerly Ideonella sakaiensis). This integrated approach reveals two metabolically distinct response layers. A carbon-source-driven response shared by all PET-exposed cells is dominated by a broad amino acid reprogramming, led by upregulation of branched-chain amino acid transport genes, enhanced serine biosynthesis, and reduced chemotaxis. A biofilm-specific layer selectively induces tripartite tricarboxylate transporter genes from three distinct genomic loci. This transcriptional feature is accompanied by a single-cell phenotype consistent with a protein-rich and saturated membrane. These results suggest that biofilm formation is not limited to enzyme retention but is associated with selective activation of transport systems, consistent with a putative role in capturing PET-derived intermediates at the polymer interface. This two-layer model separates general metabolic adaptation to PET from biofilm-specific functions and provides a framework for understanding how surface-associated bacterial physiology contributes to plastic degradation.IMPORTANCEPolyethylene terephthalate (PET) degradation in natural and engineered environments is largely mediated by surface-attached microbial communities, yet the physiological role of biofilm state during plastic degradation remains poorly understood. Using the model PET degrader Piscinibacter sakaiensis, we show that biofilm-associated cells are not simply retained near the polymer surface but exhibit a distinct metabolic program characterized by selective induction of tripartite tricarboxylate transporters. In contrast, extensive amino acid reprogramming occurs in both biofilm and planktonic PET-exposed cells, indicating that it is driven by carbon source rather than surface attachment. These findings reveal that PET degradation involves two separable physiological layers: a general metabolic response to PET-derived carbon shared across cell phenotypes, and a biofilm-specific transport response potentially linked to substrate capture at the plastic interface. This work advances our understanding of how microbial physiology is organized during plastic biodegradation and identifies transport processes as previously unrecognized components of PET-degrading biofilms.

PET biodegradation↗

PubMed2026

Single-cell transcriptomic atlas of Alzheimer's disease middle temporal gyrus reveals region, cell type, and sex specificity of gene expression with novel genetic risk for MERTK in female.

BackgroundAlzheimer's disease (AD), the most common age-related neurodegenerative disease, is closely associated with both amyloid-β plaque and neuroinflammation. Two thirds of AD patients are female, and they have a higher disease risk; women with AD have more extensive brain histological changes than men along with more severe cognitive symptoms and neurodegeneration.ObjectiveThis study aimed to determine how sex difference induces structural brain changes and molecular cell vulnerabilities in AD, with a focus on identifying sex-specific transcriptional alterations and genetic risk factors.MethodsWe performed single nucleus RNA sequencing on postmortem brains from individuals with AD and age- and sex-matched controls, focusing on the middle temporal gyrus, a cortical brain region strongly affected by the disease, and integrated single nucleus RNA sequencing results with genome-wide association study (GWAS) data using cell type-specific enrichment and generalized gene-set analysis approaches. The analysis pipeline is provided with threshold information.ResultsWe identified a selectively vulnerable subpopulation of layer 2/3 excitatory neurons that were RORB-negative and CDH9-expressing in both males and females. Disease-associated, but sex-independent, reactive astrocyte signatures were also present. In clear contrast, the microglia signatures of AD brains differed between males and females. Integrating single cell transcriptomic data with results from GWAS, we identified MERTK genetic variation as a candidate novel risk factor for AD selectively in females.ConclusionsTaken together, our single cell atlas of middle temporal gyrus revealed a unique cellular-level view of sex-specific transcriptional changes in AD, illuminating GWAS identification of sex-specific AD genes. These data serve as a rich resource for interrogation of the molecular and cellular basis of AD.

Alzheimer's disease↗