BRCA1-Associated Protein Is a Potential Prognostic Biomarker and Is Correlated With Immune Infiltration in Liver Hepatocellular Carcinoma: A Pan-Cancer Analysis

RRP8, associated with immune infiltration, is a prospective therapeutic target in hepatocellular carcinoma

BACKGROUND

Ribosomal RNA Processing 8 (RRP8) is a nucleolar Rossman fold-like methyltransferase that exhibits increased expression in many malignant tumours. However, the role of RRP8 in hepatocellular carcinoma (HCC) is still uncertain. We explored the relationships between RRP8 and prognosis and immune infiltration, as well as the putative pathological function and mechanism of RRP8 in HCC.

METHODS

Analysis of RRP8 expression across cancers was performed by using multiple databases. Associations between RRP8 expression and clinicopathological factors were further examined. Gene enrichment analysis was used to identify various putative biological activities and regulatory networks of RRP8 in HCC. The relationship between RRP8 expression and immune infiltration was confirmed by single-sample gene set enrichment analysis (ssGSEA). Univariate and multivariate Cox regression analyses were conducted to assess the impact of clinical variables on patient outcomes. Furthermore, a nomogram was constructed to estimate survival probability based on multivariate Cox regression analysis. Functional validation of RRP8 in HCC was performed with two different systems: doxycycline-inducible shRNA knockdown and CRISPR-Cas9 knockout.

RESULTS

RRP8 was markedly overexpressed in HCC clinical specimens compared to adjacent normal tissues. Further analysis demonstrated that RRP8 was directly connected to multiple clinical characteristics and strongly associated with various immune markers in HCC. Moreover, elevated RRP8 expression indicated an unfavourable prognosis. Our functional studies revealed that both knockdown and knockout of RRP8 dramatically attenuated liver cancer cells to proliferate and migrate. Knockout of RRP8 decreased the phosphorylation of MEK1/2 and β-catenin-(Y654) signalling pathway components; downregulated downstream signalling effectors, including Cyclin D1 and N-cadherin; and upregulated E-cadherin.

CONCLUSIONS

RRP8 is strongly implicated in immune infiltration and could be a potential therapeutic target in HCC.

Exploring gene-patient association to identify personalized cancer driver genes by linear neighborhood propagation

BACKGROUND

Driver genes play a vital role in the development of cancer. Identifying driver genes is critical for diagnosing and understanding cancer. However, challenges remain in identifying personalized driver genes due to tumor heterogeneity of cancer. Although many computational methods have been developed to solve this problem, few efforts have been undertaken to explore gene-patient associations to identify personalized driver genes.

RESULTS

Here we propose a method called LPDriver to identify personalized cancer driver genes by employing linear neighborhood propagation model on individual genetic data. LPDriver builds personalized gene network based on the genetic data of individual patients, extracts the gene-patient associations from the bipartite graph of the personalized gene network and utilizes a linear neighborhood propagation model to mine gene-patient associations to detect personalized driver genes. The experimental results demonstrate that as compared to the existing methods, our method shows competitive performance and can predict cancer driver genes in a more accurate way. Furthermore, these results also show that besides revealing novel driver genes that have been reported to be related with cancer, LPDriver is also able to identify personalized cancer driver genes for individual patients by their network characteristics even if the mutation data of genes are hidden.

CONCLUSIONS

LPDriver can provide an effective approach to predict personalized cancer driver genes, which could promote the diagnosis and treatment of cancer. The source code and data are freely available at https://github.com/hyr0771/LPDriver .

A novel endoplasmic reticulum stress-related lncRNA signature for prognosis prediction and immune response evaluation in Stomach adenocarcinoma

BACKGROUND

Stomach adenocarcinoma (STAD) is a significant contributor to cancer-related mortality worldwide. Although previous research has identified endoplasmic reticulum stress (ERS) as a regulator of various tumor-promoting properties of cancer cells, the impact of ERS-related long non-coding RNAs (lncRNAs) on STAD prognosis has not yet been investigated. Therefore, our study aims to develop and validate an ERS-related lncRNA signature that can accurately predict the prognosis of STAD patients.

METHODS

We collected RNA expression profiles and clinical data of STAD patients from The Cancer Genome Atlas (TCGA) and identified ERS-related genes from the Molecular Signature Database (MSigDB). Co-expression analysis enabled us to identify ERS-related lncRNAs, and we applied univariate Cox, least absolute shrinkage, and selection operator (LASSO), and multivariate Cox regression analyses to construct a predictive signature comprising of 9 ERS-related lncRNAs. We assessed the prognostic accuracy of our signature using Kaplan-Meier survival analysis, and validated our predictive signature in an independent gene expression omnibus (GEO) cohort. We also performed tumor mutational burden (TMB) and tumor immune microenvironment (TIME) analyses. Enrichment analysis was used to investigate the functions and biological processes of the signature, and we identified two distinct STAD patient subgroups through consensus clustering. Finally, we performed drug sensitivity analysis and immunologic efficacy analysis to explore further insights.

RESULTS

The 9 ERS related-lncRNAs signature demonstrated satisfactory predictive performance as an independent prognostic marker and was significantly associated with STAD clinicopathological characteristics. Furthermore, patients in the high-risk group displayed a worse STAD prognosis than those in the low-risk group. Notably, gene set enrichment analysis (GSEA) revealed significant enrichment of extracellular matrix pathways in the high-risk group, indicating their involvement in STAD progression. Additionally, the high-risk group exhibited significantly lower TMB expression levels than the low-risk group. Consensus clustering revealed two distinct STAD patient subgroups, with Cluster 1 exhibiting higher immune cell infiltration and more active immune functions. Drug sensitivity analysis suggested that the low-risk group was more responsive to oxaliplatin, epirubicinl, and other drugs.

CONCLUSION

Our study highlights the crucial regulatory roles of ERS-related lncRNAs in STAD, with significant clinical implications. The 9-lncRNA signature we have constructed represents a reliable prognostic indicator that has the potential to inform more personalized treatment decisions for STAD patients. These findings shed new light on the pathogenesis of STAD and its underlying molecular mechanisms, offering opportunities for novel therapeutic strategies to be developed for STAD patients.

BRAP silencing protects against neuronal inflammation, oxidative stress and apoptosis in cerebral ischemia-reperfusion injury by promoting PON1 expression

BACKGROUND

BRCA1 associated protein (BRAP) participates in the regulation of myocardial infarction and atherosclerosis. But the function of BRAP in cerebral ischemia-reperfusion (CIR) injury has not been elucidated yet.

METHODS

BRAP expression in PC12 cells in response to oxygen-glucose deprivation/reoxygenation (OGD/R) treatment was examined with Western blot assay. PC12 cells underwent OGD/R-treatment and were subsequently transfected with pcDNA-BRAP or sh-BRAP, followed by determination of viability, lactate dehydrogenase (LDH) production, apoptosis, inflammatory cytokine secretion, and oxidative stress marker protein levels. Paraoxonase 1 (PON1) promoter methylation was evaluated with methylation-specific PCR assay. the effect of BRAP/PON1 axis on CIR injury was investigated by rescue experiments. Additionally, sh-BRAP was injected into a middle cerebral artery occlusion (MCAO) rat model, and the changes of neurological damage were evaluated.

RESULTS

BRAP overexpression exacerbated OGD/R-induced viability reduction, LDH production, apoptosis, inflammatory cytokine secretion and oxidative stress in PC12 neuronal cells. In contrast, BRAP silencing showed the opposite results. Mechanistically, BRAP reduced PON1 expression by promoting DNA methyl transferase1 (DNMT1)-mediated PON1 promoter methylation. PON1 silencing reversed BRAP-mediated neuroprotection. Additionally, BRAP silencing alleviated CIR-induced neurological damage in MCAO rats.

CONCLUSION

BRAP silencing suppressed OGD/R-induced neuronal apoptosis, inflammation, and oxidative stress, and alleviated CIR-induced neurological damage in MCAO rats through facilitating PON1 expression.

Decoding the significant diagnostic and prognostic importance of maternal embryonic leucine zipper kinase in human cancers through deep integrative analyses

BACKGROUND

Cancer is a multifactorial disease and the second leading cause of human deaths worldwide. So far, the underlying mechanisms of cancer have not been yet fully elucidated.

METHODS

By using TCGA expression data, we determine the pathogenic roles of the maternal embryonic leucine zipper kinase (MELK) gene in various human cancers in this study. For this purpose, different online databases and tools (UALCAN, Kaplan-Meier (KM) plotter, TNMplot, GENT2, GEPIA, HPA, cBioPortal, STRING, Enrichr, TIMER, Cytoscape, DAVID, MuTarget, and CTD) were used.

RESULTS

MELK gene expression was analyzed in a total of 24 human cancers and was found notably up-regulated in all the 24 analyzed tumor tissues relative to controls. Moreover, across a few specific cancers, including kidney renal clear cell carcinoma (KIRC), stomach adenocarcinoma (STAD), lung adenocarcinoma (LUAD), and liver hepatocellular carcinoma (LIHC) patients, MELK up-regulation was observed to be correlated with the shorter survival duration and metastasis. This valuable information highlighted that MELK plays a significant role in the development and progression of these four cancers. Based on clinical variables, MELK higher expression was also found in KIRC, STAD, LUAD, and LIHC patients with different clinical variables. Gene ontology and pathway analysis outcomes showed that MELK-associated genes notably co-expressed with MELK and belongs to a variety of diverse biological processes, molecular functions, and pathways. MELK expression was also correlated with promoter methylation levels, genetic alterations, other mutant genes, tumor purity, CD8+ T, and CD+4 T immune cells infiltrations in KIRC, STAD, LUAD, and LIHC.

CONCLUSION

This pan-cancer study revealed the diagnostic and prognostic roles of MELK across four different cancers.

The construction of a prognostic model of cervical cancer based on four immune-related LncRNAs and an exploration of the correlations between the model and oxidative stress

Introduction: The immune-related lncRNAs (IRLs) are critical for the development of cervical cancer (CC), but it is still unclear how exactly ILRs contribute to CC. In this study, we aimed to examine the relationship between IRL and CC in detail. Methods: First, the RNAseq data and clinical data of CC patients were collected from The Cancer Genome Atlas (TCGA) database, along with the immune genes from the Import database. We used univariate cox and least absolute shrinkage and selection operator (lasso) to obtain IRLs for prediction after screening the variables. According to the expression levels and risk coefficients of IRLs, the riskscore were calculated. We analyzed the relationship between the model and oxidative stress. We stratified the risk model into two as the high and low-risk groups. We also evaluated the survival differences, immune cell differences, immunotherapeutic response differences, and drug sensitivity differences between the risk groups. Finally, the genes in the model were experimentally validated. Results: Based on the above analyses, we further selected four IRLs (TFAP2A.AS1, AP000911.1, AL133215.2, and LINC02078) to construct the risk model. The model was associated with oxidative-stress-related genes, especially SOD2 and OGG1. Patients in the high-risk group had a lower overall survival than those in the low-risk group. Riskscore was positively correlated with resting mast cells, neutrophils, and CD8+ T-cells. Patients in the low-risk group showed a greater sensitivity to immunosuppression therapy. In addition, we found that patients with the PIK3CA mutation were more sensitive to chemotherapeutic agents such as dasatinib, afatinib, dinaciclib and pelitinib. The function of AL133215.2 was verified, which was consistent with previous findings, and AL133215.2 exerted a pro-tumorigenic effect. We also found that AL133215.2 was closely associated with oxidative-stress-related pathways. Discussion: The results suggested that risk modeling might be useful for prognosticating patients with CC and opening up new routes for immunotherapy.

Comprehensive analysis of the oncogenic roles of vascular endothelial growth factors and their receptors in stomach adenocarcinoma

Background

Vascular endothelial growth factors (VEGFs) and their receptors (VEGFRs) play complicated oncogenic roles in multiple tumors by initiating and promoting tumor angiogenesis and lymphangiogenesis. The main goal of our study was to comprehensively investigate the oncogenic roles of VEGFs and VEGFRs in stomach adenocarcinoma (STAD).

Methods

The present study applied multiple bioinformatic tools to comprehensively explore the expression levels, prognostic values, genetic alterations and immune infiltrations of VEGFs and VEGFRs in STAD patients.

Results

We found that VEGFA, VEGFC, placenta growth factor, FLT1, KDR, FLT4, and Neuropilin 1 were overexpressed in STAD, while the expression of VEGFB and VEGFD were decreased. Survival analysis revealed that higher transcription levels of VEGF/VEGFRs were obviously correlated with worse clinical outcome in STAD patients. Additionally, high alteration frequencies of VEGFs and VEGFRs (27%) were observed in STAD patients, and alterations of VEGFs and VEGFRs improved their prognosis. The expression of VEGFs and VEGFRs was remarkably associated with immune cell infiltration and immune checkpoint expression in STAD patients.

Conclusion

Our study systematically explored the transcriptome profiles and distinct prognostic values of VEGFs and their receptors in STAD and contributed to a better understanding of the oncogenic roles of VEGF/VEGFR members in STAD.

Comprehensive pan-cancer analysis of PTGES3 and its prognostic role in hepatocellular carcinoma

Introduction

PTGES3, also known as p23, is a molecule chaperone of Hsp90 that is involved in the pathogenesis of malignant tumors. Increasing studies have shown that PTGES3 plays a nonnegligible role in tumor development. However, analysis of PTGES3 in pan-cancer has not been performed yet.

Methods

We explored the role of PTGES3 in 33 types of tumors and depicted the potentialimmune-related pathways among them. Using multiple databases includingTCGA, LinkedOmics, GDSC, and TIMER, we made a comprehensive analysis to explore whether there was an interaction between PTGES3 and prognosis, DNA methylation, copy number variation (CNV), tumor mutational burden (TMB), microsatellite instability (MSI), and tumor immune microenvironment (TME).

Results

Our study revealed that PTGES3 expression level was upregulated in most cancers. PTGES3 was also associated with a positive or negative prognosis in a variety of cancers, which was mainly associated with DNA methylation, CNV, MSI, TMB, andmismatch repair-related genes. High PTGES3 expression was related to the infiltration of Th2 subsets of CD4+ T cells and immune checkpoint-related genes in most cancers, especially in hepatocellular carcinoma (HCC). Enrichment analysis demonstrated that PTGES3 was involved in cellular processes including DNA replication and spliceosome. The relationship between PTGES3 expression and HCC progression was verified at the protein level through immune histochemical analysis.

Conclusions

Our research demonstrated theprognostic predictive value of PTGES3 in a wide range of cancers, which was alsoassociated with the process of tumor immune infiltration. As a result, it suggestedthat PTGES3 was a valuable prognostic biomarker in HCC treatment.

FOXP3 and SQSTM1/P62 correlate with prognosis and immune infiltration in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most common highly malignant tumours worldwide. FOXP3 and SQSTM1/P62 have been shown to be abnormally expressed in tumour cells, but their function in different tumours remains controversial. The present study was designed to evaluate the expression of FOXP3 and P62 in HCC and their prognostic value as well as their relationship with immune infiltration in HCC patients.

METHODS

The Gene Expression Omnibus (GEO) database and TNMplot.com platform were used to analyse the expression of FOXP3 and P62. The Cancer Genome Atlas (TCGA) database and Kaplan-Meier plotter were used to assess the impacts of FOXP3 and P62 on clinical prognosis. In addition, TCGA database was also used to examine the correlation between the expression of FOXP3 and P62 and tumour immune infiltration using the CIBERSORT algorithm. Finally, immunohistochemistry (IHC) was used to determine expression levels of FOXP3 and P62 in 89 HCC and adjacent normal liver tissues, and their effects on clinicopathological features and prognosis were verified.

RESULTS

FOXP3 expression was downregulated in HCC tissues, while P62 expression was upregulated. FOXP3 underexpression and P62 overexpression were closely related to decreased overall survival (OS) in HCC patients. Additionally, the abnormal expression of FOXP3 and P62 was closely related to the infiltration levels of 12 types of immune cells, including regulatory T cells (Tregs), M2 macrophages, M0 macrophages, and CD8 T cells. Notably, in the validation model, abnormal FOXP3 and P62 expression was significantly associated with adverse clinicopathological factors in HCC patients, including elevated α-fetoprotein (AFP) levels, poor tumour differentiation, and increased Ki67 levels. Furthermore, low FOXP3 and high P62 expression were independent risk factors for predicting OS prognosis in HCC patients.

CONCLUSION

FOXP3 and P62 have been shown to be important prognostic factors in HCC patients and are associated with immune cell infiltration in HCC. These findings suggest that FOXP3 and P62 may be valuable prognostic biomarkers and potential therapeutic targets for HCC treatment.

The role of toll-like receptors (TLRs) in pan-cancer

BACKGROUND

Toll-like receptors (TLRs) are important components of the innate and adaptive immune systems, and abnormal TLR expression has been linked to a variety of cancers. However, there was a lack of clarity on the association of TLR stimulation with the carcinogenesis of cancer. The study's goal was to analyse the clinical importance of TLRs expression at the mRNA level in pan-cancer datasets, as well as the link between TLR expression and carcinogenesis, progression, and clinical prognosis.

METHODS

The expression profile of TLRs derived from UCSC pan-cancer data was analysed in multiple dimensions, including clinical analysis, immunological subtype analysis, tumour microenvironment (TME) analysis, tumour stem cell correlation analysis, and drug sensitivity analysis. Additionally, we analyse protein-protein interactions, functional enrichment, and chromatin accessibility, as well as TLR expression in single-cell sequencing data.

RESULTS

Our multi-omics analysis results imply that TLRs may operate as a biological marker for carcinogenesis and progression, a potential target for anti-tumour therapy, and a prognostic biomarker, laying the theoretical groundwork for future translational medicine research.

CONCLUSION

TLRs are involved in the formation of malignancies and can be explored in further detail as potential prognostic indicators. Key MessagesToll-like receptors (TLRs) are key factors in the process of the innate and adaptive immune response, and their aberrant expression of TLRs have been widely reported in various cancer. However, the association between TLRs stimulation and tumorigenesis of cancer has not been well clarified.In this study, in the pan-cancer data, integrated TLR family gene expression analysis, clinical correlation analysis, immune subtype correlation analysis, tumour microenvironment correlation analysis, tumour stem cell correlation analysis, and drug sensitivity correlation analysis were performed.TLRs play an important role in the development of tumours and can be studied in depth as potential prognostic markers.

Pan-cancer analysis reveals DDX21 as a potential biomarker for the prognosis of multiple tumor types

Background

DExD-box helicase 21 (DDX21) is an essential member of the RNA helicase family. DDX21 is involved in the carcinogenesis of various malignancies, but there has been no comprehensive research on its involvement in different types of cancer.

Method

This study used TCGA, CPTAC, GTEx, GEO, FANTOM5, BioGRID, TIMER2, GEPIA2, cBioPortal, STRING, and Metascape databases and Survival ROC software to evaluate DDX21 gene expression, protein expression, immunohistochemistry, gene mutation, immune infiltration, and protein phosphorylation in 33 TCGA tumor types, as well as the prognostic relationship between DDX21 and different tumors, by survival analysis and similar gene enrichment analysis. Furthermore, Cell Counting Kit-8 (CCK-8) and Transwell studies were employed to assess the effect of DDX21 expression on lung adenocarcinoma (LUAD) cell proliferation and migration.

Result

The DDX21 gene was highly expressed in most cancers, and overexpression was associated with poor overall survival (OS) and disease-free survival (DFS). DDX21 mutations were most common in uterine corpus endometrial carcinoma (UCEC; >5%), and DDX21 expression was positively correlated with the degree of infiltration of CAF and CD8⁺ cells in several tumor types. Numerous genes were co-expressed with DDX21. Gene enrichment analysis revealed close links between DDX21, RNA metabolism, and ribosomal protein production. In vitro analysis of LUAD cells showed that DDX21 expression was positively correlated with cell proliferation and migration capacity, consistent with prior bioinformatics studies.

Conclusions

DDX21 is overexpressed in a variety of cancers, and overexpression in some cancers is associated with poor prognosis. Immune infiltration and DDX21-related gene enrichment analyses indicated that DDX21 may affect cancer development through mechanisms that regulate tumor immunity, RNA metabolism, and ribosomal protein synthesis. This pan-cancer study revealed the prognostic value and the oncogenic role of DDX21.

Identification of ACSF gene family as therapeutic targets and immune-associated biomarkers in hepatocellular carcinoma

Acyl-CoA synthetases (ACSs) are responsible for acyl-CoA synthesis from nonpolar hydrophilic fatty acids and play a vital role in many metabolic processes. As a category of ACS isozymes, members of ACS family (AACS, ACSF2-3, AASDH) participate in lipid metabolism; however, their expression patterns, regulatory mechanisms and effects in hepatocellular carcinoma (HCC) are poorly understood. Here, through evaluating the expression profiles of ACSF gene family, we found that upregulated AACS might be more significant and valuable in development and progression of HCC. Consequently, the mRNA expression levels of AACS and ACSF2 was accordantly increased in HCC. Kaplan-Meier plotter revealed that HCC patients with high level of AACS were highly related to a shorter overall survival time and relapse-free survival. Genetic alterations using cBioPortal revealed that the alteration rate of AACS were 5%. We also found that the functions of ACSF gene family were linked to several cancer-associated pathways, including long-term potentiation, phospholipase D signaling pathway and purine metabolism. TIMER database indicated that the AACS and ACSF2 had a strong relationship with the infiltration of six types of immune cells (macrophages, neutrophils, CD8+ T-cells, B-cells, CD4+ T-cells and dendritic cells). Next, Diseasemeth database revealed that the global methylation levels of ACSF2 was higher in HCC patients. In conclusion, this study firstly demonstrated that Acyl-CoA synthesis gene family, in particular, AACS, could be associated with immune microenvironment, thereby influencing the development and prognosis of patients with HCC.

A novel prognostic index of stomach adenocarcinoma based on immunogenomic landscape analysis and immunotherapy options

Stomach adenocarcinoma (STAD) is one of the most common malignant tumors worldwide. In this study, we attempted to construct a valid immune-associated gene prognostic index risk model that can predict the survival of patients with STAD and the efficacy of immune checkpoint inhibitors (ICIs) treatment. Transcriptome, clinical, and gene mutational data were obtained from the TCGA database. Immune-related genes were downloaded from the ImmPort and InnateDB databases. A total of 493 immune-related genes were identified to be enriched in functions associated with immune response, as well as in immune and tumor-related pathways. Further, 36 candidate genes related to the overall survival (OS) of STAD were obtained by weighted gene co-expression network analysis (WGCNA). Next, based on a Cox regression analysis, we constructed an immune-associated gene prognostic index (IAGPI) risk model based on eight genes, which was verified using the GEO STAD cohort. The patients were divided into two subsets according to their risk score. Patients in the low-risk group had better OS than those in the high-risk group. In the low-risk group, there were more CD8, activated memory CD4, and follicular helper T cells, and M1 macrophages, whereas monocytes, M2 macrophages, eosinophils, and neutrophils were more abundant in the high-risk group. The patients in the low-risk group were more sensitive to ICIs therapy. The IAGPI risk model can precisely predict the prognosis, reflect the tumor immune microenvironment, and predict the efficacy of ICIs therapy in patients with STAD.

MND1 functions as a potential prognostic biomarker associated with cell cycle and immune infiltration in kidney renal clear cell carcinoma

Kidney renal clear cell carcinoma (KIRC) is a common and invasive subtype of renal tumors, which has poor prognosis and high mortality. MND1 is a meiosis specific protein that participates in the progress of diverse cancers. Nonetheless, its function in KIRC was unclear. Here, TIMER, TCGA, GEO databases and IHC found MND1 expression is upregulated in KIRC, leading to poor overall survival, and MND1 can serve as an independent prognostic factor. Moreover, enrichment analysis revealed the functional relationship between MND1 and cell cycle, immune infiltration. EdU and transwell assays confirmed that MND1 knockdown surely prohibited the proliferation, migration, and invasion of KIRC cells. Additionally, immune analysis showed that MND1 displayed a strong correlation with various immune cells. Interference with MND1 significantly reduces the expression of chemokines. TCGA and GEO databases indicated that MND1 expression is significantly related to two m6A modification related gene (METTL14, IGF2BP3). Finally, the drug sensitivity analysis revealed 7 potentially sensitive drugs for KIRC patients with high MND1 expression. In conclusion, MND1 can be used as a prognostic biomarker for KIRC and provides clues regarding cell cycle, immune infiltrates and m6A. Sensitive drugs may be an effective treatment strategy for KIRC patients with high expression of MND1.

Identification of cervical squamous cell carcinoma feature genes and construction of a prognostic model based on immune-related features

As heterogeneity of cervical squamous cell carcinoma (CSCC), prognosis assessment for CSCC patients remain challenging. To develop novel prognostic strategies for CSCC patients, associated biomarkers are urgently needed. This study aimed to cluster CSCC samples from a molecular perspective. CSCC expression data sets were obtained from The Cancer Genome Atlas and based on the accessed expression profile, a co-expression network was constructed with weighted gene co-expression network analysis to form different gene modules. Tumor microenvironment was evaluated using ESTIMATE algorithm, observing that the brown module was highly associated with tumor immunity. CSCC samples were clustered into three subtypes by consensus clustering based on gene expression profiles in the module. Gene set variation analysis showed differences in immune-related pathways among the three subtypes. CIBERSORT and single-sample gene set enrichment analysis analyses showed the difference in immune cell infiltration among subtype groups. Also, Human leukocyte antigen protein expression varied considerably among subtypes. Subsequently, univariate, Lasso and multivariate Cox regression analyses were performed on the genes in the brown module and an 8-gene prognostic model was constructed. Kaplan-Meier analysis illuminated that the low-risk group manifested a favorable prognosis, and receiver operating characteristic curve showed that the model has good predictive performance. qRT-PCR was used to examine the expression status of the prognosis-associated genes. In conclusion, this study identified three types of CSCC from a molecular perspective and established an effective prognostic model for CSCC, which will provide guidance for clinical subtype identification of CSCC and treatment of patients.

Integrated analysis of anti-tumor roles of BAP1 in osteosarcoma

Background

This study aims to screen out differentially expressed genes (DEGs) regulated by BRCA1-associated protein 1 (BAP1) in osteosarcoma cells, and to analyze their biological functions.

Methods

The microarray dataset GSE23035 of BAP1-knockdown osteosarcoma cells was obtained from Gene Expression Omnibus (GEO) database, consisting of shControl, shBAP1#1 and shBAP1#2 samples. The DEGs between the BAP1-knockdown osteosarcoma cells and the untreated osteosarcoma cells were screened with limma package, and then subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Gene Set Enrichment Analysis (GSEA) was also performed for the three groups of samples. Hub genes in a protein-protein interaction (PPI) network of DEGs was filtered, and then subjected to prognostic analysis and correlation analysis with BAP1 in Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. Besides, the correlation between BAP1 and biological processes/pathways was analyzed by Gene Set Variation Analysis (GSVA) method and the correlation between BAP1 and immune infiltration by CIBERSORT and ESTIMATE methods. The roles of BAP1 in regulating proliferation and epithelial-mesenchymal transition (EMT) were validated by CCK-8 and western blot.

Results

58 upregulated DEGs and 81 downregulated DEGs were obtained with |logFC| ≥ 1 and adj.p < 0.05. Cell cycle, DNA repair, and focal adhesion were associated with BAP1 in datasets. Further, BAP1 was negatively correlated with naïve CD4 T cells infiltration. In vitro, BAP1 inhibited proliferation and EMT.

Conclusion

BAP1 might be a tumor suppressor in osteosarcoma and a promising therapeutic target.

A comprehensive analysis of gasdermin family gene as therapeutic targets in pan-cancer

Six members of the gasdermin family are involved in various biological functions in malignant tumors. The present study aimed to perform a comprehensive analysis of gasdermin family genes in pan-cancer. Raw data was acquired from the genotype-tissue expression (GTEx) and the Cancer Genome Atlas. High inter-tumor heterogeneity in the expression between paracancerous and tumor tissues was observed across cancers. Survival analysis confirmed that the risk or protective effects of gasdermin family members on prognosis depended on the cancer types. The mutation frequency appeared to be high, and the mutation group had a worse prognosis. Besides, gasdermin family genes were associated with immune infiltrate subtypes, stromal and immune cell infiltration levels, TMB, MSI, immune checkpoint gene expression, and tumor stemness scores. Moreover, gasdermin family gene expressions affected the expressions of MMR genes and methyltransferases and could predict cancer cells sensitivity to chemotherapeutic drugs. Subsequently, the findings were double-checked in LIHC and PAAD. GSEA results indicated the gasdermin family genes mainly involved in tumor metabolism and immune microenvironment remodeling related signaling pathways. In conclusion, our findings confirmed that gasdermin family genes were potential therapeutic cancer targets in pan-cancer.

A Pan-Cancer Analysis Reveals CLEC5A as a Biomarker for Cancer Immunity and Prognosis

Background

CLEC5A is a member of the C-type lectin superfamily. It can activate macrophages and lead to a series of immune-inflammation reactions. Previous studies reveal the role of CLEC5A in infection and inflammation diseases.

Method

We acquire and analyze data from The Cancer Genome Atlas (TCGA) database, Genotype-Tissue Expression (GTEx) database, and other comprehensive databases via GSCALite, cBioPortal, and TIMER 2.0 platforms or software. Single-cell sequencing analysis was performed for quantifying the tumor microenvironment of several types of cancers.

Results

CLEC5A is differentially expressed in a few cancer types, of which overexpression accompanies low overall survival of patients. DNA methylation mainly negatively correlates with CLEC5A expression. Moreover, CLEC5A is positively related to immune infiltration, including macrophages, cancer-associated fibroblasts (CAFs), and regulatory T cells (Tregs). Immune checkpoint genes are significantly associated with CLEC5A expression in diverse cancers. In addition, CLEC5A expression correlates with mismatch repair (MMR) in several cancers. Tumor mutation burden (TMB), microsatellite instability (MSI), and neoantigens show a positive association with CLEC5A expression in several cancers. Furthermore, CLEC5A in cancer correlates with signal transduction, the immune system, EMT, and apoptosis process. The drug sensitivity analysis screens out potential therapeutic agents associated with CLEC5A expression, including FR-180204, Tivozanib, OSI-930, Linifanib, AC220, VNLG/124, Bexarotene, omacetaxine mepesuccinate, narciclasine, leptomycin B, PHA-793887, LRRK2-IN-1, and CR-1-31B.

Conclusion

CLEC5A overexpresses in multiple cancers in contrast to normal tissues, and high CLEC5A expression predicts poor prognosis of patients and immune infiltration. CLEC5A is a potential prognostic biomarker of diverse cancers and a target for anti-tumor therapy.

Identification and validation of necroptosis-related prognostic gene signature and tumor immune microenvironment infiltration characterization in esophageal carcinoma

Background

Esophageal carcinoma (ESCA) is a common malignancy with a poor prognosis. Previous research has suggested that necroptosis is involved in anti-tumor immunity and promotes oncogenesis and cancer metastasis, which in turn affects tumor prognosis. However, the role of necroptosis in ESCA is unclear. This study aimed to investigate the relationships between necroptosis-related genes (NRGs) and ESCA.

Methods and results

The clinical data and gene expression profiles of ESCA patients were extracted from The Cancer Genome Atlas (TCGA), and 159 NRGs were screened from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We then identified 52 differentially expressed NRGs associated with ESCA and used them for further analysis. Gene ontology (GO) and KEGG functional enrichment analyses showed that these NRGs were mostly associated with the regulation of necroptosis, Influenza A, apoptosis, NOD-like receptor, and NF-Kappa B signaling pathway. Next, univariate and multivariate Cox regression and LASSO analysis were used to identify the correlation between NRGs and the prognosis of ESCA. We constructed a prognostic model to predict the prognosis of ESCA based on SLC25A5, PPIA, and TNFRSF10B; the model classified patients into high- and low-risk subgroups based on the patient’s risk score. Furthermore, the receiver operating characteristic (ROC) curve was plotted, and the model was affirmed to perform moderately well for prognostic predictions. In addition, Gene Expression Omnibus (GEO) datasets were selected to validate the applicability and prognostic value of our predictive model. Based on different clinical variables, we compared the risk scores between the subgroups of different clinical features. We also analyzed the predictive value of this model for drug sensitivity. Moreover, Immunohistochemical (IHC) validation experiments explored that these three NRGs were expressed significantly higher in ESCA tissues than in adjacent non-tumor tissues. In addition, a significant correlation was observed between the three NRGs and immune-cell infiltration and immune checkpoints in ESCA.

Conclusions

In summary, we successfully constructed and validated a novel necroptosis-related signature containing three genes (SLC25A5, PPIA, and TNFRSF10B) for predicting prognosis in patients with ESCA; these three genes might also play a crucial role in the progression and immune microenvironment of ESCA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12876-022-02423-6.

Identification of Potential Biomarkers for Pan-Cancer Diagnosis and Prognosis Through the Integration of Large-Scale Transcriptomic Data

Cancer is one of the leading causes of death worldwide, bringing a significant burden to human health and society. Accurate cancer diagnosis and biomarkers that can be used as robust therapeutic targets are of great importance as they facilitate early and effective therapies. Shared etiology among cancers suggests the existence of pan-cancer biomarkers, performance of which could benefit from the large sample size and the heterogeneity of the studied patients. In this study, we conducted a systematic RNA-seq study of 9,213 tumors and 723 para-cancerous tissue samples of 28 solid tumors from the Cancer Genome Atlas (TCGA) database, and 7,008 normal tissue samples from the Genotype-Tissue Expression (GTEx) database. By differential gene expression analysis, we identified 214 up-regulated and 186 downregulated differentially expressed genes (DEGs) in more than 80% of the studied tumors, respectively, and obtained 20 highly linked up- and downregulated hub genes from them. These markers have rarely been reported in multiple tumors simultaneously. We further constructed pan-cancer diagnostic models to classify tumors and para-cancerous tissues using 10 up-regulated hub genes with an AUC of 0.894. Survival analysis revealed that these hub genes were significantly associated with the overall survival of cancer patients. In addition, drug sensitivity predictions for these hub genes in a variety of tumors obtained several broad-spectrum anti-cancer drugs targeting pan-cancer. Furthermore, we predicted immunotherapy sensitivity for cancers based on tumor mutational burden (TMB) and the expression of immune checkpoint genes (ICGs), providing a theoretical basis for the treatment of tumors. In summary, we identified a set of biomarkers that were differentially expressed in multiple types of cancers, and these biomarkers can be potentially used for diagnosis and used as therapeutic targets.

Aberrant Expression of ADARB1 Facilitates Temozolomide Chemoresistance and Immune Infiltration in Glioblastoma

Chemoresistance, especially temozolomide (TMZ) resistance, is a major clinical challenge in the treatment of glioblastoma (GBM). Exploring the mechanisms of TMZ resistance could help us identify effective therapies. Adenosine deaminases acting on RNA (ADARs) are very important in RNA modification through regulating the A-to-I RNA editing. Recent studies have shown that ADARs regulate multiple neurotransmitter receptors, which have been linked with the occurrence and progress of GBM. Here, data from several bioinformatics databases demonstrated that adenosine deaminase RNA specific B1 (ADARB1), also named ADAR2, was upregulated in both GBM tissues and cells, and had the prognostic value in GBM patients. Moreover, ADARB1 was found to be involved in AKT-mediated TMZ resistance in GBM cells. The KEGG analysis of ADARB1-associated co-expressed genes showed that ADARB1 was potentially involved in the mitochondrial respiratory chain complex. TISIDB and GEPIA databases were further used to analyze the role of ADARB1 in tumor-immune system interactions in GBM. These findings deepened our understanding of the function of ADARB1 in tumorigenesis and therapeutic response in GBM.

N6-Methyladenosine Methylation Regulator RBM15 is a Potential Prognostic Biomarker and Promotes Cell Proliferation in Pancreatic Adenocarcinoma

RNA binding motif protein 15 (RBM15) is a key regulatory factor involved in N6-methyladenosine (m6A) methylation. It has been reported that RBM15 plays an important role in the progress of laryngeal squamous cell carcinoma (LSCC), promoting LSCC migration and invasion. However, the role of RBM15 in human different cancers remains unknown. This study aims to analyze the prognostic value of RBM15, and to demonstrate the correlation between RBM15 expression and tumor immunity, as well as to provide clues for further mechanism research. The results showed that RBM15 was mutated or copy number varied in 25 types of cancer. RBM15 mRNA was abnormally up-regulated across various cancers. Survival analysis suggested high expression of RBM15 was associated with poor prognosis in many cancer types. Among these, it affected patients’ overall survival (OS) in 10 cancer types, disease-free interval (DFI) in 8 cancer types, progression-free interval (PFI) in 12 cancer types and disease-specific survival (DSS) in 7 cancer types. Importantly, in pancreatic adenocarcinoma (PAAD), overexpression of RBM15 is associated with patients’ OS, DFI, PFI, or DSS. In addition, RBM15 expression was positively correlated with immune infiltrating cells in kidney renal clear cell carcinoma (KIRC), brain lower grade glioma (LGG), and PAAD. Moreover, RBM15 expression showed a strong correlation with immune checkpoint markers in PAAD. Cell counting kit-8 (CCK-8) assay showed that knockdown of RBM15 significantly inhibited the proliferation of pancreatic cancer cells. PPI analysis showed USP10, USP24, SMG1, NRAS were closely connected with RBM15 alterations. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that many biological processes (BP), cellular components (CC), molecular functions (MF), cancer related pathways including “sister chromatid cohesion”, “peptidyl-serine phosphorylation”, “cell division”, “nucleoplasm”, “nucleus”, “protein binding”, “protein serine/threonine kinase activity”, “T cell receptor signaling pathway”, “Cell cycle” were regulated by RBM15 alterations. Taken together, pan-cancer analysis of RBM15 suggested it may be served as a prognostic biomarker and immunotherapeutic target for PAAD.

ACSL4 Expression Is Associated With CD8+ T Cell Infiltration and Immune Response in Bladder Cancer

Objective

This study aimed to explore the role of ACSL4 in CD8+ T cell tumor infiltration and outcomes of bladder cancer (BLCA) patients after immunotherapy.

Methods

The correlation between ACSL4 expression and tumor infiltration of immune cells was analyzed using the Tumor Immune Estimation Resource database. The prognostic significance of ACSL4 in BLCA was analyzed using Kaplan–Meier curves. Immunohistochemistry was used to detect CD8+ T cell infiltration in tumors with high and low ACSL4 expression obtained from patients at the Fudan University Shanghai Cancer Center. The relationships between immune checkpoint genes and immune response were analyzed using The Cancer Genome Atlas and IMvigor 210 cohorts. The molecular functions, cellular components, and biological processes involving ACSL4 were explored using Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment pathway analyses.

Results

The expression level of ACSL4 was significantly correlated with the infiltration of CD8+ T cells in BLCA tumors (r = 0.192, P = 2.22e-04). Elevated ACSL4 was associated with suppressed tumor progression and better outcomes for BLCA patients. The higher expression level of ACSL4 predicted better immunotherapeutic responses and was associated with higher expression levels of core immune checkpoint genes, including CD274, CTLA4, PDCD1, and LAG3, compared with the low ACSL4 expression group.

Conclusion

This study demonstrated for the first time that elevated ACSL4 correlated significantly with CD8+ T cell infiltration and contributed to better immunotherapeutic responses in BLCA patients. Furthermore, ACSL4 serves as a novel biomarker for predicting patient outcomes after immunotherapeutic treatments, which may improve the development of individualized immunotherapy for BLCA.

YTH domain family: potential prognostic targets and immune-associated biomarkers in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common high malignancy with insidious onset, invasive fast-growing, high recurrence rate and fatality. YTH domain family plays essential roles in development of HCC. However, the biological function of YTH domain family in HCC have not been clarified. Here, through evaluating the expression profiles of YTH domain family, we found that upregulated YTHDF1 might be more significant and valuable in development and progression of HCC. There was a strong correlation between YTHDC1, YTHDF1 and YTHDF2 and pathological stage of HCC patients. Kaplan-Meier plotter revealed that HCC patients with high level of YTHDF1 and YTHDF2 were highly related to a shorter overall survival time, and low level of YTHDF1 (p = 0.0017) has an important association with a longer progression-free survival time. Genetic alterations using cBioPortal revealed that the alteration rates of YTHDF3 were the highest. We also found that the functions of YTH domain family were linked to several cancer-associated pathways, including peptidyl-serine modification, peptidyl-tyrosine modification and negative regulation of cellular component movement. TIMER database indicated that the YTH domain family had a strong relationship with the infiltration of six types of immune cells (macrophages, neutrophils, CD8+ T-cells, B-cells, CD4+ T-cells and dendritic cells). Next, Ualcan databases revealed that the global methylation levels of YTHDC1 was higher in HCC patients, while YTHDF2 was lower in HCC patients. In conclusion, our findings will enhance the understanding of YTH domain family in HCC pathology, and provide novel insights into YTH-targeted therapy for HCC patients.

Pan-Cancer Analysis Reveals FH as a Potential Prognostic and Immunological Biomarker in Lung Adenocarcinoma

Fumarate hydratase (FH) is an important enzymatic component in the tricarboxylic acid cycle. Studies have reported that FH plays an important role in hereditary leiomyomatosis and renal cell cancer (HLRCC). However, the role of FH in human different cancers remains unknown. This study is aimed at analyzing the prognostic value of FH and demonstrating the correlation between FH expression and tumor immunity. Results showed that FH was mutated or copy number varied in 27 types of cancer. FH mRNA was abnormally upregulated across various cancers. Survival analysis suggested high expression of FH was associated with poor prognosis in many cancer types, including lung adenocarcinoma (LUAD). Additionally, FH expression was associated with immune infiltration, including B cells, CD4⁺ T cells, CD8⁺ T cells, neutrophils, macrophages, and dendritic cells, especially in liver hepatocellular carcinoma (LIHC), LUAD, and lung squamous cell carcinoma (LUSC). Moreover, FH expression showed a strong correlation with immune checkpoint markers in LUAD and testicular germ cell tumors (TGCT). These results indicate that FH is an immunotherapeutic target and a potential prognostic biomarker in LUAD.

Expression, prognosis value, and immune infiltration of lncRNA ASB16-AS1 identified by pan-cancer analysis

Long non-coding RNA known as ASB16 antisense RNA1 (ASB16-AS1) has been proven to be an oncogene, and the relationship between ASB16-AS1 and immunity is still under studied. This study aims to explore the expression and prognostic potential of ASB16-AS1, and to visualize the relationship between ASB16-AS1 expression and immune infiltration in pan-cancer analysis. We clarified ASB16-AS1 expression patterns and its relationship with prognosis through multi-platform and multi-database sources. We also verified the function of ASB16-AS1 in liver hepatocellular carcinoma (LIHC). A  variety of immune cell content evaluation methods were used to mutually verify the correlation between ASB16-AS1 and immune infiltration. Finally, the relationships between ASB16-AS1 and molecular characteristics were further explored. In terms of comprehensive analysis, compared with non-tumor tissues, ASB16-AS1 was highly expressed in tumor tissues, and indicated the value of poor prognosis in multiple cancer types. Functional assays, such as counting kit-8 assay, transwell assay and scratch-wound assay verified that high ASB16-AS1 expression promoted tumor progression in LIHC. ASB16-AS1 was positively correlated with B cells, T cells CD4+ and T cells CD8+ in most cancer types, and negatively correlated with macrophages, dendritic cells and neutrophils in some cancer types. In addition, there were different interaction modes between ASB16-AS1 and molecular features, such as the relationship with oncogenic signaling pathways, showing that the high ASB16-AS1 expression was related to alterations in oncogenic signaling pathways. Our study emphasizes that ASB16-AS1 is a potential pan-cancer prognostic marker, whichs is associated with the immune infiltration in multiple cancer types.

Integrated Profiling Identifies PLOD3 as a Potential Prognostic and Immunotherapy Relevant Biomarker in Colorectal Cancer

Procollagen-Lysine,2-Oxoglutarate 5-Dioxygenase 3 (PLOD3) is related to a variety of human diseases. However, its function in Colorectal cancer (CRC) remains uncertain. PLOD3 expression was analyzed using The Cancer Genome Atlas (TCGA) pan-cancer data. DAVID was used for enrichment analysis of PLOD3-related genes. The correlation between PLOD3 expression and immune cell infiltration was evaluated. Four expression profile datasets (GSE17536, GSE39582, GSE74602, and GSE113513) from Gene Expression Omnibus, and two proteomic datasets were used as validation cohorts for assessing the diagnostic and prognostic value of PLOD3 in CRC. What's more, we performed immunohistochemistry (IHC) staining for PLOD3 in 160 paired CRC specimens and corresponding adjacent non-tumor tissues. PLOD3 was highly expressed in many tumors including CRC. PLOD3 was upregulated in advanced stage CRCs, and high PLOD3 expression was associated with poor survival. High PLOD3 expression was associated with low levels of B cells, CD4⁺ T cells, M1 macrophages, CD8⁺ T cells, and multiple immunerelated characteristics. In addition, the high PLOD3 expression group had a higher TIDE score and a lower tumor mutation burden and microsatellite instability, indicating that patients with high PLOD3 expression may be resistant to immunotherapy. Additional datasets and IHC analysis were used to validate the diagnostic and prognostic value of PLOD3 at the mRNA and protein levels in CRC. Patients with non-response to immunotherapy showed increased PLOD3 expression in an immunotherapy treated dataset. PLOD3 is a potential biomarker for CRC diagnosis and prognosis prediction. CRCs with high PLOD3 expression may be resistant to immune checkpoint therapy.

Pan-Cancer Analysis of PARP1 Alterations as Biomarkers in the Prediction of Immunotherapeutic Effects and the Association of Its Expression Levels and Immunotherapy Signatures

Background

Poly (ADP-ribose) polymerases-1 (PARP1) alterations are associated with PARP1 inhibitor resistance, regulating the function of Treg cells and PDL1 expression in tumor cells, and high PARP1 expression is significantly associated with aggressive behavior and chemotherapeutic resistance in several tumors. However, a comprehensive analysis of the predictive values of PARP1 alteration for immune checkpoint inhibitor (ICI) effectiveness in tumors remains unclear, and the associations between its expression and immunotherapy signatures also needs to be explored further.

Methods

We performed some analyses with the cBioPortal online database (https://www.cbioportal.org), TIMER2.0 (Tumor Immune Estimation Resource 2.0, http://timer.comp-genomics.org/) and TCGA database (https://xenabrowser.net or https://portal.gdc.cancer.gov/). Survival analysis was conducted using Kaplan–Meier method, and the associations between PARP1 transcription levels and immune checkpoint gene expression, the number of neoantigens, tumor mutation burden (TMB) levels, and microsatellite instability (MSI) event are analyzed by spearman correlation analysis and visualization of those mentioned above is performed using R, version 3.6.3 (http://www.r-project.org/).

Results

We found that PARP1 was altered in 1338 (2.9%) out of 45604 patients with diverse tumors, which was associated with markedly higher TMB levels in a variety of tumors (P < 0.01). Impressively, patients with PARP1 alterations in advanced tumors showed better overall survival (OS) in the ICI-treated cohort (P = 0.016). PARP1 altered group was substantially correlated with higher immune infiltrates across most tumors, including CD8+ T cells in colorectal adenocarcinoma (P = 0.0061), endometrial carcinoma (P = 0.0033), stomach cancer (P = 0.033), and cervical cancer (P = 0.026), respectively. The PARP1 altered group showed high expression in transcription (P < 0.001), and higher expression of LAG3, PDCD1, CTLA-4, and TIGIT (P < 0.05). Higher PARP1 expression was present in 27 tumor compared the corresponding normal tissues using the GTEx and TCGA databases and it had a worse OS in several tumors (P < 0.05). Further, high PARP1 expression was significantly associated with six immune cells (B cells, CD4+ T cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells) in most tumors, including colon adenocarcinoma (COAD), head and neck squamous cell carcinoma (HNSC), kidney renal clear cell carcinoma (KIRC), and liver hepatocellular carcinoma (LIHC) (P < 0.05). In particular, CD8+T cell infiltration, was also positively correlated with high PARP1 expression in bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), kidney renal papillary cell carcinoma (KIRP), brain lower grade glioma (LGG), LIHC, pancreatic adenocarcinoma (PAAD), pheochromocytoma and paraganglioma (PCPG), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), testicular germ cell tumors (TGCT), thymoma (THYM), uterine corpus endometrial carcinoma (UCEC), uveal melanoma (UVM) (P < 0.05, no data shown), and PARP1 expression was significantly positively correlated with the transcription levels of some of the 47 immune checkpoint genes, such as CD274, CTLA4, and PDCD1 in several tumors, including PAAD, LIHC, KIRC, HNSC, and BLCA (P < 0.05). A significant positive association between PARP1 expression and the number of immune neoantigen was found within COAD, KIRC, lung adenocarcinoma (LUAD), PAAD and THYM (P < 0.05), and there were also significantly positive correlations between PARP1 expression and TMB in many tumors like adrenocortical carcinoma (ACC), COAD, kidney chromophobe (KICH), LGG, LUAD, READ, skin cutaneous melanoma (SKCM) and stomach adenocarcinoma (STAD) (P < 0.05). In addition, high PARP1 expression was positively associated with microsatellite instability event in COAD, KIRP, BRCA, glioblastoma multiforme (GBM), lung squamous cell carcinoma (LUSC), LGG, READ, UCEC, SKCM and LUAD (P < 0.05).

Conclusions

Our results highlight the significance of PARP1 alterations as pan-cancer predictive biomarkers for ICI treatment, and its expression levels seem to be correlated with the status of immunotherapy-associated signatures, thus may be a promising biomarker for predicting ICI response in several tumors.

A Pan-Cancer Analysis of the Prognostic Value and Expression of Adenylate Cyclase 7 (ADCY7) in Human Tumors

BACKGROUND

The role of adenylate cyclase 7 (ADCY7) in cancer is still unclear. This study analyzed the interrelationship between the expression and immune function of ADCY7.

METHODS

ADCY7 expression in multiple human cancers was analyzed using the databases of Genotype-Tissue Expression Project (GTEx), Cancer Cell Line Encyclopedia (CCLE), and The Cancer Genome Atlas (TCGA). Correlations among ADCY7 gene expression, mismatch repair (MMR) gene expression, and DNA methyltransferase (DNMT) expression were assessed using Spearman correlation analysis. Univariate survival analysis and Kaplan-Meier (KM) curve were used to examine the effect of ADCY7 expression on prognosis. The Tumor Immune Estimation Resource (TIMER) database was used to evaluate the relationship between ADCY7 gene expression and tumor immune invasion or immune checkpoint gene (ICG) expression.

RESULTS

ADCY7 was abnormally expressed in multiple human cancers and was correlated with MMR genes and DNMT expression. Univariate survival analysis and KM curve showed that ADCY7 expression influences the overall survival (OS) of six types of cancer, disease-specific survival (DSS) of eight, and progression-free interval (PFI) of three. The high expression of ADCY7 in OS, DSS, and PFI was strongly associated with poor outcomes in patients with breast cancer and lung squamous cell carcinoma. ADCY7 expression was strongly associated with immune cell infiltration and ICG expression.

CONCLUSION

The results of this study indicated that ADCY7 may be a prognostic biomarker of tumorigenesis. The study may also provide a new perspective on the role of ADCY7 in human cancers.

Pan-Cancer Analysis Identified C1ORF112 as a Potential Biomarker for Multiple Tumor Types

C1ORF112 is an evolutionarily conserved gene across vertebrates. Over the last decade, studies have suggested that C1ORF112 may play a role in tumorigenesis. Using The Cancer Genome Atlas datasets, we explored the role of C1ORF112 across various tumor types in this study. In most tumor types, C1ORF112 expression was increased in tumor tissues compared to corresponding non-tumor tissues. In patients with certain tumor types, higher C1ORF112 expression was correlated with shorter overall survival, disease-free survival, and progression-free survival. Further analyses of C1ORF112 genetic alteration data showed that C1ORF112 amplification and mutations may have an impact on liver hepatocellular carcinoma and uterine corpus endometrial carcinoma prognosis. In cancers including lower grade glioma and adrenocortical carcinoma, C1ORF112 expression was linked to cancer-associated fibroblast infiltration. Gene Ontology analysis showed that C1ORF112 was co-expressed with genes involved in biological processes such as cell cycle and mitotic regulation. The protein interaction network demonstrated that C1ORF112 physically interacted with RAD51, DMC1, and FIGNL1, which have well characterized functions in DNA repair and cell cycle regulation. This pan-cancer study revealed the prognostic value and oncogenic role of C1ORF112 across multiple tumor types.

Identification of a five-immune gene model as an independent prognostic factor in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a common malignant tumor with a poor prognosis. We aimed to identify a new prognostic model of HCC based on differentially expressed (DE) immune genes.

Methods

The DE immune genes were identified based on an analysis of 374 cases of HCC and 50 adjacent non-tumor specimens from the Cancer Genome Atlas (TCGA) database. Univariate Cox analysis, Lasso regression, and multivariate Cox analysis were used to construct the model based on the training group. Survival analysis and the receiver operating characteristic (ROC) curves were used to evaluate model performance. The testing group and the entire group were subsequently used for validation of the model.

Results

A five-immune gene model consisted of HSPA4, ISG20L2, NDRG1, EGF, and IL17D was identified. Based on the model, the overall survival was significantly different between the high-risk and low-risk groups (P = 7.953e-06). The AUCs for the model at 1- and 3-year were 0.849 and 0.74, respectively. The reliability of the model was confirmed using the validation groups. The risk score was identified as an independent prognostic parameter and closely related to the content of immune cells from human HCC specimens.

Conclusion

We identified a five-immune gene model that can be used as an independent prognostic marker for HCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08012-2.