ECM-Receptor Regulatory Network and Its Prognostic Role in Colorectal Cancer

Comparative transcriptomic study on the ovarian cancer between chicken and human

The laying hen is the spontaneous model of ovarian tumor. A comprehensive comparison based on RNA-seq from hens and women may shed light on the molecular mechanisms of ovarian cancer. We performed next-generation sequencing of microRNA and mRNA expression profiles in 9 chicken ovarian cancers and 4 normal ovaries, which has been deposited in GSE246604. Together with 6 public datasets (GSE21706, GSE40376, GSE18520, GSE27651, GSE66957, TCGA-OV), we conducted a comparative transcriptomics study between chicken and human. In the present study, miR-451, miR-2188-5p, and miR-10b-5p were differentially expressed in normal ovaries, early- and late-stage ovarian cancers. We also disclosed 499 up-regulated genes and 1,061 down-regulated genes in chicken ovarian cancer. The molecular signals from 9 cancer hallmarks, 25 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and 369 Gene Ontology (GO) pathways exhibited abnormalities in ovarian cancer compared to normal ovaries via Gene Set Enrichment Analysis (GSEA). In the comparative analysis across species, we have uncovered the conservation of 5 KEGG and 76 GO pathways between chicken and human including the mismatch repair and ECM receptor interaction pathways. Moreover, a total of 174 genes contributed to the core enrichment for these KEGG and GO pathways were identified. Among these genes, the 22 genes were found to be associated with overall survival in patients with ovarian cancer. In general, we revealed the microRNA profiles of ovarian cancers in hens and updated the mRNA profiles previously derived from microarrays. And we also disclosed the molecular pathways and core genes of ovarian cancer shared between hens and women, which informs model animal studies and gene-targeted drug development.

Ginsenosides: an immunomodulator for the treatment of colorectal cancer

Ginsenosides, the primary bioactive ingredients derived from the root of Panax ginseng, are eagerly in demand for tumor patients as a complementary and alternative drug. Ginsenosides have increasingly become a "hot topic" in recent years due to their multifunctional role in treating colorectal cancer (CRC) and regulating tumor microenvironment (TME). Emerging experimental research on ginsenosides in the treatment and immune regulation of CRC has been published, while no review sums up its specific role in the CRC microenvironment. Therefore, this paper systematically introduces how ginsenosides affect the TME, specifically by enhancing immune response, inhibiting the activation of stromal cells, and altering the hallmarks of CRC cells. In addition, we discuss their impact on the physicochemical properties of the tumor microenvironment. Furthermore, we discuss the application of ginsenosides in clinical treatment as their efficacy in enhancing tumor patient immunity and prolonging survival. The future perspectives of ginsenoside as a complementary and alternative drug of CRC are also provided. This review hopes to open up a new horizon for the cancer treatment of Traditional Chinese Medicine monomers.

ADAMTSL2 is a potential prognostic biomarker and immunotherapeutic target for colorectal cancer: Bioinformatic analysis and experimental verification

The ADAMTS Like 2 (ADAMTSL2) mutation has been identified to be associated with different human genetic diseases. The role of ADAMTSL2 is unclear in colorectal cancer (CRC). The study investigated the expression of ADAMTSL2 in both pan cancer and CRC, using data from The Cancer Genome Atlas (TCGA) database to assess its diagnostic value. The study examined the correlation between ADAMTSL2 expression levels and clinical characteristics, as well as prognosis in CRC. The study explored potential regulatory networks involving ADAMTSL2, including its association with immune infiltration, immune checkpoint genes, tumor mutational burden (TMB) / microsatellite instability (MSI), tumor stemness index (mRNAsi), and drug sensitivity in CRC. ADAMTSL2 expression was validated using GSE71187 and quantitative real-time PCR (qRT-PCR). ADAMTSL2 was aberrantly expressed in pan cancer and CRC. An increased level of ADAMTSL2 expression in patients with CRC was significantly associated with the pathologic N stage (p < 0.001), pathologic stage (p < 0.001), age (p < 0.001), histological type (p < 0.001), and neoplasm type (p = 0.001). The high expression of ADAMTSL2 in patients with CRC was found to be significantly associated with a poorer overall survival (OS) (HR: 1.67; 95% CI: 1.18-2.38; p = 0.004), progression-free survival (PFS) (HR: 1.55; 95% CI: 1.14-2.11; p = 0.005) and disease-specific survival (DSS) (HR: 1.83; 95% CI: 1.16-2.89; p = 0.010). The expression of ADAMTSL2 in patients with CRC (p = 0.009) was identified as an independent prognostic determinant. ADAMTSL2 was associated with extracellular matrix receptor (ECM-receptor) interaction, transforming growth factor β (TGF-β) signaling pathway, and more. ADAMTSL2 expression was correlated with immune infiltration, immune checkpoint genes, TMB / MSI and mRNAsi in CRC. ADAMTSL2 expression was significantly and negatively correlated with 1-BET-762, Trametinib, and WZ3105 in CRC. ADAMTSL2 was significantly upregulated in CRC cell lines. The high expression of ADAMTSL2 is significantly correlated with lower OS and immune infiltration of CRC. ADAMTSL2 may be a potential prognostic biomarker and immunotherapeutic target for CRC patients.

Novel DNA methylation-based epigenetic signatures in colorectal cancer from peripheral blood leukocytes

Colorectal cancer (CRC) is a multifactorial disease characterized by accumulation of multiple genetic and epigenetic alterations, transforming colonic epithelial cells into adenocarcinomas. Alteration of DNA methylation (DNAm) is a promising biomarker for predicting cancer risk and prognosis, but its role in CRC tumorigenesis is inconclusive. Notably, few DNAm studies have used pre-diagnostic peripheral blood (PB) DNA, causing difficulty in postulating the underlying biologic mechanism of CRC initiation. We conducted epigenome-wide association (EWA) scans in postmenopausal women from Women's Health Initiative (WHI) with their pre-diagnostic DNAm in PB leukocytes (PBLs) to prospectively evaluate CRC development. Our site-specific DNAm analyses across the genome adjusted for DNAm-age, leukocyte heterogeneities, as well as body mass index, diabetes, and insulin resistance. We validated 20 top EWA-CpGs in 2 independent CRC tissue datasets. Also, we detected differentially methylated regions (DMRs) associated with CRC, further mapped to transcriptomic profile, and finally conducted a Gene Set Enrichment Analysis. We detected multiple novel CpGs validated across WHI and tissue datasets. In particular, 2 CpGs (B4GALNT4cg10321339, SV2Bcg18144285) had the strongest effect on CRC risk. Results from our DMR scans contained MIR663cg06007966, which was also validated in EWA analyses. Also, we detected 1 methylome region in PEG10 of Chr7 shared across datasets. Our findings reflect both novel and well-established epigenomic and transcriptomic sites in CRC, warranting further functional validations. Our study contributes to better understanding of the complex interrelated mechanisms on the methylome underlying CRC tumorigenesis and suggests novel preventive DNAm-targets in PBLs for detecting at-risk individuals for CRC development.

Colorectal cancer subtyping and immune landscape analysis based on natural killer cell-related genes

BACKGROUND AND STUDY AIMS

The prognosis of colorectal cancer (CRC) is related to natural killer (NK) cells, but the molecular subtype features of CRC based on NK cells are still unknown. This study aimed to identify NK cell-related molecular subtypes of CRC and analyze the survival status and immune landscape of patients with different subtypes.

PATIENTS/MATERIAL AND METHODS

mRNA expression data, single nucleotide variant (SNV) data, and clinical information of CRC patients were obtained from The Cancer Genome Atlas. Differentially expressed genes (DEGs) were obtained through differential analysis, and the intersection was taken with NK cell-associated genes to obtain 103 NK cell-associated CRC DEGs (NCDEGs). Based on NCDEGs, CRC samples were divided into three clusters through unsupervised clustering analysis. Survival analysis, immune analysis, Gene Set Enrichment Analysis (GSEA), and tumor mutation burden (TMB) analysis were performed. Finally, NCDEG-related small-molecule drugs were screened using the CMap database.

RESULTS

Survival analysis revealed that cluster2 had a lower survival rate than cluster1 and cluster3 (p < 0.05). Immune infiltration analysis found that the immune infiltration levels and immune checkpoint expression levels of cluster1_3 were substantially higher than those of cluster2, and the tumor purity was the opposite (p < 0.05). GSEA presented that cluster1_3 was significantly enriched in the chemokine signaling pathway, ECM receptor interaction, and antigen processing and presentation pathways (p < 0.05). The TMB of cluster1_3 was significantly higher than that of cluster2 (p < 0.05). Genes with the highest mutation rate in CRC were APC, TP53, TTN, and KRAS. Drug prediction results showed that small-molecule drugs that reverse the upregulation of NCDEGs, deoxycholic acid, dipivefrine, phenformin, and other drugs may improve the prognosis of CRC.

CONCLUSION

NK cell-associated CRC subtypes can be used to evaluate the tumor characteristics of CRC patients and provide an important reference for CRC patients.

Genetic structural analysis of different breeds and geographical groups of Fenneropenaeus chinensis reveals population diversity

Fenneropenaeus chinensis is a commercially important shrimp species cultured in China. This study investigated eight F. chinensis populations in China, including four geographical populations, three commercial breeds, and one wild population captured from the Yellow Sea. Population stratification analysis revealed that the Hebei geographical population and commercial breeding "Huanghai No. 4" were relatively independent and stable, reflecting a relatively closed breeding environment, whereas gene introgression was present between other populations. Selective signature analysis detected artificial selection for vision, growth, and disease resistance in the Hebei population. Neuronal development-related genes were detected to be under selection in the Changyi and Rizhao populations. Fertility of the Rizhao population was also investigated. Additionally, genes in the glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate pathway were involved in the high pH tolerance of the "Huanghai No. 4" population. This study provided support for the genetic mechanism of parsing economic traits and the development of molecular breeding technologies.

Integrating single-cell and spatial analysis reveals MUC1-mediated cellular crosstalk in mucinous colorectal adenocarcinoma

BACKGROUND

Mucinous colorectal adenocarcinoma (MCA) is a distinct subtype of colorectal cancer (CRC) with the most aggressive pattern, but effective treatment of MCA remains a challenge due to its vague pathological characteristics. An in-depth understanding of transcriptional dynamics at the cellular level is critical for developing specialised MCA treatment strategies.

METHODS

We integrated single-cell RNA sequencing and spatial transcriptomics data to systematically profile the MCA tumor microenvironment (TME), particularly the interactome of stromal and immune cells. In addition, a three-dimensional bioprinting technique, canonical ex vivo co-culture system, and immunofluorescence staining were further applied to validate the cellular communication networks within the TME.

RESULTS

This study identified the crucial intercellular interactions that engaged in MCA pathogenesis. We found the increased infiltration of FGF7+/THBS1+ myofibroblasts in MCA tissues with decreased expression of genes associated with leukocyte-mediated immunity and T cell activation, suggesting a crucial role of these cells in regulating the immunosuppressive TME. In addition, MS4A4A+ macrophages that exhibit M2-phenotype were enriched in the tumoral niche and high expression of MS4A4A+ was associated with poor prognosis in the cohort data. The ligand-receptor-based intercellular communication analysis revealed the tight interaction of MUC1+ malignant cells and ZEB1+ endothelial cells, providing mechanistic information for MCA angiogenesis and molecular targets for subsequent translational applications.

CONCLUSIONS

Our study provides novel insights into communications among tumour cells with stromal and immune cells that are significantly enriched in the TME during MCA progression, presenting potential prognostic biomarkers and therapeutic strategies for MCA.

KEY POINTS

Tumour microenvironment profiling of MCA is developed. MUC1+ tumour cells interplay with FGF7+/THBS1+ myofibroblasts to promote MCA development. MS4A4A+ macrophages exhibit M2 phenotype in MCA. ZEB1+ endotheliocytes engage in EndMT process in MCA.

Investigation of cellular communication and signaling pathways in tumor microenvironment for high TP53-expressing osteosarcoma cells through single-cell RNA sequencing

Osteosarcoma (OS) stands as the most prevalent primary bone cancer in children and adolescents, and its limited treatment options often result in unsatisfactory outcomes, particularly for metastatic cases. The tumor microenvironment (TME) has been recognized as a crucial determinant in OS progression. However, the intercellular dynamics between high TP53-expressing OS cells and neighboring cell types within the TME are yet to be thoroughly understood. In our study, we harnessed the single-cell RNA sequencing (scRNA-seq) technology in combination with the computational tool-Cellchat, aiming to elucidate the intercellular communication networks present within OS. Through meticulous quantitative inference and subsequent analysis of these networks, we succeeded in identifying significant signaling pathways connecting high TP53-expressing OS cells with proximate cell types, namely Macrophages, Monocytes, Endothelial Cells, and PVLs. This research brings forth a nuanced understanding of the intricate patterns and coordination involved in the TME's intercellular communication signals. These findings not only provide profound insights into the molecular mechanisms underpinning OS but also indicate potential therapeutic targets that could revolutionize treatment strategies.

The Prognosis and Immunotherapy Prediction Model of Ovarian Serous Cystadenocarcinoma Patient was Constructed Based on Cuproptosis-Related LncRNA

Cuproptosis can serve as potential prognostic predictors in patients with cancer. However, the role of this relationship in ovarian serous cystadenocarcinoma (OV) remains unclear. 376 OV tumor samples were obtained from the Cancer Genome Atlas (TCGA) database, and long non-coding RNAs (lncRNAs) related to cuproptosis were obtained through correlation analysis. The risk assessment model was further constructed by univariate Cox regression analysis and LASSO Cox regression. Bioinformatics was used to analyze the regulatory effect of relevant risk assessment models on tumor mutational burden (TMB) and immune microenvironment. We obtained 5 lncRNAs (AC025287.2, AC092718.4, AC112721.2, LINC00996, and LINC01639) and incorporated them into the Cox proportional hazards model. Kaplan-Meier (KM) curve analysis of the prognosis found that the high-risk group was associated with a poorer prognosis. The receiver operating characteristic (ROC) curve showed stronger predictive power compared to other clinicopathological features. Immune infiltration analysis showed that high-risk scores were inversely correlated with CD8+ T cells, CD4+ T cells, macrophages, NK cells, and B cells. Functional enrichment analysis found that they may act via the extracellular matrix (ECM)-interacting proteins and other pathways. We successfully constructed a reliable cuproptosis-related lncRNA model for the prognosis of OV.

Platycodin D ameliorates ammonia-induced pulmonary fibrosis by repressing TGF-β1-mediated extracellular matrix remodeling

Ammonia can induce pulmonary fibrosis in humans and animals. Platycodin D (PLD) possesses various bioactive activities including anti-fibrotic properties. In this study, we aimed to explore the activity and mechanism of PLD in pulmonary fibrosis induced by ammonia. The mouse model of ammonia-induced lung fibrosis was established, and the role of PLD was assessed by H&E and Masson's trichrome staining. The differentially expressed genes (DEGs) were identified by RNA-seq and subjected to GO and KEGG pathway analyses. BEAS-2B cells were treated with NH4 Cl alone or along with PLD. Results showed that PLD attenuated ammonia-induced pulmonary inflammation and fibrosis in vivo. The extracellular matrix (ECM)-receptor interaction pathway was predicted as a prominent pathway underlying the anti-fibrotic function of PLD. In ammonia-induced mouse models and NH4 Cl-treated BEAS-2B cells, PLD could repress the activation of the TGF-β1 pathway. By incubating lung fibroblast HFL1 cells with the conditioned medium of BEAS-2B cells treated with NH4Cl alone or along with PLD, PLD was confirmed to attenuate NH4 Cl-induced ECM deposition in HFL1 cells. Our findings demonstrate that PLD exerts a protective function in ammonia-induced pulmonary fibrosis by repressing TGF-β1-mediated ECM remodeling, suggesting the potential therapeutic value of PLD in this disease.

Identification of a 5-Gene Cuproptosis Signature Predicting the Prognosis for Colon Adenocarcinoma Based on WGCNA

Background: Colorectal cancer is a highly aggressive malignant tumor that primarily affects the digestive system. It is frequently diagnosed at an advanced stage. Cuproptosis is a copper-dependent form cell death mechanism, distinct from all other known pathways underlying cell death, tumor progression, prognosis, and immune response. Although the role of cuproptosis in colorectal cancer has been investigated over time, there is still an urgent need to explore new methods and insights to understand its potential function. Methods: The Gene Expression Omnibus and The Cancer Genome Atlas gene expression data were systematically explored to investigate the role of cuproptosis in colon adenocarcinoma. The weighted gene coexpression network analysis was used to construct a gene coexpression network and identify the critical module and cuproptosis-related genes correlated with colon adenocarcinoma prognosis. A cuproptosis-related genes prognostic signature for colon adenocarcinoma was identified and validated. To validate the identified gene signature, quantitative reverse transcription-polymerase chain reaction was performed. Cell proliferation assays were analyzed by CCK8 and cell cycle detection. In addition, reactive oxygen species assay was also analyzed. Results: Five hub cuproptosis-related genes (Dihydrolipoamide S-acetyltransferase, Cyclin-dependent kinase inhibitor 2A, ATOX1, VEGFA, and ULK1) were screened and a prognostic risk model for predicting overall survival was established based on these genes. The model was successfully tested in the validation cohort and the GEPIA database. Colon adenocarcinoma patients were categorized into high-risk and low-risk groups based on risk scores. The study revealed that patients with higher risk scores were more likely to have a poor prognosis. Moreover, Dihydrolipoamide S-acetyltransferase was a tumor suppressor gene that can induce cell death and affected the redox reactions in the colon cancer cell line. Conclusions: These findings suggest that the newly identified 5-gene signature may serve as a more reliable prognostic factor than clinical factors such as age and stage of disease. These findings offer a theoretical foundation for further investigation into potential cuproptosis-related biomarkers for predicting colon adenocarcinoma prognosis in the future.

A pan-cancer analysis of TNFAIP8L1 in human tumors

TNFAIP8L1, as a recently identified member in TNFAIP8 family, plays an important role in tumorigenesis. However, a pan-cancer analysis of TNFAIP8L1 in human tumors has not been conducted until now. The main purpose of study is to investigate TNFAIP8L1 during 33 different types of human tumors by using TCGA and GTEx. The pan-cancer analysis showed that TNFAIP8L1 was significantly over-expressed in 15 cancers and low-expressed in 9 cancers. There were distinct relations between TNFAIP8L1 expression and prognosis of patients with cancer. Furthermore, we also found that DNA methylation and RNA modification of TNFAIP8L1 were associated with many cancers. And then, we detected that TNFAIP8L1 level was positively associated with cancer-associated fibroblasts (CAFs) in many tumors. And, we obtained that TNFAIP8L1 expression was related with most of immune inhibitory and stimulatory genes in multiple types of tumors. We also found TNFAIP8L1 expression was correlated with most of chemokine, receptor, MHC, immunoinhibitor and immunostimulator gens in most of cancers. Moreover, we detected TNFAIP8L1 expression was associated with TMB and MSI in several tumors. Finally, TNFAIP8L1 gene had a significant positive association with 5 genes including BCL6B, DLL4, PCDH12, COL4A1 and DLL4 in the majority of tumors. GO enrichment and KEGG pathway analyses showed that TNFAIP8L1 in thepathogenesis of cancer may be related to "purine nucleoside binding," "purine ribonucleoside binding," "ECM-receptor interaction," etc. Our first pan-cancer study may provide a deep comprehending of TNFAIP8L1 in tumoeigenesis from different tumors.

Integrating scRNA-seq and bulk RNA-seq to characterize infiltrating cells in the colorectal cancer tumor microenvironment and construct molecular risk models

Colorectal cancer (CRC) is a malignancy that is both highly lethal and heterogeneous. Although the correlation between intra-tumoral genetic and functional heterogeneity and cancer clinical prognosis is well-established, the underlying mechanism in CRC remains inadequately understood. Utilizing scRNA-seq data from GEO database, we re-isolated distinct subsets of cells, constructed a CRC tumor-related cell differentiation trajectory, and conducted cell-cell communication analysis to investigate potential interactions across cell clusters. A prognostic model was built by integrating scRNA-seq results with TCGA bulk RNA-seq data through univariate, LASSO, and multivariate Cox regression analyses. Eleven distinct cell types were identified, with Epithelial cells, Fibroblasts, and Mast cells exhibiting significant differences between CRC and healthy controls. T cells were observed to engage in extensive interactions with other cell types. Utilizing the 741 signature genes, prognostic risk score model was constructed. Patients with high-risk scores exhibited a significant correlation with unfavorable survival outcomes, high-stage tumors, metastasis, and low responsiveness to chemotherapy. The model demonstrated a strong predictive performance across five validation cohorts. Our investigation involved an analysis of the cellular composition and interactions of infiltrates within the microenvironment, and we developed a prognostic model. This model provides valuable insights into the prognosis and therapeutic evaluation of CRC.

Heat shock protein family A member 8 is a prognostic marker for bladder cancer: Evidences based on experiments and machine learning

Heat shock protein member 8 (HSPA8) is one of the most abundant chaperones in eukaryotic cells, but its biological roles in bladder cancer (BC) are largely unclear. First, we observed that HSPA8 was abundant in both cell lines and tissues of BC, and the HSPA8-high group had poorer T stages and overall survival (OS) than the HSPA8-low group in the TCGA patients. Next, when we knocked down HSPA8 in BC cells, the growth and migration abilities were significantly decreased, the apoptosis rates were significantly increased, and the Ki67 fluorescence intensity was decreased in BC cells. Moreover, caspase 3 was significantly decreased with overexpression of HSPA8 in BC cells. After that, a machine learning prognostic model was created based on the expression of HSPA8 by applying LASSO Cox regression in TCGA and GEO patients. The model indicated that the low-risk (LR) group with BC had better tumour stages, lymphovascular invasion, and OS than the high-risk (HR) group. Additionally, the risk score was demonstrated to be an independent risk factor for the prognosis of BC by univariate and multivariate Cox analyses. Moreover, the HR group showed a greater rate of TP53 mutations and was mostly enriched in the ECM-receptor interaction pathway than the LR group. Importantly, lower CD8+ T-cell and NK cell infiltration, higher immune exclusion scores, higher expression of PD-L1 and CTLA4 and poorer immune checkpoint therapy effects were found in the HR group. These findings demonstrated how crucial HSPA8 plays a role in determining the prognosis of bladder cancer.

Synergistic Inhibitory Effect of Berberine and Low-Temperature Plasma on Non-Small-Cell Lung Cancer Cells via PI3K-AKT-Driven Signaling Axis

Low-temperature plasma (LTP) is an emerging biomedical technique that has been proposed as a potential approach for cancer therapy. Meanwhile, berberine (BER), an active ingredient extracted from various medical herbs, such as Coptischinesis, has been proven antitumor effects in a broad spectrum of cancer cells. In this study, we seek to develop a novel dual cancer therapeutic method by integrating pre-administration of BER and LTP exposure and evaluating its comprehensive antitumor effect on the human non-small-cell lung cancer (NSCLC) cell lines (A549 and H1299) in vitro. Cell viability, cell cycle, cell apoptosis, and intracellular and extracellular ROS were investigated. The results showed that cotreatment of BER and LTP significantly decreased the cell viability, arrested the cell cycle in the S phase, promoted cell apoptosis, and increased intracellular and extracellular ROS. Additionally, RNA Sequencing (RNA-Seq) technology was used to explore potential mechanisms. The differentially expressed genes among different treatment groups of NSCLC cells were analyzed and were mainly enriched in the phosphatidylinositol-3-kinase (PI3K)/protein kinase B (AKT) signaling pathway. Moreover, cotreatment of BER and LTP notably depressed the total protein expression level of PI3K and AKT with immunoblotting. In conclusion, BER and LTP have a synergistic inhibitory effect on NSCLC cells via the PI3K-AKT signaling pathway, which could provide a promising strategy for supplementary therapy in the anti-NSCLC battle.

The integrated single-cell analysis developed an immunogenic cell death signature to predict lung adenocarcinoma prognosis and immunotherapy

BACKGROUND

Research on immunogenic cell death (ICD) in lung adenocarcinoma (LUAD) has been relatively limited. This study aims to create ICD-related signatures for accurate survival prognosis prediction in LUAD patients, addressing the challenge of lacking reliable early prognostic indicators for this type of cancer.

METHODS

Using single-cell RNA sequencing (scRNA-seq) analysis, ICD activity in cells was calculated by AUCell algorithm, divided into high- and low-ICD groups according to median values, and key ICD regulatory genes were identified through differential analysis, and these genes were integrated into TCGA data to construct prognostic signatures using LASSO and COX regression analysis, and multi-dimensional analysis of ICD-related signatures in terms of prognosis, immunotherapy, tumor microenvironment (TME), and mutational landscape.

RESULTS

The constructed signature reveals a pronounced disparity in prognosis between the high- and low-risk groups of LUAD patients. The statistical discrepancies in survival times among LUAD patients from both the TCGA and GEO databases further corroborate this observation. Additionally, heightened levels of immune cell infiltration expression are evidenced in the low-risk group, suggesting a potential benefit from immunotherapeutic interventions for these patients. The expression levels of pivotal risk-associated genes in tissue samples were assessed utilizing qRT-PCR, thereby unveiling PITX3 as a plausible therapeutic target in the context of LUAD.

CONCLUSIONS

Our constructed ICD-related signatures provide help in predicting the prognosis and immunotherapy of LUAD patients, and to some extent guide the clinical treatment of LUAD patients.

Hyperuricemia and coronary heart disease: The mediating role of blood pressure and thrombospondin 3

BACKGROUND & AIMS

Although hyperuricemia is a known risk factor for coronary heart disease (CHD), little is known about the role of blood pressure in mediating this association. The purpose of this study is to investigate the role of blood pressure-related indicators and Thrombospondin 3 (THBS3) in the association between hyperuricemia and CHD.

METHODS AND RESULTS

Our observational epidemiology study included 593 CHD cases and 760 controls from a residential stable sample. We also chose 43 new CHD patients and 43 controls to test the expression levels of THBS3 using ELISA kits. We used logistic regression models and mediating effect analysis to investigate the relationships between hyperuricemia and CHD, as well as the mediating role of blood pressure-related indicators and THBS3. In the general population (OR: 2.001 [95% CI: 1.528-2.622]), male population (OR: 1.591 [95% CI: 1.119-2.262]), and female population (OR: 2.813 [95% CI: 1.836-4.310]), hyperuricemia is an independent risk factor for CHD. In general, average systolic blood pressure (SBP) and average pulse pressure difference (PPD) mediated 3.35% and 4.59%, respectively, of the association between hyperuricemia and CHD, and 6.60% and 6.60% in women. However, in the male population, we have not yet found that blood pressure-related indicators had a significant mediating effect. Meanwhile, we found that THBS3 mediated 19.23% of the association between hyperuricemia and CHD.

CONCLUSIONS

Average SBP, PPD, and THBS3 all play a role in the association of hyperuricemia and CHD. In the female population, similar mediating results in blood pressure-related indicators were observed.

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.

A Review of IsomiRs in Colorectal Cancer

As advancements in sequencing technology rapidly continue to develop, a new classification of microRNAs has occurred with the discovery of isomiRs, which are relatively common microRNAs with sequence variations compared to their established template microRNAs. This review article seeks to compile all known information about isomiRs in colorectal cancer (CRC), which has not, to our knowledge, been gathered previously to any great extent. A brief overview is given of the history of microRNAs, their implications in colon cancer, the canonical pathway of biogenesis and isomiR classification. This is followed by a comprehensive review of the literature that is available on microRNA isoforms in CRC. The information on isomiRs presented herein shows that isomiRs hold great promise for translation into new diagnostics and therapeutics in clinical medicine.

Exploring the prognostic function of TMB-related prognostic signature in patients with colon cancer

Tumor mutation burden (TMB) level is identified as a useful predictor in multiple tumors including colon adenocarcinoma (COAD). However, the function of TMB related genes has not been explored previously. In this study, we obtained patients' expression and clinical data from The Cancer Genome Atlas (TCGA) and the National Center for Biotechnology Information (NCBI). TMB genes were screened and subjected to differential expression analysis. Univariate Cox and LASSO analyses were utilized to construct the prognostic signature. The efficiency of the signature was tested by using a receiver operating characteristic (ROC) curve. A nomogram was further plotted to assess the overall survival (OS) time of patients with COAD. In addition, we compared the predictive performance of our signature with other four published signatures. Functional analyses indicated that patients in the low-risk group have obviously different enrichment of tumor related pathways and tumor infiltrating immune cells from that of high-risk patients. Our findings suggested that the ten genes' prognostic signature could exert undeniable prognostic functions in patients with COAD, which might provide significant clues for the development of personalized management of these patients.

The m6A/m5C/m1A regulator genes signature reveals the prognosis and is related with immune microenvironment for hepatocellular carcinoma

BACKGROUND

RNA methylation is a crucial in many biological functions, and its aberrant regulation is associated with cancer progression. N6-Methyladenosine (m6A), 5-Methylcytosine (m5C), N1-methyladenosine (m1A) are common modifications of RNA methylation. However, the effect of methylation of m6A/m5C/m1A in hepatocellular carcinoma (HCC) remains unclear.

METHOD

The transcriptome datasets, clinic information, and mutational data of 48 m6A/m5C/m1A regulator genes were acquired from the TCGA database, and the prognostic hazard model was established by univariate and Least absolute shrinkage and selection operator (Lasso) regression. The multivariate regression was performed to determine whether the risk score was an independent prognostic indicator. Kaplan-Meier survival analysis and ROC curve analysis were used to evaluate the predictive ability of the risk model. Decision curve analysis（DCA）analysis was conducted to estimate the clinical utility of the risk model. We further analyzed the association between risk score and functional enrichment, tumor immune microenvironment, and somatic mutation.

RESULT

The four-gene (YTHDF1, YBX1, TRMT10C, TRMT61A) risk signature was constructed. The high-risk group had shorter overall survival (OS) than the low-risk group. Univariate and multivariate regression analysis indicated that risk score was an independent prognostic indicator. Risk scores in male group, T3 + T4 group and Stage III + IV group were higher in female group, T1 + T2 group and stage I + II group. The AUC values for 1-, 2-, and 3-year OS in the TCGA dataset were 0.764, 0.693, and 0.689, respectively. DCA analysis showed that the risk score had a higher clinical net benefit in 1- and 2-year OS than other clinical features.The risk score was positively related to some immune cell infiltration and most immune checkpoints.

CONCLUSION

We developed a novel m6A/m5C/m1A regulator genes' prognostic model, which could be applied as a latent prognostic tool for HCC and might guide the choice of immunotherapies.

Leveraging various extracellular matrix levels to assess prognosis and sensitivity to immunotherapy in patients with ovarian cancer

Background

Ovarian cancer (OC) is the fifth leading cause of cancer-related deaths among women. Late diagnosis and heterogeneous treatment result in a poor prognosis for patients with OC. Therefore, we aimed to develop new biomarkers to predict accurate prognoses and provide references for individualized treatment strategies.

Methods

We constructed a co-expression network applying the "WGCNA" package and identified the extracellular matrix-associated gene modules. We figured out the best model and generated the extracellular matrix score (ECMS). The ECMS' ability to predict accurate OC patients' prognoses and responses to immunotherapy was evaluated.

Results

The ECMS was an independent prognostic factor in the training [hazard ratio (HR) = 3.132 (2.068-4.744), p< 0.001] and testing sets [HR = 5.514 (2.084-14.586), p< 0.001]. The receiver operating characteristic curve (ROC) analysis showed that the AUC values for 1, 3, and 5 years were 0.528, 0.594, and 0.67 for the training set, respectively, and 0.571, 0.635, and 0.684 for the testing set, respectively. It was found that the high ECMS group had shorter overall survival than the low ECMS group [HR = 2 (1.53-2.61), p< 0.001 in the training set; HR = 1.62 (1.06-2.47), p = 0.021 in the testing set; HR = 1.39 (1.05-1.86), p = 0.022 in the training set]. The ROC values of the ECMS model for predicting immune response were 0.566 (training set) and 0.572 (testing set). The response rate to immunotherapy was higher in patients with low ECMS.

Conclusion

We created an ECMS model to predict the prognosis and immunotherapeutic benefits in OC patients and provided references for individualized treatment of OC patients.

Microarray profiling identifies hsa_circ_0082003 as a novel tumor promoter for papillary thyroid carcinoma

BACKGROUND

Circular RNAs (circRNAs) are non-coding RNAs that have essential regulatory roles in the development of various tumors. This study explored whether circRNAs are involved in the progression of papillary thyroid carcinoma (PTC).

METHODS

Differentially expressed circRNAs (DECs) in four pairs of PTC and matched normal thyroid tissues were screened using a circRNA microarray. The potential functions of dysregulated circRNAs were predicted by bioinformatic analyses. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) was used to determine hsa_circ_0082003 expression in 80 pairs of PTC and matched normal thyroid tissues. Cell counting kit-8, colony formation, wound healing, and Transwell assays were performed to evaluate the biological functions of hsa_circ_0082003 in PTC cells. The role of hsa_circ_0082003 in PTC tumorigenesis in vivo was validated in nude mice.

RESULTS

In total, 3150 DECs (2317 upregulated and 833 downregulated) were identified. Pathway enrichment analyses indicated that the dysregulated circRNAs may play roles in PTC development. RT-qPCR validation demonstrated that hsa_circ_0082003 expression was significantly increased in PTC tissues and correlated with poor clinicopathological parameters. Receiver operating characteristic curve analysis showed that hsa_circ_0082003 had good performance for diagnosing PTC and judging whether it was accompanied by lymph node metastasis. Knockdown of hsa_circ_0082003 inhibited PTC cell proliferation, migration, and invasion. Tumor formation assays in vivo showed that downregulation of hsa_circ_0082003 significantly suppressed the growth of PTC.

CONCLUSION

Hsa_circ_0082003 may serve as a novel diagnostic biomarker and potential therapeutic target for PTC.

Prognostic risk analysis related to radioresistance genes in colorectal cancer

Background

Radiotherapy (RT) is one of the most important treatments for patients with colorectal cancer (CRC). Radioresistance is the crucial cause of poor therapeutic outcomes in colorectal cancer. However, the underlying mechanism of radioresistance in colorectal cancer is still poorly defined. Herein we established a radioresistant colorectal cancer cell line and performed transcriptomics analyses to search for the underlying genes that contribute to radioresistance and investigate its association with the prognosis of CRC patients.

Methods

The radioresistant cell line was developed from the parental HCT116 cell by a stepwise increased dose of irradiation. Differential gene analysis was performed using cellular transcriptome data to identify genes associated with radioresistance, from which extracellular matrix (ECM) and cell adhesion-related genes were screened. Survival data from a CRC cohort in the TCGA database were used for further model gene screening and validation. The correlation between the risk score model and tumor microenvironment, clinical phenotype, drug treatment sensitivity, and tumor mutation status were also investigated.

Results

A total of 493 different expression genes were identified from the radioresistant and wild-type cell line, of which 94 genes were associated with ECM and cell adhesion-related genes. The five model genes TNFRSF13C, CD36, ANGPTL4, LAMB3, and SERPINA1 were identified for CRC radioresistance via screening using the best model. A ROC curve indicated that the AUC of the resulting prognostic model (based on the 5-gene risk score and other clinical parameters, including age, sex, and tumor stages) was 0.79, 0.77, and 0.78 at 1, 2, and 3 years, respectively. The calibration curve showed high agreement between the risk score prediction and actual survival probability. The immune microenvironment, drug treatment sensitivity, and tumor mutation status significantly differed between the high- and low-risk groups.

Conclusions

The risk score model built with five radioresistance genes in this study, including TNFRSF13C, CD36, ANGPTL4, LAMB3, and SERPINA1, showed favorable performance in prognosis prediction after radiotherapy for CRC.

DNA methylation-mediated silencing of Neuronatin promotes hepatocellular carcinoma proliferation through the PI3K-Akt signaling pathway

AIMS

To explore the methylation status, function, and underlying mechanism of the imprinted gene Neuronatin (NNAT) in hepatocellular carcinoma (HCC) progression.

MAIN METHODS

Immunohistochemistry (IHC) was performed to evaluate the expression of NNAT in HCC samples. Bisulfite genomic sequencing PCR (BSP) was applied to examine the methylation status of the NNAT promoter. In addition, colony formation, 5-Ethynyl-20-deoxyuridine (EdU) assays and subcutaneous xenograft nude models were used to explore the roles of NNAT in HCC cell proliferation. Furthermore, RNA-seq and phospho-specific protein microarray assays were conducted to illustrate the underlying mechanism by which NNAT regulates HCC progression.

KEY FINDINGS

NNAT was obviously downregulated in HCC tissues, and its expression level was closely associated with tumor growth and patient prognosis. The downregulation of NNAT in HCC was induced by hypermethylation of CpG islands in the promoter region, and hypermethylation was correlated with overall survival of HCC. Moreover, the enforced expression of NNAT significantly inhibited HCC cell proliferation in vitro and in vivo. Transcriptome analysis showed that the alteration of NNAT expression was mainly related to dysregulation of the PI3K-Akt signaling pathway. Finally, phospho-specific antibody microarray detection further revealed that overexpressed NNAT can increase the phosphorylation levels of LKB1, Met, and elF4E and decrease the phosphorylation levels of PTEN, which are all involved in the PI3K-Akt signaling pathway.

SIGNIFICANCE

Our research provides new insights into the epigenetic regulation of imprinted genes in tumorigenesis and implies that the imprinted gene NNAT may act as a prognostic biomarker and tumor suppressor in HCC.

Gastrointestinal pan-cancer landscape of tumor matrix heterogeneity identifies biologically distinct matrix stiffness subtypes predicting prognosis and chemotherapy efficacy

Gastrointestinal (GI) cancers are a heterogeneous group of primary solid tumors, arising in GI tract from the esophagus to rectum. Matrix stiffness (MS) is a critical physical factor for cancer progression; however, its importance in tumor progression remains to be comprehensively recognized. Herein, we conducted a comprehensive pan-cancer analysis of MS subtypes across seven GI-cancer types. Using unsupervised clustering based on literature-derived MS-specific pathway signatures, the GI-tumor samples were divided into three MS subtypes, termed as the Soft, Mixed and Stiff. Then, distinct prognoses, biological features, tumor microenvironments and mutation landscapes among three MS subtypes were revealed. The Stiff tumor subtype was associated with the poorest prognosis, the most malignant biological behaviors, and the immunosuppressive tumor stromal microenvironment. Furthermore, multiple machine learning algorithms were used to develop an 11-gene MS-signature to identify the MS subtypes of GI-caner and predict chemotherapy sensitivity, which were further validated in two external GI-cancer cohorts. This novel MS-based classification on GI-cancers could enhance our understanding of the important role of MS in tumor progression, and may have implications for the optimization of individualized cancer management.

TGF-β Pathways Stratify Colorectal Cancer into Two Subtypes with Distinct Cartilage Oligomeric Matrix Protein (COMP) Expression-Related Characteristics

BACKGROUND

Colorectal cancers (CRCs) continue to be the leading cause of cancer-related deaths worldwide. The exact landscape of the molecular features of TGF-β pathway-inducing CRCs remains uncharacterized.

METHODS

Unsupervised hierarchical clustering was performed to stratify samples into two clusters based on the differences in TGF-β pathways. Weighted gene co-expression network analysis was applied to identify the key gene modules mediating the different characteristics between two subtypes. An algorithm integrating the least absolute shrinkage and selection operator (LASSO), XGBoost, and random forest regression was performed to narrow down the candidate genes. Further bioinformatic analyses were performed focusing on COMP-related immune infiltration and functions.

RESULTS

The integrated machine learning algorithm identified COMP as the hub gene, which exhibited a significant predictive value for two subtypes with an area under the curve (AUC) value equaling 0.91. Further bioinformatic analysis revealed that COMP was significantly upregulated in various cancers, especially in advanced CRCs, and regulated the immune infiltration, especially M2 macrophages and cancer-associated fibroblasts in CRCs.

CONCLUSIONS

Comprehensive immune analysis and experimental validation demonstrate that COMP is a reliable signature for subtype prediction. Our results could provide a new point for TGFβ-targeted anticancer drugs and contribute to guiding clinical decision making for CRC patients.

CpG Site-Based Signature Predicts Survival of Colorectal Cancer

BACKGROUND

A critical unmet medical need in clinical management of colorectal cancer (CRC) pivots around lack of noninvasive and or minimally invasive techniques for early diagnosis and prognostic prediction of clinical outcomes. Because DNA methylation can capture the regulatory landscape of tumors and can be measured in body fluids, it provides unparalleled opportunities for the discovery of early diagnostic and prognostics markers predictive of clinical outcomes. Here we investigated use of DNA methylation for the discovery of potential clinically actionable diagnostic and prognostic markers for predicting survival in CRC.

METHODS

We analyzed DNA methylation patterns between tumor and control samples to discover signatures of CpG sites and genes associated with CRC and predictive of survival. We conducted functional analysis to identify molecular networks and signaling pathways driving clinical outcomes.

RESULTS

We discovered a signature of aberrantly methylated genes associated with CRC and a signature of thirteen (13) CpG sites predictive of survival. We discovered molecular networks and signaling pathways enriched for CpG sites likely to drive clinical outcomes.

CONCLUSIONS

The investigation revealed that CpG sites can predict survival in CRC and that DNA methylation can capture the regulatory state of tumors through aberrantly methylated molecular networks and signaling pathways.

Construction and Evaluation of a Risk Score Model for Lymph Node Metastasis-Associated Circadian Clock Genes in Esophageal Squamous Carcinoma

Background: Studies suggested that circadian clock genes (CCGs) in human esophageal squamous carcinoma (ESCC) samples are dysregulated. However, the relevance of CCGs to lymph node metastasis (LNM) and prognosis of ESCC remains unclear. Methods: The differentially expressed genes (DEGs) between normal and ESCC samples in The Cancer Genome Atlas database (TCGA) database were intersected with the genes associated with LNM (LNMGs) in ESCC samples and 300 CCGs to obtain the differentially expressed LNM-associated CCGs (DE-LNM-CCGs). The risk model was constructed by Cox regression analysis in the TCGA-ESCC training set, and the accuracy of the risk model was verified by risk profile and overall survival profile. Furthermore, differences of 23 immune cells, 13 immune functions, and immune checkpoint molecules between the high- and low-risk groups were assessed using the single-sample gene set enrichment analysis (ssGSEA) algorithm. Gene set enrichment analysis (GSEA) was conducted to investigate the functional differences between low- and high-risk groups. Finally, we validated the mRNA expression levels of prognostic model genes by quantitative real-time polymerase chain reaction (qRT-PCR). Results: A total of six DE-LNM-CCGs were identified in TCGA-ESCC. TP53 and NAGLU were selected by Cox regression analysis to construct the risk model. Risk profile plots, overall survival plots, and validation results of the risk model in the validation set indicated that the constructed risk model was reliable. The result of ssGSEA showed that the percentages of activated B cells, activated dendritic cells, effector memory CD8 T cells, immune function in neutrophils, plasmacytoid dendritic cells, T cell co-inhibition, and Type 17 T helper cells were different between the high- and low-risk groups. In addition, the expression of CD274, PDCD1, TNFRSF18, and TNFRSF9 was dysregulated between the high- and low-risk groups. GSEA revealed that the high-risk group was associated with cell differentiation, oxidative phosphorylation, and steroid biosynthesis pathways, while the low-risk group was associated with chromosome, ECM–receptor interaction, and other pathways. Finally, qRT-PCR results showed that the mRNA expression levels of two prognostic genes were consistent with TCGA. Conclusion: In conclusion, the risk model constructed based on TP53 and NAGLU could accurately predict the prognosis.

Development of necroptosis-related gene signature to predict the prognosis of colon adenocarcinoma

Colon adenocarcinoma (COAD) is a common malignancy and has a high mortality rate. However, the current tumor node metastasis (TNM) staging system is inadequate for prognostic assessment of COAD patients. Therefore, there is an urgent need to identify reliable biomarkers for the prognosis COAD patients. The aberrant expression of necroptosis-related genes (NRGs) is reported to be associated with tumorigenesis and metastasis. In the present work, we compared the expression profiles of NRGs between COAD patients and normal individuals. Based on seven differentially expressed NRGs, a risk score was defined to predict the prognosis of COAD patients. The validation results from both training and independent external cohorts demonstrated that the risk score is able to distinguish the high and low risk COAD patients with higher accuracies, and is independent of the other clinical factors. To facilitate its clinical use, by integrating the proposed risk score, a nomogram was built to predict the risk of individual COAD patients. The C-index of the nomogram is 0.75, indicating the reliability of the nomogram in predicting survival rates. Furthermore, two candidate drugs, namely dapsone and xanthohumol, were screed out and validated by molecular docking, which hold the potential for the treatment of COAD. These results will provide novel clues for the diagnosis and treatment of COAD.

Development of a novel hypoxia-immune–related LncRNA risk signature for predicting the prognosis and immunotherapy response of colorectal cancer

Background

Colorectal cancer (CRC) is one of the most common digestive system tumors worldwide. Hypoxia and immunity are closely related in CRC; however, the role of hypoxia-immune–related lncRNAs in CRC prognosis is unknown.

Methods

Data used in the current study were sourced from the Gene Expression Omnibus and The Cancer Genome Atlas (TCGA) databases. CRC patients were divided into low- and high-hypoxia groups using the single-sample gene set enrichment analysis (ssGSEA) algorithm and into low- and high-immune groups using the Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE) algorithm. Differentially expressed lncRNAs (DElncRNAs) between low- and high-hypoxia groups, low- and high-immune groups, and tumor and control samples were identified using the limma package. Hypoxia-immune–related lncRNAs were obtained by intersecting these DElncRNAs. A hypoxia-immune–related lncRNA risk signature was developed using univariate Cox regression and least absolute shrinkage and selection operator (LASSO) analyses. The tumor microenvironments in the low- and high-risk groups were evaluated using ssGSEA, ESTIMATE, and the expression of immune checkpoints. The therapeutic response in the two groups was assessed using TIDE, IPS, and IC50. A ceRNA network based on signature lncRNAs was constructed. Finally, we used RT-qPCR to verify the expression of hypoxia-immune–related lncRNA signatures in normal and cancer tissues.

Results

Using differential expression analysis, and univariate Cox and LASSO regression analyses, ZNF667-AS1, LINC01354, LINC00996, DANCR, CECR7, and LINC01116 were selected to construct a hypoxia-immune–related lncRNA signature. The performance of the risk signature in predicting CRC prognosis was validated in internal and external datasets, as evidenced by receiver operating characteristic curves. In addition, we observed significant differences in the tumor microenvironment and immunotherapy response between low- and high-risk groups and constructed a CECR7–miRNA–mRNA regulatory network in CRC. Furthermore, RT-qPCR results confirmed that the expression patterns of the six lncRNA signatures were consistent with those in TCGA-CRC cohort.

Conclusion

Our study identified six hypoxia-immune–related lncRNAs for predicting CRC survival and sensitivity to immunotherapy. These findings may enrich our understanding of CRC and help improve CRC treatment. However, large-scale long-term follow-up studies are required for verification.

Prognostic and tumor immunity implication of inflammatory bowel disease-associated genes in colorectal cancer

BACKGROUND

Epidemiologic studies continue to emphasize that increasing patients with inflammatory bowel disease (IBD) develop to colorectal cancer (CRC). Although the function and mechanisms of IBD-associated genes (IBDGs) in CRC tumorigenesis have been extensively researched, the implications of IBDGs in the prognosis value and tumor immunity of CRC remain unclear.

RESULTS

In this study, the expression, pathological stages and prognostic value of IBDGs in CRC were systematically analyzed, and 7 prognostic genes including CDH1, CCL11, HLA-DRA, NOS2, NAT2, TIMP1 and TP53 were screened through LASSO-Cox regression analysis. Then, a prognostic signature was established based on the 7 prognostic genes, and the model exhibited a good ability in risk stratification of CRC patients. Subsequent results showed that the genetic alterations of the 7 prognostic genes exhibited more significant and extensive influence on immune cells infiltration in colon adenocarcinoma than that in rectal adenocarcinoma. Meanwhile, immune cells infiltration also showed a significant difference between low-risk group and high-risk group. What's more, 7 prognostic genes-based risk stratification was associated with microsatellite instability, and its prognostic characteristics were significantly negatively correlated with mismatch repair genes.

CONCLUSIONS

This study provided a promising insight that the 7 IBDGs could be used as valuable biomarkers for prognostic diagnosis and personalized immunotherapy of CRC patients.

ExhauFS: exhaustive search-based feature selection for classification and survival regression

Feature selection is one of the main techniques used to prevent overfitting in machine learning applications. The most straightforward approach for feature selection is an exhaustive search: one can go over all possible feature combinations and pick up the model with the highest accuracy. This method together with its optimizations were actively used in biomedical research, however, publicly available implementation is missing. We present ExhauFS-the user-friendly command-line implementation of the exhaustive search approach for classification and survival regression. Aside from tool description, we included three application examples in the manuscript to comprehensively review the implemented functionality. First, we executed ExhauFS on a toy cervical cancer dataset to illustrate basic concepts. Then, multi-cohort microarray breast cancer datasets were used to construct gene signatures for 5-year recurrence classification. The vast majority of signatures constructed by ExhauFS passed 0.65 threshold of sensitivity and specificity on all datasets, including the validation one. Moreover, a number of gene signatures demonstrated reliable performance on independent RNA-seq dataset without any coefficient re-tuning, i.e., turned out to be cross-platform. Finally, Cox survival regression models were used to fit isomiR signatures for overall survival prediction for patients with colorectal cancer. Similarly to the previous example, the major part of models passed the pre-defined concordance index threshold 0.65 on all datasets. In both real-world scenarios (breast and colorectal cancer datasets), ExhauFS was benchmarked against state-of-the-art feature selection models, including L₁-regularized sparse models. In case of breast cancer, we were unable to construct reliable cross-platform classifiers using alternative feature selection approaches. In case of colorectal cancer not a single model passed the same 0.65 threshold. Source codes and documentation of ExhauFS are available on GitHub: https://github.com/s-a-nersisyan/ExhauFS.