Comprehensive analysis of genes based on chr1p/19q co-deletion reveals a robust 4-gene prognostic signature for lower grade glioma

Glycolysis‑related lncRNA may be associated with prognosis and immune activity in grade II‑III glioma

Glucose metabolism, as a novel theory to explain tumor cell behavior, has been intensively studied in various tumors. The present study explored the long non-coding RNAs (lncRNAs) related to glycolysis in grade II-III glioma, aiming to provide a promising target for further research. Pearson correlation analysis was used to identify glycolysis-related lncRNAs. Univariate/multivariate Cox regression analysis and the Least Absolute Shrinkage and Selection Operator algorithm were applied to identify glycolysis-related lncRNAs to construct a prognosis prediction model. Subsequently, multi-dimensional evaluations were used to verify whether the risk model could predict the prognosis and survival rate of patients with grade II-III glioma. Finally, it was verified by functional experiments. The present study finally identified seven glycolysis-related lncRNAs (CRNDE, AC022034.1, RHOQ-AS1, AL159169.2, AL133215.2, AC007098.1 and LINC02587) to construct a prognosis prediction model. The present study further investigated the underlying immune microenvironment, somatic landscape and functional enrichment pathways. Additionally, individualized immunotherapeutic strategies and candidate compounds were identified to guide clinical treatment. The experimental results demonstrated that CRNDE could increase the proliferation of SHG-44 cells. In conclusion, a large sample of human grade II-III glioma in The Cancer Genome Atlas database was used to construct a risk model using glycolysis-related lncRNAs to predict the prognosis of patients with grade II-III glioma.

Pan-Cancer Interrogation of MUTYH Variants Reveals Biallelic Inactivation and Defective Base Excision Repair Across a Spectrum of Solid Tumors

PURPOSE

Biallelic germline pathogenic variants of the base excision repair (BER) pathway gene MUTYH predispose to colorectal cancer (CRC) and other cancers. The possible association of heterozygous variants with broader cancer susceptibility remains uncertain. This study investigated the prevalence and consequences of pathogenic MUTYH variants and MUTYH loss of heterozygosity (LOH) in a large pan-cancer analysis.

MATERIALS AND METHODS

Data from 354,366 solid tumor biopsies that were sequenced as part of routine clinical care were analyzed using a validated algorithm to distinguish germline from somatic MUTYH variants.

RESULTS

Biallelic germline pathogenic MUTYH variants were identified in 119 tissue biopsies. Most were CRCs and showed increased tumor mutational burden (TMB) and a mutational signature consistent with defective BER (COSMIC Signature SBS18). Germline heterozygous pathogenic variants were identified in 5,991 biopsies and their prevalence was modestly elevated in some cancer types. About 12% of these cancers (738 samples: including adrenal gland cancers, pancreatic islet cell tumors, nonglioma CNS tumors, GI stromal tumors, and thyroid cancers) showed somatic LOH for MUTYH, higher rates of chromosome 1p loss (where MUTYH is located), elevated genomic LOH, and higher COSMIC SBS18 signature scores, consistent with BER deficiency.

CONCLUSION

This analysis of MUTYH alterations in a large set of solid cancers suggests that in addition to the established role of biallelic pathogenic MUTYH variants in cancer predisposition, a broader range of cancers may possibly arise in MUTYH heterozygotes via a mechanism involving somatic LOH at the MUTYH locus and defective BER. However, the effect is modest and requires confirmation in additional studies before being clinically actionable.

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.

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.

Prognostic Value of Immune-Related lncRNA SBF2-AS1 in Diffuse Lower-Grade Glioma

LncRNA SET-binding factor 2 (SBF2) antisense RNA1 (SBF2-AS1) has been proven to play an oncogenic role in various types of tumors, but the prognostic role of SBF2-AS1 in tumors, especially in diffuse lower-grade glioma (LGG), is still unclear. Here, we aimed to investigate the prognostic value of SBF2-AS1 in LGG. The LGG expression profiles from The Cancer Genome Atlas (TCGA, n = 524) and Chinese Glioma Genome Atlas (CGGA, n = 431) were mined by Kaplan-Meier analysis, Cox regression analysis, Chi-square test and GSEA analysis. Through Kaplan-Meier analysis, we found the prognosis of LGG patients with high expression of SBF2-AS1 were worse than that of patients with low expression (Log Rank P < 0.001). Cox analysis showed SBF2-AS1 was an independent prognostic factor for poorer overall survival in LGG (P < 0.05). SBF2-AS1 was found to be significantly related to IDH mutation status and SBF2-AS1 was highly expressed in IDH wildtype group. GSEA analysis obtained a total of 126 GO terms and 6 KEGG pathways that were significantly enriched in SBF2-AS1 high expression phenotype (NOM P value < 0.05). We found these 126 GO terms and KEGG pathways were mainly related to immunity. In conclusion, lncRNA SBF2-AS1 expression is an immune-related lncRNA associated with unfavorable overall survival in LGG. SBF2-AS1 could be a reliable prognostic biomarker for patients with LGG.

Development of MR-based preoperative nomograms predicting DNA copy number subtype in lower grade gliomas with prognostic implication

OBJECTIVES

We aimed to determine the value of MR-based preoperative nomograms in predicting DNA copy number (CN) subtype in lower grade glioma (LGG) patients.

METHODS

The overall survival (OS) data were analyzed. MRI data of 170 subjects were retrospectively analyzed. The correlation was explored by univariate and multivariate regression analysis.

RESULTS

CN2 subtype was associated with shortest median OS (CN2 subtype vs. others: 46.8 vs. 221.7 months, p < 0.05). The time-dependent receiver operating characteristic for the CN2 subtype was 0.80 (95% CI: 0.74-0.85) for survival at 1 year, 0.80 (95% CI: 0.75-0.85) for survival at 2 years, and 0.77 (95% CI: 0.73-0.83) for survival at 3 years. On multivariate analysis, hemorrhage (OR: 0.118; p < 0.001; 95% CI: 0.037-0.376), poorly defined margin (OR: 4.592; p < 0.001; 95% CI: 1.965-10.730), extranodular growth (OR: 0.247; p = 0.006; 95% CI: 0.091-0.671), and volume ≥ 60 cm³ (OR: 4.734.256; p < 0.001; 95% CI: 2.051-10.924) were associated with CN1 subtype (AUC: 0.781). Proportion CE tumor (OR: 5.905; p = 0.007; 95% CI: 1.622-21.493), extranodular growth (OR: 9.047; p = 0.001; 95% CI: 2.349-34.846), width ≥ median (OR: 0.231; p = 0.049; 95% CI: 0.054-0.998), and depth ≥ median (OR: 0.192; p = 0.023; 95% CI: 0.046-0.799) were associated with CN2 subtype (AUC: 0.854). Necrosis/cystic (OR: 6.128; p = 0.007; 95% CI: 1.635-22.968), hemorrhage (OR: 5.752; p = 0.002; 95% CI: 1.953-16.942), poorly defined margin (OR: 0.164; p < 0.001; 95% CI: 0.063-0.427), and volume ≥ median (OR: 4.422; p < 0.001; 95% CI: 1.925-10.160) were associated with CN3 subtype (AUC: 0.808). All three nomograms showed good discrimination and calibration. Decision curve analysis supported that all nomograms were clinically useful. The average accuracy of the tenfold cross-validation was 0.680 (CN1), 0.794 (CN2), and 0.894 (CN3), respectively.

CONCLUSIONS

The shortest OS was observed in patients with CN2 subtype. This preliminary radiogenomics analysis revealed that the MR-based preoperative nomograms provide individualized prediction of DNA copy number subtype in LGG patients.

KEY POINTS

• This preliminary radiogenomics analysis of LGG revealed that the MR-based preoperative nomograms provide individualized prediction of DNA copy number subtype in LGG patients. • The AUC for the ROC curve was 0.781 for CN1 subtype, 0.854 for CN2 subtype, and 0.808 for CN3 subtype. Decision curve analysis supported that all nomograms were clinically useful. • The sensitivity was 0.779 (CN1), 0.731 (CN2), and 0.851 (CN3), respectively. The specificity was 0.664 (CN1), 0.872 (CN2), and 0.625 (CN3), respectively. And the accuracy was 0.717 (CN1), 0.849 (CN2), and 0.692 (CN3), respectively.

Identification and Validation of an Energy Metabolism-Related lncRNA-mRNA Signature for Lower-Grade Glioma

Energy metabolic processes play important roles for tumor malignancy, indicating that related protein-coding genes and regulatory upstream genes (such as long noncoding RNAs (lncRNAs)) may represent potential biomarkers for prognostic prediction. This study will develop a new energy metabolism-related lncRNA-mRNA prognostic signature for lower-grade glioma (LGG) patients. A GSE4290 dataset obtained from Gene Expression Omnibus was used for screening the differentially expressed genes (DEGs) and lncRNAs (DELs). The Cancer Genome Atlas (TCGA) dataset was used as the prognosis training set, while the Chinese Glioma Genome Atlas (CGGA) was for the validation set. Energy metabolism-related genes were collected from the Molecular Signatures Database (MsigDB), and a coexpression network was established between energy metabolism-related DEGs and DELs to identify energy metabolism-related DELs. Least absolute shrinkage and selection operator (LASSO) analysis was performed to filter the prognostic signature which underwent survival analysis and nomogram construction. A total of 1613 DEGs and 37 DELs were identified between LGG and normal brain tissues. One hundred and ten DEGs were overlapped with energy metabolism-related genes. Twenty-seven DELs could coexpress with 67 metabolism-related DEGs. LASSO regression analysis showed that 9 genes in the coexpression network were the optimal signature and used to construct the risk score. Kaplan-Meier curve analysis showed that patients with a high risk score had significantly worse OS than those with a low risk score (TCGA: HR = 3.192, 95%CI = 2.182‐4.670; CGGA: HR = 1.922, 95%CI = 1.431‐2.583). The predictive accuracy of the risk score was also high according to the AUC of the ROC curve (TCGA: 0.827; CGGA: 0.806). Multivariate Cox regression analyses revealed age, IDH1 mutation, and risk score as independent prognostic factors, and thus, a prognostic nomogram was established based on these three variables. The excellent prognostic performance of the nomogram was confirmed by calibration and discrimination analyses. In conclusion, our findings provided a new biomarker for the stratification of LGG patients with poor prognosis.

Prognostic model and nomogram construction based on autophagy signatures in lower grade glioma

The median survival time of lower grade glioma (LGG) tumors spans a wide range of 2-10 years and is highly dependent on the molecular characteristics and tumor location. Currently, there is no prognostic predictor for these tumors based on autophagy-related (ATG) genes. A prognostic risk score model based on the most significant seven ATG genes was established for LGG. These seven genes, including GRID2, FOXO1, MYC, PTK6, IKBKE, BIRC5, and TP73, have been screened as potentially therapeutic targets. The Kaplan-Meier survival curve analyses validated that patients with high or low risk scores had significantly different overall survival. Following the multivariate Cox regression and area under the ROC curve (AUC) analysis, a final prognostic model based on age, World Health Organization grade, 1p19q-codeletion status, and ATG risk score was performed as an independent prognostic indicator (training set: p = 4.09E-05, AUC = 0.901; validation set-1: p = .00069, AUC = 0.808; validation set-2: p = .0376, AUC = 0.830). Subsequently, a prognostic nomogram was constructed for individualized survival prediction. The calibration plots showed excellent predict efficiency between probability and actual overall survival. In this study, we provided several potential biomarkers for further developing potentially therapeutic targets of LGG. We also established a prognostic model and nomogram to improve the clinical glioma management and assist individualized survival prediction.

Identification and validation of a three-gene signature as a candidate prognostic biomarker for lower grade glioma

Background

Lower grade glioma (LGG) are a heterogeneous tumor that may develop into high-grade malignant glioma seriously shortens patient survival time. The clinical prognostic biomarker of lower-grade glioma is still lacking. The aim of our study is to explore novel biomarkers for LGG that contribute to distinguish potential malignancy in low-grade glioma, to guide clinical adoption of more rational and effective treatments.

Methods

The RNA-seq data for LGG was downloaded from UCSC Xena and the Chinese Glioma Genome Atlas (CGGA). By a robust likelihood-based survival model, least absolute shrinkage and selection operator regression and multivariate Cox regression analysis, we developed a three-gene signature and established a risk score to predict the prognosis of patient with LGG. The three-gene signature was an independent survival predictor compared to other clinical parameters. Based on the signature related risk score system, stratified survival analysis was performed in patients with different age group, gender and pathologic grade. The prognostic signature was validated in the CGGA dataset. Finally, weighted gene co-expression network analysis (WGCNA) was carried out to find the co-expression genes related to the member of the signature and enrichment analysis of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were conducted for those co-expression network. To prove the efficiency of the model, time-dependent receiver operating characteristic curves of our model and other models are constructed.

Results

In this study, a three-gene signature (WEE1, CRTAC1, SEMA4G) was constructed. Based on the model, the risk score of each patient was calculated with LGG (low-risk vs. high-risk, hazard ratio (HR) = 0.198 (95% CI [0.120-0.325])) and patients in the high-risk group had significantly poorer survival results than those in the low-risk group. Furthermore, the model was validated in the CGGA dataset. Lastly, by WGCNA, we constructed the co-expression network of the three genes and conducted the enrichment of GO and KEGG. Our study identified a three-gene model that showed satisfactory performance in predicting the 1-, 3- and 5-year survival of LGG patients compared to other models and may be a promising independent biomarker of LGG.