A Prognostic Signature for Lower Grade Gliomas Based on Expression of Long Non-Coding RNAs

The aggrephagy-related gene TUBA1B influences clinical outcomes in glioma patients by regulating the cell cycle

Background

Gliomas are common primary malignant brain tumors, with glioblastoma (GBM) being the most aggressive subtype. GBM is characterized by high recurrence rates and treatment resistance, leading to poor patient outcomes. Current prognostic models have limited predictive power, underscoring the need to elucidate underlying mechanisms and identify novel biomarkers to improve therapeutic strategies and prognostic models.

Methods

Gene expression and clinical data for GBM and LGG were obtained from the TCGA and CGGA database, while single-cell sequencing data from GSE167960 were selected from the GEO database. Molecular characteristics of gliomas were revealed through normalization, consensus clustering analysis, immune scoring, cell infiltration analysis, and pathway analysis. TUBA1B, identified as a key gene through machine learning, was incorporated into a nomogram model using multivariate Cox regression. Its functions were validated through qRT-PCR, in vitro functional assays, and mouse xenograft models. All data analyses and statistics were performed using R software.

Results

Consensus clustering of the TCGA glioma dataset identified two aggrephagy subtypes (C1 and C2), with C2 showing worse survival outcomes and higher immune infiltration. TUBA1B was identified as an independent prognostic marker, with high expression associated with upregulated cell cycle pathways and alterations in the immune microenvironment. TUBA1B was shown to influence glioma cell proliferation, migration, invasion, and autophagy, impacting tumor progression and treatment response through intercellular communication and metabolic pathways.

Conclusion

The study demonstrates that high TUBA1B expression is closely associated with glioma malignancy and poor prognosis, making it a potential therapeutic target.

DSN1 may predict poor prognosis of lower-grade glioma patients and be a potential target for immunotherapy

DSN1 has been previously found to be positively correlated with various cancers. However, the effect of DSN1 or its methylation on the prognosis, molecular characteristics, and immune cell infiltration of low-grade glioma (LGG) has not yet been studied. We obtained 1046 LGG samples from the The Cancer Genome Atlas, The Chinese Glioma Genome Atlas (CGGA) microarray, and CGGA RNA-Seq databases. Bioinformatic methods (gene set enrichment analysis (GSEA), chi-square test, multivariate), and laboratory validation were used to investigate DSN1 in LGG. The expression levels of DSN1 mRNA and protein in LGG were substantially higher than those in normal brain tissue, and their expression was negatively regulated by methylation. The survival time of patients with low expression of DSN1 and cg12601032 hypermethylation was considerably prolonged. DSN1 was a risk factor, and of good diagnostic and prognostic value for LGG. Importantly, the expression of DSN1 is related to many types of tumor-infiltrating immune cells and has a positive correlation with PDL1. DSN1 promoted the activation of multiple cancer-related pathways, such as the cell cycle. Additionally, knockdown of DSN1 substantially inhibited the proliferation and invasion of LGG cells. To the best of our knowledge, this study is the first comprehensive analysis of the mechanism of DSN1 leading to poor prognosis of LGG, which provides a new perspective for revealing the pathogenesis of LGG. DSN1 or its methylation has diagnostic value for the prognosis of glioma, and may become a new biological target of anti-tumor immunotherapy.

Machine learning-based identification of histone deacetylase-associated prognostic factors and prognostic modeling for low-grade glioma

BACKGROUND

Low-grade glioma (LGG) is a slow-growing but invasive tumor that affects brain function. Histone deacetylases (HDACs) play a critical role in gene regulation and tumor progression. This study aims to develop a prognostic model based on HDAC-related genes to aid in risk stratification and predict therapeutic responses.

METHODS

Expression data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) were analyzed to identify an optimal HDAC-related risk signature from 73 genes using 10 machine learning algorithms. Patients were stratified into high- and low-risk groups based on the median risk score. Prognostic accuracy was evaluated using Kaplan-Meier survival analysis and receiver operating characteristic (ROC) curves. Functional enrichment analyses, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA), were performed to explore pathways linked to the gene signature. Immune infiltration and tumor microenvironment characteristics were assessed using Single Sample Gene Set Enrichment Analysis (ssGSEA) and ESTIMATE algorithm. SubMap was applied to predict responsiveness to immune checkpoint inhibitors, and chemotherapeutic sensitivity was analyzed via the Genomics of Drug Sensitivity in Cancer (GDSC) database.

RESULTS

A prognostic model consisting of four HDAC-related genes-SP140, BAZ1B, SP100, and SIRT1-was identified. This signature displayed strong prognostic accuracy, achieving a C-index of 0.945. Individuals with LGG were systematically divided into high-risk and low-risk cohorts based on the median risk value, enabling more precise risk stratification. The survival prognosis was significantly worse in the high-risk cohort compared to the low-risk group, highlighting distinct survival trajectories. Notably, the two cohorts exhibited marked shifts in immune checkpoint gene transcriptional profiles and immune cell infiltration maps, underscoring fundamental biological differences that contribute to these differing prognoses.

CONCLUSION

We developed an HDAC-related four-gene prognostic model that correlates with survival, immune landscape, and therapeutic response in LGG patients. This model may guide personalized treatment strategies and improve prognostic accuracy, warranting further validation in clinical settings.

A comprehensive analysis of Transcribed Ultra Conserved Regions uncovers important regulatory functions of novel non-coding transcripts in gliomas

Transcribed Ultra-Conserved Regions (TUCRs) represent a severely understudied class of putative non-coding RNAs (ncRNAs) that are 100% conserved across multiple species. We performed the first-ever analysis of TUCRs in glioblastoma (GBM) and low-grade gliomas (LGG). We leveraged large human datasets to identify the genomic locations, chromatin accessibility, transcription, differential expression, correlation with survival, and predicted functions of all 481 TUCRs, and identified TUCRs that are relevant to glioma biology. Of these, we investigated the expression, function, and mechanism of action of the most highly upregulated intergenic TUCR, uc.110, identifying it as a new tumor enhancer. Uc.110 was highly overexpressed in GBM and LGG, where it promoted malignancy and tumor growth. Uc.110 activated the WNT pathway by upregulating the expression of membrane frizzled-related protein (MFRP), by sponging the tumor suppressor microRNA miR-544. This pioneering study shows important roles for TUCRs in gliomas and provides an extensive database and novel methods for future TUCR research.

Radiomics Analysis of Quantitative Maps from Synthetic MRI for Predicting Grades and Molecular Subtypes of Diffuse Gliomas

PURPOSE

To investigate the feasibility of using radiomics analysis of quantitative maps from synthetic MRI to preoperatively predict diffuse glioma grades, isocitrate dehydrogenase (IDH) subtypes, and 1p/19q codeletion status.

METHODS

Data from 124 patients with diffuse glioma were used for analysis (n = 87 for training, n = 37 for testing). Quantitative T1, T2, and proton density (PD) maps were obtained using synthetic MRI. Enhancing tumour (ET), non-enhancing tumour and necrosis (NET), and peritumoral edema (PE) regions were segmented followed by manual fine-tuning. Features were extracted using PyRadiomics and then selected using Levene/T, BorutaShap and maximum relevance minimum redundancy algorithms. A support vector machine was adopted for classification. Receiver operating characteristic curve analysis and integrated discrimination improvement analysis were implemented to compare the performance of different radiomics models.

RESULTS

Radiomics models constructed using features from multiple tumour subregions (ET + NET + PE) in the combined maps (T1 + T2 + PD) achieved the highest AUC in all three prediction tasks, among which the AUC for differentiating lower-grade and high-grade diffuse gliomas, predicting IDH mutation status and predicting 1p/19q codeletion status were 0.92, 0.95 and 0.86 respectively. Compared with those constructed on individual T1, T2, and PD maps, the discriminant ability of radiomics models constructed on the combined maps separately increased by 11, 17 and 10% in predicting glioma grades, 35, 52 and 19% in predicting IDH mutation status, and 16, 15 and 14% in predicting 1p/19q codeletion status (p < 0.05).

CONCLUSION

Radiomics analysis of quantitative maps from synthetic MRI provides a new quantitative imaging tool for the preoperative prediction of grades and molecular subtypes in diffuse gliomas.

Prognostic value of four immune-related genes in lower-grade gliomas: a biomarker discovery study

Introduction

The tumor microenvironment and IRGs are highly correlated with tumor occurrence, progression, and prognosis. However, their roles in grade II and III gliomas, termed LGGs in this study, remain to be fully elucidated. Our research aims to develop immune-related features for risk stratification and prognosis prediction in LGG.

Methods

Using the ssGSEA method, we assessed the immune characteristics of the LGG population. We conducted differential analysis using LGG samples from the TCGA database and normal samples from GTEx, identifying 412 differentially expressed immune-related genes (DEIRGs). Subsequently, we utilized univariate Cox, LASSO, and multivariate Cox regression analyses to establish both a gene predictive model and a nomogram predictive model.

Results

Here, we found that the ESTIMATE score, immune score and stromal score of high-immunity, high-grade and isocitrate dehydrogenase (IDH) wild-type glioma were higher than those of the corresponding group, and the tumor purity was lower. Higher ESTIMATE scores, stromal scores and immune scores indicated a poor prognosis in patients with LGG. Our four-gene prognostic model demonstrated superior accuracy compared to other molecular features. Validation using the CGGA as a testing set and the combined TCGA and CGGA cohort confirmed its robust prognostic value. Additionally, a nomogram integrating the prognostic model and clinical variables showed enhanced predictive capability.

Discussion

Our study highlights the prognostic significance of the identified four DEIRGs (KLRC3, MR1, PDIA2, and RFXAP) in LGG patients. The predictive model and nomogram developed herein offer valuable tools for personalized treatment strategies in LGG. Future research should focus on further validating these findings and exploring the functional roles of these DEIRGs within the LGG tumor microenvironment.

The Potential Significance of the EMILIN3 Gene in Augmenting the Aggressiveness of Low-Grade Gliomas is Noteworthy

Purpose

Low-grade gliomas (LGG) are common brain tumors with high mortality rates. Cancer cell invasion is a significant factor in tumor metastasis. Novel biomarkers are urgently needed to predict LGG prognosis effectively.

Methods

The data for LGG were obtained from the Bioinformatics database. A consensus clustering analysis was performed to identify molecular subtypes linked with invasion in LGG. Differential expression analysis was performed to identify differentially expressed genes (DEGs) between the identified clusters. Enrichment analyses were then conducted to explore the function for DEGs. Prognostic signatures were placed, and their predictive power was assessed. Furthermore, the invasion-related prognostic signature was validated using the CGGA dataset. Subsequently, clinical specimens were procured in order to validate the expression levels of the distinct genes examined in this research, and to further explore the impact of these genes on the glioma cell line LN229 and HS-683.

Results

Two invasion-related molecular subtypes of LGG were identified, and we sifted 163 DEGs between them. The enrichment analyses indicated that DEGs are mainly related to pattern specification process. Subsequently, 10 signature genes (IGF2BP2, SRY, CHI3L1, IGF2BP3, MEOX2, ABCC3, HOXC4, OTP, METTL7B, and EMILIN3) were sifted out to construct a risk model. Besides, the survival (OS) in the high-risk group was lower. The performance of the risk model was verified. Furthermore, a highly reliable nomogram was generated. Cellular experiments revealed the ability to promote cell viability, value-addedness, migratory ability, invasive ability, and colony-forming ability of the glioma cell line LN229 and HS-683. The qRT-PCR analysis of clinical glioma samples showed that these 10 genes were expressed at higher levels in high-grade gliomas than in low-grade gliomas, suggesting that these genes are associated with poor prognosis of gliomas.

Conclusion

Our study sifted out ten invasion-related biomarkers of LGG, providing a reference for treatments and prognostic prediction in LGG.

A first comprehensive analysis of Transcribed Ultra Conserved Regions uncovers important regulatory functions of novel non-coding transcripts in gliomas

Transcribed Ultra-Conserved Regions (TUCRs) represent a severely understudied class of putative non-coding RNAs (ncRNAs) that are 100% conserved across multiple species. We performed the first-ever analysis of TUCRs in glioblastoma (GBM) and low-grade gliomas (LGG). We leveraged large human datasets to identify the genomic locations, chromatin accessibility, transcription, differential expression, correlation with survival, and predicted functions of all 481 TUCRs, and identified TUCRs that are relevant to glioma biology. Of these, we investigated the expression, function, and mechanism of action of the most highly upregulated intergenic TUCR, uc.110, identifying it as a new oncogene. Uc.110 was highly overexpressed in GBM and LGG, where it promoted malignancy and tumor growth. Uc.110 activated the WNT pathway by upregulating the expression of membrane frizzled-related protein (MFRP), by sponging the tumor suppressor microRNA miR-544. This pioneering study shows important roles for TUCRs in gliomas and provides an extensive database and novel methods for future TUCR research.

LncRNA NDUFA6-DT: A Comprehensive Analysis of a Potential LncRNA Biomarker and Its Regulatory Mechanisms in Gliomas

Gliomas are the most prevalent primary malignant tumors affecting the brain, with high recurrence and mortality rates. Accurate diagnoses and effective treatment challenges persist, emphasizing the need for identifying new biomarkers to guide clinical decisions. Long noncoding RNAs (lncRNAs) hold potential as diagnostic and therapeutic biomarkers in cancer. However, only a limited subset of lncRNAs in gliomas have been explored. Therefore, this study aims to identify lncRNA signatures applicable to patients with gliomas across all grades and explore their clinical significance and potential biological mechanisms. Data used in this study were obtained from TCGA, CGGA, and GEO datasets to identify key lncRNA signatures in gliomas through differential and survival analyses and machine learning algorithms. We examined their associations with the clinical characteristics, gene mutations, diagnosis, and prognosis of gliomas. Functional enrichment analysis was employed to elucidate the potential biological mechanisms associated with these significant lncRNA signatures. We explored competing endogenous RNA (ceRNA) regulatory networks. We found that NDUFA6-DT emerged as a significant lncRNA signature in gliomas, with reduced NDUFA6-DT expression associated with a worse prognosis in gliomas. Nomogram analysis incorporating NDUFA6-DT expression levels exhibited excellent prognostic and predictive capabilities. Functional annotation suggested that NDUFA6-DT might influence immunological responses and synaptic transmission, potentially modifying glioma initiation and progression. The associated ceRNA network revealed the possible presence of the NDUFA6-DT-miR-455-3p-YWHAH/YWHAG axis in low-grade glioma (LGG) and glioblastoma multiforme (GBM), regulating the PI3K-AKT signaling pathway and influencing glioma cell survival and apoptosis. We believe that NDUFA6-DT is a novel lncRNA linked to glioma diagnosis and prognosis, potentially becoming a pivotal biomarker for glioma.

Identification of a novel disulfideptosis-related gene signature for prognostic implication in lower-grade gliomas

Lower-grade gliomas (GBMLGG) are common, fatal, and difficult-to-treat cancers. The current treatment choices have impressive efficacy constraints. As a result, the development of effective treatments and the identification of new therapeutic targets are urgent requirements. Disulfide metabolism is the cause of the non-apoptotic programmed cell death known as disulfideptosis, which was only recently discovered. The mRNA expression data and related clinical information of GBMLGG patients downloaded from public databases were used in this study to investigate the prognostic significance of genes involved in disulfideptosis. In the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohort, our findings showed that many disulfidptosis-related genes were expressed differently in normal and GBMLGG tissues. It was discovered that IQ motif-containing GTPase-activating protein 1 (IQGAP1) is a key gene that influences the outcome of GBMLGG. Besides, a nomogram model was built to foresee the visualization of GBMLGG patients. In addition, in vivo and in vitro validation of IQGAP1's cancer-promoting function was done. In conclusion, we discovered a gene signature associated with disulfideptosis that can effectively predict OS in GBMLGG patients. As a result, treating disulfideptosis may be a viable alternative for GBMLGG patients.

APOBEC3C is a novel target for the immune treatment of lower-grade gliomas

BACKGROUND

Apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) type 3C (A3C) has been identified as a cancer molecular biomarker in the past decade. However, the practical role of A3C in lower-grade gliomas (LGGs) in improving the clinical outcome remains unclear. This study aims to discuss the function of A3C in immunotherapy in LGGs.

METHODS

The RNA-Sequencing (RNA-seq) and corresponding clinical data were extracted from UCSC Xena and the results were verified in the Chinese Glioma Genome Atlas (CGGA). Weighted gene co-expression network analysis (WGCNA) was used for screening A3C-related genes. Comprehensive bioinformation analyses were performed and multiple levels of expression, survival rate, and biological functions were assessed to explore the functions of A3C.

RESULTS

A3C expression was significantly higher in LGGs than in normal tissues but lower than in glioblastoma (GBM), indicating its role as an independent prognosis predictor for LGGs. Twenty-eight A3C-related genes were found with WGCNA for unsupervised clustering analysis and three modification patterns with different outcomes and immune cell infiltration were identified. A3C and the A3C score were also correlated with immune cell infiltration and the expression of immune checkpoints. In addition, the A3C score was correlated with increased sensitivity to chemotherapy. Single-cell RNA (scRNA) analysis indicated that A3C most probably expresses on immune cells, such as T cells, B cells and macrophage.

CONCLUSIONS

A3C is an immune-related prognostic biomarker in LGGs. Developing drugs to block A3C could enhance the efficiency of immunotherapy and improve disease survival.Abbreviation: A3C: Apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) type 3C; LGGs: lower-grade gliomas; CGGA: Chinese Glioma Genome Atlas; WGCNA: Weighted gene co-expression network analysis; scRNA: Single-cell RNA; HGG: higher-grade glioma; OS: overall survival; TME: tumor microenvironment; KM: Kaplan-Meier; PFI: progression-free interval; IDH: isocitrate dehydrogenase; ROC: receiver operating characteristic; GS: gene significance; MM: module membership; TIMER: Tumor IMmune Estimation Resource; GSVA: gene set variation analysis; ssGSEA: single-sample gene-set enrichment analysis; PCA: principal component analysis; AUC: area under ROC curve; HAVCR2: hepatitis A virus cellular receptor 2; PDCD1: programmed cell death 1; PDCD1LG2: PDCD1 ligand 2; PTPRC: protein tyrosine phosphatase receptor type C; ACC: Adrenocortical carcinoma; BLCA: Bladder Urothelial Carcinoma;BRCA: Breast invasive carcinoma; CESC: Cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOLCholangiocarcinoma; COADColon adenocarcinoma; DLBC: Lymphoid Neoplasm Diffuse Large B-cell Lymphoma; ESCA: Esophageal carcinoma; GBM: Glioblastoma multiforme; HNSC: Head and Neck squamous cell carcinoma; KICH: Kidney Chromophobe; KIRC: Kidney renal clear cell carcinoma; KIRP: Kidney renal papillary cell carcinoma; LAML: Acute Myeloid Leukemia; LGG: Brain Lower Grade Glioma; LIHC: Liver hepatocellular carcinoma; LUAD: Lung adenocarcinoma; LUSC: Lung squamous cell carcinoma; MESO: Mesothelioma; OV: Ovarian serous cystadenocarcinoma; PAAD: Pancreatic adenocarcinoma; PCPG: Pheochromocytoma and Paraganglioma; PRAD: Prostate adenocarcinoma; READ: Rectum adenocarcinoma; SARC: Sarcoma; SKCM: Skin Cutaneous Melanoma; STAD: Stomach adenocarcinoma; TGCT: Testicular Germ Cell Tumors; THCA: Thyroid carcinoma; THYM: Thymoma; UCEC: Uterine Corpus Endometrial Carcinoma; UCS: Uterine Carcinosarcoma; UVM: Uveal Melanoma.

An inflammatory gene-related prognostic risk score model for prognosis and immune infiltration in glioblastoma

This study aimed to screen for key genes related to the prognosis of patients with glioblastoma (GBM). First, bioinformatics analysis was performed based on databases such as TCGA and MSigDB. Inflammatory-related genes were obtained from the MSigDB database. The TCGA-tumor samples were divided into cluster A and B groups based on consensus clustering. Multivariate Cox regression was applied to construct the risk score model of inflammatory-related genes based on the TCGA database. Second, to understand the effects of model characteristic genes on GBM cells, U-87 MG cells were used for knockdown experiments, which are important means for studying gene function. PLAUR is an unfavorable prognostic biomarker for patients with glioma. Therefore, the model characteristic gene PLAUR was selected for knockdown experiments. The prognosis of cluster A was significantly better than that of cluster B. The verification results also demonstrate that the risk score could predict overall survival. Although the immune cells in cluster B and high-risk groups increased, no matching survival advantage was observed. It may be that stromal activation inhibits the antitumor effect of immune cells. PLAUR knockdown inhibits tumor cell proliferation, migration, and invasion, and promoted tumor cell apoptosis. In conclusion, a prognostic prediction model for GBM composed of inflammatory-related genes was successfully constructed. Increased immune cell expression may be linked to a poor prognosis for GBM, as stromal activation decreased the antitumor activity of immune cells in cluster B and high-risk groups. PLAUR may play an important role in tumor cell proliferation, migration, invasion, and apoptosis.

Intra-operative applications of augmented reality in glioma surgery: a systematic review

Background

Augmented reality (AR) is increasingly being explored in neurosurgical practice. By visualizing patient-specific, three-dimensional (3D) models in real time, surgeons can improve their spatial understanding of complex anatomy and pathology, thereby optimizing intra-operative navigation, localization, and resection. Here, we aimed to capture applications of AR in glioma surgery, their current status and future potential.

Methods

A systematic review of the literature was conducted. This adhered to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guideline. PubMed, Embase, and Scopus electronic databases were queried from inception to October 10, 2022. Leveraging the Population, Intervention, Comparison, Outcomes, and Study design (PICOS) framework, study eligibility was evaluated in the qualitative synthesis. Data regarding AR workflow, surgical application, and associated outcomes were then extracted. The quality of evidence was additionally examined, using hierarchical classes of evidence in neurosurgery.

Results

The search returned 77 articles. Forty were subject to title and abstract screening, while 25 proceeded to full text screening. Of these, 22 articles met eligibility criteria and were included in the final review. During abstraction, studies were classified as "development" or "intervention" based on primary aims. Overall, AR was qualitatively advantageous, due to enhanced visualization of gliomas and critical structures, frequently aiding in maximal safe resection. Non-rigid applications were also useful in disclosing and compensating for intra-operative brain shift. Irrespective, there was high variance in registration methods and measurements, which considerably impacted projection accuracy. Most studies were of low-level evidence, yielding heterogeneous results.

Conclusions

AR has increasing potential for glioma surgery, with capacity to positively influence the onco-functional balance. However, technical and design limitations are readily apparent. The field must consider the importance of consistency and replicability, as well as the level of evidence, to effectively converge on standard approaches that maximize patient benefit.

Mutation-derived, genomic instability-associated lncRNAs are prognostic markers in gliomas

Background

Gliomas are the most commonly-detected malignant tumors of the brain. They contain abundant long non-coding RNAs (lncRNAs), which are valuable cancer biomarkers. LncRNAs may be involved in genomic instability; however, their specific role and mechanism in gliomas remains unclear. LncRNAs that are related to genomic instability have not been reported in gliomas.

Methods

The transcriptome data from The Cancer Genome Atlas (TCGA) database were analyzed. The co-expression network of genomic instability-related lncRNAs and mRNA was established, and the model of genomic instability-related lncRNA was identified by univariate Cox regression and LASSO analyses. Based on the median risk score obtained in the training set, we divided the samples into high-risk and low-risk groups and proved the survival prediction ability of genomic instability-related lncRNA signatures. The results were verified in the external data set. Finally, a real-time quantitative polymerase chain reaction assay was performed to validate the signature.

Results

The signatures of 17 lncRNAs (LINC01579, AL022344.1, AC025171.5, LINC01116, MIR155HG, AC131097.3, LINC00906, CYTOR, AC015540.1, SLC25A21.AS1, H19, AL133415.1, SNHG18, FOXD3.AS1, LINC02593, AL354919.2 and CRNDE) related to genomic instability were identified. In the internal data set and Gene Expression Omnibus (GEO) external data set, the low-risk group showed better survival than the high-risk group (P < 0.001). In addition, this feature was identified as an independent risk factor, showing its independent prognostic value with different clinical stratifications. The majority of patients in the low-risk group had isocitrate dehydrogenase 1 (IDH1) mutations. The expression levels of these lncRNAs were significantly higher in glioblastoma cell lines than in normal cells.

Conclusions

Our study shows that the signature of 17 lncRNAs related to genomic instability has prognostic value for gliomas and could provide a potential therapeutic method for glioblastoma.

TAF1D Functions as a Novel Biomarker in Osteosarcoma

Background: The most frequent primary bone cancer in teenagers, osteosarcoma (OS), is particularly aggressive with a high mortality rate. Methods: By combining public databases, OS and non-cancer samples were obtained. The Wilcoxon test and standardized mean difference (SMD) were utilized to evaluate the mRNA expression level of TATA-box binding protein associated factor, RNA polymerase 1 subunit D (TAF1D). The potential of TAF1D to discriminate OS samples from non-cancer samples was revealed by summary receiver operating characteristic curve (sROC). To investigate the prognostic significance, Kaplan‒Meier curve and univariate Cox analysis were performed. Immunohistochemistry (IHC) was used to determine the TAF1D protein expression level. ESTIMATE algorithm and TIMER2.0 database were used to reveal the association between TAF1D expression and the immune microenvironment. Enrichment analysis and potential drug prediction were performed to clarify the underlying molecular mechanisms and possible therapeutic directions of TAF1D. Ultimately, the transcription factors (TFs) and the TAF1D binding site were predicted based on the Cistrome and JASPAR databases. Results: TAF1D was upregulated in OS at the mRNA and protein levels and possessed robust discriminatory power. TAF1D upregulation was suggestive of worse prognosis and enhancement of tumor purity in OS patients. The cell cycle was the most significantly enriched pathway, and NU.1025 was considered to be the potential target agent. Finally, MYC was identified as a TF that regulates the expression of TAF1D. Conclusions: Altogether, TAF1D has the potential to serve as a biological marker and therapeutic target in OS, which could offer new perspectives for OS treatment.

MAGOH is correlated with poor prognosis and is essential for cell proliferation in lower-grade glioma

OBJECTIVE

Mago-nashi homolog (MAGOH) has been shown to play a pivotal part in various tumors. However, its specific contribution in lower-grade glioma (LGG) is still unknown.

METHODS

Pan-cancer analysis was implemented to inspect the expression characteristics and prognostic significance of MAGOH in multiple tumors. The associations between MAGOH expression patterns and the pathological features of LGG were analyzed, as were the connections between MAGOH expression and the clinical traits, prognosis, biological activities, immune features, genomic variations, and responses to treatment in LGG. Additionally, in vitro studies were performed to detect the expression levels and biomedical functions of MAGOH in LGG.

RESULTS

Abnormally increased levels of MAGOH expression were connected with adverse prognosis in patients with several types of tumors, including LGG. Importantly, we found that levels of MAGOH expression were independent prognostic biomarker of patients with LGG. Increased MAGOH expression was also highly associated with several immune-related markers, immune cell infiltration, immune checkpoint genes (ICPGs), gene mutations, and responses to chemotherapy in patients with LGG. In vitro studies ascertained that abnormally increased MAGOH was essential for cell proliferation in LGG.

CONCLUSION

MAGOH is a valid predictive biomarker in LGG and may become a novel therapeutic target in these patients.

DUSP10 is a novel immune-related biomarker connected with survival and cellular proliferation in lower-grade glioma

OBJECTIVE

The role of dual-specificity phosphatase 10 (DUSP10) has been investigated in several types of cancer. Nevertheless, the underlying function of DUSP10 in lower-grade glioma (LGG) remains undetermined.

METHODS

We entirely determined the expression features and prognostic significance of DUSP10 in numerous tumors by implementing a pan-cancer analysis. Adjacently, we thoroughly inspected the correlation between DUSP10 expression and clinicopathologic features, prognosis, biological processes, immune traits, gene variations, and treatment responses based on the expression features in LGG. In vitro studies were conducted to detect the underlying functions of DUSP10 in LGG.

RESULTS

Unconventionally boosted DUSP10 expression and higher DUSP10 expression correlated with poorer prognosis were discovered in various tumors, including LGG. Fortunately, DUSP10 expression was proven to be an independent prognostic indicator of patients with LGG. Additionally, DUSP10 expression was tightly linked to the immune modulation, gene mutations, and response to immunotherapy/chemotherapy in LGG patients. In vitro studies illustrated that the DUSP10 was abnormally increased and pivotal for cell proliferation in LGG.

CONCLUSIONS

Collectively, we verified that DUSP10 was an independent prognostic indicator and may become a novelty target of targeted therapy of LGG.

A Novel Prognostic Pyroptosis-Related Gene Signature Correlates to Oxidative Stress and Immune-Related Features in Gliomas

Gliomas are highly invasive and aggressive tumors having the highest incidence rate of brain cancer. Identifying effective prognostic and potential therapeutic targets is necessitated. The relationship of pyroptosis, a form of programmed cellular death, with gliomas remains elusive. We constructed and validated a prognostic model for gliomas using pyroptosis-related genes. Differentially expressed pyroptosis-related genes were screened using the "limma" package. Based on LASSO-Cox regression, nine significant genes including CASP1, CASP3, CASP6, IL32, MKI67, MYD88, PRTN3, NOS1, and VIM were employed to construct a prognostic model in the TCGA cohort; the results were validated in the CGGA cohort. According to the median risk score, the patients were classified into two risk groups, namely, high- and low-risk groups. Patients at high risk had worse prognoses relative to those at low risk evidenced by the Kaplan-Meier curve analysis. The two groups exhibited differences in immune cell infiltration and TMB scores, with high immune checkpoint levels, TMB scores, and immune cell infiltration levels in the high-risk group. KEGG and GO analyses suggested enrichment in immune-related pathways. Furthermore, we found that the genes in our signature strongly correlated with oxidative stress-related pathways and the subgroups exhibited different ssGSEA scores. Some small molecules targeted the genes in the model, and we verified their drug sensitivities between the risk groups. The scRNA-seq dataset, GSE138794, was processed using the "Seurat" package to assess the level of risk gene expression in specific cell types. Finally, the MYD88 level was lowered in the U87 glioma cell line using si-RNA constructs. Cellular proliferation was impaired, and fewer pyroptosis-related cytokines were released upon exposure to LPS. In summary, we built a pyroptosis-related gene model that accurately classified glioma patients into high- and low-risk groups. The findings suggest that the signature may be an effective prognostic predictive tool for gliomas.

Inefficient splicing of long non-coding RNAs is associated with higher transcript complexity in human and mouse

Recent reports show that long non-coding RNAs (lncRNAs) have inefficient splicing and fewer alternative splice variants than mRNAs. Here, we have explored the efficiency of lncRNAs and mRNAs in producing various splice variants, given the number of exons in humans and mice. Intriguingly, lncRNAs produce more splice variants per exon, referred to as Transcript Complexity, than mRNAs. Most lncRNA splice variants are the product of the alternative last exon and exon skipping. LncRNAs and mRNAs with higher transcript complexity have shorter intron lengths. Longer exon length and GC/AG at 5'/3' splice sites are associated with higher transcript complexity in lncRNAs. Lastly, our results indicate that inefficient splicing of lncRNAs may facilitate multiple introns splicing and, thus, more spliced products per exon.

Identification of cuproptosis-related lncRNAs for prognosis and immunotherapy in glioma

Glioma is a highly invasive primary brain tumour, making it challenging to accurately predict prognosis for glioma patients. Cuproptosis is a recently discovered cell death attracting significant attention in the tumour field. Whether cuproptosis-related genes have prognostic predictive value has not been clarified. In this study, uni-/multi-variate Cox and Lasso regression analyses were applied to construct a risk model based on cuproptosis-related lncRNAs using TCGA and CGGA cohorts. A nomogram was constructed to quantify individual risk, including clinical and genic characteristics and risk. GO and KEGG analyses were used to define functional enrichment of DEGs. Tumour mutation burden (TMB) and immune checkpoint analyses were performed to evaluate potential responses to ICI therapy. Ten prognostic lncRNAs were obtained from Cox regression. Based on the median risk score, patients were divided into high- and low-risk groups. Either for grade 2-3 or for grade 4, glioma patients with high-risk exhibited significant poorer prognoses. The risk was an independent risk factor associated with overall survival. The high-risk group was functionally associated with immune responses and cancer-related pathways. The high-risk group was associated with higher TMB scores. The expression levels of many immune checkpoints in the high-risk group were significantly higher than those in the low-risk group. Differentiated immune pathways were primarily enriched in the IFN response, immune checkpoint and T-cell co-stimulation pathways. In conclusion, we established a risk model based on cuproptosis-related lncRNAs showing excellent prognostic prediction ability but also indicating the immuno-microenvironment status of glioma.

Latest updates on cellular and molecular biomarkers of gliomas

Gliomas are the most common central nervous system malignancies, compromising almost 80% of all brain tumors and is associated with significant mortality. The classification of gliomas has shifted from basic histological perspective to one that is based on molecular biomarkers. Treatment of this type of tumors consists currently of surgery, chemotherapy and radiation therapy. During the past years, there was a limited development of effective glioma diagnostics and therapeutics due to multiple factors including the presence of blood-brain barrier and the heterogeneity of this type of tumors. Currently, it is necessary to highlight the advantage of molecular diagnosis of gliomas to develop patient targeted therapies based on multiple oncogenic pathway. In this review, we will evaluate the development of cellular and molecular biomarkers for the diagnosis of gliomas and the impact of these diagnostic tools for better tailored and targeted therapies.

The cuproptosis-related signature associated with the tumor environment and prognosis of patients with glioma

Background

Copper ions are essential for cellular physiology. Cuproptosis is a novel method of copper-dependent cell death, and the cuproptosis-based signature for glioma remains less studied.

Methods

Several glioma datasets with clinicopathological information were collected from TCGA, GEO and CGGA. Robust Multichip Average (RMA) algorithm was used for background correction and normalization, cuproptosis-related genes (CRGs) were then collected. The TCGA-glioma cohort was clustered using ConsensusClusterPlus. Univariate Cox regression analysis and the Random Survival Forest model were performed on the differentially expressed genes to identify prognostic genes. The cuproptosis-signature was constructed by calculating CuproptosisScore using Multivariate Cox regression analysis. Differences in terms of genomic mutation, tumor microenvironment, and enrichment pathways were evaluated between high- or low-CuproptosisScore. Furthermore, drug response prediction was carried out utilizing pRRophetic.

Results

Two subclusters based on CRGs were identified. Patients in cluster2 had better clinical outcomes. The cuproptosis-signature was constructed based on CuproptosisScore. Patients with higher CuproptosisScore had higher WHO grades and worse prognosis, while patients with lower grades were more likely to develop IDH mutations or MGMT methylation. Univariate and Multivariate Cox regression analysis demonstrated CuproptosisScore was an independent prognostic factor. The accuracy of the signature in prognostic prediction was further confirmed in 11 external validation datasets. In groups with high-CuproptosisScore, PIK3CA, MUC16, NF1, TTN, TP53, PTEN, and EGFR showed high mutation frequency. IDH1, TP53, ATRX, CIC, and FUBP1 demonstrated high mutation frequency in low-CuproptosisScore group. The level of immune infiltration increased as CuproptosisScore increased. SubMap analysis revealed patients with high-CuproptosisScore may respond to anti-PD-1 therapy. The IC50 values of Bexarotene, Bicalutamide, Bortezomib, and Cytarabine were lower in the high-CuproptosisScore group than those in the low-CuproptosisScore group. Finally, the importance of IGFBP2 in TCGA-glioma cohort was confirmed.

Conclusion

The current study revealed the novel cuproptosis-based signature might help predict the prognosis, biological features, and appropriate treatment for patients with glioma.

Identification of a novel cuproptosis-related gene signature and integrative analyses in patients with lower-grade gliomas

Background

Cuproptosis is a newly discovered unique non-apoptotic programmed cell death distinguished from known death mechanisms like ferroptosis, pyroptosis, and necroptosis. However, the prognostic value of cuproptosis and the correlation between cuproptosis and the tumor microenvironment (TME) in lower-grade gliomas (LGGs) remain unknown.

Methods

In this study, we systematically investigated the genetic and transcriptional variation, prognostic value, and expression patterns of cuproptosis-related genes (CRGs). The CRG score was applied to quantify the cuproptosis subtypes. We then evaluated their values in the TME, prognostic prediction, and therapeutic responses in LGG. Lastly, we collected five paired LGG and matched normal adjacent tissue samples from Sun Yat-sen University Cancer Center (SYSUCC) to verify the expression of signature genes by quantitative real-time PCR (qRT-PCR) and Western blotting (WB).

Results

Two distinct cuproptosis-related clusters were identified using consensus unsupervised clustering analysis. The correlation between multilayer CRG alterations with clinical characteristics, prognosis, and TME cell infiltration were observed. Then, a well-performed cuproptosis-related risk model (CRG score) was developed to predict LGG patients' prognosis, which was evaluated and validated in two external cohorts. We classified patients into high- and low-risk groups according to the CRG score and found that patients in the low-risk group showed significantly higher survival possibilities than those in the high-risk group (P<0.001). A high CRG score implies higher TME scores, more significant TME cell infiltration, and increased mutation burden. Meanwhile, the CRG score was significantly correlated with the cancer stem cell index, chemoradiotherapy sensitivity-related genes and immune checkpoint genes, and chemotherapeutic sensitivity, indicating the association with CRGs and treatment responses. Univariate and multivariate Cox regression analyses revealed that the CRG score was an independent prognostic predictor for LGG patients. Subsequently, a highly accurate predictive model was established for facilitating the clinical application of the CRG score, showing good predictive ability and calibration. Additionally, crucial CRGs were further validated by qRT-PCR and WB.

Conclusion

Collectively, we demonstrated a comprehensive overview of CRG profiles in LGG and established a novel risk model for LGG patients' therapy status and prognosis. Our findings highlight the potential clinical implications of CRGs, suggesting that cuproptosis may be the potential therapeutic target for patients with LGG.

Next-Generation Grade and Survival Expression Biomarkers of Human Gliomas Based on Algorithmically Reconstructed Molecular Pathways

In gliomas, expression of certain marker genes is strongly associated with survival and tumor type and often exceeds histological assessments. Using a human interactome model, we algorithmically reconstructed 7494 new-type molecular pathways that are centered each on an individual protein. Each single-gene expression and gene-centric pathway activation was tested as a survival and tumor grade biomarker in gliomas and their diagnostic subgroups (IDH mutant or wild type, IDH mutant with 1p/19q co-deletion, MGMT promoter methylated or unmethylated), including the three major molecular subtypes of glioblastoma (proneural, mesenchymal, classical). We used three datasets from The Cancer Genome Atlas and the Chinese Glioma Genome Atlas, which in total include 527 glioblastoma and 1097 low grade glioma profiles. We identified 2724 such gene and 2418 pathway survival biomarkers out of total 17,717 genes and 7494 pathways analyzed. We then assessed tumor grade and molecular subtype biomarkers and with the threshold of AUC > 0.7 identified 1322/982 gene biomarkers and 472/537 pathway biomarkers. This suggests roughly two times greater efficacy of the reconstructed pathway approach compared to gene biomarkers. Thus, we conclude that activation levels of algorithmically reconstructed gene-centric pathways are a potent class of new-generation diagnostic and prognostic biomarkers for gliomas.

A TP53 Related Immune Prognostic Model for the Prediction of Clinical Outcomes and Therapeutic Responses in Lung Adenocarcinoma

TP53 is the most frequently mutated gene in lung adenocarcinoma (LUAD). The tumor immune microenvironment (TIM) is considered a vital factor that influences tumor progression and survival rate. The influence of TP53 mutation on TIM in LUAD has not been fully studied. Here we systematically investigated the relationship and potential mechanisms between TP53 mutation status and immune response in LUAD. We constructed an immune prognostic model (IPM) using immune associated genes, which were expressed differentially between the TP53 mutant and wild type LUAD patients. We discovered that TP53 mutations were significantly associated with 5 immune related biological processes. Thirty-six immune genes were expressed differentially between TP53 mutant and wild type LUAD patients. An IPM was constructed using 3 immune genes to differentiate the prognostic survival in LUAD. The high-risk LUAD group displayed significantly higher proportions of dendritic cell resting, T cell CD4 memory resting and mast cell resting, and significantly low proportions of dendritic cell activated, T cell CD4 memory activated, and mast cell activated. Moreover, IPM was found to be an independent clinical feature and can be used to predict immunotherapy responses. In summary, we constructed and validated an IPM using 3 immune related genes, which provides a better understanding of the mechanism from an immunological perspectives.

Construction of a Risk Model to Predict the Prognosis and Immunotherapy of Low-Grade Glioma Ground on 7 Ferroptosis-Related Genes

Purpose

Ferroptosis is closely associated with tumors. The purpose of this study was to investigate the correlation between ferroptosis and prognosis of low grade glioma (LGG) via construction and verification of a risk model.

Patients and Methods

The data of LGG were downloaded from public databases. Through LASSO analysis of characteristic genes, a gene signature was constructed. Patients into were divided two groups based on risk score. Subsequently, survival, clinical phenotype, functional enrichment, immune cell infiltration and somatic mutation analysis were performed. In addition, whether ferroptosis-related genes (FRGs) signature can predict the patient's response to anti-PD-1/PD-L1 immunotherapy was also investigated.

Results

FRGs signature had strong prognostic assessment ability, and high risk score was associated with poor overall survival (OS) of LGG. The high risk score group had higher degree of immune cell infiltration, stronger stromal activity, higher immune score, and high expression of immune checkpoint. In low risk score group anti-PD-1/PD-L1 immunotherapy has significant therapeutic advantages and clinical response. Genes and frequency of somatic mutations and clinical phenotypes in the high and low risk score groups were significantly different.

Conclusion

A prognostic model based on 7 FRGs can be used to predict the prognosis and immunotherapeutic response of LGG.

APOLLO: An accurate and independently validated prediction model of lower-grade gliomas overall survival and a comparative study of model performance

BACKGROUND

Virtually few accurate and robust prediction models of lower-grade gliomas (LGG) survival exist that may aid physicians in making clinical decisions. We aimed to develop a prognostic prediction model of LGG by incorporating demographic, clinical and transcriptional biomarkers with either main effects or gene-gene interactions.

METHODS

Based on gene expression profiles of 1,420 LGG patients from six independent cohorts comprising both European and Asian populations, we proposed a 3-D analysis strategy to develop and validate an Accurate Prediction mOdel of Lower-grade gLiomas Overall survival (APOLLO). We further conducted decision curve analysis to assess the net benefit (NB) of identifying true positives and the net reduction (NR) of unnecessary interventions. Finally, we compared the performance of APOLLO and the existing prediction models by the first systematic review.

FINDINGS

APOLLO possessed an excellent discriminative ability to identify patients at high mortality risk. Compared to those with less than the 20th percentile of APOLLO risk score, patients with more than the 90th percentile of APOLLO risk score had significantly worse overall survival (HR=54·18, 95% CI: 34·73-84·52, P=2·66 × 10-69). Further, APOLLO can accurately predict both 36- and 60-month survival in six independent cohorts with a pooled AUC36-month=0·901 (95% CI: 0·879-0·923), AUC60-month=0·843 (95% CI: 0·815-0·871) and C-index=0·818 (95% CI: 0·800-0·835). Moreover, APOLLO offered an effective screening strategy for detecting LGG patients susceptible to death (NB36-month=0·166, NR36-month=40·1% and NB60-month=0·258, NR60-month=19·2%). The systematic comparisons revealed APOLLO outperformed the existing models in accuracy and robustness.

INTERPRETATION

APOLLO has the demonstrated feasibility and utility of predicting LGG survival (http://bigdata.njmu.edu.cn/APOLLO).

FUNDING

National Key Research and Development Program of China (2016YFE0204900); Natural Science Foundation of Jiangsu Province (BK20191354); National Natural Science Foundation of China (81973142 and 82103946); China Postdoctoral Science Foundation (2020M681671); National Institutes of Health (CA209414, CA249096, CA092824 and ES000002).

The Value of H2BC12 for Predicting Poor Survival Outcomes in Patients With WHO Grade II and III Gliomas

Purpose: Glioma is a common primary malignant brain tumor. Grade II (GII) gliomas are prone to develop into anaplastic grade III (GIII) gliomas, which indicate a higher malignancy and poorer survival outcome. This study aimed to satisfy the increasing demand for novel sensitive biomarkers and potential therapeutic targets in the treatment of GII and GIII gliomas.

Methods: A TCGA dataset was used to investigate the expression of H2BC12 mRNA in GII and GIII gliomas and its relation to clinical pathologic characteristics. Glioma tissues were collected to verify results from the TCGA dataset, and H2BC12 mRNA was detected by RT-qPCR. ROC analysis was employed to evaluate the classification power for GII and GIII. The significance of H2BC12 mRNA GII and GIII gliomas was also investigated. In addition, H2BC12 expression-related pathways were enriched by gene set enrichment analysis (GSEA). DNA methylation level and mutation of H2BC12 were analyzed by the UALCAN and CBioPortal databases, respectively.

Results: Based on the sample data from multiple databases and RT-qPCR, higher expression of H2BC12 mRNA was found in GII and GIII glioma tissue compared to normal tissue, which was consistent with a trend with our clinical specimen. H2BC12 mRNA had a better power in distinguishing between GII and GIII and yielded an AUC of 0.706 with a sensitivity of 76.9% and specificity of 81.8%. Meanwhile, high H2BC12 levels were associated with IDH status, 1p/19q codeletion, primary therapy outcome, and the histological type of gliomas. Moreover, the overall survival (OS), disease-specific survival (DSS), and progress-free interval (PFI) of GII glioma patients with higher levels of H2BC12 were shorter than those of patients with lower levels as well as GIII patients. In the multivariate analysis, a high H2BC12 level was an independent predictor for poor survival outcomes of gliomas. The Wnt or PI3K-AKT signaling pathways, DNA repair, cellular senescence, and DNA double-strand break repair were differentially activated in phenotypes that were positively associated with H2BC12. H2BC12 DNA methylation was high in TP53 nonmutant patients, and no H2BC12 mutation was observed in gliomas patients.

Conclusion: H2BC12 is a promising biomarker for the diagnosis and prognosis of patients with WHO grade II and III gliomas.

A 1p/19q Codeletion-Associated Immune Signature for Predicting Lower Grade Glioma Prognosis

Lower grade gliomas (LGGs) with codeletion of chromosomal arms 1p and 19q (1p/19 codeletion) have a favorable outcome. However, its overall survival (OS) varies. Here, we established an immune signature associated with 1p/19q codeletion for accurate prediction of prognosis of LGGs. The Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA) databases with RNA sequencing and corresponding clinical data were dichotomized into training group and testing group. The immune-related differentially expressed genes (DEGs) associated with 1p/19q codeletion were screened using Cox proportional hazards regression analyses. A prognostic signature was established using dataset from CGGA and tested in TCGA database. Subsequently, we explored the correlation between the prognostic signature and immune response. Thirteen immune genes associated with 1p/19q codeletion were used to construct a prognostic signature. The 1-, 3-, 5-year survival rates of the low-risk group were approximately 97%, 89%, and 79%, while those of the high-risk group were 81%, 50% and 34%, respectively, in the training group. The nomogram which comprised age, WHO grade, primary or recurrent types, 1p/19q codeletion status and risk score provided accurate prediction for the survival rate of glioma. DEGs that were highly expressed in the high-risk group clustered with many immune-related pathways. Immune checkpoints including TIM3, PD1, PDL1, CTLA4, TIGIT, MIR155HG, and CD48 were correlated with the risk score. VAV3 and TNFRFSF11B were found to be candidate immune checkpoints associated with prognosis. The 1p/19q codeletion-associated immune signature provides accurate prediction of OS. VAV3 and TNFRFSF11B are novel immune checkpoints.

Testing lncRNAs signature as clinical stage–related prognostic markers in gastric cancer progression using TCGA database

LncRNA expression can be conducive to gastric cancer (GC) prognosis. The objective of this study is to ascertain five specific lncRNAs involved in tumor progression of GC and their role as prognostic markers to diagnose clinical stage-wise GC. High-throughput RNA sequencing data were obtained from The Cancer Genome Atlas (TCGA) database and performed genome-wide lncRNA expression analysis using edgeR package, Bioconductor.org , and R-statistical computing to analyze differentially expressed lncRNA analysis. Cutoff parameters were FDR < 0.05 and |Log2FC| > 2. Total 351 tumor samples with differentially expressed lncRNAs were divided into group-1 lncRNAs such as AC019117.2 and LINC00941, and group-2 lncRNAs such as LINC02410, AC012317.2, and AC141273.1 by 2:1. The Spearman correlation coefficients ( p < 0.05) and correlation test function (cor.test ()) were performed for lncRNAs as per clinical stage. Cytoscape software was used to construct lncRNA–mRNA interaction networks. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway ( p < 0.05) analysis were conducted using the clusterProfiler package. Kaplan–Meier survival analysis was performed to determine the overall survival of patients based on the expression of five lncRNAs in different clinical stages of GC. AC019117.2 and LINC00941 of group 1 inferred a positive correlation with clinical stages of stage I to stage IV, and their expressions were higher in tumor tissues than normal tissues. On the contrary, LINC02410, AC012317.2, and AC141273.1 of group 2 exhibited a negative correlation with clinical stage, and they exhibited more expression in normal tissues compared to tumor tissues. GO and KEGG pathway analysis reported that AC019117.2 may interact with LINC00941 via ITGA3 and trophoblast glycoprotein (TPBG) to foster tumor progression. Tumor-specific group-1 lncRNAs were conducive to the poor overall survival and exhibited a positive correlation with the clinical stages of stage I to stage IV in GC as per the lncRNA–mRNA networking analysis. These five lncRNAs could be considered as clinically useful lncRNA-based prognostic markers to predict clinical stage-wise GC progression.

Transcriptional Networks Identify BRPF1 as a Potential Drug Target Based on Inflammatory Signature in Primary Lower-Grade Gliomas

Gliomas are the most common tumors of the central nervous system and are classified into grades I-IV based on their histological characteristics. Lower-grade gliomas (LGG) can be divided into grade II diffuse low-grade gliomas and grade III moderate gliomas and have a relatively good prognosis. However, LGG often develops into high-grade glioma within a few years. This study aimed to construct and identify the prognostic value of an inflammatory signature and discover potential drug targets for primary LGG. We first screened differentially expressed genes in primary LGG (TCGA) compared with normal brain tissue (GTEx) that overlapped with inflammation-related genes from MSigDB. After survival analysis, nine genes were selected to construct an inflammatory signature. LGG patients with a high inflammatory signature score had a poor prognosis, and the inflammatory signature was a strong independent prognostic factor in both the training cohort (TCGA) and validation cohort (CGGA). Compared with the low-inflammatory signature group, differentially expressed genes in the high-inflammatory signature group were mainly enriched in immune-related signaling pathways, which is consistent with the distribution of immune cells in the high- and low-inflammatory signature groups. Integrating driver genes, upregulated genes and drug targets data, bromodomain and PHD finger-containing protein 1 (BRPF1) was selected as a potential drug target. Inhibition of BRPF1 function or knockdown of BRPF1 expression attenuated glioma cell proliferation and colony formation.

Genome-wide methylomic analyses identify prognostic epigenetic signature in lower grade glioma

Glioma is the most malignant and aggressive type of brain tumour with high heterogeneity and mortality. Although some clinicopathological factors have been identified as prognostic biomarkers, the individual variants and risk stratification in patients with lower grade glioma (LGG) have not been fully elucidated. The primary aim of this study was to identify an efficient DNA methylation combination biomarker for risk stratification and prognosis in LGG. We conducted a retrospective cohort study by analysing whole genome DNA methylation data of 646 patients with LGG from the TCGA and GEO database. Cox proportional hazard analysis was carried out to screen and construct biomarker model that predicted overall survival (OS). The Kaplan‐Meier survival curves and time‐dependent ROC were constructed to prove the efficiency of the signature. Then, another independent cohort was used to further validate the finding. A two‐CpG site DNA methylation signature was identified by multivariate Cox proportional hazard analysis. Further analysis indicated that the signature was an independent survival predictor from other clinical factors and exhibited higher predictive accuracy compared with known biomarkers. This signature was significantly correlated with immune‐checkpoint blockade, immunotherapy‐related signatures and ferroptosis regulator genes. The expression pattern and functional analysis showed that these two genes corresponding with two methylation sites contained in the model were correlated with immune infiltration level, and involved in MAPK and Rap1 signalling pathway. The signature may contribute to improve the risk stratification of patients and provide a more accurate assessment for precision medicine in the clinic.

SP1-induced lncRNA ZFPM2 antisense RNA 1 (ZFPM2-AS1) aggravates glioma progression via the miR-515-5p/Superoxide dismutase 2 (SOD2) axis

Glioma is a common life-threatening tumor with high malignancy and high invasiveness. LncRNA ZFPM2 antisense RNA 1 (ZFPM2-AS1) was confirmed to be implicated in numerous tumors, while its biological function and mechanism have not been thoroughly understood in glioma. The gene expression was measured by RT-qPCR. Cell proliferation, cell cycle, and cell apoptosis of glioma cells were validated by CCK-8, colony formation, flow cytometry and TUNEL assays. The effect of ZFPM2-AS1 on tumor growth was verified by in vivo assay. The exploration on ZFPM2-AS1-mediated mechanism was carried out via ChIP, luciferase reporter, and RIP assays. In the present study, ZFPM2-AS1 was demonstrated as a highly-expressed lncRNA in glioma tissues and cells. ZFPM2-AS1 silencing suppressed cell proliferation and cell cycle, but facilitated cell apoptosis. In addition, the inhibitive effect of silenced ZFPM2-AS1 was also observed in tumor growth. Furthermore, we found that SP1 interacted with ZFPM2-AS1 promoter to transcriptionally activate ZFPM2-AS1 expression. Moreover, ZFPM2-AS1 was identified as a competing endogenous RNA (ceRNA) for miR-515-5p to target SOD2. Rescue assays verified that SOD2 overexpression partially abolished the suppressive impact of ZFPM2-AS1 silencing on glioma cell growth. In conclusion, this study corroborated the regulatory mechanism of SP1/ZFPM2-AS1/miR-515-5p/SOD2 axis in glioma, indicating that targeting ZFPM2-AS1 might be an effective way to treat glioma.

Development and Validation of an Immune-Related Gene Pairs Signature in Grade II/III Glioma

BACKGROUND

Gliomas are prevalent primary intracerebral malignant tumors. Increasing evidence indicates an association between the immune signature and Grade II/III glioma prognosis. Thus, we aimed to develop an immune-related gene pair (IRGP) signature that can be used as a prognostic tool in Grade II/III glioma.

METHODS

The gene expression levels and clinical information of Grade II/III glioma patients were collected from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases. The TCGA data were randomly divided into a training cohort (n = 249) and a validation cohort (n = 162), and a CGGA dataset served as an external validation group (n = 605). IRGPs significantly associated with prognosis were selected by Cox regression. Gene set enrichment analysis and filtration were performed with the IRGPs.

RESULTS

Within a set of 1991 immune genes, 8 IRGPs including 15 unique genes that significantly affect survival constituted a gene signature. In the validation datasets, the IRGP signature significantly stratified patients with Grade II/III glioma into low- and high-risk groups (P < 0.001), and the IRGP index was found to be an independent prognostic factor through univariate and multivariate analyses (P < 0.05). Additionally, 26 functional pathways were identified through the intersection of Gene Set Enrichment Analysis (GSEA) and Gene Ontology (GO) enrichment analysis.

CONCLUSION

The IRGP signature demonstrated good prognostic value for Grade II/III gliomas, which may provide new insights into individual treatment for glioma patients. The IRGPs might function through the identified 26 functional pathways.

Identification and Potential Mechanisms of a 7-MicroRNA Signature That Predicts Prognosis in Patients with Lower-Grade Glioma

Background

Lower-grade glioma is an intracranial cancer that may develop into glioblastoma with high mortality. The main objective of our study is to develop microRNA for LGG patients which will provide novel prognostic biomarkers along with therapeutic targets.

Methods

Clinicopathological data of LGG patients and their RNA expression profile were downloaded through The Cancer Genome Atlas Relevant expression profiles of RNA, and clinicopathological data of the LGG patients had been extracted from the database of "The Cancer Genome Atlas." Differential expression analysis had been conducted for identification of the differentially expressed microRNAs as well as mRNAs in LGG samples and normal ones. ROC curves and K-M plots were plotted to confirm performance and for predictive accuracy. For the confirmation of microRNAs as an independent prognostic factor, an independent prognosis analysis was conducted. Moreover, target differentially expressed genes of these identified prognostic microRNAs that were extracted and protein-protein interaction networks were developed. Moreover, the biological functions of signature were determined through Genome Ontology analysis, genome pathway analysis, and Kyoto Encyclopedia of Genes.

Results

7-microRNA signature was identified that has the ability of categorization of individuals with LGG into high- and low-risk groups on the basis of significant difference in survival during training and testing cohorts (P < 0.001). The 7-microRNA signature had appeared to be robust in predictive accuracy (all AUC> 0.65). It was also approved with multivariate Cox regression along with some traditional clinical practices that we can use 7-microRNA signature for therapeutic purposes as a self-regulating predictive OS factor (P < 0.001). KEGG and Gene Ontology (GO) analyses reported that 7-microRNAs had mainly developed in important pathways related with glioma, e.g., the "cAMP signaling pathway," "glutamatergic synapses," and "calcium signaling pathway".

Conclusion

A newly discovered 7-microRNA signature could be a potential target for the diagnosis and treatment for LGG patients.

Identification of biomarkers and ceRNA network in glioblastoma through bioinformatic analysis and evaluation of potential prognostic values

Background

Glioblastoma (GBM) is one of the most common and malignant primary brain tumors in adults, with high mortality rates and limited treatment. Based on bioinformatic analyses, this study aimed to identify biomarkers and relevant molecular pathways that may serve as potential targets for the treatment of GBM.

Methods

Expression profiles were downloaded from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database; nine GBM samples and three normal samples were extracted from the GSE104267 dataset. Differentially-expressed messenger RNA (mRNA) and long non-coding RNA (lncRNA) were screened from the preprocessed dataset. The clusterProfiler package in R was used to perform a biological process (BP) analysis of gene ontology (GO), and a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed separately in upregulated and downregulated groups. A competing endogenous RNA (ceRNA) network was constructed using Cytoscape. Based on data downloaded from The Cancer Genome Atlas (TCGA), Kaplan-Meier (K-M) survival curves were established. Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed to evaluate IL10RB antisense RNA 1 (IL10RB-AS1) expression in GBM tissue compared with that in normal brain tissue.

Results

A total of 253 differentially-expressed genes (DEGs) were obtained. Based on BP and KEGG enrichment annotation analyses, 11 lncRNA-related pathways were identified through function prediction analysis. A competing endogenous RNA (ceRNA) subnetwork, including 21 nodes and 29 regulatory pairs, was then constructed. Based on the clinical data of GBM in TCGA, one survival-related DEG, IL10RB-AS1, was identified using the log-rank statistical test. K-M survival curves of IL10RB-AS1 and expression levels of IL10RB-AS1 in both GBM and normal brain tissue were obtained.

Conclusions

Through the combination of bioinformatic analyses, one survival-related differentially-expressed lncRNA, IL10RB-AS1, was identified. This, along with several related signaling pathways and ceRNA systems that were elucidated in GBM have potential prognostic value and might offer new possibilities for the treatment of GBM.

Comprehensive Analysis of Inflammatory Response-Related Genes, and Prognosis and Immune Infiltration in Patients With Low-Grade Glioma

Background: Although low-grade glioma (LGG) has a good prognosis, it is prone to malignant transformation into high-grade glioma. It has been confirmed that the characteristics of inflammatory factors and immune microenvironment are closely related to the occurrence and development of tumors. It is necessary to clarify the role of inflammatory genes and immune infiltration in LGG.

Methods: We downloaded the transcriptome gene expression data and corresponding clinical data of LGG patients from the TCGA and GTEX databases to screen prognosis-related differentially expressed inflammatory genes with the difference analysis and single-factor Cox regression analysis. The prognostic risk model was constructed by LASSO Cox regression analysis, which enables us to compare the overall survival rate of high- and low-risk groups in the model by Kaplan–Meier analysis and subsequently draw the risk curve and survival status diagram. We analyzed the accuracy of the prediction model via ROC curves and performed GSEA enrichment analysis. The ssGSEA algorithm was used to calculate the score of immune cell infiltration and the activity of immune-related pathways. The CellMiner database was used to study drug sensitivity.

Results: In this study, 3 genes (CALCRL, MMP14, and SELL) were selected from 9 prognosis-related differential inflammation genes through LASSO Cox regression analysis to construct a prognostic risk model. Further analysis showed that the risk score was negatively correlated with the prognosis, and the ROC curve showed that the accuracy of the model was better. The age, grade, and risk score can be used as independent prognostic factors (p < 0.001). GSEA analysis confirmed that 6 immune-related pathways were enriched in the high-risk group. We found that the degree of infiltration of 12 immune cell subpopulations and the scores of 13 immune functions and pathways in the high-risk group were significantly increased by applying the ssGSEA method (p < 0.05). Finally, we explored the relationship between the genes in the model and the susceptibility of drugs.

Conclusion: This study analyzed the correlation between the inflammation-related risk model and the immune microenvironment. It is expected to provide a reference for the screening of LGG prognostic markers and the evaluation of immune response.

Immunogenomic Analyses of the Prognostic Predictive Model for Patients With Renal Cancer

Background

Renal cell carcinoma (RCC) is associated with poor prognostic outcomes. The current stratifying system does not predict prognostic outcomes and therapeutic benefits precisely for RCC patients. Here, we aim to construct an immune prognostic predictive model to assist clinician to predict RCC prognosis.

Methods

Herein, an immune prognostic signature was developed, and its predictive ability was confirmed in the kidney renal clear cell carcinoma (KIRC) cohorts based on The Cancer Genome Atlas (TCGA) dataset. Several immunogenomic analyses were conducted to investigate the correlations between immune risk scores and immune cell infiltrations, immune checkpoints, cancer genotypes, tumor mutational burden, and responses to chemotherapy and immunotherapy.

Results

The immune prognostic signature contained 14 immune-associated genes and was found to be an independent prognostic factor for KIRC. Furthermore, the immune risk score was established as a novel marker for predicting the overall survival outcomes for RCC. The risk score was correlated with some significant immunophenotypic factors, including T cell infiltration, antitumor immunity, antitumor response, oncogenic pathways, and immunotherapeutic and chemotherapeutic response.

Conclusions

The immune prognostic, predictive model can be effectively and efficiently used in the prediction of survival outcomes and immunotherapeutic responses of RCC patients.

Mechanisms of Long Non-Coding RNAs in Biological Characteristics and Aerobic Glycolysis of Glioma

Glioma is the most common and aggressive tumor of the central nervous system. The uncontrolled proliferation, cellular heterogeneity, and diffusive capacity of glioma cells contribute to a very poor prognosis of patients with high grade glioma. Compared to normal cells, cancer cells exhibit a higher rate of glucose uptake, which is accompanied with the metabolic switch from oxidative phosphorylation to aerobic glycolysis. The metabolic reprogramming of cancer cell supports excessive cell proliferation, which are frequently mediated by the activation of oncogenes or the perturbations of tumor suppressor genes. Recently, a growing body of evidence has started to reveal that long noncoding RNAs (lncRNAs) are implicated in a wide spectrum of biological processes in glioma, including malignant phenotypes and aerobic glycolysis. However, the mechanisms of diverse lncRNAs in the initiation and progression of gliomas remain to be fully unveiled. In this review, we summarized the diverse roles of lncRNAs in shaping the biological features and aerobic glycolysis of glioma. The thorough understanding of lncRNAs in glioma biology provides opportunities for developing diagnostic biomarkers and novel therapeutic strategies targeting gliomas.

Establishment of an Immune-Related Gene Signature for Risk Stratification for Patients with Glioma

Glioma is a frequently seen primary malignant intracranial tumor, characterized by poor prognosis. The study is aimed at constructing a prognostic model for risk stratification in patients suffering from glioma. Weighted gene coexpression network analysis (WGCNA), integrated transcriptome analysis, and combining immune-related genes (IRGs) were used to identify core differentially expressed IRGs (DE IRGs). Subsequently, univariate and multivariate Cox regression analyses were utilized to establish an immune-related risk score (IRRS) model for risk stratification for glioma patients. Furthermore, a nomogram was developed for predicting glioma patients' overall survival (OS). The turquoise module (cor = 0.67; P < 0.001) and its genes (n = 1092) were significantly pertinent to glioma progression. Ultimately, multivariate Cox regression analysis constructed an IRRS model based on VEGFA, SOCS3, SPP1, and TGFB2 core DE IRGs, with a C-index of 0.811 (95% CI: 0.786-0.836). Then, Kaplan-Meier (KM) survival curves revealed that patients presenting high risk had a dismal outcome (P < 0.0001). Also, this IRRS model was found to be an independent prognostic indicator of gliomas' survival prediction, with HR of 1.89 (95% CI: 1.252-2.85) and 2.17 (95% CI: 1.493-3.14) in the Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) datasets, respectively. We established the IRRS prognostic model, capable of effectively stratifying glioma population, convenient for decision-making in clinical practice.

Promoting Prognostic Model Application: A Review Based on Gliomas

Malignant neoplasms are characterized by poor therapeutic efficacy, high recurrence rate, and extensive metastasis, leading to short survival. Previous methods for grouping prognostic risks are based on anatomic, clinical, and pathological features that exhibit lower distinguishing capability compared with genetic signatures. The update of sequencing techniques and machine learning promotes the genetic panels-based prognostic model development, especially the RNA-panel models. Gliomas harbor the most malignant features and the poorest survival among all tumors. Currently, numerous glioma prognostic models have been reported. We systematically reviewed all 138 machine-learning-based genetic models and proposed novel criteria in assessing their quality. Besides, the biological and clinical significance of some highly overlapped glioma markers in these models were discussed. This study screened out markers with strong prognostic potential and 27 models presenting high quality. Conclusively, we comprehensively reviewed 138 prognostic models combined with glioma genetic panels and presented novel criteria for the development and assessment of clinically important prognostic models. This will guide the genetic models in cancers from laboratory-based research studies to clinical applications and improve glioma patient prognostic management.

Promoting Prognostic Model Application: A Review Based on Gliomas

Malignant neoplasms are characterized by poor therapeutic efficacy, high recurrence rate, and extensive metastasis, leading to short survival. Previous methods for grouping prognostic risks are based on anatomic, clinical, and pathological features that exhibit lower distinguishing capability compared with genetic signatures. The update of sequencing techniques and machine learning promotes the genetic panels-based prognostic model development, especially the RNA-panel models. Gliomas harbor the most malignant features and the poorest survival among all tumors. Currently, numerous glioma prognostic models have been reported. We systematically reviewed all 138 machine-learning-based genetic models and proposed novel criteria in assessing their quality. Besides, the biological and clinical significance of some highly overlapped glioma markers in these models were discussed. This study screened out markers with strong prognostic potential and 27 models presenting high quality. Conclusively, we comprehensively reviewed 138 prognostic models combined with glioma genetic panels and presented novel criteria for the development and assessment of clinically important prognostic models. This will guide the genetic models in cancers from laboratory-based research studies to clinical applications and improve glioma patient prognostic management.

LINC02381 contributes to cell proliferation and hinders cell apoptosis in glioma by transcriptionally enhancing CBX5

Glioma, featured with high incidence and low survival rate, is the most common type of primary brain tumor, severely affecting human life worldwide. LINC02381 is an interesting lncRNA functioning as oncogenic lncRNA in some cancers but as tumor-suppressor in others, but no report demonstrates its association with and function in glioma. Intriguingly, we found in a bioinformatics website LncRNADisease that LINC02381 was closely related to malignant glioma, so this study aimed to figure out the expression and function of LINC02381 in glioma. By RT-qPCR, we confirmed LINC02381 upregulation in glioma cells. Functional experiments demonstrated that LINC02381 knockdown repressed glioma cell proliferation and induced apoptosis. Boinformatics tools and RT-qPCR revealed the positive correlation between LINC02381 and CBX5 in glioma cells. More importantly, we confirmed that LINC02381 could interact and work synergistically with CEBPβ to bind to CBX5 promoter and activate CBX5 transcriptionally. Additionally, rescue experiments indicated that CBX5 up-regulation reversed the decline in cell proliferation and the augment in cell apoptosis caused by LINC02381 knockdown. To conclude, LINC02381 could facilitate CBX5 transcription via interaction with CEBPβ, thus exerting its oncogenic role in glioma cells, which could contribute to better understanding of glioma.

Identification and Validation of an 11-Ferroptosis Related Gene Signature and Its Correlation With Immune Checkpoint Molecules in Glioma

Background

Glioma is the most common primary malignant brain tumor with significant mortality and morbidity. Ferroptosis, a novel form of programmed cell death (PCD), is critically involved in tumorigenesis, progression and metastatic processes.

Methods

We revealed the relationship between ferroptosis-related genes and glioma by analyzing the mRNA expression profiles from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), GSE16011, and the Repository of Molecular Brain Neoplasia Data (REMBRANDT) datasets. The least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to construct a ferroptosis-associated gene signature in the TCGA cohort. Glioma patients from the CGGA, GSE16011, and REMBRANDT cohorts were used to validate the efficacy of the signature. Receiver operating characteristic (ROC) curve analysis was applied to measure the predictive performance of the risk score for overall survival (OS). Univariate and multivariate Cox regression analyses of the 11-gene signature were performed to determine whether the ability of the prognostic signature in predicting OS was independent. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to identify the potential biological functions and pathways of the signature. Subsequently, we performed single sample gene set enrichment analysis (ssGSEA) to explore the correlation between risk scores and immune status. Finally, seven putative small molecule drugs were predicted by Connectivity Map.

Results

The 11-gene signature was identified to divide patients into two risk groups. ROC curve analysis indicated the 11-gene signature as a potential diagnostic factor in glioma patients. Multivariate Cox regression analyses showed that the risk score was an independent predictive factor for overall survival. Functional analysis revealed that genes were enriched in iron-related molecular functions and immune-related biological processes. The results of ssGSEA indicated that the 11-gene signature was correlated with the initiation and progression of glioma. The small molecule drugs we selected showed significant potential to be used as putative drugs.

Conclusion

we identified a novel ferroptosis-related gene signature for prognostic prediction in glioma patients and revealed the relationship between ferroptosis-related genes and immune checkpoint molecules.

Immune-Related Gene SERPINE1 Is a Novel Biomarker for Diffuse Lower-Grade Gliomas via Large-Scale Analysis

Background

Glioma is one of the highly fatal primary tumors in the central nervous system. As a major component of tumor microenvironment (TME), immune cell has been proved to play a critical role in the progression and prognosis of the diffuse lower-grade gliomas (LGGs). This study aims to screen the key immune-related factors of LGGs by investigating the TCGA database.

Methods

The RNA-sequencing data of 508 LGG patients were downloaded in the TCGA database. ESTIMATE algorithm was utilized to calculate the stromal, immune, and ESTIMATE scores, based on which, the differentially expressed genes (DEGs) were analyzed by using “limma” package. Cox regression analysis and the cytoHubba plugin of Cytoscape software were subsequently applied to screen the survival-related genes and hub genes, the intersection of which led to the identification of SERPINE1 that played key roles in the LGGs. The expression patterns, clinical features, and regulatory mechanisms of SERPINE1 in the LGGs were further analyzed by data mining of the TCGA database. What’s more, the above analyses of SERPINE1 were further validated in the LGG cohort from the CGGA database.

Result

We found that stromal and immune cell infiltrations were strongly related to the prognosis and malignancy of the LGGs. A total of 54 survival-related genes and 46 hub genes were screened out in the DEGs, within which SERPINE1 was identified to be significantly overexpressed in the LGG samples compared with the normal tissues. Moreover, the upregulation of SERPINE1 was more pronounced in the gliomas of WHO grade III and IDH wild type, and its expression was correlated with poor prognosis in the LGG patients. The independent prognostic value of SERPINE1 in the LGG patients was also confirmed by Cox regression analysis. In terms of the functions of SERPINE1, the results of enrichment analysis indicated that SERPINE1 was mainly enriched in the immune‐related biological processes and signaling pathways. Furthermore, it was closely associated with infiltrations of immune cells in the LGG microenvironment and acted synergistically with PD1, PD-L1, PD-L2.

Conclusion

These findings proved that SERPINE1 could serve as a prognostic biomarker and potential immunotherapy target of LGGs.

A novel tailored immune gene pairs signature for overall survival prediction in lower-grade gliomas

•

The Immune-related gene pairs (IRGPs) pronostic signature for LGG is correlated with immune cells infiltration.

•

WGCNA presented a gene set correlating with immune cells infiltration and genes co-expression relationships were visualized.

•

The nomogram constrcted by three IRGPs and clinical factors is a novel tailored tool for individual-level prediction.

Lower-grade gliomas (LGGs) have a good prognosis with a wide range of overall survival (OS) outcomes. An accurate prognostic system can better predict survival time. An RNA-Sequencing (RNA-seq) prognostic signature showed a better predictive power than clinical predictor models. A signature constructed using gene pairs can transcend changes from biological heterogeneity, technical biases, and different measurement platforms. RNA-seq coupled with corresponding clinical information were extracted from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). Immune-related gene pairs (IRGPs) were used to establish a prognostic signature through univariate and multivariate Cox proportional hazards regression. Weighted gene co-expression network analysis (WGCNA) was used to evaluate module eigengenes correlating with immune cell infiltration and to construct gene co-expression networks. Samples in the training and testing cohorts were dichotomized into high- and low-risk groups. Risk score was identified as an independent predictor, and exhibited a closed relationship with prognosis. WGCNA presented a gene set that was positively correlated with age, WHO grade, isocitrate dehydrogenase (IDH) mutation status, 1p/19 codeletion, risk score, and immune cell infiltrations (CD4 T cells, B cells, dendritic cells, and macrophages). A nomogram comprising of age, WHO grade, 1p/19q codeletion, and three gene pairs (BIRC5|SSTR2, BMP2|TNFRSF12A, and NRG3|TGFB2) was established as a tool for predicting OS. The IPGPs signature, which is associated with immune cell infiltration, is a novel tailored tool for individual-level prediction.

Identification and validation of a five-lncRNA prognostic signature related to Glioma using bioinformatics analysis

Background

To accurately predict the prognosis of glioma patients.

Methods

A total of 541 samples from the TCGA cohort, 181 observations from the CGGA database and 91 samples from our cohort were included in our study. Long non-coding RNAs (LncRNAs) associated with glioma WHO grade were evaluated by weighted gene co-expression network analysis (WGCNA). Five lncRNA features were selected out to construct prognostic signatures based on the Cox regression model.

Results

By weighted gene co-expression network analysis (WGCNA), 14 lncRNAs related to glioma grade were identified. Using univariate and multivariate Cox analysis, five lncRNAs (CYTOR, MIR155HG, LINC00641, AC120036.4 and PWAR6) were selected to develop the prognostic signature. The Kaplan-Meier curve depicted that the patients in high risk group had poor prognosis in all cohorts. The areas under the receiver operating characteristic curve of the signature in predicting the survival of glioma patients at 1, 3, and 5 years were 0.84, 0.92, 0.90 in the CGGA cohort; 0.8, 0.85 and 0.77 in the TCGA set and 0.72, 0.90 and 0.86 in our own cohort. Multivariate Cox analysis demonstrated that the five-lncRNA signature was an independent prognostic indicator in the three sets (CGGA set: HR = 2.002, p < 0.001; TCGA set: HR = 1.243, p = 0.007; Our cohort: HR = 4.457, p = 0.008, respectively). A nomogram including the lncRNAs signature and clinical covariates was constructed and demonstrated high predictive accuracy in predicting 1-, 3- and 5-year survival probability of glioma patients.

Conclusion

We established a five-lncRNA signature as a potentially reliable tool for survival prediction of glioma patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-07972-9.

A Five Immune-Related lncRNA Signature as a Prognostic Target for Glioblastoma

Background: A variety of regulatory approaches including immune modulation have been explored as approaches to either eradicate antitumor response or induce suppressive mechanism in the glioblastoma microenvironment. Thus, the study of immune-related long noncoding RNA (lncRNA) signature is of great value in the diagnosis, treatment, and prognosis of glioblastoma.

Methods: Glioblastoma samples with lncRNA sequencing and corresponding clinical data were acquired from the Cancer Genome Atlas (TCGA) database. Immune-lncRNAs co-expression networks were built to identify immune-related lncRNAs via Pearson correlation. Based on the median risk score acquired in the training set, we divided the samples into high- and low-risk groups and demonstrate the survival prediction ability of the immune-related lncRNA signature. Both principal component analysis (PCA) and gene set enrichment analysis (GSEA) were used for immune state analysis.

Results: A cohort of 151 glioblastoma samples and 730 immune-related genes were acquired in this study. A five immune-related lncRNA signature (AC046143.1, AC021054.1, AC080112.1, MIR222HG, and PRKCQ-AS1) was identified. Compared with patients in the high-risk group, patients in the low-risk group showed a longer overall survival (OS) in the training, validation, and entire TCGA set (p = 1.931e-05, p = 1.706e-02, and p = 3.397e-06, respectively). Additionally, the survival prediction ability of this lncRNA signature was independent of known clinical factors and molecular features. The area under the ROC curve (AUC) and stratified analyses were further performed to verify its optimal survival predictive potency. Of note, the high-and low-risk groups exhibited significantly distinct immune state according to the PCA and GSEA analyses.

Conclusions: Our study proposes that a five immune-related lncRNA signature can be utilized as a latent indicator of prognosis and potential therapeutic approach for glioblastoma.

A novel Foxp3-related immune prognostic signature for glioblastoma multiforme based on immunogenomic profiling

Foxp3⁺ regulatory T cells (Treg) play an important part in the glioma immunosuppressive microenvironment. This study analyzed the effect of Foxsp3 on the immune microenvironment and constructed a Foxp3-related immune prognostic signature (IPS)for predicting prognosis in glioblastoma multiforme (GBM). Immunohistochemistry (IHC) staining for Foxp3 was performed in 72 high-grade glioma specimens. RNA-seq data from 152 GBM samples were obtained from The Cancer Genome Atlas database (TCGA) and divided into two groups, Foxp3 High (Foxp3_H) and Foxp3 Low (Foxp3_L), based on Foxp3 expression. We systematically analyzed the influence of Foxp3 on the immune microenvironment. Least Absolute Shrinkage and Selection Operator (LASSO) Cox analysis was conducted for immune-related genes that were differentially expressed between Foxp3_H and Foxp3_L GBM patients. We found a differential expression of Foxp3 in high-grade glioma tissues. The presence of Foxp3 was significantly associated with poor OS. From the four-gene IPS developed, GBM patients were stratified into low-risk and high-risk groups in both the training set and validation sets. Furthermore, we developed a novel nomogram to evaluate the overall survival in GBM patients. This study offers innovative insights into the GBM immune microenvironment and these findings contribute to individualized treatment and improvement in the prognosis for GBM patients.

Integrative Analysis of Neuregulin Family Members-Related Tumor Microenvironment for Predicting the Prognosis in Gliomas

Gliomas, including brain lower grade glioma (LGG) and glioblastoma multiforme (GBM), are the most common primary brain tumors in the central nervous system. Neuregulin (NRG) family proteins belong to the epidermal growth factor (EGF) family of extracellular ligands and they play an essential role in both the central and peripheral nervous systems. However, roles of NRGs in gliomas, especially their effects on prognosis, still remain to be elucidated. In this study, we obtained raw counts of RNA-sequencing data and corresponding clinical information from 510 LGG and 153 GBM samples from The Cancer Genome Atlas (TCGA) database. We analyzed the association of NRG1-4 expression levels with tumor immune microenvironment in LGG and GBM. GSVA (Gene Set Variation Analysis) was performed to determine the prognostic difference of NRGs gene set between LGG and GBM. ROC (receiver operating characteristic) curve and the nomogram model were constructed to estimate the prognostic value of NRGs in LGG and GBM. The results demonstrated that NRG1-4 were differentially expressed in LGG and GBM in comparison to normal tissue. Immune score analysis revealed that NRG1-4 were significantly related to the tumor immune microenvironment and remarkably correlated with immune cell infiltration. The investigation of roles of m⁶A (N6-methyladenosine, m⁶A)-related genes in gliomas revealed that NRGs were prominently involved in m⁶A RNA modification. GSVA score showed that NRG family members are more associated with prognosis in LGG compared with GBM. Prognostic analysis showed that NRG3 and NRG1 can serve as potential independent biomarkers in LGG and GBM, respectively. Moreover, GDSC drug sensitivity analysis revealed that NRG1 was more correlated with drug response compared with other NRG subtypes. Based on these public databases, we preliminarily identified the relationship between NRG family members and tumor immune microenvironment, and the prognostic value of NRGs in gliomas. In conclusion, our study provides comprehensive roles of NRG family members in gliomas, supporting modulation of NRG signaling in the management of glioma.