SEARCH · PubMed Health

Results for “Single-cell mRNA”

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

PubMed2025

Droplet-Based Single-Cell 3' mRNA Sequencing of Marburg Virus-Infected Samples.

Single-cell technologies are continually evolving with emerging methods that are gradually uncovering the central DNA-RNA-protein dogma. Single-cell RNA sequencing is one arm of a multi-omic approach that achieves an astounding level of granularity to reveal the complexity of virus-host interactions at the transcriptomic level. Cell tropism, virus replication, pathogenesis, and gene expression changes mediated by the virus and the host's immune response to infection are just some areas of study that are gaining better clarity due to the high-resolution analysis afforded by the technology.We describe a single-cell sequencing protocol for Marburg virus infection in vivo using nonhuman primate blood and the 10× Chromium Next GEM single-cell genomics methodology. Working with pathogens of high consequence is logistically complicated, requiring containment in biosafety level (BSL)-4 laboratories and harsh inactivation procedures before samples can safely be removed to lower biosafety conditions. We provide procedural insight into sample isolation and processing conducted in BSL-4 and describe the requirements for safe sample removal without jeopardizing quality for down-stream sequencing and analysis in BSL-2 conditions. Characterization of complicated biological processes mediated by high-containment pathogens, typically restricted to analogous model systems, e.g., minigenome, can be achieved using live virus.

Animals↗

PubMed2026

Multimodal Analysis Reveals Aberrant Expression of SUMO2 and Its Significant Association With Key Mechanisms of Metabolic Pathways in Hepatocellular Carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related deaths worldwide. However, the role of small ubiquitin-like modifier 2 (SUMO2), a core member of the small ubiquitin-like modifier (SUMO) family, regarding its expression patterns and metabolism-related functions in HCC remains inadequately understood. METHODS: A multidimensional analytical framework was applied, integrating immunohistochemistry (153 HCC vs. 21 non-HCC samples), proteomics (159 paired samples), bulk transcriptomics (3240 HCC vs. 2267 non-HCC samples), single-cell RNA sequencing (RNA-seq) (10 HCC vs. 8 non-HCC samples), spatial transcriptomics, and external CRISPR/Cas9 functional genomics data. Systematic analyses included standardized mean difference (SMD), pathway enrichment, pseudotime trajectory inference, in silico knockout, cell-cell communication, metabolic flux scoring, immune infiltration, clinical correlation, drug sensitivity prediction, and molecular docking. RESULTS: At the protein level, immunohistochemistry (nuclear positivity) and external proteomic data collectively demonstrated consistent SUMO2 overexpression in HCC. Consistent upregulation was also observed at the mRNA level across large-scale cohorts. Single-cell RNA-seq and spatial transcriptomics localized SUMO2 enrichment to malignant hepatocytes and tumor-dominant regions. CRISPR-mediated SUMO2 knockout suppressed proliferation in multiple HCC cell lines. Mechanistically, high SUMO2 expression was significantly associated with metabolic reprogramming involving glycolysis/gluconeogenesis, pyruvate metabolism, and the tricarboxylic acid cycle. SUMO2-high malignant hepatocyte subpopulations exhibited enhanced activity of the macrophage migration inhibitory factor signaling axis and enhanced iron-sensor interactions. Further, the immune infiltration analysis revealed a negative correlation between SUMO2 expression and M1 macrophages and a positive correlation between follicular helper T cells and regulatory T cells. Clinically, elevated SUMO2 levels were found to be associated with adverse prognostic features. Furthermore, high SUMO2 expression was associated with increased sensitivity to dasatinib, and molecular docking simulations predicted potential binding between SUMO2 and dasatinib, with a Vina score of -8.5 kcal/mol. CONCLUSIONS: SUMO2 is aberrantly expressed at the protein, mRNA, single-cell, and spatial transcriptomic levels in HCC and is significantly associated with metabolic reprogramming and altered migration inhibitory factor (MIF)-mediated intercellular communication, suggesting its potential as a novel biomarker for diagnosis and treatment.

Humans↗

PubMed2026

Dose-dependent IFN programs in myeloid cells after mRNA and adenovirus COVID-19 vaccination.

BACKGROUNDThe SARS-CoV-2 pandemic provided a rare opportunity to study how human immune responses develop to a novel viral antigen delivered through different vaccine platforms. However, to date, no study has directly compared immune responses to all 3 FDA-approved COVID-19 vaccines at single-cell multiomic resolution.METHODSWe longitudinally profiled SARS-CoV-2-naive adults (n = 31) vaccinated with BNT162b2, mRNA-1273, or Ad26.COV2.S, integrating plasma cytokines, antibody titers, and single-cell multiomic data (DOGMA-Seq).RESULTSWe discovered a distinct, transient IFN program termed ISG-dim, which emerged specifically 1-2 days after the first mRNA dose in approximately 10% of myeloid cells. This state was characterized by ISGF3 complex activation and its target genes (e.g., MX1, MX2, DDX58), with transcriptional and epigenetic profiles distinct from the robust IFN program observed after mRNA boosting or a single Ad26.COV2.S dose (ISG-high). In vitro stimulation of human monocytes showed that IFN-α alone recapitulates ISG-dim, whereas both IFN-α and IFN-γ are required for ISG-high.CONCLUSIONThese findings define dose-dependent IFN programming in human myeloid cells and highlight mechanistic differences between priming and boosting, with implications for optimizing vaccine platform choice, dose scheduling, and formulation.FUNDINGNIH grants AI142086, U19 AI135972, U01 AI165452, U01 AI165452, R01 AI160706, and P30 AG067988.

Humans↗

PubMed2026

MX1 promotes gastric cancer cell migration via inhibiting ANXA2 ubiquitination and degradation.

Gastric cancer (GC) is a globally lethal malignancy, with invasion and metastasis driving treatment failure and poor prognosis. MX dynamin like GTPase 1 (MX1) shows tumor-specific functional heterogeneity, while its expression, biological functions and molecular mechanisms in GC remain unclear. Here, we explored MX1's clinical significance and its regulatory mechanism in GC cell migration. We integrated public databases and institutional paired clinical samples for bioinformatics analysis of MX1's correlation with clinical outcomes, and verified its pro-migratory effect via Transwell and wound healing assays. Co-immunoprecipitation/mass spectrometry (Co-IP/MS), immunofluorescence and ubiquitination assays were used to identify MX1-interacting proteins and dissect the underlying mechanism, and the Genomics of Drug Sensitivity in Cancer database was applied for chemosensitivity analysis. MX1 was aberrantly upregulated in GC tissues and served as an independent prognostic biomarker, with high expression associated with shortened overall, first-progression and post-progression survival. MX1 promoted GC cell migration and epithelial-mesenchymal transition pathway enrichment, and directly bound Annexin A2 (ANXA2) in the cytoplasm; both were co-enriched in endothelial and epithelial cells by single-cell sequencing. MX1 dose-dependently upregulated ANXA2 protein (without affecting its mRNA) by inhibiting NEDD4L/TRIM65-mediated ANXA2 ubiquitination and degradation, enhancing ANXA2 stability. Additionally, high MX1 expression correlated with increased paclitaxel sensitivity in GC patients based on database analysis, and CCK-8 assays confirmed that MX1 overexpression significantly reduced the paclitaxel IC50 in gastric cancer cells, supporting its potential as a predictive biomarker for paclitaxel efficacy. This study demonstrates that MX1 promotes GC cell migration by suppressing ANXA2 ubiquitination and degradation, highlighting the critical role of the MX1-ANXA2 axis in GC progression. These findings provide novel molecular targets and theoretical support for GC prognostic evaluation, individualized chemotherapy and targeted therapy.

ANXA2↗

PubMed2026

Transcriptome changes in circulating immune cells of critical COVID-19 patients predict a specific metabolic and epigenetic imprint.

BACKGROUND: The progression to critical COVID-19 arises predominantly from a dysregulated host immune response although the underlying regulatory mechanisms still remain partially elusive. This limits a prompt prediction of the disease progression, reduces the therapeutic options and restrains our understanding of “long COVID”. METHODS: Here, we analyzed the transcriptome of peripheral blood mononuclear cells (PBMCs) collected from COVID-19 patients experiencing different degrees of the disease (mild and critical), and control patients enrolled in the clinical trial COntAGIouS as well as independent bulk RNA-seq, single-cell RNA-seq and proteomic datasets. RESULTS: In critical COVID-19 patients, the integrative analysis of transcriptomic data revealed an altered regulatory network involving microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and coding genes that control mRNA translation-related genes, epigenetics, and metabolism. In parallel, we observed an upregulation of tRNA aminoacylation genes in critical COVID-19 patients by the analysis of either bulk or single-cell RNA-seq data from publicly available independent cohorts. Additionally, we found increased expression of coding genes enriched for the cognate amino acids (glycine, alanine, isoleucine and tyrosine), all related to protein localization, post-translational modifications, and cell metabolism in our cohort. Similar alterations in amino acid frequency were found in an independent proteomic dataset. CONCLUSIONS: Collectively, our findings indicate a broad perturbation of the gene expression landscape that characterizes the aberrant host immune response in critical COVID-19 patients and is potentially coordinated by miRNA and tRNA metabolism alterations. TRIAL REGISTRATION: COntAGIouS, NCT04327570. Registered 26 March 2020, https://clinicaltrials.gov/ct2/show/NCT04327570 .

Female↗

PubMed2025

An advanced cytosine base editor enabled the generation of cattle with a stop codon in the β-lactoglobulin gene.

β-Lactoglobulin (BLG) is an allergen present in milk that can induce an acute immune response in certain individuals. The successful use of cytosine base editors (CBEs) can introduce stop codons into premature mRNA, thereby generating animals with disrupted genes that negatively regulate target traits. In this study, we employed a CBE system to target the major milk allergen BLG in bovine embryos, mammary epithelial cells, and live cattle. First, the precise single-base editing of the BLG gene in bovine embryos was achieved by designing an effective sgRNA to induce a c.61C > T substitution in the coding region, converting codon 21Gln (p.21Gln) to a premature stop codon. Sanger sequencing revealed an editing efficiency of 83.3% (20 out of 24 embryos), including two homozygous edits. Second, a bovine mammary epithelial cell line harboring BLG edits was constructed using the same CBE system. Sequencing showed that the designed sgRNA1 enabled the simultaneous conversion of three consecutive cytosines (c.59-61CCC > TTT) to thymines. At position c.61, single-cell clones exhibited monoallelic or biallelic editing (BLGc.61C > T), with monoallelic edits at positions c.59 and c.60 (CC > TT). Gene expression analysis confirmed that the BLGc.61C > T mutation effectively suppressed BLG expression at both the mRNA and protein levels, even in monoallelically edited cells. Finally, we successfully generated a heterozygous BLGc.61C > T single-base-edited dairy cow that despite its heterozygosity, showed significantly reduced BLG expression in the mammary epithelial cells and milk. Collectively, this study demonstrates the feasibility of using CBEs to disrupt BLG expression in dairy cows and provides a foundation for application in generating hypoallergenic dairy products.

Animals↗

PubMed2025

Insights into KIF11 pathogenesis in microcephaly-lymphedema-chorioretinopathy syndrome from a lymphatic perspective.

Pathogenic variants in kinesin KIF11 underlie microcephaly-lymphedema-chorioretinopathy (MLC) syndrome. Although well known for regulating spindle dynamics ensuring successful cell division, the association of KIF11 (encoding EG5) with development of the lymphatic system and how KIF11 pathogenic variants lead to lymphatic dysfunction and lymphedema remain unknown. Using patient-derived lymphoblastoid cells, we demonstrated that patients with MLC carrying pathogenic stop-gain variants in KIF11 have reduced mRNA and protein levels. Lymphoscintigraphy showed reduced tracer absorption, and intestinal lymphangiectasia was detected in one patient, pointing to impairment of lymphatic function caused by KIF11 haploinsufficiency. We revealed that KIF11 is expressed in early human and mouse development with the lymphatic markers VEGFR3, podoplanin, and PROX1. In zebrafish, single-cell RNA-Seq identified kif11 specifically expressed in endothelial precursors. In human lymphatic endothelial cells, EG5 inhibition with ispinesib reduced VEGFC-driven AKT phosphorylation, migration, and spheroid sprouting. KIF11 knockdown reduced PROX1 and VEGFR3 expression, providing for the first time to our knowledge a link between KIF11 and drivers of lymphangiogenesis and lymphatic identity.

Humans↗

PubMed2026

TET2 promotes monocyte inflammatory activation in asthma via ALKBH5-m6A regulation and PI3K signaling: evidence from m6A-SNP and single-cell analyses.

Asthma is a complex inflammatory airway disease with strong genetic determinants, yet the functional relevance of most asthma-associated non-coding variants remains unclear. Emerging evidence suggests that N6-methyladenosine (m6A) modification may serve as a critical epitranscriptomic link between genetic variation and immune regulation. In this study, we aimed to systematically identify functionally relevant m6A-regulated genes in asthma by integrating large-scale GWAS data, m6A-SNP annotations, and single-cell transcriptomic analyses, and to investigate their roles in monocyte-driven airway inflammation. We identified TET2 as a key m6A-regulated gene associated with both asthma and lung function, which was selectively upregulated in monocytes during asthma and accompanied by activation of inflammatory and PI3K signaling pathways. Mechanistic experiments further demonstrated that inflammatory stimulation induced ALKBH5 expression, reduced m6A modification of TET2 mRNA, and increased TET2 protein levels, thereby promoting PI3K/AKT signaling and pro-inflammatory cytokine production, whereas inhibition of TET2 or ALKBH5 attenuated these effects. Collectively, these findings demonstrate that ALKBH5-mediated m6A regulation of TET2 enhances PI3K/AKT signaling in monocytes, thereby promoting inflammatory responses in asthma. Our study establishes TET2 as a key m6A-regulated gene linking genetic susceptibility to monocyte-driven inflammation, and highlights the ALKBH5-m6A-TET2 axis as a potential therapeutic target for modulating aberrant immune responses in asthma.

Humans↗

PubMed2026

Increased IL4I1 expression predicts poor survival and modulates the immune microenvironment in acute myeloid leukemia.

BACKGROUND: The immunometabolic enzyme Interleukin-4-induced-1 (IL4I1) is implicated in cancer pathogenesis, yet its specific function and clinical relevance in acute myeloid leukemia (AML) remain unclear. METHODS: Comparative analysis of IL4I1 mRNA levels between AML patients and normal controls was performed using the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. The Kaplan–Meier survival analysis was conducted to evaluate the prognostic value of IL4I1. Functional insights were derived from analyses of differentially expressed genes (DEGs), Gene Set Enrichment Analysis (GSEA), and Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Immune infiltration was evaluated using the ssGSEA, ESTIMATE, quanTIseq and single-cell RNA sequencing (scRNA-seq) analysis. Finally, in vitro and in vivo functional experiments were perfromed to explore the impact of IL4I1 on AML progression and immunoregulation. RESULTS: IL4I1 expression was significantly elevated in AML compared to normal controls (p = 0.0004) and associated with poorer overall survival (p = 0.003). Bioinformatic analysis revealed that IL4I1 was linked to immune-related pathways—including humoral immune response, leukocyte interactions, and chemokine signaling—and to cellular amino acid metabolism. Its expression correlated with immune cell infiltration and checkpoint molecule expression. Experimentally, IL4I1 promoted leukemia cell proliferation in vitro and in vivo (p < 0.05). Furthermore, silencing IL4I1 suppressed M2 macrophage polarization and reduced secretion of inflammatory factors (p < 0.05). CONCLUSIONS: IL4I1 may serve as a potential biomarker for poor prognosis and an attractive target for immune-based therapeutic interventions in AML.

Humans↗

PubMed2025

Integrative multi-omics analysis proposes a metabolic classification of gliomas: distinct metabolic states, immune infiltration, and prognosis.

BACKGROUND: The tumor microenvironment (TME) of glioma harbors diverse cell types; however, cell metabolic heterogeneity remains to be explored. This study aims to characterize the metabolic features of different cell types in the TME by integrating multiple datasets, including genomics, bulk and single-cell transcriptomics, and metabolomics. METHODS: Unsupervised machine learning was used to construct an energy metabolic classifier based on the metabolic pathways identified from bulk RNA-seq of gliomas in the TCGA dataset. The classifier was externally validated using multiple datasets, including genomics, bulk RNA-seq, snRNA-seq, and the metabolomics data. Furthermore, metabolic heterogeneity associated with the classifier was further characterized at single-cell resolution. RESULTS: The energy metabolism-based classifier stratified patients into two prognostic clusters: patients in cluster 1 were characterized by high pathway activity of glycolysis, the pentose phosphate pathway (PPP), and fatty acid oxidation (FAO), whereas patients in cluster 2 exhibited higher activity in glutaminolysis. This metabolic classifier revealed both intratumoral and intertumoral metabolic heterogeneity, and the complexity was further validated by the metabolomics profiling and snRNA-seq data from the CPTAC dataset. Notably, OSMR, highly expressed in cluster 1, showed significant co-expression with key glycolytic enzyme genes. The OSM/OSMR/JAK1/STAT3 axis potently drives malignant progression of glioma cells, specially enhancing their invasive and migratory capabilities. Single-cell resolution analyses demonstrated that tumor metabolic heterogeneity is primarily driven by malignant cells rather than non-malignant components, while tumor microenvironment (TME) factors were also found to modulate malignant cell metabolism. Significantly, glycolytic activity in glioma cells increased during the phenotypic transition from PN (proneural) to MES (mesenchymal), with cluster 1 metabolic phenotypes predominating in the tumor core. Compared to cluster 2, cluster 1 patients exhibited higher mRNA expression of immunosuppressive checkpoint genes, which correlated with pronounced immunosuppression in the TME. Furthermore, various immune cells demonstrated distinct metabolic preferences at single-cell resolution. CONCLUSIONS: This study developed an energy metabolic-based classifier for gliomas with prognostic and therapeutic potential. Metabolic reprogramming was linked with the PN-to-MES transition of glioma cells and immunosuppression in the tumor microenvironment. Multi-omics data, especially snRNA-seq, offered insights into metabolism heterogeneity at single-cell resolution, enabling personalized treatment strategies.

Humans↗

PubMed2026

Exploration and experimental verification of triaptosis-related prognostic genes and cells in gastric cancer.

BACKGROUND: Triaptosis is a recently characterized form of programmed cell death with unclear implications in cancer. This study aimed to investigate the prognostic significance and biological relevance of triaptosis in gastric cancer (GC). METHODS: Transcriptomic and clinical data from TCGA-STAD and GSE62254, and single-cell RNA sequencing data from GSE183904 were analyzed. Triaptosis-related gene (TRG) scores were calculated using single-sample gene set enrichment analysis. Differentially expressed genes identified in TRG-score and GC-versus-normal comparisons underwent functional enrichment, Cox regression, and least absolute shrinkage and selection operator regression to develop an externally validated signature. Immune profiles, pathway activity, somatic mutations, tumor mutational burden (TMB), predicted drug sensitivity, and clinical features were compared by risk group. Single-cell analyses assessed TRG activity, prognostic gene expression, cell-cell communication, and pseudotime. Reverse transcription-quantitative PCR and Western blotting assessed mRNA expression and protein levels, respectively. RESULTS: A TRG-based prognostic model comprising ASPN, GRB14, and VTN was developed and externally validated, effectively distinguishing patients into two distinct risk groups with notably different survival outcomes. mRNA expression of all three genes and their protein levels were significantly higher in SGC-7901 cells than in GES-1 cells. High-risk patients had higher stromal scores and distinct immune profiles; 15 immune cell types differed between groups. Single-cell analysis revealed fibroblasts and pericytes among high-TRG-active cell types. Prognostic genes were significantly overexpressed in fibroblasts, which also showed high TRG activity. Fibroblasts demonstrated enhanced communication with pericytes, whereas tumor-derived fibroblasts showed weaker communication with macrophages, indicating immune microenvironment remodeling. CONCLUSION: The three-gene prognostic signature predicted GC prognosis and was associated with distinct immune and genomic features, suggesting potential value for risk stratification and personalized treatment.

Humans↗

PubMed2025

Founder Homozygous Nonsense CREB3 Variant and Variable-Onset Retinal Degeneration.

IMPORTANCE: Uncovering the genetic basis of inherited retinal diseases (IRDs) can enhance both diagnostic accuracy and the development of targeted treatment strategies. OBJECTIVE: To evaluate the association between a homozygous nonsense variant in CREB3 with IRDs. DESIGN, SETTING, AND PARTICIPANTS: Thirteen patients with a clinical diagnosis of retinitis pigmentosa or cone-rod degeneration were analyzed by whole-genome sequencing (WGS) and whole-exome sequencing (WES). Clinically, patients presented with 2 main phenotypes, rod-cone and cone-rod dystrophies, demonstrating variable electrophysiological and fundoscopic findings. Expression analysis was performed on patient-derived skin fibroblasts using the reverse transcription-polymerase chain reaction and Western blot analysis, and by interrogating previously published retinal single-cell RNA sequence data. Immunohistochemistry staining was performed on wild-type mouse retinal sections using an anti-CREB3 antibody. Patients with variable phenotypes of IRDs were recruited from 3 medical centers in Israel and Italy. Ophthalmologists clinically diagnosed patients at the relevant medical centers and referred them for genetic screening. WES and WGS were performed at different national and international centers, and the findings of the previously unreported gene were shared between investigators. EXPOSURES: CREB3 and IRDs. MAIN OUTCOMES AND MEASURES: The main outcome was evidence supporting an association between CREB3 and IRD. Measures included WES, WGS, and immunohistochemistry staining. RESULTS: A founder homozygous nonsense variant in CREB3 (c.881G>A, p.Trp294*) was identified in 13 patients from 4 unrelated families; 12 descendent from North-African Jewish origins and 1 from Italian origins. All patients manifested retinal degeneration with varying ages at onset. In patient-derived fibroblasts, the variant mRNA transcript generated a truncated CREB3 protein. Expression analysis and immunohistochemistry staining revealed CREB3 RNA and protein expression in various retinal cell types, indicating its vital role in photoreceptor function. CONCLUSIONS AND RELEVANCE: This study found an association between CREB3 and IRDs. CREB3 was previously shown to be upregulated following ultraviolet radiation. This might contribute to the extensive clinical variability observed in this relatively large cohort of homozygous patients with the same truncated variant.

Humans↗

PubMed2026

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α), 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+  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↗

PubMed2026

A single-nucleus transcriptome atlas of soybean anthers.

Anther development is crucial for plant sexual reproduction. However, a high-resolution, cell-type-specific transcriptomic atlas of this process is lacking for the legume crop soybean (Glycine max). Here, we construct a comprehensive transcriptional atlas of developing soybean anthers using single-nucleus RNA sequencing (snRNA-seq). We identify and characterize nine distinct cell types spanning both somatic and reproductive lineages. Our analysis reveals robust transcriptional continuity across anther developmental stages and dynamic reprogramming during key transitions. Notably, the shift from diploid meiocytes to haploid unicellular microspores is marked by the induction of previously inactive genes, despite an overall reduction in transcript abundance. Subsequently, within bicellular microspores, generative and vegetative cell lineages exhibit sharply divergent transcriptional programs: generative cells specialize in mRNA export and turnover, whereas vegetative cells up-regulate translational machinery. Evolutionary analysis further indicates that generative-cell-specific genes are subject to more relaxed purifying selection compared to those specific to vegetative cells. Functional validation using mutants generated by CRISPR/Cas9-mediated genome editing and EMS mutagenesis reveals the essential roles of OSD1A and PKSA in pollen development and fertility. This high-resolution atlas provides fundamental insights into the transcriptional regulation of soybean anther development and serves as a valuable resource for manipulating male fertility to advance hybrid breeding programs. The data are available at https://databases.genedenovo.com/pollen.

Glycine max↗

PubMed2025

Enhancing and accelerating cell type deconvolution of large-scale spatial transcriptomics slices with dual network model.

MOTIVATION: Cell type deconvolution deciphers spatial distribution of mRNA transcripts at single cell level by integrating single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics data to infer mixture of cell types of spots in slices. Current algorithms are criticized for neglecting connection between scRNA-seq and spatial transcriptomics data, as well as time-consuming, hampering their application to large-scale datasets. RESULTS: In this study, we propose a joint learning nonnegative matrix factorization algorithm for fast cell type deconvolution (aka jMF2D), which integrates scRNA-seq and spatial transcriptomics data with network models. To bridge scRNA-seq and spatial transcriptomics data, jMF2D jointly learns cell type similarity network to enhance quality of signatures of cell types, thereby promoting accuracy and efficiency of deconvolution. Experiments demonstrate that jMF2D outperforms state-of-the-art baselines in terms of accuracy by saving about 90% running time on various datasets generated by different platforms. Furthermore, it can also facilitates the identification of spatial domains and bio-marker genes, providing an efficient and effective model for analyzing spatial transcriptomics data. AVAILABILITY AND IMPLEMENTATION: The software is coded using python, and is free available for academic https://github.com/xkmaxidian/jMF2D.

Algorithms↗

PubMed2025

Targeting RECQL4 in hepatocellular carcinoma: from prognosis to therapeutic potential.

OBJECTIVE: The aim of this study is to assess the clinical utility of RecQ Like Helicase 4 (RECQL4) as a prognostic marker in hepatocellular carcinoma (HCC) and investigate its associations with various biological processes, angiogenesis-related factors, immune cell infiltration, immune checkpoints, and drug sensitivity. METHODS: RECQL4 expression was analyzed across a range of cancer types utilizing data from the TCGA database. Disparities in RECQL4 expression levels between normal and malignant tissues were evaluated, alongside an analysis of progression-free interval (PFI), disease-specific survival (DSS), and overall survival (OS) curves. Exploration of pertinent pathways, immune cell infiltration, single-cell RNA-seq data, and drug sensitivity was conducted employing The Cancer Genome Atlas (TCGA) and Tumor Immune Single-Cell Hub (TISCH) databases. Furthermore, validation of in-silico results was validated through qPCR, Western blotting, CCK-8 assay, EdU assay, clonogenic assay, wound-healing assay, and transwell assay. RESULTS: In HCC, RECQL4 was highly expressed and associated with poorer prognosis (p < 0.05). It positively correlated with pathways related to MYC targets, DNA replication, PI3K/AKT/mTOR signaling, DNA repair mechanisms, and the G2/M checkpoint (R > 0.24, p < 0.001). RECQL4 also showed significant correlations with angiogenesis-related genes, including PTK2 (R > 0.4, p < 0.05), suggesting a potential role in angiogenesis regulation. Immune analysis indicated that RECQL4 was associated with immune cell types such as T helper 2 cells, NK CD56bright cells, and follicular helper T cells, suggesting a positive relationship with their infiltration. High RECQL4 expression was also linked to increased sensitivity to drugs including Sorafenib, 5-Fluorouracil, Cisplatin, and Doxorubicin. Cellular experiments showed that RECQL4 expression at the mRNA and protein levels were significantly higher in HCC cell lines Hep3B and Huh7 compared to the normal liver cell line MHA. Moreover, RECQL4 knockdown resulted in reduced proliferation and migration in HCC cell lines (p < 0.05). CONCLUSIONS: RECQL4 shows promise as a biomarker for predicting recurrence and survival in HCC and may affect angiogenesis regulation. Its expression also appears to impact sensitivity to drugs such as Sorafenib, 5-Fluorouracil, Cisplatin, and Doxorubicin. Furthermore, silencing RECQL4 significantly inhibits HCC cell line proliferation and migration.

Humans↗

PubMed2025

Putative function and prognostic molecular marker of mast cells in colorectal cancer.

BACKGROUND: The increased demand for markers for colorectal cancer (CRC) highlights the importance of investigating immune cells involved in CRC progression. This study aims to dissect the mast cells in CRC, characterize the role of mast cells in CRC development, coordinate molecular communication between mast cells and malignant cells, and construct and validate a prognostic classification model based on mast cell markers. METHODS: Single-cell transcriptome data of CRC patients were extracted from GSE146771 for cell classification and annotation. The malignant cells were identified by copykat and the communication between mast cells and malignant cells was analyzed by CellChat. Least absolute shrinkage and selection operator (LASSO) regression analysis and Cox regression analysis of mast cell markers were performed in the TCGA-COAD cohort to construct a prognostic classification model. qRT-PCR was performed to detect the mRNA expression of the molecules in the classification model in P815 and MC-9 cells. The co-culture experiment of MC38 and P815 cells were performed in 12-well transwell dish. Wound healing assay and Transwell assay were performed to detect cell migration and invasion. RESULTS: 10,186 high-quality cells in GSE146771 were annotated to 9 cell types. Six markers in mast cells (HDC, GATA2, ASAH1, BTBD19, TIMP1, FAM110A) were selected to construct a classification model. The high-risk score defined showed high infiltration of immunosuppressive cells, including endothelial cells, CAFs, Tregs and high angiogenesis and epithelial-mesenchymal transition (EMT) activities. In the model, HDC were abnormally low expressed in P815 cells, while BTBD19, FAM110A, GATA2, ASAH1 and TIMP1 showed excessive expression in P815 cells. Knockdown of GATA2 in the co-culture system of P815 and MC38 cells blocked cell migration and invasion. CONCLUSION: This study identified the cell types within CRC, elaborated the cellular functions of mast cells in CRC development and their molecular communication to coordinate malignant cells, and highlighted the molecular components and biological features that constitute promising prognostic classification model.

Mast Cells↗

PubMed2026

NPLOC4 Constructs Tumor Immunosuppressive Microenvironment in Pan-cancer and Hepatocellular Carcinoma.

INTRODUCTION: NPLOC4 (nuclear protein localization 4 homolog) is mainly involved in DNA damage, cell cycle, and ubiquitination promotion. Nonetheless, the role of NPLOC4 in the tumor immune microenvironment (TIME) and its potential as a promising tumor therapeutic target remains unclear. METHODS: Therefore, analyses of NPLOC4 mRNA and protein expression, RNA subcellular localization, and patient prognosis associated with NPLOC4 expression were conducted across multiple tumor types. Additionally, the correlations between NPLOC4 and immune cells, non-immune cells, and immune molecules within the tumor immune microenvironment (TIME) were investigated. These analyses utilized data from various public resources, including the Genotype-Tissue Expression (GTEx) project, The Cancer Genome Atlas (TCGA), Cancer Cell Line Encyclopedia (CCLE), The Human Protein Atlas (HPA), Clinical Proteomic Tumor Analysis Consortium (CPTAC), TIMER2.0, KM-Plotter, The University of Alabama at Birmingham Cancer Data Analysis Portal (UALCAN), and Tumor Immune Single-cell Hub 2 (TISCH2). Subsequently, we utilized hepatocellular carcinoma (HCC) patients' cancer and adjacent tissues plus tumor cell lines to verify the differential RNA and protein expression of NPLOC4 via qRT-PCR and immunohistochemistry (IHC). Then, the relationship of NPLOC4 expression level with immune infiltration score, infiltration of effector immune cells, suppressive immune cells, and several vital immune checkpoints was analyzed in HCC immune microenvironment. Furthermore, the distribution of expression of NPLOC4 in various cells in the HCC microenvironment was determined through single-cell sequencing analysis. RESULTS: We discovered that NPLOC4 was up-regulated in a variety of tumors and was correlated with poor prognosis. NPLOC4 not only had the potential as a tumor prognostic marker and therapeutic target but also was strongly linked to immune cells, immune checkpoints, and immune-related molecules and pathways in HCC immune microenvironment. CONCLUSION: In summary, NPLOC4 may serve as a promising target for immunotherapy.

Humans↗