Objective To explore the prognostic value of glycolysis-related genes in colorectal cancer (CRC) patients and formulate a novel glycolysis-related prognostic risk model. Methods Single-cell and bulk transcriptomic data of CRC patients, along with clinical information, were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Glycolysis scores for each sample were calculated using single-sample Gene Set Enrichment Analysis (ssGSEA). Kaplan-Meier survival curves were generated to analyze the relationship between glycolysis scores and overall survival. Novel glycolysis-related subgroups were defined among the cell type with the highest glycolysis scores. Gene enrichment analysis, metabolic activity assessment, and univariate Cox regression were performed to explore the biological functions and prognostic impact of these subgroups. A prognostic risk model was built and validated based on genes significantly affecting the prognosis. Gene Set Enrichment Analysis (GSEA) was conducted to explore differences in biological processes between high- and low-risk groups. Differences in immune microenvironment and drug sensitivity between these groups were assessed using R packages. Potential targeted agents for prognostic risk genes were predicted using the Enrichr database. Results Tumor tissues showed significantly higher glycolysis scores than normal tissues, which was associated with a poor prognosis in CRC patients. The highest glycolysis score was observed in epithelial cells, within which we defined eight novel glycolysis-related cell subpopulations. Specifically, the P4HA1⁺ epithelial cell subpopulation was associated with a poor prognosis. Based on signature genes of this subpopulation, a six-gene prognostic risk model was formulated. GSEA revealed significant biological differences between high- and low-risk groups. Immune microenvironment analysis demonstrated that the high-risk group had increased infiltration of macrophages and tumor-associated fibroblasts, along with evident immune exclusion and suppression, while the low-risk group exhibited higher levels of B cell and T cell infiltration. Drug sensitivity analysis indicated that high-risk patients were more sensitive to Abiraterone, while low-risk patients responded to Cisplatin. Additionally, Valproic acid was predicted as a potential targeted agent. Conclusion High glycolytic activity is associated with a poor prognosis in CRC patients. The novel glycolysis-related prognostic risk model formulated in this study offers significant potential for enhancing the diagnosis and treatment of CRC.
Humans ; Colorectal Neoplasms/pathology* ; Glycolysis/genetics* ; Prognosis ; Transcriptome ; Tumor Microenvironment/genetics* ; Gene Expression Profiling ; Single-Cell Analysis ; Gene Expression Regulation, Neoplastic ; Male ; Female ; Kaplan-Meier Estimate

Objective To explore macrophage subpopulations in ischemic stroke (IS) by using single-cell RNA sequencing (scRNA-seq) data analysis and High-Dimensional Weighted Gene Co-Expression Network Analysis (hdWGCNA). Methods Based on single-cell sequencing data, transcriptomic information for different cell types was obtained, and macrophages were selected for subpopulation identification. hdWGCNA, cell-cell communication, and pseudotime trajectory analysis were used to explore the characteristics of macrophage subpopulations following IS. Key genes related to IS were identified using microarray data and validated for diagnostic potential through Receiver Operating Characteristic (ROC) analysis. Gene Set Enrichment Analysis (GSEA) was conducted to investigate the potential functions of these genes. Results The scRNA-seq data analysis revealed significant changes in macrophage subpopulation composition after IS. A specific macrophage subpopulation enriched in the stroke group was identified and designated as MCAO-specific macrophages (MSM). Pseudotime trajectory analysis indicated that MSM cells were in an intermediate stage of macrophage differentiation. Cell-cell communication analysis uncovered complex interactions between MSM cells and other cells, with the CCL6-CCR1 signaling axis potentially playing a crucial role in neuroinflammation. Two gene modules associated with MSM were identified via hdWGCNA, significantly enriched in pathways related to NOD-like receptors and antigen processing. By integrating differentially expressed MSM genes with conventional transcriptomic data, three IS-related hub genes were identified: Arg1, CLEC4D, and CLEC4E. Conclusion This study reveals the characteristics and functions of macrophage subpopulations following IS and identifies three hub genes with potential diagnostic value, providing novel insights into the pathological mechanisms of IS.
Macrophages/metabolism* ; Computational Biology/methods* ; Single-Cell Analysis/methods* ; Transcriptome ; Ischemic Stroke/metabolism* ; Animals ; Gene Regulatory Networks ; Gene Expression Profiling ; Humans ; Male

OBJECTIVES: To elucidate the underlying mechanisms of macrophage-mediated inflammation and tissue injury in diabetic foot ulcer (DFU). METHODS: Skin tissue samples were collected from patients with DFU and with non-DFU. A total of 79 272 high-quality cell transcriptomes were obtained using single-cell RNA sequencing. An unbiased clustering approach was employed to identify cell subpopulations. Seurat functions were used to identify differentially expressed genes between DFU and non-DFU groups, and gene ontology (GO) enrichment analysis was used to reveal gene function. Furthermore, cell-cell communication network construction and ligand-receptor interaction analysis were performed to reveal the mechanisms underlying cellular interactions and signaling regulation in the DFU microenvironment from multiple perspectives. RESULTS: The results revealed a significant expansion of myeloid cells in DFU tissues, alongside a marked reduction in structural cells such as endothelial cells, epithelial cells, and smooth muscle cells. Major cell types underwent functional reprogramming, characterized by immune activation and impaired tissue remodeling. Specifically, macrophages in DFU skin tissues exhibited a shift toward a pro-inflammatory M1 phenotype, with upregulation of genes associated with inflammation and oxidative stress. Cell communication analysis further demonstrated that M1 macrophages served as both primary signal receivers and influencers in the COMPLEMENT pathway mediated communication network, and as key signal senders and mediators in the secreted phosphoprotein 1 (SPP1) pathway mediated communication network, actively shaping the inflammatory microenvironment. Key ligand-receptor interactions driving macrophage signaling were identified, including C3-(ITGAM+ITGB2) and SPP1-CD44. CONCLUSIONS This study establishes a comprehensive single-cell atlas of DFU, revealing the role of macrophage-driven cellular networks in chronic inflammation and impaired healing. These findings may offer potential novel therapeutic targets for DFU treatment.
Humans ; Macrophages/immunology* ; Diabetic Foot/pathology* ; Single-Cell Analysis ; Transcriptome ; Gene Expression Profiling ; Inflammation ; Skin ; Cell Communication ; Signal Transduction ; Cellular Reprogramming

Recent data suggest that vascular endothelial growth factor receptor inhibitor (VEGFRi) can enhance the anti-tumor activity of the anti-programmed cell death-1 (anti-PD-1) antibody in colorectal cancer (CRC) with microsatellite stability (MSS). However, the comparison between this combination and standard third-line VEGFRi treatment is not performed, and reliable biomarkers are still lacking. We retrospectively enrolled MSS CRC patients receiving anti-PD-1 antibody plus VEGFRi (combination group, n=54) or VEGFRi alone (VEGFRi group, n=32), and their efficacy and safety were evaluated. We additionally examined the immune characteristics of the MSS CRC tumor microenvironment (TME) through single-cell and spatial transcriptomic data, and an MSS CRC immune cell-related signature (MCICRS) that can be used to predict the clinical outcomes of MSS CRC patients receiving immunotherapy was developed and validated in our in-house cohort. Compared with VEGFRi alone, the combination of anti-PD-1 antibody and VEGFRi exhibited a prolonged survival benefit (median progression-free survival: 4.4 vs. 2.0 months, P=0.0024; median overall survival: 10.2 vs. 5.2 months, P=0.0038) and a similar adverse event incidence. Through single-cell and spatial transcriptomic analysis, we determined ten MSS CRC-enriched immune cell types and their spatial distribution, including naive CD4⁺ T, regulatory CD4⁺ T, CD4⁺ Th17, exhausted CD8⁺ T, cytotoxic CD8⁺ T, proliferated CD8⁺ T, natural killer (NK) cells, plasma, and classical and intermediate monocytes. Based on a systemic meta-analysis and ten machine learning algorithms, we obtained MCICRS, an independent risk factor for the prognosis of MSS CRC patients. Further analyses demonstrated that the low-MCICRS group presented a higher immune cell infiltration and immune-related pathway activation, and hence a significant relation with the superior efficacy of pan-cancer immunotherapy. More importantly, the predictive value of MCICRS in MSS CRC patients receiving immunotherapy was also validated with an in-house cohort. Anti-PD-1 antibody combined with VEGFRi presented an improved clinical benefit in MSS CRC with manageable toxicity. MCICRS could serve as a robust and promising tool to predict clinical outcomes for individual MSS CRC patients receiving immunotherapy.
Humans ; Colorectal Neoplasms/drug therapy* ; Male ; Female ; Immunotherapy ; Middle Aged ; Aged ; Tumor Microenvironment/immunology* ; Retrospective Studies ; Microsatellite Instability ; Transcriptome ; Single-Cell Analysis ; Programmed Cell Death 1 Receptor/immunology* ; Gene Expression Profiling ; Immune Checkpoint Inhibitors/therapeutic use* ; Adult ; Receptors, Vascular Endothelial Growth Factor/antagonists & inhibitors*

Imprinted genes play a key role in regulating mammalian placental and embryonic development. Here, we generated glutaminyl-peptide cyclotransferase-knockout (Qpct^-/-) mice utilizing the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) platform and identified Qpct as a novel anti-angiogenic factor in regulating mouse placentation. Compared with Qpct^+/+ mice, placentae and embryos (Qpct^-/+ and Qpct^-/-) showed significant overgrowth at embryonic Day 12.5 (E12.5), E15.5, and E18.5. Using single-cell transcriptome analysis of 32 309 cells from Qpct^+/+ and Qpct^-/- mouse placentae, we identified 13 cell clusters via single-nucleus RNA sequencing (snRNA-seq) (8880 Qpct^+/+ and 13 577 Qpct^-/- cells) and 20 cell clusters via single-cell RNA sequencing (scRNA-seq) (6567 Qpct^+/+ and 3285 Qpct^-/- cells). Furthermore, we observed a global up-regulation of pro-angiogenic genes in the Qpct^-/- background. Immunohistochemistry assays revealed a notable increase in the number of blood vessels in the decidual and labyrinthine layers of E15.5 Qpct^-/+ and Qpct^-/- mice. Moreover, the elevation of multiple pairs of ligand-receptor interactions was observed in decidual cells, endothelial cells, and macrophages, promoting angiogenesis and inflammatory response. Our findings indicate that loss of maternal Qpct leads to altered phenotypic characteristics of placentae and embryos and promotes angiogenesis in murine placentae.
Animals ; Female ; Pregnancy ; Mice ; Placentation/genetics* ; Single-Cell Analysis ; Gene Expression Profiling ; Mice, Knockout ; Transcriptome ; Placenta/blood supply* ; Neovascularization, Pathologic/genetics* ; Genomic Imprinting ; Single-Cell Gene Expression Analysis ; Angiogenesis

Nasopharyngeal carcinoma（NPC） is a distinct type of head and neck cancer closely associated with Epstein-Barr virus（EBV） infection and exhibits significant geographic variations in its incidence. Despite recent advancements in radiotherapy techniques and precision medicine for NPC, the overall survival rate remains unsatisfactory due to tumor metastasis, recurrence, and drug resistance. Single-cell RNA sequencing（scRNA-seq） is an emerging technology that allows for the analysis of gene expression at single-cell resolution, providing a clearer understanding of tumor cell subpopulations, the evolutionary trajectory of tumor cells, and the functional roles and interactions of cells within the tumor microenvironment. This provides new ideas for the development of precision medicine in NPC. Here, we review the applications of scRNA-seq in exploring the mechanisms of NPC pathogenesis, tumor heterogeneity, the tumor microenvironment, drug resistance, and therapeutic response.
Humans ; Nasopharyngeal Neoplasms/genetics* ; Tumor Microenvironment ; Nasopharyngeal Carcinoma ; Single-Cell Analysis ; Sequence Analysis, RNA ; Precision Medicine ; Drug Resistance, Neoplasm ; Epstein-Barr Virus Infections ; Herpesvirus 4, Human ; Single-Cell Gene Expression Analysis

The oxytocin receptor (OXTR) has garnered increasing attention for its role in regulating both mature behaviors and brain development. It has been established that OXTR mediates a range of effects that are region-specific or period-specific. However, the current studies of OXTR expression patterns in mice only provide limited help due to limitations in resolution. Therefore, our objective was to generate a comprehensive, high-resolution spatiotemporal expression map of Oxtr mRNA across the entire developing mouse brain. We applied RNAscope in situ hybridization to investigate the spatiotemporal expression pattern of Oxtr in the brains of male mice at six distinct postnatal developmental stages (P7, P14, P21, P28, P42, P56). We provide detailed descriptions of Oxtr expression patterns in key brain regions, including the cortex, basal forebrain, hippocampus, and amygdaloid complex, with a focus on the precise localization of Oxtr⁺ cells and the variance of expression between different neurons. Furthermore, we identified some neuronal populations with high Oxtr expression levels that have been little studied, including glutamatergic neurons in the ventral dentate gyrus, Vgat⁺Oxtr⁺ cells in the basal forebrain, and GABAergic neurons in layers 4/5 of the cortex. Our study provides a novel perspective for understanding the distribution of Oxtr and encourages further investigations into its functions.
Animals ; Receptors, Oxytocin/metabolism* ; Male ; Brain/growth & development* ; Mice ; Mice, Inbred C57BL ; Neurons/metabolism* ; Single-Cell Analysis ; Gene Expression Regulation, Developmental ; RNA, Messenger/metabolism* ; Animals, Newborn

Dysregulated RNA splicing is a well-recognized characteristic of colorectal cancer (CRC); however, its intricacies remain obscure, partly due to challenges in profiling full-length transcript variants at the single-cell level. Here, we employ high-depth long-read scRNA-seq to define the full-length transcriptome of colorectal epithelial cells in 12 CRC patients, revealing extensive isoform diversities and splicing alterations. Cancer cells exhibited increased transcript complexity, with widespread 3'-UTR shortening and reduced intron retention. Distinct splicing regulation patterns were observed between intrinsic-consensus molecular subtypes (iCMS), with iCMS3 displaying even higher splicing factor activities and more pronounced 3'-UTR shortening. Furthermore, we revealed substantial shifts in isoform usage that result in alterations of protein sequences from the same gene with distinct carcinogenic effects during tumorigenesis of CRC. Allele-specific expression analysis revealed dominant mutant allele expression in key oncogenes and tumor suppressors. Moreover, mutated PPIG was linked to widespread splicing dysregulation, and functional validation experiments confirmed its critical role in modulating RNA splicing and tumor-associated processes. Our findings highlight the transcriptomic plasticity in CRC and suggest novel candidate targets for splicing-based therapeutic strategies.
Humans ; Colorectal Neoplasms/metabolism* ; RNA Isoforms/metabolism* ; Single-Cell Analysis ; RNA Splicing ; Gene Expression Regulation, Neoplastic ; RNA, Neoplasm/metabolism* ; Transcriptome

Ancient traditional Chinese medicine (TCM) doctrine says "The superior doctor prevents illnesses," pointing out preventative medicine as the ultimate goal for medical care. TCM recognizes that genetic predisposition and environmental and lifestyle influences contribute to diseases. It divides people into eight constitutions in addition to one normal/healthy kind. People with one of the eight subhealth constitutions are prone to develop different kinds of corresponding illnesses. The goal for this type of categorization is to help people take preemptive measures to prevent or delay disease onset. As the peripheral immune system through surveying the body, it can capture information from essentially all organs and reflect anomalies occurring in each organ. Thus, the detailed profiling of the peripheral immune-system function can generally reflect a person's overall heath state. In this study, we performed the single-cell RNA sequencing (scRNA-seq) of peripheral blood mononuclear cells (PBMCs) from individuals with Tanshi (phlegm dampness) constitution. They were prone to develop metabolic disorders including diabetes. scRNA-seq revealed greatly reduced mucosal-associated invariable T cell content and heightened TNFα-NFκB, JAK-STAT, and interferon signaling. These findings indicated heightened chronic inflammation, as well as increased hypoxia/apoptosis responses, likely resulting from frequent sleep apnea that Tanshi individuals experienced. Altogether, this pilot study demonstrated effectiveness in using scRNA-seq to reveal molecular-immunological bases for constitution categorization, thereby substantiating that preventative medicine originated from TCM.
Humans ; Leukocytes, Mononuclear/metabolism* ; Male ; Female ; Gene Expression Profiling ; Single-Cell Analysis ; Middle Aged ; Adult ; Medicine, Chinese Traditional ; Transcriptome ; Single-Cell Gene Expression Analysis

OBJECTIVE: To identify prognostic genes associated with lysosome-dependent cell death (LDCD) in patients with gastric cancer (GC). METHODS: Differentially expressed genes (DEGs) were identified using The Cancer Genome Atlas - Stomach Adenocarcinoma. Weighted gene co-expression network analysis was performed to identify the key module genes associated with LDCD score. Candidate genes were identified by DEGs and key module genes. Univariate Cox regression analysis, and least absolute shrinkage and selection operator regression and multivariate Cox regression analyses were performed for the selection of prognostic genes, and risk module was established. Subsequently, key cells were identified in the single-cell dataset (GSE183904), and prognostic gene expression was analyzed. Cell proliferation and migration were assessed using the Cell Counting Kit-8 assay and the wound healing assay. RESULTS: A total of 4,465 DEGs, 95 candidate genes, and 4 prognostic genes, including C19orf59, BATF2, TNFAIP2, and TNFSF18, were identified in the analysis. Receiver operating characteristic curves indicated the excellent predictive power of the risk model. Three key cell types (B cells, chief cells, and endothelial/pericyte cells) were identified in the GSE183904 dataset. C19orf59 and TNFAIP2 exhibited predominant expression in macrophage species, whereas TNFAIP2 evolved over time in endothelial/pericyte cells and chief cells. Functional experiments confirmed that interfering with C19orf59 inhibited proliferation and migration in GC cells. CONCLUSION C19orf59, BATF2, TNFAIP2, and TNFSF18 are prognostic genes associated with LDCD in GC. Furthermore, the risk model established in this study showed robust predictive power.
Stomach Neoplasms/pathology* ; Humans ; Prognosis ; Lysosomes/physiology* ; RNA-Seq ; Cell Death ; Single-Cell Analysis ; Gene Expression Regulation, Neoplastic ; Cell Proliferation ; Single-Cell Gene Expression Analysis