SEARCH · PubMed Health

Results for “Single-cell RNA sequencing”

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

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

At least 73 records · Page 4Linked to original sources

PubMed2026

Identification of molecular subtypes in clear cell renal cell carcinoma based on chromatin regulators and tumor immune microenvironment profiling.

In the histological classification of renal cell carcinoma, clear cell renal cell carcinoma (ccRCC) accounts for the highest proportion and is the most common subtype. Despite advances in management, it continues to be associated with considerable incidence and mortality. Although surgery and systemic therapies are available, their efficacy is constrained by pronounced intratumoral heterogeneity and treatment resistance. Identifying robust biomarkers and clarifying the underlying biological mechanisms are therefore essential to improving diagnosis, risk stratification and therapeutic decision-making. In this work, we identified two ccRCC molecular subtypes displaying divergent chromatin regulator (CR) profiles and different clinical prognoses. Using the genes differentially expressed between these subgroups, we constructed a CR-related score (CRS) that effectively stratified patients according to survival. More analysis concluded that the low expression of CR was more linked with the immune-activated tumors, which encompassed the immune pathway enrichment, as well as the elevation of numerous immune cell subtypes. Moreover, elevated CRS was associated with improved immunotherapy responsiveness. Drug-sensitivity analyses nominated several candidate agents, and SMARCD3 knockdown in 786-O cells inhibited proliferation and migration and reduced sensitivity to masitinib. Collectively, these findings support the prognostic and therapeutic relevance of CR-related states in ccRCC and provide a framework for future experimental validation of chromatin-regulated tumor-immune interactions.

Humans↗

PubMed2026

Genome-wide association study reveals two novel genetic loci associated with chronic lung allograft dysfunction.

BACKGROUND: Chronic lung allograft dysfunction (CLAD) leads to declining respiratory function and high mortality, representing the main barrier to long-term survival in lung transplantation (LT). We performed the first genome-wide association study (GWAS) investigating donor's and recipient's genetic factors associated with CLAD. METHOD: We genotyped 392 donor-recipient pairs from the multicentric Cohort in Lung Transplantation. We tested 4.5 million SNPs for association with CLAD using multivariable logistic regression models corrected for age, sex, initial disease and genetic ancestry. Three levels of explanatory variables were separately considered to conduct GWAS: donors-only, recipients-only, and donor-recipient mismatches. We also ran HLA-centric analyses using the same models. RESULTS: Our analysis confirmed the deleterious impact of HLA allelic and epitopic mismatches on CLAD risk, mostly driven by class I HLA (p=0.004). No significant associations with CLAD were found for donors' genotypes or donor-recipient non-HLA mismatches. We highlighted two independent recipient's loci associated with CLAD, including one protective signal (0.39 in CLAD vs 0.66 in non-CLAD recipients, p-value=5.05×10-7, q-value=0.017, OR=0.35) encompassing the PLXDC2 gene, and one risk signal (0.66 in CLAD vs 0.38 in non-CLAD recipients, p-value=9.86×10-7, q-value=0.017, OR=2.83) encompassing the ZNF518A/BLNK genes. These non-coding SNPs are putative regulatory variants of gene expression. Importantly, our single-cell RNA-sequencing showed a down-regulation of PLXDC2 in fibroblasts and lung epithelium in CLAD vs healthy controls. CONCLUSION: This first LT GWAS revealed two candidate loci from the recipient's genome, both biologically relevant for CLAD pathogenesis. Our study calls for larger LT genomic initiatives to increase power for signal discovery.

Humans↗

PubMed2026

A machine learning-derived and functionally validated circadian rhythm signature predicts clinical outcomes and in silico drug sensitivity in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) displays considerable heterogeneity in clinical outcomes, highlighting the need for reliable prognostic biomarkers. While the aberrant expression of circadian rhythm-related genes has been implicated in cancer pathogenesis, its comprehensive role in CRC progression and predicted therapeutic vulnerabilities remains inadequately characterized. METHODS: Bulk and single-cell RNA-sequencing data were integrated from multiple CRC cohorts. A circadian rhythm signature (CRS) was developed through machine learning algorithms and validated for prognostic value. Comprehensive analyses of tumor microenvironment, genomic alterations, and drug sensitivity were performed. Furthermore, the biological function of the core gene, BHLHE40, was validated in CRC cell lines through CCK-8, EdU, and wound healing assays. RESULTS: Single-cell analysis demonstrated an elevated expression signature of circadian rhythm-related genes in dendritic cells. The optimized CRS, comprising 14 circadian rhythm-related genes, successfully categorized patients into high- and low-risk groups. Patients with a high CRS showed markedly poorer overall survival and computationally inferred immunosuppressive features, including reduced CD8+ T cell infiltration and increased M2 macrophage polarization. Genomic analysis revealed enhanced mutation burden in TP53 and alterations in RTK-RAS/WNT pathways. Notably, in vitro assays confirmed that BHLHE40 is significantly overexpressed in CRC cells. Knockdown of BHLHE40 markedly inhibited tumor cell proliferation and migration. Drug sensitivity profiling identified bexarotene and SMER-3 as potential therapeutic options for high-CRS patients. A nomogram integrating CRS with clinical parameters demonstrated superior predictive accuracy for 1-, 3-, and 5-year survival. CONCLUSIONS: The CRS represents a promising prognostic biomarker that reflects tumor immune status and genomic features, providing valuable insights for personalized treatment strategies in CRC.

Circadian rhythm↗

PubMed2026

Single-cell transcriptome revealed the aberrant keratinocytes activation in antigen presentation in atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD), a common chronic inflammatory skin disease, has been extensively studied using single-cell genomics. However, keratinocytes, as key effector cells in AD, have underlying mechanisms remain incompletely understood and require further investigation. METHODS: We integrated single-cell transcriptomic data from skin tissues of healthy controls, chronic active AD patients, spontaneously healed AD (SHAD) patients, and an ovalbumin-induced AD mouse model. The study particularly emphasized the gene expression and cellular dynamics of keratinocytes across the different groups, as well as their interactions with immune cells. RESULTS: Compared to healthy controls, we observed significant changes in the keratinocyte transcriptome, cellular state, and keratinocyte-immune cell ligand-receptor interactions in AD skin, particularly the marked activation of genes involved in antigen processing and presentation. Interestingly, such gene activation was not observed in keratinocytes from the ovalbumin-induced AD mouse model, despite its phenotype closely resembling human AD. Furthermore, in SHAD, we identified a recovery of both the ligand-receptor interaction patterns and antigen processing and presentation genes, accompanied by a notable shift in the transcriptome. This involved a significant downregulation of genes related to cytoplasmic transcription and oxidative phosphorylation. Notably, this pattern was not observed in the self-healing mouse model following the removal of ovalbumin stimulation. CONCLUSION: Our results suggest that the persistent activation of antigen processing and presentation pathways in keratinocytes may be a key driver of chronic inflammation in AD. Therefore, redirecting anti-allergic therapeutic strategies from solely targeting immune cells to targeting of keratinocyte-mediated antigen presentation may offer a more effective approach. Furthermore, we raise concerns about the use of ovalbumin-induced mouse models to recapitulate human chronic AD, as the underlying mechanisms may differ significantly.

Dermatitis, Atopic↗

PubMed2026

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase↗

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

The establishment of prostate-specific, SKP2 humanized mice by CRISPR knock-in method reveals neoplastic initiation and microenvironmental reprogramming.

Genetic inactivation of SKP2 has been shown to effectively prevent cancer initiation and block tumorigenesis. However, direct in vivo evidence for SKP2 on cancer initiation and prostatic microenvironment is still lacking and a SKP2 humanized mouse model is critical for developing prostate cancer immunoprevention approaches through targeting SKP2. We therefore have established a prostate-specific human SKP2 knock-in mouse model driven by an endogenous mouse probasin promoter. Overexpression of hSKP2 induces PIN and low-grade carcinoma. RNA-sequencing analysis revealed significant gene expression alterations in EMT, extracellular matrix, and interferon signaling. Single-cell deconvolution showed an increase of fibroblast population and a decrease of CD8+ T cell and B cell populations. Consistent with these results from the SKP2 humanized mouse, SKP2 protein is overexpressed in human prostatic hyperplasia, PIN and prostate adenocarcinoma compared to normal prostate tissues. Overexpression of SKP2 markedly increased cell migration and invasion and induced the gene expression of EMT and interferon pathways. Inhibition of SKP2 signaling by Flavokawain A and C1 reverses EMT and affects EMT and interferon-related gene expression. In addition, paired prostate organoids were derived from SKP2 humanized and wild-type mice for drug screening and validated by known SKP2 inhibitors, Flavokawain A and C1. Both of which selectively decreased viability and altered the morphologies of organoids of hSKP2 knock-in rather than wild-type mice. Our studies provide a well-characterized prostate-specific hSKP2 knock-in mouse model and offer new mechanistic insights for understanding the oncogenic role of SKP2 in shaping the prostatic microenvironment during early carcinogenesis.

Animals↗

PubMed2026

Full-length single-cell spatial transcriptomics reveals spatial and cell-type-specific transcript isoforms in the primate brain.

The primate brain exhibits complex RNA alternative splicing heterogeneity crucial for functional complexity, yet systematic spatial isoform characterization has been lacking. We developed Fullscope-seq, a full-length single-molecule large field-of-view spatial transcriptomics sequencing method at single-cell resolution, based on programmed concatenation cDNA for multiple long-read sequencing platforms. Applying Fullscope-seq to the macaque brain, we uncovered thousands of genes exhibiting differential transcript usage (DTU) across cortical layers, cell types and brain regions. Fullscope-seq resolved hundreds of major isoform switches across distinct brain regions and identified DTUs between superficial and deep cortical layers. Cortical layer-specific DTUs showed cell-composition dependence, whereas regional DTUs were regulated according to both cellular composition and spatial contexts. These isoform variations showed substantial enrichment for neuropsychiatric disorder-associated genes and were conserved across platforms and species. Our study establishes a scalable framework for spatial isoform analysis and provides a resource for understanding transcriptomic diversity in complex tissues.

Animals↗

PubMed2026

Integrated bioinformatics analysis reveals cross-talking hub genes and therapeutic agents between sepsis and acute myocardial infarction.

BACKGROUND: Sepsis and acute myocardial infarction (AMI) are two significant diseases that may share overlapping etiological mechanisms. This study aims to systematically identify core genes common to both conditions and to explore their potential as therapeutic targets and drug candidates through an integrative analysis of clinical data and bioinformatics. METHODS: The AMI dataset was obtained from the GEO database, and RNA sequencing data were collected from blood samples of patients with sepsis at our hospital. Common genes were identified using differential expression gene analysis (DEG) and weighted gene co-expression network analysis (WGCNA). Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were performed. A protein-protein interaction (PPI) network was constructed, and hub genes were identified using the MCC/Degree algorithm. Diagnostic value was assessed via receiver operating characteristic curve analysis. Immune infiltration patterns, single-cell sequencing data, and molecular docking simulations were employed to evaluate immune relevance and identify potential therapeutic compounds. RESULTS: A total of 417 genes were identified between sepsis and AMI, with enrichment analysis revealing significant involvement in inflammatory responses. Three hub genes-JAK2, MYD88, and TIMP1-were selected for further investigation. ROC curves confirmed their strong diagnostic performance for both diseases. Immune infiltration analysis showed that these core genes were significantly correlated with the infiltration levels of various immune cell types. Molecular docking indicated that quercetin exhibited stable binding affinity with the proteins encoded by these genes. qPCR validation further confirmed the upregulation of these three genes, supporting the anti-inflammatory effects of quercetin as a potential targeted therapy. CONCLUSION: JAK2, MYD88, and TIMP1 were identified as shared core genes in sepsis and AMI. These genes not only serve as potential diagnostic biomarkers but also offer novel targets for developing common therapeutic strategies for both conditions. Furthermore, quercetin emerges as a promising candidate for targeted treatment.

Humans↗

PubMed2026

Receptor-defined targeting of a genomically unique melanoma-enriched noncanonical antigen.

Effective T cell-based immunotherapies require functional receptors that can be engineered and redeployed to recognize tumor-restricted antigens. Noncanonical peptides arising from transcription outside annotated protein-coding regions expand the antigenic landscape of cancer; however, systematic strategies to biologically prioritize and functionally validate such targets remain underdeveloped. Here, we integrated de novo transcript analysis, exon-resolved quantification, RNA in situ hybridization, and immunopeptidomics to identify melanoma-associated noncanonical transcripts and advance candidates through receptor-level validation. Among three recurrent melanoma-associated transcripts, EVA003 emerged as a lead target based on its distinct repeat-enriched genomic architecture, consistent tumor-enriched exon-level expression across independent datasets, and a genomically unique immunogenic core sequence. We demonstrate endogenous presentation of EVA003-derived peptides on HLA-A*03:01 and detect specific reactivity in patient-derived tumor-infiltrating lymphocytes. Single-cell transcriptomic profiling identified a dominant peptide-reactive clonotype, enabling isolation of a naturally occurring T cell receptor. Transfer of this receptor into healthy donor T cells conferred antigen-dependent activation and cytotoxicity against both peptide-pulsed targets and melanoma cells expressing EVA003 endogenously. Together, these findings establish a biologically informed strategy for prioritizing noncanonical tumor antigens and demonstrate that genomically unique, tumor-enriched noncanonical peptides can be presented to molecularly defined receptors capable of mediating cancer cell killing. These findings support the integration of prioritized noncanonical antigens into engineered T cell therapeutic strategies.

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

Cell-type signatures of Alzheimer's disease shared across population groups.

Genomic studies at single-cell resolution have identified several cell types associated with clinical and pathological traits in Alzheimer's disease1-9, but have not examined associations that are shared across populations. To bridge this gap, here we use single-nucleus RNA sequencing and assay for transposase-accessible chromatin with sequencing to profile cortical and subcortical regions in post-mortem brain-tissue samples from Latin, white (excluding Latin) and African American (excluding Latin) individuals. Using discrete and continuous dissections of molecular programs, we identify cell-type-specific clusters associated with Alzheimer's disease in a region-specific manner across all three population groups, including microglial (GPNMB+ and CD74+ subgroups), astrocytic (SERPINH1+, CD44+ and WIF1+ subgroups) and neuronal (SST+ GABAergic and superficial-layer glutamatergic) signatures. We also report continuous gene-expression factors in astrocytes and oligodendrocytes that are not captured by discrete cluster assignments, but which show strong associations with disease phenotypes; these factors are enriched for genes associated with annotated functions such as lipid processing and neurotransmitter reuptake. Finally, we find that molecular programs reveal six distinct subgroups of individuals with cognitive impairment that span all three populations, are not captured by neuropathology, and are instead distinguished by molecular signatures that are not universally present but are nonetheless associated with ante-mortem impairment. Overall, our study identifies key cell types and gene programs implicated in Alzheimer's disease that are shared across population groups, and underscores how representative sampling can capture both shared signatures and disease heterogeneity, thereby enabling better prioritization of key cell types for further investigation.

Female↗

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

Integrative quantum and systems biology of cancer: From molecular fluctuations to ecological outcomes.

This review treats cancer as a multiscale adaptive system, asks what the framework must predict to be worth adopting, and separates at each scale what the evidence establishes from what is proposed. It is an expert narrative synthesis, not a systematic review, and states the limits of that design. Proton transfer and tautomeric shifts contribute to spontaneous mispairing but do not license claims of directed or non-random mutation: replication timing, three-dimensional chromatin organization, sequence context and known mutagenic processes explain most mutational heterogeneity, leaving any quantum contribution as a residual against that baseline. The Waddington quasi-potential is bounded: outside detailed balance the dynamics are not gradient-derivable and require a probability-flux term. Hysteresis, rate-limited bimodality and return to state after perturbation distinguish an attractor from a transcriptomic cluster. Single-cell karyotype and live-imaging evidence supports whole-genome doubling as an unstable intermediate of heterogeneous origin and context-dependent consequence, not a uniform adaptive strategy. Systems and synthetic biology, virtual cells and digital twins are assessed against benchmarks, not promise. Tissue-scale ecology is reported with the spatial measurements now quantifying it, including evidence that stromal niche construction is not uniformly tumor-supporting. RNA modification is a layer in its own right, showing that the interpretation of a regulatory signal, not its magnitude, is biologically decisive. A dedicated section states the framework's commitments, the observable and evidence at each scale, and what would falsify them, asking what this adds to somatic mutation theory with clonal evolution and plasticity.

Neoplasms↗

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

Single-nucleus profiling of postmortem diffuse midline gliomas identifies mitochondrial biogenesis as a resistance mechanism to imipridone therapy.

BACKGROUND: Imipridone ONC201 is the first FDA-approved therapy for H3K27-altered diffuse midline glioma; however, clinical responses remain limited. Defining tumor-intrinsic determinants and microenvironmental, extrinsic factors that shape sensitivity or resistance to imipridones will identify actionable therapeutic opportunities and inform improved clinical strategies. METHODS: To identify mechanisms of imipridone resistance, we obtained postmortem brain tissue from DMG patients who had received imipridones and/or standard care. Single-nucleus RNA and open-chromatin sequencing were performed on N = 22 cases. Immunofluorescence-based myeloid phenotyping was performed on N = 46 cases. Mitochondrial copy-number analysis was performed on N = 19 cases. Validation of imipridone sensitivity, its effect on mitochondrial density, and its synergy with inhibition of mitochondrial biogenesis were assessed in DMG primary cells. RESULTS: We established a single-cell RNA/open-chromatin atlas from postmortem DMG cases and found imipridone treatment resulting in regressed mesenchymal transition, reduced myeloid-derived suppressive cells, and reversed aberrant H3K27-altered enhancer activity. Resistant tumors showed increased mitochondrial density, turnover, and membrane potential. Mitochondrial biogenesis and PPARGC1A emerged as resistance biomarkers and actionable targets. CONCLUSIONS: These studies implicate mitochondrial biogenesis as a biomarker of imipridone resistance and a focus for the development of combinatorial strategies to provide effective therapeutic options for a challenging pediatric brain tumor.

Humans↗

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↗