Single-cell transcriptomic profiling unravels the adenoma-initiation role of protein tyrosine kinases during colorectal tumorigenesis

Systematic Analysis of Human Colorectal Cancer scRNA-seq Revealed Limited Pro-tumoral IL-17 Production Potential in Gamma Delta T Cells

Gamma delta (γδ) T cells play a crucial role in anti-tumor immunity due to their cytotoxic properties. However, the role and extent of γδ T cells in production of pro-tumorigenic interleukin- 17 (IL-17) within the tumor microenvironment (TME) of colorectal cancer (CRC) remains controversial. In this study, we re-analyzed nine published human CRC whole-tissue single-cell RNA sequencing (scRNA-seq) datasets, identifying 18,483 γδ T cells out of 951,785 total cells, in the neoplastic or adjacent normal tissue of 165 human CRC patients. Our results confirm that tumor-infiltrating γδ T cells exhibit high cytotoxicity-related transcription in both tumor and adjacent normal tissues, but critically, none of the γδ T cell clusters showed IL-17 production potential. We also identified various γδ T cell subsets, including Teff, TRM, Tpex, and Tex, and noted an increased expression of cytotoxic molecules in tumor-infiltrating γδ T cells compared to their normal area counterparts. Our work demonstrates that γδ T cells in CRC primarily function as cytotoxic effector cells rather than IL-17 producers, mitigating the concerns about their potential pro-tumorigenic roles in CRC, highlighting the importance of accurately characterizing these cells for cancer immunotherapy research and the unneglectable cross-species discrepancy between the mouse and human immune system in the study of cancer immunology.

Using a pan-cancer atlas to investigate tumour associated macrophages as regulators of immunotherapy response

The paradigm for macrophage characterization has evolved from the simple M1/M2 dichotomy to a more complex model that encompasses the broad spectrum of macrophage phenotypic diversity, due to differences in ontogeny and/or local stimuli. We currently lack an in-depth pan-cancer single cell RNA-seq (scRNAseq) atlas of tumour-associated macrophages (TAMs) that fully captures this complexity. In addition, an increased understanding of macrophage diversity could help to explain the variable responses of cancer patients to immunotherapy. Our atlas includes well established macrophage subsets as well as a number of additional ones. We associate macrophage composition with tumour phenotype and show macrophage subsets can vary between primary and metastatic tumours growing in sites like the liver. We also examine macrophage-T cell functional cross talk and identify two subsets of TAMs associated with T cell activation. Analysis of TAM signatures in a large cohort of immune checkpoint inhibitor-treated patients (CPI1000 + ) identify multiple TAM subsets associated with response, including the presence of a subset of TAMs that upregulate collagen-related genes. Finally, we demonstrate the utility of our data as a resource and reference atlas for mapping of novel macrophage datasets using projection. Overall, these advances represent an important step in both macrophage classification and overcoming resistance to immunotherapies in cancer.

Recent advances in hematopoietic cell kinase in cancer progression: Mechanisms and inhibitors

Hematopoietic cell kinase (Hck), a non-receptor tyrosine kinase belonging to the Src kinase family, is intricately linked to the pathogenesis of numerous human diseases, with a particularly pronounced association with cancer. Hck not only directly impacts the proliferation, migration, and apoptosis of cancer cells but also interacts with JAK/STAT, MEK/ERK, PI3K/AKT, CXCL12/CXCR4, and other pathways. Hck also influences the tumor microenvironment to facilitate the onset and progression of cancer. This paper delves into the functional role and regulatory mechanisms of Hck in various solid tumors. Additionally, it explores the implications of Hck in hematological malignancies. The review culminates with a summary of the current research status of Hck inhibitors, the majority of which are in the pre-clinical phase of investigation. Notably, these inhibitors are predominantly utilized in the therapeutic management of leukemia, with their combinatorial potential indicating promising avenues for future research. In conclusion, this review underscores the significance of the mechanism of Hck in solid tumors. This insight is crucial for comprehending the current research trends regarding Hck: targeted therapy against Hck shows great promise in both diagnosis and treatment of malignant tumors. Further investigation into the role of Hck in cancer, coupled with the development of specific inhibitors, has the potential to revolutionize approaches to cancer treatment.

Tumor suppressive role of the epigenetic master regulator BRD3 in colorectal cancer

Bromodomain and extraterminal domain (BET) family proteins are epigenetic master regulators of gene expression via recognition of acetylated histones and recruitment of transcription factors and co-activators to chromatin. Hence, BET family proteins have emerged as promising therapeutic targets in cancer. In this study, we examined the functional role of bromodomain containing 3 (BRD3), a BET family protein, in colorectal cancer (CRC). In vitro and vivo analyses using BRD3-knockdown or BRD3-overexpressing CRC cells showed that BRD3 suppressed tumor growth and cell cycle G1/S transition and induced p21 expression. Clinical analysis of CRC datasets from our hospital or The Cancer Genome Atlas revealed that BET family genes, including BRD3, were overexpressed in tumor tissues. In immunohistochemical analyses, BRD3 was observed mainly in the nucleus of CRC cells. According to single-cell RNA sequencing in untreated CRC tissues, BRD3 was highly expressed in malignant epithelial cells, and cell cycle checkpoint-related pathways were enriched in the epithelial cells with high BRD3 expression. Spatial transcriptomic and single-cell RNA sequencing analyses of CRC tissues showed that BRD3 expression was positively associated with high p21 expression. Furthermore, overexpression of BRD3 combined with knockdown of, a driver gene in the BRD family, showed strong inhibition of CRC cells in vitro. In conclusion, we demonstrated a novel tumor suppressive role of BRD3 that inhibits tumor growth by cell cycle inhibition in part via induction of p21 expression. BRD3 activation might be a novel therapeutic approach for CRC.

Spatial and single-cell colocalisation analysis reveals MDK-mediated immunosuppressive environment with regulatory T cells in colorectal carcinogenesis

BACKGROUND

Cell-cell interaction factors that facilitate the progression of adenoma to sporadic colorectal cancer (CRC) remain unclear, thereby hindering patient survival.

METHODS

We performed spatial transcriptomics on five early CRC cases, which included adenoma and carcinoma, and one advanced CRC. To elucidate cell-cell interactions within the tumour microenvironment (TME), we investigated the colocalisation network at single-cell resolution using a deep generative model for colocalisation analysis, combined with a single-cell transcriptome, and assessed the clinical significance in CRC patients.

FINDINGS

CRC cells colocalised with regulatory T cells (Tregs) at the adenoma-carcinoma interface. At early-stage carcinogenesis, cell-cell interaction inference between colocalised adenoma and cancer epithelial cells and Tregs based on the spatial distribution of single cells highlighted midkine (MDK) as a prominent signalling molecule sent from tumour epithelial cells to Tregs. Interaction between MDK-high CRC cells and SPP1+ macrophages and stromal cells proved to be the mechanism underlying immunosuppression in the TME. Additionally, we identified syndecan4 (SDC4) as a receptor for MDK associated with Treg colocalisation. Finally, clinical analysis using CRC datasets indicated that increased MDK/SDC4 levels correlated with poor overall survival in CRC patients.

INTERPRETATION

MDK is involved in the immune tolerance shown by Tregs to tumour growth. MDK-mediated formation of the TME could be a potential target for early diagnosis and treatment of CRC.

FUNDING

Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Science Research; OITA Cancer Research Foundation; AMED under Grant Number; Japan Science and Technology Agency (JST); Takeda Science Foundation; The Princess Takamatsu Cancer Research Fund.

Multiple-matrix metabolomics analysis for the distinct detection of colorectal cancer and adenoma

OBJECTIVES

Although colorectal cancer (CRC) is the leading cause of cancer-related morbidity and mortality, current diagnostic tests for early-stage CRC and colorectal adenoma (CRA) are suboptimal. Therefore, there is an urgent need to explore less invasive screening procedures for CRC and CRA diagnosis.

METHODS

Untargeted gas chromatography-mass spectrometry (GC-MS) metabolic profiling approach was applied to identify candidate metabolites. We performed metabolomics profiling on plasma samples from 412 subjects including 200 CRC patients, 160 CRA patients and 52 normal controls (NC). Among these patients, 45 CRC patients, 152 CRA patients and 50 normal controls had their fecal samples tested simultaneously.

RESULTS

Differential metabolites were screened in the adenoma-carcinoma sequence. Three diagnostic models were further developed to identify cancer group, cancer stage, and cancer microsatellite status using those significant metabolites. The three-metabolite-only classifiers used to distinguish the cancer group always keeps the area under the receiver operating characteristic curve (AUC) greater than 0.7. The AUC performance of the classifiers applied to discriminate CRC stage is generally greater than 0.8, and the classifiers used to distinguish microsatellite status of CRC is greater than 0.9.

CONCLUSION

This finding highlights potential early-driver metabolites in CRA and early-stage CRC. We also find potential metabolic markers for discriminating the microsatellite state of CRC. Our study and diagnostic model have potential applications for non-invasive CRC and CRA detection.

Deciphering the molecular and clinical characteristics of TREM2, HCST, and TYROBP in cancer immunity: A comprehensive pan-cancer study

Background

Hematopoietic cell signal transducer (HCST) and tyrosine kinase-binding protein (TYROBP) are triggering receptors expressed on myeloid cells 2 (TREM2), which are pivotal in the immune response to disease. Despite growing evidence underscoring the significance of TREM2, HCST, and TYROBP in certain forms of tumorigenesis, a comprehensive pan-cancer analysis of these proteins is lacking.

Methods

Multiple databases were synthesized to investigate the relationship between TREM2, HCST, TYROBP, and various cancer types. These include prognosis, methylation, regulation by long non-coding RNAs and transcription factors, immune signatures, pathway activity, microsatellite instability (MSI), tumor mutational burden (TMB), single-cell transcriptome profiling, and drug sensitivity.

Results

TREM2, HCST, and TYROBP displayed extensive somatic changes across numerous tumors, and their mRNA expression and methylation levels influenced patient outcomes across multiple cancer types. long non-coding RNA (lncRNA) -messenger RNA (mRNA) and TF-mRNA regulatory networks involving TREM2, HCST, and TYROBP were identified, with lncRNA MEG3 and the transcription factor SIP1 emerging as potential key regulators. Further immune analyses indicated that TREM2, HCST, and TYROBP play critical roles in immune-related pathways and macrophage differentiation, and may be significantly associated with TGF-β and SMAD9. Furthermore, the expression of TREM2, HCST, and TYROBP correlated with the immunotherapy markers TMB and MSI, and influenced sensitivity to immune-targeted drugs, thereby indicating their potential as predictors of immunotherapy outcomes.

Conclusion

This study offers valuable insights into the roles of TREM2, HCST, and TYROBP in tumor immunotherapy, suggesting their potential as prognostic markers and therapeutic targets for various cancers.

Construction of an immune predictive model and identification of TRIP6 as a prognostic marker and therapeutic target of CRC by integration of single-cell and bulk RNA-seq data

BACKGROUND

Investigations elucidating the complex immunological mechanisms involved in colorectal cancer (CRC) and accurately predicting patient outcomes via bulk RNA-Seq analysis have been notably limited. This study aimed to identify the immune status of CRC patients, construct a prognostic model, and identify prognostic signatures via bulk RNA sequencing (RNA-seq) and single-cell RNA-seq (scRNA-seq).

METHODS

The scRNA-seq data of CRC were downloaded from Gene Expression Omnibus (GEO). The UCSC Xena database was used to obtain bulk RNA-seq data. Differentially expressed gene (DEG), functional enrichment, and random forest analyses were conducted in order to identify core genes associated with colorectal cancer (CRC) that were relevant to prognosis. A molecular immune prediction model was developed using logistic regression after screening features using the least absolute shrinkage and selection operator (LASSO). The differences in immune cell infiltration, mutation, chemotherapeutic drug sensitivity, cellular senescence, and communication between patients who were at high and low risk of CRC according to the predictive model were investigated. The prognostic genes that were closely associated with CRC were identified by random survival forest (RSF) analysis. The expression levels and clinical significance of the hub genes were analyzed in vitro. The LoVo cell line was employed to ascertain the biological role of thyroid hormone receptor-interacting protein 6 (TRIP6).

RESULTS

A total of seven main cell subtypes were identified by scRNA-seq analysis. A molecular immune predictive model was constructed based on the risk scores. The risk score was significantly associated with OS, stage, mutation burden, immune cell infiltration, response to immunotherapy, key pathways, and cell-cell communication. The functions of the six hub genes were determined and further utilized to establish a regulatory network. Our findings unequivocally confirmed that TRIP6 upregulation was verified in the CRC samples. After knocking down TRIP6, cell proliferation, migration, and invasion of LoVo cells were inhibited, and apoptosis was promoted.

CONCLUSIONS

The molecular predictive model reliably distinguished the immune status of CRC patients. We further revealed that TRIP6 may act as an oncogene in CRC, making it a promising candidate for targeted therapy and as a prognostic marker for CRC.

Increased circulating regulatory T cells and decreased follicular T helper cells are associated with colorectal carcinogenesis

Objective

Colorectal cancer (CRC) is the third most prevalent cancer worldwide and is associated with high morbidity and mortality rates. Colorectal carcinogenesis occurs via the conventional adenoma-to-carcinoma and serrated pathways. Conventional T helper (Th) and innate lymphoid cells (ILCs) play vital roles in maintaining intestinal homeostasis. However, the contribution of these two major lymphoid cell populations and their associated cytokines to CRC development is unclear. Therefore, we aimed to analyze peripheral lymphocyte profiles during colorectal carcinogenesis.

Methods

We collected 86 blood samples concurrently, and pathologists confirmed the presence of various pathological conditions (i.e., HPs, adenoma, and carcinoma) using hematoxylin and eosin staining. Ten healthy donors were recruited as healthy controls (HCs) from the physical examination center. We performed flow cytometry on peripheral blood mononuclear cells collected from patients with various pathological conditions and the HCs, and cytokines (interleukin-2, interleukin-4, interleukin-5, interleukin-13, interleukin-17A, interleukin-17F, interleukin-22, interferon-γ, and tumor necrosis factor-α) were quantified. We also analyzed the published single-cell RNA sequence data derived from tissue samples from different stages of colorectal carcinogenesis.

Results

The cytokine response in peripheral CD4+ T cells was upregulated during the carcinoma process. The frequency of peripheral regulatory T cells (Tregs) increased in the adenoma and carcinoma stages. While the T follicular helper (Tfh) cell proportion was downregulated in the adenoma and carcinoma processes. Thus, Th cell subsets, especially Tregs and Tfh cells, were involved in colonic diseases. Moreover, the immunological profile characteristics in the HPs were clarified.

Conclusion

We comprehensively analyzed circulating ILCs and adaptive T-cell lymphocyte subtypes in colorectal carcinoma progression. Our results show the immunological profile characteristics and support the involvement of Th subsets, especially Treg and Tfh cell populations, in colonic diseases. These findings significantly enhance our understanding of the immune mechanisms underlying CRC and its precancerous lesions. Further investigation of the Treg and Tfh cells' function in colorectal disease development will provide potential therapeutic targets for monitoring and preventing CRC development.

The oncogenic mechanisms of the Janus kinase-signal transducer and activator of transcription pathway in digestive tract tumors

Digestive tract tumors are heterogeneous and involve the dysregulation of multiple signaling pathways. The Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathway plays a notable role in the oncogenesis of digestive tract tumors. Typically activated by pro-inflammatory cytokines, it regulates important biological processes, such as cell growth, differentiation, apoptosis, immune responses, and inflammation. The aberrant activation of this pathway manifests in different forms, including mutations in JAKs, overexpression of cytokine receptors, and sustained STAT activation, and contributes to promoting the malignant characteristics of cancer cells, including uncontrolled proliferation, resistance to apoptosis, enhanced invasion and metastasis, angiogenesis, acquisition of stem-like properties, and drug resistance. Numerous studies have shown that aberrant activation of the JAK-STAT pathway is closely related to the development and progression of digestive tract tumors, contributing to tumor survival, angiogenesis, changes in the tumor microenvironment, and even immune escape processes. In addition, this signaling pathway also affects the sensitivity of digestive tract tumors to chemotherapy and targeted therapy. Therefore, it is crucial to comprehensively understand the oncogenic mechanisms underlying the JAK-STAT pathway in order to develop effective therapeutic strategies against digestive tract tumors. Currently, several JAK-STAT inhibitors are undergoing clinical and preclinical trials as potential treatments for various human diseases. However, further investigation is required to determine the role of this pathway, as well as the effectiveness and safety of its inhibitors, especially in the context of digestive tract tumors. In this review, we provide an overview of the structure, classic activation, and negative regulation of the JAK-STAT pathway. Furthermore, we discuss the pathogenic mechanisms of JAK-STAT signaling in different digestive tract tumors, with the aim of identifying potential novel therapeutic targets. Video Abstract.

Gains of 20q11.21 in human pluripotent stem cells: Insights from cancer research

The genetic abnormalities observed in hPSC cultures worldwide have been suggested to pose an important hurdle in their safe use in regenerative medicine due to the possibility of oncogenic transformation by mutant cells in the patient posttransplantation. One of the best-characterized genetic lesions in hPSCs is the gain of 20q11.21, found in 20% of hPSC lines worldwide, and strikingly, also amplified in 20% of human cancers. In this review, we have curated the existing knowledge on the incidence of this mutation in hPSCs and cancer, explored the significance of chromosome 20q11.21 amplification in cancer progression, and reviewed the oncogenic role of the genes in the smallest common region of gain, to shed light on the significance of this mutation in hPSC-based cell therapy. Lastly, we discuss the state-of-the-art strategies devised to detect aneuploidies in hPSC cultures, avoid genetic changes in vitro cultures of hPSCs, and strategies to eliminate genetically abnormal cells from culture.

Development and validation of a nomogram predictive model for colorectal adenoma with low-grade intraepithelial neoplasia using routine laboratory tests: A single-center case-control study in China

Background

Colorectal cancer (CRC) is the third most common cancer in the world and has a high mortality rate. Colorectal adenoma (CRA) is precancerous lesions of CRC. The purpose of the present study was to construct a nomogram predictive model for CRA with low-grade intraepithelial neoplasia (LGIN) in order to identify high-risk individuals, facilitating early diagnosis and treatment, and ultimately reducing the incidence of CRC.

Methods

We conducted a single-center case-control study. Based on the results of colonoscopy and pathology, 320 participants were divided into the CRA group and the control group, the demographic and laboratory test data were collected. A development cohort (n = 223) was used for identifying the risk factors for CRA with LGIN and to develop a predictive model, followed by an internal validation. An independent validation cohort (n = 97) was used for external validation. Receiver operating characteristic curve, calibration plot and decision curve analysis were used to evaluate discrimination ability, accuracy and clinical practicability of the model.

Results

Four predictors, namely sex, age, albumin and monocyte count, were included in the predictive model. In the development cohort, internal validation and external validation cohort, the area under the curve (AUC) of this risk predictive model were 0.946 (95%CI: 0.919-0.973), 0.909 (95 % CI: 0.869-0.940) and 0.928 (95%CI: 0.876-0.980), respectively, which demonstrated the model had a good discrimination ability. The calibration plots showed a good agreement and the decision curve analysis (DCA) suggested the predictive model had a high clinical net benefit.

Conclusion

The nomogram model exhibited good performance in predicting CRA with LGIN, which can aid in the early detection of high-risk patients, improve early treatment, and ultimately reduce the incidence of CRC.

Tumor-infiltrating gamma delta T-cells reveal exhausted subsets with remarkable heterogeneity in colorectal cancer

The γδT-cells recognize infected or transformed cells. However, unlike αβT-cells, γδT-cells are innate-like immune cells, with no major histocompatibility complex restriction requirements. γδT-cells are the main population of intestinal intraepithelial lymphocytes (IELs) and are associated with the antitumor immune response, particularly in colorectal cancer (CRC). Although CD8+ T-cells exhibit dysfunction and even exhaustion in the tumor microenvironment (TME), which contributes to tumor immune escape, whether the same applies to tumor-infiltrating (TI)-γδT-cells is not completely understood. Here, we sought to investigate the expression pattern of inhibitory receptors and functional state of TI-γδT-cells, and reveal the features of exhausted TI-γδT-cells in the CRC TME. We demonstrated that TI-γδT-cells exhibited exhaustion phenotypes and displayed more severe functional exhaustion than TI-CD8+ T-cells or NK-cells in the TME of CRC. In addition, scRNA-seq analysis of TI-γδT-cells revealed three exhausted subsets with remarkable heterogeneity. The presence of three heterogeneous exhausted γδT-cell (Tex) populations, including Texprog , Textran and Texterm were further confirmed by flow cytometry, on the basis of PD-1 and TIM-3 expression. Finally, we revealed that c-Maf not only contributed to γδT-cell exhaustion via upregulation of inhibitory receptors, but also involved in the exhaustion of CD8+ T and NK-cells. c-Maf may also be an important contributor to γδT-cell exhaustion in CRC patients. These findings indicated that TI-γδT-cells exhibit phenotypic and functional exhaustion in the CRC TME. The revealed features of exhausted TI-γδT-cells may provide help for understanding the mechanisms and the association of γδT-cell exhaustion with tumor development and pathogenesis.

Integrative analysis reveals a four-gene signature for predicting survival and immunotherapy response in colon cancer patients using bulk and single-cell RNA-seq data

Background

Colon cancer (CC) ranks as one of the leading causes of cancer-related mortality globally. Single-cell transcriptome sequencing (scRNA-seq) offers precise gene expression data for distinct cell types. This study aimed to utilize scRNA-seq and bulk transcriptome sequencing (bulk RNA-seq) data from CC samples to develop a novel prognostic model.

Methods

scRNA-seq data was downloaded from the GSE161277 database. R packages including "Seurat", "Harmony", and "singleR" were employed to categorize eight major cell types within normal and tumor tissues. By comparing tumor and normal samples, differentially expressed genes (DEGs) across these major cell types were identified. Gene Ontology (GO) enrichment analyses of DEGs for each cell type were conducted using "Metascape". DEGs-based signature construction involved Cox regression and least absolute shrinkage operator (LASSO) analyses, performed on The Cancer Genome Atlas (TCGA) training cohort. Validation occurred in the GSE39582 and GSE33382 datasets. The expression pattern of prognostic genes was verified using spatial transcriptome sequencing (ST-seq) data. Ultimately, an established prognostic nomogram based on the gene signature and age was established and calibrated. Sensitivity to chemotherapeutic drugs was predicted with the "oncoPredict" R package.

Results

Using scRNA-Seq data, we examined 33,213 cells, categorizing them into eight cell types within normal and tumor samples. GO enrichment analysis revealed various cancer-related pathways across DEGs in these cell types. Among the 55 DEGs identified via univariate Cox regression, four independent prognostic genes emerged: PTPN6, CXCL13, SPINK4, and NPDC1. Expression validation through ST-seq confirmed PTPN6 and CXCL13 predominance in immune cells, while SPINK4 and NPDC1 were relatively epithelial cell-specific. Creating a four-gene prognostic signature, Kaplan-Meier survival analyses emphasized higher risk scores correlating with unfavorable prognoses, confirmed across training and validation cohorts. The risk score emerged as an independent prognostic factor, supported by a reliable nomogram. Intriguingly, drug sensitivity analysis unveiled contrasting anti-cancer drug responses in the two risk groups, suggesting significant clinical implications.

Conclusion

We developed a novel prognostic four-gene risk model, and these genes may act as potential therapeutic targets for CC.

Identification of GUCA2A and COL3A1 as prognostic biomarkers in colorectal cancer by integrating analysis of RNA-Seq data and qRT-PCR validation

By 2030, it is anticipated that there will be 2.2 million new instances of colorectal cancer worldwide, along with 1.1 million yearly deaths. Therefore, it is critical to develop novel biomarkers that could help in CRC early detection. We performed an integrated analysis of four RNA-Seq data sets and TCGA datasets in this study to find novel biomarkers for diagnostic, prediction, and as potential therapeutic for this malignancy, as well as to determine the molecular mechanisms of CRC carcinogenesis. Four RNA-Seq datasets of colorectal cancer were downloaded from the Sequence Read Archive (SRA) database. The metaSeq package was used to integrate differentially expressed genes (DEGs). The protein-protein interaction (PPI) network of the DEGs was constructed using the string platform, and hub genes were identified using the cytoscape software. The gene ontology and KEGG pathway enrichment analysis were performed using enrichR package. Gene diagnostic sensitivity and its association to clinicopathological characteristics were demonstrated by statistical approaches. By using qRT-PCR, GUCA2A and COL3A1 were examined in colon cancer and rectal cancer. We identified 5037 differentially expressed genes, including (4752 upregulated, 285 downregulated) across the studies between CRC and normal tissues. Gene ontology and KEGG pathway analyses showed that the highest proportion of up-regulated DEGs was involved in RNA binding and RNA transport. Integral component of plasma membrane and mineral absorption pathways were identified as containing down-regulated DEGs. Similar expression patterns for GUCA2A and COL3A1 were seen in qRT-PCR and integrated RNA-Seq analysis. Additionally, this study demonstrated that GUCA2A and COL3A1 may play a significant role in the development of CRC.

Construction of a prognostic model for colorectal adenocarcinoma based on Zn transport-related genes identified by single-cell sequencing and weighted co-expression network analysis

Background

Colorectal cancer (CRC) is one of the most prevalent malignancies and the third most lethal cancer globally. The most reported histological subtype of CRC is colon adenocarcinoma (COAD). The zinc transport pathway is critically involved in various tumors, and its anti-tumor effect may be through improving immune function. However, the Zn transport pathway in COAD has not been reported.

Methods

The determination of Zn transport-related genes in COAD was carried out through single-cell analysis of the GSE 161277 obtained from the GEO dataset. Subsequently, a weighted co-expression network analysis of the TCGA cohort was performed. Then, the prognostic model was conducted utilizing univariate Cox regression and least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Functional enrichment, immune microenvironment, and survival analyses were also carried out. Consensus clustering analysis was utilized to verify the validity of the prognostic model and explore the immune microenvironment. Ultimately, cell experiments, including CCK-8,transwell and scratch assays, were performed to identify the function of LRRC59 in COAD.

Results

According to the Zn transport-related prognostic model, the individuals with COAD in TCGA and GEO databases were classified into high- and low-risk groups. The group with low risk had a comparatively more favorable prognosis. Two groups had significant variations in the immune infiltration, MHC, and the expression of genes related to the immune checkpoint. The cell experiments indicated that the proliferation, migration, and invasion of the HCT-116, DLD-1, and RKO cell lines were considerably increased after LRRC59 knockdown. It proved that LRRC59 was indeed a protective factor for COAD.

Conclusion

A prognostic model for COAD was developed using zinc transport-related genes. This model can efficiently assess the immune microenvironment and prognosis of individuals with COAD. Subsequently, the function of LRRC59 in COAD was validated via cell experiments, highlighting its potential as a biomarker.

Chronic exposure to UVB induces nephritis and gut microbiota dysbiosis in mice based on the integration of renal transcriptome profiles and 16S rRNA sequencing data

Ultraviolet (UV) is a common and abundant environmental factor that affects daily life. Although the effects of UV radiation on the skin have been extensively reported, studies on the influence of UV radiation on internal organs are still limited. This study aimed to evaluate the influence of UVB exposure on the kidney of mice and to investigate the possible mechanism. In the present study, histopathology changes, oxidative stress, and inflammatory response were used to evaluate the kidney and colon injury induced by UVB exposure. The results showed that the 14-week chronic skin exposure to UVB triggers a kidney injury response characterized by macrophage infiltration, elevated oxidative stress as well as inflammatory and injury markers. The RNA sequencing demonstrated that chronic UVB exposure could alter the kidney transcriptomic profile distinguished by the regulation of genes involved in the Notch signaling pathway, JAK-STAT signaling pathway, and ECM-receptor interaction. Besides, chronic UVB exposure also resulted in gut dysbiosis, manifested as colon macrophage infiltration, stimulated inflammatory responses, impaired barrier integrity, and microbiota structural and functional disorders. The Spearman analysis results further revealed a strong correlation between gut microbiota and kidney injury. In conclusion, skin chronic exposure to UVB causes nephritis and gut microbiota dysbiosis in mice, and these findings provide new insight into the underlying risks of chronic UVB exposure to human wellness.

Single-cell RNA-Seq and bulk RNA-Seq reveal reliable diagnostic and prognostic biomarkers for CRC

PURPOSE

The potential role of epithelium-specific genes through the adenoma-carcinoma sequence in the development of colorectal cancer (CRC) remains unknown. Therefore, we integrated single-cell RNA sequencing and bulk RNA sequencing data to select diagnosis and prognosis biomarkers for CRC.

METHODS

The CRC scRNA-seq dataset was used to describe the cellular landscape of normal intestinal mucosa, adenoma and CRC and to further select epithelium-specific clusters. Differentially expressed genes (DEGs) of epithelium-specific clusters were identified between intestinal lesion and normal mucosa in the scRNA-seq data throughout the adenoma-carcinoma sequence. Diagnostic biomarkers and prognostic biomarker (the risk score) for CRC were selected in the bulk RNA-seq dataset based on DEGs shared by the adenoma epithelium-specific cluster and the CRC epithelium-specific cluster (shared-DEGs).

RESULTS

Among the 1063 shared-DEGs, we selected 38 gene expression biomarkers and 3 methylation biomarkers that had promising diagnostic power in plasma. Multivariate Cox regression identified  174 shared-DEGs as prognostic genes for CRC. We combined 1000 times LASSO-Cox regression and two-way stepwise regression to select 10 prognostic shared-DEGs to construct the risk score in the CRC meta-dataset. In the external validation dataset, the 1- and 5-year AUCs of the risk score were higher than those of stage, the pyroptosis-related genes (PRG) score and the cuproptosis-related genes (CRG) score. In addition, the risk score was closely associated with the immune infiltration of CRC.

CONCLUSION

The combined analysis of the scRNA-seq dataset and the bulk RNA-seq dataset in this study provides reliable biomarkers for the diagnosis and prognosis of CRC.

Integrated single-cell and bulk RNA sequencing analysis identifies a prognostic signature related to ferroptosis dependence in colorectal cancer

Ferroptosis is an iron-dependent form of cell death induced by lipid oxidation with an essential role in diseases, including cancer. Since prognostic value of ferroptosis-dependent related genes (FDRGs) in colorectal cancer (CRC) remains unclear, we explored the significance of FDRGs in CRC through comprehensive single-cell analysis. We downloaded the GSE161277 dataset for single-cell analyses and calculated the ferroptosis-dependent gene score (FerrScore) for each cell type. According to each cell type-specific median FerrScore, we categorized the cells into low- and high-ferroptosis groups. By analyzing differentially-expressed genes across the two groups, we identified FDRGs. We further screened these prognosis-related genes used to develop a prognostic signature and calculated its correlation with immune infiltration. We also compared immune checkpoint gene efficacy among different risk groups, and qRT-PCR was performed in colorectal normal and cancer cell lines to explore whether the signature genes could be used as clinical prognostic indicators. In total, 523 FDRGs were identified. A prognostic signature including five signature genes was constructed, and patients were divided into two risk groups. The high-risk group had poor survival rates and displayed high levels of immune infiltration. Our newly developed ferroptosis-based prognostic signature possessed a high predictive ability for CRC.

Dynamic heterogeneity of colorectal cancer during progression revealed clinical risk-associated cell types and regulations in single-cell resolution and spatial context

Background

Tumor heterogeneity is contributed by tumor cells and the microenvironment. Dynamics of tumor heterogeneity during colorectal cancer (CRC) progression have not been elucidated.

Methods

Eight single-cell RNA sequencing (scRNA-seq) data sets of CRC were included. Milo was utilized to reveal the differential abundance of cell clusters during progression. The differentiation trajectory was imputed by using the Palantir algorithm and metabolic states were assessed by using scMetabolism. Three spatial transcription sequencing (ST-seq) data sets of CRC were used to validate cell-type abundances and colocalization. Cancer-associated regulatory hubs were defined as communication networks affecting tumor biological behaviors. Finally, quantitative reverse transcription polymerase chain reaction and immunohistochemistry staining were performed for validation.

Results

TM4SF1⁺, SOX4⁺, and MKI67⁺ tumor cells; CXCL12⁺ cancer-associated fibroblasts; CD4⁺ resident memory T cells; Treg; IgA⁺ plasma cells; and several myeloid subsets were enriched in stage IV CRC, most of which were associated with overall survival of patients. Trajectory analysis indicated that tumor cells from patients with advanced-stage CRC were less differentiated, when metabolic heterogeneity showed a highest metabolic signature in terminal states of stromal cells, T cells, and myeloid cells. Moreover, ST-seq validated cell-type abundance in a spatial context and also revealed the correlation of immune infiltration between tertiary lymphoid structures and tumors followed by validation in our cohort. Importantly, analysis of cancer-associated regulatory hubs revealed a cascade of activated pathways including leukocyte apoptotic process, MAPK pathway, myeloid leukocyte differentiation, and angiogenesis during CRC progression.

Conclusions

Tumor heterogeneity was dynamic during progression, with the enrichment of immunosuppressive Treg, myeloid cells, and fibrotic cells. The differential state of tumor cells was associated with cancer staging. Assessment of cancer-associated regulatory hubs suggested impaired antitumor immunity and increased metastatic ability during CRC progression.

Insights on hematopoietic cell kinase: An oncogenic player in human cancer

Hematopoietic cell kinase (Hck) is a member of the Src family and is expressed in hematopoietic cells. By regulating multiple signaling pathways, HCK can interact with multiple receptors to regulate signaling events involved in cell adhesion, proliferation, migration, invasion, apoptosis, and angiogenesis. However, aberrant expression of Hck in various hematopoietic cells and solid tumors plays a crucial role in tumor-related properties, including cell proliferation and epithelial-mesenchymal transition. In addition, Hck signaling regulates the function of immune cells such as macrophages, contributing to an immunosuppressive tumor microenvironment. The clinical success of various kinase inhibitors targeting the Src kinase family has validated the efficacy of targeting Src, and therapies with highly selective Hck kinase inhibitors are in clinical trials. This article reviews Hck inhibition as an emerging cancer treatment strategy, focusing on the expressions and functions of Hck in tumors and its impact on the tumor microenvironment. It also explores preclinical and clinical pharmacological strategies for Hck targeting to shed light on Hck-targeted tumor therapy.

Dissecting the heterogeneities of the tumor microenvironment between metastatic and nonmetastatic primary colorectal cancer patients by single-cell RNA sequencing

AIMS

One of the main factors hampering the long-term prognosis of colorectal cancer (CRC) patients is distant metastasis. However, the driving factors of CRC metastasis have not been clarified at the single-cell level, which limits the in-depth study of accurate prediction and prevention of CRC metastasis to improve the prognosis.

MATERIALS AND METHODS

Heterogeneities in the tumor microenvironment (TME) between metastatic and nonmetastatic CRC were investigated by single-cell RNA (scRNA) sequencing data. In detail, 50,462 single cells from 20 primary CRC samples, including 40,910 cells from nonmetastatic CRC (M0 group) and 9552 cells from metastatic CRC (M1 group), were systematically analyzed in this study.

KEY FINDINGS

Based on the single-cell atlas, we revealed that cancer cells and fibroblasts accounted for relatively high proportions in metastatic CRC compared with nonmetastatic CRC. Moreover, two specific cancer cell subtypes (FGGY⁺SLC6A6⁺ and IGFBP3⁺KLK7⁺ cancer cells) and three specific fibroblast subtypes (ADAMTS6⁺CAPG⁺, PIM1⁺SGK1⁺ and CA9⁺UPP1⁺ fibroblasts) in metastatic CRC were identified. The functional and differentiation characteristics of these specific cell subclusters were elucidated by enrichment and trajectory analyses.

SIGNIFICANCE

These results provide fundamental knowledge for future in-depth research to screen effective methods and drugs to predict and prevent CRC metastasis to improve prognosis.

Low expression of the dynamic network markers FOS/JUN in pre-deteriorated epithelial cells is associated with the progression of colorectal adenoma to carcinoma

BACKGROUND

Deterioration of normal intestinal epithelial cells is crucial for colorectal tumorigenesis. However, the process of epithelial cell deterioration and molecular networks that contribute to this process remain unclear.

METHODS

Single-cell data and clinical information were downloaded from the Gene Expression Omnibus (GEO) database. We used the recently proposed dynamic network biomarker (DNB) method to identify the critical stage of epithelial cell deterioration. Data analysis and visualization were performed using R and Cytoscape software. In addition, Single-Cell rEgulatory Network Inference and Clustering (SCENIC) analysis was used to identify potential transcription factors, and CellChat analysis was conducted to evaluate possible interactions among cell populations. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set variation analysis (GSVA) analyses were also performed.

RESULTS

The trajectory of epithelial cell deterioration in adenoma to carcinoma progression was delineated, and the subpopulation of pre-deteriorated epithelial cells during colorectal cancer (CRC) initialization was identified at the single-cell level. Additionally, FOS/JUN were identified as biomarkers for pre-deteriorated epithelial cell subpopulations in CRC. Notably, FOS/JUN triggered low expression of P53-regulated downstream pro-apoptotic genes and high expression of anti-apoptotic genes through suppression of P53 expression, which in turn inhibited P53-induced apoptosis. Furthermore, malignant epithelial cells contributed to the progression of pre-deteriorated epithelial cells through the GDF signaling pathway.

CONCLUSIONS

We demonstrated the trajectory of epithelial cell deterioration and used DNB to characterize pre-deteriorated epithelial cells at the single-cell level. The expression of DNB-neighboring genes and cellular communication were triggered by DNB genes, which may be involved in epithelial cell deterioration. The DNB genes FOS/JUN provide new insights into early intervention in CRC.

Diabetes mellitus induces a novel inflammatory network involving cancer progression: Insights from bioinformatic analysis and in vitro validation

Background

Patients with diabetes mellitus (DM) have a higher incidence of malignant tumors than people without diabetes, but the underlying molecular mechanisms are still unclear.

Methods

To investigate the link between DM and cancer, we screened publicly available databases for diabetes and cancer-related genes (DCRGs) and constructed a diabetes-based cancer-associated inflammation network (DCIN). We integrated seven DCRGs into the DCIN and analyzed their role in different tumors from various perspectives. We also investigated drug sensitivity and single-cell sequencing data in colon adenocarcinoma as an example. In addition, we performed in vitro experiments to verify the expression of DCRGs and the arachidonic acid metabolic pathway.

Results

Seven identified DCRGs, including PPARG, MMP9, CTNNB1, TNF, TGFB1, PTGS2, and HIF1A, were integrated to construct a DCIN. The bioinformatics analysis showed that the expression of the seven DCRGs in different tumors was significantly different, which had varied effects on diverse perspectives. Single-cell sequencing analyzed in colon cancer showed that the activity of the DCRGs was highest in Macrophage and the lowest in B cells among all cell types in adenoma and carcinoma tissue. In vitro experiments showed that the DCRGs verified by western bolt and PEG2 verified by ELISA were all highly expressed in COAD epithelial cells stimulated by high glucose.

Conclusion

This study, for the first time, constructed a DCIN, which provides novel insights into the underlying mechanism of how DM increases tumor occurrence and development. Although further research is required, our results offer clues for new potential therapeutic strategies to prevent and treat malignant tumors.

Comprehensive genomics analysis of aging related gene signature to predict the prognosis and drug resistance of colon adenocarcinoma

Background: Colon adenocarcinoma (COAD) is a heterogeneous tumor and senescence is crucial in the occurrence of cancer. This study aimed to identify senescence-based subtypes and construct a prognostic signature to predict the prognosis and guide immunotherapy or chemotherapy decisions for COAD patients. Methods: Based on the single-cell RNA sequencing (scRNA-seq) data of 13 samples from the Gene Expression Omnibus (GEO) database, we assessed cellular senescence characteristics. Transcriptome data, copy number variations (CNVs) and single nucleotide variations (SNVs) data were obtained from The Cancer Genome Atlas (TCGA) database. GSE39582 and GSE17537 were used for validation. Senescence subtypes were identified using unsupervised consensus clustering analysis, and a prognostic signature was developed using univariate Cox analysis and least absolute shrinkage and selection operator (LASSO). Response of risk groups to chemotherapy was predicted using the half-maximal inhibitory concentration (IC50) values. We further analyzed the relationship between risk gene expression and methylation level. The prediction performance was assessed by nomogram. Results: Senescence-related pathways were highly enriched in malignant cells and bulk RNA-seq verified cellular senescence. Three senescence subtypes were identified, in which patients in clust3 had poorest prognosis and higher T stage, accompanied with higher tumor mutation burden (TMB) and mutations, activated inflammatory response, more immune cell infiltration, and higher immune escape tendency. A senescence-based signature using 11 genes (MFNG, GPRC5B, TNNT1, CCL22, NOXA1, PABPC1L, PCOLCE2, MID2, CPA3, HSPA1A, and CALB1) was established, and accurately predicted a lower prognosis in high risk patients. Its robustness was validated by external cohort. Low risk patients were more sensitive to small molecule drugs including Erlotinib, Sunitinib, MG-132, CGP-082996, AZ628, Sorafenib, VX-680, and Z-LLNle-CHO. Risk score was an independent prognostic factor and nomogram confirmed its reliability. Four risk genes (CALB1, CPA3, NOXA1, and TNNT1) had significant positive correlation with their methylation level, while six genes (CCL22, GPRC5B, HSPA1A, MFNG, PABPC1L, and PCOLCE2) were negatively correlated with their methylation level. Conclusion: This study provides novel understanding of heterogeneity in COAD from the perspective of senescence, and develops signatures for prognosis prediction in COAD.

Targeting cancer stem cells and signalling pathways through phytochemicals: A promising approach against colorectal cancer

BACKGROUND

Cancer stem cells (CSCs) are strongly associated with high tumourigenicity, chemotherapy or radiotherapy resistance, and metastasis and recurrence, particularly in colorectal cancer (CRC). Therefore, targeting CSCs may be a promising approach. Recently, discovery and research on phytochemicals that effectively target colorectal CSCs have been gaining popularity because of their broad safety profile and multi-target and multi-pathway modes of action.

PURPOSE

This review aimed to elucidate and summarise the effects and mechanisms of phytochemicals with potential anti-CSC agents that could contribute to the better management of CRC.

METHODS

We reviewed PubMed, EMBASE, Web of Science, Ovid, ScienceDirect and China National Knowledge Infrastructure databases from the original publication date to March 2022 to review the mechanisms by which phytochemicals inhibit CRC progression by targeting CSCs and their key signalling pathways. Phytochemicals were classified and summarised based on the mechanisms of action.

RESULTS

We observed that phytochemicals could affect the biological properties of colorectal CSCs. Phytochemicals significantly inhibit self-renewal, migration, invasion, colony formation, and chemoresistance and induce apoptosis and differentiation of CSCs by regulating the Wnt/β-catenin pathway (e.g., diallyl trisulfide and genistein), the phosphatidylinositol-3-kinase/protein kinase B/mammalian target of rapamycin pathway (e.g., caffeic acid and piperlongumine), the neurogenic locus notch homolog protein pathway (e.g., honokiol, quercetin, and α-mangostin), the Janus kinase-signal transducer and activator of transcription pathway (e.g., curcumin, morin, and ursolic acid), and other key signalling pathways. It is worth noting that several phytochemicals, such as resveratrol, silibinin, evodiamine, and thymoquinone, highlight multi-target and multi-pathway effects in restraining the malignant biological behaviour of CSCs.

CONCLUSIONS

This review demonstrates the potential of targeted therapies for colorectal CSCs using phytochemicals. Phytochemicals could serve as novel therapeutic agents for CRC and aid in drug development.

Silence of a dependence receptor CSF1R in colorectal cancer cells activates tumor-associated macrophages

BACKGROUND

Colony-stimulating factor 1 receptor (CSF1R), a classic tyrosine kinase receptor, has been identified as a proto-oncogene in multiple cancers. The CSF1/CSF1R axis is essential for the survival and differentiation of M2-phenotype tumor-associated macrophages (M2 TAMs). However, we found here that the CSF1R expression was abnormally down-regulated in colorectal cancer (CRC), and its biological functions and underlying mechanisms have become elusive in CRC progression.

METHODS

The expression of class III receptor tyrosine kinases in CRC and normal intestinal mucosa was accessed using The Cancer Genome Atlas and Gene Expression Omnibus datasets and was further validated by our tested cohort. CSF1R was reconstructed in CRC cells to identify its biological functions in vitro and in vivo. We compared CSF1R expression and methylation differences between CRC cells and macrophages. Furthermore, a co-culture system was used to mimic a competitive mechanism between CSF1R-overexpressed CRC cells and M2-like macrophages. We utilized a CSF1R inhibitor PLX3397 to ablate M2 TAMs and evaluated its efficacy on CRC treatment in animal models.

RESULTS

We found here that the CSF1R is silenced in CRC, and the reintroduced expression of the receptor in CRC cells can be cleaved by caspases and constrain tumor growth in vitro and in vivo, functioning as a tumor suppressor gene. We further identified CSF1R as a novel dependence receptor, which has the potential to act as either a tumor suppressor gene or an oncogene, depending on its activated state. In CRC tumors, CSF1R expression is enriched in TAMs, and its expression is associated with poor prognosis in patients ith CRC. In a co-culture system, CRC cells expressing CSF1R compete with M2-like macrophages for CSF1R ligands, resulting in a decrease in CSF1R activation and cell proliferation in macrophages. Blocking CSF1R by PLX3397 could deplete M2 TAMs and augments CD8⁺ T cell infiltration, effectively inhibiting tumor growth and metastasis and improving responses to chemotherapy and immunotherapy.

CONCLUSION

Our findings revealed that CSF1R is a novel identified dependence receptor silenced in CRC. The silence abalienates its ligands to stimulate CSF1R expressed on M2 TAMs, which is an appealing therapeutic target for M2 TAM depletion and CRC treatment.

Applications of human organoids in the personalized treatment for digestive diseases

Digestive system diseases arise primarily through the interplay of genetic and environmental influences; there is an urgent need in elucidating the pathogenic mechanisms of these diseases and deploy personalized treatments. Traditional and long-established model systems rarely reproduce either tissue complexity or human physiology faithfully; these shortcomings underscore the need for better models. Organoids represent a promising research model, helping us gain a more profound understanding of the digestive organs; this model can also be used to provide patients with precise and individualized treatment and to build rapid in vitro test models for drug screening or gene/cell therapy, linking basic research with clinical treatment. Over the past few decades, the use of organoids has led to an advanced understanding of the composition of each digestive organ and has facilitated disease modeling, chemotherapy dose prediction, CRISPR-Cas9 genetic intervention, high-throughput drug screening, and identification of SARS-CoV-2 targets, pathogenic infection. However, the existing organoids of the digestive system mainly include the epithelial system. In order to reveal the pathogenic mechanism of digestive diseases, it is necessary to establish a completer and more physiological organoid model. Combining organoids and advanced techniques to test individualized treatments of different formulations is a promising approach that requires further exploration. This review highlights the advancements in the field of organoid technology from the perspectives of disease modeling and personalized therapy.

Loss of MHC-I antigen presentation correlated with immune checkpoint blockade tolerance in MAPK inhibitor-resistant melanoma

Immune checkpoint blockade and MAPK-targeted combined therapy is a promising regimen for advanced melanoma patients. However, the clinical benefit from this combo regimen remains limited, especially in patients who acquired resistance to MAPK-targeted therapy. Here, we systematically characterized the immune landscape during MAPK-targeted therapy in patients and mouse melanoma models. We observed that both the abundance of tumor-infiltrated T cells and the expression of immune-related genes were upregulated in the drug-responsive period, but downregulated in the resistance period, implying that acquired drug resistance dampens the antitumor immune response. Further transcriptomic dissection indicated that loss of MHC-I antigen presentation on tumor cells plays a critical role in the reduction of T cell infiltration during drug resistance. Survival analysis demonstrates that loss of antigen presentation and reduction of T-cell infiltration during acquired drug resistance are associated with poorer clinical response and prognosis of anti-PD-1 therapy in melanoma patients. In addition, we identified that alterations in the MAPK inhibitor resistance-related oncogenic signaling pathway closely correlated with deficiency of MHC-I antigen presentation, including activation of the PI3K-mTOR, MAPK, and Wnt pathways. In conclusion, our research illuminates that decreased infiltration of T cells is associated with acquired drug resistance during MAPK-targeted therapy, which may underlie the cross-resistance to immune checkpoint blockade.

Single-Cell RNA Sequencing Reveals the Role of Epithelial Cell Marker Genes in Predicting the Prognosis of Colorectal Cancer Patients

Single-cell RNA sequencing (scRNA-seq) is increasingly used in studies on gastrointestinal cancers. This study investigated the prognostic value of epithelial cell-associated biomarkers in colorectal cancer (CRC) using scRNA-seq data. We downloaded and analysed scRNA-seq data from four CRC samples from the Gene Expression Omnibus (GEO), and we identified marker genes of malignant epithelial cells (MECs) using CRC transcriptome and clinical data downloaded from The Cancer Genome Atlas (TCGA) and GEO as training and validation cohorts, respectively. In the TCGA training cohort, weighted gene correlation network analysis, univariate Cox proportional hazard model (Cox) analysis, and least absolute shrinkage and selection operator regression analysis were performed on the marker genes of MEC subsets to identify a signature of nine prognostic MEC-related genes (MECRGs) and calculate a risk score based on the signature. CRC patients were divided into high- and low-risk groups according to the median risk score. We found that the MECRG risk score was significantly correlated with the clinical features and overall survival of CRC patients, and that CRC patients in the high-risk group showed a significantly shorter survival time. The univariate and multivariate Cox regression analyses showed that the MECRG risk score can serve as an independent prognostic factor for CRC patients. Gene set enrichment analysis revealed that the MECRG signature genes are involved in fatty acid metabolism, p53 signalling, and other pathways. To increase the clinical application value, we constructed a MECRG nomogram by combining the MECRG risk score with other independent prognostic factors. The validity of the nomogram is based on receiver operating characteristics and calibration curves. The MECRG signature and nomogram models were well validated in the GEO dataset. In conclusion, we established an epithelial cell marker gene-based risk assessment model based on scRNA-seq analysis of CRC samples for predicting the prognosis of CRC patients.

Exploring the Interplay between Metabolism and Tumor Microenvironment Based on Four Major Metabolism Pathways in Colon Adenocarcinoma

Tumor metabolism plays a critical role in tumor progression. However, the interaction between metabolism and tumor microenvironment (TME) has not been comprehensively revealed in colon adenocarcinoma (COAD). We used unsupervised consensus clustering to establish three molecular subtypes (clusters) based on the enrichment score of four major metabolism pathways in TCGA-COAD dataset. GSE17536 was used as a validation dataset. Single-cell RNA sequencing data (GSE161277) was employed to further verify the reliability of subtyping and characterize the correlation between metabolism and TME. Three clusters were identified and they performed distinct prognosis and molecular features. Clust3 had the worst overall survival and the highest enrichment score of glycolysis. 86 differentially expressed genes (DEGs) were identified, in which 11 DEGs were associated with favorable prognosis and 75 DEGs were associated with poor prognosis. Striking correlations were observed between hypoxia and glycolysis, clust3 and hypoxia, and clust3 and angiogenesis (P < 0.001).We constructed a molecular subtyping system which was effective and reliable for predicting COAD prognosis. The 86 identified key DEGs may be greatly involved in COAD progression, and they provide new perspectives and directions for further understanding the mechanism of metabolism in promoting aggressive phenotype by interacting with TME.

Sodium Tanshinone IIA Sulfonate as a Potent IDO1/TDO2 Dual Inhibitor Enhances Anti-PD1 Therapy for Colorectal Cancer in Mice

Although the antitumor efficacy of immune checkpoint blockade (ICB) has been proved in colorectal cancer (CRC), the results are unsatisfactory, presumably owing to the presence of tryptophan metabolism enzymes indoleamine 2,3-dioxygenase 1 (IDO1) and tryptophan 2,3-dioxygenase 2 (TDO2). However, only a few dual inhibitors for IDO1 and TDO2 have been reported. Here, we discovered that sodium tanshinone IIA sulfonate (STS), a sulfonate derived from tanshinone IIA (TSN), reduced the enzymatic activities of IDO1 and TDO2 with a half inhibitory concentration (IC₅₀) of less than 10 μM using enzymatic assays for natural product screening. In IDO1- or TDO2- overexpressing cell lines, STS decreased kynurenine (kyn) synthesis. STS also reduced the percentage of forkhead box P3 (FOXP3) T cells in lymphocytes from the mouse spleen cocultured with CT26. In vivo, STS suppressed tumor growth and enhanced the antitumor effect of the programmed cell death 1 (PD1) antibody. Compared with anti-PD1 (α-PD1) monotherapy, combined with STS had lower level of plasma kynurenine. Immunofluorescence assay suggested that STS decreased the number of FOXP3+ T cells and increased the number of CD8+ T cells in tumors. Flow cytometry analysis of immune cells in tumor tissues demonstrated an increase in the percentage of tumor-infiltrating CD8+ T cells. According to our findings, STS acts as an immunotherapy agent in CRC by inhibiting both IDO1 and TDO2.