Development of a membrane lipid metabolism-based signature to predict overall survival for personalized medicine in ccRCC patients

Integration of bioinformatics and machine learning strategies identifies APM-related gene signatures to predict clinical outcomes and therapeutic responses for breast cancer patients

BACKGROUND

Tumor antigenicity and efficiency of antigen presentation jointly influence tumor immunogenicity, which largely determines the effectiveness of immune checkpoint blockade (ICB). However, the role of altered antigen processing and presentation machinery (APM) in breast cancer (BRCA) has not been fully elucidated.

METHODS

A series of bioinformatic analyses and machine learning strategies were performed to construct APM-related gene signatures to guide personalized treatment for BRCA patients. A single-sample gene set enrichment analysis (ssGSEA) algorithm and weighted gene co-expression network analysis (WGCNA) were combined to screen for BRCA-specific APM-related genes. The non-negative matrix factorization (NMF) algorithm was used to divide the cohort into different clusters and the fgsea algorithm was applied to investigate the altered signaling pathways. Random survival forest (RSF) and the least absolute shrinkage and selection operator (Lasso) Cox regression analysis were combined to construct an APM-related risk score (APMrs) signature to predict overall survival. Furthermore, a nomogram and decision tree were generated to improve predictive accuracy and risk stratification for individual patients. Based on Tumor Immune Dysfunction and Exclusion (TIDE) method, random forest (RF) and Lasso logistic regression model were combined to establish an APM-related immunotherapeutic response score (APMis). Finally, immune infiltration, immunomodulators, mutational patterns, and potentially applicable drugs were comprehensively analyzed in different APM-related risk groups. IHC staining was used to assess the expression of APM-related genes in clinical samples.

RESULTS

In this study, APMrs and APMis showed favorable performances in risk stratification and therapeutic prediction for BRCA patients. APMrs exhibited more powerful prognostic capacity and accurate survival prediction compared to conventional clinicopathological features. APMrs was closely associated with distinct mutational patterns, immune cell infiltration and immunomodulators expression. Furthermore, the two APM-related gene signatures were independently validated in external cohorts with prognosis or immunotherapeutic responses. Potential applicable drugs and targets were mined in the APMrs-high group. APM-related genes were further validated in our in-house samples.

CONCLUSION

The APM-related gene signatures established in our study could improve the personalized assessment of survival risk and guide ICB decision-making for BRCA patients.

Identification of two novel lipid metabolism-related long non-coding RNAs (SNHG17 and LINC00837) as potential signatures for osteosarcoma prognosis and precise treatment

OBJECTIVE

Dysregulated lipid metabolism enhances the development and advancement of many cancers, including osteosarcoma (OS); however, the underlying mechanisms are still largely unknown. Therefore, this investigation aimed to elucidate novel potential lipid metabolism-related long non-coding RNAs (lncRNAs) that regulate OS development and provide novel signatures for its prognosis and precise treatment.

MATERIALS AND METHODS

The GEO datasets (GSE12865 and GSE16091) were downloaded and analyzed using R software packages. Immunohistochemistry (IHC) was used to evaluate protein levels in OS tissues while real-time qPCR was used to measure lncRNA levels, and MTT assays were used to assess OS cell viability.

RESULTS

Two lipid metabolism-associated lncRNAs (LM-lncRNAs), small nucleolar RNA host gene 17 (SNHG17) and LINC00837, were identified as efficient and independent prognostic indicators for OS. In addition, further experiments confirmed that SNHG17 and LINC00837 were significantly elevated in OS tissues and cells than para-cancerous counterparts. Knockdown of SNHG17 and LINC00837 synergistically suppressed the viability of OS cells, whereas overexpression of the two lncRNAs promoted OS cell proliferation. Moreover, bioinformatics analysis was conducted to construct six novel SNHG17-microRNA-mRNA competing endogenous RNA (ceRNA) networks, and three lipid metabolism-associated genes (MIF, VDAC2, and CSNK2A2) were found to be abnormally upregulated in OS tissues, suggesting that they were potential effector genes of SNHG17.

CONCLUSION

In summary, SNHG17 and LINC00837 were found to promote OS cell malignancy, suggesting their use as ideal biomarkers for OS prognosis and treatment.

A novel fatty-acid metabolism-based classification for triple negative breast cancer

BACKGROUND

The high heterogeneity of triple negative breast cancer (TNBC) is the main clinical challenge for individualized therapy. Considering that fatty acid metabolism (FAM) plays an indispensable role in tumorigenesis and development of TNBC, we proposed a novel FAM-based classification to characterize the tumor microenvironment immune profiles and heterogeneous for TNBC.

METHODS

Weighted gene correlation network analysis (WGCNA) was performed to identify FAM-related genes from 221 TNBC samples in Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) dataset. Then, non-negative matrix factorization (NMF) clustering analysis was applied to determine FAM clusters based on the prognostic FAM-related genes, which chosen from the univariate/multivariate Cox regression model and the least absolute shrinkage and selection operator (LASSO) regression algorithm. Then, a FAM scoring scheme was constructed to further quantify FAM features of individual TNBC patient based on the prognostic differentially expressed genes (DEGs) between different FAM clusters. Systematically analyses were performed to evaluate the correlation between the FAM scoring system (FS) with survival outcomes, genomic characteristics, tumor microenvironment (TME) features and immunotherapeutic response for TNBC, which were further validated in the Cancer Genome Atlas (TCGA) and GSE58812 datasets. Moreover, the expression level and clinical significancy of the selected FS gene signatures were further validated in our cohort.

RESULTS

1860 FAM-genes were screened out using WGCNA. Three distinct FAM clusters were determined by NMF clustering analysis, which allowed to distinguish different groups of patients with distinct clinical outcomes and tumor microenvironment (TME) features. Then, prognostic gene signatures based on the DEGs between different FAM clusters were identified using univariate Cox regression analysis and Lasso regression algorithm. A FAM scoring scheme was constructed, which could divide TNBC patients into high and low-FS subgroups. Low FS subgroup, characterized by better prognosis and abundance with effective immune infiltration. While patients with higher FS were featured with poorer survival and lack of effective immune infiltration. In addition, two independent immunotherapy cohorts (Imvigor210 and GSE78220) confirmed that patients with lower FS demonstrated significant therapeutic advantages from anti-PD-1/PD-L1 immunotherapy and durable clinical benefits. Further analyses in our cohort found that the differential expression of CXCL13, FBP1 and PLCL2 were significantly associated with clinical outcomes of TNBC samples.

CONCLUSIONS

This study revealed FAM plays an indispensable role in formation of TNBC heterogeneity and TME diversity. The novel FAM-based classification could provide a promising prognostic predictor and guide more effective immunotherapy strategies for TNBC.

Bioinformatics Analysis and Structure of Gastric Cancer Prognosis Model Based on Lipid Metabolism and Immune Microenvironment

OBJECTIVES

The reprogramming of lipid metabolism is a new trait of cancers. However, the role of lipid metabolism in the tumor immune microenvironment (TIME) and the prognosis of gastric cancer remains unclear.

METHODS

Consensus clustering was applied to identify novel subgroups. ESTIMATE, TIMER, and MCPcounter algorithms were used to determine the TIME of the subgroups. The underlying mechanisms were elucidated using functional analysis. The prognostic model was established using the LASSO algorithm and multivariate Cox regression analysis.

RESULTS

Three molecular subgroups with significantly different survival were identified. The subgroup with relatively low lipid metabolic expression had a lower immune score and immune cells. The differentially expressed genes (DEGs) were concentrated in immune biological processes and cell migration via GO and KEGG analyses. GSEA analysis showed that the subgroups were mainly enriched in arachidonic acid metabolism. Gastric cancer survival can be predicted using risk models based on lipid metabolism genes.

CONCLUSIONS

The TIME of gastric cancer patients is related to the expression of lipid metabolism genes and could be used to predict cancer prognosis accurately.

A novel mTOR-associated gene signature for predicting prognosis and evaluating tumor immune microenvironment in lung adenocarcinoma

BACKGROUND

The mechanistic target of rapamycin (mTOR) was proven to have great impact on apoptosis, cell proliferation, autophagy, and many other fundamental cellular processes; moreover, it closely correlates with tumor occurrence and development. However, few studies have constructed signatures based on mTOR-associated genes to assess multiple indicators of prognosis in lung adenocarcinoma (LUAD) patients.

METHODS

mTOR-associated gene sets, whole mRNA expression matrices, and clinical information of LUAD patients in training and validation cohorts were obtained from multiple public databases. Multiple methods were used to screen candidate genes, construct signatures, validate internally and externally, and conduct further studies: differentially expressed gene analysis, LASSO Cox regression analysis, Cox regression analysis, risk factor analysis, nomogram analysis, functional enrichment analysis, analyses in tumor immune microenvironment, and therapy.

RESULTS

A prognostic signature containing 8 genes (LDHA, SLA, WNT7A, PLK1, CCT6A, BTG2, TXNRD1, and DDIT4) was constructed. It performed well in both internal and external validation. Subsequent analysis found that the prognostic signature was of great significance in evaluating the tumor immune microenvironment and could guide the treatment of patients with LUAD to a certain extent.

CONCLUSION

The constructed mTOR-associated gene signature accurately predicted the prognostic pattern of patients with LUAD and is expected to be extremely useful in guiding LUAD therapy.

Expression of Lipid-Metabolism Genes Is Correlated With Immune Microenvironment and Predicts Prognosis in Osteosarcoma

Objectives

Osteosarcoma was the most popular primary malignant tumor in children and adolescent, and the 5-year survival of osteosarcoma patients gained no substantial improvement over the past 35 years. This study aims to explore the role of lipid metabolism in the development and diagnosis of osteosarcoma.

Methods

Clinical information and corresponding RNA data of osteosarcoma patients were downloaded from TRGET and GEO databases. Consensus clustering was performed to identify new molecular subgroups. ESTIMATE, TIMER and ssGSEA analyses were applied to determinate the tumor immune microenvironment (TIME) and immune status of the identified subgroups. Functional analyses including GO, KEGG, GSVA and GSEA analyses were conducted to elucidate the underlying mechanisms. Prognostic risk model was constructed using LASSO algorithm and multivariate Cox regression analysis.

Results

Two molecular subgroups with significantly different survival were identified. Better prognosis was associated with high immune score, low tumor purity, high abundance of immune infiltrating cells and relatively high immune status. GO and KEGG analyses revealed that the DEGs between the two subgroups were mainly enriched in immune- and bone remodeling-associated pathways. GSVA and GSEA analyses indicated that, lipid catabolism downregulation and lipid hydroxylation upregulation may impede the bone remodeling and development of immune system. Risk model based on lipid metabolism related genes (LMRGs) showed potent potential for survival prediction in osteosarcoma. Nomogram integrating risk model and clinical characteristics could predict the prognosis of osteosarcoma patients accurately.

Conclusion

Expression of lipid-metabolism genes is correlated with immune microenvironment of osteosarcoma patients and could be applied to predict the prognosis of in osteosarcoma accurately.

Utility of a metabolic-associated nomogram to predict the recurrence-free survival of stage I cervical cancer

Aims: To identify metabolism-associated genes (MAGs) that serve as biomarkers to predict prognosis associated with recurrence-free survival (RFS) for stage I cervical cancer (CC). Patients & methods: By analyzing the Gene Expression Omnibus (GEO) database for 258 cases of stage I CC via univariate Cox analysis, LASSO and multivariate Cox regression analysis, we unveiled 11 MAGs as a signature that was also validated using Kaplan-Meier and receiver operating characteristic analyses. In addition, a metabolism-related nomogram was developed. Results: High accuracy of this signature for prediction was observed (area under the curve at 1, 3 and 5 years was 0.964, 0.929 and 0.852 for the internal dataset and 0.759, 0.719 and 0.757 for the external dataset). The high-risk score group displayed markedly worse RFS than did the low-risk score group. The indicators performed well in our nomogram. Conclusions: We identified a novel signature as a biomarker for predicting prognosis and a nomogram to facilitate the individual management of stage I CC patients.

Development and Validation of Epigenetic Signature Predict Survival for Patients with Laryngeal Squamous Cell Carcinoma

Establishing epigenetic signature to improve the accuracy of survival prediction and optimize therapeutic strategies for laryngeal squamous cell carcinoma (LSCC) by a genome-wide integrated analysis of methylation and the transcriptome. LSCC DNA methylation datasets and RNA sequencing datasets were acquired from the Cancer Genome Atlas (TCGA). MethylMix was applied to detect DNA methylation-driven genes (MDGs), which developed an epigenetic signature. The predictive accuracy and clinical value of the epigenetic signature were evaluated by receiver operating characteristic and decision curve analysis, and compared with tumor-node-metastasis (TNM) stage system. In addition, prognostic value of the epigenetic signature was validated by external Gene Expression Omnibus (GEO) database. According to five MDGs of epigenetic signature, the candidate small molecules for LSCC were screen out by the CMap database. A total of 88 DNA MDGs were identified, five of which (MAGEB2, SUSD1, ZNF382, ZNF418, and ZNF732) were chosen to construct an epigenetic signature. The epigenetic signature can effectively divide patients into high-risk and low-risk group, with the area under curve (AUC) of 0.8 (5-year overall survival [OS]) and AUC of 0.745 (3-year OS). Stratification analysis affirmed that the epigenetic signature was still a significant statistical prognostic model in subsets of patients with different clinical variables. Multivariate Cox regression analysis indicated that the efficacy of epigenetic signature appears independent of other clinicopathological characteristics. In terms of predictive capacity and clinical usefulness, the epigenetic signature was superior to traditional TNM stage. In addition, the epigenetic signature was confirmed in external LSCC cohorts from GEO. Finally, CMap matched the 10 most significant small molecules as promising therapeutic drugs to reverse the LSCC gene expression. An epigenetic signature, with five DNA MDGs, was identified and validated in LSCC patients by integrating multidimensional genomic data, which may offer novel research directions and prospects for individualized treatment of patients with LSCC.

Identification and validation of methylation-driven genes prognostic signature for recurrence of laryngeal squamous cell carcinoma by integrated bioinformatics analysis

Background

Recurrence remains a major obstacle to long-term survival of laryngeal squamous cell carcinoma (LSCC). We conducted a genome-wide integrated analysis of methylation and the transcriptome to establish methylation-driven genes prognostic signature (MDGPS) to precisely predict recurrence probability and optimize therapeutic strategies for LSCC.

Methods

LSCC DNA methylation datasets and RNA sequencing (RNA-seq) dataset were acquired from the Cancer Genome Atlas (TCGA). MethylMix was applied to detect DNA methylation-driven genes (MDGs). By univariate and multivariate Cox regression analyses, five genes of DNA MDGs was developed a recurrence-free survival (RFS)-related MDGPS. The predictive accuracy and clinical value of the MDGPS were evaluated by receiver operating characteristic (ROC) and decision curve analysis (DCA), and compared with TNM stage system. Additionally, prognostic value of MDGPS was validated by external Gene Expression Omnibus (GEO) database. According to 5 MDGs, the candidate small molecules for LSCC were screen out by the CMap database. To strengthen the bioinformatics analysis results, 30 pairs of clinical samples were evaluated by digoxigenin-labeled chromogenic in situ hybridization (CISH).

Results

A total of 88 DNA MDGs were identified, and five RFS-related MDGs (LINC01354, CCDC8, PHYHD1, MAGEB2 and ZNF732) were chosen to construct a MDGPS. The MDGPS can effectively divide patients into high-risk and low-risk group, with the area under curve (AUC) of 0.738 (5-year RFS) and AUC of 0.74 (3-year RFS). Stratification analysis affirmed that the MDGPS was still a significant statistical prognostic model in subsets of patients with different clinical variables. Multivariate Cox regression analysis indicated the efficacy of MDGPS appears independent of other clinicopathological characteristics. In terms of predictive capacity and clinical usefulness, the MDGPS was superior to traditional TNM stage. Additionally, the MDGPS was confirmed in external LSCC cohorts from GEO. CMap matched the 9 most significant small molecules as promising therapeutic drugs to reverse the LSCC gene expression. Finally, CISH analysis in 30 LSCC tissues and paired adjacent normal tissues revealed that MAGEB2 has significantly higher expression of LSCC compared to adjacent non-neoplastic tissues; LINC01354, CCDC8, PHYHD1, and ZNF732 have significantly lower expression of LSCC compared to adjacent non-neoplastic tissues, which were in line with bioinformatics analysis results.

Conclusion

A MDGPS, with five DNA MDGs, was identified and validated in LSCC patients by combining transcriptome and methylation datasets analysis. Compared TNM stage alone, it generates more accurate estimations of the recurrence prediction and maybe offer novel research directions and prospects for individualized treatment of patients with LSCC.

Mast cell-based molecular subtypes and signature associated with clinical outcome in early-stage lung adenocarcinoma

Mast cells are a major component of the immune microenvironment in tumour tissues and modulate tumour progression by releasing pro‐tumorigenic and antitumorigenic molecules. Regarding the impact of mast cells on the outcomes of patients with lung adenocarcinoma (LUAD) patient, several published studies have shown contradictory results. Here, we aimed at elucidating the role of mast cells in early‐stage LUAD. We found that high mast cell abundance was correlated with prolonged survival in early‐stage LUAD patients. The mast cell‐related gene signature and gene mutation data sets were used to stratify early‐stage LUAD patients into two molecular subtypes (subtype 1 and subtype 2). The neural network‐based framework constructed with the mast cell‐related signature showed high accuracy in predicting response to immunotherapy. Importantly, the prognostic mast cell‐related signature predicted the survival probability and the potential relationship between TP53 mutation, c‐MYC activation and mast cell activities. The meta‐analysis confirmed the prognostic value of the mast cell‐related gene signature. In summary, this study might improve our understanding of the role of mast cells in early‐stage LUAD and aid in the development of immunotherapy and personalized treatments for early‐stage LUAD patients.

A novel epigenetic signature for overall survival prediction in patients with breast cancer

BACKGROUND

Breast cancer is the most common malignancy in female patients worldwide. Because of its heterogeneity in terms of prognosis and therapeutic response, biomarkers with the potential to predict survival or assist in making treatment decisions in breast cancer patients are essential for an individualised therapy. Epigenetic alterations in the genome of the cancer cells, such as changes in DNA methylation pattern, could be a novel marker with an important role in the initiation and progression of breast cancer.

METHOD

DNA methylation and RNA-seq datasets from The Cancer Genome Atlas (TCGA) were analysed using the Least Absolute Shrinkage and Selection Operator (LASSO) Cox model. Applying gene ontology (GO) and single sample gene set enrichment analysis (ssGSEA) an epigenetic signature associated with the survival of breast cancer patients was constructed that yields the best discrimination between tumour and normal breast tissue. A predictive nomogram was built for the optimal strategy to distinguish between high- and low-risk cases.

RESULTS

The combination of mRNA-expression and of DNA methylation datasets yielded a 13-gene epigenetic signature that identified subset of breast cancer patients with low overall survival. This high-risk group of tumor cases was marked by upregulation of known cancer-related pathways (e.g. mTOR signalling). Subgroup analysis indicated that this epigenetic signature could distinguish high and low-risk patients also in different molecular or histological tumour subtypes (by Her2-, EGFR- or ER expression or different tumour grades). Using Gene Expression Omnibus (GEO) the 13-gene signature was confirmed in four external breast cancer cohorts.

CONCLUSION

An epigenetic signature was discovered that effectively stratifies breast cancer patients into low and high-risk groups. Since its efficiency appears independent of other known classifiers (such as staging, histology, metastasis status, receptor status), it has a high potential to further improve likely individualised therapy in breast cancer.