YTHDF1 Is a Potential Pan-Cancer Biomarker for Prognosis and Immunotherapy

Bioinformatics Analysis Revealing the Correlation between NF-κB Signaling Pathway and Immune Infiltration in Gastric Cancer

Although the emerging of immunotherapy conferred a new landscape of gastric cancer (GC) treatment, its response rate was of significant individual differences. Insight into GC immune microenviroment may contribute to breaking the dilemma. To this end, the enrichment score of NF-κB signaling pathway was calculated in each GC sample from The Cancer Genome Atlas (TCGA) via ssGSEA algorithm, and its association with immune infiltration was estimated. Based on NF-κB-related genes, a risk score was established and its involvement in immune infiltration, tumor mutational burden (TMB), and N6-methyladenosine (M6A) modification was analyzed in GC. The results showed that NF-κB signaling pathway promoted the infiltration of immune cells in GC. In addition, GC samples were divided into low- and high-risk groups according to a seven-gene (CARD11, CCL21, GADD45B, LBP, RELB, TRAF1, and VCAM1) risk score. Although the high-risk group displayed high immune infiltration and high expression of M6A regulatory genes, it remains in an immunosuppressive microenviroment and whereby suffers a poorer outcome. Of note, most of hub genes were related to immune infiltration and could serve as an independent prognostic biomarker. Conclusively, our study emphasized the crucial role of NF-κB signaling pathway in GC immune microenviroment and provided several candidate genes that may participate in immune infiltration.

EXTL3 could serve as a potential biomarker of prognosis and immunotherapy for prostate cancer and its potential mechanisms

Background

Exostosin like glycosyltransferase 3 (EXTL3) had been reported to be associated with immune deficiency and play prognostic roles in various cancers. However, little is known about the associations between EXTL3 and prostate cancer (PCa). Hence, this article was designed to clarify their associations.

Methods

All original data were downloaded from The Cancer Genome Atlas (TCGA) database. Gene set enrichment analysis (GSEA) and CellMiner database was utilized, respectively, to identify EXTL3-related signaling pathways and drugs. We explored the relationships between EXTL3 expression and immunity to further evaluate the involvement of EXTL3 in response to immunotherapies. LncRNA/RBP/EXTL3 mRNA networks were also identified for its potential mechanism.

Results

Compared with normal prostate samples, EXTL3 was poorly expressed in PCa samples not only in mRNA expression levels, but also in protein expression levels, with worse overall survival (P < 0.05) and this gene could be an independent prognostic biomarker for PCa (both P < 0.05). EXTL3 was revealed to be markedly linked with seven signaling pathways in PCa by GSEA, including calcium, chemokine, ERBB, JAK STAT, MAPK, WNT, oxidative phosphorylation pathways. EXTL3 expression was also revealed to be significantly associated with MSI, immune cells, immune checkpoint molecules, tumor microenvironment and immune cells infiltration. We further predicted immune responses of EXTL3 gene to immunotherapies by TIDE database and the IMvigor210 cohort. A total of six LncRNA/RBP/EXTL3 mRNA networks were eventually identified for its potential mechanisms.

Conclusions

EXTL3 could serve as a potential biomarker of prognosis and immunotherapy for PCa and six LncRNA/RBP/EXTL3 mRNA networks were also identified for its potential mechanisms.

GRAMD1A Is a Biomarker of Kidney Renal Clear Cell Carcinoma and Is Associated with Immune Infiltration in the Tumour Microenvironment

Background

GRAM structural domain-containing protein 1A (GRAMD1A) is upregulated in a variety of human cancer tissues and is closely associated with tumourigenesis and progression.

Methods

Patient RNA-sequencing data and clinicopathological information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. The expression of GRAMD1A in kidney cancer cell lines and KIRC patients was analysed by quantitative polymerase chain reaction (qPCR). Receiver Operator Characteristic (ROC) curves, nomograms, Kaplan-Meier analysis, forest plots, and COX analysis were used to assess the diagnostic and prognostic value of GRAMD1A in KIRC, and gene set enrichment analysis (GSEA) was used to explore its potential signalling pathways. In addition, the Sangerbox website, Kaplan-Meier plotter database, and TISIDB and TIMER databases were used to further analyse the correlation of GRAMD1A with microsatellite instability (MSI), tumour mutational burden (TMB), immune checkpoint genes, and tumour-infiltrating lymphocytes (TILs).

Results

GRAMD1A was significantly highly expressed in KIRC and associated with shorter overall survival and relapse-free survival (P < 0.05). The AUC value of the ROC curve to identify KIRC and normal renal tissues was 0.942. Forest plot and COX analysis visualized that GRAMD1A could be an independent prognostic factor in KIRC patients (P < 0.01), and nomograms to determine the overall survival (OS) of KIRC patients also showed good efficacy (C-index: 0.776). Moreover, Spearman correlation analysis showed a positive correlation between GRAMD1A and MSI, TMB (P < 0.01). On the other hand, GRAMD1A was also found to be closely associated with immune checkpoint genes. Meanwhile, patients with KIRC with high GRAMD1A expression had a relatively low hazard ratio (HR) of death when B lymphocytes, natural killer T cells, CD4+ T lymphocytes, CD8+ T lymphocytes, and macrophages were enriched in the tumour microenvironment (TME), and a greater HR of death when regulatory T lymphocytes with tumour-specific immunosuppressive effects were significantly enriched. Last, GSEA shows that GRAMD1A is closely associated with the regulation of energy metabolism in KIRC.

Conclusions

GRAMD1A is a promising diagnostic and prognostic biomarker for patients with KIRC, and its biological function correlates to some extent with immune infiltration in TME.

Liquid-liquid phase separation in tumor biology

Liquid–liquid phase separation (LLPS) is a novel principle for explaining the precise spatial and temporal regulation in living cells. LLPS compartmentalizes proteins and nucleic acids into micron-scale, liquid-like, membraneless bodies with specific functions, which were recently termed biomolecular condensates. Biomolecular condensates are executors underlying the intracellular spatiotemporal coordination of various biological activities, including chromatin organization, genomic stability, DNA damage response and repair, transcription, and signal transduction. Dysregulation of these cellular processes is a key event in the initiation and/or evolution of cancer, and emerging evidence has linked the formation and regulation of LLPS to malignant transformations in tumor biology. In this review, we comprehensively summarize the detailed mechanisms of biomolecular condensate formation and biophysical function and review the recent major advances toward elucidating the multiple mechanisms involved in cancer cell pathology driven by aberrant LLPS. In addition, we discuss the therapeutic perspectives of LLPS in cancer research and the most recently developed drug candidates targeting LLPS modulation that can be used to combat tumorigenesis.

The Potential Value of m6A RNA Methylation in the Development of Cancers Focus on Malignant Glioma

N6-methyladenosine (m6A) RNA methylation is an epigenetic modification that has emerged in the last few years and has received increasing attention as the most abundant internal RNA modification in eukaryotic cells. m6A modifications affect multiple aspects of RNA metabolism, and m6A methylation has been shown to play a critical role in the progression of multiple cancers through a variety of mechanisms. This review summarizes the mechanisms by which m6A RNA methylation induced peripheral cancer cell progression and its potential role in the infiltration of immune cell of the glioblastoma microenvironment and novel immunotherapy. Assessing the pattern of m6A modification in glioblastoma will contribute to improving our understanding of microenvironmental infiltration and novel immunotherapies, and help in developing immunotherapeutic strategies.

AFF3 is a novel prognostic biomarker and a potential target for immunotherapy in gastric cancer

Background

Gastric cancer (GC) is one of the most common cancers worldwide with a poor prognosis. The tumor microenvironment (TME) serves a pivotal role in affecting the prognosis and efficacy of immunotherapy. Given the poor prognosis of GC patients and the limitation of immunotherapy, we urged to identify new prognostic and immunotherapeutic biomarkers.

Methods

The transcriptome data were downloaded from the TCGA, GEO, and GEPIA databases, and performed differential analysis of AFF3 in tumor samples and normal samples. The UALCAN, Kaplan–Meier plotter and GEPIA databases were employed to assess the correlation of AFF3 with clinicopathological characteristics and prognosis. The potential mechanism of AFF3 was explored by the GO and KEGG enrichment. The potential role of AFF3 on tumor‐infiltrating immune cells (TIICs) was explored by TIMER2.0 and TISIDB. TIMER2.0 and SangerBox3.0 databases were, respectively, used to determine the correlation of AFF3 with immune checkpoint (ICs), tumor mutational burden (TMB), and microsatellite instability (MSI) in GC.

Results

We found significant downregulation of AFF3 in GC tissues as compared with normal tissues. However, GC patients having a higher expression of AFF3 were found to have worse clinicopathological characteristics and prognosis. Moreover, the GO enrichment analysis illustrated that AFF3 might regulate the immune cells in the TME. In addition, the AFF3 was positively correlated with TIICs, ICs, TMB, and MSI.

Conclusion

Here, we conclude that AFF3 may be a promising potential marker for the diagnosis and prognosis of GC patients, and may influence response to ICIs by affecting TIICs and ICs expression in the TME.

Pan-Cancer Analysis Shows That ALKBH5 Is a Potential Prognostic and Immunotherapeutic Biomarker for Multiple Cancer Types Including Gliomas

Background

AlkB homolog 5 (ALKBH5) is a N⁶-methyladenosine (m⁶A) demethylase associated with the development, growth, and progression of multiple cancer types. However, the biological role of ALKBH5 has not been investigated in pan-cancer datasets. Therefore, in this study, comprehensive bioinformatics analysis of pan-cancer datasets was performed to determine the mechanisms through which ALKBH5 regulates tumorigenesis.

Methods

Online websites and databases such as NCBI, UCSC, CCLE, HPA, TIMER2, GEPIA2, cBioPortal, UALCAN, STRING, SangerBox, ImmuCellAl, xCell, and GenePattern were used to extract data of ALKBH5 in multiple cancers. The pan-cancer patient datasets were analyzed to determine the relationship between ALKBH5 expression, genetic alterations, methylation status, and tumor immunity. Targetscan, miRWalk, miRDB, miRabel, LncBase databases and Cytoscape tool were used to identify microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) that regulate expression of ALKBH5 and construct the lncRNA-miRNA-ALKBH5 network. In vitro CCK-8, wound healing, Transwell and M2 macrophage infiltration assays as well as in vivo xenograft animal experiments were performed to determine the biological functions of ALKBH5 in glioma cells.

Results

The pan-cancer analysis showed that ALKBH5 was upregulated in several solid tumors. ALKBH5 expression significantly correlated with the prognosis of cancer patients. Genetic alterations including duplications and deep mutations of the ALKBH5 gene were identified in several cancer types. Alterations in the ALKBH5 gene correlated with tumor prognosis. GO and KEGG enrichment analyses showed that ALKBH5-related genes were enriched in the inflammatory, metabolic, and immune signaling pathways in glioma. ALKBH5 expression correlated with the expression of immune checkpoint (ICP) genes, and influenced sensitivity to immunotherapy. We constructed a lncRNA-miRNA network that regulates ALKBH5 expression in tumor development and progression. In vitro and in vivo experiments showed that ALKBH5 promoted proliferation, migration, and invasion of glioma cells and recruited the M2 macrophage to glioma cells.

Conclusions

ALKBH5 was overexpressed in multiple cancer types and promoted the development and progression of cancers through several mechanisms including regulation of the tumor-infiltration of immune cells. Our study shows that ALKBH5 is a promising prognostic and immunotherapeutic biomarker in some malignant tumors.

Targeting the RNA m6A modification for cancer immunotherapy

N⁶-methyladenosine (m⁶A) is the most abundant epigenetic modification of RNA, and its dysregulation drives aberrant transcription and translation programs that promote cancer occurrence and progression. Although defective gene regulation resulting from m⁶A often affects oncogenic and tumor-suppressing networks, m⁶A can also modulate tumor immunogenicity and immune cells involved in anti-tumor responses. Understanding this counterintuitive concept can aid the design of new drugs that target m⁶A to potentially improve the outcomes of cancer immunotherapies. Here, we provide an up-to-date and comprehensive overview of how m⁶A modifications intrinsically affect immune cells and how alterations in tumor cell m⁶A modifications extrinsically affect immune cell responses in the tumor microenvironment (TME). We also review strategies for modulating endogenous anti-tumor immunity and discuss the challenge of reshaping the TME. Strategies include: combining specific and efficient inhibitors against m⁶A regulators with immune checkpoint blockers; generating an effective programmable m⁶A gene-editing system that enables efficient manipulation of individual m⁶A sites; establishing an effective m⁶A modification system to enhance anti-tumor immune responses in T cells or natural killer cells; and using nanoparticles that specifically target tumor-associated macrophages (TAMs) to deliver messenger RNA or small interfering RNA of m⁶A-related molecules that repolarize TAMs, enabling them to remodel the TME. The goal of this review is to help the field understand how m⁶A modifications intrinsically and extrinsically shape immune responses in the TME so that better cancer immunotherapy can be designed and developed.

Comprehensive Analyses of the Immunological and Prognostic Roles of an IQGAP3AR/let-7c-5p/IQGAP3 Axis in Different Types of Human Cancer

IQ motif containing GTPase-activating protein 3 (IQGAP3) is a member of the Rho family of guanosine-5′-triphosphatases (GTPases). IQGAP3 plays a crucial part in the development and progression of several types of cancer. However, the prognostic, upstream-regulatory, and immunological roles of IQGAP3 in human cancer types are not known. We found that IQGAP3 expression was increased in different types of human cancer. The high expression of IQGAP3 was correlated with tumor stage, lymph node metastasis, and a poor prognosis in diverse types of human cancer. The DNA methylation of IQGAP3 was highly and negatively correlated with IQGAP3 expression in diverse cancer types. High DNA methylation in IQGAP3 was correlated with better overall survival in human cancer types. High mRNA expression of IQGAP3 was associated with tumor mutational burden, microsatellite instability, immune cell infiltration, and immune modulators. Analyses of signaling pathway enrichment showed that IQGAP3 was involved in the cell cycle. IQGAP3 expression was associated with sensitivity to a wide array of drugs in cancer cells lines. We revealed that polypyrimidine tract–binding protein 1 (PTBP1) and an IQGAP3-associated lncRNA (IQGAP3AR)/let-7c-5p axis were potential regulations for IQGAP3 expression. We provided the first evidence to show that an IQGAP3AR/let-7c-5p/IQGAP3 axis has indispensable roles in the progression and immune response in different types of human cancer.

FOXM1/lncRNA TYMSOS/miR-214-3p–Mediated High Expression of NCAPG Correlates With Poor Prognosis and Cell Proliferation in Non–Small Cell Lung Carcinoma

Lung cancer is the most common cancer with high mortality. Increasing evidence has demonstrated that nonstructural maintenance of chromosomes condensin I complex subunit G (NCAPG) plays a crucial role in the progression of human cancers. However, the biological function and underlying mechanism of NCAPG in non–small cell lung cancer (NSCLC) are still unclear. Here, we utilized diverse public databases to analyze the expression of NCAPG in pan-cancer. We found that NCAPG was highly expressed in various human cancers, especially in NSCLC. NCAPG expression was significantly positively correlated with poor clinical-pathological features, poor prognosis, tumor mutational burden, DNA microsatellite instability, and immune cell infiltration in NSCLC. In addition, our results showed that depletion of NCAPG significantly inhibited NSCLC cell proliferation, migration, and self-renewal abilities, yet these could be reversed by adding microRNA (miRNA)-214-3p. Knockdown of long noncoding RNA (lncRNA) thymidylate synthetase opposite strand (TYMSOS) also inhibits the NSCLC cell proliferation, migration, and self-renewal abilities. In summary, our findings demonstrated that the crucial roles of the FOXM1/lncRNA-TYMSOS/miRNA-214-3p/NCAPG axis in NSCLC may shed light on how NCAPG may act as a therapeutic target for NSCLC.

Pan-cancer analysis of N4-acetylcytidine adaptor THUMPD1 as a predictor for prognosis and immunotherapy

Background: THUMPD1 is a specific RNA adaptor that assists acetylation of mRNA and production of N4-acetylcytidine (ac4C). However, it remains unclear whether THUMPD1 plays a part in tumorigenesis and therapeutic efficacy. Here, we analyzed the expression profiles and prognostic value of THUMPD1 in pan-cancer and gained insights into the correlation between THUMPD1 expression level and immunotherapy efficacy.

Methods: Gene expression pattern and its correlation with prognosis, immune cell infiltration in pan-cancer were obtained from Genotype-Tissue Expression (GTEx), Cancer Cell Line Encyclopedia (CCLE) and The Cancer Genome Atlas (TCGA) databases, with Kaplan–Meier method and Spearman correlation analysis used. Western blotting and immunofluorescence on clinical samples were performed to validate our database-derived results. Correlation between THUMPD1 expression level and immunotherapy responses was also explored, based on clinical cohorts receiving programmed cell death protein 1 ligand (PD-L1) antibody therapy. Finally, gene set enrichment analysis (GSEA) was performed to show the possible tumorigenic mechanism.

Results: THUMPD1 was highly expressed in most cancer types, and this elevated expression indicated poor or improved prognosis for different cancers. In kidney renal clear cell carcinoma (KIRC) and rectum adenocarcinoma (READ), patients with higher THUMPD1 expression exhibited a better prognosis, while liver hepatocellular carcinoma (LIHC) patients had worse prognosis. Besides, THUMPD1 was significantly associated with immune cell infiltration, tumor mutational burden (TMB), microsatellite instability (MSI), immune checkpoints and neoantigen in many cancer types. Further, more clinical advantages and therapeutic responses were observed in patients with high THUMPD1 expression.

Conclusions: THUMPD1 may serve as a novel predictor to evaluate cancer prognosis and immune therapy efficacy in diverse cancer types.

Comprehensive Analysis of Regulatory Factors and Immune-Associated Patterns to Decipher Common and BRCA1/2 Mutation-Type-Specific Critical Regulation in Breast Cancer

Background: BRCA1/2 mutations are closely related to high lifetime risk of breast cancer (BC). The objective of this study was to identify the genes, regulators, and immune-associated patterns underlying disease pathology in BC with BRCA1/2 somatic mutations and their associations with clinical traits. Methods: RNA sequencing data and clinical information from The Cancer Genome Atlas (TCGA; N = 36 BRCA1-mutant BC; N = 49 BRCA2-mutant BC; and N = 117 BRCA1/2-wild-type BC samples) were used for discovery, which included consensus network analysis, function enrichment, and analysis of hub genes; other TCGA data (N = 117 triple-negative BC) and two Gene Expression Omnibus database expression profiles were used as validation cohorts. Results: Consensus network analysis helped to identify specific co-expressed modules that showed positive correlations with tumor stage, number of positive lymph nodes, and margin status in BRCA1/2-mutant BC but lacking correlations in BRCA1/2-wild-type BC. Functional enrichment suggested potential mechanisms in BRCA1/2 carriers that could regulate the cell cycle, immune response, cellular metabolic processes, and cell migration, via enriched pathways including p53 and JAK-STAT signaling. Consensus network analysis identified the specific and common carcinogenic mechanisms involving BRCA mutations. Regulators cross-linking these modules include E2F or IRF transcription factor family, associated with cell cycle or immune response regulation module, respectively. Eight hub genes, including ISG15, BUB1, and TTK, were upregulated in several BRCA1/2-mutant BC datasets and showed prognostic value in BC. Furthermore, their genetic expression was related to higher levels of immune infiltration in BRCA1/2-mutant BC, which manifested as recruitment of T helper cells (Th1 cells), follicular helper T cells, and regulatory T cells, and T cell exhaustion. Moreover, important indicators for evaluation of BC immunotherapy, tumor mutational burden and neoantigen load also positively correlated with expression of some hub genes. Conclusion: We constructed a BRCA1/2 mutation-type-specific co-expressed gene network with related transcription factors and immune-associated patterns that could regulate and influence tumor metastasis and immune microenvironment, providing novel insights into the pathological process of this disease and the corresponding BRCA mutations.

Pan-cancer analysis combined with experiments predicts CTHRC1 as a therapeutic target for human cancers

BACKGROUND

The function of collagen triple helix repeat containing 1 (CTHRC1) as an oncogene has been reported in a growing number of publications. Bioinformatics methods represent a beneficial approach to examine the mechanism and function of the CTHRC1 gene in the disease process of cancers from a pan-cancer perspective.

METHODS

In this study, using the online databases UCSC, NCBI, HPA, TIMER2, Oncomine, GEPIA, UALCAN, cBioPortal, COSMIC, MEXPRESS, STRING, CCLE, LinkedOmics, GTEx, TCGA, CGGA, and SangerBox, we focused on the relationship between CTHRC1 and tumorigenesis, progression, methylation, immunity, and prognosis. qPCR was used to detect CTHRC1 expression in glioma tissues and cell lines.

RESULTS

The pan-cancer analysis showed that CTHRC1 was overexpressed in most tumors, and a significant correlation was observed between CTHRC1 expression and the prognosis of patients with cancer. CTHRC1 genetic alterations occur in diverse tumors and are associated with tumor progression. Levels of CTHRC1 promoter methylation were decreased in most cancer tissues compared with normal tissues. In addition, CTHRC1 coordinated the activity of ICP genes through diverse signal transduction pathways, was also associated with immune cell infiltration and the tumor microenvironment, and potentially represented a promising immunotherapy target. We identified CTHRC1-related genes across cancers using the GEPIA2 tool. The single-gene GO analysis of CTHRC1 across cancers showed that it was involved in some signaling pathways and biological processes, such as the Wnt signaling pathway, cell migration, and positive regulation of protein binding. The expression and function of CTHRC1 were also further verified in glioma tissues and cell lines.

CONCLUSIONS

CTHRC1 is overexpressed in various cancer types and functions as an important oncogene that may promote tumorigenesis and development through different mechanisms. CTHRC1 may represent an important therapeutic target for human cancers.

Elevated Ras related GTP binding B (RRAGB) expression predicts poor overall survival and constructs a prognostic nomogram for colon adenocarcinoma

Currently, no articles have explored the roles of RRAGB gene in the occurrence and development of cancer. By means of The Cancer Genome Atlas (TCGA) data mining, we found that this gene might be a novel prognostic predictor for colon adenocarcinoma (COAD). Hence, this article was carried out to explore its roles in COAD and associations with immunity. RRAGB single-gene expression matrix and corresponding clinical information were extracted from TCGA database. Univariate/multivariate cox regression analyses and gene set enrichment analysis (GSEA) were utilized to identify independent prognostic factors and RRAGB related pathways, respectively. Relationships between RRAGB and immunity were also analyzed. Boxplot and K-M survival analysis indicated that RRAGB was not only differently expressed in COAD (P < 0.05), but also significantly associated with overall survival (OS; P < 0.05). Univariate and multivariate Cox hazard regression analyses indicated that RRAGB could serve as an independent prognostic factor for COAD (both P < 0.05). GSEA identified five signaling pathways significantly enriched in the high-RRAGB expression phenotype. Moreover, a RRAGB-based nomogram was successfully constructed and displayed a satisfactory performance. In addition, RRAGB expression was found to be significantly associated with microsatellite instability (MSI), tumor mutational burden (TMB) and immunity. Our results revealed that RRAGB could be a prognostic biomarker for COAD in terms of OS and markedly related to MSI, TMB, and immunity. We also constructed an RRAGB-based nomogram with a satisfactory performance. Further researches should be carried out to validate our findings.