OBJECTIVES: Metastasis is the primary cause of death in osteosarcoma, and current clinical treatments remain limited. BRD4, a key epigenetic regulator, has shown therapeutic promise in various cancers through its inhibition. However, the mechanistic role of BRD4 in osteosarcoma remains poorly understood. This study aims to elucidate the molecular mechanisms by which BRD4 regulate osteosarcoma progression and to explore novel therapeutic strategies. METHODS: Immunofluorescence was used to assess BRD4 expression levels in a tissue microarray containing 80 osteosarcoma samples from different patients. The Gene Expression Omnibus (GEO) dataset (GSE42352, containing survival data from 88 osteosarcoma patients) was downloaded to perform Kaplan-Meier survival analysis based on BRD4 gene expression levels. In vivo, an orthotopic intramedullary osteosarcoma model was established using HOS cells in C57 mice, followed by treatment with varying doses of the BRD4 inhibitor (+)-JQ1. Micro-CT, 3D reconstruction of bone tissue, and HE staining were employed to evaluate pathological changes in bone and intestinal lymph nodes. In vitro, cell viability was measured using the methyl thiazolyl tetrazolium (MTT) assay, while colony formation and Transwell assays assessed proliferative and invasive capacities. Chromatin-bound BRD4 was analyzed via co-immunoprecipitation combined with mass spectrometry (Co-IP/MS), and O-GlcNAc glycosylation sites and glycan chains of BRD4 were identified using Co-IP with Nano-LC MS/MS. Real-time PCR and Western blotting were used to analyze the relative mRNA and protein expression levels of target genes, respectively. RESULTS: BRD4 was positively expressed in 61.25% (49/80) of osteosarcoma tissues. Patients with high BRD4 expression exhibited significantly shorter survival times (P<0.05). In the orthotopic mouse model, intervention with (+)-JQ1, a potent and commonly used BETi, significantly inhibited tumor growth in vivo and reduced bone destruction (P<0.05). (+)-JQ1 treatment significantly suppressed the proliferation (P<0.001), invasion (P<0.001), and migration (P<0.05) of HOS cells. In osteosarcoma cells, BRD4 exhibited O-GlcNAc modifications at both N- and C- C-termini, particularly at Thr73, which is essential for protein stability. This modification also contributed to the activation of the EGFR tyrosine kinase inhibitor resistance pathway (KEGG Pathway: hsa01521). (+)-JQ1 treatment displaced BRD4 from enhancers and downregulated the transcription of pathway-related genes, such as EGFR and PDGFC, thereby suppressing the malignant behavior of osteosarcoma cells. CONCLUSIONS BRD4 promotes osteosarcoma progression via O-GlcNAc modification at Thr73 and plays a crucial role in tumor growth and metastasis.
Osteosarcoma/drug therapy* ; Humans ; Transcription Factors/metabolism* ; Animals ; Cell Cycle Proteins ; Mice ; Bone Neoplasms/drug therapy* ; Azepines/pharmacology* ; Cell Line, Tumor ; Cell Proliferation/drug effects* ; Triazoles/pharmacology* ; Mice, Inbred C57BL ; Nuclear Proteins/metabolism* ; Gene Expression Regulation, Neoplastic ; Male ; Bromodomain Containing Proteins

OBJECTIVES: Osteoarthritis (OA) and sarcopenia are significant health concerns in the elderly, substantially impacting their daily activities and quality of life. However, the relationship between them remains poorly understood. This study aims to uncover common biomarkers and pathways associated with both OA and sarcopenia. METHODS: Gene expression profiles related to OA and sarcopenia were retrieved from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between disease and control groups were identified using R software. Common DEGs were extracted via Venn diagram analysis. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to identify biological processes and pathways associated with shared DEGs. Protein-protein interaction (PPI) networks were constructed, and candidate hub genes were ranked using the maximal clique centrality (MCC) algorithm. Further validation of hub gene expression was performed using 2 independent datasets. Receiver operating characteristic (ROC) curve analysis was used to evaluate the predictive value of key genes for OA and sarcopenia. Mouse models of OA and sarcopenia were established. Hematoxylin-eosin and Safranin O/Fast Green staining were used to validate the OA model. The sarcopenia model was validated via rotarod testing and quadriceps muscle mass measurement. Real-time reverse transcription PCR (real-time RT-PCR) was employed to assess the mRNA expression levels of candidate key genes in both models. Gene set enrichment analysis (GSEA) was conducted to identify pathways associated with the selected shared key genes in both diseases. RESULTS: A total of 89 common DEGs were identified in the gene expression profiles of OA and sarcopenia, including 76 upregulated and 13 downregulated genes. These 89 DEGs were significantly enriched in protein digestion and absorption, the PI3K-Akt signaling pathway, and extracellular matrix-receptor interaction. PPI network analysis and MCC algorithm analysis of the 89 common DEGs identified the top 17 candidate hub genes. Based on the differential expression analysis of these 17 candidate hub genes in the validation datasets, AEBP1 and COL8A2 were ultimately selected as the common key genes for both diseases, both of which showed a significant upregulation trend in the disease groups (all P<0.05). The value of area under the curve (AUC) for AEBP1 and COL8A2 in the OA and sarcopenia datasets were all greater than 0.7, indicating that both genes have potential value in predicting OA and sarcopenia. Real-time RT-PCR results showed that the mRNA expression levels of AEBP1 and COL8A2 were significantly upregulated in the disease groups (all P<0.05), consistent with the results observed in the bioinformatics analysis. GSEA revealed that AEBP1 and COL8A2 were closely related to extracellular matrix-receptor interaction, ribosome, and oxidative phosphorylation in OA and sarcopenia. CONCLUSIONS AEBP1 and COL8A2 have the potential to serve as common biomarkers for OA and sarcopenia. The extracellular matrix-receptor interaction pathway may represent a potential target for the prevention and treatment of both OA and sarcopenia.
Sarcopenia/genetics* ; Osteoarthritis/genetics* ; Computational Biology/methods* ; Humans ; Protein Interaction Maps/genetics* ; Animals ; Mice ; Gene Expression Profiling ; Gene Ontology ; Transcriptome ; Male ; Signal Transduction/genetics* ; Gene Regulatory Networks

OBJECTIVES: Polycystic ovary syndrome (PCOS) is a common endocrine disorder that affects women's health. This study aims to investigate gene and transcription factor (TF) expression differences between PCOS patients and healthy individuals using bioinformatics approaches, and to verify the function of key transcription factors, with the goal of providing new insights into the pathogenesis of PCOS. METHODS: Differentially expressed genes (DEGs) and differentially expressed transcription factors (DETFs) between PCOS patients and controls were identified from the RNA sequencing dataset GSE168404 using bioinformatics methods. Functional enrichment analysis was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The expression and function of core transcription factors were further validated in ovarian tissues of PCOS model mice and control mice using Western blotting and reverse transcription quantitative polymerase chain reaction (RT-qPCR). RESULTS: A total of 332 DEGs were identified between PCOS patients and controls, including 259 upregulated and 73 downregulated genes in the PCOS group. 19 DETFs were further screened, of which 16 were upregulated and 3 were downregulated in PCOS. The upregulated DETFs (including TFCP2L1, DACH1, ESR2, AFF3, SMAD9, ZNF331, HOPX,ATOH8, HIF3α, DPF3, HOXC4, HES1, ID1, JDP2, SOX4, and ID3) were primarily associated with lipid metabolism, development, and cell adhesion. Protein and mRNA expression analysis in PCOS model mice revealed significantly decreased levels of hypoxia-inducible factor (HIF) 1α and HIF2α, and significantly increased expression of HIF3α compared to control mice (all P<0.001). CONCLUSIONS Significant differences in gene and TF expression exist between PCOS patients and healthy individuals. HIF-3α may play a crucial role in PCOS and could serve as a novel biomarker for diagnosis and a potential therapeutic target.
Polycystic Ovary Syndrome/metabolism* ; Female ; Humans ; Animals ; Mice ; Transcription Factors/metabolism* ; Computational Biology ; Gene Expression Profiling ; Adult

OBJECTIVES: Multiple myeloma (MM) is a highly heterogeneous hematologic malignancy, with disease progression driven by cytogenetic abnormalities and a complex bone marrow microenvironment. This study aims to construct a prognostic model for MM based on transcriptomic data and lipid metabolism related genes (LRGs), and to identify potential drug targets for high-risk patients to support clinical decision-making. METHODS: In this study, 2 transcriptomic datasets covering 985 newly diagnosed MM patients were retrieved from the Gene Expression Omnibus (GEO) database. Univariate Cox regression and 101 machine learning algorithms were used for gene selection. An LRG-based prognostic model was constructed using Stepwise Cox (both directions) and random survival forest (RSF) algorithms. The association between the prognostic score and clinical events was evaluated, and model performance was assessed using time-dependent receiver operating characteristic (ROC) curves and the C-index. The added predictive value of combining prognostic scores with clinical variables and staging systems was also analyzed. Differentially expressed genes between high- and low-risk groups were identified using limma and clusterProfiler and subjected to pathway enrichment analysis. Drug sensitivity analysis was conducted using the Genomics of Drug Sensitivity in Cancer (GDSC) database and oncoPredict to identify potential therapeutic targets for high-risk patients. The functional role of key LRGs in the model was validated via in vitro cell experiments. RESULTS: An LRG-based prognostic model (LRG17) was successfully developed using transcriptomic data and machine learning. The model demonstrated robust predictive performance, with area under the curve (AUC) values of 0.962, 0.912, and 0.842 for 3-, 5-, and 7-year survival, respectively. Patients were stratified into high- and low-risk groups, with high-risk patients showing significantly shorter overall survival (OS) and event-free survival (EFS) (both P<0.001) and worse clinical profiles (e.g., lower albumin, higher β2-microglobulin and lactate dehydrogenase levels). Enrichment analysis revealed that high-risk patients were significantly enriched for pathways related to chromosome segregation and mitosis, whereas low-risk patients were enriched for immune response and immune cell activation pathways. Drug screening suggested that AURKA inhibitor BMS-754807 and FGFR3 inhibitor I-BET-762 may be more effective in high-risk patients. Functional assays demonstrated that silencing of key LRG PLA2G4A significantly inhibited cell viability and induced apoptosis. CONCLUSIONS LRGs serve as promising biomarkers for prognosis prediction and risk stratification in MM. The overexpression of chromosomal instability-related and high-risk genetic event-associated genes in high-risk patients may explain their poorer outcomes. Given the observed resistance to bortezomib and lenalidomide in high-risk patients, combination therapies involving BMS-754807 or I-BET-762 may represent effective alternatives.
Humans ; Multiple Myeloma/mortality* ; Prognosis ; Lipid Metabolism/genetics* ; Transcriptome ; Machine Learning ; Male ; Female ; Gene Expression Profiling ; Algorithms

OBJECTIVES: Increasing detection of low-risk papillary thyroid carcinoma (PTC) is associated with overdiagnosis and overtreatment. N6-methyladenosine (m⁶A)-mediated microRNA (miRNA) dysregulation plays a critical role in tumor metastasis and progression. However, the functional role of m⁶A-miRNAs in PTC remains unclear. This study aims to elucidate the regulatory mechanism of m⁶A-miR-139-5p expression in PTC, determine its association with PTC metastasis, and evaluate its potential as a diagnostic biomarker for PTC metastasis, thereby providing experimental evidence for precision diagnosis and therapy. METHODS: Expression profiles of m⁶A-miRNAs were compared between the The Cancer Genome Atlas (TCGA) and GSE130512 cohorts to identify metastasis-associated candidates. Clinical specimens from 13 metastasis and 18 non-metastasis PTC patients were analyzed to assess m⁶A-miR-139-5p expression and its correlation with metastasis. Functional experiments were conducted to investigate the effect of fat mass and obesity-associated protein (FTO) on pri-miR-139 methylation and processing, clarifying its regulatory role in miR-139-5p expression. In TPC-1 cells, MTT assays were performed to evaluate whether miR-139-5p overexpression could counteract FTO-mediated cell proliferation. Transwell invasion assays were used to determine the impact of miR-139-5p on PTC cell invasion, exploring whether it functions through the ZEB1/E-cadherin axis. RESULTS: By comparing TCGA and GSE130512 cohorts, it was found that circulating m⁶A-miR-139-5p could serve as a biological indicator for detecting PTC metastasis. Detection of 13 metastatic and 18 non-metastatic clinical specimens showed that FTO inhibited the processing of pri-miR-139 by reducing its methylation level, leading to the dysregulation of miR-139-5p in PTC (P<0.05). In TPC-1 cells, MTT assay showed that overexpression of miR-139-5p could partially reverse FTO overexpression-mediated cell proliferation (P<0.05). In addition, miR-139-5p inhibited the invasive ability of PTC cells by targeting the ZEB1/E-cadherin axis, while FTO overexpression could partially weaken this inhibitory effect. CONCLUSIONS Circulating miR-139-5p can be a potential marker for evaluating PTC metastasis. FTO affects the expression and function of miR-139-5p by regulating m⁶A modification of pri-miR-139, but its clinical value needs further verification.
Humans ; MicroRNAs/metabolism* ; Thyroid Cancer, Papillary/metabolism* ; Alpha-Ketoglutarate-Dependent Dioxygenase FTO/metabolism* ; Thyroid Neoplasms/metabolism* ; Cell Line, Tumor ; Neoplasm Metastasis ; Adenosine/genetics* ; Gene Expression Regulation, Neoplastic ; Female ; Male ; Cadherins/metabolism* ; Cell Proliferation ; Zinc Finger E-box-Binding Homeobox 1/genetics*

OBJECTIVES: Bladder cancer is a common malignancy with high incidence and poor prognosis. N⁶-methyladenosine (m⁶A) modification is widely involved in diverse physiological processes, among which the m⁶A recognition protein YTH N⁶-methyladenosine RNA binding protein F2 (YTHDF2) plays a crucial role in bladder cancer progression. This study aims to elucidate the molecular mechanism by which O-linked N-acetylglucosamine (O-GlcNAc) modification of YTHDF2 regulates its downstream target, period circadian regulator 1 (PER1), thereby promoting bladder cancer cell proliferation. METHODS: Expression of YTHDF2 in bladder cancer was predicted using The Cancer Genome Atlas (TCGA). Twenty paired bladder cancer and adjacent normal tissues were collected at the clinical level. Normal bladder epithelial cells (SV-HUC-1) and bladder cancer cell lines (T24, 5637, EJ-1, SW780, BIU-87) were examined by quantitative real-time PCR (RT-qPCR), Western blotting, and immunohistochemistry for expression of YTHDF2, PER1, and proliferation-related proteins [proliferating cell nuclear antigen (PCNA), minichromosome maintenance complex component 2 (MCM2), Cyclin D1]. YTHDF2 was silenced in 5637 and SW780 cells, and cell proliferation was assessed by Cell Counting Kit-8 (CCK-8), colony formation, and EdU assays. Bioinformatics was used to predict glycosylation sites of YTHDF2, and immunoprecipitation (IP) was performed to detect O-GlcNAc modification levels of YTHDF2 in tissues and cells. Bladder cancer cells were treated with DMSO, OSMI-1 (O-GlcNAc inhibitor), or Thiamet G (O-GlcNAc activator), followed by cycloheximide (CHX), to assess YTHDF2 ubiquitination by IP. YTHDF2 knockdown and Thiamet G treatment were further used to evaluate PER1 mRNA stability, PER1 m⁶A modification, and cell proliferation. TCGA was used to predict PER1 expression in tissues; SRAMP predicted potential PER1 m⁶A sites. Methylated RNA immunoprecipitation (MeRIP) assays measured PER1 m⁶A modification. Finally, the effects of knocking down YTHDF2 and PER1 on 5637 and SW780 cell proliferation were assessed. RESULTS: YTHDF2 expression was significantly upregulated in bladder cancer tissues compared with adjacent tissues (mRNA: 2.5-fold; protein: 2-fold), which O-GlcNAc modification levels increased 3.5-fold (P<0.001). YTHDF2 was upregulated in bladder cancer cell lines, and its knockdown suppressed cell viability (P<0.001), downregulated PCNA, MCM2, and CyclinD1 (all P<0.05), reduced colony numbers 3-fold (P<0.01), and inhibited proliferation. YTHDF2 exhibited elevated O-GlcNAc modification in cancer cells. OSMI-1 reduced YTHDF2 protein stability (P<0.01) and enhanced ubiquitination, while Thiamet G exerted opposite effects (P<0.001). Thiamet G reversed the proliferation-suppressive effects of YTHDF2 knockdown, promoting cell proliferation (P<0.01) and upregulating PCNA, MCM2, and CyclinD1 (all P<0.05). Mechanistically, YTHDF2 targeted PER1 via m⁶A recognition, promoting PER1 mRNA degradation. Rescue experiments showed that PER1 knockdown reversed the inhibitory effect of YTHDF2 knockdown on cell proliferation, upregulated PCNA, MCM2, and Cyclin D1 (all P<0.05), and promoted bladder cancer cell proliferation (P<0.001). CONCLUSIONS O-GlcNAc modification YTHDF2 promotes bladder cancer development by downregulating the tumor suppressor gene PER1 through m⁶A-mediated post-transcriptional regulation.
Humans ; Urinary Bladder Neoplasms/metabolism* ; RNA-Binding Proteins/genetics* ; Cell Proliferation ; Cell Line, Tumor ; Disease Progression ; Acetylglucosamine/metabolism* ; Adenosine/metabolism* ; Gene Expression Regulation, Neoplastic ; Genes, Tumor Suppressor

Drugs and pesticide residues in broiler feed can compromise the therapeutic and production benefits of antibiotic (ANT) application and affect gene expression. In this study, we analyzed the expression of 13 key pancreatic genes and blood physiology parameters after administering one maximum residue limit of herbicide glyphosate (GLY), two ANTs, and one anticoccidial drug (AD). A total of 260 Ross 308 broilers aged 1‍-‍40 d were divided into the following four groups of 65 birds each: control group, which was fed the main diet (MD), and three experimental groups, which were fed MD supplemented with GLY, GLY+ANTs (enrofloxacin and colistin methanesulfonate), and GLY+AD (ammonium maduramicin), respectively. The results showed that the addition of GLY, GLY+ANTs, and GLY+AD caused significant changes in the expression of several genes of physiological and economic importance. In particular, genes related to inflammation and apoptosis (interleukin 6 (IL6), prostaglandin-endoperoxide synthase 2 (PTGS2), and caspase 6 (CASP6)) were downregulated by up to 99.1%, and those related to antioxidant protection (catalase (CAT), superoxide dismutase 1 (SOD1) and peroxiredoxin 6 (PRDX6)) by up to 98.6%, compared to controls. There was also a significant decline in the values of immunological characteristics in the blood serum observed in the experimental groups, and certain changes in gene expression were concordant with changes in the functioning of the pancreas and blood. The changes revealed in gene expression and blood indices in response to GLY, ANTs, and AD provide insights into the possible mechanisms of action of these agents at the molecular level. Specifically, these changes may be indicative of physiological mechanisms to overcome the negative effects of GLY, GLY+ANTs, and GLY+AD in broilers.
Animals ; Glyphosate ; Glycine/administration & dosage* ; Chickens/blood* ; Pancreas/metabolism* ; Anti-Bacterial Agents/pharmacology* ; Animal Feed ; Gene Expression/drug effects* ; Herbicides

Given that ovarian stimulation is vital for assisted reproductive technology (ART) and results in elevated serum estrogen levels, exploring the impact of elevated estrogen exposure on oocytes and embryos is necessary. We investigated the effects of various ovarian stimulation treatments on oocyte and embryo morphology and gene expression using a mouse model and estrogen-treated mouse embryonic stem cells (mESCs). Female C57BL/6J mice were subjected to two types of conventional ovarian stimulation and ovarian hyperstimulation; mice treated with only normal saline served as controls. Hyperstimulation resulted in high serum estrogen levels, enlarged ovaries, an increased number of aberrant oocytes, and decreased embryo formation. The messenger RNA (mRNA)‍-sequencing of oocytes revealed the dysregulated expression of lysine-specific demethylase 6b (Kdm6b), which may be a key factor indicating hyperstimulation-induced aberrant oocytes and embryos. In vitro, Kdm6b expression was downregulated in mESCs treated with high-dose estrogen; treatment with an estrogen receptor antagonist could reverse this downregulated expression level. Furthermore, treatment with high-dose estrogen resulted in the upregulated expression of histone H3 lysine 27 trimethylation (H3K27me3) and phosphorylated H2A histone family member X (γ-H2AX). Notably, knockdown of Kdm6b and high estrogen levels hindered the formation of embryoid bodies, with a concomitant increase in the expression of H3K27me3 and γ-H2AX. Collectively, our findings revealed that hyperstimulation-induced high-dose estrogen could impair the demethylation of H3K27me3 by reducing Kdm6b expression. Accordingly, Kdm6b could be a promising marker for clinically predicting ART outcomes in patients with ovarian hyperstimulation syndrome.
Female ; Mice ; Demethylation/drug effects* ; Embryonic Stem Cells ; Estrogens/administration & dosage* ; Gene Expression/drug effects* ; Histones/metabolism* ; Jumonji Domain-Containing Histone Demethylases/metabolism* ; Mice, Inbred C57BL ; Oocytes ; Ovary/drug effects* ; Reproductive Techniques, Assisted ; Animals

BACKGROUND: PANoptosis has the features of pyroptosis, apoptosis, and necroptosis. Numerous studies have confirmed the diverse roles of various types of cell death in acute liver failure (ALF), but limited attention has been given to the crosstalk among them. In this study, we aimed to explore the role of PANoptosis in ALF and uncover new targets for its prevention or treatment. METHODS: Three ALF-related datasets (GSE14668, GSE62029, and GSE74000) were downloaded from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs). Hub genes were identified through intersecting DEGs, genes obtained from weighted gene co-expression network analysis (WGCNA), and genes related to PANoptosis. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), protein‍‒‍protein interaction (PPI) analyses and gene set enrichment analysis (GSEA) were performed to determine functional roles. Verification was performed using an ALF mouse model. RESULTS: Our results showed that expression of seven hub genes (B-cell lymphoma-2-modifying factor (BMF), B-cell lymphoma-2-interacting protein 3-like (BNIP3L), Caspase-1 (CASP1), receptor-interacting protein kinase 3 (RIPK3), uveal autoantigen with coiled-coil domains and ankyrin repeats protein (UACA), uncoordinated-5 homolog B receptor (UNC5B), and Z-DNA-binding protein 1 (ZBP1)) was up-regulated in liver samples of patients. However, in the ALF mouse model, the expression of BNIP3L, RIPK3, phosphorylated RIPK3 (P-RIPK3), UACA, and cleaved caspase-1 was up-regulated, while the expression of CASP1 and UNC5B was down-regulated. The expression of ZBP1 and BMF increased only during the development of ALF, and there was no significant change in the end stage. Immunofluorescence of mouse liver tissue showed that macrophages expressed all seven markers. Western blot results showed that pyroptosis, apoptosis, and necroptosis were always involved in lipopolysaccharide (LPS)/ d-galactosamine (d-gal)‍-induced ALF mice. The ALF cell model showed that bone marrow-derived macrophages (BMDMs) form PANoptosomes after LPS stimulation. CONCLUSIONS: Our results suggest that PANoptosis of macrophages promotes the development of ALF. The seven new ALF biomarkers identified and validated in this study may contribute to further investigation of diagnostic markers or novel therapeutic targets of ALF.
Animals ; Liver Failure, Acute/genetics* ; Computational Biology ; Mice ; Pyroptosis/genetics* ; Humans ; Protein Interaction Maps ; Apoptosis/genetics* ; Necroptosis/genetics* ; Gene Regulatory Networks ; Gene Ontology ; Gene Expression Profiling ; Disease Models, Animal

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*