Hypermethylation of the RSK4 promoter associated with BRAF V600E promotes papillary thyroid carcinoma

Analysis of the potential biological value of pyruvate dehydrogenase E1 subunit β in human cancer

BACKGROUND

The pyruvate dehydrogenase E1 subunit β (PDHB) gene which regulates energy metabolism is located in mitochondria. However, few studies have elucidated the role and mechanism of PDHB in different cancers.

AIM

To comprehensive pan-cancer analysis of PDHB was performed based on bioinformatics approaches to explore its tumor diagnostic and prognostic value and tumor immune relevance in cancer. In vitro experiments were performed to examine the biological regulation of PDHB in liver cancer.

METHODS

Pan-cancer data related to PDHB were obtained from the Cancer Genome Atlas (TCGA) database. Analysis of the gene expression profiles of PDHB was based on TCGA and Genotype Tissue Expression Dataset databases. Cox regression analysis and Kaplan-Meier methods were used to assess the correlation between PDHB expression and survival prognosis in cancer patients. The correlation between PDHB and receiver operating characteristic diagnostic curve, clinicopathological staging, somatic mutation, tumor mutation burden (TMB), microsatellite instability (MSI), DNA methylation, and drug susceptibility in pan-cancer was also analyzed. Various algorithms were used to analyze the correlation between PDHB and immune cell infiltration and tumor chemotaxis environment, as well as the co-expression analysis of PDHB and immune checkpoint (ICP) genes. The expression and functional phenotype of PDHB in single tumor cells were studied by single-cell sequencing, and the functional enrichment analysis of PDHB-related genes was performed. The study also validated the level of mRNA or protein expression of PDHB in several cancers. Finally, in vitro experiments verified the regulatory effect of PDHB on the proliferation, migration, and invasion of liver cancer.

RESULTS

PDHB was significantly and differently expressed in most cancers. PDHB was significantly associated with prognosis in patients with a wide range of cancers, including kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, breast invasive carcinoma, and brain lower grade glioma. In some cancers, PDHB expression was clearly associated with gene mutations, clinicopathological stages, and expression of TMB, MSI, and ICP genes. The expression of PDHB was closely related to the infiltration of multiple immune cells in the immune microenvironment and the regulation of tumor chemotaxis environment. In addition, single-cell sequencing results showed that PDHB correlated with different biological phenotypes of multiple cancer single cells. This study further demonstrated that down-regulation of PDHB expression inhibited the proliferation, migration, and invasion functions of hepatoma cells.

CONCLUSION

As a member of pan-cancer, PDHB may be a novel cancer marker with potential value in diagnosing cancer, predicting prognosis, and in targeted therapy.

Prognostic and Diagnostic Values of Semaphorin 5B and Its Correlation With Tumor-Infiltrating Immune Cells in Kidney Renal Clear-Cell Carcinoma

Background: Semaphorin 5B (SEMA5B) has been described to be involved in the development and progression of cancer. However, the potential diagnostic and prognosis roles and its correlation with tumor-infiltrating immune cells in KIRC have not been clearly reported yet. Methods: The mRNA level of SEMA5B was analyzed via the TCGA and GTEx database as well as the CCLE dataset and verified by GSE53757 and GSE40435 datasets. Meanwhile, the protein level of SEMA5B was analyzed by CPTAC and validated by HPA. The diagnostic value of SEMA5B was analyzed according to the TCGA database and validated by GSE53757, GSE46699, and GSE11024 + GSE46699 datasets. Then, the survival analysis was conducted using GEPIA2. R software (v3.6.3) was applied to investigate the relevance between SEMA5B and immune checkpoints and m6A RNA methylation regulator expression. The correlation between SEMA5B and MMRs and DNMT expression and tumor-infiltrating immune cells was explored via TIMER2. Co-expressed genes of SEMA5B were assessed by cBioPortal, and enrichment analysis was conducted by Metascape. The methylation analysis was conducted with MEXPRESS and MethSurv online tools. Gene set enrichment analysis (GSEA) was applied to annotate the biological function of SEMA5B. Results: SEMA5B was significantly upregulated at both the mRNA and protein levels in KIRC. Further analysis demonstrated that the mRNA expression of SEMA5B was significantly correlated with gender, age, T stage, pathologic stage, and histologic grade. High levels of SEMA5B were found to be a favorable prognostic factor and novel diagnostic biomarker for KIRC. SEMA5B expression was shown to be significantly associated with the abundance of immune cells in KIRC. Also, SEMA5B expression was significantly correlated with the abundance of MMR genes, DNMTs, and m6A regulators in KIRC. Enrichment analysis indicated that the co-expressed genes may involve in crosslinking in the extracellular matrix (ECM). GSEA disclosed that SYSTEMIC_LUPUS_ERYTHEMATOSUS and NABA_ECM_REGULATORS were prominently enriched in the SEMA5B low-expression phenotype. Finally, the methylation analysis demonstrated a correlation between hypermethylation of the SEMA5B gene and a poor prognosis in KIRC. Conclusion: Increased SEMA5B expression correlated with immune cell infiltration, which can be served as a favorable prognostic factor and a novel diagnostic biomarker for KIRC.

RAP1GAP Functions as a Tumor Suppressor Gene and Is Regulated by DNA Methylation in Differentiated Thyroid Cancer

Inactivation of tumor suppressor genes, such as RAP1GAP, by hypermethylation of their regulatory region can give rise to thyroid tumors. The aim of this study was to investigate the expression of the RAP1GAP gene and the DNA methylation patterns of its CpG74a, CpG74b, and CpG24 in an Iranian population with differentiated thyroid cancer (DTC). In this study, 160 individuals who underwent thyroidectomy in the Tehran Erfan Hospital between 2018 and 2020 were selected. DNA methylation patterns of selected CpG islands (CpG74a, CpG74b, and CpG24) were determined using methylation-specific PCR. The mRNA expression and protein level of -RAP1GAP were also evaluated. SW1736 and B-CPAP cells were treated with 5-aza-2'-deoxycytidine (5-Aza) to demethylate these regions. The hypermethylation rates of CpG74a and CpG24 in DTC samples were significantly higher than in the control. The mRNA expression and protein level of -RAP1GAP were significantly decreased in the DTC group. In the DTC group, hypermethylation in CpG74a was correlated with decreasing RAP1GAP expression (R2: 0.34; p = 0.043). CpG74a with a specificity of 86.4% has significant prediction power to distinguish between DTC and normal thyroid tissues. Additionally, hypermethylation of CpG74a was significantly associated with higher tumor stages (stage III-IV: 77%; stage I-II: 23%; p = 0.012). Increasing expression of RAP1GAP after demethylation with 15 µM of 5-Aza was observed in both cell lines. These results indicate that DNA hypermethylation in CpG74a can be considered as an epigenetic biomarker in DTC.

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

Background

BRCA1-associated protein (BRAP) is a critical gene that regulates inflammation-related signaling pathway and affects patients’ prognosis in esophageal squamous cell carcinoma (ESCC). However, its roles in different cancers remain largely unknown.

Methods

BRAP expression in human pan-cancer was analyzed via the Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA) database. Pearson correlation analysis was used to evaluate the association between BRAP expression with mismatch repair (MMR) gene mutation and DNA methyltransferase. We evaluated the influence of BRAP on clinical prognosis by univariate survival analysis. Moreover, the correlation between BRAP and tumor immune infiltration was analyzed via the Tumor Immune Evaluation Resource (TIMER) database. Pearson correlation analysis was used to investigate the correlation between BRAP expression and immune checkpoint genes expression.

Results

BRAP is abnormally overexpressed and significantly correlated with MMR gene mutation level and DNA methyltransferase expression in human pan-cancer. Univariate survival analysis showed that BRAP was significant with patients’ overall survival (OS) in six cancer types, disease-free interval (DFI) in three cancer types, and progression-free interval (PFI) in two cancer types. Remarkably, increased BRAP expression was strongly correlated with patients’ poor prognosis in liver hepatocellular carcinoma (LIHC), whether OS (P < 0.0001, hazard ratio (HR) = 1.1), DFI (P = 0.00099, HR = 1.06), or PFI (P = 0.00025, HR = 1.07). Moreover, a positive relationship was found between BRAP expression and immune infiltrating cells including B cell, CD4 + T cell, CD8 + T cell, dendritic cell, macrophage cell, and neutrophil cell in colon adenocarcinoma (COAD), kidney renal clear cell carcinoma (KIRC), and LIHC. Additionally, BRAP expression showed strong correlations with immune checkpoint genes in LIHC.

Conclusion

BRAP expression is increased in human pan-cancer samples compared with normal tissues. Overexpression of BRAP is correlated with poor prognosis and immune infiltration in multiple cancers, especially in LIHC. These findings suggest that BRAP may be used as a potential molecular biomarker for determining prognosis and immune infiltration in LIHC.

NFE2L2 Is a Potential Prognostic Biomarker and Is Correlated with Immune Infiltration in Brain Lower Grade Glioma: A Pan-Cancer Analysis

Nuclear factor, erythroid 2 like 2 (NFE2L2, NRF2) is a transcription factor that regulates various antioxidant enzymes. It plays a vital physiological role in regulating oxidative stress and inflammatory response. However, the roles of NFE2L2 in human cancers are still unclear. Our study is aimed at analyzing the prognostic value of NFE2L2 in pan-cancer and at revealing the relationship between NFE2L2 expression and tumor immunity. The present study revealed that NFE2L2 was abnormally expressed and significantly correlated with mismatch repair (MMR) gene mutation levels and DNA methyltransferase expression in human pan-cancer. In particular, pan-cancer survival analysis indicated that NFE2L2 expression was associated with adverse outcomes-overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI)-in adrenocortical carcinoma (ACC), brain lower grade glioma (LGG), and pancreatic adenocarcinoma (PAAD) patients. A positive relationship was also found between NFE2L2 expression and immune infiltration, including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells, especially in breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), kidney renal clear cell carcinoma (KIRC), LGG, liver hepatocellular carcinoma (LIHC), and prostate adenocarcinoma (PRAD). Additionally, NFE2L2 expression was positively correlated with the immune score and the expression of immune checkpoint markers in LGG. In conclusion, these results indicate that transcription factor NFE2L2 is a potential prognostic biomarker and is correlated with immune infiltration in LGG.