Identification of BGN and THBS2 as metastasis-specific biomarkers and poor survival key regulators in human colon cancer by integrated analysis

Integrating bulk RNA-seq and scRNA-seq analyses revealed the function and clinical value of thrombospondins in colon cancer

Background

Acting as mediators in cell-matrix and cell-cell communication, matricellular proteins play a crucial role in cancer progression. Thrombospondins (TSPs), a type of matricellular glycoproteins, are key regulators in cancer biology with multifaceted roles. Although TSPs have been implicated in anti-tumor immunity and epithelial-mesenchymal transition (EMT) in several malignancies, their specific roles to colon cancer remain elusive. Addressing this knowledge gap is essential, as understanding the function of TSPs in colon cancer could identify new therapeutic targets and prognostic markers.

Methods

Analyzing 1981 samples from 10 high-throughput datasets, including six bulk RNA-seq, three scRNA-seq, and one spatial transcriptome dataset, our study investigated the prognostic relevance, risk stratification value, immune heterogeneity, and cellular origin of TSPs, as well as their influence on cancer-associated fibroblasts (CAFs). Utilizing survival analysis, unsupervised clustering, and functional enrichment, along with multiple correlation analyses of the tumor-microenvironment (TME) via Gene Set Variation Analysis (GSVA), spatial localization, Monocle2, and CellPhoneDB, we provided insights into the clinical and cellular implications of TSPs.

Results

First, we observed significant upregulation of THBS2 and COMP in colon cancer, both of which displayed significant prognostic value. Additionally, we detected a significant positive correlation between TSPs and immune cells, as well as marker genes of EMT. Second, based on TSPs expression, patients were divided into two clusters with distinct prognoses: the high TSPs expression group (TSPs-H) was characterized by pronounced immune and stromal cell infiltration, and notably elevated T-cell exhaustion scores. Subsequently, we found that THBS2 and COMP may be associated with the differentiation of CAFs into pan-iCAFs and pan-dCAFs, which are known for their heightened matrix remodeling activities. Moreover, THBS2 enhanced CAFs communication with vascular endothelial cells and monocyte-macrophages. CAFs expressing THBS2 (THBS2+ CAFs) demonstrated higher scores across multiple signaling pathways, including angiogenic, EMT, Hedgehog, Notch, Wnt, and TGF-β, when compared to THBS2- CAFs. These observations suggest that THBS2 may be associated with stronger pro-carcinogenic activity in CAFs.

Conclusions

This study revealed the crucial role of TSPs and the significant correlation between THBS2 and CAFs interactions in colon cancer progression, providing valuable insights for targeting TSPs to mitigate cancer progression.

Plasma circulating tumor DNA unveils the efficacy of PD-1 inhibitors and chemotherapy in advanced gastric cancer

Programmed Death Receptor 1 (PD-1) inhibitors, when combined with chemotherapy, have exhibited notable effectiveness in enhancing the survival outcomes of patients afflicted with advanced gastric cancer. However, it is important to acknowledge that not all patients derive substantial benefits from this therapeutic approach, highlighting the crucial necessity of identifying efficacious biomarkers to inform immunotherapy interventions. In this study, we sought to investigate the predictive utility of circulating tumor DNA (ctDNA) as a biomarker in a cohort of 30 patients diagnosed with advanced gastric cancer, all of whom underwent first-line treatment involving PD-1 inhibitor administration alongside chemotherapy. We procured peripheral blood samples both at baseline and following the completion of two treatment cycles. Additionally, baseline tissue specimens were collected for the purpose of genomic alteration assessment, employing both 47-gene and 737-gene next-generation sequencing panels for plasma and tumor tissue, respectively. We delineated a ctDNA response as the eradication of maximum variant allele frequencies relative to baseline levels. Notably, the objective response rate among individuals exhibiting a ctDNA response proved significantly superior in comparison to non-responders (P = 0.0073). Furthermore, patients who manifested a ctDNA response experienced markedly prolonged progression-free survival (PFS) and overall survival (OS) when juxtaposed with those devoid of a ctDNA response (median PFS: 15.6 vs. 6.0 months, P = 0.003; median OS: not reached [NR] vs. 9.0 months, P = 0.011). In summation, patients with advanced gastric cancer receiving first-line treatment with PD-1 inhibitors and chemotherapy, dynamic changes in ctDNA can serve as a potential biomarker for predicting treatment efficacy and long-term outcomes.

Unlocking the potential of Escherichia coli K-12: A novel approach for malignancy reduction in colorectal cancer through gene expression modulation

Studies have noted the association between Escherichia coli K-12 (E. coli K-12) and the reduction of malignancy in colorectal cancer (CRC). However, the molecular mechanisms underlying this relationship have not been thoroughly explored. The aim of this study was to identify the genes influenced by E. coli K-12 and their connection to CRC. We identified the genes affected by E. coli K-12 using the GSE50040 dataset. Additionally, we investigated the relationship between the expression of genes affected by E. coli K-12 and CRC using the cancer genome atlas data. The association between the expression of E. coli K-12-affected genes and patient prognosis was investigated using clinical data. Pathways related to CRC and E. coli K-12-related genes were analyzed using the Enrichr tool. Furthermore, we employed a protein-protein interaction (PPI) network to identify hub genes associated with both E. coli K-12 and CRC. To validate our findings, we conducted RT-qPCR analysis on CRC samples and adjacent normal tissue. The results of GSE50040 showed that E. coli K-12 could change the expression of many genes related to CRC in colorectal cell lines. The results showed that E. coli K-12 reduces the expression of several genes linked to the main pathways used by cancer cells, such as the metastasis, WNT, cell proliferation pathway, and mTORC1. It was demonstrated that elevated BGN, FJX1, and LZTS1 expression is linked to a bad prognosis in patients and that E. coli K-12 may be able to lower this expression. Also, based on the PPI network, genes such as KLF4 and CXCL3 were identified as hub genes related to genes affected by E. coli K-12. When KLF4 and CXCL3 expression levels in cancer samples were compared to nearby normal tissue, a significant change in these genes' expression levels was found in CRC. Our findings demonstrated the potential relationship between oncogene genes and genes impacted by E. coli K-12. Also, our findings demonstrated that E. coli K-12 may regulate the expression of genes linked to a high death rate. In summary, the results of this study suggest that E. coli K-12 can be regarded as a significant probiotic with the potential to mitigate the risk of CRC development.

Investigating gene signatures associated with immunity in colon adenocarcinoma to predict the immunotherapy effectiveness using NFM and WGCNA algorithms

Colon adenocarcinoma (COAD), a frequently encountered and highly lethal malignancy of the digestive system, has been the focus of intensive research regarding its prognosis. The intricate immune microenvironment plays a pivotal role in the pathological progression of COAD; nevertheless, the underlying molecular mechanisms remain incompletely understood. This study aims to explore the immune gene expression patterns in COAD, construct a robust prognostic model, and delve into the molecular mechanisms and potential therapeutic targets for COAD liver metastasis, thereby providing critical support for individualized treatment strategies and prognostic evaluation. Initially, we curated a comprehensive dataset by screening 2600 immune-related genes (IRGs) from the ImmPort and InnateDB databases, successfully obtaining a rich data resource. Subsequently, the COAD patient cohort was classified using the non-negative matrix factorization (NMF) algorithm, enabling accurate categorization. Continuing on, utilizing the weighted gene co-expression network analysis (WGCNA) method, we analyzed the top 5000 genes with the smallest p-values among the differentially expressed genes (DEGs) between immune subtypes. Through this rigorous screening process, we identified the gene modules with the strongest correlation to the COAD subpopulation, and the intersection of genes in these modules with DEGs (COAD vs COAD vs Normal colon tissue) is referred to as Differentially Expressed Immune Genes Associated with COAD (DEIGRC). Employing diverse bioinformatics methodologies, we successfully developed a prognostic model (DPM) consisting of six genes derived from the DEIGRC, which was further validated across multiple independent datasets. Not only does this predictive model accurately forecast the prognosis of COAD patients, but it also provides valuable insights for formulating personalized treatment regimens. Within the constructed DPM, we observed a downregulation of CALB2 expression levels in COAD tissues, whereas NOXA1, KDF1, LARS2, GSR, and TIMP1 exhibited upregulated expression levels. These genes likely play indispensable roles in the initiation and progression of COAD and thus represent potential therapeutic targets for patient management. Furthermore, our investigation into the molecular mechanisms and therapeutic targets for COAD liver metastasis revealed associations with relevant processes such as fat digestion and absorption, cancer gene protein polysaccharides, and nitrogen metabolism. Consequently, genes including CAV1, ANXA1, CPS1, EDNRA, and GC emerge as promising candidates as therapeutic targets for COAD liver metastasis, thereby providing crucial insights for future clinical practices and drug development. In summary, this study uncovers the immune gene expression patterns in COAD, establishes a robust prognostic model, and elucidates the molecular mechanisms and potential therapeutic targets for COAD liver metastasis, thereby possessing significant theoretical and clinical implications. These findings are anticipated to offer substantial support for both the treatment and prognosis management of COAD patients.

Arenobufagin inhibits lung metastasis of colorectal cancer by targeting c-MYC/Nrf2 axis

BACKGROUND

Colorectal cancer (CRC) is one of the commonest cancers worldwide. Metastasis is the most common cause of death in patients with CRC. Arenobufagin is an active component of bufadienolides, extracted from toad skin and parotid venom. Arenobufagin reportedly inhibits epithelial-to-mesenchymal transition (EMT) and metastasis in various cancers. However, the mechanism through which arenobufagin inhibits CRC metastasis remains unclear.

PURPOSE

This study aimed to elucidate the molecular mechanisms by which arenobufagin inhibits CRC metastasis.

METHODS

Wound-healing and transwell assays were used to assess the migration and invasion of CRC cells. The expression of nuclear factor erythroid-2-related factor 2 (Nrf2) in the CRC tissues was assessed using immunohistochemistry. The protein expression levels of c-MYC and Nrf2 were detected by immunoblotting. A mouse model of lung metastasis was used to study the effects of arenobufagin on CRC lung metastasis in vivo.

RESULTS

Arenobufagin observably inhibited the migration and invasion of CRC cells by downregulating c-MYC and inactivating the Nrf2 signaling pathway. Pretreatment with the Nrf2 inhibitor brusatol markedly enhanced arenobufagin-mediated inhibition of migration and invasion, whereas pretreatment with the Nrf2 agonist tert‑butylhydroquinone significantly attenuated arenobufagin-mediated inhibition of migration and invasion of CRC cells. Furthermore, Nrf2 knockdown with short hairpin RNA enhanced the arenobufagin-induced inhibition of the migration and invasion of CRC cells. Importantly, c-MYC acts as an upstream modulator of Nrf2 in CRC cells. c-MYC knockdown markedly enhanced arenobufagin-mediated inhibition of the Nrf2 signaling pathway, cell migration, and invasion. Arenobufagin inhibited CRC lung metastasis in vivo. Together, these findings provide evidence that interruption of the c-MYC/Nrf2 signaling pathway is crucial for arenobufagin-inhibited cell metastasis in CRC.

CONCLUSIONS

Collectively, our findings show that arenobufagin could be used as a potential anticancer agent against CRC metastasis. The arenobufagin-targeted c-MYC/Nrf2 signaling pathway may be a novel chemotherapeutic strategy for treating CRC.

Identification of BGN positive fibroblasts as a driving factor for colorectal cancer and development of its related prognostic model combined with machine learning

BACKGROUND

Numerous studies have indicated that cancer-associated fibroblasts (CAFs) play a crucial role in the progression of colorectal cancer (CRC). However, there are still many unknowns regarding the exact role of CAF subtypes in CRC.

METHODS

The data for this study were obtained from bulk, single-cell, and spatial transcriptomic sequencing data. Bioinformatics analysis, in vitro experiments, and machine learning methods were employed to investigate the functional characteristics of CAF subtypes and construct prognostic models.

RESULTS

Our study demonstrates that Biglycan (BGN) positive cancer-associated fibroblasts (BGN + Fib) serve as a driver in colorectal cancer (CRC). The proportion of BGN + Fib increases gradually with the progression of CRC, and high infiltration of BGN + Fib is associated with poor prognosis in terms of overall survival (OS) and recurrence-free survival (RFS) in CRC. Downregulation of BGN expression in cancer-associated fibroblasts (CAFs) significantly reduces migration and proliferation of CRC cells. Among 101 combinations of 10 machine learning algorithms, the StepCox[both] + plsRcox combination was utilized to develop a BGN + Fib derived risk signature (BGNFRS). BGNFRS was identified as an independent adverse prognostic factor for CRC OS and RFS, outperforming 92 previously published risk signatures. A Nomogram model constructed based on BGNFRS and clinical-pathological features proved to be a valuable tool for predicting CRC prognosis.

CONCLUSION

In summary, our study identified BGN + Fib as drivers of CRC, and the derived BGNFRS was effective in predicting the OS and RFS of CRC patients.

THBS1 and THBS2 Enhance the In Vitro Proliferation, Adhesion, Migration and Invasion of Intrahepatic Cholangiocarcinoma Cells

In intrahepatic cholangiocarcinoma (iCCA), thrombospondin 1 (THBS1) and 2 (THBS2) are soluble mediators released in the tumor microenvironment (TME) that contribute to the metastatic spreading of iCCA cells via a lymphatic network by the trans-differentiation of vascular endothelial cells to a lymphatic-like phenotype. To study the direct role of THBS1 and THBS2 on the iCCA cells, well-established epithelial (HuCCT-1) and mesenchymal (CCLP1) iCCA cell lines were subjected to recombinant human THBS1 and THBS2 (rhTHBS1, rhTHBS2) for cellular function assays. Cell growth, cell adhesion, migration, and invasion were all enhanced in both CCLP1 and HuCCT-1 cells by the treatment with either rhTHBS1 or rhTHBS2, although they showed some variability in their intensity of speeding up cellular processes. rhTHBS2 was more intense in inducing invasiveness and in committing the HuCCT-1 cells to a mesenchymal-like phenotype and was therefore a stronger enhancer of the malignant behavior of iCCA cells compared to rhTHBS1. Our data extend the role of THBS1 and THBS2, which are not only able to hinder the vascular network and promote tumor-associated lymphangiogenesis but also exacerbate the malignant behavior of the iCCA cells.

Identification of a basement membrane gene signature for predicting prognosis and estimating the tumor immune microenvironment in prostate cancer

Basement membrane plays an important role in tumor invasion and metastasis, which is closely related to prognosis. However, the prognostic value and biology of basement membrane genes (BMGs) in prostate cancer (PCa) remain unknown. In the TCGA training set, we used differentially expressed gene analysis, protein-protein interaction networks, univariate and multivariate Cox regression, and least absolute shrinkage and selection operator regression to construct a basement membrane-related risk model (BMRM) and validated its effectiveness in the MSKCC validation set. Furthermore, the accurate nomogram was constructed to improve clinical applicability. Patients with PCa were divided into high-risk and low-risk groups according to the optimal cut-off value of the basement membrane-related risk score (BMRS). It was found that BMRS was significantly associated with RFS, T-stage, Gleason score, and tumor microenvironmental characteristics in PCa patients. Further analysis showed that the model grouping was closely related to tumor immune microenvironment characteristics, immune checkpoint inhibitors, and chemotherapeutic drug sensitivity. In this study, we developed a new BMGs-based prognostic model to determine the prognostic value of BMGs in PCa. Finally, we confirmed that THBS2, a key gene in BMRM, may be an important link in the occurrence and progression of PCa. This study provides a novel perspective to assess the prognosis of PCa patients and provides clues for the selection of future personalized treatment regimens.

Cdh11: Roles in different diseases and potential value in disease diagnosis and treatment

Cadherin is a homophilic, Ca2+-dependent cell adhesion glycoprotein that mediates cell-cell adhesion. Among them, Cadherin-11 (CDH11), as a classical cadherin, participates in and influences many crucial aspects of human growth and development. Furthermore, The involvement of CDH11 has been identified in an increasing number of diseases, primarily including various tumorous diseases, fibrotic diseases, autoimmune diseases, neurodevelopmental disorders, and more. In various tumorous diseases, CDH11 acts not only as a tumor suppressor but can also promote migration and invasion of certain tumors through various mechanisms. Likewise, in non-tumorous diseases, CDH11 remains a pivotal factor in disease progression. In this context, we summarize the specific functionalities and mechanisms of CDH11 in various diseases, aiming to gain a more comprehensive understanding of the potential value of CDH11 in disease diagnosis and treatment. This endeavor seeks to provide more effective diagnostic and therapeutic strategies for clinical management across diverse diseases.

RUNX1-Regulated Signaling Pathways in Ovarian Cancer

Ovarian cancer is the leading cause of gynecological death worldwide, and its poor prognosis and high mortality seriously affect the life of ovarian cancer patients. Runt-related transcription factor 1 (RUNX1) has been widely studied in hematological diseases and plays an important role in the occurrence and development of hematological diseases. In recent years, studies have reported the roles of RUNX1 in solid tumors, including the significantly increased expression of RUNX1 in ovarian cancer. In ovarian cancer, the dysregulation of the RUNX1 signaling pathway has been implicated in tumor progression, metastasis, and response to therapy. At the same time, the decreased expression of RUNX1 in ovarian cancer can significantly improve the sensitivity of clinical chemotherapy and provide theoretical support for the subsequent diagnosis and treatment target of ovarian cancer, providing prognosis and treatment options to patients with ovarian cancer. However, the role of RUNX1 in ovarian cancer remains unclear. Therefore, this article reviews the relationship between RUNX1 and the occurrence and development of ovarian cancer, as well as the closely regulated signaling pathways, to provide some inspiration and theoretical support for future research on RUNX1 in ovarian cancer and other diseases.

The Landscape of Small Leucine-Rich Proteoglycan Impact on Cancer Pathogenesis with a Focus on Biglycan and Lumican

Cancer development is a multifactorial procedure that involves changes in the cell microenvironment and specific modulations in cell functions. A tumor microenvironment contains tumor cells, non-malignant cells, blood vessels, cells of the immune system, stromal cells, and the extracellular matrix (ECM). The small leucine-rich proteoglycans (SLRPs) are a family of nineteen proteoglycans, which are ubiquitously expressed among mammalian tissues and especially abundant in the ECM. SLRPs are divided into five canonical classes (classes I-III, containing fourteen members) and non-canonical classes (classes IV-V, including five members) based on their amino-acid structural sequence, chromosomal organization, and functional properties. Variations in both the protein core structure and glycosylation status lead to SLRP-specific interactions with cell membrane receptors, cytokines, growth factors, and structural ECM molecules. SLRPs have been implicated in the regulation of cancer growth, motility, and invasion, as well as in cancer-associated inflammation and autophagy, highlighting their crucial role in the processes of carcinogenesis. Except for the class I SLRP decorin, to which an anti-tumorigenic role has been attributed, other SLPRs' roles have not been fully clarified. This review will focus on the functions of the class I and II SLRP members biglycan and lumican, which are correlated to various aspects of cancer development.

Identification of hub genes and pathways in lung metastatic colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is one of the most prevalent types of malignant tumours. Metastasis is the leading cause of cancer-related mortality, with lung metastases accounting for 32.9% of all metastatic CRCs. However, since the biological mechanism of lung metastatic CRC is poorly understood, limited therapeutic targets are available. In the present study, we aimed to identify the key genes and molecular processes involved in CRC lung metastasis.

METHODS

The differentially expressed genes (DEGs) between primary and lung metastatic CRC patients were obtained from the Gene Expression Omnibus (GEO) database via the GEO2R tool. The enriched biological processes and pathways modulated by the DEGs were determined with Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Reactome Gene Sets analyses. The search tool Retrieval of Interacting Genes (STRING) and Cytoscape were used to construct a protein-protein interaction (PPI) network among DEGs.

RESULTS

The DEGs were enriched in surfactant metabolism, cell-cell communication and chemokine signaling pathways. The defined hub genes were included CLU, SFTPD, CCL18, SPP1, APOE, BGN and MMP3. Among them, CLU, SFTPD and CCL18 might be associated with the specific lung tropism metastasis in CRC. In addition, the expression and prognostic values of the hub genes in CRC patients were verified in database of The Cancer Genome Atlas (TCGA) and GEO. Moreover, the protein levels of the hub genes were detected in primary and lung metastatic CRC cells, serum or tissues. Furthermore, SFTPD was confirmed to facilitate cellular proliferation and lung metastasis in CRC.

CONCLUSION

This bioinformatics study may provide a better understanding of the candidate therapeutic targets and molecular mechanisms for CRC lung metastasis.

Unusual Suspects: Bone and Cartilage ECM Proteins as Carcinoma Facilitators

The extracellular matrix (ECM) is the complex three-dimensional network of fibrous proteins and proteoglycans that constitutes an essential part of every tissue to provide support for normal tissue homeostasis. Tissue specificity of the ECM in its topology and structure supports unique biochemical and mechanical properties of each organ. Cancers, like normal tissues, require the ECM to maintain multiple processes governing tumor development, progression and spread. A large body of experimental and clinical evidence has now accumulated to demonstrate essential roles of numerous ECM components in all cancer types. Latest findings also suggest that multiple tumor types express, and use to their advantage, atypical ECM components that are not found in the cancer tissue of origin. However, the understanding of cancer-specific expression patterns of these ECM proteins and their exact roles in selected tumor types is still sketchy. In this review, we summarize the latest data on the aberrant expression of bone and cartilage ECM proteins in epithelial cancers and their specific functions in the pathogenesis of carcinomas and discuss future directions in exploring the utility of this selective group of ECM components as future drug targets.