Immune cell infiltration-based signature for prognosis and immunogenomic analysis in breast cancer

T cell-related circRNA pairs to predict prognosis of patients with esophageal squamous cell carcinoma

The prognosis for esophageal squamous cell carcinoma (ESCC), a prevalent and aggressive form of cancer, remains poor despite advancements in treatment options. Addressing the gap in comprehensive prognostic information derived from circRNA expression profiles for ESCC, our study aimed to establish a linkage between circRNA expressions and ESCC prognosis. To achieve this, we first developed an optimized prognostic model named T cell-related risk score (TRRS), which integrates T cell-associated features with machine learning algorithms. In parallel, we re-analyzed existing RNA-seq datasets to redefine the expression profiles of circRNAs and mRNAs. Utilizing the TRRS as a foundational "bridge," we identified circRNAs correlated with TRRS, leading to the development of a novel circRNA pair-based prognostic model, the TCRS, which is independent of specific expression levels. Further investigations uncovered two circRNAs, circNLK(5,6,7).1 and circRC3H1(2).1, with potential functional significance. These findings underscore the utility of these risk scores as tools for predicting overall survival and identifying potential therapeutic targets for ESCC patients.

HLF is a promising prognostic, immunological, and therapeutic biomarker in human tumors

Despite past research linking HLF mutations to cancer development, no pan-cancer analyses of HLF have been published. As a result, we utilized multiple databases to illustrate the potential roles of HLF in diverse types of cancers. Several databases were used to assess HLF expression in the TCGA cancer samples. Additional assessments were undertaken to investigate the relationship between HLF and overall survival, immune cell infiltration, genetic alterations, promoter methylation, and protein-protein interaction. HLF's putative roles and the relationship between HLF expression and drug reactivity were investigated. HLF expression was shown to be lower in tumor tissues from a variety of malignancies when compared to normal tissues. There was a substantial link found between HLF expression and patient survival, genetic mutations, and immunological infiltration. HLF influenced the pathways of apoptosis, cell cycle, EMT, and PI3K/AKT signaling. Abnormal expression of HLF lowered sensitivity to numerous anti-tumor drugs and small compounds. According to our findings, reduced HLF expression drives cancer growth, and it has the potential to be identified as a vital biomarker for use in prognosis, immunotherapy, and targeted treatment of a range of malignancies.

Integration of Radiomics and Immune-Related Genes Signatures for Predicting Axillary Lymph Node Metastasis in Breast Cancer

BACKGROUND

To develop a radiogenomics nomogram for predicting axillary lymph node (ALN) metastasis in breast cancer and reveal underlying associations between radiomics features and biological pathways.

MATERIALS AND METHODS

This study included 1062 breast cancer patients, 90 patients with both DCE-MRI and gene expression data. The optimal immune-related genes and radiomics features associated with ALN metastasis were firstly calculated, and corresponding feature signatures were constructed to further validate their performances in predicting ALN metastasis. The radiogenomics nomogram for predicting the risk of ALN metastasis was established by integrating radiomics signature, immune-related genes (IRG) signature, and critical clinicopathological factors. Gene modules associated with key radiomics features were identified by weighted gene co-expression network analysis (WGCNA) and submitted to functional enrichment analysis. Gene set variation analysis (GSVA) and correlation analysis were performed to investigate the associations between radiomics features and biological pathways.

RESULTS

The radiogenomics nomogram showed promising predictive power for predicting ALN metastasis, with AUCs of 0.973 and 0.928 in the training and testing groups, respectively. WGCNA and functional enrichment analysis revealed that gene modules associated with key radiomics features were mainly enriched in breast cancer metastasis-related pathways, such as focal adhesion, ECM-receptor interaction, and cell adhesion molecules. GSVA also identified pathway activities associated with radiomics features such as glycogen synthesis, integration of energy metabolism.

CONCLUSION

The radiogenomics nomogram can serve as an effective tool to predict the risk of ALN metastasis. This study provides further evidence that radiomics phenotypes may be driven by biological pathways related to breast cancer metastasis.

Identification of immune-related tumor antigens and immune subtypes in osteosarcoma

Purpose

The development of tumor vaccines has become a hot topic in immunotherapy for osteosarcoma (OS); however, more tumor antigens with stronger immunogenicity need to be identified.

Methods

We downloaded six sets of gene expression profile data from online databases. The overexpressed genes were analyzed, intersected, and used to calculate the immune infiltration abundance in the TARGET OS dataset based on their expression matrix. Potential tumor antigen genes were identified based on whether they exhibited a high correlation with the antigen-presenting cells (APCs). A total of 1330 immune-related genes (IRGs) from the ImmPort website were retrieved based on their expression, and the Consensus Cluster method was used to obtain immune subtypes of the OS samples. Prognosis, immune microenvironment, and sensitivity to drugs were compared among the immune subtypes.

Results

In total, 680 genes were overexpressed in at least two datasets, of which TREM2, TNFRSF12A, and THY1 were positively correlated with different APCs. Based on the expression matrix of 1330 IRGs in TARGET-OS, two immune subtypes, IS1 and IS2, were identified. The prognosis of the IS1 subtype was better than that of IS2, the expression of immune checkpoint (ICP)-related genes was higher in patients with the IS1 subtype, and immune cell infiltration and sensitivity to 16 drugs were generally higher in IS1 subtype patients.

Conclusion

We identified three APC-correlated genes that can be considered to code for potential novel tumor antigens for OS vaccines. Two immune subtypes in patients with OS were identified to implement personalized treatments using mRNA vaccines.

High expression of CCNB2 is an independent predictive poor prognostic biomarker and correlates with immune infiltrates in breast carcinoma

Background

Cyclin B2 (CCNB2) is associated with cell cycle progression, acting as a cell cycle checkpoint in progression of G2/M transition. In many cancer patients, it has been observed that overexpression of CCNB2 enhances tumor invasiveness and leads to adverse prognosis. However, the association of CCNB2 with the tumor microenvironment remains unclear. Therefore, it is necessary to clarify the associations of CCNB2 with the immune status and prognosis of breast carcinoma (BRCA).

Methods

Gene expression and clinical data for BRCA were obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases, followed by association analyses of CCNB2 expression with prognosis, immune cell infiltration, and immune checkpoints. This study further performed drug sensitivity analysis and constructed a prognostic nomogram for CCNB2.

Results

3619 differentially expressed genes were identified in BRCA, including CCNB2 that emerged as a key gene in the network. High CCNB2 expression correlated with poor prognosis. Functional analysis demonstrated enrichment of CCNB2 co-expressed genes with the cell cycle, cancer progression, cell energy, and immune pathways. Microsatellite instability and tumor mutation burden analyses indicated CCNB2 as a candidate immunotherapy target. Tumor-infiltrating myeloid-derived suppressor cells, regulatory T cells, and T helper 2 cells were associated with CCNB2-related tumor progression and metastasis. CCNB2 expression positively correlated with immune checkpoints, indicating that high CCNB2 expression might facilitate tumor immune escape. Tumors with high CCNB2 expression showed sensitivity to phosphoinositide 3-kinase-protein kinase B-mammalian target of rapamycin and cyclin-dependent kinase (CDK) 4/6 inhibitors, and the nomogram had good prognostic predictive ability for patients with BRCA.

Conclusions

CCNB2 may play a crucial role in tumorigenesis and serve as an independent prognostic biomarker associated with tumor microenvironment, tumor immune infiltration and immunotherapy in BRCA.

Applying image features of proximal paracancerous tissues in predicting prognosis of patients with hepatocellular carcinoma

BACKGROUND

Most of the methods using digital pathological image for predicting Hepatocellular carcinoma (HCC) prognosis have not considered paracancerous tissue microenvironment (PTME), which are potentially important for tumour initiation and metastasis. This study aimed to identify roles of image features of PTME in predicting prognosis and tumour recurrence of HCC patients.

METHODS

We collected whole slide images (WSIs) of 146 HCC patients from Sun Yat-sen Memorial Hospital (SYSM dataset). For each WSI, five types of regions of interests (ROIs) in PTME and tumours were manually annotated. These ROIs were used to construct a Lasso Cox survival model for predicting the prognosis of HCC patients. To make the model broadly useful, we established a deep learning method to automatically segment WSIs, and further used it to construct a prognosis prediction model. This model was tested by the samples of 225 HCC patients from the Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC).

RESULTS

In predicting prognosis of the HCC patients, using the image features of manually annotated ROIs in PTME achieved C-index 0.668 in the SYSM testing dataset, which is higher than the C-index 0.648 reached by the model only using image features of tumours. Integrating ROIs of PTME and tumours achieved C-index 0.693 in the SYSM testing dataset. The model using automatically segmented ROIs of PTME and tumours achieved C-index of 0.665 (95% CI: 0.556-0.774) in the TCGA-LIHC samples, which is better than the widely used methods, WSISA (0.567), DeepGraphSurv (0.593), and SeTranSurv (0.642). Finally, we found the Texture SumAverage Skew HV on immune cell infiltration and Texture related features on desmoplastic reaction are the most important features of PTME in predicting HCC prognosis. We additionally used the model in prediction HCC recurrence for patients from SYSM-training, SYSM-testing, and TCGA-LIHC datasets, indicating the important roles of PTME in the prediction.

CONCLUSIONS

Our results indicate image features of PTME is critical for improving the prognosis prediction of HCC. Moreover, the image features related with immune cell infiltration and desmoplastic reaction of PTME are the most important factors associated with prognosis of HCC.

Multiomics-based classifier to decipher immune landscape of uveal melanoma and predict patient outcomes

Uveal melanoma (UVM) prognosis and the possibilities for targeted therapy depend on a thorough understanding of immune infiltration features and the analysis of genomic and immune signatures. Leveraging multi-omics data from The Cancer Genome Atlas and GEO datasets, we employed an unsupervised clustering algorithm to categorize UVM into immune-related subgroups. Subsequent multi-omics analysis revealed two distinct UVM subtypes, each characterized by unique genomic mutations and immune microenvironment disparities. The aggressive UMCS2 subtype exhibited higher TNM stage and poorer survival, marked by elevated metabolism and increased immune infiltration. However, UMCS2 displayed heightened tumor mutational burden and immune dysfunction, leading to reduced responsiveness to immunotherapy. Importantly, these subtypes demonstrated differential sensitivity to targeted drugs due to significant variances in metabolic and immune environments, with UMCS2 displaying lower sensitivity. We developed a robust, subtype-specific marker-based risk scoring system. This system's diagnostic accuracy was validated through ROC curves, decision curve analysis, and calibration curves, all yielding satisfactory results. Additionally, cell experiments identified the pivotal function of HTR2B, the most crucial factor in this risk model. Knocking down HTR2B significantly reduced the activity, proliferation, and invasion ability of the UVM cell line. These findings underscored the impact of gene and immune microenvironment alterations in driving distinct molecular subtypes, emphasizing the need for precise treatment strategies. The molecular subtyping-based risk assessment system not only aids in predicting patient prognosis but also guides the identification of populations suitable for combined treatment. Molecules represented by HTR2B in the model may serve as effective therapeutic targets for UVM.Communicated by Ramaswamy H. Sarma.

Construction of immune-related molecular diagnostic and predictive models of hepatocellular carcinoma based on machine learning

Background

To exploit hepatocellular carcinoma (HCC) diagnostic substances, we identify potential predictive markers based on machine learning and to explore the significance of immune cell infiltration in this pathology.

Method

Three HCC gene expression datasets were used for weighted gene co-expression network analysis (WGCNA) and differential expression analysis. Least Absolute Shrinkage and Selection Operator (LASSO) and Random Forest were applied to identify candidate biomarkers. The diagnostic value of HCC diagnostic gene biomarkers was further assessed by the area under the ROC curve observed in the validation dataset. CIBERSORT was used to analyze 22 immune cell fractions from HCC patients and to analyze their correlation with diagnostic markers. In addition, the prognostic value of the markers and the sensitivity of the drugs were analyzed.

Result

WGCNA and differential expression analysis were used to screen 396 distinct gene signatures in HCC tissues. They were mostly engaged in cytoplasmic fusion and the cell division cycle, according to gene enrichment analyses. Five genes were shown to have a high diagnostic value for use as diagnostic biomarkers for HCC, including EFHD1 (AUC = 0.77), KIF4A (AUC = 0.97), UBE2C (AUC = 0.96), SMYD3 (AUC = 0.91), and MCM7 (AUC = 0.93). T cells, NK cells, macrophages, and dendritic cells were found to be related to diagnostic markers in HCC tissues by immune cell infiltration analysis, indicating that these cells are intimately linked to the onset and spread of HCC. Concurrently, these five genes and their constructed models have considerable prognostic value.

Conclusion

These five genes (EFHD1, KIF4A, UBE2C, SMYD3, and MCM7) may serve as new candidate molecular markers for HCC, providing new insights for future diagnosis, prognosis, and molecular therapy of HCC.

Identification and characterization of CLEC11A and its derived immune signature in gastric cancer

Introduction

C-type lectin domain family 11 member A (CLEC11A) was characterized as a growth factor that mainly regulates hematopoietic function and differentiation of bone cells. However, the involvement of CLEC11A in gastric cancer (GC) is not well understood.

Methods

Transcriptomic data and clinical information pertaining to GC were obtained and analyzed from publicly available databases. The relationships between CLEC11A and prognoses, genetic alterations, tumor microenvironment (TME), and therapeutic responses in GC patients were analyzed by bioinformatics methods. A CLEC11A-derived immune signature was developed and validated, and its mutational landscapes, immunological characteristics as well as drug sensitivities were explored. A nomogram was established by combining CLEC11A-derived immune signature and clinical factors. The expression and carcinogenic effects of CLEC11A in GC were verified by qRT-PCR, cell migration, invasion, cell cycle analysis, and in vivo model analysis. Myeloid-derived suppressor cells (MDSCs), regulatory T cells (Tregs), M2 macrophages, and T cells in tumor samples extracted from mice were analyzed utilizing flow cytometry analysis.

Results

CLEC11A was over-expressed in GC, and the elevated CLEC11A expression indicated an unfavorable prognosis in GC patients. CLEC11A was involved in genomic alterations and associated with the TME in GC. Moreover, elevated CLEC11A was found to reduce the benefit of immunotherapy according to immunophenoscore (IPS) and the tumor immune dysfunction, exclusion (TIDE). After validation, the CLEC11A-derived immune signature demonstrated a consistent ability to predict the survival outcomes in GC patients. A nomogram that quantifies survival probability was constructed to improve the accuracy of prognosis prediction in GC patients. Using shRNA to suppress the expression of CLEC11A led to significant inhibitions of cell cycle progression, migration, and invasion, as well as a marked reduction of in vivo tumor growth. Moreover, the flow cytometry assay showed that the knock-down of CLEC11A increased the infiltration of cytotoxic CD8+ T cells and helper CD4+ T into tumors while decreasing the percentage of M2 macrophages, MDSCs, and Tregs.

Conclusion

Collectively, our findings revealed that CLEC11A could be a prognostic and immunological biomarker in GC, and CLEC11A-derived immune signature might serve as a new option for clinicians to predict outcomes and formulate personalized treatment plans for GC patients.

Identification and validation of methylation-CpG prognostic signature for prognosis of hepatocellular carcinoma

Epigenetic biomarkers help predict the prognosis of cancer patients and evaluating the clinical outcome of immunization therapy. In this study, we present a personalized gene methylation-CpG signature to enhance the accuracy of survival prediction for individuals with hepatocellular carcinoma (HCC). Utilizing RNA sequencing and methylation datasets from GEO as well as TCGA, we conducted single sample GSEA (ssGSEA), WGCNA, as well as Cox regression. Through these analyses, we identified 175 oxidative stress and immune-related genes along with 4 CpG loci that are associated with the prognosis of HCC. Subsequently, we constructed a prognostic signature for HCC utilizing these 4 CpG sites, referred to as the HCC Prognostic Signature of Methylation-CpG sites (HPSM). Further investigation revealed an enrichment of immune-related signal pathways in the HPSM-low group, which demonstrated a positive correlation with better survival among HCC patients. Moreover, the methylation of the CpG sites in HPSM was found to be closely linked to drug sensitivity. In vitro experiments tentatively confirmed that promoter methylation regulated the expression of BMPER, one of the CpG sites within HPSM. The expression of BMPER was significantly correlated with cell death in the oxidative stress pathway, and overexpression of BMPER effectively inhibited HCC cell proliferation. Consequently, our findings suggest that HPSM is an independent predictive factor and holds promise for accurately predicting the prognosis of HCC patients.

The tertiary lymphoid structure-related signature identified PTGDS in regulating PD-L1 and promoting the proliferation and migration of glioblastoma

Background

Tertiary lymphoid structure (TLS) is a unique organ that carries out tumor cell elimination at tumor sites. It is continuously stimulated by inflammatory tumor signals and has been found to augment immunotherapy response. However, the detailed mechanisms behind it still need to be defined.

Methods

To explore and grasp the whole picture of TLS from a pan-cancer view, we collected nine TLS-related genes from previous studies. We performed a comprehensive analysis of 9637 samples across 33 tumor types accessed from The Cancer Genome Atlas (TCGA) database. EdU, Transwell, and flow cytometry were performed on the feature gene PTGDS in U251 cells. The regulatory role of PTGDS on PD-L1 expression and macrophage polarization was verified.

Results

Alteration analysis showed that mutations of TLS-related genes were widespread and relatively high. Clustering analysis based on the expression of these nine genes obtained two distinct clusters, with high EIF1AY and PTGDS in cluster 2 and better overall survival in cluster 1. To distinguish the two clusters, we utilized six machine learning algorithms and filtrated EIF1AY, PTGDS, SKAP1, and RBP5 as the characteristic genes, among which the former two genes were proved to be hazardous. PTGDS was found to regulate PD-L1 expression and also promoted the proliferation and migration of U251 cells. The knockdown of PTGDS could reduce the migration of macrophages and inhibit the polarization of macrophages into M2-phenotype. In addition, we established a TLS score to demonstrate patients' TLS activity. The low TLS-score group overlapped with cluster 1 and displayed a better prognosis. Besides, the low TLS-score group was related to better immunotherapy responses. The HE staining of histopathological sections confirmed that the low TLS-score group exhibited higher infiltration of immune cells.

Conclusion

This study reveals broad molecular, tumorigenic, and immunogenic signatures for further functional and therapeutic studies of tertiary lymphoid structure. The TLS score we established effectively predicted immunotherapy response and patients' survival. Its future application and combination await more research.

KNSTRN Is a Prognostic Biomarker That Is Correlated with Immune Infiltration in Breast Cancer and Promotes Cell Cycle and Proliferation

Kinetochore-localized astrin/SPAG5-binding protein (KNSTRN) promotes the progression of bladder cancer and lung adenocarcinoma. However, its expression and biological function in breast cancer remain largely unknown. Therefore, this study aimed to analyze KNSTRN expression, prognoses, correlation with immune infiltration, expression-associated genes, and regulated signaling pathways to characterize its role in regulating the cell cycle using both bioinformatics and in vitro functional experiments. Analyses of The Cancer Genome Atlas, Gene Expression Omnibus, TIMER, and The Human Protein Atlas databases revealed a significant upregulation of KNSTRN transcript and protein levels in breast cancer. Kaplan-Meier survival analyses demonstrated a significant association between high expression of KNSTRN and poor overall survival, relapse-free survival, post-progression survival, and distant metastases-free survival in patients with breast cancer. Furthermore, multivariate Cox regression analyses confirmed that KNSTRN is an independent prognostic factor for breast cancer. Immune infiltration analysis indicated a positive correlation between KNSTRN expression and T regulatory cell infiltration while showing a negative correlation with Tgd and natural killer cell infiltration. Gene set enrichment analysis along with single-cell transcriptome data analysis suggested that KNSTRN promoted cell cycle progression by regulating the expression of key cell cycle proteins. The overexpression and silencing of KNSTRN in vitro, respectively, promoted and inhibited the proliferation of breast cancer cells. The overexpression of KNSTRN enhanced the expression of key cell cycle regulators, including CDK4, CDK6, and cyclin D3, thereby accelerating the G1/S phase transition and leading to aberrant proliferation of breast cancer cells. In conclusion, our study demonstrates that KNSTRN functions as an oncogene in breast cancer by regulating immune response, promoting G1/S transition, and facilitating breast cancer cell proliferation. Moreover, KNSTRN has potential as a molecular biomarker for diagnostic and prognostic prediction in breast cancer.

Identification of genes and pathways associated with menopausal status in breast cancer patients using two algorithms

BACKGROUND

Menopausal status has a known relationship with the levels of estrogen, progesterone, and other sex hormones, potentially influencing the activity of ER, PR, and many other signaling pathways involved in the initiation and progression of breast cancer. However, the differences between premenopausal and postmenopausal breast cancer patients at the molecular level are unclear.

METHODS

We retrieved eight datasets from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) associated with menopausal status in breast cancer patients were identified using the MAMA and LIMMA methods. Based on these validated DEGs, we performed Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein-protein interaction (PPI) networks were constructed. We used DrugBank data to investigate which of these validated DEGs are targetable. Survival analysis was performed to explore the influence of these genes on breast cancer patient prognosis.

RESULTS

We identified 762 DEGs associated with menopausal status in breast cancer patients. PPI network analysis indicated that these genes are primarily involved in pathways such as the cell cycle, oocyte meiosis and progesterone-mediated oocyte maturation pathways. Notably, several genes played roles in multiple signaling pathways and were associated with patient survival. These genes were also observed to be targetable according to the DrugBank database.

CONCLUSION

We identified DEGs associated with menopausal status in breast cancer patients. The association of these genes with several key pathways may promote understanding of the complex characterizations of breast cancer. Our findings offer valuable insights for developing new therapeutic strategies tailored to the menopausal status of breast cancer patients.

9S1R nullomer peptide induces mitochondrial pathology, metabolic suppression, and enhanced immune cell infiltration, in triple-negative breast cancer mouse model

Nullomers are the shortest strings of absent amino acid (aa) sequences in a species or group of species. Primes are those nullomers that have not been detected in the genome of any species. 9S1R is a 5-aa peptide prime sequence attached to 5-arginine aa, used to treat triple negative breast cancer (TNBC) in an in vivo mouse model. This unique peptide, administered with a trehalose carrier (9S1R-NulloPT), offers enhanced solubility and exhibits distinct anti-cancer effects against TNBC. In our study, we investigated the effect of 9S1R-NulloPT on tumor growth, metabolism, metastatic burden, tumor immune-microenvironment (TME), and transcriptome of aggressive mouse TNBC tumors. Notably, treated mice had smaller tumors in the initial phase of the treatment, as compared to untreated control, and diminished in vivo and ex vivo bioluminescence at later-stages - indicative of metabolically quiescent, dying tumors. The treatment also caused changes in TME with increased infiltration of immune cells and altered tumor transcriptome, with 365 upregulated genes and 710 downregulated genes. Consistent with in vitro data, downregulated genes were enriched in cellular metabolic processes (179), specifically mitochondrial TCA cycle/oxidative phosphorylation (44), and translation machinery/ribosome biogenesis (45). The upregulated genes were associated with the developmental (13), ECM organization (12) and focal adhesion pathways (7). In conclusion, our study demonstrates that 9S1R-NulloPT effectively reduced tumor growth during its initial phase, altering the TME and tumor transcriptome. The treatment induced mitochondrial pathology which led to a metabolic deceleration in tumors, aligning with in vitro observations.

MIRS: An AI scoring system for predicting the prognosis and therapy of breast cancer

Tumor-infiltrating immune cells (TIICs) and metastasis are crucial characteristics for tumorigenesis. However, the potential role of their combination in breast cancer (BRCA) remains elusive. Herein, on the basis of quantifying TIICs and tumor metastasis together, we established a precise prognostic scoring system named metastatic and immunogenomic risk score (MIRS) using a neural network model. MIRS showed better performance when compared with other published signatures. MIRS stratifies patients into a high risk subtype (MIRShigh) and a low risk subtype (MIRSlow). The MIRShigh patients exhibit significantly lower survival rate compared with MIRSlow patients (P < 0.0001 ), higher response to chemotherapy, but lower response to immunotherapy. Conversely, higher infiltration level of TIICs and significantly prolonged survival (P = 0.029 ) are observed in MIRSlow patients, indicating sensitive response in immunotherapy. This work presents a promising indicator to guide treatment options of the BRCA population and provides a predicted webtool that is almost universally applicable to BRCA patients.

Identification of disulfidptosis-related subtypes, characterization of tumor microenvironment infiltration, and development of a prognosis model in colorectal cancer

BACKGROUND

Colorectal cancer is the second leading cause of cancer-related deaths, which imposes a significant societal burden. Regular screening and emerging molecular tumor markers have important implications for detecting the progression and development of colorectal cancer. Disulfidptosis is a newly defined type of programmed cell death triggered by abnormal accumulation of disulfide compounds in cells that stimulate disulfide stress. Currently, there is no relevant discussion on this mechanism and colorectal cancer.

METHODS

We classified the disulfidptosis-related subtypes of colorectal cancer using bioinformatics methods. Through secondary clustering of differentially expressed genes between subtypes, we identified characteristic genes of the disulfidptosis subtype, constructed a prognostic model, and searched for potential biomarkers through clinical validation.

RESULTS

Using disulfidptosis-related genes collected from the literature, we classified colorectal cancer patients from public databases into three subtypes. The differentially expressed genes between subtypes were clustered into three gene subtypes, and eight characteristic genes were screened to construct a prognostic model.

CONCLUSION

The disulfidptosis mechanism has important value in the classification of colorectal cancer patients, and characteristic genes selected based on this mechanism can serve as a new potential biological marker for colorectal cancer.

Identification of a novel competing endogenous RNA network and candidate drugs associated with ferroptosis in aldosterone-producing adenomas

Aldosterone-producing adenoma (APA), characterized by unilaterally excessive aldosterone production, is a common cause of primary aldosteronism. Ferroptosis, a recently raised iron-dependent mode of programmed cell death, has been involved in the development and therapy of various diseases. This study obtained datasets of the mRNA and lncRNA expression profiles for APA and adjacent adrenal gland (AAG) from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and lncRNAs (DE lncRNAs) associated with ferroptosis were identified. Enrichment analyses indicated 89 ferroptosis-related DEGs were primarily enriched in ROS related processes and ferroptosis. Two physical cores, and one combined core were identified in the protein-protein interaction (PPI). DEGs and clinical traits were used in conjunction to screen eight hub genes from two hub modules and 89 DEGs. A competitive endogenous RNA (ceRNA) network was constructed via co-express analysis. Thereafter, molecular docking was used to identify potential targets. Two active compounds, QL-X-138 and MK-1775, bound to AURKA and DUOX1, respectively, with the lowest binding energies. Molecular dynamics simulation verified the stability of the two complexes. In summary, our studies identified eight hub genes and a novel ceRNA regulatory network associated with ferroptosis, wherein QL-X-138 and MK-1775 were considered to be potential drugs for treating APA.

Nullomer peptide increases immune cell infiltration and reduces tumor metabolism in triple negative breast cancer mouse model

Background

Nullomers are the shortest strings of absent amino acid (aa) sequences in a species or group of species. Primes are those nullomers that have not been detected in the genome of any species. 9S1R is a 5-aa peptide derived from a prime sequence that is tagged with 5 arginine aa, used to treat triple negative breast cancer (TNBC) in an in vivo TNBC mouse model. 9S1R is administered in trehalose (9S1R-NulloPT), which enhances solubility and exhibits some independent effects against tumor growth and is thus an important component in the drug preparation.

Method

We examined the effect of 9S1R-NulloPT on tumor growth, metabolism, metastatic burden, necrosis, tumor immune microenvironment, and the transcriptome of aggressive mouse TNBC tumors.

Results

The peptide-treated mice had smaller tumors in the initial phase of the treatment, as compared to the untreated control, and reduced in vivo bioluminescence at later stages, which is indicative of metabolically inactive tumors. A decrease in ex vivo bioluminescence was also observed in the excised tumors of treated mice, but not in the secondary metastasis in the lungs. The treatment also caused changes in tumor immune microenvironment with increased infiltration of immune cells and margin inflammation. The treatment upregulated 365 genes and downregulated 710 genes in tumors compared to the untreated group. Consistent with in vitro findings in breast cancer cell lines, downregulated genes in the treated TNBC tumors include Cellular Metabolic Process Related genes (179), specifically mitochondrial genes associated with TCA cycle/oxidative phosphorylation (44), and translation machinery/ribosome biogenesis genes (45). Among upregulated genes, the Developmental Pathway (13), ECM Organization (12) and Focal Adhesion Related Pathways (7) were noteworthy. We also present data from a pilot study using a bilateral BC mouse model, which supports our findings.

Conclusion

In conclusion, although 9S1R-NulloPT was moderate at reducing the tumor volume, it altered the tumor immune microenvironment as well as the tumor transcriptome, rendering tumors metabolically less active by downregulating the mitochondrial function and ribosome biogenesis. This corroborates previously published in vitro findings.

Regional Variation in the Tumor Microenvironment, Immune Escape and Prognostic Factors in Breast Cancer in Sub-Saharan Africa

The low overall survival rates of patients with breast cancer in sub-Saharan Africa (SSA) are driven by regionally differing tumor biology, advanced tumor stages at diagnosis, and limited access to therapy. However, it is not known whether regional differences in the composition of the tumor microenvironment (TME) exist and affect patients' prognosis. In this international, multicentre cohort study, 1,237 formalin-fixed, paraffin-embedded breast cancer samples, including samples of the "African Breast Cancer-Disparities in Outcomes (ABC-DO) Study," were analyzed. The immune cell phenotypes, their spatial distribution in the TME, and immune escape mechanisms of breast cancer samples from SSA and Germany (n = 117) were investigated using histomorphology, conventional and multiplex IHC, and RNA expression analysis. The data revealed no regional differences in the number of tumor-infiltrating lymphocytes (TIL) in the 1,237 SSA breast cancer samples, while the distribution of TILs in different breast cancer IHC subtypes showed regional diversity, particularly when compared with German samples. Higher TIL densities were associated with better survival in the SSA cohort (n = 400), but regional differences concerning the predictive value of TILs existed. High numbers of CD163+ macrophages and CD3+CD8+ T cells accompanied by reduced cytotoxicity, altered IL10 and IFNγ levels and downregulation of MHC class I components were predominantly detected in breast cancer samples from Western SSA. Features of nonimmunogenic breast cancer phenotypes were associated with reduced patient survival (n = 131). We therefore conclude that regional diversity in the distribution of breast cancer subtypes, TME composition, and immune escape mechanisms should be considered for therapy decisions in SSA and the design of personalized therapies. See related Spotlight by Bergin et al., p. 705.

Bioinformatics-Based Identification of CircRNA-MicroRNA-mRNA Network for Calcific Aortic Valve Disease

Background

Calcific aortic valve disease (CAVD) is the most common native valve disease. Valvular interstitial cell (VIC) osteogenic differentiation and valvular endothelial cell (VEC) dysfunction are key steps in CAVD progression. Circular RNA (circRNAs) is involved in regulating osteogenic differentiation with mesenchymal cells and is associated with multiple disease progression, but the function of circRNAs in CAVD remains unknown. Here, we aimed to investigate the effect and potential significance of circRNA-miRNA-mRNA networks in CAVD.

Methods

Two mRNA datasets, one miRNA dataset, and one circRNA dataset of CAVD downloaded from GEO were used to identify DE-circRNAs, DE-miRNAs, and DE-mRNAs. Based on the online website prediction function, the common mRNAs (FmRNAs) for constructing circRNA-miRNA-mRNA networks were identified. GO and KEGG enrichment analyses were performed on FmRNAs. In addition, hub genes were identified by PPI networks. Based on the expression of each data set, the circRNA-miRNA-hub gene network was constructed by Cytoscape (version 3.6.1).

Results

32 DE-circRNAs, 206 DE-miRNAs, and 2170 DE-mRNAs were identified. Fifty-nine FmRNAs were obtained by intersection. The KEGG pathway analysis of FmRNAs was enriched in pathways in cancer, JAK-STAT signaling pathway, cell cycle, and MAPK signaling pathway. Meanwhile, transcription, nucleolus, and protein homodimerization activity were significantly enriched in GO analysis. Eight hub genes were identified based on the PPI network. Three possible regulatory networks in CAVD disease were obtained based on the biological functions of circRNAs including: hsa_circ_0026817-hsa-miR-211-5p-CACNA1C, hsa_circ_0007215-hsa-miR-1252-5p-MECP2, and hsa_circ_0007215-hsa-miR-1343-3p- RBL1.

Conclusion

The present bionformatics analysis suggests the functional effect for the circRNA-miRNA-mRNA network in CAVD pathogenesis and provides new targets for therapeutics.

The transcriptional risk scores for kidney renal clear cell carcinoma using XGBoost and multiple omics data

Most kidney cancers are kidney renal clear cell carcinoma (KIRC) that is a main cause of cancer-related deaths. Polygenic risk score (PRS) is a weighted linear combination of phenotypic related alleles on the genome that can be used to assess KIRC risk. However, standalone SNP data as input to the PRS model may not provide satisfactory result. Therefore, Transcriptional risk scores (TRS) based on multi-omics data and machine learning models were proposed to assess the risk of KIRC. First, we collected four types of multi-omics data (DNA methylation, miRNA, mRNA and lncRNA) of KIRC patients from the TCGA database. Subsequently, a novel TRS method utilizing multiple omics data and XGBoost model was developed. Finally, we performed prevalence analysis and prognosis prediction to evaluate the utility of the TRS generated by our method. Our TRS methods exhibited better predictive performance than the linear models and other machine learning models. Furthermore, the prediction accuracy of combined TRS model was higher than that of single-omics TRS model. The KM curves showed that TRS was a valid prognostic indicator for cancer staging. Our proposed method extended the current definition of TRS from standalone SNP data to multi-omics data and was superior to the linear models and other machine learning models, which may provide a useful implement for diagnostic and prognostic prediction of KIRC.

A generalized non-linear model predicting efficacy of neoadjuvant therapy in HER2+ breast cancer

Neoadjuvant therapy (NAT) is currently recommended to patients with human epidermal growth factor receptor 2-positive breast cancer (HER2+ BC) that typically exhibit a poor prognosis. The tumor immune microenvironment profoundly affects the efficacy of NAT. However, the correlation between tumor-infiltrating lymphocytes or their specific subpopulations and the response to NAT in HER2+ BC remains largely unknown. In our study, the immune infiltration status of 295 patients was classified as "immune-rich" or "immune-poor" phenotypes. The "immune-rich" phenotype was significantly positively related to pathological complete response (pCR). Ten genes were correlated with both pCR and the immune phenotype based on the results of spline and logistic regression. We constructed a generalized non-linear model combining linear and non-linear gene effects and successfully validated its predictive power using an internal and external validation set (AUC = 0.819, 0.797; respectively) and a clinical set (accuracy = 0.75).

Single-Cell Transcriptome Analysis Revealed Heterogeneity and Identified Novel Therapeutic Targets for Breast Cancer Subtypes

Breast cancer (BC) is a heterogeneous disease, which is primarily classified according to hormone receptors and HER2 expression. Despite the many advances in BC diagnosis and management, the identification of novel actionable therapeutic targets expressed by cancerous cells has always been a daunting task due to the large heterogeneity of the disease and the presence of non-cancerous cells (i.e., immune cells and stromal cells) within the tumor microenvironment. In the current study, we employed computational algorithms to decipher the cellular composition of estrogen receptor-positive (ER⁺), HER2⁺, ER⁺HER2⁺, and triple-negative BC (TNBC) subtypes from a total of 49,899 single cells' publicly available transcriptomic data derived from 26 BC patients. Restricting the analysis to EPCAM⁺Lin^- tumor epithelial cells, we identified the enriched gene sets in each BC molecular subtype. Integration of single-cell transcriptomic with CRISPR-Cas9 functional screen data identified 13 potential therapeutic targets for ER⁺, 44 potential therapeutic targets for HER2⁺, and 29 potential therapeutic targets for TNBC. Interestingly, several of the identified therapeutic targets outperformed the current standard of care for each BC subtype. Given the aggressive nature and lack of targeted therapies for TNBC, elevated expression of ENO1, FDPS, CCT6A, TUBB2A, and PGK1 predicted worse relapse-free survival (RFS) in basal BC (n = 442), while elevated expression of ENO1, FDPS, CCT6A, and PGK1 was observed in the most aggressive BLIS TNBC subtype. Mechanistically, targeted depletion of ENO1 and FDPS halted TNBC cell proliferation, colony formation, and organoid tumor growth under 3-dimensional conditions and increased cell death, suggesting their potential use as novel therapeutic targets for TNBC. Differential expression and gene set enrichment analysis in TNBC revealed enrichment in the cycle and mitosis functional categories in FDPS^high, while ENO1^high was associated with numerous functional categories, including cell cycle, glycolysis, and ATP metabolic processes. Taken together, our data are the first to unravel the unique gene signatures and to identify novel dependencies and therapeutic vulnerabilities for each BC molecular subtype, thus setting the foundation for the future development of more effective targeted therapies for BC.

A computational framework of routine test data for the cost-effective chronic disease prediction

Chronic diseases, because of insidious onset and long latent period, have become the major global disease burden. However, the current chronic disease diagnosis methods based on genetic markers or imaging analysis are challenging to promote completely due to high costs and cannot reach universality and popularization. This study analyzed massive data from routine blood and biochemical test of 32 448 patients and developed a novel framework for cost-effective chronic disease prediction with high accuracy (AUC 87.32%). Based on the best-performing XGBoost algorithm, 20 classification models were further constructed for 17 types of chronic diseases, including 9 types of cancers, 5 types of cardiovascular diseases and 3 types of mental illness. The highest accuracy of the model was 90.13% for cardia cancer, and the lowest was 76.38% for rectal cancer. The model interpretation with the SHAP algorithm showed that CREA, R-CV, GLU and NEUT% might be important indices to identify the most chronic diseases. PDW and R-CV are also discovered to be crucial indices in classifying the three types of chronic diseases (cardiovascular disease, cancer and mental illness). In addition, R-CV has a higher specificity for cancer, ALP for cardiovascular disease and GLU for mental illness. The association between chronic diseases was further revealed. At last, we build a user-friendly explainable machine-learning-based clinical decision support system (DisPioneer: http://bioinfor.imu.edu.cn/dispioneer) to assist in predicting, classifying and treating chronic diseases. This cost-effective work with simple blood tests will benefit more people and motivate clinical implementation and further investigation of chronic diseases prevention and surveillance program.

NAD+ biosynthesis metabolism predicts prognosis and indicates immune microenvironment for breast cancer

The growing evidence implies that tumor cells need to increase NAD+ levels by upregulating NAD+ biosynthesis to satisfy their growth demand. NAD+ biosynthesis metabolism is implicated in tumor progression. Breast cancer (BC) is the most common malignant malignancy in the world. Nevertheless, the prognostic significance of NAD+ biosynthesis and its relationship with the tumor immune microenvironment in breast cancer still need further investigation. In this study, we obtained the mRNA expression data and clinical information of BC samples from public databases and calculated the level of NAD+ biosynthesis activity by single-sample gene set enrichment analysis (ssGSEA). We then explored the relationship between the NAD+ biosynthesis score, infiltrating immune cells, prognosis significance, immunogenicity and immune checkpoint molecules. The results demonstrated that patients with high NAD+ biosynthetic score displayed poor prognosis, high immune infiltration, high immunogenicity, elevated PD-L1 expression, and might more benefit from immunotherapy. Taken together, our studies not only deepened the understanding of NAD+ biosynthesis metabolism of breast cancer but also provided new insights into personalized treatment strategies and immunological therapies to improve the outcomes of breast cancer patients.

Risk scoring based on DNA methylation-driven related DEGs for colorectal cancer prognosis with systematic insights

Colorectal cancer is a common malignant tumor of the digestive tract. Despite advances in diagnostic techniques and medications. Its prognosis remains challenging. DNA methylation-driven related circulating tumor cells have attracted enormous interest in diagnosing owing to their non-invasive nature and early recognition properties. However, the mechanism through which risk biomarkers act remains elusive. Here, we designed a risk model based on differentially expressed genes, DNA methylation, robust, and survival-related factors in the framework of Cox regression. The model has satisfactory performance and is independently verified by an external and isolated dataset in terms of C-index value, ROC, and tROC. The model was applied to Colorectal cancer patients who were subsequently divided into high- and low-risk groups. Functional annotations, genomic alterations, tumor immune environment, and drug sensitivity were analyzed. We observed that up-regulated genes are associated with epithelial cell differentiation and MAPK signaling pathways. The down-regulated genes are related to IL-7 signaling and apoptosis-induced DNA fragmentation. Interestingly, the immune system was inhibited in high-risk groups. High-frequency mutation genes tend to co-occur. High-risk score patients are related to copy number amplification events. To address the challenges, we suggested eleven and twenty-one drugs that are sensitive to low- and high-risk patients. Finally, an artificial neural network was provided to evaluate the immunotherapeutic efficiency. Taken together, the findings demonstrated that our risk score model is robust and reliable for evaluating the prognosis with novel diagnostic and treatment targets. It also yields benefits for the treatment and provides unique insights into developing therapeutic strategies.

Profiling specific cell populations within the inflammatory tumor microenvironment by oscillating-gradient diffusion-weighted MRI

BACKGROUND

The inflammatory tumor microenvironment (TME) is formed by various immune cells, being closely associated with tumorigenesis. Especially, the interaction between tumor-infiltrating T-cells and macrophages has a crucial impact on tumor progression and metastatic spread. The purpose of this study was to investigate whether oscillating-gradient diffusion-weighted MRI (OGSE-DWI) enables a cell size-based discrimination between different cell populations of the TME.

METHODS

Sine-shaped OGSE-DWI was combined with the Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion (IMPULSED) approach to measure microscale diffusion distances, here relating to cell sizes. The accuracy of IMPULSED-derived cell radii was evaluated using in vitro spheroid models, consisting of either pure cancer cells, macrophages, or T-cells. Subsequently, in vivo experiments aimed to assess changes within the TME and its specific immune cell composition in syngeneic murine breast cancer models with divergent degrees of malignancy (4T1, 67NR) during tumor progression, clodronate liposome-mediated depletion of macrophages, and immune checkpoint inhibitor (ICI) treatment. Ex vivo analysis of IMPULSED-derived cell radii was conducted by immunohistochemical wheat germ agglutinin staining of cell membranes, while intratumoral immune cell composition was analyzed by CD3 and F4/80 co-staining.

RESULTS

OGSE-DWI detected mean cell radii of 8.8±1.3 µm for 4T1, 8.2±1.4 µm for 67NR, 13.0±1.7 for macrophage, and 3.8±1.8 µm for T-cell spheroids. While T-cell infiltration during progression of 4T1 tumors was observed by decreasing mean cell radii from 9.7±1.0 to 5.0±1.5 µm, increasing amount of intratumoral macrophages during progression of 67NR tumors resulted in increasing mean cell radii from 8.9±1.2 to 12.5±1.1 µm. After macrophage depletion, mean cell radii decreased from 6.3±1.7 to 4.4±0.5 µm. T-cell infiltration after ICI treatment was captured by decreasing mean cell radii in both tumor models, with more pronounced effects in the 67NR tumor model.

CONCLUSIONS

OGSE-DWI provides a versatile tool for non-invasive profiling of the inflammatory TME by assessing the dominating cell type T-cells or macrophages.

Glioblastoma glycolytic signature predicts unfavorable prognosis, immunological heterogeneity, and ENO1 promotes microglia M2 polarization and cancer cell malignancy

Glioblastomas are the most malignant brain tumors, whose progress was promoted by aberrate aerobic glycolysis. The immune environment was highly engaged in glioblastoma formation, while its interaction with aerobic glycolysis remained unclear. Herein, we build a 7-gene Glycolytic Score (GS) by Elastic Net in the training set and two independent validating sets. The GS predicted malignant features and poor survival with good performances. Immune functional analyses and Cibersort calculation identified depressed T cells, B cells, natural killer cells immunity, and high immunosuppressive cell infiltration in the high-GS group. Also, high expressions of the immune-escape genes were discovered. Subsequently, the single-cell analyses validated the glycolysis-related immunosuppression. The functional results manifested the high-GS neoplastic cells' association with T cells, NK cells, and macrophage function regulation. The intercellular cross-talk showed strong associations between high-GS neoplastic cells and M2 macrophages/microglia in several immunological pathways. We finally confirmed that ENO1, the key gene of the GS, promoted M2 microglia polarization and glioblastoma cell malignant behaviors via immunofluorescence, clone formation, CCK8, and transwell rescue experiments. These results indicated the interactions between cancerous glycolysis and immunosuppression and glycolysis' role in promoting glioblastoma progression. Conclusively, we built a robust model and discovered strong interaction between GS and immune, shedding light on prognosis management improvement and therapeutic strategies development for glioblastoma patients.

MARCO is a potential prognostic and immunotherapy biomarker

BACKGROUND

Macrophage receptor with collagenous structure (MARCO), a novel immune checkpoint expressed on tumor-associated macrophages, has antitumor therapeutic properties. However, the association between MARCO and patient prognosis, immune infiltration, and ICI immunotherapy needs to be studied urgently.

METHODS

MARCO distribution in cancer tissues was investigated using the TCGA and GTEx databases. The PrognoScan and KM Plotter databases was used to assess the MARCO prognosis. TIMER2.0, GEPIA, cBioPortal, and GSEA all confirmed the link between MARCO and immune infiltration, mutation profile, and enrichment pathway analysis. Data visualization was implemented by R language.

RESULTS

In general, MARCO had a substantial impact on the prognosis of cancer patients and was expressed differently in cancer and adjacent normal tissues. High expression of MARCO was associated with poorer OS in bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), lung squamous cell carcinoma (LUSC), colon adenocarcinoma (COAD), and prostate adenocarcinoma (PRAD). However, high expression of MARCO had a better PFI in brain lower-grade glioma (LGG) and skin cutaneous melanoma (SKCM). We discovered that MARCO expression was lowest in pancreatic adenocarcinoma (PAAD) and rectum adenocarcinoma (READ) stage 1, BLCA stage 2, LUSC and stomach adenocarcinoma (STAD) stage 3, and liver hepatocellular carcinoma (LIHC) stage 4. Subsequently, we analyzed the correlation between MARCO and 47 immune checkpoints and observed that MARCO was positively connected with CD80, CD86, and leukocyte-associated immunoglobulin-like receptor 1(LAIR1) in most cancers. In COAD, MARCO has the most microsatellite instability (MSI). In addition, we discovered that high expression of MARCO patients had a better prognosis after immune checkpoint inhibitor (ICI) treatment in SKCM. Finally, GSEA revealed a significant correlation between MARCO and TNF/NFκB signaling, KRAS signaling, PI3K/AKT/mTOR pathway, IL-6-STAT3 signaling, TGFβ pathway, and p53 pathway.

CONCLUSION

This study comprehensively investigated the relationship between MARCO and clinical prognosis, immune infiltration, and ICI immunotherapy in various cancers. We demonstrated the potential of MARCO as an emerging biomarker, exploring new avenues for future tumor immunotherapy.

Prognostic Significance of the CXCLs and Its Impact on the Immune Microenvironment in Ovarian Cancer

The chemokine (C-X-C motif) ligand (CXCL) family in tumor tissue is closely related to tumor growth, metastasis, and survival. However, the differential expression profile and prognostic value of the CXCLs in ovarian cancer (OC) have not been elucidated. Therefore, we studied the expression levels and mutations of CXCLs in OC patient in TCGA and various public databases. The expression differences of CXCLs in OC cancer tissues and normal tissues were compared through the Gene Expression Profiling Interactive Analysis (GEPIA) database. The effect of CXCLs on OC prognosis was analyzed using the Kaplan-Meier curves in GEPIA database. The impact of CXCLs on immune infiltration and clinicopathological outcomes in OC was assessed using the TIMER algorithm. Compared with normal tissues, we found that eight CXCLs were significantly differentially expressed in OC. The expression levels of CXCL9 (P = 0.0201), CXCL11 (P = 0.0385), and CXCL13 (P = 0.0288) were significantly associated with tumor stage. CXCL13 was the only gene that significantly affected both disease-free survival (DFS) and overall survival (OS) in OC, and higher CXCL13 transcript levels implied longer DFS and OS. Although there was no significant impact on DFS, CXCL10 (P = 0.0079) and CXCL11 (P = 0.0011) expression levels had a significant effect on OS in OC. At the same time, CXCLs were significantly associated with several immune-infiltrating cells in OC tissues. The CXCLs were significantly associated with one or more immune-infiltrating cells in OC tissue. CXCL13 was differentially expressed in OC and significantly affected the prognosis of patients and was a potential marker of OC prognosis.

Prediction of immune subtypes and overall survival in lung squamous cell carcinoma

OBJECTIVE

Lung squamous cell carcinoma (LUSC), one of the most common subtypes of lung cancer, is a leading cause of cancer-caused deaths in the world. It has been well demonstrated that a deep understanding of the tumor environment in cancer would be helpful to predict the prognosis of patients. This study aimed to evaluate the tumor environment in LUSC, and to construct an efficient prognosis model involved in specific subtypes.

METHODS

Four expression files were downloaded from the Gene Expression Omnibus (GEO) database. Three datasets (GSE19188, GSE2088, GSE6044) were considered as the testing group and the other dataset (GSE11969) was used as the validation group. By performing LUSC immune subtype consensus clustering (CC), LUSC patients were separated into two immune subtypes comprising subtype 1 (S1) and subtype 2 (S2). Weighted gene co-expression network (WGCNA) and least absolute shrinkage and selection operator (LASSO) were performed to identify and narrow down the key genes among subtype 1 related genes that were closely related to the overall survival (OS) of LUSC patients. Using immune subtype related genes, a prognostic model was also constructed to predict the OS of LUSC patients.

RESULTS

It showed that LUSC patients in the S1 immune subtype exhibited a better OS than in the S2 immune subtype. WGCNA and LASSO analyses screened out important immune subtype related genes in specific modules that were closely associated with LUSC prognosis, followed by construction of the prognostic model. Both the testing datasets and validation dataset confirmed that the prognostic model could be efficiently used to predict the OS of LUSC patients in subtype 1.

CONCLUSION

We explored the tumor environment in LUSC and established a risk prognostic model that might have the potential to be applied in clinical practice.

Cell Pair Algorithm-Based Immune Infiltrating Cell Signature for Improving Outcomes and Treatment Responses in Patients with Hepatocellular Carcinoma

BACKGROUND

Immune interactions play important roles in the regulation of T cells' cytotoxic function, further impacting the anti-tumor efficacy of immunotherapy. A comprehensive analysis of immune cell types in HCC and immune-cell-related signatures predicting prognosis and monitoring immunotherapy efficacy is still absent.

METHODS

More than 1,300 hepatocellular carcinomas (HCC) patients were collected from public databases and included in the present study. The ssGSEA algorithm was applied to calculate the infiltration level of 28 immunocyte subpopulations. A cell pair algorithm was applied to construct an immune-cell-related prognostic index (ICRPI). Survival analyses were performed to measure the survival difference across ICRPI risk groups. Spearman's correlation analyses were used for the relevance assessment. A Wilcoxon test was used to measure the expression level's differences.

RESULTS

In this study, 28 immune subpopulations were retrieved, and 374 immune cell pairs (ICPs) were established, 38 of which were picked out by the least absolute shrinkage and selection operator (LASSO) algorithm. By using the selected ICPs, the ICRPI was constructed and validated to play crucial roles in survival stratification and dynamic monitoring of immunotherapy effect. We also explored several candidate drugs targeting ICRPI. A composite ICRPI and clinical prognostic index (ICPI) was then constructed, which achieved a more accurate estimation of HCC's survival and is a better choice for prognosis predictions in HCC.

CONCLUSIONS

In conclusion, we constructed and validated ICRPI based on the cell pair algorithm in this study, which might provide some novel insights for increasing the survival estimation and clinical response to immune therapy for individual HCC patients and contribute to the personalized precision immunotherapy strategy of HCC.

Novel exosome-related risk signature as prognostic biomarkers in glioblastoma

Exosomes are progressively being detected as an indicator for the diagnosis and prognosis of cancer in clinical settings. Many clinical trials have confirmed the impact of exosomes on tumor growth, particularly in anti-tumor immunity and immunosuppression of exosomes. Therefore, we developed a risk score based on genes found in glioblastoma-derived exosomes. In this study, we used the TCGA dataset as the training queue and GSE13041, GSE43378, GSE4412, and CGGA datasets as the external validation queue. Based on machine algorithms and bioinformatics methods, an exosome-generalized risk score was established. We found that the risk score could independently predict the prognosis of patients with glioma, and there were significant differences in the outcomes of patients in the high- and low-risk groups. Univariate and multivariate analyses showed that risk score is a valid predictive biomarker for gliomas. Two immunotherapy datasets, IMvigor210 and GSE78220, were obtained from previous studies. A high-risk score showed a significant association with multiple immunomodulators that could act on cancer immune evasion. The exosome-related risk score could predict the effectiveness of anti-PD-1 immunotherapy. Moreover, we compared the sensitivity of patients with high- and low-risk scores to various anti-cancer drugs and found that patients with high-risk scores had better responses to a variety of anti-cancer drugs. The risk-scoring model established in this study provides a useful tool to predict the total survival time of patients with glioma and guide immunotherapy.

Identification of an Immune-Related Gene Signature Associated with Prognosis and Tumor Microenvironment in Esophageal Cancer

Background

Esophageal cancer (EC) is a common malignant tumor of the digestive system with high mortality and morbidity. Current evidence suggests that immune cells and molecules regulate the initiation and progression of EC. Accordingly, it is necessary to identify immune-related genes (IRGs) affecting the biological behaviors and microenvironmental characteristics of EC.

Methods

Bioinformatics methods, including differential expression analysis, Cox regression, and immune infiltration prediction, were conducted using R software to analyze the Gene Expression Omnibus (GEO) dataset. The Cancer Genome Atlas (TCGA) cohort was used to validate the prognostic signature. Patients were stratified into high- and low-risk groups for further analyses, including functional enrichment, immune infiltration, checkpoint relevance, clinicopathological characteristics, and therapeutic sensitivity analyses.

Results

A prognostic signature was established based on 21 IRGs (S100A7, S100A7A, LCN1, CR2, STAT4, GAST, ANGPTL5, TRAV39, F2RL2, PGLYRP3, KLRD1, TRIM36, PDGFA, SLPI, PCSK2, APLN, TICAM1, ITPR3, MAPK9, GATA4, and PLAU). Compared with high-risk patients, better overall survival rates and clinicopathological characteristics were found in low-risk patients. The areas under the curve of the two cohorts were 0.885 and 0.718, respectively. Higher proportions of resting CD4⁺ memory T lymphocytes, M2 macrophages, and resting dendritic cells and lower proportions of follicular helper T lymphocytes, plasma cells, and neutrophils were found in the high-risk tumors. Moreover, the high-risk group showed higher expression of CD44 and TNFSF4, lower expression of PDCD1 and CD40, and higher TIDE scores, suggesting they may respond poorly to immunotherapy. High-risk patients responded better to chemotherapeutic agents such as docetaxel, doxorubicin, and gemcitabine. Furthermore, IRGs associated with tumor progression, including PDGFA, ITPR3, SLPI, TICAM1, and GATA4, were identified.

Conclusion

Our immune-related signature yielded reliable value in evaluating the prognosis, microenvironmental characteristics, and therapeutic sensitivity of EC and may help with the precise treatment of this patient population.

Vemurafenib inhibits immune escape biomarker BCL2A1 by targeting PI3K/AKT signaling pathway to suppress breast cancer

Objectives

To investigate the role of immune escape encoding genes on the prognosis of BC, and to predict the novel targeting agents.

Methods

Human immune genes and immune escape encoding genes were obtained from the IMMPORT database and the previous study. Sample information and clinical data on BC were obtained from the TCGA and GTEX databases. Obtaining differentially expressed protein data from cBioportal database. To construct a risk score model by lasso analysis, and nomogram was used to predict score core. GSCA, TIMER and CELLMINER databases were used for immune and drug susceptibility correlation analyses. Cell experiments were verified by MTT, Western blotting, and RT-qPCR.

Results

We found prognostic models consisting of eleven immune escape related protein-coding genes with ROC curves that performed well in the ontology data (AUC for TCGA is 0.672) and the external data (AUC for GSE20685 is 0.663 and for GES42568 is 0.706). Five core prognostic models are related to survival (EIF4EBP1, BCL2A1, NDRG1, ERRFI1 and BRD4) were summarized, and a nomogram was constructed to validate a C-index of 0.695, which was superior to other prognostic models. Relevant drugs targeting core genes were identified based on drug sensitivity analysis, and found that Vemurafenib downregulates the PI3K-AKT pathway and BCL2A1 protein in BC, as confirmed by external data and cellular assays.

Conclusions

Briefly, our work establishes and validates an 11-immune escape risk model, and five core prognostic factors that are mined deeply from this model, and elucidates in detail that Vemurafenib suppresses breast cancer by targeting the PI3K/AKT signaling pathway to inhibit the immune escape biomarker BCL2A1, confirms the validity of the prognostic model, and provides corresponding targeted agents to guide individualized treatment of BC patients.

EIF4G1 Is a Potential Prognostic Biomarker of Breast Cancer

BACKGROUND

Breast cancer (BRCA) is one of the most common cancers in women worldwide and a leading cause of death from malignancy. This study was designed to identify a novel biomarker for prognosticating the survival of BRCA patients.

METHODS

The prognostic potential of eukaryotic translation initiation factor 4 gamma 1 (EIF4G1) was assessed using RNA sequencing (RNA-seq) data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) as training cohort and validation set, respectively. The functional enrichment analysis of differentially expressed genes (DEGs) was performed. The relationship between EIF4G1 and tumor microenvironment (TME) was analyzed. Immunotherapy responses were explored by the immunophenoscores (IPS) and tumor immune dysfunction and exclusion (TIDE) score. The Connectivity Map (CMap) was used to discover potentially effective therapeutic molecules against BRCA. Immunohistochemistry (IHC) was applied to compare the protein levels of EIF4G1 in normal and cancer tissues and to verify the prognostic value of EIF4G1.

RESULTS

BRCA patients with increased expression of EIF4G1 had a shorter overall survival (OS) in all cohorts and results from IHC. EIF4G1-related genes were mainly involved in DNA replication, BRCA metastasis, and the MAPK signaling pathway. Infiltration levels of CD4⁺-activated memory T cells, macrophages M0, macrophages M1, and neutrophils were higher in the EIF4G1 high-expression group than those in the EIF4G1 low-expression group. EIF4G1 was positively correlated with T cell exhaustion. Lower IPS was revealed in high EIF4G1 expression patients. Five potential groups of drugs against BRCA were identified.

CONCLUSION

EIF4G1 might regulate the TME and affect BRCA metastasis, and it is a potential prognostic biomarker and therapeutic target for BRCA.

A prognostic signature consisting of metabolism-related genes and SLC17A4 serves as a potential biomarker of immunotherapeutic prediction in prostate cancer

Background

Prostate cancer (PCa), a prevalent malignant cancer in males worldwide, screening for patients might benefit more from immuno-/chemo-therapy remained inadequate and challenging due to the heterogeneity of PCa patients. Thus, the study aimed to explore the metabolic (Meta) characteristics and develop a metabolism-based signature to predict the prognosis and immuno-/chemo-therapy response for PCa patients.

Methods

Differentially expressed genes were screened among 2577 metabolism-associated genes. Univariate Cox analysis and random forest algorithms was used for features screening. Multivariate Cox regression analysis was conducted to construct a prognostic Meta-model based on all combinations of metabolism-related features. Then the correlation between MetaScore and tumor was deeply explored from prognostic, genomic variant, functional and immunological perspectives, and chemo-/immuno-therapy response. Multiple algorithms were applied to estimate the immunotherapeutic responses of two MeteScore groups. Further in vitro functional experiments were performed using PCa cells to validate the association between the expression of hub gene SLC17A4 which is one of the model component genes and tumor progression. GDSC database was employed to determine the sensitivity of chemotherapy drugs.

Results

Two metabolism-related clusters presented different features in overall survival (OS). A metabolic model was developed weighted by the estimated regression coefficients in the multivariate Cox regression analysis (0.5154*GAS2 + 0.395*SLC17A4 - 0.1211*NTM + 0.2939*GC). This Meta-scoring system highlights the relationship between the metabolic profiles and genomic alterations, gene pathways, functional annotation, and tumor microenvironment including stromal, immune cells, and immune checkpoint in PCa. Low MetaScore is correlated with increased mutation burden and microsatellite instability, indicating a superior response to immunotherapy. Several medications that might improve patients` prognosis in the MetaScore group were identified. Additionally, our cellular experiments suggested knock-down of SLC17A4 contributes to inhibiting invasion, colony formation, and proliferation in PCa cells in vitro.

Conclusions

Our study supports the metabolism-based four-gene signature as a novel and robust model for predicting prognosis, and chemo-/immuno-therapy response in PCa patients. The potential mechanisms for metabolism-associated genes in PCa oncogenesis and progression were further determined.

Construction of a novel prognostic signature based on the composition of tumor-infiltrating immune cells in clear cell renal cell carcinoma

Emerging evidence has uncovered that tumor-infiltrating immune cells (TIICs) play significant roles in regulating the tumorigenesis and progression of clear cell renal cell carcinoma (ccRCC). However, the exact composition of TIICs and their prognostic values in ccRCC have not been well defined. A total of 534 ccRCC samples with survival information and TIIC data from The Cancer Genome Atlas (TCGA) dataset were included in our research. The ImmuCellAI tool was employed to estimate the abundance of 24 TIICs and further survival analysis explored the prognostic values of TIICs in ccRCC. In addition, the expression levels of immunosuppressive molecules (PDL1, PD1, LAG3, and CTLA4) in the high- and low-risk groups were explored. Various subtypes of TIICs had distinct infiltrating features and most TIICs exhibited dysregulated abundance between normal and tumor tissues. Moreover, specific kinds of TIICs had encouraging prognostic values in ccRCC. Further analysis constructed a 4-TIICs signature to evaluate the prognosis of ccRCC patients. Cox regression analyses confirmed the independent prognostic role of the signature in ccRCC. Moreover, immunosuppressive molecules, including PD1, LAG3, and CTLA4, were significantly upregulated in the high-risk group and predicted poor prognosis. However, PDL1 was not changed between high- and low-risk groups and could not predict poor prognosis. To sum up, our research explored the landscape of TIICs in ccRCC and established a novel 4-TIIC prognostic signature, which could effectively predict the prognosis for patients with ccRCC. Based on this signature, we also concluded that PDL1 may not predict prognosis in ccRCC.

Immune-infiltrating signature-based classification reveals CD103+CD39+ T cells associate with colorectal cancer prognosis and response to immunotherapy

BACKGROUND

Current stratification systems for tumor prognostic prediction and immunotherapeutic efficacy evaluation are less satisfying in colorectal cancer (CRC). As infiltrating immune cells in tumor microenvironment (TME) played a key role in tumor progression and responses to immune checkpoint blockade (ICB) therapy, we want to construct an immune-related scoring system with detailed immune profiles to stratify CRC patients.

METHODS

We developed a scoring system based on immune-related signatures and validated its ability to predict prognosis and immunotherapeutic outcomes in CRC. CD45+ cells from CRC patients were sorted to investigate detailed immune profiles of the stratification system using mass cytometry. A single-cell RNA sequencing dataset was used to analyze transcriptomic profiles.

RESULTS

We constructed an immune-related signature score (IRScore) based on 54 recurrence-free survival (RFS)-related immune signatures to stratify CRC patients. We revealed that IRScore was positively correlated with RFS and favorable outcomes in ICB treatment. Moreover, we depicted a detailed immune profile in TME using mass cytometry and identified that CD103+CD39+ T cells, characterized by an exhaustive, cytotoxic and proliferative phenotype, were enriched in CRC patients with high IRScore. As a beneficial immune signature, CD103+CD39+ T cells could predict prognosis and responses to ICB therapy in CRC.

CONCLUSIONS

All the analyses above revealed that IRScore could be a valuable tool for predicting prognosis and facilitating the development of new therapeutic strategies in CRC, and CD103+CD39+ T cells were one of defined immune signatures in IRScore, which might be a key factor for antitumor immunity.

Cuproptosis Combined with lncRNAs Predicts the Prognosis and Immune Microenvironment of Breast Cancer

Breast cancer (BC), the most common cancer in women, is caused by the uncontrolled proliferation of mammary epithelial cells under the action of a variety of carcinogenic factors. Cuproptosis-related targets have been found to be closely associated with breast cancer development. TCGA obtained 1226 tumor samples, 1073 clinical data, and 37 lncRNAs during univariate Cox multivariate analysis. We used nonnegative matrix factoring (NMF) agglomeration to spot thirty-three potential molecular subsets with totally different cuproptosis-related lncRNA expression patterns. The least absolute shrinkage and selection operator (LASSO) formula and variable Cox multivariate analysis were not used to construct the best prognostic model. The variations in neoplasm mutation burden and factor gene ontology (GO) and gene set enrichment analysis (GSEA) within the high- and low-risk teams were analyzed, and therefore, the potential mechanism of the development of carcinoma was analyzed. We created a prognostic profile consisting of nineteen cuproptosis-related genes (NFE2L2, LIPT1, LIPT2, DLD, etc.) and their connected targets. The correlation between tumor mutational burden (TMB) and clinical manifestations of tumors demonstrates the importance of high- and low-expression bunch data on the incidence of clinical manifestations of tumors. The area under the curve (AUC) shows moderate prophetic power for copper mortality. GO enrichment analysis showed that immunorelated responses were enriched. Correlation analysis of immune cells showed that pathology could play an important role in the prevalence and prognosis of tumors, and there were variations in immune cells between the probable and low-risk groups. Our study suggests that the prognostic characteristic genes associated with cuproptosis can be used as new biomarkers to predict the prognosis of breast cancer patients. In addition, we found that immunotherapy may play a key role in breast cancer treatment regimens. Levels of immune-associated cells and pathways vary significantly among risk groups of breast cancer patients.

Comprehensive analysis of PTPN gene family revealing PTPN7 as a novel biomarker for immuno-hot tumors in breast cancer

Background: The non-receptor protein tyrosine phosphatase (PTPN) gene family has been considered to be involved in the oncogenesis and development of multiple cancers. However, its prognostic utility and immunological relevance in breast cancer (BrCa) have not been clarified.

Methods: A transcriptional level interpretation of the expressions and prognostic values was analyzed using the data from The Cancer Genome Atlas (TCGA) cohort. In addition, GO and DAVID pinpoint the functional enrichment of PTPNs. Moreover, the immune correlations of PTPN7 in BrCa and pan-cancer were further investigated based on the TCGA cohort and were testified using the in-house and the Gene Expression Omnibus (GEO) cohorts.

Results: For systematic analysis of the PTPN family, we found that the expression levels of PTPN1, PTPN6, PTPN7, PTPN18, PTPN20, and PTPN22 was promoted in tumor tissues while comparing with paraneoplastic tissues during our study. We further investigated their functions and protein-protein interactions (PPI), and these results strongly suggested that PTPN family was associated with protein dephosphorylation. Next, we performed an immunological relevance analysis and found that PTPN7 was correlated with immune infiltration, suggesting a stronger association of PTPN7 with immuno-hot tumors in BrCa. In addition, results from the in-house cohort confirmed the positive correlation between PTPN7 and PD-L1. The pan-cancer analysis revealed that PTPN7 was related to PD-L1 and CTLA-4 expression in almost all cancer types. Finally, the predictive value of PTPN7 for immunotherapy was significant in two independent GEO cohorts.

Conclusion: In conclusion, this is the first extensive research on the correlation between PTPN family expression and immune characterization in BrCa. As results, PTPN7 expression is associated with immuno-hot tumors and could be a promising predictive biomarker for immunotherapy in not only BrCa but multiple cancers.

Systematic pan-cancer analysis identifies APOC1 as an immunological biomarker which regulates macrophage polarization and promotes tumor metastasis

Apolipoprotein C1 (APOC1) has been found to play an essential part in proliferation and metastasis of numerous cancers, but related mechanism has not been elucidated, especially its function and role in tumor immunity. Through systematic pan-cancer analysis, we identified that APOC1 was closely associated with the infiltration of various immune cells in multiple cancers. Besides, APOC1 was significantly co-expressed with the immune checkpoints, major histocompatibility complex (MHC) molecules, chemokines and other immune-related genes. Furthermore, single-cell sequencing analysis suggested that the vast majority of APOC1 was expressed in macrophages or tumor-associated macrophages (TAMs). Additionally, the expression of APOC1 was significantly related to the prognosis of different cancers. Since APOC1 was most significantly abnormally expressed in renal cell cancer (RCC), subsequent experiments were carried out in RCC to explore the role of APOC1 in tumor immunity. The expression of APOC1 was significantly elevated in the tumor and serum of RCC patients. Besides, APOC1 was mainly expressed in the macrophage and it was closely related to the immune cell infiltration of RCC. Co-culture with RCC cells could induce the generation of TAMs with M2 phenotype which be blocked by silencing APOC1. The expression of APOC1 was elevated in the M2 or TAMs and APOC1 promoted M2 polarization of macrophages through interacting with CD163 and CD206. Furthermore, macrophages overexpressing APOC1 promoted the metastasis of RCC cells via secreting CCL5. Together, these data indicate that APOC1 is an immunological biomarker which regulates macrophage polarization and promotes tumor metastasis.

Pan-Cancer Analysis Reveals the Relation between TRMT112 and Tumor Microenvironment

Dysregulated epigenetic modifications play a critical role in cancer development where TRMT112 is a member of the transfer RNA (tRNA) methyltransferase family. Till now, no studies have revealed the linkage between TRMT112 expression and diverse types of tumors. Based on TCGA data, we first probed into the relation between TRMT112 and prognosis and the potential role of TRMT112 in tumor microenvironment across 33 types of tumor. TRMT112 presented with increased expression in most cancers, which was significantly prognostic. Furthermore, TRMT112 was associated with tumor-associated fibroblasts in a variety of cancers. Additionally, a positive relationship was identified between TRMT112 expression and multiple tumor-related immune infiltrations, such as dendritic cells, CD8+ T cells, macrophages, CD4+ T cells, neutrophils, and B cells in lung adenocarcinoma and breast invasive carcinoma. In summary, our results suggest that TRMT112 might be a potential prognostic predictor of cancers and involved in regulating multiple cancer-related immune responses to some extent.

Spindle and Kinetochore-associated Family Genes are Prognostic and Predictive Biomarkers in Hepatocellular Carcinoma

BACKGROUND AND AIMS

Hepatocellular carcinoma (HCC) is one of the most frequent malignant tumors. Spindle and kinetochore-associated (SKA) family genes are essential for the maintenance of the metaphase plate and spindle checkpoint silencing during mitosis. Recent studies have indicated that dysregulation of SKA family genes induces tumorigenesis, tumor progression, and chemoresistance via modulation of cell cycle and DNA replication. However, the differential transcription of SKAs in the context of HCC and its prognostic significance has not been demonstrated.

METHODS

Bioinformatics analyses were performed using TCGA, ONCOMINE, HCCDB, Kaplan-Meier plotter, STRING, GEPIA databases. qRT-PCR, western blot, and functional assays were utilized for in vitro experiments.

RESULTS

We found remarkable upregulation of transcripts of SKA family genes in HCC samples compared with normal liver samples on bioinformatics analyses and in vitro validation. Interaction analysis and enrichment analysis showed that SKA family members were mainly related to microtubule motor activity, mitosis, and cell cycle. Immuno-infiltration analysis showed a correlation of all SKA family genes with various immune cell subsets, especially T helper 2 (Th2) cells. Transcriptional levels of SKA family members were positively associated with histologic grade, T stage, and α-fetoprotein in HCC patients. Receiver operating characteristic curve analysis demonstrated a strong predictive ability of SKA1/2/3 for HCC. Increased expression of these SKAs was associated with unfavorable overall survival, progression-free survival, and disease-specific survival. On Cox proportional hazards regression analyses, SKA1 upregulation and pathological staging were independent predictors of overall survival and disease-specific survival of HCC patients. Finally, clinical tissue microarray validation and in vitro functional assays revealed SKA1 acts an important regulatory role in tumor malignant behavior.

CONCLUSIONS

SKA family members may potentially serve as diagnostic and prognostic markers in the context of HCC. The correlation between SKAs and immune cell infiltration provides a promising research direction for SKA-targeted immunotherapeutics for HCC.

The construction of transcriptional risk scores for breast cancer based on lightGBM and multiple omics data

BACKGROUND

Polygenic risk score (PRS) can evaluate the individual-level genetic risk of breast cancer. However, standalone single nucleotide polymorphisms (SNP) data used for PRS may not provide satisfactory prediction accuracy. Additionally, current PRS models based on linear regression have insufficient power to leverage non-linear effects from thousands of associated SNPs. Here, we proposed a transcriptional risk score (TRS) based on multiple omics data to estimate the risk of breast cancer.

METHODS

The multiple omics data and clinical data of breast invasive carcinoma (BRCA) were collected from the cancer genome atlas (TCGA) and the gene expression omnibus (GEO). First, we developed a novel TRS model for BRCA utilizing single omic data and LightGBM algorithm. Subsequently, we built a combination model of TRS derived from each omic data to further improve the prediction accuracy. Finally, we performed association analysis and prognosis prediction to evaluate the utility of the TRS generated by our method.

RESULTS

The proposed TRS model achieved better predictive performance than the linear models and other ML methods in single omic dataset. An independent validation dataset also verified the effectiveness of our model. Moreover, the combination of the TRS can efficiently strengthen prediction accuracy. The analysis of prevalence and the associations of the TRS with phenotypes including case-control and cancer stage indicated that the risk of breast cancer increases with the increases of TRS. The survival analysis also suggested that TRS for the cancer stage is an effective prognostic metric of breast cancer patients.

CONCLUSIONS

Our proposed TRS model expanded the current definition of PRS from standalone SNP data to multiple omics data and outperformed the linear models, which may provide a powerful tool for diagnostic and prognostic prediction of breast cancer.

Identification of Potential Biomarkers of Platelet RNA in Glioblastoma by Bioinformatics Analysis

Objective

Glioblastoma is one of the most common and fatal malignancies in adults. Current treatment is still not optimistic. Glioblastoma (GBM) transports RNA to platelets in the blood system via microvesicles, suggesting that platelet RNA can be a potential diagnostic and therapeutic target. The roles of specific platelet RNAs in treatment of GBM are not well understood.

Methods

Platelet RNA profiling of 8 GBM and 12 normal samples were downloaded from the GEO database. Differentially expressed genes (DEGs) were identified between tumors and normal samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to elucidate the functions of up- and downregulated genes. miRNA was predicted by miRTarBase, TargetScan, and miRDB databases. circBase and circBank were used for circRNA prediction. ceRNA (circRNA-mRNA-miRNA) network was constructed to investigate the potential interactions.

Results

22 genes were upregulated and 9 genes were downregulated. There are only two genes (CCR7 and FAM102A) that connect to miRNAs (hsa-let-7a-5p, hsa-miR-1-3p). We assessed the overall survival rates by Kaplan-Meier plotter, and relative expression of GBM and subtypes for overlapped mRNA (CCR7 and FAM102A) were evaluated, and further, we obtained circRNAs (has-circ-0015164, hsa-circ-0003243) by circBank and circBase and bind sites through the CSCD database. Finally, a ceRNA network (circRNA-mRNA-miRNA) was constructed based on 2 miRNAs, 2 mRNAs, and 2 circRNAs by Cytoscape. This study focused on potential mRNA and ceRNA biomarkers to targeted treatment of GBM and provided ideas for clinical treatment through the combination of hematology and oncology.

Conclusion

The findings of this study contribute to better understand the relationship between GBM and the blood system (platelets) and might lay a solid foundation for improving GBM molecule and gene diagnosis and prognosis.

Crosstalk of three novel types of programmed cell death defines distinct microenvironment characterization and pharmacogenomic landscape in breast cancer

Background

Prior studies have highlighted that novel programmed cell death (PCD) modalities, including ferroptosis, pyroptosis, and necroptosis, are correlated with tumor progression and antitumor immunity. Nonetheless, comprehensive analysis of tumor microenvironment (TME) profiles mediated by the crosstalk of distinct PCD forms has not been conducted in breast cancer (BC).

Methods

Here, we curated 34 identified PCD-associated genes (PCDAGs) and applied the consensus clustering algorithm to establish PCD-mediated tumor patterns in BC. Subsequently, based on prognostic differentially expressed genes extracted from distinct PCD-mediated patterns, we applied the LASSO algorithm to construct CD_Score. Furthermore, the correlation analysis between CD_Score and TME features, molecular subtypes, clinicopathological characteristics, drug response, and immunotherapeutic efficacy was performed.

Results

Three distinct PCD-clusters were determined among 2,038 BC samples, which did not only display different clinical outcomes but highly correlated to the established immunological tumor phenotypes: “desert,” “excluded,” and “inflamed” immune profiles. Based on the CD_Score derived from the PCD-related gene signature, BC patients could be stratified into CD_Score-low and -high group, of which the former displayed satisfactory survival outcome and enhanced immune infiltration. Further exploration identified that the CD_Score-high group significantly correlated with elevated neoantigen load and higher mutation frequency in SMGs (e.g., TP53 and MAP3K1) and reduced expression of immune checkpoint proteins.

Conclusions

This research is the first to emphasize the close relationship between distinct cell death modalities and the diversity and complexity of immune infiltration in TME. We established the CD_Score, which could help enhance our cognition of TME features and facilitate the clinical application of immunotherapy.

Hallmark-guided subtypes of hepatocellular carcinoma for the identification of immune-related gene classifiers in the prediction of prognosis, treatment efficacy, and drug candidates

Hepatocellular carcinoma (HCC), accounting for ~90% of all primary liver cancer, is a prevalent malignancy worldwide. The intratumor heterogeneity of its causative etiology, histology, molecular landscape, and immune phenotype makes it difficult to precisely recognize individuals with high mortality risk or tumor-intrinsic treatment resistance, especially immunotherapy. Herein, we comprehensively evaluated the activities of cancer hallmark gene sets and their correlations with the prognosis of HCC patients using gene set variation analysis (GSVA) and identified two HCC subtypes with distinct prognostic outcomes. Based on these subtypes, seven immune-related genes (TMPRSS6, SPP1, S100A9, EPO, BIRC5, PLXNA1, and CDK4) were used to construct a novel prognostic gene signature [hallmark-guided subtypes-based immunologic signature (HGSIS)] via multiple statistical approaches. The HGSIS-integrated nomogram suggested an enhanced predictive performance. Interestingly, oncogenic hallmark pathways were significantly enriched in the high-risk group and positively associated with the risk score. Distinct mutational landscapes and immune profiles were observed between different risk groups. Moreover, immunophenoscore (IPS) and tumor immune dysfunction and exclusion (TIDE) analysis showed different sensitivities of HGSIS risk groups for immune therapy efficacy, and the pRRophetic algorithm indicated distinguishable responses for targeted/chemotherapies in different groups. KIF2C was picked out as the key target concerning HGSIS, and the top 10 small molecules were predicted to bind to the active site of KIF2C via molecular docking, which might be further used for candidate drug discovery of HCC. Taken together, our study offers novel insights for clinically significant subtype recognition, and the proposed signature may be a helpful guide for clinicians to improve the treatment regimens.

Oncogenic signaling pathway-related long non-coding RNAs for predicting prognosis and immunotherapy response in breast cancer

Background

The clinical outcomes of breast cancer (BC) are unpredictable due to the high level of heterogeneity and complex immune status of the tumor microenvironment (TME). When set up, multiple long non-coding RNA (lncRNA) signatures tended to be employed to appraise the prognosis of BC. Nevertheless, predicting immunotherapy responses in BC is still essential. LncRNAs play pivotal roles in cancer development through diverse oncogenic signal pathways. Hence, we attempted to construct an oncogenic signal pathway–based lncRNA signature for forecasting prognosis and immunotherapy response by providing reliable signatures.

Methods

We preliminarily retrieved RNA sequencing (RNA-seq) data from The Cancer Genome Atlas (TCGA) database and extracted lncRNA profiles by matching them with GENCODE. Following this, Gene Set Variation Analysis (GSVA) was used to identify the lncRNAs closely associated with 10 oncogenic signaling pathways from the TCGA-BRCA (breast-invasive carcinoma) cohort and was further screened by the least absolute shrinkage and selection operator Cox regression model. Next, an lncRNA signature (OncoSig) was established through the expression level of the final 29 selected lncRNAs. To examine survival differences in the stratification described by the OncoSig, the Kaplan–Meier (KM) survival curve with the log-rank test was operated on four independent cohorts (n = 936). Subsequently, multiple Cox regression was used to investigate the independence of the OncoSig as a prognostic factor. With the concordance index (C-index), the time-dependent receiver operating characteristic was employed to assess the performance of the OncoSig compared to other publicly available lncRNA signatures for BC. In addition, biological differences between the high- and low-risk groups, as portrayed by the OncoSig, were analyzed on the basis of statistical tests. Immune cell infiltration was investigated using gene set enrichment analysis (GSEA) and deconvolution tools (including CIBERSORT and ESTIMATE). The combined effect of the Oncosig and immune checkpoint genes on prognosis and immunotherapy was elucidated through the KM survival curve. Ultimately, a pan-cancer analysis was conducted to attest to the prevalence of the OncoSig.

Results

The OncoSig score stratified BC patients into high- and low-risk groups, where the latter manifested a significantly higher survival rate and immune cell infiltration when compared to the former. A multivariate analysis suggested that OncoSig is an independent prognosis predictor for BC patients. In addition, compared to the other four publicly available lncRNA signatures, OncoSig exhibited superior predictive performance (AUC = 0.787, mean C-index = 0.714). The analyses of the OncoSig and immune checkpoint genes clarified that a lower OncoSig score meant significantly longer survival and improved response to immunotherapy. In addition to BC, a high OncoSig score in several other cancers was negatively correlated with survival and immune cell infiltration.

Conclusions

Our study established a trustworthy and discriminable prognostic signature for BC patients with similar clinical profiles, thus providing a new perspective in the evaluation of immunotherapy responses. More importantly, this finding can be generalized to be applicable to the vast majority of human cancers.