In recent years, development in high‐throughput sequencing technologies has experienced an increasing application of statistics, pattern recognition, and machine learning in bioinformatics analyses. SangeBox platform to meet different scientific demands. The new version of Sangs is a widely used tool among many researchers, which encourages us to continuously improve the plerBox 2 (http://vip.sangerbox.com) and extends and optimizes the functions of interactive graphics and analysis of clinical bioinformatics data. We introduced novel analytical tools such as random forests and support vector machines, as well as corresponding plotting functions. At the same time, we also optimized the performance of the platform and fixed known problems to allow users to perform data analyses more quickly and efficiently. SangerBox 2 improved the speed of analysis, reduced resource required for computer performance, and provided more analysis methods, greatly promoting the research efficiency.

Deep learning models have been developed for various predictions in glioma; yet, they were constrained by manual segmentation, task-specific design, or a lack of biological interpretation. Herein, we aimed to develop an end-to-end multi-task deep learning (MDL) pipeline that can simultaneously predict molecular alterations and histological grade (auxiliary tasks), as well as prognosis (primary task) in gliomas. Further, we aimed to provide the biological mechanisms underlying the model’s predictions. We collected multiscale data including baseline MRI images from 2776 glioma patients across two private (FAHZU and HPPH, n = 1931) and three public datasets (TCGA, n = 213; UCSF, n = 410; and EGD, n = 222). We trained and internally validated the MDL model using our private datasets, and externally validated it using the three public datasets. We used the model-predicted deep prognosis score (DPS) to stratify patients into low-DPS and high-DPS subtypes. Additionally, a radio-multiomics analysis was conducted to elucidate the biological basis of the DPS. In the external validation cohorts, the MDL model achieved average areas under the curve of 0.892–0.903, 0.710–0.894, and 0.850–0.879 for predicting IDH mutation status, 1p/19q co-deletion status, and tumor grade, respectively. Moreover, the MDL model yielded a C-index of 0.723 in the TCGA and 0.671 in the UCSF for the prediction of overall survival. The DPS exhibits significant correlations with activated oncogenic pathways, immune infiltration patterns, specific protein expression, DNA methylation, tumor mutation burden, and tumor-stroma ratio. Accordingly, our work presents an accurate and biologically meaningful tool for predicting molecular subtypes, tumor grade, and survival outcomes in gliomas, which provides personalized clinical decision-making in a global and non-invasive manner.

Background: To develop and validate a conventional MRI-based radiomic model for predicting prognosis in patients with IDH wild-type glioblastoma (GBM) and reveal the biological underpinning of the radiomic phenotypes. Methods: A total of 801 adult patients (training set, N = 471; internal validation set, N = 239; external validation set, N = 91) diagnosed with IDH wild-type GBM were included. A 20-feature radiomic risk score (Radscore) was built for overall survival (OS) prediction by univariate prognostic analysis and least absolute shrinkage and selection operator (LASSO) Cox regression in the training set. GSEA and WGCNA were applied to identify the intersectional pathways underlying the prognostic radiomic features in a radiogenomic analysis set with paired MRI and RNA-seq data (N = 132). The biological meaning of the conventional MRI sequences was revealed using a Mantel test. Results: Radscore was demonstrated to be an independent prognostic factor (P < 0.001). Incorporating the Radscore into a clinical model resulted in a radiomic-clinical nomogram predicting survival better than either the Radscore model or the clinical model alone, with better calibration and classification accuracy (a total net reclassification improvement of 0.403, P < 0.001). Three pathway categories (proliferation, DNA damage response, and immune response) were significantly correlated with the prognostic radiomic phenotypes. Conclusion: Our findings indicated that the prognostic radiomic phenotypes derived from conventional MRI are driven by distinct pathways involved in proliferation, DNA damage response, and immunity of IDH wild-type GBM.

Background: Isocitrate dehydrogenase (IDH) mutation and 1p19q codeletion status have been identified as significant markers for therapy and prognosis in lower-grade glioma (LGG). The current study aimed to construct a combined machine learning-based model for predicting the molecular subtypes of LGG, including (1) IDH wild-type astrocytoma (IDHwt), (2) IDH mutant and 1p19q non-codeleted astrocytoma (IDHmut-noncodel), and (3) IDH-mutant and 1p19q codeleted oligodendroglioma (IDHmut-codel), based on multiparametric magnetic resonance imaging (MRI) radiomics, qualitative features, and clinical factors. Methods: A total of 335 patients with LGG (WHO grade II/III) were retrospectively enrolled. The sum of 5,929 radiomics features were extracted from multiparametric MRI. Selected robust, non-redundant, and relevant features were used to construct a random forest model based on a training cohort (n = 269) and evaluated on a testing cohort (n = 66). Meanwhile, preoperative MRIs of all patients were scored in accordance with Visually Accessible Rembrandt Images (VASARI) annotations and T2-fluid attenuated inversion recovery (T2-FLAIR) mismatch sign. By combining radiomics features, qualitative features (VASARI annotations and T2-FLAIR mismatch signs), and clinical factors, a combined prediction model for the molecular subtypes of LGG was built. Results: The 17-feature radiomics model achieved area under the curve (AUC) values of 0.6557, 0.6830, and 0.7579 for IDHwt, IDHmut-noncodel, and IDHmut-codel, respectively, in the testing cohort. Incorporating qualitative features and clinical factors into the radiomics model resulted in improved AUCs of 0.8623, 0.8056, and 0.8036 for IDHwt, IDHmut-noncodel, and IDHmut-codel, with balanced accuracies of 0.8924, 0.8066, and 0.8095, respectively. Conclusion: The combined machine learning algorithm can provide a method to non-invasively predict the molecular subtypes of LGG preoperatively with excellent predictive performance.

Brain network dynamics not only endow the brain with flexible coordination for various cognitive processes but also with a huge potential of neuroplasticity for development, skill learning, and after cerebral injury. Diffusive and progressive glioma infiltration triggers the neuroplasticity for functional compensation, which is an outstanding pathophysiological model for the investigation of network reorganization underlying neuroplasticity. In this study, we employed dynamic conditional correlation to construct framewise language networks and investigated dynamic reorganizations in 83 patients with left hemispheric gliomas involving language networks (40 patients without aphasia and 43 patients with aphasia). We found that, in healthy controls (HCs) and patients, the language network dynamics in resting state clustered into 4 temporal-reoccurring states. Language deficits-severity-dependent topological abnormalities of dFCs were observed. Compared with HCs, suboptimal language network dynamics were observed for those patients without aphasia, while more severe network disruptions were observed for those patients with aphasia. Machine learning-based dFC-linguistics prediction analyses showed that dFCs of the 4 states significantly predicted individual patients’ language scores. These findings shed light on our understanding of metaplasticity in glioma. Glioma-induced language network reorganizations were investigated under a dynamic “meta-networking” (network of networks) framework. In healthy controls and patients with glioma, the framewise language network dynamics in resting-state robustly clustered into 4 temporal-reoccurring states. The spatial but not temporal language deficits-severity-dependent abnormalities of dFCs were observed in patients with left hemispheric gliomas involving language network. Language network dynamics significantly predicted individual patients’ language scores.