Objective:To analyze the clinicopathological factors affecting the prognosis of patients after curative resection of lung metastases (LMs) from colorectal cancer (CRC) and to construct a clinical risk scoring (CRS) model.Methods:This study retrospectively collected clinicopathological data and follow-up information on 132 patients who underwent radical resection of LMs from CRC at the Cancer Hospital of the Chinese Academy of Medical Sciences between January 2010 and December 2020. We analyzed the clinicopathological factors influencing patient prognosis using univariate and multivariate Cox proportional hazards regression models, and we developed a risk stratification model for prognostic prediction.Results:The median follow-up duration for the cohort of 132 patients was 54.2 months. During this period, 61 patients (46.2%) experienced recurrence or distant metastasis, resulting in a 5-year DFS rate of 54.1%. Additionally, 33 patients (25.0%) died, corresponding to a 5-year overall survival (OS) rate of 76.7%. Univariate Cox proportional hazards regression model analysis indicated that ten clinicopathological factors were significantly associated with OS (all P＜0.05). These factors include the total number of lymph nodes (LNs) dissected from the primary tumor (PT) <16, the number of negative LNs from the PT <13, pN(+) of the PT, logarithmic odds of positive lymph nodes (LODDS) of the PT ≥-1.1, lymph nodes ratio (LNR) of the PT ≥0.02, preoperative carcinoembryonic antigen (CEA) level before LMs resection ≥10 ng/ml, the presence of hilar/mediastinal LN metastasis, the number of LMs ≥2, the maximum diameter of LMs ≥2.5 cm, and the necessity for hilar/mediastinal lymphadenectomy. Multivariate Cox proportional hazards regression analysis identified the number of negative LNs <13 ( HR=3.01, 95% CI: 1.28-7.03, P=0.011), pN(+) of the PT ( HR=5.04, 95% CI: 1.51-16.84, P=0.009), preoperative CEA level before LMs resection ≥10 ng/ml ( HR=5.39, 95% CI: 1.80-16.19, P=0.003), the number of LMs ≥2 ( HR=2.47, 95% CI: 1.09-5.60, P=0.030), and the necessity for hilar/mediastinal lymphadenectomy ( HR=2.74, 95% CI: 1.15-6.52, P=0.023) as independent prognostic risk factors. Patients were categorized based on independent risk factors, revealing statistically significant differences in OS across the groups with CRS scores of ≤2, 3～4, and ≥5 ( P＜0.001). Conclusions:Independent risk factors associated with LMs from CRC patients include the number of negative LNs <13, pN(+) of the PT, preoperative CEA level before LMs resection ≥10 ng/ml, the number of LMs ≥2, and the necessity for hilar/mediastinal lymphadenectomy. Patients scoring 3 or higher on the CRS model may warrant cautious assessment for the appropriateness of direct surgical treatment.

Objective:To construct and validate a prognostic prediction model for patients with ovarian metastases from colorectal cancer after radical resection.Methods:A retrospective case series analysis was conducted on the clinical and pathological data of 81 patients with colorectal cancer and ovarian metastases who underwent radical resection for ovarian metastases at the Department of Colorectal Surgery, Cancer Hospital, Chinese Academy of Medical Sciences, between January 2014 and December 2023. The patients were all female, with an age ( M(IQR)) of 49(13) years (range: 22 to 79 years). The primary tumor was located in the colon in 60 cases (74.1%) and in the rectum in 21 cases (25.9%). Univariate and multivariate Cox regression analyses were used to identify independent risk factors affecting prognosis. A risk scoring system was constructed, and patients were assigned to high-risk and low-risk groups based on their risk scores. The predictive performance of the scoring system was assessed, and 5-fold cross-validation was performed to evaluate the model′s stability on the internal dataset. Results:Among the 81 patients with ovarian metastases, a high proportion had T4 stage (58 cases, 71.6%), lymph node positivity (68 cases, 84.0%), and colon cancer (60 cases, 74.1%). Preoperative imaging suggested unilateral ovarian metastasis in 15 patients (23.4%), but pathological examination after bilateral oophorectomy confirmed bilateral ovarian metastases. Among the 17 patients who initially underwent unilateral oophorectomy, 11 developed contralateral ovarian metastases at varying times postoperatively. Univariate Cox proportional hazards regression analysis revealed that positive lymph node ratio ( HR=2.68,95% CI:1.41 to 5.09, P=0.003), N stage ( HR=2.07,95% CI:1.08 to 3.95, P=0.028),maximum diameter of metastatic tumors ( HR=2.27,95% CI:1.04 to 4.96, P=0.040),and peritoneal metastasis or ascites at the time of ovarian metastasis ( HR=2.04,95% CI:1.02 to 4.08, P=0.043) were significantly associated with overall survival in patients with ovarian metastasis from colorectal cancer. Multivariate regression analysis identified that positive lymph node ratio ( HR=3.34,95% CI:1.08 to 10.34, P=0.037) and maximum diameter of metastatic tumors ( HR=2.65,95% CI:1.19 to 5.88, P=0.017) were independent prognostic factors for overall survival following radical oophorectomy in patients with ovarian metastasis from colorectal cancer. Based on the regression coefficients from the multivariate analysis for variables (ovarian metastatic tumor diameter ≥6 cm, positive lymph node ratio ≥0.3,and presence of peritoneal metastasis or ascites), a risk scoring system was developed. Using the optimal cutoff value (154 points) for the risk score,patients were divided into high-risk (19 cases) and low-risk (62 cases) groups. Kaplan-Meier survival curves demonstrated that the high-risk group had significantly lower median overall survival (27 months) and median disease-free survival (22 months) compared to the low-risk group (median overall survival 90 months,median disease-free survival not reached; both P<0.01). Receiver operating characteristic curve analysis showed that the area under the curve(AUC) for predicting 1-,3-,and 5-year overall survival was 0.731(95% CI:0.563 to 0.899), 0.703(95% CI:0.573 to 0.833), and 0.776(95% CI: 0.657 to 0.894), respectively. The AUC for predicting 1-,3-, and 5-year disease-free survival was 0.724(95% CI:0.397 to 0.993),0.710(95% CI:0.514 to 0.906),and 0.688(95% CI:0.478 to 0.898),respectively,indicating good performance of the model.The decision curve analysis showed that the model has good clinical net benefit and the results of the 5-fold cross-validation showed that the model demonstrated stability in the internal dataset. Conclusions:When performing radical resection for ovarian metastasis from colorectal cancer,bilateral oophorectomy should be considered to minimize the risk of postoperative recurrence. Patients with ovarian metastasis from colorectal cancer,characterized by a metastatic tumor diameter ≥6 cm,a positive lymph node ratio ≥0.3,and the presence of peritoneal metastasis or ascites, tend to have a poorer prognosis. Based on these findings,a clinical prognostic scoring system for radical resection of ovarian metastasis from colorectal cancer has been developed to stratify patients into different risk groups and may assist in postoperative risk assessment and management.

Objective:To explore the clinicopathological features of rectal neuroendocrine tumor (R-NET) G2, identify prognostic factors, and summarize treatment experience.Methods:The clinical data of patients diagnosed with R-NET G2 by pathological diagnosis admitted to Cancer Hospital of the Chinese Academy of Medical Sciences from January 2003 to September 2023 were retrospectively analyzed. The Fisher's exact test and Kaplan-Meier curves were performed to analyze the association between pathological features and prognosis.Results:A total of 22 patients were enrolled in this study and 21 patients were followed up for a period of 6-98 months with a median follow-up time of 42 months. 5 patients died due to tumor progression during the follow-up period. The 1-, 3-, and 5-year cancer-specific survival (CSS) of the whole group were 100.0%, 92.9%, and 69.6%, respectively. Of the 22 patients, 20 underwent surgical treatment, of which 15 underwent postoperative adjuvant therapy; 2 underwent medical treatment for liver and bone multiple metastases. The 5-year survival rates of patients with tumours ≥2 cm in length, T2-3 stage, lymph node metastasis, and distant metastasis (57.1%, 68.8%, 66.7%, and 63.6%, respectively) were shorter than those of patients with tumours <2 cm in length, T1 stage, no lymph node metastasis, and no distant metastasis (all 100.0%, P＜0.001). In addition, patients with liver metastases had larger primary tumor diameters and higher T-stages compared with those without distant metastasis ( P＜0.05). Conclusions:R-NET G2 has a high degree of malignancy compared with G1 and a high propensity for metastasis. Clinicians should formulate appropriate diagnostic and treatment strategies based on factors such as tumor size, depth of invasion, lymph node status, presence of distant metastasis, and the location and extent of distant metastasis.

Objective:To analyze the clinicopathological factors affecting the prognosis of patients after curative resection of lung metastases (LMs) from colorectal cancer (CRC) and to construct a clinical risk scoring (CRS) model.Methods:This study retrospectively collected clinicopathological data and follow-up information on 132 patients who underwent radical resection of LMs from CRC at the Cancer Hospital of the Chinese Academy of Medical Sciences between January 2010 and December 2020. We analyzed the clinicopathological factors influencing patient prognosis using univariate and multivariate Cox proportional hazards regression models, and we developed a risk stratification model for prognostic prediction.Results:The median follow-up duration for the cohort of 132 patients was 54.2 months. During this period, 61 patients (46.2%) experienced recurrence or distant metastasis, resulting in a 5-year DFS rate of 54.1%. Additionally, 33 patients (25.0%) died, corresponding to a 5-year overall survival (OS) rate of 76.7%. Univariate Cox proportional hazards regression model analysis indicated that ten clinicopathological factors were significantly associated with OS (all P＜0.05). These factors include the total number of lymph nodes (LNs) dissected from the primary tumor (PT) <16, the number of negative LNs from the PT <13, pN(+) of the PT, logarithmic odds of positive lymph nodes (LODDS) of the PT ≥-1.1, lymph nodes ratio (LNR) of the PT ≥0.02, preoperative carcinoembryonic antigen (CEA) level before LMs resection ≥10 ng/ml, the presence of hilar/mediastinal LN metastasis, the number of LMs ≥2, the maximum diameter of LMs ≥2.5 cm, and the necessity for hilar/mediastinal lymphadenectomy. Multivariate Cox proportional hazards regression analysis identified the number of negative LNs <13 ( HR=3.01, 95% CI: 1.28-7.03, P=0.011), pN(+) of the PT ( HR=5.04, 95% CI: 1.51-16.84, P=0.009), preoperative CEA level before LMs resection ≥10 ng/ml ( HR=5.39, 95% CI: 1.80-16.19, P=0.003), the number of LMs ≥2 ( HR=2.47, 95% CI: 1.09-5.60, P=0.030), and the necessity for hilar/mediastinal lymphadenectomy ( HR=2.74, 95% CI: 1.15-6.52, P=0.023) as independent prognostic risk factors. Patients were categorized based on independent risk factors, revealing statistically significant differences in OS across the groups with CRS scores of ≤2, 3～4, and ≥5 ( P＜0.001). Conclusions:Independent risk factors associated with LMs from CRC patients include the number of negative LNs <13, pN(+) of the PT, preoperative CEA level before LMs resection ≥10 ng/ml, the number of LMs ≥2, and the necessity for hilar/mediastinal lymphadenectomy. Patients scoring 3 or higher on the CRS model may warrant cautious assessment for the appropriateness of direct surgical treatment.

Objective:To construct and validate a prognostic prediction model for patients with ovarian metastases from colorectal cancer after radical resection.Methods:A retrospective case series analysis was conducted on the clinical and pathological data of 81 patients with colorectal cancer and ovarian metastases who underwent radical resection for ovarian metastases at the Department of Colorectal Surgery, Cancer Hospital, Chinese Academy of Medical Sciences, between January 2014 and December 2023. The patients were all female, with an age ( M(IQR)) of 49(13) years (range: 22 to 79 years). The primary tumor was located in the colon in 60 cases (74.1%) and in the rectum in 21 cases (25.9%). Univariate and multivariate Cox regression analyses were used to identify independent risk factors affecting prognosis. A risk scoring system was constructed, and patients were assigned to high-risk and low-risk groups based on their risk scores. The predictive performance of the scoring system was assessed, and 5-fold cross-validation was performed to evaluate the model′s stability on the internal dataset. Results:Among the 81 patients with ovarian metastases, a high proportion had T4 stage (58 cases, 71.6%), lymph node positivity (68 cases, 84.0%), and colon cancer (60 cases, 74.1%). Preoperative imaging suggested unilateral ovarian metastasis in 15 patients (23.4%), but pathological examination after bilateral oophorectomy confirmed bilateral ovarian metastases. Among the 17 patients who initially underwent unilateral oophorectomy, 11 developed contralateral ovarian metastases at varying times postoperatively. Univariate Cox proportional hazards regression analysis revealed that positive lymph node ratio ( HR=2.68,95% CI:1.41 to 5.09, P=0.003), N stage ( HR=2.07,95% CI:1.08 to 3.95, P=0.028),maximum diameter of metastatic tumors ( HR=2.27,95% CI:1.04 to 4.96, P=0.040),and peritoneal metastasis or ascites at the time of ovarian metastasis ( HR=2.04,95% CI:1.02 to 4.08, P=0.043) were significantly associated with overall survival in patients with ovarian metastasis from colorectal cancer. Multivariate regression analysis identified that positive lymph node ratio ( HR=3.34,95% CI:1.08 to 10.34, P=0.037) and maximum diameter of metastatic tumors ( HR=2.65,95% CI:1.19 to 5.88, P=0.017) were independent prognostic factors for overall survival following radical oophorectomy in patients with ovarian metastasis from colorectal cancer. Based on the regression coefficients from the multivariate analysis for variables (ovarian metastatic tumor diameter ≥6 cm, positive lymph node ratio ≥0.3,and presence of peritoneal metastasis or ascites), a risk scoring system was developed. Using the optimal cutoff value (154 points) for the risk score,patients were divided into high-risk (19 cases) and low-risk (62 cases) groups. Kaplan-Meier survival curves demonstrated that the high-risk group had significantly lower median overall survival (27 months) and median disease-free survival (22 months) compared to the low-risk group (median overall survival 90 months,median disease-free survival not reached; both P<0.01). Receiver operating characteristic curve analysis showed that the area under the curve(AUC) for predicting 1-,3-,and 5-year overall survival was 0.731(95% CI:0.563 to 0.899), 0.703(95% CI:0.573 to 0.833), and 0.776(95% CI: 0.657 to 0.894), respectively. The AUC for predicting 1-,3-, and 5-year disease-free survival was 0.724(95% CI:0.397 to 0.993),0.710(95% CI:0.514 to 0.906),and 0.688(95% CI:0.478 to 0.898),respectively,indicating good performance of the model.The decision curve analysis showed that the model has good clinical net benefit and the results of the 5-fold cross-validation showed that the model demonstrated stability in the internal dataset. Conclusions:When performing radical resection for ovarian metastasis from colorectal cancer,bilateral oophorectomy should be considered to minimize the risk of postoperative recurrence. Patients with ovarian metastasis from colorectal cancer,characterized by a metastatic tumor diameter ≥6 cm,a positive lymph node ratio ≥0.3,and the presence of peritoneal metastasis or ascites, tend to have a poorer prognosis. Based on these findings,a clinical prognostic scoring system for radical resection of ovarian metastasis from colorectal cancer has been developed to stratify patients into different risk groups and may assist in postoperative risk assessment and management.

Objective:To explore the clinicopathological features of rectal neuroendocrine tumor (R-NET) G2, identify prognostic factors, and summarize treatment experience.Methods:The clinical data of patients diagnosed with R-NET G2 by pathological diagnosis admitted to Cancer Hospital of the Chinese Academy of Medical Sciences from January 2003 to September 2023 were retrospectively analyzed. The Fisher's exact test and Kaplan-Meier curves were performed to analyze the association between pathological features and prognosis.Results:A total of 22 patients were enrolled in this study and 21 patients were followed up for a period of 6-98 months with a median follow-up time of 42 months. 5 patients died due to tumor progression during the follow-up period. The 1-, 3-, and 5-year cancer-specific survival (CSS) of the whole group were 100.0%, 92.9%, and 69.6%, respectively. Of the 22 patients, 20 underwent surgical treatment, of which 15 underwent postoperative adjuvant therapy; 2 underwent medical treatment for liver and bone multiple metastases. The 5-year survival rates of patients with tumours ≥2 cm in length, T2-3 stage, lymph node metastasis, and distant metastasis (57.1%, 68.8%, 66.7%, and 63.6%, respectively) were shorter than those of patients with tumours <2 cm in length, T1 stage, no lymph node metastasis, and no distant metastasis (all 100.0%, P＜0.001). In addition, patients with liver metastases had larger primary tumor diameters and higher T-stages compared with those without distant metastasis ( P＜0.05). Conclusions:R-NET G2 has a high degree of malignancy compared with G1 and a high propensity for metastasis. Clinicians should formulate appropriate diagnostic and treatment strategies based on factors such as tumor size, depth of invasion, lymph node status, presence of distant metastasis, and the location and extent of distant metastasis.

Objective: To explore mental health status of vocational college students and its association with physical exercise and exercise motivation, so as to provide data support for the promotion of physical exercise, exercise motivation theory and the protection of mental health among vocational college students. Methods: From June 2 to July 2, 2023, a convenience sampling method was used to select 1 763 college students from three vocational colleges in Sichuan, Chongqing, and Guizhou for the survey. Physical Activity Rating Scale (PARS-3), Motives for Physical Activities Measure-Revised (MPAM-R), Patient Health Questionnaire-9 (PHQ-9) and Flourishing Scale (FS) were used to assess physical exercise, exercise motivation, depressive symptoms and flourishing levels of vocational college students. Multiple linear regression analysis was conducted to examine the relationship between physical exercise, exercise motivation, depression as the negative indicator of mental health and flourishing as the positive indicator of mental health. Results: There were negative correlations between physical exercise, exercise motivation and their five factors with depression scores among vocational college students ( r=-0.162--0.133, P <0.01). Physical exercise, health motivation and pleasure motivation were positively correlated with flourishing and its eight factors among vocational college students ( r= 0.054 -0.099, P <0.05). Depression scores and flourishing scores varied significantly across students with different levels of physical exercise and exercise motivation ( F=11.18-15.69, 2.80-16.05, P <0.05). Multiple linear regression analysis indicated that physical exercise, health motivation, and total exercise motivation scores were negatively correlated with depressive symptoms ( β = -0.02 , -0.11, -0.35), and physical exercise, health motivation, pleasure motivation, and total exercise motivation scores were positively correlated with flourishing ( β =0.31, 0.32, 0.50, 0.30) ( P <0.05). Conclusions Physical exercise and exercise motivation are significantly associated with mental health of vocational college students. Enhancing physical exercise and increasing exercise motivation in aspects of health and pleasure can promote mental health of vocational college students.

Oral squamous cell carcinoma (OSCC) escape from the immune system is mediated through several immunosuppressive phenotypes that are critical to the initiation and progression of tumors. As a hallmark of cancer, DNA damage repair is closely related to changes in the immunophenotypes of tumor cells. Although flap endonuclease-1 (FEN1), a pivotal DNA-related enzyme is involved in DNA base excision repair to maintain the stability of the cell genome, the correlation between FEN1 and tumor immunity has been unexplored. In the current study, by analyzing the clinicopathological characteristics of FEN1, we demonstrated that FEN1 overexpressed and that an inhibitory immune microenvironment was established in OSCC. In addition, we found that downregulating FEN1 inhibited the growth of OSCC tumors. In vitro studies provided evidence that FEN1 knockdown inhibited the biological behaviors of OSCC and caused DNA damage. Performing multiplex immunohistochemistry (mIHC), we directly observed that the acquisition of critical immunosuppressive phenotypes was correlated with the expression of FEN1. More importantly, FEN1 directly or indirectly regulated two typical immunosuppressive phenotype-related proteins human leukocyte antigen (HLA-DR) and programmed death receptor ligand 1 (PD-L1), through the interferon-gamma (IFN-γ)/janus kinase (JAK)/signal transducer and activator transcription 1 (STAT1) pathway. Our study highlights a new perspective on FEN1 action for the first time, providing theoretical evidence that it may be a potential immunotherapy target for OSCC.
Humans ; Carcinoma, Squamous Cell/pathology* ; DNA ; Down-Regulation ; Flap Endonucleases/metabolism* ; Head and Neck Neoplasms ; Interferon-gamma/metabolism* ; Mouth Neoplasms/pathology* ; Phenotype ; Squamous Cell Carcinoma of Head and Neck ; Tumor Microenvironment ; Janus Kinases/metabolism*

OBJECTIVE: To examine the anti-inflammatory effects and potential mechanisms of polypeptide from Moschus (PPM) in lipopolysaccharide (LPS)-induced THP-1 macrophages and BALB/c mice. METHODS: The polypeptide was extracted from Moschus and analyzed by high-performance liquid chromatography and sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). Subsequently, LPS was used to induce inflammation in THP-1 macrophages and BALB/c mice. In LPS-treated or untreated THP-1 macrophages, cell viability was observed by cell counting kit 8 and lactate dehydrogenase release assays; the proinflammatory cytokines and reactive oxygen species (ROS) were measured by enzyme-linked immunosorbent assay and flow cytometry, respectively; and protein and mRNA levels were measured by Western blot and real-time quantitative polymerase chain reaction (qRT-PCR), respectively. In LPS-induced BALB/c mice, the proinflammatory cytokines were measured, and lung histology and cytokines were observed by hematoxylin and eosin (HE) and immunohistochemical (IHC) staining, respectively. RESULTS: The SDS-PAGE results suggested that the molecular weight of purified PPM was in the range of 10-26 kD. In vitro, PPM reduced the production of interleukin 1β (IL-1β), IL-18, tumor necrosis factor α (TNF-α), IL-6 and ROS in LPS-induced THP-1 macrophages (P<0.01). Western blot analysis demonstrated that PPM inhibited LPS-induced nuclear factor κB (NF-κB) pathway and thioredoxin interacting protein (TXNIP)/nucleotide-binding oligomerization domain, leucine-rich repeat and pyrin domain containing 3 (NLRP3) inflammasome pathway by reducing protein expression of phospho-NF-κB p65, phospho-inhibitors of NF-κB (Iκ Bs) kinase α/β (IKKα/β), TXNIP, NLRP3, apoptosis-associated speck-like protein containing a caspase recruitment domain (ASC), and pro-caspase-1 (P<0.05 or P<0.01). In addition, qRT-PCR revealed the inhibitory effects of PPM on the mRNA levels of TXNIP, NLRP3, ASC, and caspase-1 (P<0.05 or P<0.01). Furthermore, in LPS-induced BALB/c mice, PPM reduced TNF-α and IL-6 levels in serum (P<0.05 or P<0.01), decreased IL-1β and IL-18 levels in the lungs (P<0.01) and alleviated pathological injury to the lungs. CONCLUSION PPM could attenuate LPS-induced inflammation by inhibiting the NF-κB-ROS/NLRP3 pathway, and may be a novel potential candidate drug for treating inflammation and inflammation-related diseases.

Background::With functionally heterogeneous cells, tumors comprise a complex ecosystem to promote tumor adaptability and evolution under strong selective pressure from the given microenvironment. Diversifying tumor cells or intra-tumor heterogeneity is essential for tumor growth, invasion, and immune evasion. However, no reliable method to classify tumor cell subtypes is yet available. In this study, we introduced the single-cell sequencing combined with copy number characteristics to identify the types of tumor cells in microsatellite stable (MSS) colorectal cancer (CRC).Methods::To characterize the somatic copy number alteration (SCNA) of MSS CRC in a single cell profile, we analyzed 26 tissue samples from 19 Korean patients (GSE132465, the Samsung Medical Center [SMC] dataset) and then verified our findings with 15 tissue samples from five Belgian patients (GSE144735, the Katholieke Universiteit Leuven 3 [KUL3] dataset). The Cancer Genome Atlas (TCGA) cohort, GSE39582 cohort, and National Cancer Center (NCC) cohort (24 MSS CRC patients were enrolled in this study between March 2017 and October 2017) were used to validate the clinical features of prognostic signatures.Results::We employed single cell RNA-sequencing data to identify three types of tumor cells in MSS CRC by their SCNA characteristics. Among these three types of tumor cells, C1 and C3 had a higher SCNA burden; C1 had significant chromosome 13 and 20 amplification, whereas C3 was the polar opposite of C1, which exhibited deletion in chromosome 13 and 20. The three types of tumor cells exhibited various functions in the tumor microenvironment and harbored different mutations. C1 and C2 were linked to the immune response and hypoxia, respectively, while C3 was critical for cell adhesion activity and tumor angiogenesis. Additionally, one gene ( OLFM4) was identified as epithelium-specific biomarker of better prognosis of CRC (TCGA cohort: P = 0.0110; GSE39582 cohort: P= 0.0098; NCC cohort: P= 0.0360). Conclusions::On the basis of copy number characteristics, we illustrated tumor heterogeneity in MSS CRC and identified three types of tumor cells with distinct roles in tumor microenvironment. By understanding heterogeneity in the intricate tumor microenvironment, we gained an insight into the mechanisms of tumor evolution, which may support the development of therapeutic strategies.