Profiling of Tumor Microenvironment Components Identifies Five Stroma-Related Genes with Prognostic Implications in Colorectal Cancer

Machine learning-based model for CD4+ conventional T cell genes to predict survival and immune responses in colorectal cancer

Globally, CRC ranks as a principal cause of mortality, with projections indicating a substantial rise in both incidence and mortality by the year 2040. The immunological responses to cancer heavily rely on the function of CD4Tconv. Despite this critical role, prognostic studies on CRC-related CD4Tconv remain insufficient. In this investigation, transcriptomic and clinical data were sourced from TCGA and GEO. Initially, we pinpointed CD4TGs using single-cell datasets. Prognostic genes were then isolated through univariate Cox regression analysis. Building upon this, 101 machine learning algorithms were employed to devise a novel risk assessment framework, which underwent rigorous validation using Kaplan-Meier survival analysis, univariate and multivariate Cox regression, time-dependent ROC curves, nomograms, and calibration plots. Furthermore, GSEA facilitated the examination of these genes' potential roles. The RS derived from this model was also analyzed for its implications in the TME, and its potential utility in immunotherapy and chemotherapy contexts. A novel prognostic model was developed, utilizing eight CD4TGs that are significantly linked to the outcomes of patients with CRC. This model's RS showcased remarkable predictive reliability for the overall survival rates of CRC patients and strongly correlated with malignancy levels. RS serves as an autonomous prognostic indicator, capable of accurately forecasting patient prognoses. Based on the median value of RS, patients were categorized into subgroups of high and low risk. The subgroup with higher risk demonstrated increased immune infiltration and heightened activity of genes associated with immunity. This investigation's establishment of a CD4TGs risk model introduces novel biomarkers for the clinical evaluation of CRC risks. These biomarkers may enhance therapeutic approaches and, in turn, elevate the clinical outcomes for patients with CRC by facilitating an integrated treatment strategy.

An epithelial-mesenchymal transition-related mRNA signature associated with the prognosis, immune infiltration and therapeutic response of colon adenocarcinoma

Background: Epithelial-mesenchymal transition (EMT) is closely associated with cancer cell metastasis. Colon adenocarcinoma (COAD) is one of the most common malignancies in the world, and its metastasis leading to poor prognosis remains a challenge for clinicians. The purpose of this study was to explore the prognostic value of EMT-related genes (EMTRGs) by bioinformatics analysis and to develop a new EMTRGs prognostic signature for COAD. Methods: The TCGA-COAD dataset was downloaded from the TCGA portal as the training cohort, and the GSE17538 and GSE29621 datasets were obtained from the GEO database as the validation cohort. The best EMTRGs prognostic signature was constructed by differential expression analysis, Cox, and LASSO regression analysis. Gene set enrichment analysis (GSEA) is used to reveal pathways that are enriched in high-risk and low-risk groups. Differences in tumor immune cell levels were analyzed using microenvironmental cell population counter and single sample gene set enrichment analysis. Subclass mapping analysis and Genomics of Drug Sensitivity in Cancer were applied for prediction of immunotherapy response and chemotherapy response, respectively. Results: A total of 77 differentially expressed EMTRGs were identified in the TCGA-COAD cohort, and they were significantly associated with functions and pathways related to cancer cell metastasis, proliferation, and apoptosis. We constructed EMTRGs prognostic signature with COMP, MYL9, PCOLCE2, SCG2, and TIMP1 as new COAD prognostic biomarkers. The high-risk group had a poorer prognosis with enhanced immune cell infiltration. The GSEA demonstrated that the high-risk group was involved in "ECM Receptor Interaction," "WNT Signaling Pathway" and "Colorectal Cancer." Furthermore, patients with high risk scores may respond to anti-CTLA4 therapy and may be more resistant to targeted therapy agents BI 2536 and ABT-888. Conclusion: Together, we developed a new EMTRGs prognostic signature that can be an independent prognostic factor for COAD. This study has guiding implications for individualized counseling and treatment of COAD patients.

Identification of immune subtypes to guide immunotherapy and targeted therapy in clear cell renal cell carcinoma

Accumulating pieces of evidence suggested that immunotypes may indicate the overall immune landscape in the tumor microenvironment, which were closely related to therapeutic response. The purpose of this study was to classify and define the immune subtypes of clear cell renal cell carcinoma (ccRCC), so as to authenticate the potential immune subtypes that respond to immunotherapy. Transcriptome expression profile and mutation profile data of ccRCC, as well as clinical characteristics used in this study were obtained from TCGA database. There were significant differences in the infiltration of immune cells, immune checkpoints, and antigens between ccRCC and para-cancerous tissues. According to immune components, patients with ccRCC were divided into three immune subtypes, with different clinical and molecular characteristics. Compared with other subtypes, IS2 showed cold immune phenotype, and was associated with better survival. IS1 represented complex immune populations and was associated with poor overall survival (OS) and progression free survival (PFS). Further analysis indicated that expression of immune checkpoints also differed among the three subtypes, and was abnormally up-regulated in IS3. Pathway enrichment analysis indicated that the mTOR signaling pathway was abnormally enriched in IS3, while the TGF_BETA, ANGIOGENESIS and receptor tyrosine kinase signaling pathways were abnormally enriched in IS2. Furthermore, there was an abnormal enrichment of the epithelial-to-mesenchymal transition (EMT) signaling pathway in IS1, which may be associated with a higher rate of metastasis. Finally, SCG2 was screened as a specific antigen of ccRCC, which was not only related to poor prognosis, but also significantly associated with immune cells and immune checkpoints. In conclusion, the immune subtypes of ccRCC may provide new insights into the tumor biology and the precise clinical management of this disease.

Classification of colon adenocarcinoma based on immunological characterizations: Implications for prognosis and immunotherapy

Accurate immune molecular typing is pivotal for screening out patients with colon adenocarcinoma (COAD) who may benefit from immunotherapy and whose tumor microenvironment (TME) was needed for reprogramming to beneficial immune-mediated responses. However, little is known about the immune characteristic of COAD. Here, by calculating the enrichment score of immune characteristics in three online COAD datasets (TCGA-COAD, GSE39582, and GSE17538), we identified 17 prognostic-related immune characteristics that overlapped in at least two datasets. We determined that COADs could be stratified into three immune subtypes (IS1–IS3), based on consensus clustering of these 17 immune characteristics. Each of the three ISs was associated with distinct clinicopathological characteristics, genetic aberrations, tumor-infiltrating immune cell composition, immunophenotyping (immune “hot” and immune “cold”), and cytokine profiles, as well as different clinical outcomes and immunotherapy/therapeutic response. Patients with the IS1 tumor had high immune infiltration but immunosuppressive phenotype, IS3 tumor is an immune “hot” phenotype, whereas those with the IS2 tumor had an immune “cold” phenotype. We further verified the distinct immune phenotype of IS1 and IS3 by an in-house COAD cohort. We propose that the immune subtyping can be utilized to identify COAD patients who will be affected by the tumor immune microenvironment. Furthermore, the ISs may provide a guide for personalized cancer immunotherapy and for tumor prognosis.

Identification of Angiogenesis-Related Prognostic Biomarkers Associated With Immune Cell Infiltration in Breast Cancer

Background: This study aimed to explore the prognostic value of angiogenesis-related genes (ARGs) and their association with immune cell infiltration (ICI) in breast cancer (BC).

Methods: Transcriptome data of BC were obtained from the TCGA and GEO databases. Differentially expressed ARGs were identified by the limma package. The identification of key genes and construction of the risk score model were performed by univariate and multivariate Cox regression algorithms. The prognostic value of the risk score was assessed by ROC curves and nomogram. GO, KEGG pathway, and GSEA were used to investigate the biological functions of differentially expressed genes (DEGs), and CIBERSORT, ssGSEA, and xCell algorithms were performed to estimate the ICI in high-risk and low-risk groups. The correlations between prognostic biomarkers and differentially distributed immune cells were assessed. Moreover, a ceRNA regulatory network based on prognostic biomarkers was constructed and visualized by Cytoscape software.

Results: A total of 18 differentially expressed ARGs were identified between tumor and adjacent normal tissue samples. TNFSF12, SCG2, COL4A3, and TNNI3 were identified as key prognostic genes by univariate and multivariate Cox regression analyses. The risk score model was further constructed based on the four-gene signature and validated in GSE7390 and GSE88770 datasets. ROC curves and nomogram indicated that the risk score had good accuracy for determining BC patient survival. Biological function analysis showed that DEGs in high- and low-risk groups had a high enrichment in immune-related biological processes and signaling pathways. Moreover, significantly different ICIs were found between high- and low-risk groups, such as memory B cells, CD8⁺ T cells, resting memory CD4⁺ T cells, follicular helper T cells, regulatory T cells, monocytes, M2 macrophages, and neutrophils, and each prognostic biomarker was significantly correlated with one or more immune cell types.

Conclusion: The current study identified novel prognostic ARGs and developed a prognostic model for predicting survival in patients with BC. Furthermore, this study indicated that ICI may act as a bond between angiogenesis and BC. These findings enhance our understanding of angiogenesis in BC and provide novel guidance on developing therapeutic targets for BC patients.

SCG2 is a Prognostic Biomarker Associated With Immune Infiltration and Macrophage Polarization in Colorectal Cancer

Colorectal cancer (CRC) is the second most lethal malignancy around the world. Limited efficacy of immunotherapy creates an urgent need for development of novel treatment targets. Secretogranin II (SCG2) is a member of the chromogranin family of acidic secretory proteins, has a role in tumor microenvironment (TME) of lung adenocarcinoma and bladder cancer. Besides, SCG2 is a stroma-related gene in CRC, its potential function in regulating tumor immune infiltration of CRC needs to be fully elucidated. In this study, we used western blot, real-time PCR, immunofluorescence and public databases to evaluate SCG2 expression levels and distribution. Survival analysis and functional enrichment analysis were performed. We examined TME and tumor infiltrating immune cells using ESTIMATE and CIBERSORT algorithm. The results showed that SCG2 expression was significantly decreased in CRC tumor tissues, and differentially distributed between tumor and adjacent normal tissues. SCG2 was an independent prognostic predictor in CRC. High expression of SCG2 correlated with poor survival and advanced clinical stage in CRC patients. SCG2 might regulate multiple tumor- and immune-related pathways in CRC, influence tumor immunity by regulating infiltration of immune cells and macrophage polarization in CRC.

Circ_0025039 acts an oncogenic role in the progression of non-small cell lung cancer through miR-636-dependent regulation of CORO1C

Non-small cell lung cancer remains the leading cause of cancer-related death worldwide. Circular RNA plays vital roles in NSCLC progression. This study is designed to reveal the role of circ_0025039 in NSCLC cell malignancy. The RNA expression of circ_0025039, microRNA-636 (miR-636), and coronin 1C was detected by quantitative real-time polymerase chain reaction. Protein expression was checked by Western blot analysis or immunohistochemistry assay. Cell proliferation, migration, invasion, tube formation ability, sphere formation capacity, and apoptosis were investigated by cell counting kit-8, 5-Ethynyl-29-deoxyuridine, transwell assay, tube formation assay, sphere formation assay, and flow cytometry analysis, respectively. Mouse model assay was conducted to reveal the effect of circ_0025039 silencing on tumor formation in vivo. The interaction between miR-636 and circ_0025039 or CORO1C was identified through dual-luciferase reporter and RNA pull-down assays. The expression of circ_0025039 and CORO1C was significantly increased, while miR-636 was decreased in NSCLC tissues and cells compared with controls. Circ_0025039 depletion repressed NSCLC cell proliferation, migration, invasion, tube-forming capacity, and sphere formation ability, but induced cell apoptosis. The neoplasm formation was repressed after circ_0025039 silencing. Additionally, circ_0025039 acted as a sponge for miR-636, which was found to target CORO1C. Importantly, the contribution of circ_0025039 to NSCLC progression was mediated by miR-636/CORO1C axis. Circ_0025039 silencing repressed NSCLC malignant progression by reducing CORO1C expression through miR-636, showing the possibility of circ_0025039 as a therapeutic target for NSCLC.

Comparison of Gene Expression Profiles of Signet Ring Cell Carcinoma and Poorly Cohesive Carcinoma in Early Gastric Cancer

Introduction:

Although signet ring cell carcinoma (SRC) is a subtype of poorly cohesive carcinoma (PC), the differences in the biological behavior between the 2 morphologically similar carcinomas have not been fully studied. Therefore, we performed transcriptome analysis to evaluate the differences of genetic expressions between SRC and PC.

Methods:

The study group consisted of patients with SRC or PC pathology from patients with early gastric cancer (EGC) whose depth of invasion was localized in the mucosal layer. A total of 18 patients were enrolled. The patients were divided into 3 groups based on their histologic type and lymph node (LN) status. Group 1 consisted of patients with PC and positive LN metastasis, Group 2 consisted of patients with PC without LN metastasis, and Group 3 consisted of patients with SRC without LN metastasis. Transcriptome analysis was performed using the nCounter PanCancer Progression Panel Kit.

Results:

The expression of 77 genes in Group 1 was altered compared to that in normal tissues. The expression of 49 and 13 genes in Groups 2 and 3, respectively, was altered when compared to that in normal tissues. Groups 1 and 2 showed similar genetic expressions. However, Group 3 showed numerous differences in gene expression including Roundabout4 (Robo4) compared to the other groups, especially Group 1.

Conclusion:

Our data suggest that gene expression patterns were different between SRC and PC and expression of ROBO4 may play an important role in the prognosis of SRC and PC type of EGC.

Transcriptomics based multi-dimensional characterization and drug screen in esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) remains one of the deadly cancer types. Comprehensively dissecting the molecular characterization and the heterogeneity of ESCC paves the way for developing more promising therapeutics.

METHODS

Expression profiles of multiple ESCC datasets were integrated. ATAC-seq and RNA-seq were combined to reveal the chromatin accessibility features. A prognosis-related subtype classifier (PrSC) was constructed, and its association with the tumor microenvironment (TME) and immunotherapy was assessed. The key gene signature was validated in clinical samples. Based on the TME heterogeneity of ESCC patients, potential subtype-specific therapeutic agents were screened.

FINDINGS

The common differentially expressed genes (cDEGs) in ESCC were identified. Up-regulated genes (HEATR1, TIMELESS, DTL, GINS1, RUVBL1, and ECT2) were found highly important in ESCC cell survival. The expression alterations of PRIM2, HPGD, NELL2, and TFAP2B were associated with chromatin accessibility changes. PrSC was a robust scoring tool that was not only associated with the prognosis of ESCC patients, but also could reflect the TME heterogeneity. TNS1^high fibroblasts were associated with immune exclusion. TG-101348 and Vinorelbine were identified as potential subtype-specific therapeutic agents. Besides, the application of PrSC into two immunotherapy cohorts indicated its potential value in assessing treatment response to immunotherapy.

INTERPRETATION

Our study depicted the multi-dimensional characterization of ESCC, established a robust scoring tool for the prognosis assessment, highlighted the role of TNS1^high fibroblasts in TME, and identified potential drugs for clinical use.

FUNDING

A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

Secretogranin II impairs tumor growth and angiogenesis by promoting degradation of hypoxia-inducible factor-1α in colorectal cancer

Distant metastasis is a major cause of death in patients with colorectal cancer (CRC) but the management of advanced and metastatic CRC still remains problematic due to the distinct molecular alterations during tumor progression. Tumor angiogenesis is a key step in tumor growth, invasion and metastasis. However, the signaling pathways involved in angiogenesis are poorly understood. The results of the present study showed that secretogranin II (SCG2) was significantly downregulated in malignant CRC tissues, and higher expression of SCG2 was correlated with longer disease‐free survival and overall survival of CRC patients. The results of an animal study showed that ectopic expression of SCG2 significantly inhibited CRC tumor growth by disrupting angiogenesis. Furthermore, the inhibition of expression of vascular endothelial growth factor (VEGF) by SCG2 and rescue of VEGF effectively blocked SCG2‐induced inhibition of angiogenesis. Investigations into the underlying mechanism suggested that SCG2 promoted degradation of hypoxia‐inducible factor (HIF)‐1α by interacting with the von Hippel–Lindau tumor suppressor in CRC cells. Blocking of degradation of HIF‐1α effectively attenuated the SCG2‐mediated decrease in expression of VEGF in CRC cells. Collectively, these results demonstrated that treatment with SCG2 effectively inhibited CRC tumor growth by disrupting the activities of HIF‐1α/VEGF, thereby clarifying the anti‐tumor and anti‐angiogenesis roles of SCG2 in CRC, while providing a novel therapeutic target and a potential prognostic marker of disease progression.

Weighted Gene Co-expression Network Analysis Identifies CALD1 as a Biomarker Related to M2 Macrophages Infiltration in Stage III and IV Mismatch Repair-Proficient Colorectal Carcinoma

Immunotherapy has achieved efficacy for advanced colorectal cancer (CRC) patients with a mismatch-repair-deficient (dMMR) subtype. However, little immunotherapy efficacy was observed in patients with the mismatch repair-proficient (pMMR) subtype, and hence, identifying new immune therapeutic targets is imperative for those patients. In this study, transcriptome data of stage III/IV CRC patients were retrieved from the Gene Expression Omnibus database. The CIBERSORT algorithm was used to quantify immune cellular compositions, and the results revealed that M2 macrophage fractions were higher in pMMR patients as compared with those with the dMMR subtype; moreover, pMMR patients with higher M2 macrophage fractions experienced shorter overall survival (OS). Subsequently, weighted gene co-expression network analysis and protein–protein interaction network analysis identified six hub genes related to M2 macrophage infiltrations in pMMR CRC patients: CALD1, COL6A1, COL1A2, TIMP3, DCN, and SPARC. Univariate and multivariate Cox regression analyses then determined CALD1 as the independent prognostic biomarker for OS. CALD1 was upregulated specifically the in CMS4 CRC subtype, and single-sample Gene Set Enrichment Analysis (ssGSEA) revealed that CALD1 was significantly correlated with angiogenesis and TGF-β signaling gene sets enrichment scores in stage III/IV pMMR CRC samples. The Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm and correlation analysis revealed that CALD1 was significantly associated with multiple immune and stromal components in a tumor microenvironment. In addition, GSEA demonstrated that high expression of CALD1 was significantly correlated with antigen processing and presentation, chemokine signaling, leukocyte transendothelial migration, vascular smooth muscle contraction, cytokine–cytokine receptor interaction, cell adhesion molecules, focal adhesion, MAPK, and TGF-beta signaling pathways. Furthermore, the proliferation, invasion, and migration abilities of cancer cells were suppressed after reducing CALD1 expression in CRC cell lines. Taken together, multiple bioinformatics analyses and cell-level assays demonstrated that CALD1 could serve as a prognostic biomarker and a prospective therapeutic target for stage III/IV pMMR CRCs.

Knockdown of hsa_circ_0037658 inhibits the progression of osteoarthritis via inducing autophagy

Osteoarthritis (OA) is a chronic musculoskeletal degeneration disease that can result in chronic pain and functional disability. Circular RNAs (CirRNAs) are known to be involved in OA. It was reported that hsa_circ_0037658 was notably upregulated in OA tissues; however, the biological role of hsa_circ_0037658 in OA remains unclear. To investigate the function of hsa_circ_0037658 in OA, CHON-001 cells were treated with IL-1β. The effect of hsa_circ_0037658 knockdown on cell growth was tested by CCK-8 and immunofluorescence staining. In addition, the correlation between hsa_circ_0037658 and autophagy was explored by LC3 staining and western blot. The results indicated that hsa_circ_0037658 was significantly upregulated in IL-1β-treated CHON-001 cells. The silencing of hsa_circ_0037658 could protect CHON-001 cell injury against IL-1β. Moreover, hsa_circ_0037658 shRNA reversed IL-1β-induced cell growth inhibition via inducing cell autophagy. Furthermore, knockdown of hsa_circ_0037658 notably alleviated the symptom of OA in vivo. To sum up, knockdown of hsa_circ_0037658 suppressed the progression of OA via inducing autophagy. Thus, hsa_circ_0037658 might serve as a potential target for the treatment of OA.