Search PubMed⌕ Search

SEARCH · Search PubMed

Search Search PubMed

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.

11 recordsLinked to original sources
PubMed2026

Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.

Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature‑supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR = 0.52) and its potential regulation of risk factors IL2RA (OR = 0.46) and HLA-DR (OR = 0.40). Conversely, IL2RA (OR = 1.42), HLA-DR (OR = 1.88), and MIF (OR = 1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+ HLA-DR+ CD74+ monocytes and CD4+ IL2RA+ T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.

Humans↗
PubMed2026

Decoding tumor immune microenvironment heterogeneity by single-cell and spatial multi-omics: From immunotherapy resistance to translational biomarkers.

Immune checkpoint blockade has transformed cancer therapy, yet primary and acquired resistance remain major clinical challenges. Increasing evidence indicates that immunotherapy resistance cannot be fully explained by tumor-intrinsic alterations or conventional biomarkers such as PD-L1 expression, tumor mutational burden, or microsatellite instability. Instead, therapeutic response is shaped by the tumor immune microenvironment (TIME) as a heterogeneous, spatially organized, and dynamically evolving ecosystem. Single-cell omics has revealed diverse immune and stromal cell states, including progenitor and terminally exhausted T cells, suppressive myeloid programs, B-cell/TLS-associated immune-reactive states, and CAF-mediated exclusion phenotypes. Spatial transcriptomics, spatial proteomics, and imaging-based approaches further demonstrate that these cell states assemble into distinct immune niches, including immune-inflamed, T-cell-excluded, myeloid-suppressive, metabolic/hypoxic, and TLS-associated niches. These spatial ecosystems determine whether antitumor immune cells can access malignant cells, receive antigen-presenting support, or become restrained by stromal, vascular, metabolic, and myeloid barriers. In this review, we summarize how single-cell and spatial multi-omics redefine TIME heterogeneity in immunotherapy resistance, highlight ligand-receptor communication networks linking cell states to spatial immune dysfunction, and discuss emerging translational biomarkers for patient stratification. We further propose that future immunotherapy biomarkers should evolve from static single-marker assays toward longitudinal, spatially resolved, and interpretable multi-omics models that guide precision combination immunotherapy.

Humans↗
PubMed2026

Single-cell RNA sequencing provides further insights into the immunostimulatory action of freeze-dried Lactiplantibacillus plantarum on Penaeus vannamei shrimp.

Immunostimulation through dietary interventions opened new avenues in developing disease control and prevention tools for shrimp aquaculture. We have previously shown that feeding with freeze-dried Lactiplantibacillus plantarum (LAB) increased disease resistance of Penaeus vannamei against both Vibrio parahaemolyticus and white spot syndrome virus (WSSV) based on bulk RNA sequencing of shrimp gills. This tissue participates in ion transport and serves as a first line of defense against environmental stressors and pathogenic infections. However, characterization of their cell composition and functions remains limited. Here, we implemented a single-cell RNA sequencing approach to further gather insights into how feeding with freeze-dried LAB modulates host immunity which may not be evident with bulk RNA sequencing approach. A total of five clusters with unique transcriptional signatures were identified, corresponding to pillar cells, septal cells, and sessile hemocytes. Pseudo-bulk analyses at global- and cluster-levels showed differential expression of genes related to host immunity and metabolism. We further revealed how overall transcriptomic changes are not exclusively caused by gene expression changes but may also be driven by cell population dynamics. This study highlighted how single-cell RNA sequencing approach may shed light on the mechanisms of action of immunostimulants which may be masked in bulk transcriptome analyses.

Animals↗
PubMed2026

Single-cell transcriptomics reveals heterogeneous stress responses and Mg2+-mediated survival mechanisms in Lactobacillus delbrueckii subsp. bulgaricus during freeze-drying and storage.

Maintaining the viability of lactic acid bacteria during dehydration and subsequent storage remains a significant challenge. Here, we employed single-cell RNA sequencing to reveal the heterogeneous stress responses of Lactobacillus delbrueckii subsp. bulgaricus, identifying seven distinct transcriptional clusters across the liquid culture, freeze-drying, and storage phases. The dominant clusters in the freeze-drying and storage were not completely consistent, showing significant functional differentiation. Genomic stability may be important for survival during freeze-drying and storage, while intracellular energy homeostasis appears important for viability during storage. The magnesium transporter mgtB was highly expressed in clusters tolerant to freeze-drying and storage, suggesting a critical role for Mg2+ homeostasis. Further experimental validation confirmed that Mg2+ treatment significantly bolstered stress resistance, increasing immediate post-freeze-drying survival by over 2-fold (up to 92.90%) and post-storage survival by over 5-fold (up to 5.98%). Proteomic data indicated that Mg2+ supplementation correlated with the maintenance of several biological functions potentially relevant to bacterial survival during freeze-drying and storage, including DNA repair, translation, and central carbon metabolism. These findings provide a map of microbial stress resistance through population heterogeneity and offer a potential strategy that may be adapted for enhancing the stability of other industrial lactic acid bacteria products.

Freeze Drying↗
PubMed2026

Endocrine-disrupting chemical-induced gene networks confer coronary heart disease risk revealed by causal inference and single-cell analyses.

BACKGROUND: Endocrine-disrupting chemicals (EDCs) are linked to coronary heart disease (CHD), but underlying mechanisms remain unclear. We aimed to identify EDC-related genes and evaluate their causal roles in CHD. METHODS: We curated EDC-related genes from a compound-gene interaction database and integrated them with CHD genome-wide association study (GWAS) summary statistics and tissue-specific expression quantitative trait loci (eQTL) data. Two-sample Mendelian randomization (MR) and Bayesian colocalization were applied to infer causality. Functional enrichment, single-cell RNA sequencing of human coronary arteries, and EDC-gene networks were further analyzed. RESULTS: After FDR correction, 39 genes were significantly associated with CHD risk via MR. Four genes-ZNF827, FCHO1, IPO9 (protective), and RPL13 (risk-increasing)-showed strong colocalization (PPH4 > 0.9). Pathway and single-cell analyses of coronary artery tissue indicated that vascular and immune pathways mediate these effects. An interaction network highlighted associations between specific EDCs and candidate genes implicated in CHD susceptibility. CONCLUSION: This integrative genomic study provides evidence that EDCs influence CHD susceptibility through distinct gene networks, revealing potential mechanisms and molecular targets for prevention and therapy.

Humans↗
PubMed2026

Mapping the immune-genetic architecture of Epstein-Barr virus-related phenotypes and multiple sclerosis through a single-cell genetic framework for target prioritization and pharmacologic hypothesis generation.

BACKGROUND: Multiple sclerosis (MS) is a severe neuroinflammatory disease causing substantial long-term disability. Strong epidemiologic evidence links Epstein-Barr virus (EBV) exposure with MS risk, but genetic evidence for immune target prioritization in EBV-related phenotypes remains limited. METHODS: We integrated single-cell cis-eQTL data from 14 immune cell types with GWASs of an EBV-related clinical phenotype and MS using a single-cell Mendelian randomization framework with colocalization analyses. Candidate eGenes were evaluated in independent cohorts. For multi-SNP instruments, we performed heterogeneity, pleiotropy, MR-Egger, weighted median, mode-based, and MR-PRESSO sensitivity analyses. We also conducted phenome-wide association analyses and queried DrugBank to annotate candidate compounds targeting prioritized genes. RESULTS: We prioritized 43 immune-cell-specific candidate eGenes with convergent genetic support, including 6 for the EBV-related phenotype and 37 for MS. SERPINB1 in NK cells was associated with increased risk of the EBV-related phenotype, whereas HLA-G was associated with decreased risk. For MS, APOM and MSH5 showed protective associations, while AHI1 showed cell-type-dependent, bidirectional associations across immune lineages. Colocalization and independent cohort evaluation supported these findings. Among FDR-significant multi-SNP associations, MR-Egger intercept tests did not indicate directional pleiotropy, although a small subset showed heterogeneity or MR-PRESSO signals. Phenome-wide analyses identified no significant adverse phenotypic associations among evaluable genes at the prespecified threshold. DrugBank annotation nominated sodium nitroprusside, fasudil, artenimol, and choline as hypothesis-generating compounds for experimental follow-up. CONCLUSIONS: This study provides a single-cell genetic framework for prioritizing immune-cell-specific candidate targets for EBV-related phenotypes and MS, and nominates genetically supported targets and pharmacologic hypotheses for experimental investigation.

Humans↗
PubMed2026

Deciphering CD8+ T cell exhaustion in human cancers through single-cell and spatial transcriptomics.

Exhausted CD8+ T cells (Tex) within the tumor microenvironment (TME) represents a critical barrier limiting anti-tumor immune responses. Tex cells are characterized by upregulated inhibitory immune checkpoint receptors, reduced cytotoxicity, and functional heterogeneity. Their genomic features and regulatory networks remain poorly defined, and only a minority of patients respond to immune checkpoint blockade (ICB) therapy. Single-cell RNA sequencing (scRNA-seq), through high-resolution transcriptomic profiling, has revealed diverse Tex subpopulations, identified subpopulation-specific marker genes and regulatory pathways. Spatial transcriptomics has further mapped the spatial distribution of Tex and their interaction networks with immune cells, tumor cells, and stromal cells, elucidating the impact of spatial heterogeneity on Tex functionality. Current studies indicate that the exhausted state of Tex is dynamic and modifiable, with functional differences among subpopulations closely associated with tumor progression and therapeutic response. However, the genomic characteristics, epigenetic regulation, and spatial interaction mechanisms of Tex require further exploration. This review summarizes recent advances in high-resolution omics technologies for precisely dissecting Tex heterogeneity, functional features, and interactions with other cells. It emphasizes the central value of optimizing Tex-targeted tumor immunotherapy strategies, providing theoretical foundations and directional guidance for developing more effective anti-tumor immunotherapies.

Humans↗
PubMed2026

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

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

Adenosine Deaminase↗
PubMed2026

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↗
PubMed2026

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

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

Animals↗
PubMed2026

Single-cell 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↗