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 19 recordsLinked to original sources

PubMed2026

Immune-Like Malignant Epithelial Programs Shape Tumor-Immune Interactions and Inform Prognostic Stratification in Lung Adenocarcinoma.

Lung adenocarcinoma (LUAD) is characterized by marked cellular heterogeneity, yet how malignant epithelial states contribute to immune regulation and clinical outcomes remains incompletely defined. We integrated single-cell RNA-sequencing data to map the cellular landscape of LUAD and identify malignant epithelial cells based on inferred copy-number alterations. Epithelial states were further examined through trajectory inference, transcription factor analysis, and cell-cell communication profiling. Single-cell-derived genes were subsequently integrated with TCGA and independent GEO cohorts to construct and validate a machine learning-based prognostic signature. Malignant epithelial cells displayed distinct functional programs, including an immune-like state associated with genomic instability, immune-related transcriptional activity, tumor-immune communication, and patient outcomes. The resulting immune-like malignant epithelial cell signature (IMEC-Sig) consistently stratified survival across multiple cohorts. Low IMEC-Sig scores were accompanied by greater immune infiltration, higher immune checkpoint expression, and increased immunophenoscore, whereas high scores were linked to a comparatively immunosuppressive phenotype. Pan-cancer analyses further identified KRT8 as a gene associated with unfavorable prognosis, and functional experiments showed that KRT8 silencing suppressed proliferation, migration, invasion, and colony formation in LUAD cells. Together, these findings connect malignant epithelial heterogeneity with the immune context and clinical outcomes, support IMEC-Sig as a biologically informed prognostic tool, and nominate KRT8 as a potential therapeutic target in LUAD.

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

Ligand-based directed differentiation to produce granulosa-like cells expressing steroidogenic enzyme genes.

The ovarian granulosa cells are responsible for producing hormones and supporting oocytes through maturation and meiotic resumption. There is a need to generate granulosa-like cells (GLCs) from human induced pluripotent stem cells (hiPSCs) to better model human gonadal development and to test the effects of exogenous or pharmaceutical compounds on the ovary. Here we report a rapid ligand-based protocol for differentiating hiPSCs into cells that express markers of the transient developmental lineages and steroidogenic pathway genes. Single-cell RNA-sequencing (scRNA-seq) analysis identified canonical granulosa cell genes were expressed in a subset of cells and identified new genes of interest that were significantly associated with computationally modeled pseudotime. HSD17B1 was expressed in resulting GLCs but at low levels, suggesting an immature granulosa cell phenotype. The GLCs were produced using a simple culture method that could be augmented for granulosa cell functions such as sustaining oocyte growth. Producing GLCs through protocols such as this one is a first step toward designing large-scale ovarian endocrinology assays and developing personalized cell-based fertility and hormone restoration technologies in the future. This rapid protocol produced cells that express steroidogenic enzyme genes etoc blurb. Kubo and colleagues present a 5-day rapid protocol to generate immature granulosa-like cells from hiPSCs. Cells differentiated with inhibition of DKK1, a WNT signaling target gene, expressed gonadal ridge markers and FOXL2 transcripts and protein. Additionally, steroidogenic enzyme genes were expressed. A small population of differentiated cells were identified as expressing early-stage granulosa cell genes by single-cell RNA-seq.

Female↗

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

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

Beta cell-derived cholecystokinin drives obesity-associated pancreatic adenocarcinoma development.

Pancreatic endocrine-exocrine crosstalk plays a key role in normal physiology and disease and can be altered by host metabolic states, such as obesity. Classically, endocrine islet beta (β) cell secretion of insulin is thought to promote the development of obesity-associated pancreatic adenocarcinoma (PDAC), an exocrine cell-derived tumor. Here, we show that β cell expression of the peptide hormone cholecystokinin (CCK) is necessary and sufficient for obesity-associated PDAC progression in mice and that CCK expression - rather than insulin - correlates strongly with enhanced tumorigenesis. Single-cell RNA-sequencing, in silico latent-space archetypal and trajectory analysis, and experimental lineage tracing in vivo reveal that obesity induces the expansion of postnatal immature β cells, which adapt to express CCK via stress-responsive JNK/cJun signaling. Finally, obesity perturbs CCK-dependent peri-islet exocrine cell transcriptional states and enhances islet-proximal tumor formation. These results define endocrine-exocrine CCK signaling as a bona fide driver of obesity-associated PDAC development and uncover avenues to target the endocrine pancreas to subvert exocrine tumorigenesis.

Animals↗

PubMed2026

shinyDeepGxP: a user-friendly R shiny app for predicting surface protein abundance from scRNA-seq expression using deep learning in blood cells.

MOTIVATION: Understanding accurate immune cell heterogeneity and function in single-cell datasets requires access to protein-level information, which is often unavailable due to experimental limitations. RESULTS: We present shinyDeepGxP, an interactive web application featuring our deep learning model, DeepGxP, for predicting surface protein abundance from single-cell RNA-sequencing (scRNA-seq) data. This platform makes DeepGxP accessible to researchers without programming skills. Users can upload scRNA-seq count matrices and use "Predict Protein" to predict the abundance of 224 biologically relevant surface proteins. shinyDeepGxP provides visualizations to help identify distinct cell populations based on predicted protein profiles. Moreover, users can choose "Explore Model" to reveal key RNA predictors and their associated biological pathways for each protein. Overall, shinyDeepGxP is a user-friendly, freely available web tool that provides protein-level detail for RNA-only single-cell datasets, enabling multimodal discovery without additional experiments. AVAILABILITY AND IMPLEMENTATION: shinyDeepGxP can be launched on https://shiny.crc.pitt.edu/deepgxp/.

Journal Article↗

PubMed2026

REACTOR: REgulon Activity analysis and Comparison Tool for single-cell transcriptOmics Research.

SUMMARY: We introduce REACTOR, a computational tool designed to detect differential activity of transcriptional regulators and their target genes (regulons) in single-cell RNA-sequencing data. It expands the currently available framework for regulon analysis by introducing a robust statistical test to detect differential regulon activity between conditions, such as disease versus control, with multiple replicates. By contrasting different conditions, REACTOR enables identification of key condition- and cell type-specific regulons. To demonstrate the use of REACTOR, we illustrate its performance in a publicly available COVID-19 dataset. AVAILABILITY: REACTOR R-package together with an implementation vignette are available at https://www.github.com/elolab/REACTOR.

Regulon↗

PubMed2026

Elevated intron retention implicates neuroinflammation in brains of individuals with alcohol use disorder.

Intron retention, a form of alternative RNA splicing, can occur as part of normal gene regulation or result from disruption of the splicing machinery. Retained introns can potentially form double-stranded RNA, activating innate immune sensors and inflammation. This mechanism has been implicated in cancer but has not been studied in neuropsychiatric diseases like alcohol use disorder. We systematically analysed transcriptome-wide intron retention events in post-mortem brain tissue from 142 individuals (66 with alcohol use disorder and 76 controls), encompassing 320 region-specific samples from the superior frontal cortex, nucleus accumbens, central nucleus and basolateral amygdala. Analyses were adjusted for demographic, technical and biological covariates. Validation was performed in alcohol-preferring (P) rats using long-read sequencing. In complementary experiments, immunofluorescent staining was used to detect double-stranded RNA in rat brain tissue, while single-cell RNA-sequencing was performed to test activation of double-stranded RNA-sensing pathways in human brains. Brains from individuals with alcohol use disorder showed significantly higher total intron retention compared with controls, independent of age, with females showing greater increases than males. A total of 368 introns were positively associated with alcohol use disorder, and these introns were significantly longer and had weaker splice acceptor sites compared with non-associated introns. Genes harbouring these intron retention events were enriched in Purkinje neurons, visual cortex neurons and oligodendrocytes. Computational predictions indicated these long introns could form duplex RNA structures. Increased double-stranded RNA was confirmed experimentally in multiple brain regions of alcohol-consuming rats, where it co-localized primarily with neuronal nuclei and dendrites. In individuals with alcohol use disorder, we found that multiple pathways including double-stranded RNA responses, neuroinflammation, interferon and NF-κB signalling, adaptive immunity and apoptosis were activated. In addition, NeuN-positive neuronal counts significantly decreased in both the prefrontal and visual cortices. Furthermore, single-cell analysis demonstrated upregulation of TICAM1, the target of double-stranded RNA sensor TLR3, in oligodendrocytes, as well as widespread activation of downstream inflammatory pathways across glial and neuronal cell types. These findings provide the first evidence that chronic alcohol consumption promotes an overall increase of intron retention in the brain and is associated with the presence of double-stranded RNA. Furthermore, the double-stranded RNA may contribute to neuronal loss and brain pathology by activating a neuroinflammatory response.

alcohol use disorder↗

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

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

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

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

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↗