Identification of key genes for predicting colorectal cancer prognosis by integrated bioinformatics analysis

Identification of PANoptosis-related subtypes, construction of a prognosis signature, and tumor microenvironment landscape of hepatocellular carcinoma using bioinformatic analysis and experimental verification

Background

Hepatocellular carcinoma (HCC) is one of the most lethal malignancies worldwide. PANoptosis is a recently unveiled programmed cell death pathway, Nonetheless, the precise implications of PANoptosis within the context of HCC remain incompletely elucidated.

Methods

We conducted a comprehensive bioinformatics analysis to evaluate both the expression and mutation patterns of PANoptosis-related genes (PRGs). We categorized HCC into two clusters and identified differentially expressed PANoptosis-related genes (DEPRGs). Next, a PANoptosis risk model was constructed using LASSO and multivariate Cox regression analyses. The relationship between PRGs, risk genes, the risk model, and the immune microenvironment was studies. In addition, drug sensitivity between high- and low-risk groups was examined. The expression profiles of these four risk genes were elucidate by qRT-PCR or immunohistochemical (IHC). Furthermore, the effect of CTSC knock down on HCC cell behavior was verified using in vitro experiments.

Results

We constructed a prognostic signature of four DEPRGs (CTSC, CDCA8, G6PD, and CXCL9). Receiver operating characteristic curve analyses underscored the superior prognostic capacity of this signature in assessing the outcomes of HCC patients. Subsequently, patients were stratified based on their risk scores, which revealed that the low-risk group had better prognosis than those in the high-risk group. High-risk group displayed a lower Stromal Score, Immune Score, ESTIMATE score, and higher cancer stem cell content, tumor mutation burden (TMB) values. Furthermore, a correlation was noted between the risk model and the sensitivity to 56 chemotherapeutic agents, as well as immunotherapy efficacy, in patient with. These findings provide valuable guidance for personalized clinical treatment strategies. The qRT-PCR analysis revealed that upregulated expression of CTSC, CDCA8, and G6PD, whereas downregulated expression of CXCL9 in HCC compared with adjacent tumor tissue and normal liver cell lines. The knockdown of CTSC significantly reduced both HCC cell proliferation and migration.

Conclusion

Our study underscores the promise of PANoptosis-based molecular clustering and prognostic signatures in predicting patient survival and discerning the intricacies of the tumor microenvironment within the context of HCC. These insights hold the potential to advance our comprehension of the therapeutic contribution of PANoptosis plays in HCC and pave the way for generating more efficacious treatment strategies.

Glucose metabolism-based signature predicts prognosis and immunotherapy strategies for colon adenocarcinoma

BACKGROUND

The global prevalence and metastasis rates of colon adenocarcinoma (COAD) are high, and therapeutic success is limited. Although previous research has primarily explored changes in gene phenotypes, the incidence rate of COAD remains unchanged. Metabolic reprogramming is a crucial aspect of cancer research and therapy. The present study aims to develop cluster and polygenic risk prediction models for COAD based on glucose metabolism pathways to assess the survival status of patients and potentially identify novel immunotherapy strategies and related therapeutic targets.

METHODS

COAD-specific data (including clinicopathological information and gene expression profiles) were sourced from The Cancer Genome Atlas (TCGA) and two Gene Expression Omnibus (GEO) datasets (GSE33113 and GSE39582). Gene sets related to glucose metabolism were obtained from the MSigDB database. The Gene Set Variation Analysis (GSVA) method was utilized to calculate pathway scores for glucose metabolism. The hclust function in R, part of the Pheatmap package, was used to establish a clustering system. The mutation characteristics of identified clusters were assessed via MOVICS software, and differentially expressed genes (DEGs) were filtered using limma software. Signature analysis was performed using the least absolute shrinkage and selection operator (LASSO) method. Survival curves, survival receiver operating characteristic (ROC) curves and multivariate Cox regression were analyzed to assess the efficacy and accuracy of the signature for prognostic prediction. The pRRophetic program was employed to predict drug sensitivity, with data sourced from the Genomics of Drug Sensitivity in Cancer (GDSC) database.

RESULTS

Four COAD subgroups (i.e., C1, C2, C3 and C4) were identified based on glucose metabolism, with the C4 group having higher survival rates. These four clusters were bifurcated into a new Clust2 system (C1 + C2 + C3 and C4). In total, 2175 DEGs were obtained (C1 + C2 + C3 vs. C4), from which 139 prognosis-related genes were identified. ROC curves predicting 1-, 3- and 5-year survival based on a signature containing nine genes showed an area under the curve greater than 0.7. Meanwhile, the study also found this feature to be an important predictor of prognosis in COAD and accordingly assessed the risk score, with higher risk scores being associated with a worse prognosis. The high-risk and low-risk groups responded differently to immunotherapy and chemotherapeutic agents, and there were differences in functional enrichment pathways.

CONCLUSIONS

This unique signature based on glucose metabolism may potentially provide a basis for predicting patient prognosis, biological characteristics and more effective immunotherapy strategies for COAD.

Identification of novel T cell proliferation patterns, potential biomarkers and therapeutic drugs in colorectal cancer

Background: T cells are crucial components of antitumor immunity. A list of genes associated with T cell proliferation was recently identified; however, the impact of T cell proliferation-related genes (TRGs) on the prognosis and therapeutic responses of patients with colorectal cancer (CRC) remains unclear. Methods: 33 TRG expression information and clinical information of patients with CRC gathered from multiple datasets were subjected to bioinformatic analysis. Consensus clustering was used to determine the molecular subtypes associated with T cell proliferation. Utilizing the Lasso-Cox regression, a predictive signature was created and verified in external cohorts. A tumor immune environment analysis was conducted, and potential biomarkers and therapeutic drugs were identified and confirmed via in vitro and in vivo studies. Results: CRC patients were separated into two TRG clusters, and differentially expressed genes (DEGs) were identified. Patient information was divided into three different gene clusters, and the determined molecular subtypes were linked to patient survival, immune cells, and immune functions. Prognosis-associated DEGs in the three gene clusters were used to evaluate the risk score, and a predictive signature was developed. The ability of the risk score to predict patient survival and treatment response has been successfully validated using multiple datasets. To discover more possible biomarkers for CRC, the weighted gene co-expression network analysis algorithm was utilized to screen key TRG variations between groups with high- and low-risk. CDK1, BATF, IL1RN, and ITM2A were screened out as key TRGs, and the expression of key TRGs was confirmed using real-time reverse transcription polymerase chain reaction. According to the key TRGs, 7,8-benzoflavone was identified as the most significant drug molecule, and MTT, colony formation, wound healing, transwell assays, and in vivo experiments indicated that 7,8-benzoflavone significantly suppressed the proliferation and migration of CRC cells. Conclusion: T cell proliferation-based molecular subtypes and predictive signatures can be utilized to anticipate patient results, immunological landscape, and treatment response in CRC. Novel biomarker candidates and potential therapeutic drugs for CRC were identified and verified using in vitro and in vivo tests.

Bioinformatics screening of colorectal-cancer causing molecular signatures through gene expression profiles to discover therapeutic targets and candidate agents

BACKGROUND

Detection of appropriate receptor proteins and drug agents are equally important in the case of drug discovery and development for any disease. In this study, an attempt was made to explore colorectal cancer (CRC) causing molecular signatures as receptors and drug agents as inhibitors by using integrated statistics and bioinformatics approaches.

METHODS

To identify the important genes that are involved in the initiation and progression of CRC, four microarray datasets (GSE9348, GSE110224, GSE23878, and GSE35279) and an RNA_Seq profiles (GSE50760) were downloaded from the Gene Expression Omnibus database. The datasets were analyzed by a statistical r-package of LIMMA to identify common differentially expressed genes (cDEGs). The key genes (KGs) of cDEGs were detected by using the five topological measures in the protein-protein interaction network analysis. Then we performed in-silico validation for CRC-causing KGs by using different web-tools and independent databases. We also disclosed the transcriptional and post-transcriptional regulatory factors of KGs by interaction network analysis of KGs with transcription factors (TFs) and micro-RNAs. Finally, we suggested our proposed KGs-guided computationally more effective candidate drug molecules compared to other published drugs by cross-validation with the state-of-the-art alternatives of top-ranked independent receptor proteins.

RESULTS

We identified 50 common differentially expressed genes (cDEGs) from five gene expression profile datasets, where 31 cDEGs were downregulated, and the rest 19 were up-regulated. Then we identified 11 cDEGs (CXCL8, CEMIP, MMP7, CA4, ADH1C, GUCA2A, GUCA2B, ZG16, CLCA4, MS4A12 and CLDN1) as the KGs. Different pertinent bioinformatic analyses (box plot, survival probability curves, DNA methylation, correlation with immune infiltration levels, diseases-KGs interaction, GO and KEGG pathways) based on independent databases directly or indirectly showed that these KGs are significantly associated with CRC progression. We also detected four TFs proteins (FOXC1, YY1, GATA2 and NFKB) and eight microRNAs (hsa-mir-16-5p, hsa-mir-195-5p, hsa-mir-203a-3p, hsa-mir-34a-5p, hsa-mir-107, hsa-mir-27a-3p, hsa-mir-429, and hsa-mir-335-5p) as the key transcriptional and post-transcriptional regulators of KGs. Finally, our proposed 15 molecular signatures including 11 KGs and 4 key TFs-proteins guided 9 small molecules (Cyclosporin A, Manzamine A, Cardidigin, Staurosporine, Benzo[A]Pyrene, Sitosterol, Nocardiopsis Sp, Troglitazone, and Riccardin D) were recommended as the top-ranked candidate therapeutic agents for the treatment against CRC.

CONCLUSION

The findings of this study recommended that our proposed target proteins and agents might be considered as the potential diagnostic, prognostic and therapeutic signatures for CRC.

Identification of cuproptosis-based molecular subtypes, construction of prognostic signature and characterization of immune landscape in colon cancer

Background

Cuproptosis is a newly discovered form of cell death induced by targeting lipoacylated proteins involved in the tricarboxylic acid cycle. However, the roles of cuproptosis-related genes (CRGs) in the clinical outcomes and immune landscape of colon cancer remain unknown.

Methods

We performed bioinformatics analysis of the expression data of 13 CRGs identified from a previous study and clinical information of patients with colon cancer obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases. Colon cancer cases were divided into two CRG clusters and prognosis-related differentially expressed genes. Patient data were separated into three corresponding distinct gene clusters, and the relationships between the risk score, patient prognosis, and immune landscape were analyzed. The identified molecular subtypes correlated with patient survival, immune cells, and immune functions. A prognostic signature based on five genes was identified, and the patients were divided into high- and low-risk groups based on the calculated risk score. A nomogram model for predicting patient survival was developed based on the risk score and other clinical features.

Results

The high-risk group showed a worse prognosis, and the risk score was related to immune cell abundance, microsatellite instability, cancer stem cell index, checkpoint expression, immune escape, and response to chemotherapeutic drugs and immunotherapy. Findings related to the risk score were validated in the imvigor210 cohort of patients with metastatic urothelial cancer treated with anti-programmed cell death ligand 1.

Conclusion

We demonstrated the potential of cuproptosis-based molecular subtypes and prognostic signatures for predicting patient survival and the tumor microenvironment in colon cancer. Our findings may improve the understanding of the role of cuproptosis in colon cancer and lead to the development of more effective treatment strategies.

A novel oxidative stress- and ferroptosis-related gene prognostic signature for distinguishing cold and hot tumors in colorectal cancer

Oxidative stress and ferroptosis exhibit crosstalk in many types of human diseases, including malignant tumors. We aimed to develop an oxidative stress- and ferroptosis-related gene (OFRG) prognostic signature to predict the prognosis and therapeutic response in patients with colorectal cancer (CRC). Thirty-four insertion genes between oxidative stress-related genes and ferroptosis-related genes were identified as OFRGs. We then performed bioinformatics analysis of the expression profiles of 34 OFRGs and clinical information of patients obtained from multiple datasets. Patients with CRC were divided into three OFRG clusters, and differentially expressed genes (DEGs) between clusters were identified. OFRG clusters correlated with patient survival and immune cell infiltration. Prognosis-related DEGs in three clusters were used to calculate the risk score, and a prognostic signature was constructed according to the risk score. In this study, patients in the low-risk group had better prognosis, higher immune cell infiltration levels, and better responses to fluorouracil-based chemotherapy and immune checkpoint blockade therapy than high-risk patients; these results were successfully validated with multiple independent datasets. Thus, low-risk CRC could be defined as hot tumors and high-risk CRC could be defined as cold tumors. To further identify potential biomarkers for CRC, the expression levels of five signature genes in CRC and adjacent normal tissues were further verified via an in vitro experiment. In conclusion, we identified 34 OFRGs and constructed an OFRG-related prognostic signature, which showed excellent performance in predicting survival and therapeutic responses for patients with CRC. This could help to distinguish cold and hot tumors in CRC, and the results might be helpful for precise treatment protocols in clinical practice.

PANoptosis-based molecular clustering and prognostic signature predicts patient survival and immune landscape in colon cancer

PANoptosis is a newly-discovered cell death pathway that involves crosstalk and co-ordination between pyroptosis, apoptosis, and necroptosis processes. However, the roles of PANoptosis-related genes (PRGs) in prognosis and immune landscape of colon cancer remain widely unknown. Here, we performed a bioinformatics analysis of expression data of nineteen PRGs identified from previous studies and clinical data of colon cancer patients obtained from TCGA and GEO databases. Colon cancer cases were divided into two PRG clusters, and prognosis-related differentially expressed genes (PRDEGs) were identified. The patient data were then separated into two corresponding distinct gene clusters, and the relationship between the risk score, patient prognosis, and immune landscape was analyzed. The identified PRGs and gene clusters correlated with patient survival and immune system and cancer-related biological processes and pathways. A prognosis signature based on seven genes was identified, and patients were divided into high-risk and low-risk groups based on the calculated risk score. A nomogram model for prediction of patient survival was also developed based on the risk score and other clinical features. Accordingly, the high-risk group showed worse prognosis, and the risk score was related to immune cell abundance, cancer stem cell (CSC) index, checkpoint expression, and response to immunotherapy and chemotherapeutic drugs. Results of quantitative real-time polymerase chain reaction (qRT-PCR) showed that LGR5 and VSIG4 were differentially expressed between normal and colon cancer samples. In conclusion, we demonstrated the potential of PANoptosis-based molecular clustering and prognostic signatures for prediction of patient survival and tumor microenvironment (TME) in colon cancer. Our findings may improve our understanding of the role of PANoptosis in colon cancer, and enable the development of more effective treatment strategies.

A network-based pharmacological investigation to identify the mechanistic regulatory pathway of andrographolide against colorectal cancer

Traditional cancer treatments have posed numerous obstacles, including toxicity, multiple drug resistance, and financial cost. On the contrary, bioactive phytochemicals used in complementary alternative medicine have recently increased attention due to their potential to modulate a wide range of molecular mechanisms with a less toxic effect. Therefore, we investigated the potential regulatory mechanisms of andrographolide to treat colorectal cancer (CRC) using a network pharmacology approach. Target genes of andrographolide were retrieved from public databases (PharmMapper, Swiss target prediction, Targetnet, STITCH, and SuperPred), while targets related to CRC were retrieved from disease databases (Genecards and DisGeNet) and expression datasets (GSE32323 and GSE8671) were retrieved from gene expression omnibus (GEO). Protein-protein interaction networks (PPI) were generated using STRING and Cytoscape, and hub genes were identified by topology analysis and MCODE. Annotation of target proteins was performed using Gene Ontology (GO) database DAVID and signaling pathway enrichment analysis using the Kyoto Encyclopedia and Genome Database (KEGG). Survival and molecular docking analysis for the hub genes revealed three genes (PDGFRA, PTGS2, and MMP9) were involved in the overall survival of CRC patients, and the top three genes with the lowest binding energy include PDGFRA, MET, and MAPK1. MET gene upregulation and PDGFRA and PTGS2 gene downregulation are associated with the survival of CRC patients, as revealed by box plots and correlation analysis. In conclusion, this study has provided the first scientific evidence to support the use of andrographolide to inhibit cellular proliferation, migration, and growth, and induce apoptosis by targeting the hub genes (PDGFRA, PTGS2, MMP9, MAPK1, and MET) involved in CRC migration and invasion.

Innovative signature establishment using lymphangiogenesis-related lncRNA pairs to predict prognosis of hepatocellular carcinoma

Aims

Hepatocellular carcinoma (HCC) remains a major tumoral burden globally, and its heterogeneity encumbers prognostic prediction. The lymphangiogenesis-related long non-coding RNAs (lrlncRNAs) reported to be implicated in immune response regulation show potential importance in predicting the prognostic and therapeutic outcome. Hence, this study aims to establish a lrlncRNA pairs-based signature not requiring specific expression levels of transcripts, which displays promising clinical practicality and satisfactory predictive capability.

Main methods

Transcriptomic and clinical information of the Liver Hepatocellular Carcinoma (LIHC) project retrieved from the TCGA portal were used to find differently expressed lrlncRNA (DElrlncRNA) via analysis performed between lymphangiogenesis-related genes (lr-genes) and lncRNAs(lrlncRNA), and to ultimately construct the signature based on lrlncRNA pairs screened out via Lasso and Cox regression analyses. Akaike information criterion (AIC) values were computed to find the cut-off point optimum for high-risk and low-risk group allocation. The signature then underwent trials in terms of its predictive value for survival, clinicopathological features, immune cells infiltration in tumoral microenvironment, selected checkpoint biomarkers and chemosensitivity.

Key findings

A novel lymphangiogenesis-related lncRNA pair signature was established using nine lrlncRNA pairs identified and significantly related to overall survival, clinicopathological features, immune cells infiltration and susceptibility to chemotherapy. Moreover, the signature efficacy was verified in acknowledged clinicopathological subgroups and partially validated by qRT-PCR assay in various human HCC cell lines.

Significance

The novel lrlncRNA-pairs based signature was shown to effectively and independently estimate HCC prognosis and help screen patients suitable for anti-tumor immunotherapy and chemotherapy.

Identification of Six Genes as Diagnostic Markers for Colorectal Cancer Detection by Integrating Multiple Expression Profiles

Background

Many studies have demonstrated the promising utility of DNA methylation and miRNA as biomarkers for colorectal cancer (CRC) early detection. However, mRNA is rarely reported. This study aimed to identify novel fecal-based mRNA signatures.

Methods

The differentially expressed genes (DEGs) were first determined between CRCs and matched normal samples by integrating multiple datasets. Then, Least Absolute Shrinkage and Selection Operator (LASSO) regression was used to reduce the number of candidates of aberrantly expressed genes. Next, the potential functions were investigated for the candidate signatures and their ability to detect CRC and pan-cancers was comprehensively evaluated.

Results

We identified 1841 common DEGs in two independent datasets. Functional enrichment analysis revealed they were mainly related to extracellular structure, biosynthesis, and cell adhesion. The CRC classifier was established based on six genes screened by LASSO regression. Sensitivity, specificity, and area under the ROC curve (AUC) for CRC detection were 79.30%, 80.40%, and 0.85 (0.76–0.92) in the training set, and these indexes achieved 93.20%, 41.80%, and 0.73 (0.65–0.83) in the testing set. For validation set, the sensitivity, specificity, and AUC were 98.90%, 98.00%, and 0.97 (0.94–0.99). The average sensitivities exceeded 90.00% for CRCs with different clinical features. For adenomas detection, the sensitivity and specificity were 74.50% and 64.00%. Besides, the six genes obtained an average AUC of 0.855 for pan-cancer detection.

Conclusion

The six-gene signatures showed ability to detect CRC and pan-cancer samples, which could be served as potential diagnostic markers.

A Diagnostic Model Using Exosomal Genes for Colorectal Cancer

Colorectal cancer (CRC) is a leading cause of cancer-related deaths worldwide. Exosomes have great potential as liquid biopsy specimens due to their presence and stability in body fluids. However, the function and diagnostic values of exosomal genes in CRC are poorly understood. In the present study, exosomal data of CRC and healthy samples from the exoRBase 2.0 and Gene Expression Omnibus (GEO) databases were used, and 38 common exosomal genes were identified. Through the least absolute shrinkage and selection operator (Lasso) analysis, support vector machine recursive feature elimination (SVM-RFE) analysis, and logistic regression analysis, a diagnostic model of the training set was constructed based on 6 exosomal genes. The diagnostic model was internally validated in the test and exoRBase 2.0 database and externally validated in the GEO database. In addition, the co-expression analysis was used to cluster co-expression modules, and the enrichment analysis was performed on module genes. Then a protein–protein interaction and competing endogenous RNA network were constructed and 10 hub genes were identified using module genes. In conclusion, the results provided a comprehensive understanding of the functions of exosomal genes in CRC as well as a diagnostic model related to exosomal genes.

Machine Learning-Based Identification of Colon Cancer Candidate Diagnostics Genes

Simple Summary

We developed a predictive approach using different machine learning methods to identify a number of genes that can potentially serve as novel diagnostic colon cancer biomarkers.

Abstract

Background: Colorectal cancer (CRC) is the third leading cause of cancer-related death and the fourth most commonly diagnosed cancer worldwide. Due to a lack of diagnostic biomarkers and understanding of the underlying molecular mechanisms, CRC’s mortality rate continues to grow. CRC occurrence and progression are dynamic processes. The expression levels of specific molecules vary at various stages of CRC, rendering its early detection and diagnosis challenging and the need for identifying accurate and meaningful CRC biomarkers more pressing. The advances in high-throughput sequencing technologies have been used to explore novel gene expression, targeted treatments, and colon cancer pathogenesis. Such approaches are routinely being applied and result in large datasets whose analysis is increasingly becoming dependent on machine learning (ML) algorithms that have been demonstrated to be computationally efficient platforms for the identification of variables across such high-dimensional datasets. Methods: We developed a novel ML-based experimental design to study CRC gene associations. Six different machine learning methods were employed as classifiers to identify genes that can be used as diagnostics for CRC using gene expression and clinical datasets. The accuracy, sensitivity, specificity, F1 score, and area under receiver operating characteristic (AUROC) curve were derived to explore the differentially expressed genes (DEGs) for CRC diagnosis. Gene ontology enrichment analyses of these DEGs were performed and predicted gene signatures were linked with miRNAs. Results: We evaluated six machine learning classification methods (Adaboost, ExtraTrees, logistic regression, naïve Bayes classifier, random forest, and XGBoost) across different combinations of training and test datasets over GEO datasets. The accuracy and the AUROC of each combination of training and test data with different algorithms were used as comparison metrics. Random forest (RF) models consistently performed better than other models. In total, 34 genes were identified and used for pathway and gene set enrichment analysis. Further mapping of the 34 genes with miRNA identified interesting miRNA hubs genes. Conclusions: We identified 34 genes with high accuracy that can be used as a diagnostics panel for CRC.

Predictors of colorectal cancer survival using cox regression and random survival forests models based on gene expression data

Understanding and identifying the markers and clinical information that are associated with colorectal cancer (CRC) patient survival is needed for early detection and diagnosis. In this work, we aimed to build a simple model using Cox proportional hazards (PH) and random survival forest (RSF) and find a robust signature for predicting CRC overall survival. We used stepwise regression to develop Cox PH model to analyse 54 common differentially expressed genes from three mutations. RSF is applied using log-rank and log-rank-score based on 5000 survival trees, and therefore, variables important obtained to find the genes that are most influential for CRC survival. We compared the predictive performance of the Cox PH model and RSF for early CRC detection and diagnosis. The results indicate that SLC9A8, IER5, ARSJ, ANKRD27, and PIPOX genes were significantly associated with the CRC overall survival. In addition, age, sex, and stages are also affecting the CRC overall survival. The RSF model using log-rank is better than log-rank-score, while log-rank-score needed more trees to stabilize. Overall, the imputation of missing values enhanced the model’s predictive performance. In addition, Cox PH predictive performance was better than RSF.

Construction of a novel immune-related lncRNA signature and its potential to predict the immune status of patients with hepatocellular carcinoma

Background

The accuracy of existing biomarkers for predicting the prognosis of hepatocellular carcinoma (HCC) is not satisfactory. It is necessary to explore biomarkers that can accurately predict the prognosis of HCC.

Methods

In this study, original transcriptome data were downloaded from The Cancer Genome Atlas (TCGA) database. Immune-related long noncoding ribonucleic acids (irlncRNAs) were identified by coexpression analysis, and differentially expressed irlncRNA (DEirlncRNA) pairs were distinguished by univariate analysis. In addition, the least absolute shrinkage and selection operator (LASSO) penalized regression was modified. Next, the cutoff point was determined based on the area under the curve (AUC) and Akaike information criterion (AIC) values of the 5-year receiver operating characteristic (ROC) curve to establish an optimal model for identifying high-risk and low-risk groups of HCC patients. The model was then reassessed in terms of clinicopathological features, survival rate, tumor-infiltrating immune cells, immunosuppressive markers, and chemotherapy efficacy.

Results

A total of 1009 pairs of DEirlncRNAs were recognized in this study, 30 of these pairs were included in the Cox regression model for subsequent analysis. After regrouping according to the cutoff point, we could more effectively identify factors such as aggressive clinicopathological features, poor survival outcomes, specific immune cell infiltration status of tumors, high expression level of immunosuppressive biomarkers, and low sensitivity to chemotherapy drugs in HCC patients.

Conclusions

The nonspecific expression level signature involved with irlncRNAs shows promising clinical value in predicting the prognosis of HCC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-09059-x.

A Network-Based Methodology to Identify Subnetwork Markers for Diagnosis and Prognosis of Colorectal Cancer

The development of reliable methods for identification of robust biomarkers for complex diseases is critical for disease diagnosis and prognosis efforts. Integrating multi-omics data with protein-protein interaction (PPI) networks to investigate diseases may help better understand disease characteristics at the molecular level. In this study, we developed and tested a novel network-based method to detect subnetwork markers for patients with colorectal cancer (CRC). We performed an integrated omics analysis using whole-genome gene expression profiling and copy number alterations (CNAs) datasets followed by building a gene interaction network for the significantly altered genes. We then clustered the constructed gene network into subnetworks and assigned a score for each significant subnetwork. We developed a support vector machine (SVM) classifier using these scores as feature values and tested the methodology in independent CRC transcriptomic datasets. The network analysis resulted in 15 subnetwork markers that revealed several hub genes that may play a significant role in colorectal cancer, including PTP4A3, FGFR2, PTX3, AURKA, FEN1, INHBA, and YES1. The 15-subnetwork classifier displayed over 98 percent accuracy in detecting patients with CRC. In comparison to individual gene biomarkers, subnetwork markers based on integrated multi-omics and network analyses may lead to better disease classification, diagnosis, and prognosis.

Using Immune-Related Long Non-coding Ribonucleic Acids to Develop a Novel Prognosis Signature and Predict the Immune Landscape of Colon Cancer

Purpose: This study aimed to construct a novel signature to predict the survival of patients with colon cancer and the associated immune landscape, based on immune-related long noncoding ribonucleic acids (irlncRNAs).

Methods: Expression profiles of irlncRNAs in 457 patients with colon cancer were retrieved from the TCGA database (https://portal.gdc.cancer.gov). Differentially expressed (DE) irlncRNAs were identified and irlncRNA pairs were recognized using Lasso regression and Cox regression analyses. Akaike information criterion (AIC) values of receiver operating characteristic (ROC) curve were calculated to identify the ideal cut-off point for dividing patients into two groups and constructing the prognosis signature. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed to validate the expression of LINC02195 and SCARNA9 in colon cancer.

Results: We identified 22 irlncRNA pairs and patients were divided into high-risk and low-risk groups based on the calculated risk score using these 22 irlncRNA pairs. The irlncRNA pairs were significantly related to patient survival. Low-risk patients had a significantly longer survival time than high-risk patients (p < 0.001). The area under the curve of the signature to predict 5-year survival was 0.951. The risk score correlated with tumor stage, infiltration depth, lymph node metastasis, and distant metastasis. The risk score remained significant after univariate and multivariate Cox regression analyses. A nomogram model to predict patient survival was developed based on the results of Cox regression analysis. Immune cell infiltration status, expression of some immune checkpoint genes, and sensitivity to chemotherapeutics were also related to the risk score. The results of qRT-PCR revealed that LINC02195 and SCARNA9 were significantly upregulated in colon cancer tissues.

Conclusion: The constructed prognosis signature showed remarkable efficiency in predicting patient survival, immune cell infiltration status, expression of immune checkpoint genes, and sensitivity to chemotherapeutics.

Novel diagnostic and prognostic biomarkers of colorectal cancer: Capable to overcome the heterogeneity-specific barrier and valid for global applications

INTRODUCTION

The heterogeneity-specific nature of the available colorectal cancer (CRC) biomarkers is significantly contributing to the cancer-associated high mortality rate worldwide. Hence, this study was initiated to investigate a system of novel CRC biomarkers that could commonly be employed to the CRC patients and helpful to overcome the heterogenetic-specific barrier.

METHODS

Initially, CRC-related hub genes were extracted through PubMed based literature mining. A protein-protein interaction (PPI) network of the extracted hub genes was constructed and analyzed to identify few more closely CRC-related hub genes (real hub genes). Later, a comprehensive bioinformatics approach was applied to uncover the diagnostic and prognostic role of the identified real hub genes in CRC patients of various clinicopathological features.

RESULTS

Out of 210 collected hub genes, in total 6 genes (CXCL12, CXCL8, AGT, GNB1, GNG4, and CXCL1) were identified as the real hub genes. We further revealed that all the six real hub genes were significantly dysregulated in colon adenocarcinoma (COAD) patients of various clinicopathological features including different races, cancer stages, genders, age groups, and body weights. Additionally, the dysregulation of real hub genes has shown different abnormal correlations with many other parameters including promoter methylation, overall survival (OS), genetic alterations and copy number variations (CNVs), and CD8+T immune cells level. Finally, we identified a potential miRNA and various chemotherapeutic drugs via miRNA, and real hub genes drug interaction network that could be used in the treatment of CRC by regulating the expression of real hub genes.

CONCLUSION

In conclusion, we have identified six real hub genes as potential biomarkers of CRC patients that could help to overcome the heterogenetic-specific barrier across different clinicopathological features.

Identification and Validation of a Novel Immune-Related lncRNA Signature for Bladder Cancer

Purpose

We aimed to construct an immune-related long noncoding ribonucleic acids (irlncRNA) signature to evaluate the prognosis of patients without specific expression level of these irlncRNA.

Methods

The raw transcriptome data were downloaded from The Cancer Genome Atlas (TCGA), irlncRNAs were filtered out using an online immune related gene database and coexpression analysis, differently expressed irlncRNA (DEirlncRNA) pairs were identified by univariate analysis. The areas under curve (AUC) were compared and the Akaike information criterion (AIC) values of receiver operating curve (ROC) was counted, the most optimal model was constructed to divide bladder cancer patients into high- and low-risk groups usingõ the cut-off point of ROC. Then, we evaluated them from multiple perspectives, such as survival time, clinic-pathological characteristics, immune-related cells infiltrating, chemotherapeutics efficacy and immune checkpoint inhibitors.

Results

14 DEirlncRNA pairs were included in this signature. Patients in high-risk groups demonstrated apparent shorter survival time, more aggressive clinic-pathological characteristics, different immune-related cells infiltrating status, lower chemotherapeutics efficacy.

Conclusion

The irlncRNA signature demonstrated a promising prediction value for bladder cancer patients and was important in guiding clinical treatment.

Identification of the miRNA signature and key genes in colorectal cancer lymph node metastasis

BACKGROUND

Because its metastasis to the lymph nodes are closely related to poor prognosis, miRNAs and mRNAs can serve as biomarkers for the diagnosis, prognosis, and therapy of colorectal cancer (CRC). This study aimed to identify novel gene signatures in the lymph node metastasis of CRC.

METHODS

GSE56350, GSE70574, and GSE95109 datasets were downloaded from the Gene Expression Omnibus (GEO) database, while data from 569 colorectal cancer cases were also downloaded from The Cancer Genome Atlas (TCGA) database. Differentially expressed miRNAs (DE-miRNAs) were calculated using R programming language (Version 3.6.3), while gene ontology and enrichment analysis of target mRNAs were performed using FunRich ( http://www.funrich.org ). Furthermore, the mRNA-miRNA network was constructed using Cytoscape software (Version 3.8.0). Gene expression levels were verified using the GEO datasets. Similarly, quantitative real-time PCR (qPCR) was used to examine expression profiles from 20 paired non-metastatic and metastatic lymph node tissue samples obtained from patients with CRC.

RESULTS

In total, five DE-miRNAs were selected, and 34 mRNAs were identified after filtering the results. Moreover, two key miRNAs (hsa-miR-99a, hsa-miR-100) and one gene (heparan sulfate-glucosamine 3-sulfotransferase 2 [HS3ST2]) were identified. The GEO datasets analysis and qPCR results showed that the expression of key miRNA and genes were consistent with that obtained from the bioinformatic analysis. A novel miRNA-mRNA network capable of predicting the prognosis and confirmed experimentally, hsa-miR-99a-HS3ST2-hsa-miR-100, was found after expression analysis in metastasized lymph node tissue from CRC samples.

CONCLUSION

In summary, miRNAs and genes with potential as biomarkers were found and a novel miRNA-mRNA network was established for CRC lymph node metastasis by systematic bioinformatic analysis and experimental validation. This network may be used as a potential biomarker in the development of lymph node metastatic CRC.

Identification of Prognostic miRNAs Associated With Immune Cell Tumor Infiltration Predictive of Clinical Outcomes in Patients With Non-Small Cell Lung Cancer

Background

A detailed means of prognostic stratification in patients with non-small cell lung cancer (NSCLC) is urgently needed to support individualized treatment plans. Recently, microRNAs (miRNAs) have been used as biomarkers due to their previously reported prognostic roles in cancer. This study aimed to construct an immune-related miRNA signature that effectively predicts NSCLC patient prognosis.

Methods

The miRNAs and mRNA expression and mutation data of NSCLC was obtained from The Cancer Genome Atlas (TCGA). Immune-associated miRNAs were identified using immune scores calculated by the ESTIMATE algorithm. LASSO-penalized multivariate survival models were using for development of a tumor immune-related miRNA signature (TIM-Sig), which was evaluated in several public cohorts from the Gene Expression Omnibus (GEO) and the CellMiner database. The miRTarBase was used for constructing the miRNA-target interactions.

Results

The TIM-Sig, including 10 immune-related miRNAs, was constructed and successfully predicted overall survival (OS) in the validation cohorts. TIM-Sig score negatively correlated with CD8+ T cell infiltration, IFN-γ expression, CYT activity, and tumor mutation burden. The correlation between TIM-Sig score and genomic mutation and cancer chemotherapeutics was also evaluated. A miRNA-target network of 10 miRNAs in TIM-Sig was constructed. Further analysis revealed that these target genes showed prognostic value in both lung squamous cell carcinoma and adenocarcinoma.

Conclusions

We concluded that the immune-related miRNAs demonstrated a potential value in clinical prognosis.

Identification of Hub Genes in Different Stages of Colorectal Cancer through an Integrated Bioinformatics Approach

Colorectal cancer (CRC) is the third most common cancer that contributes to cancer-related morbidity. However, the differential expression of genes in different phases of CRC is largely unknown. Moreover, very little is known about the role of stress-survival pathways in CRC. We sought to discover the hub genes and identify their roles in several key pathways, including oxidative stress and apoptosis in the different stages of CRC. To identify the hub genes that may be involved in the different stages of CRC, gene expression datasets were obtained from the gene expression omnibus (GEO) database. The differentially expressed genes (DEGs) common among the different datasets for each group were obtained using the robust rank aggregation method. Then, gene enrichment analysis was carried out with Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases. Finally, the protein-protein interaction networks were constructed using the Cytoscape software. We identified 40 hub genes and performed enrichment analysis for each group. We also used the Oncomine database to identify the DEGs related to stress-survival and apoptosis pathways involved in different stages of CRC. In conclusion, the hub genes were found to be enriched in several key pathways, including the cell cycle and p53 signaling pathway. Some of the hub genes were also reported in the stress-survival and apoptosis pathways. The hub DEGs revealed from our study may be used as biomarkers and may explain CRC development and progression mechanisms.

The Role of AKR1B10 in Physiology and Pathophysiology

AKR1B10 is a human nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reductase belonging to the aldo-keto reductase (AKR) 1B subfamily. It catalyzes the reduction of aldehydes, some ketones and quinones, and interacts with acetyl-CoA carboxylase and heat shock protein 90α. The enzyme is highly expressed in epithelial cells of the stomach and intestine, but down-regulated in gastrointestinal cancers and inflammatory bowel diseases. In contrast, AKR1B10 expression is low in other tissues, where the enzyme is upregulated in cancers, as well as in non-alcoholic fatty liver disease and several skin diseases. In addition, the enzyme's expression is elevated in cancer cells resistant to clinical anti-cancer drugs. Thus, growing evidence supports AKR1B10 as a potential target for diagnosing and treating these diseases. Herein, we reviewed the literature on the roles of AKR1B10 in a healthy gastrointestinal tract, the development and progression of cancers and acquired chemoresistance, in addition to its gene regulation, functions, and inhibitors.

Immune-Related lncRNA to Construct Novel Signature and Predict the Immune Landscape of Human Hepatocellular Carcinoma

The signature composed of immune-related long noncoding ribonucleic acids (irlncRNAs) with no requirement of specific expression level seems to be valuable in predicting the survival of patients with hepatocellular carcinoma (HCC). Here, we retrieved raw transcriptome data from The Cancer Genome Atlas (TCGA), identified irlncRNAs by co-expression analysis, and recognized differently expressed irlncRNA (DEirlncRNA) pairs using univariate analysis. In addition, we modified Lasso penalized regression. Then, we compared the areas under curve, counted the Akaike information criterion (AIC) values of 5-year receiver operating characteristic curve, and identified the cut-off point to set up an optimal model for distinguishing the high- or low-disease-risk groups among patients with HCC. We then reevaluated them from the viewpoints of survival, clinic-pathological characteristics, tumor-infiltrating immune cells, chemotherapeutics efficacy, and immunosuppressed biomarkers. 36 DEirlncRNA pairs were identified, 12 of which were included in a Cox regression model. After regrouping the patients by the cut-off point, we could more effectively differentiate between them based on unfavorable survival outcome, aggressive clinic-pathological characteristics, specific tumor immune infiltration status, low chemotherapeutics sensitivity, and highly expressed immunosuppressed biomarkers. The signature established by paring irlncRNA regardless of expression levels showed a promising clinical prediction value.

Prognostic Value of SLC4A4 and its Correlation with Immune Infiltration in Colon Adenocarcinoma

BACKGROUND SLC4A4 is differentially expressed in a variety of tumors, but its significance in colon adenocarcinoma has not been determined. MATERIAL AND METHODS Transcriptomes of two cohorts, GSE41258 and GSE32323, contained in The Cancer Genome Atlas (TCGA) were analysed to determine differences in SLC4A4 expression between tumor and normal tissue and their correlations with overall survival. The relationships between SLC4A4 expression and clinical characteristics were determined by COX regression analysis and logistic regression analysis, and correlations of SLC4A4 levels with tumor infiltrating immune cells (TIICs) and genes with high mutation frequency were evaluated by Pearson correlation analysis. Molecular functions and signaling pathways that might be affected by changes in SLC4A4 expression were determined by gene set enrichment analysis (GSEA). The overall distribution of TIICs was determined by two web servers: tumor immune estimation resource (TIMER) and CIBERSORT. RESULTS SLC4A4 expression was lower in colon adenocarcinoma than in normal colon tissue, suggesting that SLC4A4 was associated with poor prognosis. Reduced SLC4A4 expression was also associated with lymph node invasion and distant metastasis and was moderately correlated with increased expression of MUC4 and SMAD4, two genes with high mutation frequency in colon adenocarcinoma. GSEA indicated that changes in SLC4A4 expression affects several biological processes, including mismatch repair, base excision repair, and DNA replication. Eight TIICs in the tumor microenvironment differed significantly in groups with low and high expression of SLC4A4. CONCLUSIONS SLC4A4 may be a novel biomarker predicting prognosis in patients with colon adenocarcinoma. TIICs differed significantly in samples with higher and lower expression of SLC4A4.