1.Construction and validation of a prognostic model for colon cancer based on anoikis-related genes
Tao ZHANG ; Ziyao LI ; Yingying SUN ; Boyang LI ; Zhao WANG ; Zhifu YANG
Cancer Research and Clinic 2025;37(1):55-63
Objective:To construct and validate a prognostic model of colon cancer based on differentially expressed anoikis-related genes, and to preliminarily investigate the relationship between anoikis-related genes and the tumor immune microenvironment of colon cancer.Methods:A total of 472 cancer tissues samples of patients with colon cancer, RNA sequencing data and clinical data of 41 normal tissues samples were downloaded from the Cancer Genome Atlas (TCGA) database between the establishment time and July in 2024. A total of 919 genes related to anoikis were screened out from GeneCards database, and the common genes were selected from the RNA sequencing gene datasets of colon cancer and normal colon tissues in the TCGA database, among which the differentially expressed anoikis-related genes of colon cancer and normal colon tissues were screened out based on P < 0.05. Furthermore, genes related to the prognosis of 446 colon cancer patients with prognostic data in the TCGA database were screened by using univariate Cox proportional risk model; the genes with P < 0.05 were further screened out and a colon cancer prognosis model was constructed by using LASSO-Cox proportional risk model. The risk score of the above 446 colon cancer patients in the TCGA database was calculated according to the prognostic model, and the patients were divided into high-risk (≥ median value) group and low-risk (< median value) group according to the median risk score, and the overall survival of the 2 groups was analyzed by using the Kaplan-Meier method. The risk score based on R software-based time ROC program package was used to predict 1-year, 2-year, 3-year overall survival therapeutic efficacy of colon cancer patients in the TCGA database. According to the median risk score of colon cancer patients in the TCGA database, the patients in the International Cancer Genome Consortium (ICGC) database were divided into high-risk group and low-risk group. Kaplan-Meier method and receiver operating characteristic (ROC) curve were used to verify the predictive effect of the prognostic model. The differentially expressed genes between low-risk group and high-risk group stratified by prognostic model risk score in the TCGA database were used to perform single sample gene set enrichment analysis (ssGESA) of immune cells and immune function by using R software related programs. The differences in risk scores of patients with different immunophenotypes (including inflammator response type, wound healing type, interferon gamma dominant type and lymphocyte depletion type) were compared; and correlation analysis of infiltration and risk scores between immune cells and stromal cells in tumor microenvironment was made. Based on the tumor immune function and rejection (TIDE) database, the relationship between the prognostic model risk score and programmed death receptor ligand 1 (PD-L1) gene expression level was analyzed. Results:Based on anoikis-related genes in the GeneCards database, 236 differentially expressed anoikis-related genes between colon cancer tissues and normal tissues were obtained from the TCGA database. LASSO Cox regression was applied to establish a prognostic model constructed by 7 differentially expressed anoikis-related genes in cancer tissues and normal colon tissues related to the prognosis of colon cancer. Risk score = 0.366×TIMP1-0.404×NAT1+0.207×LTB4R2+0.075×INHBB+0.140×CD36-0.109×MMP3+2.994×OFCC1. The median risk score of 446 colon cancer patients in the TCGA database was 1.754 719 545. Survival analysis showed that the overall survival of colon cancer patients in high-risk group of the TCGA database was worse than that in low-risk group ( P < 0.001); ROC curve analysis showed that the area under the curve for predicting 1-year, 2-year and 3-year overall survival of patients in the TCGA database based on the prognostic model risk score was 0.705, 0.731 and 0.723, respectively. Kaplan-Meier method analysis showed that in the ICGC database, the overall survival of colon cancer patients in high-risk group was worse than that in low-risk group ( P = 0.041); ROC curve analysis showed that the area under the curve of prognostic model risk score for predicting 1-year and 2-year overall survival of colon cancer patients in the ICGC database was 0.663 and 0.966, respectively. ssGESA analysis showed that macrophage level in high-risk group was higher than that in low-risk group, helper T (Th) 1 cell and Th2 cell levels in high-risk group were lower than those in low-risk group (all P < 0.01). In terms of immune function, the cell killing activity and histocompatibility complex Ⅰ level in high-risk group were lower than those in low-risk group, and type Ⅱ interferon response score in high-risk group was higher than that in low-risk group (all P < 0.05). The analysis of immunophenotype showed that the risk score of inflammatory response type was higher than that of wound healing type ( P < 0.05), and there was no statistically significant difference in risk score between the other 2 types (all P > 0.05). Risk score was positively correlated with stromal cell infiltration score ( R = 0.340, P < 0.001) and immune cell infiltration score ( R = 0.148, P < 0.05); the expression level of PD-L1 in high-risk group was higher than that in low-risk group in the TCGA database ( P = 0.048), and the expression level of PD-L1 was positively correlated with risk score ( R = 0.130, P = 0.009). Conclusions:A prognostic model of colon cancer constructed by anoikis-related genes can better predict the prognosis of colon cancer patients, and anoikis-related genes may play an important role in tumor immunity of colon cancer.
2.Anti-central-fatigue effect of maca via mitochondrial biogenesis via the AMPK/SIRT1/PGC-1α pathway in rats
Wenhuan YAO ; Wen ZHOU ; Yaxuan LI ; Ziyao LI ; Jing ZHANG ; Shibo LYU ; Hui LI
Chinese Journal of Comparative Medicine 2025;35(7):36-43
Objective To examine the anti-central-fatigue function of maca and its underlying mechanism.Methods Forty Sprague-Dawley rats were divided randomly into a negative control group,model control group,and low-,medium-,and high-dose maca groups(0.6,1.2,and 2.4 g/kg·body weight).Rats in all groups except the negative control group were subjected to multi-factor stimulation,including cold-water swimming,sleep deprivation,restraining,and tail-clamping,to establish central fatigue rat models.Rats in the low-,medium-,and high-dose maca groups received 0.6,1.2,or 2.4 g/kg maca,respectively,by gavage for 35 days.Behavioral testing was carried out using the Morris water-maze,sucrose-preference,and tail-suspension tests.Markers of oxidative stress in the hippocampus,including superoxide dismutase(SOD),malondialdehyde(MDA),and catalase(CAT),were detected using test kits.Proteins connected with the AMP-activated protein kinase(AMPK)/sirtuin 1(SIRT1)/peroxisome proliferator-activated receptor γ coactivator 1-α(PGC-1α)signaling pathway in the hippocampus were detected by Western blot.Results Rats in the low-,medium-,and high-dose maca groups spent significantly more time in the target quadrant compared with the model control group(P<0.05 or P<0.01),but there was no significant dose-effect relationship.Rats in the medium-and high-dose maca groups showed decreased escape latency(P<0.05),increased time crossing the platform location(P<0.05),increased sucrose preference(P<0.05),decreased tail suspension time(P<0.05),increased the activities of CAT(P<0.01)and SOD(P<0.05),and decreased MDA content(P<0.01).Rats in the low-,medium-,and high-dose maca groups also showed significantly increased protein expression levels of AMPK and nuclear respiratory factor 1(P<0.01 or P<0.05),but no significant dose-effect relationship was observed.Rats in the medium-and high-dose maca groups showed increased protein expression of PGC-1α(P<0.05 or P<0.01),and rats in the high-dose maca group showed increased protein expression of SIRT1 and mitochondrial transcription factor A(P<0.05 or P<0.01).Conclusions Maca can improve the indicators of central fatigue in rats,determined by behavioral testing and oxidative stress-related factors.The underlying mechanism may be related to its regulatory effects on the AMPK/SIRT1/PGC-1α signaling pathway.
3.Work Memory Impairment in Patients with Chronic Fatigue Syndrome and Spleen Deficiency
Tian ZHOU ; Yunhe ZHANG ; Ziyao WU ; Sitong FENG ; Yanzhe NING ; Hongxiao JIA
World Science and Technology-Modernization of Traditional Chinese Medicine 2025;27(11):3148-3156
Objective To characterize working memory performance in patients with chronic fatigue syndrome(CFS)and spleen-deficiency syndrome and to examine its associations with clinical symptoms by Sternberg working memory task(SWMT).Methods 31 CFS patients meeting both CDC-1994 criteria and consensus criteria for spleen-deficiency pattern were recruited from outpatient clinics and universities from September 2022 and June 2025.31 healthy controls were also recruited based on age,sex,and education.All subjects completed the SWMT.Group differences were analyzed.Within the CFS cohort,reaction time(RT)was correlated with scores on the checklist individual strength(CIS),36-item short-form health survey(sf-36),and fatigue scale-14(FS-14).Mediation was examined.Results RT lengthened with increasing memory load in both groups.CFS patients displayed slower RTs than controls in the baseline and 6-digit set(P<0.05).The 3-digit RT difference,though not significant(P>0.05),yielded a medium effect size(r=0.36).Accuracy did not differ between two groups.Among CFS patients,3-digit RT correlated positively with CIS total and the 4 sub-scale scores.6-digit RT correlated with the SF-36 health-transition dimension(r=0.396,P=0.027).CIS and FS-14 scores directly impaired SF-36 social functioning without working-memory mediating.Conclusion CFS patients with spleen-deficiency exhibit slowed processing speed rather than capacity loss.The close link between working-memory slowing and fatigue suggests a distinct neural basis.These results support the traditional concept"the spleen stores Yi"and integrate TCM pattern differentiation with modern cognitive neuroscience in CFS.
4.Brain functional magnetic resonance imaging changes in patients with depression after acupuncture treatment: ALE meta-analysis
Jing TANG ; Juan YU ; Ziyao GENG ; Yating YUE ; Jingxia ZHANG
Chinese Journal of Behavioral Medicine and Brain Science 2025;34(11):1005-1011
Objective:Through the meta-analysis of activation likelihood estimation (ALE), the functional magnetic resonance imaging research evidence of acupuncture in the treatment of depression was systematically integrated to analyze the effect of acupuncture on brain activity in patients with depression and systematically reveal its potential neural mechanism.Methods:The literatures which used functional magnetic resonance imaging to explore brain image changes of patients with depression after receiving acupuncture intervention were retrieved from CNKI, VIP database, Wanfang database, PubMed, Embase, EBSCO host MEDLINE and Web of Science database. The retrieval time is from the establishment of the database to March 2025. Data processing was performed using the ALE method by the Ginger ALE 2.3 software. The coordinates of the subjects included in the study were extracted, and the brain regions with abnormal spontaneous brain activity in patients with depression were integrated and analyzed.Results:A total of 28 articles involving 770 subjects were included. Depressed patients receiving acupuncture exhibited increased brain activity in the bilateral medial frontal gyrus and bilateral anterior cingulate cortex. Exploratory analysis revealed that depressed patients receiving acupuncture exhibited reduced activation in the left posterior cerebellar lobe (MNI: x, y, z = -6, -68, -21; x, y, z = -9, -81, -9) compared to healthy controls (voxel sum=5 824 mm 3,P<0.05). Conclusion:Patients with depression who received acupuncture intervention exhibited relatively consistent activation in brain regions compared to depressed patients who did not receive acupuncture and healthy individuals. Compared to healthy individuals, patients with depression who received acupuncture showed reduced activation in the left posterior cerebellum. It is suggested that acupuncture intervention may play a therapeutic role by regulating the functional state of emotion-related brain regions and brain networks.
5.Changes in S100A8/9 and NLRP3/Caspase-1/interleukin-1β pathway in kidney-aging rats induced by D-galactose
Dandan FENG ; Ying ZHOU ; Ziyao PANG ; Yueqin CAI ; Chu CHEN ; Jianyou ZHANG ; Dejun WANG
Acta Laboratorium Animalis Scientia Sinica 2025;33(6):823-835
Objective To investigate changes in the pro-inflammatory mediator S100A8/9 and NLRP3/Caspase-1/IL-1β pathway in a rat kidney-aging model induced by D-galactose.Methods Twelve SD rats were divided into control and D-galactose groups,and injected subcutaneously in the back of the neck with D-galactose(150 mg/kg)to establish a rat model of kidney aging.Kidney samples were collected under anesthesia after 8 weeks.Kidneys were stained for senescence-associated beta-galactosidase(SA-β-Gal),mRNA expression levels of the aging-related genes p21,p16,and p53 were detected by quantitative reverse transcription-polymerase chain reaction(qRT-PCR),and histopathological changes were observed by hematoxylin-eosin(HE)and Masson staining.Serum urea nitrogen and creatinine,and catalase(CAT),glutathione peroxidase(GSH-PX),superoxide dismutase(SOD),and malondialdehyde(MDA)levels in the kidney tissues were detected.Reactive oxygen species(ROS)were detected by dihydroethdium staining and protein expression levels of collagen Ⅲ,α-smooth muscle actin(α-SMA),Protein expression of S100A8/9 was detected by immunofluorescence,and transforming growth factor(TGF)-β1 levels in kidney tissues and key factors in the NLRP3/Caspase-1/IL-1β inflammatory pathway were detected by Western Blot.A renal senescence model using HK-2 cells was constructed using H2O2 in vitro,and expression levels of the senescence proteins p21 and p16 and mRNA expression levels of the inflammatory factors IL-18 and tumor necrosis factor-α(TNF-α)were detected.Cell senescence was observed by SA-β-Gal staining.The effects of the S100A8/9 inhibitor paquinimod on expression levels of S100A8/9 and NLRP3/Caspase-1/IL-1β pathway-related proteins in the aging model were also detected.Results mRNA levels of the aging genes p21,p16,and p53 in kidney tissues were significantly increased in rats in the D-galactose group compared with the control group(P<0.01),and SA-β-Gal staining showed a significant increase in senescent cells(P<0.01).Serum blood urea nitrogen and creatinine levels increased(P<0.05),CAT,GSH-PX,and SOD activities decreased(P<0.01),while MDA activity increased in the D-galactose group(P<0.01).Collagen Ⅲ,α-SMA,and TGFβ1 expression and the ROS content in tissues increased(P<0.05).Glomeruli were atrophied or absent in the D-galactose group,the lumens of the renal sacs and renal tubules were enlarged,the nuclei were deeply stained and constricted,and numerous collagen fibers were deposited.Levels of S100A8 and S100A9 protein(P<0.01),as well as NLRP3,Caspase-1,and IL-1β increased(P<0.05).Paquinimod alleviated HK-2 cell senescence and decreased expression levels of the senescence proteins p21 and p16,and mRNA levels of the inflammatory factors IL-18 and TNF-α(P<0.05,P<0.01).The number of senile cells was also decreased,shown by SA-β-Gal staining(P<0.01).Paquinimod also inhibited the protein expression of S100A8 and S100A9(P<0.01)and NLRP3,Caspase-1,and IL-1β(P<0.05 or P<0.01).Conclusions S100A8/9 participates in the chronic inflammatory response by activating the NLRP3/Caspase-1/IL-1β pathway,thereby promoting D-galactose-induced renal aging.
6.Brain functional magnetic resonance imaging changes in patients with depression after acupuncture treatment: ALE meta-analysis
Jing TANG ; Juan YU ; Ziyao GENG ; Yating YUE ; Jingxia ZHANG
Chinese Journal of Behavioral Medicine and Brain Science 2025;34(11):1005-1011
Objective:Through the meta-analysis of activation likelihood estimation (ALE), the functional magnetic resonance imaging research evidence of acupuncture in the treatment of depression was systematically integrated to analyze the effect of acupuncture on brain activity in patients with depression and systematically reveal its potential neural mechanism.Methods:The literatures which used functional magnetic resonance imaging to explore brain image changes of patients with depression after receiving acupuncture intervention were retrieved from CNKI, VIP database, Wanfang database, PubMed, Embase, EBSCO host MEDLINE and Web of Science database. The retrieval time is from the establishment of the database to March 2025. Data processing was performed using the ALE method by the Ginger ALE 2.3 software. The coordinates of the subjects included in the study were extracted, and the brain regions with abnormal spontaneous brain activity in patients with depression were integrated and analyzed.Results:A total of 28 articles involving 770 subjects were included. Depressed patients receiving acupuncture exhibited increased brain activity in the bilateral medial frontal gyrus and bilateral anterior cingulate cortex. Exploratory analysis revealed that depressed patients receiving acupuncture exhibited reduced activation in the left posterior cerebellar lobe (MNI: x, y, z = -6, -68, -21; x, y, z = -9, -81, -9) compared to healthy controls (voxel sum=5 824 mm 3,P<0.05). Conclusion:Patients with depression who received acupuncture intervention exhibited relatively consistent activation in brain regions compared to depressed patients who did not receive acupuncture and healthy individuals. Compared to healthy individuals, patients with depression who received acupuncture showed reduced activation in the left posterior cerebellum. It is suggested that acupuncture intervention may play a therapeutic role by regulating the functional state of emotion-related brain regions and brain networks.
7.Anti-central-fatigue effect of maca via mitochondrial biogenesis via the AMPK/SIRT1/PGC-1α pathway in rats
Wenhuan YAO ; Wen ZHOU ; Yaxuan LI ; Ziyao LI ; Jing ZHANG ; Shibo LYU ; Hui LI
Chinese Journal of Comparative Medicine 2025;35(7):36-43
Objective To examine the anti-central-fatigue function of maca and its underlying mechanism.Methods Forty Sprague-Dawley rats were divided randomly into a negative control group,model control group,and low-,medium-,and high-dose maca groups(0.6,1.2,and 2.4 g/kg·body weight).Rats in all groups except the negative control group were subjected to multi-factor stimulation,including cold-water swimming,sleep deprivation,restraining,and tail-clamping,to establish central fatigue rat models.Rats in the low-,medium-,and high-dose maca groups received 0.6,1.2,or 2.4 g/kg maca,respectively,by gavage for 35 days.Behavioral testing was carried out using the Morris water-maze,sucrose-preference,and tail-suspension tests.Markers of oxidative stress in the hippocampus,including superoxide dismutase(SOD),malondialdehyde(MDA),and catalase(CAT),were detected using test kits.Proteins connected with the AMP-activated protein kinase(AMPK)/sirtuin 1(SIRT1)/peroxisome proliferator-activated receptor γ coactivator 1-α(PGC-1α)signaling pathway in the hippocampus were detected by Western blot.Results Rats in the low-,medium-,and high-dose maca groups spent significantly more time in the target quadrant compared with the model control group(P<0.05 or P<0.01),but there was no significant dose-effect relationship.Rats in the medium-and high-dose maca groups showed decreased escape latency(P<0.05),increased time crossing the platform location(P<0.05),increased sucrose preference(P<0.05),decreased tail suspension time(P<0.05),increased the activities of CAT(P<0.01)and SOD(P<0.05),and decreased MDA content(P<0.01).Rats in the low-,medium-,and high-dose maca groups also showed significantly increased protein expression levels of AMPK and nuclear respiratory factor 1(P<0.01 or P<0.05),but no significant dose-effect relationship was observed.Rats in the medium-and high-dose maca groups showed increased protein expression of PGC-1α(P<0.05 or P<0.01),and rats in the high-dose maca group showed increased protein expression of SIRT1 and mitochondrial transcription factor A(P<0.05 or P<0.01).Conclusions Maca can improve the indicators of central fatigue in rats,determined by behavioral testing and oxidative stress-related factors.The underlying mechanism may be related to its regulatory effects on the AMPK/SIRT1/PGC-1α signaling pathway.
8.Work Memory Impairment in Patients with Chronic Fatigue Syndrome and Spleen Deficiency
Tian ZHOU ; Yunhe ZHANG ; Ziyao WU ; Sitong FENG ; Yanzhe NING ; Hongxiao JIA
World Science and Technology-Modernization of Traditional Chinese Medicine 2025;27(11):3148-3156
Objective To characterize working memory performance in patients with chronic fatigue syndrome(CFS)and spleen-deficiency syndrome and to examine its associations with clinical symptoms by Sternberg working memory task(SWMT).Methods 31 CFS patients meeting both CDC-1994 criteria and consensus criteria for spleen-deficiency pattern were recruited from outpatient clinics and universities from September 2022 and June 2025.31 healthy controls were also recruited based on age,sex,and education.All subjects completed the SWMT.Group differences were analyzed.Within the CFS cohort,reaction time(RT)was correlated with scores on the checklist individual strength(CIS),36-item short-form health survey(sf-36),and fatigue scale-14(FS-14).Mediation was examined.Results RT lengthened with increasing memory load in both groups.CFS patients displayed slower RTs than controls in the baseline and 6-digit set(P<0.05).The 3-digit RT difference,though not significant(P>0.05),yielded a medium effect size(r=0.36).Accuracy did not differ between two groups.Among CFS patients,3-digit RT correlated positively with CIS total and the 4 sub-scale scores.6-digit RT correlated with the SF-36 health-transition dimension(r=0.396,P=0.027).CIS and FS-14 scores directly impaired SF-36 social functioning without working-memory mediating.Conclusion CFS patients with spleen-deficiency exhibit slowed processing speed rather than capacity loss.The close link between working-memory slowing and fatigue suggests a distinct neural basis.These results support the traditional concept"the spleen stores Yi"and integrate TCM pattern differentiation with modern cognitive neuroscience in CFS.
9.Changes in S100A8/9 and NLRP3/Caspase-1/interleukin-1β pathway in kidney-aging rats induced by D-galactose
Dandan FENG ; Ying ZHOU ; Ziyao PANG ; Yueqin CAI ; Chu CHEN ; Jianyou ZHANG ; Dejun WANG
Acta Laboratorium Animalis Scientia Sinica 2025;33(6):823-835
Objective To investigate changes in the pro-inflammatory mediator S100A8/9 and NLRP3/Caspase-1/IL-1β pathway in a rat kidney-aging model induced by D-galactose.Methods Twelve SD rats were divided into control and D-galactose groups,and injected subcutaneously in the back of the neck with D-galactose(150 mg/kg)to establish a rat model of kidney aging.Kidney samples were collected under anesthesia after 8 weeks.Kidneys were stained for senescence-associated beta-galactosidase(SA-β-Gal),mRNA expression levels of the aging-related genes p21,p16,and p53 were detected by quantitative reverse transcription-polymerase chain reaction(qRT-PCR),and histopathological changes were observed by hematoxylin-eosin(HE)and Masson staining.Serum urea nitrogen and creatinine,and catalase(CAT),glutathione peroxidase(GSH-PX),superoxide dismutase(SOD),and malondialdehyde(MDA)levels in the kidney tissues were detected.Reactive oxygen species(ROS)were detected by dihydroethdium staining and protein expression levels of collagen Ⅲ,α-smooth muscle actin(α-SMA),Protein expression of S100A8/9 was detected by immunofluorescence,and transforming growth factor(TGF)-β1 levels in kidney tissues and key factors in the NLRP3/Caspase-1/IL-1β inflammatory pathway were detected by Western Blot.A renal senescence model using HK-2 cells was constructed using H2O2 in vitro,and expression levels of the senescence proteins p21 and p16 and mRNA expression levels of the inflammatory factors IL-18 and tumor necrosis factor-α(TNF-α)were detected.Cell senescence was observed by SA-β-Gal staining.The effects of the S100A8/9 inhibitor paquinimod on expression levels of S100A8/9 and NLRP3/Caspase-1/IL-1β pathway-related proteins in the aging model were also detected.Results mRNA levels of the aging genes p21,p16,and p53 in kidney tissues were significantly increased in rats in the D-galactose group compared with the control group(P<0.01),and SA-β-Gal staining showed a significant increase in senescent cells(P<0.01).Serum blood urea nitrogen and creatinine levels increased(P<0.05),CAT,GSH-PX,and SOD activities decreased(P<0.01),while MDA activity increased in the D-galactose group(P<0.01).Collagen Ⅲ,α-SMA,and TGFβ1 expression and the ROS content in tissues increased(P<0.05).Glomeruli were atrophied or absent in the D-galactose group,the lumens of the renal sacs and renal tubules were enlarged,the nuclei were deeply stained and constricted,and numerous collagen fibers were deposited.Levels of S100A8 and S100A9 protein(P<0.01),as well as NLRP3,Caspase-1,and IL-1β increased(P<0.05).Paquinimod alleviated HK-2 cell senescence and decreased expression levels of the senescence proteins p21 and p16,and mRNA levels of the inflammatory factors IL-18 and TNF-α(P<0.05,P<0.01).The number of senile cells was also decreased,shown by SA-β-Gal staining(P<0.01).Paquinimod also inhibited the protein expression of S100A8 and S100A9(P<0.01)and NLRP3,Caspase-1,and IL-1β(P<0.05 or P<0.01).Conclusions S100A8/9 participates in the chronic inflammatory response by activating the NLRP3/Caspase-1/IL-1β pathway,thereby promoting D-galactose-induced renal aging.
10.Construction and validation of a prognostic model for colon cancer based on anoikis-related genes
Tao ZHANG ; Ziyao LI ; Yingying SUN ; Boyang LI ; Zhao WANG ; Zhifu YANG
Cancer Research and Clinic 2025;37(1):55-63
Objective:To construct and validate a prognostic model of colon cancer based on differentially expressed anoikis-related genes, and to preliminarily investigate the relationship between anoikis-related genes and the tumor immune microenvironment of colon cancer.Methods:A total of 472 cancer tissues samples of patients with colon cancer, RNA sequencing data and clinical data of 41 normal tissues samples were downloaded from the Cancer Genome Atlas (TCGA) database between the establishment time and July in 2024. A total of 919 genes related to anoikis were screened out from GeneCards database, and the common genes were selected from the RNA sequencing gene datasets of colon cancer and normal colon tissues in the TCGA database, among which the differentially expressed anoikis-related genes of colon cancer and normal colon tissues were screened out based on P < 0.05. Furthermore, genes related to the prognosis of 446 colon cancer patients with prognostic data in the TCGA database were screened by using univariate Cox proportional risk model; the genes with P < 0.05 were further screened out and a colon cancer prognosis model was constructed by using LASSO-Cox proportional risk model. The risk score of the above 446 colon cancer patients in the TCGA database was calculated according to the prognostic model, and the patients were divided into high-risk (≥ median value) group and low-risk (< median value) group according to the median risk score, and the overall survival of the 2 groups was analyzed by using the Kaplan-Meier method. The risk score based on R software-based time ROC program package was used to predict 1-year, 2-year, 3-year overall survival therapeutic efficacy of colon cancer patients in the TCGA database. According to the median risk score of colon cancer patients in the TCGA database, the patients in the International Cancer Genome Consortium (ICGC) database were divided into high-risk group and low-risk group. Kaplan-Meier method and receiver operating characteristic (ROC) curve were used to verify the predictive effect of the prognostic model. The differentially expressed genes between low-risk group and high-risk group stratified by prognostic model risk score in the TCGA database were used to perform single sample gene set enrichment analysis (ssGESA) of immune cells and immune function by using R software related programs. The differences in risk scores of patients with different immunophenotypes (including inflammator response type, wound healing type, interferon gamma dominant type and lymphocyte depletion type) were compared; and correlation analysis of infiltration and risk scores between immune cells and stromal cells in tumor microenvironment was made. Based on the tumor immune function and rejection (TIDE) database, the relationship between the prognostic model risk score and programmed death receptor ligand 1 (PD-L1) gene expression level was analyzed. Results:Based on anoikis-related genes in the GeneCards database, 236 differentially expressed anoikis-related genes between colon cancer tissues and normal tissues were obtained from the TCGA database. LASSO Cox regression was applied to establish a prognostic model constructed by 7 differentially expressed anoikis-related genes in cancer tissues and normal colon tissues related to the prognosis of colon cancer. Risk score = 0.366×TIMP1-0.404×NAT1+0.207×LTB4R2+0.075×INHBB+0.140×CD36-0.109×MMP3+2.994×OFCC1. The median risk score of 446 colon cancer patients in the TCGA database was 1.754 719 545. Survival analysis showed that the overall survival of colon cancer patients in high-risk group of the TCGA database was worse than that in low-risk group ( P < 0.001); ROC curve analysis showed that the area under the curve for predicting 1-year, 2-year and 3-year overall survival of patients in the TCGA database based on the prognostic model risk score was 0.705, 0.731 and 0.723, respectively. Kaplan-Meier method analysis showed that in the ICGC database, the overall survival of colon cancer patients in high-risk group was worse than that in low-risk group ( P = 0.041); ROC curve analysis showed that the area under the curve of prognostic model risk score for predicting 1-year and 2-year overall survival of colon cancer patients in the ICGC database was 0.663 and 0.966, respectively. ssGESA analysis showed that macrophage level in high-risk group was higher than that in low-risk group, helper T (Th) 1 cell and Th2 cell levels in high-risk group were lower than those in low-risk group (all P < 0.01). In terms of immune function, the cell killing activity and histocompatibility complex Ⅰ level in high-risk group were lower than those in low-risk group, and type Ⅱ interferon response score in high-risk group was higher than that in low-risk group (all P < 0.05). The analysis of immunophenotype showed that the risk score of inflammatory response type was higher than that of wound healing type ( P < 0.05), and there was no statistically significant difference in risk score between the other 2 types (all P > 0.05). Risk score was positively correlated with stromal cell infiltration score ( R = 0.340, P < 0.001) and immune cell infiltration score ( R = 0.148, P < 0.05); the expression level of PD-L1 in high-risk group was higher than that in low-risk group in the TCGA database ( P = 0.048), and the expression level of PD-L1 was positively correlated with risk score ( R = 0.130, P = 0.009). Conclusions:A prognostic model of colon cancer constructed by anoikis-related genes can better predict the prognosis of colon cancer patients, and anoikis-related genes may play an important role in tumor immunity of colon cancer.

Result Analysis
Print
Save
E-mail