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 37 records · Page 2Linked to original sources

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

BACH1 orchestrates macrophage state transitions to coordinate regenerative inflammation.

Efficient tissue regeneration requires the precise coordination of inflammatory and regenerative programs, principally mediated by monocyte-derived macrophages. However, the transcriptional wiring and epigenomic processes behind complex macrophage subtype specification and transition between the different states are not known. Here we have identified the transcriptional repressor BACH1 as a critical, cell-intrinsic regulator of monocyte-derived macrophage specification during skeletal muscle regeneration. Using a myeloid-specific BACH1 knockout mouse model, we demonstrate that BACH1 deficiency disrupts the temporal coordination of monocyte-to-macrophage differentiation, leading to aberrant macrophage subsets with concurrent opposing pro- and anti-inflammatory features. Single-cell RNA-sequencing profiling reveals that BACH1 controls a core transcriptional network, including Nfkb1, Cebpb, and interferon signaling, governing inflammatory resolution and functional macrophage specialization. Mechanistically, BACH1 loss accelerates macrophage differentiation but also affects its core cellular identity, resulting in sustained, rather than declining inflammatory programs including upregulation of Il1b and thus, defective tissue remodeling. These immune alterations compromise the paracrine landscape during regenerative inflammation and impair muscle stem cell differentiation. Our findings establish BACH1 as a molecular tuner or controller that integrates early innate immune signaling with regenerative output, positioning it as a central node linking transcriptional control, immune fate decisions, and tissue repair.

Animals↗

PubMed2026

TCGA-based identification of prognostic biomarkers and candidate traditional Chinese medicine compounds in papillary thyroid carcinoma: An observational study.

This study aimed to identify prognostic genes associated with papillary thyroid carcinoma (PTC) and explore candidate traditional Chinese medicine (TCM) compounds using integrated bioinformatics and molecular docking. In this observational study, PTC gene expression profiles and clinical data were obtained from The Cancer Genome Atlas. Differentially expressed genes were screened using differential-expression sequencing (DESeq2), followed by protein-protein interaction network analysis to identify hub genes. Their expression, diagnostic value, immune relevance, prognostic significance, protein-level validation, and single-cell distribution were assessed using gene expression profiling interactive analysis, receiver operating characteristic analysis, immune infiltration analysis, Kaplan-Meier survival analysis, the human protein atlas, and single-cell RNA-sequencing data. Candidate TCM compounds were predicted using symptom mapping (SymMap) and the TCM Systems Pharmacology Database and Analysis Platform, and molecular docking was performed to evaluate potential ligand-target interactions. Five hub genes, colony-stimulating factor 2, apolipoprotein E, fibronectin 1 (FN1), collagen type I alpha 1 chain (COL1A1), and intercellular adhesion molecule 1, were identified and found to be significantly upregulated in PTC tissues, with diagnostic value in receiver operating characteristic analysis. Immune infiltration analysis showed associations with macrophages, dendritic cells, and T helper 1 cells, whereas single-cell analysis demonstrated heterogeneous expression across immune and stromal cell populations, including fibroblasts. Higher FN1 and COL1A1 expression was associated with poorer outcomes. Immunohistochemistry supported the expression patterns, while single-cell analysis provided exploratory cell-type-level context for the cellular distribution of selected genes. Ginseng and Smilax glabra were predicted as common candidate TCMs, and docking suggested favorable binding between their active compounds and selected hub targets. Colony-stimulating factor 2, apolipoprotein E, FN1, COL1A1, and intercellular adhesion molecule 1 may be biologically relevant hub genes in PTC, while FN1 and COL1A1 may have prognostic value. Predicted TCM compounds provide preliminary computational evidence for possible compound-target interactions, requiring experimental and clinical validation.

Female↗

PubMed2026

Multiomics approaches reveal direct NF-κB p65 target genes in pancreatic islets during cytokine exposure and in type 1 diabetes.

Autoimmune diseases, including Type 1 diabetes (T1D), are often characterized by overactive inflammatory signaling pathways. The proinflammatory cytokine interleukin-1β (IL-1β) elicits global gene expression changes in islet β-cells which overlap with islets obtained from human donors with T1D. The direct transcriptional link between NF-κB subunit p65 and target genes involved with autoimmune events was investigated. We used a multiomics approach including bulk RNA-sequencing (RNA-Seq), single-cell RNA-sequencing (scRNA-Seq), and chromatin immunoprecipitation coupled to deep sequencing (ChIP-Seq), alongside molecular docking simulations, and transcriptional assays. Through the various experimental modalities, we identified early response genes driven by IL-1β that were differentially expressed in pancreatic islets from human T1D donors and also conserved across mouse, rat, and human tissues. ChIP-Seq revealed genes that are direct genomic targets of the NF-κB p65 transcription factor. Moreover, regions that gained RNA polymerase II binding following cellular exposure to IL-1β were identified, complementing the early response gene profile induced by β-cell exposure to IL-1β. Molecular docking simulations predicted that mutations reducing p65 transcriptional capacity do not alter DNA binding ability. These findings clearly show that IL-1β signaling in pancreatic β-cells directs p65 to specific genomic regions congruent with increased gene expression relevant to T1D in β-cell lines as well as mouse and human islets exposed to cytokines. Islets from human donors with T1D express genes identified as direct p65 targets using unbiased approaches, implicating heightened NF-κB activity as a critical component of autoimmune disease etiology.NEW & NOTEWORTHY Using multiple Seq-based approaches, this study identified genes expressed in human pancreatic tissue from donors with Type 1 diabetes that are regulated acutely by exposure to the cytokine interleukin-1beta. The NF-kB transcription factor p65 (RelA) was determined via ChIP-Seq to be a major control node regulating this immediate early response. These collective datasets are consistent with a paradigm of overactive NF-kB signaling as a critical component of autoimmunity in both rodents and humans.

Humans↗

PubMed2026

BET inhibition blunts antibody production and macrophage-mediated fibrosis to restore lung function in murine cGVHD.

Chronic graft-versus-host disease (cGVHD) significantly contributes to late mortality after allogeneic stem cell transplantation, with bronchiolitis obliterans syndrome (BOS) being a particularly lethal and treatment-resistant complication despite available therapies. Bromodomain and extraterminal (BET) proteins are epigenetic readers driving inflammatory transcriptional programs across multiple cell types. We hypothesized that BET inhibition would suppress inflammatory T and B cells and decrease macrophage polarization to a profibrotic phenotype, alleviating disease. In an established BOS cGVHD model, BET inhibition reduced germinal center (GC) formation and responses through a reduction of the CXCL13:CXCR5 axis and inflammatory T follicular helper/GC B cells in the spleen, along with a reduction in plasma cell infiltration within the lung. Mice with cGVHD had elevated pathogenic immunoglobulin G1 (IgG1) and IgM levels, both in circulation and deposited on lung tissue, which were attenuated under BET inhibition. Single-cell RNA-sequencing analysis revealed distinct cell states in the BOS lung vs control. In cGVHD mice, gene set enrichment analysis revealed the upregulation of profibrotic Arginase1 and Tgfb1 expression in alveolar macrophages (AM) and interstitial macrophages (IM), which was significantly reduced with BET inhibition. Furthermore, BET inhibition targeted lung-infiltrating M2 macrophages through the selective depletion of CD206+FcγR+ IM and AM, ultimately resulting in reduced collagen deposition and improved lung function. Our findings reveal a previously unrecognized mechanistic axis of BET regulation during cGVHD fibrosis and highlight BET inhibition as a promising therapeutic strategy.

Animals↗

PubMed2026

Matrine Alleviates Sepsis-Induced Acute Lung Injury by Reinforcing NQO1/SLC7A11/GPX4-Associated Anti-Ferroptotic Defenses and Attenuating NF-κB-Driven Inflammation.

BACKGROUND: Sepsis triggers dysregulated systemic inflammation and multiple-organ dysfunction, with the lungs being particularly susceptible to injury. Sepsis-induced acute respiratory distress syndrome (ARDS) is mainly driven by TLR4/NF-κB-mediated hyperinflammation and alveolar macrophage activation. Matrine, a bioactive alkaloid derived from Sophora flavescens, has been reported to modulate redox homeostasis and ferroptosis-associated lipid peroxidation. However, the target-specific mechanisms underlying its effects on ferroptosis and inflammatory signaling in sepsis-induced acute lung injury (SALI) remain incompletely understood. PURPOSE: This study aimed to evaluate the therapeutic effects of matrine in a cecal ligation and puncture (CLP)-induced SALI model and to determine whether its protective effects involve reinforcement of NQO1/SLC7A11/GPX4-associated anti-ferroptotic defenses and suppression of NF-κB-driven inflammation. METHODS: We analyzed the single-cell RNA-sequencing (scRNA-seq) dataset GSE273924 to characterize CD45-enriched pulmonary immune-cell subsets in sham mice and mice with intratracheal Escherichia coli-induced pneumonia. Network pharmacology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to predict Kushen (KS)-related targets and pathways associated with SALI. Differential expression analysis and weighted gene co-expression network analysis (WGCNA) of GSE245013 were used to identify candidate targets. Matrine-NQO1 binding and intracellular target engagement were evaluated using molecular docking, molecular dynamics simulations, surface plasmon resonance (SPR), and the cellular thermal shift assay (CETSA). The therapeutic effects of matrine were assessed in mice with CLP-induced SALI and in lipopolysaccharide (LPS)-stimulated MH-S cells. Lung histopathology, inflammatory cytokine production, target protein expression, ferroptosis-associated indicators, and NF-κB activation were evaluated using molecular, biochemical, and histological assays. The functional contribution of NQO1 was further examined using the NQO1 inhibitor ES936. RESULTS: scRNA-seq analysis of GSE273924 revealed substantial remodeling of the CD45-enriched pulmonary immune-cell landscape in mice with intratracheal E. coli-induced pneumonia, including macrophage transcriptional programs associated with ferroptosis and inflammatory signaling. Integrated network pharmacology and bioinformatics analyses prioritized NQO1 as a candidate target of matrine and identified NF-κB signaling as a potentially relevant pathway. Molecular docking, molecular dynamics simulations, SPR, and CETSA supported matrine-NQO1 binding and intracellular target engagement. Functionally, matrine improved survival, attenuated lung injury, reinforced NQO1/SLC7A11/GPX4-associated anti-ferroptotic defenses, and suppressed NF-κB activation in CLP mice. Similar protective effects were observed in LPS-stimulated MH-S cells. ES936 partially attenuated the matrine-mediated improvements in cell viability, redox homeostasis, ferroptosis-associated indicators, and NF-κB p65 phosphorylation, supporting a functional contribution of NQO1 to the protective effects of matrine. CONCLUSION: Matrine alleviates SALI by reinforcing NQO1/SLC7A11/GPX4-associated anti-ferroptotic defenses and attenuating NF-κB-driven inflammation.

Animals↗

PubMed2025

Intrinsic changes in cell differentiation and identity drive impaired wound healing in aged female murine skin.

Cellular and molecular mechanisms that drive a perturbed wound microenvironment and impaired healing in aged skin have not been fully delineated. To obtain a comprehensive understanding of cell-intrinsic changes acquired during ageing that impact early responses to injury, we performed single-cell RNA sequencing in young and aged intact female murine skin and wounds 3 days post-injury. We observed that substantial changes in the mean proportional distribution and transcriptomic state of skin resident subpopulations in aged, but not young, tissues accompany a global increase in basal inflammation. This is driven by an altered signalling environment leading to impaired keratinocyte differentiation, loss of fibroblast identity and defective macrophage function. Further, we show that ageing-induced changes in skin resident cells persist after injury, resulting in increased expression of senescence-related genes in wound fibroblasts and aberrant monocyte-to-macrophage transitioning coupled to an enhanced inflammatory signature and defective intercellular signalling in comparison to wounds in young mice. In summary, our data highlights a contribution of both cell-intrinsic changes and an altered tissue microenvironment to poor wound healing responses in aged mice.

Animals↗

PubMed2026

Malignant epithelial states drive immune dysfunction in ampulla of Vater carcinoma.

BACKGROUND: Ampulla of Vater (AoV) carcinoma is a rare malignancy arising at the junction of intestinal and pancreatobiliary epithelium. Its heterogeneous clinical behavior and histological diversity have hindered therapeutic advances, and the cellular basis of this heterogeneity remains unclear. We aimed to construct a single-cell transcriptomic atlas of AoV carcinoma, with a focus on identifying epithelial subtypes and their interactions with the tumor microenvironment (TME). METHODS: We performed single-cell RNA sequencing on eight primary AoV tumors and four matched normal tissues. Comprehensive clustering and transcriptomic analyses identified cell-type composition, epithelial heterogeneity, and tumor-immune interactions. Findings were validated using deconvolution of bulk RNA-seq data from 62 AoV carcinoma patients. Results Malignant epithelial cells were categorized into four distinct subtypes: Int-Wnt, PB-KRAS, Int-Hypoxia, and Cycling stage. PB-KRAS cells exhibited stem-like transcriptional programs and high genomic instability. Deconvolution analysis of bulk RNA-seq data from the independent AoV cohort revealed that enrichment of the PB-KRAS subtype correlated with tumor recurrence and poor survival. Our immune profiling analysis discovered a significant association between PB-KRAS subtype and GZMK+ CD8+ T cells, which are in a pre-dysfunctional state, alongside SPP1+ macrophages exhibiting immunosuppressive traits. Spatial transcriptome data further supports the immunosuppressive natures of TME around PB-KRAS subtype malignant epithelial cells in AoV carcinoma. CONCLUSIONS: Our study presents a single-cell atlas of AoV carcinoma, highlighting the molecular diversity of malignant epithelium and its association with the immune microenvironment. The PB-KRAS subtype emerges as a stem-like, immunosuppressive tumor state associated with poor prognosis, providing insights for future therapeutic targeting.

Ampulla of Vater carcinoma↗

PubMed2026

Machine learning-integrated multi-omics risk prediction for pulmonary fungal infection in COPD and lung cancer: a transcriptomic and immune profiling study.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) and lung cancer are major risk factors for invasive pulmonary fungal infection (IPFI), carrying an attributable mortality of 30%-80%. Their coexistence further amplifies immunosuppression, while current diagnostic criteria remain inadequate for early risk identification. METHODS: Transcriptomic data from the GEO dataset GSE296912 (scRNA-seq; 12,078 cells from normal and COPD lung tissue) and The Cancer Genome Atlas (TCGA)-lung adenocarcinoma (LUAD) bulk RNA-seq cohort (539 tumor and 59 normal samples) underwent differential expression and cross-omics integration analysis. Five machine learning models were constructed: logistic regression, SVM, random forest, XGBoost, and LASSO. Candidate genes were validated by qRT-PCR in A549 cells and THP-1-derived macrophages stimulated with heat-inactivated Aspergillus fumigatus conidia, a protocol selected to ensure BSL-2 biosafety compliance and isolate PAMP-mediated innate immune signaling. Model performance was evaluated using 5-fold stratified cross-validation with AUC, calibration curves, and decision curve analysis. RESULTS: Single-cell transcriptomic analysis of 12,078 cells identified 14 distinct cell populations, with marked myeloid expansion and immune dysregulation in COPD lung tissue. Cross-omics integration with TCGA-LUAD data identified 1,145 shared genes (79 immune-related), converging on NF-κB, TLR4, and cytokine receptor signaling. The random forest model achieved excellent discriminative performance (5-fold CV AUC = 0.988), with Treg infiltration, TLR4, and MMP9 as the top predictors. qRT-PCR confirmed significant upregulation of all five candidate genes (DEFB4A, S100A8, IL-8, MMP9, and TLR4) in both A549 and THP-1 cells following fungal stimulation. CONCLUSION: This multi-omics machine learning model integrating scRNA-seq and TCGA transcriptomic data demonstrates excellent discriminative performance (AUC = 0.988), with mechanistic convergence of NF-κB, TLR4, and oncogenic signaling pathways identified across shared immune gene signatures. In vitro qRT-PCR validation confirms the biological relevance of five key antifungal immune genes, providing a transcriptomic foundation for future prospective IPFI risk stratification in patients with COPD and lung cancer.

TLR4↗

PubMed2026

Follicular Lymphoma Transformation is Characterized by Cytokine-associated Remodeling of Stromal and Macrophage Compartments.

Across cancer, one of the most frequent examples of histologic transformation is the evolution of follicular lymphoma (FL) to an aggressive large cell lymphoma. Despite recent progress, understanding of the molecular and cellular underpinnings of transformation remains incomplete. Here, we dissect the interplay of tumor and microenvironment cell populations across transformation through a multimodal investigation of 95 FL and transformed FL (tFL) samples, including single-cell and bulk RNA-sequencing alongside spatial transcriptomics and proteomics, and validate findings across independent FL-tFL pairs. Upon transformation, fibroblasts and GPNMB+ macrophages increase while lymph-node organizing follicular dendritic and CCL21+ fibroblastic reticular cells were lost, resulting in an altered spatial distribution of cytokines that impacts T cell infiltration and macrophage differentiation and function. Secreted stromal and macrophage signals were further evident by non-invasive plasma proteomics. Taken together, our data reveal expansion of macrophages and fibroblasts as key features of transformation with potential diagnostic and therapeutic implications.

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

Unveiling tumor heterogeneity by single cell RNA-sequencing: From basic considerations to clinical applications.

Tumor heterogeneity-encompassing diverse cellular phenotypes, genomic alterations, and microenvironmental contexts-is a principal barrier to effective cancer therapy. Single-cell RNA sequencing (scRNA-seq) has transformed our ability to resolve this complexity by capturing transcriptomes at single-cell resolution. Here, we review the technical foundations required for high-quality scRNA-seq studies. We then trace the evolution of scRNA-seq platforms from manual micromanipulation to high-throughput systems, and describe the computational pipelines that enable reliable data interpretation. The application of scRNA-seq is exemplarily shown in the context of lung cancer, where single-cell profiling has revealed (i) the clonal and sub-clonal architecture of tumors, (ii) extensive remodeling of the immune microenvironment, iii) key mechanisms underlying resistance to targeted agents and immune-checkpoint blockade, and (iv) the dynamics of neo-antigen-specific T-cell responses. Integrating machine-learning techniques-such as deep-learning classifiers and graph-based models-with single-cell transcriptomic data has markedly sped up biomarker discovery, produced more accurate risk-stratification scores, and enabled the generation of patient-specific therapeutic predictions. We surveyed the major trial registry ClinicalTrials.gov and identified ∼380 ongoing or completed studies that explicitly incorporate scRNA-seq as a correlative or pharmacodynamic endpoint. Overall, the analysis shows that scRNA-seq becomes an increasingly important component of modern trials, providing high-resolution cellular and molecular readouts that complement conventional imaging and bulk-omics endpoints. While key challenges remain, ranging from costs, scalability and need for rigorous validation before routine clinical deployment, ongoing technological advances continue to expand the potential of scRNA-seq as a cornerstone of precision medicine.

Humans↗

PubMed2025

Colorectal Liver Metastasis Pathomics Model: Integrating Single-Cell and Spatial Transcriptome Analysis With Pathomics for Predicting Liver Metastasis in Colorectal Cancer.

The liver is the primary target organ for hematologic metastasis of colorectal cancer (CRC), and CRC liver metastasis (CRLM) often precludes radical resection, making it the leading cause of death in patients with CRC. To improve the identification and prediction of liver metastasis risk, we identified a cell type of liver metastasis--triggering malignant cells (LMTMCs) through integrating single-cell RNA sequencing and spatial transcriptome analysis. Multiomics cell communication analysis indicated that the interaction between fibroblasts and LMTMCs through the COL1A1-CD44/SDC4 and LAMA4-CD44 signaling axes could promote CRLM. By applying the one-class logistic regression algorithm, we developed a CRLM scoring system in the bulk RNA-sequencing data according to the abundance of LMTMCs in each individual. Using the grouping labels derived from the CRLM scoring system in the bulk data and the corresponding whole-slide images without any manual annotations at the region or pixel level, processed via slide-level weakly supervised learning, a deep-learning model based on the ResNet18 architecture, called Colorectal Liver Metastasis Pathomics Model, was developed to predict the risk of liver metastasis in patients with CRC. The Colorectal Liver Metastasis Pathomics Model achieved an area under the curve of 0.84 at the internal test set of The Cancer Genome Atlas-CRC histology images. In the external independent validation sets, namely the Affiliated Hospital of Southwest Medical University and the Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University cohorts, the areas under the curve were 0.89 and 0.72, respectively, indicating effective classification performances. This study provided new insights and tools for the early identification of CRLM and demonstrated the potential of combining multiomics with deep learning-based pathomics in cancer research.

Humans↗

PubMed2026

An integrated single-cell and spatial transcriptomic atlas of thyroid cancer progression identifies prognostic fibroblast subpopulations.

Although well-differentiated thyroid carcinoma (WDTC) is characterized by a robust treatment response, aggressive subtypes, such as anaplastic thyroid carcinoma (ATC), remain highly lethal. To understand thyroid cancer evolution in both children and adults, we analyzed single-cell transcriptomes of 423,733 cells from 81 samples and spatially resolved key tumor and microenvironment populations across 28 tumors with spatial transcriptomics, including rare and unique composite WDTC/ATC tumors and pediatric diffuse sclerosing thyroid carcinomas. Additionally, we identified gene signatures of stromal cell populations in 5 large thyroid cancer bulk RNA-sequencing cohorts. Through this multi-institutional effort, we defined a population of POSTN+ myofibroblast cancer-associated fibroblasts (myCAFs) that are intimately associated with invasive tumor cells and correlate with poor prognosis, lymph node metastasis, and disease progression in thyroid carcinoma. We also revealed a population of inflammatory CAFs that are distant to tumor cells and are found in the inflammatory stromal microenvironment of autoimmune thyroiditis. Together, our study provides spatial profiling of thyroid cancer evolution in samples with mixed WDTC/ATC histopathology and identifies a prognostic myCAF subtype with potential clinical utility in predicting aggressive disease in both children and adults.

Humans↗

PubMed2026

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

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

Gene expression profiling↗

PubMed2026

OscillomeR infers ultradian oscillations and targets of the Hes family.

The Hes family, basic-helix-loop-helix transcription factors and downstream effectors of Notch signaling, regulate the fate choices of pancreatic progenitors, muscle stem cells, neuronal progenitors, and presomitic mesoderm cells. Bioluminescence imaging (BLI) has revealed ultradian oscillatory dynamics of Hes-family members Hes1, Hes5, and Hes7. However, identifying which of the Hes target genes also oscillate remains challenging due to the time-consuming and costly nature of tracking individual target genes using BLI. Here, we propose OscillomeR, a computational framework that reconstructs ultradian oscillations from RNA-sequencing data to identify oscillatory target genes at high throughput. OscillomeR predicts thousands of oscillatory genes in synchronized or unsynchronized cell types, identifying both known and novel Hes-family targets. It also captures the dynamic rewiring of gene-regulatory networks during cell differentiation. Overall, OscillomeR is an effective tool for elucidating the functions of oscillatory transcription factors at the genomic scale.

Hes family↗

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

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↗