Search PubMedSearch

SEARCH · Search PubMed

Results for “single-cell transcriptome sequencing”

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

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

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

Integration of single-cell transcriptomics and genomic mutation analysis identifies an immunotherapy-resistant tumor subcluster and validates ARNTL2 as a malignant driver in lung adenocarcinoma.

BACKGROUND: Immunotherapy resistance in lung adenocarcinoma (LUAD) remains a critical clinical challenge, and the mechanisms underlying resistance-associated intratumoral heterogeneity are poorly characterized. METHODS: We performed single-cell RNA sequencing of LUAD patients receiving neoadjuvant immunotherapy (responders vs. non-responders), integrating inferCNV, GSVA, and differential expression analyses. Cluster-specific genes were validated across seven independent cohorts (TCGA-LUAD, GSE13213, GSE26939, GSE29016, GSE30219, GSE31210, GSE42127). A multi-algorithm machine learning framework was used to construct a prognostic model, and the immune microenvironment was characterized using TCIA scoring, seven infiltration algorithms, and ESTIMATE. ARNTL2 function was assessed by CCK-8 and Transwell assays in A549 and H1299 cells. RESULTS: Non-responders showed significant enrichment of epithelial cells, depletion of cytotoxic T/NK cells, and elevated copy number variation burden versus responders (p < 0.0001). A resistance-enriched malignant subcluster (Cluster 2) exhibited hyperproliferative and metabolic reprogramming signatures with upregulated KRT17, S100A2, and CST6, which showed tumor-specific overexpression, adverse prognostic value, and genomic amplification across cohorts. CoxBoost combined with survivalSVM achieved optimal predictive performance (C-index = 0.686), yielding robust risk stratification (HR: 2.54-10.51, all p < 0.05). Low-risk patients showed greater immune infiltration and higher TCIA immunophenoscores. ARNTL2 was an independent prognostic factor (HR: 2.07-4.64) strongly correlated with risk score (r = 0.69), and its knockdown suppressed proliferation and invasion in both LUAD cell lines (all p < 0.05). CONCLUSION: This study identifies a resistance-associated malignant subcluster in LUAD, constructs a validated CoxBoost + survivalSVM prognostic model with robust immune stratification, and establishes ARNTL2 as a core oncogenic driver and therapeutic target.

ARNTL2

Integrated single-cell and bulk transcriptomic analysis identifies a novel senescent fibroblast subtype associated with poor prognosis in acral melanoma.

BACKGROUND: Acral melanoma (AM) exhibits significant intratumoral heterogeneity, but its tumor microenvironment (TME) and immune regulation remain unclear. This study aims to dissect TME heterogeneity and establish a prognostic model based on key cell subpopulations. METHODS: We collected AM single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). Unsupervised clustering, CellChat, and Scissor analysis were performed to characterize cellular heterogeneity, cell-cell communication, and prognosis-related cell subpopulations. Kaplan-Meier analysis was used to assess the prognostic value of key genes, which were further validated by multiplex immunohistochemistry (mIHC). RESULTS: In AM, Mel_C2, C7, and C9 with high SEMA6A and KIT expression were strongly linked to poor prognosis. We further identified a senescent fibroblast subpopulation (sCAF_CDKN2A) characterized by high fibroblast senescence signature (FSS) scores. Integrating Scissor analysis of fibroblast subtypes with bulk prognostic data, we identified COL3A1, VCAN, and KIT as prognosis-associated genes upregulated in poor-outcome-related fibroblast subsets. Cell-cell communication analysis revealed that sCAF_CDKN2A engages in an immunosuppressive network, interacting with regulatory T cells (Tregs) via MIF signaling and receiving signals from exhausted CD8+ T cells through PPIA-BSG interactions. Using transcription factor expression patterns from these fibroblast subtypes, we constructed a prognostic model that effectively stratified patients into distinct risk groups with significant differences in overall survival (OS). mIHC confirmed significantly higher protein levels of SEMA6A and COL3A1 in tumor tissues compared to matched normal tissues. CONCLUSIONS: We established a novel prognostic model for AM and identified sCAF_CDKN2A as an immunosuppressive senescent fibroblast subpopulation driving poor prognosis.

Acral melanoma

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

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

Multilevel genomic, transcriptomic, and epidemiologic evidence linking diabetic retinopathy to Alzheimer disease.

BACKGROUND: Diabetic retinopathy (DR) and Alzheimer disease (AD) share metabolic and vascular dysfunctions, but the extent to which they reflect overlapping genetic susceptibility and neurovascular-metabolic regulatory pathways remains unclear. We combined multi-omics analyses with population-based data to examine the genetic convergence, cellular pathways, and longitudinal association between DR and AD. METHODS: We performed a two-sample Mendelian randomisation (MR) to estimate the association between genetically predicted DR liability and AD risk. We used Bayesian colocalisation analysis to identify shared genomic loci, and summary-data-based MR (SMR) to detect expression-mediated genes jointly associated with DR and AD. We analysed single-cell RNA sequencing data to characterise shared cellular features and related biological pathways. We also conducted an MR-based mediation analysis to explore whether lipid-related, metabolic, or inflammatory traits mediated the observed DR-AD association, and a longitudinal analysis of the UK Biobank cohort to assess the association between DR and incident AD. RESULTS: With the MR analysis, we found that genetically predicted liability to DR was associated with a modest increase in AD risk. Colocalisation analysis supported a shared genetic signal. We identified three genes with shared expression-mediated associations across DR and AD through SMR. Functional enrichment analyses revealed partially overlapping neurovascular and metabolic pathways. Using MR-based mediation analysis, we found no significant intermediary traits linking DR and AD. Findings from the UK Biobank cohort were directionally consistent with the genetic analyses. CONCLUSIONS: Genetic liability to DR is associated with an increased risk of AD and is accompanied by shared expression-mediated effects and convergent neurovascular-metabolic pathways. These findings support the possibility that DR may serve as a clinically accessible indicator of increased neurodegenerative vulnerability.

Humans

Proteomics combined with single-cell sequencing reveals key genes and computational lead compound related to ligamentum flavum hypertrophy, lactate metabolism and lactate modification.

Ligamentum flavum hypertrophy (LFH) is a hallmark pathological feature of lumbar spinal stenosis; however, its underlying molecular mechanisms remain incompletely understood. Lactate metabolism and related lactylation modifications have emerged as critical links between cellular metabolism and epigenetic regulation, with established roles in various fibrotic and inflammatory diseases. Nevertheless, the specific contribution of lactylation to LFH pathogenesis remains unexplored. In this study, we integrated proteomic profiling of ligamentum flavum tissues with single-cell transcriptomic data to identify differentially expressed proteins associated with LFH. Cross-referencing these genes with genes involved in lactate metabolism and lactylation yielded 16 candidate genes. Through functional enrichment analysis, protein-protein interaction network construction, and GraphBAN model prediction, we identified five hub genes (NDUFS2, HMOX1, SPR, FABP5, and PFKP) and two potential lead compounds (ZINC000014879975 and ZINC000242437513). Molecular docking analysis confirmed favorable binding affinities between these compounds, suggesting that they may serve as potential lead compounds worthy of further experimental investigation. Single-cell analysis further revealed that macrophages occupy a central position in the LFH microenvironment, resulting in pronounced metabolic reprogramming and remodeling of intercellular communication networks, particularly via the MIF-CD74/CD44 axis, under pathological conditions.

Proteomics

Translating single-cell RNA sequencing into monocyte direct leukocyte subpopulation-transcript abundance assay ratio-based biomarkers (IFI27/PSAP or IFI27/CTSS) for clinical detection of viral infection.

A rapid method for triaging febrile patients by aetiology (e.g., viral or bacterial infection) using gene expression in peripheral blood (PB) is an intensively researched area. However, gene expression in blood represents a composite sum of gene expression of all the component cell types present in the sample. As a result, numerous genes are measured in most proposed signatures. Herein, we propose a simple ratio-based biomarker (RBB) called direct leukocyte subpopulation-transcript abundance assay (DIRECT LS-TA) that recapitulates gene expressions of a single cell type in PB (i.e., monocytes). Based on single-cell RNA sequencing (scRNAseq) data and bulk expression data, IFI27 and SIGLEC1 are found as interferon-stimulated genes (ISGs) predominantly expressed by monocytes. The DIRECT LS-TA method can use a simple ratio of two genes measured in PB as an RBB to represent the target gene expression in monocytes without the need for monocyte purification. Both scRNAseq and bulk RNA sequencing datasets were used to evaluate the correlation between ISG expression in monocytes and PB, with a particular focus on monocyte expression of IFI27. An iceberg plot of bulk transcriptome data was used to identify genes that were predominantly expressed by monocytes in PB. DIRECT LS-TA RBBs of the three genes (IFI27, IFI44L and SIGLEC1) were evaluated by group-wise comparison, receiver operating characteristic and meta-analysis. In addition, the conventional interferon (IFN) score was evaluated for comparison of diagnostic performance. In viral infection datasets, DIRECT LS-TA of IFI27 (IFI27/PSAP or IFI27/CTSS) was most intensely activated (p value by t test <1e-9) and had the best area under the curve (0.94) among the three potential monocyte ISGs analysed. DIRECT LS-TA SIGLEC1 was also another monocyte biomarker but showed a lower activation (p<9e-5). IFI27/PSAP showed better diagnostic performance than the conventional IFN score. On the other hand, IFI44L was not a predominant monocyte expression gene. DIRECT LS-TA of IFI27 (IFI27/PSAP or IFI27/CTSS) measured in PB was the best biomarker of viral infection and IFN activation among ISGs predominantly expressed by monocytes. It performed even better than the conventional IFN score which required quantification of eight genes. The results suggest that DIRECT LS-TA of IFI27 is a monocyte-informative biomarker which is easy to determine in PB without the need for cell sorting.

Humans

Strategies for mosaic variant calling in brain disorders.

The human brain is a genomic mosaic, where postzygotic mutations arising from embryogenesis to senescence drive diverse neurodevelopmental and neurodegenerative diseases. Because of numerous sequencing artifacts at ultralow variant allele frequencies (VAFs), detecting these variants remains a significant analytical challenge. This review focuses on single-nucleotide variants and small indels, summarizing current strategies for aligning sampling methods, including bulk, laser capture microdissection, and single-cell genomics, with the expected clonal architecture of the brain. It emphasizes that mosaic detection sensitivity is fundamentally constrained by sequencing depth, since even the most advanced algorithms cannot identify variants not physically represented in the sequencing library. The review further recommends the selection of variant calling algorithms based on validated VAF detection performance, matching tools like MuTect2 and MosaicForecast to their optimal performance ranges. Furthermore, we discuss how multitissue sampling, as emphasized by the SMaHT project, addresses the matched-control dilemma and supports accurate variant classification via cross-tissue VAF gradients. Integrating these established pipelines with multiomics modalities, including transcriptomic and epigenetic data, could advance the field toward a functional understanding of how the somatic genome impacts human brain health and disease.

Humans

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

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

Animals

Single-cell sequencing reveals synovial fluid &#x3b3;&#x3b4; T-cell expansion in equine experimental osteoarthritis.

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

equine

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

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

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&#x2009;014 cells) from adult patients with burn injuries within postburn day 17 (n&#x2009;=&#x2009;4) and healthy donors (n&#x2009;=&#x2009;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&#xef;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&#x2009;=&#x2009;158; AUC 0.81) and intensive care unit COVID-19 (n&#x2009;=&#x2009;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

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

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

Atherosclerosis

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

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

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

Animals

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