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

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

Objective： To investigate the relationship between the lipid metabolism-related gene sterol carrier protein 2(SCP2) and prostate cancer (PCa) from a multi-omics perspective using single-cell transcriptomes combined with Mendelian randomization. Methods： Single-cell transcriptome data of benign and malignant prostate tissues were obtained from GSE120716, GSE157703 and GSE141445 datasets, respectively. Integration, quality control and annotation were performed on the data to categorize the epithelial cells into high and low SCP2 expression groups, followed by further differential and trajectory analyses. Single nucleotide polymorphism (SNP) data for SCP2 expression quantitative trait loci (eQTL) were subsequently downloaded from Genotype-Tissue Expression (GTEx) and investigated from the PCa Society Cancer-Related Genomic Alteration Panel for the Investigation of Cancer-Related Alterations (PRACTICAL) to obtain PCa outcome data for Mendelian randomization analysis to validate the causal relationship between SCP2 and PCa. Results： High SCP2-expressing epithelial cells had higher energy metabolism and proliferation capacity with low immunotherapy response and metastatic tendency. Trajectory analysis showed that epithelial cells with high SCP2 expression may have a higher degree of malignancy, and SCP2 may be a key marker gene for differentiation of malignant epithelial cells in the prostate. Further Mendelian randomization results showed a significant causal relationship between SCP2 and PCa development (OR=1.045, 95% CI: 1.010 －1.083, P=0.011). Conclusion： By combining single-cell transcriptome and Mendelian randomization, the role of the lipid metabolism-related gene SCP2 in PCa development has been confirmed, and new targets and therapeutic directions for PCa treatment have been provided.
Humans ; Prostatic Neoplasms/genetics* ; Male ; Mendelian Randomization Analysis ; Polymorphism, Single Nucleotide ; Quantitative Trait Loci ; Single-Cell Analysis ; Sequence Analysis, RNA ; Carrier Proteins/genetics* ; Transcriptome ; Lipid Metabolism

The spatio-temporal heterogeneity of breast cell subsets forms the fundamental biological basis for physiological development and pathological progression, including tumorigenesis; however, its complex regulatory mechanisms are not yet fully elucidated. With its high-resolution capabilities, single-cell RNA sequencing (scRNA-seq) technology offers a powerful tool for dissecting this cellular heterogeneity. This technology enables the construction of high-precision breast cell atlases, the accurate identification of distinct cell subsets, and the reconstruction of differentiation trajectories from stem/progenitor cells to functional epithelial cells. By resolving the transcriptional regulatory networks that govern cell fate determination, intercellular communication patterns, and dynamic microenvironmental interactions, scRNA-seq has unveiled the molecular foundations of breast development and provided new perspectives on the pathogenesis of related diseases such as breast cancer and macromastia. Furthermore, scRNA-seq demonstrates significant potential for discovering early molecular markers of disease, deciphering tumor heterogeneity, and elucidating mechanisms of therapeutic resistance. The continued application of scRNA-seq for dissecting breast cell heterogeneity, combined with its integration with multi-modal data such as spatial omics, promises to provide critical evidence and new insights for revealing the molecular mechanisms of breast development-related diseases and for formulating precision therapeutic strategies.
Humans ; Single-Cell Analysis/methods* ; Female ; Breast Neoplasms/pathology* ; Sequence Analysis, RNA/methods* ; Breast/cytology*

OBJECTIVES: Keloids are fibrotic skin disorders characterized by excessive collagen deposition and a high recurrence rate, closely associated with inflammatory mediators. However, existing epidemiological studies are limited by confounding factors and reverse causality, making it difficult to establish causation. This study aims to investigate the causal relationship between circulating cytokines and keloids using Mendelian randomization analysis. METHODS: Significant single nucleotide polymorphisms (SNPs) associated with circulating cytokines (exposures) and keloids (outcomes) were extracted from genome-wide association study (GWAS) summary datasets. Eligible SNPs were selected as instrumental variables (IVs). Exposure data were derived from a cytokine GWAS including 8 293 Finnish participants, and outcome data from a keloid GWAS based on the UK Biobank. The inverse-variance weighted (IVW) method served as the primary analytical approach to estimate causal effects, supplemented by weighted median (WME), MR-Egger regression, and other sensitivity analyses. Horizontal pleiotropy was assessed using MR-Egger regression and the MR pleiotropy residual sum and outlier (MR-PRESSO) test, while Cochran's Q test evaluated heterogeneity. Leave-one-out analysis was used to verify robustness and consistency. A reverse MR analysis was also conducted, with keloid as the exposure and cytokines as outcomes, to rule out reverse causation. RESULTS: IVW analysis identified significant positive causal associations between two cytokines and keloids-macrophage migration inhibitory factor (MIF) [odds ratio (OR)=2.081, 95% confidence interval (CI) 1.219 to 3.552, P=0.007] and monocyte chemoattractant protein-1 (MCP-1) (OR=1.673, 95% CI 1.036 to 2.701, P=0.035). Conversely, stem cell factor (SCF) showed a negative causal relationship with keloids (OR=0.518, 95% CI 0.269 to 0.998, P=0.049). Results from the MR-Egger and weighted median analyses were consistent with IVW findings. No evidence of horizontal pleiotropy was observed (P>0.05). Except for interleukin-6 (P=0.014), no heterogeneity was detected in other cytokines. Leave-one-out analysis further confirmed the robustness of the causal associations. In reverse MR analysis, keloids were causally related only to β-nerve growth factor (beta-NGF) (OR=1.048, 95% CI 1.002 to 1.095, P=0.039), with no heterogeneity or pleiotropy detected in most cytokines (P>0.05). CONCLUSIONS MIF and MCP-1 exhibit positive causal associations with keloid formation, while SCF shows a negative causal relationship. These findings provide new evidence for the causal involvement of inflammatory cytokines in keloid pathogenesis and offer potential molecular targets for developing novel keloid therapies.
Humans ; Keloid/blood* ; Mendelian Randomization Analysis ; Cytokines/genetics* ; Polymorphism, Single Nucleotide ; Genome-Wide Association Study ; Chemokine CCL2/genetics* ; Interleukin-6/genetics* ; Macrophage Migration-Inhibitory Factors/genetics* ; Male ; Stem Cell Factor/blood* ; Female ; Intramolecular Oxidoreductases

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

Meniscus injuries are widespread and the available treatments do not offer enough healing potential. Here, we provide critical support for using pigs as a biological model for meniscal degeneration and the development of cutting-edge therapies in orthopedics. We present a single-cell transcriptome atlas of the meniscus, consisting of cell clusters corresponding to four major cell types: chondrocytes, endothelial cells, smooth muscle cells, and immune cells. Five distinct chondrocyte subclusters (CH0‒CH4) were annotated, of which only one was widespread in both the red and white zones, indicating a major difference in the cellular makeup of the zones. Subclusters distinct to the white zone appear responsible for cartilage-specific matrix deposition and protection against adverse microenvironmental factors, while those in the red zone exhibit characteristics of mesenchymal stem cells and are more likely to proliferate and migrate. Additionally, they induce remodeling actions in other chondrocyte subclusters and promote the proliferation and maturation of endothelial cells, inducing healing and vascularization processes. Considering that they have substantial remodeling capabilities, these subclusters should be of great interest for tissue engineering studies. We also show that the cellular makeup of the pig meniscus is comparable to that of humans, which supports the use of pigs as a model in orthopedic therapy development.
Animals ; Swine ; Chondrocytes/physiology* ; Single-Cell Analysis ; Meniscus/blood supply* ; Endothelial Cells/cytology* ; Transcriptome ; Mesenchymal Stem Cells/cytology* ; Neovascularization, Physiologic ; Cell Proliferation