Alternations and Applications of the Structural and Functional Connectome in Gliomas: A Mini-Review.

In the central nervous system, gliomas are the most common, but complex primary tumors. Genome-based molecular and clinical studies have revealed different classifications and subtypes of gliomas. Neuroradiological approaches have non-invasively provided a macroscopic view for surgical resection and therapeutic effects. The connectome is a structural map of a physical object, the brain, which raises issues of spatial scale and definition, and it is calculated through diffusion magnetic resonance imaging (MRI) and functional MRI. In this study, we reviewed the basic principles and attributes of the structural and functional connectome, followed by the alternations of connectomes and their influences on glioma. To extend the applications of connectome, we demonstrated that a series of multi-center projects still need to be conducted to systemically investigate the connectome and the structural-functional coupling of glioma. Additionally, the brain-computer interface based on accurate connectome could provide more precise structural and functional data, which are significant for surgery and postoperative recovery. Besides, integrating the data from different sources, including connectome and other omics information, and their processing with artificial intelligence, together with validated biological and clinical findings will be significant for the development of a personalized surgical strategy.

Introduction

In the central nervous system (CNS), gliomas are the most common, but complex primary tumors. An average annual incidence rate of glioma is 6.7 per 100,000 individuals in China (Jiang et al., 2016; Zhao et al., 2021), which is similar to that in the United States (Ostrom et al., 2020). Genome-based molecular and clinical studies have revealed different classifications and subtypes of gliomas. Cancer-related molecular characteristics could be used to indicate the type of treatment modality, such as surgery and adjuvant therapies, as well as clinical outcomes. The 2021 World Health Organization (WHO) classification of CNS tumors (Wen and Packer, 2021) divided gliomas into 6 different types, including adult-type diffuse gliomas, pediatric-type diffuse low-grade gliomas (LGG), pediatric-type diffuse high-grade gliomas, circumscribed astrocytic gliomas, glioneuronal and neuronal tumors, and ependymomas. The classification and grading of gliomas combine clinical, pathological, histological, and especially molecular features. Molecular profiles, such as isocitrate dehydrogenase (IDH)-mutant, 1p19q co-deletion, the X-linked alpha-thalassemia mental retardation (ATRX)-mutant, Telomerase reverse transcriptase (TERT) promoter mutant, and O6-methylguanine–DNA methyltransferase (MGMT) promoter methylation, provide a microscopic view for predicting tumor biological behavior and patients’ outcomes. Besides, neuroradiological approaches non-invasively provide a macroscopic view for surgical resection and improving therapeutic effects.

In clinical practice, computed tomography (CT), magnetic resonance imaging (MRI), biphasic 2-deoxy-2-[18 F]fluoro-D-glucose positron emission tomography/CT (18F-FDG PET/CT), and neuronavigation are the main imaging methods to assess cancer metabolism, cancer progression, conditions of nerve fibers, and cortical and subcortical functions.

Along with the development of artificial intelligence (AI), quantitative data from CT and ordinary MR images, such as T1, T2, and T1-contrast-enhanced, and fluid-attenuated inversion recovery (FLAIR) images could be extracted and processed by machine learning or deep learning algorithms (Lambin et al., 2012; Chang et al., 2019). State-of-the-art neural networks could select shape, texture, and gray level of tumors for predicting tumor classification (Lee et al., 2020; Chakrabarty et al., 2021), segmentation (Tang et al., 2020; Rudie et al., 2021), grading (Naser and Deen, 2020), and even molecular parameters (Bangalore Yogananda et al., 2020; Choi et al., 2021).

Studies (Lundervold and Lundervold, 2019; Tomaszewski and Gillies, 2021) have investigated the imaging features, as well as their potential relationships with pathological and molecular characteristics. However, in clinical work, as some special gliomas are located in some functional areas or important subcortical areas, their physical properties with surrounding brain tissues and functional effects on both hemispheres are worthy of further investigation. Thus, special imaging sequences, including structural and functional MRI pulse sequences, can be exploited to determine the functional connectivity of tumors and other tissues, which are called brain tumor connectome.

The brain’s anatomical connectivity, or connectome, was defined as the map of neural connections in the brain (Hagmann, 2005). The connections emphasized the fact that the huge brain neuronal communication capacity and computational power critically relied on subtle and incredibly complex connectivity architecture through diffusion MRI (dMRI) (Mori and Barker, 1999; Sporns et al., 2005; Assaf and Pasternak, 2008), which was designed to assess the trajectory of protons in white matter fiber bundles, illustrating structural connectivity in the brain. Apart from dMRI, functional MRI (fMRI) (Ogawa et al., 1992; Logothetis, 2008) indirectly reflects the activity of electrical signals of neurons in a certain brain area by detecting changes in blood oxygen levels, thereby reflecting the functional connections of neurons. By analyzing the functional connections through task-fMRI or resting-state fMRI (rs-fMRI), we could construct the functional connectivity atlas, followed by establishment of the brain connectome that is composed of functional and structural connectivity. In general, as a key component of brain multi-omics, connectome, depending on improvement of imaging hardware and software (Wang et al., 2013), could be used to explore neuropsychiatric activities and their changes (Sporns, 2011). In the present study, we reviewed the structural–functional connectome attributes in normal human brains and their abnormalities and instabilities in glioma, and highlighted their potential applications in glioma clinical management.

Structural Connectome and the Alternations in Gliomas

Structural Connectome

According to the different degrees of free-water molecules in neurons, dMRI could distinguish the white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF) by imposing gradients of magnetic fields and controlling the time intervals through different b-values (Cohen and Assaf, 2002). The myelinated axon of a neuron in one region of the brain extends to another region following a particular anatomic course or trajectory. Thus, the graph theory (Euler, 1741; Semmel et al., 2021), proposed in the 18th century, could provide a robust mathematical and practical foundation for structural connectome. Specifically, after acquiring ordinary MRI-T1, -T2, and diffusion sequences, as well as preprocessing the images including space alignment, image denoise, and field correction in a standard workflow, we could estimate the response functions of WM, GM, and CSF. Then, the response functions were applied to deconvolve the orientations of fibers and construct the tractography by various atlases (Tzourio-Mazoyer et al., 2002; Desikan et al., 2006; Fan et al., 2016; Glasser et al., 2016). A brain structural network graph would be architected by measuring the strengths and directions of neurons (Figure 1). Numerous studies (Hagmann, 2005; He et al., 2007; Toga et al., 2012; Liao et al., 2017) analyzed the dMRI tractography and network topology of the brain and revealed the long-range axonal network and the hierarchical organization of the functional inter-cortical connectivity. A small-world network architecture (Watts and Strogatz, 1998) has mainly demonstrated a low-cost, but high clustered network organization. Hub neuronal nodes as the centralities in the network could be connected together to enable an efficient and robust signaling in the network. A rich club is defined as a set of high-degree nodes that is more densely interconnected than predicted on the basis of the node degrees alone (de Reus and van den Heuvel, 2013), suggesting the theoretical basement of anatomical organization for data transmission. It is also a conservative network architecture across various species (Watts and Strogatz, 1998; Oh et al., 2014; Beul et al., 2015) and spatial scales (Das, 2020). Previous studies (Collin and van den Heuvel, 2013; van den Heuvel et al., 2015) demonstrated that seeking the differences of brain connectome between preterm and term infant brains is a vital step for understanding the developmental impairment caused by preterm.

FIGURE 1

Workflow of diffusion MRI (dMRI) and resting-state functional MRI (fMRI). After image acquisition, images obtained by dMRI and fMRI were initially preprocessed by denoise, bias correction, motion correction, alignment, etc. Then, in dMRI, fiber orientation distribution (FOD), and fractional anisotropy (FA) were calculated for tractography, followed by structural connectome evaluation. In fMRI, time series sequences were used to assess the brain functional activities. Besides, seed-based approaches were used for the analysis of regions of interests (ROIs) and the graph theory was employed for the whole-brain functional connectome analysis.

Tractography constructed by high-resolution MRI systems provides the WM connections of cortico-cortical and cortico-subcortical spatial regions and their anatomic relationships with normal functions of awareness, cognition, behaviors, etc. Thus, abnormalities in structural connectome affected by tumors could reflect the biological behaviors, such as tumor progression and invasion.

Structural Connectome Alternations in Gliomas

In gliomas, a number of scholars (Venkataramani et al., 2019, 2021; Venkatesh et al., 2019) have reported a relevant direct synaptic communication and electrical integration between tumors cells and neurons, providing the potential explanations of tumors, involving the neuronal signaling pathways. Glutamatergic synapses could medicate connections between presynaptic neurons and postsynaptic glioma cells and form a positive current, promoting calcium-related invasiveness of tumor cells and glioma proliferation until it would be inhibited by glutamate receptor antagonists (Venkataramani et al., 2019). Meanwhile, glioma interconnections via gap junctions (Osswald et al., 2015) and membrane tubes (Jung et al., 2017; Weil et al., 2017) could construct a self-repairing and resistant network, disturbing the normal brain network, as well as promoting tumor progression and therapeutic resistance.

Therefore, implementation of dMRI and establishment of tractography might detect the accumulation of the abnormal cells. As we mentioned above, hub neuronal regions play a central role in the brain small-world network and can be clustered into a rich club. Glioma patients, however, have more hub-related connections in the contralateral hemisphere (Yu et al., 2016; Douw et al., 2019; Yuan T. et al., 2020), while a slight difference in the ipsilateral hemisphere was found. These contralesional connections suggested the decompensation of structure and function and the particular redistribution of connectivity (Crossley et al., 2014; Stam, 2014). The cascade network failure theory (Buldyrev et al., 2010; van den Heuvel and Sporns, 2019), indicating the overloading of functional hubs and the increase of connectivity, has been studied in Alzheimer’s disease (Jones et al., 2016) and schizophrenia (van den Heuvel and Sporns, 2019). More specifically, if a hub bears the throughput beyond its threshold, it could lead to catastrophic system collapse. For the whole brain level, network analysis showed decreased global and local connectivity efficiency in the ipsilesional hemisphere (Fekonja et al., 2021); Jutten et al. (2021) investigated fractional anisotropy (FA) of dMRI and achieved similar results of structural connectivity and systemic microstructural WM dis-integrity, especially in the IDH-wild-type (IDH-wt) group due to the worse biological attributes (Eckel-Passow et al., 2015). In addition, patients with IDH-wt and a greater lesion volume have significantly lower connectivity efficiency and cognitive capability (Kesler et al., 2017). In other words, functional connectivity was traditionally considered based on brain anatomy and the structural connectivity. Functional aberrations in cognition or perception might be caused by impairment of structural integrity. Wei et al. (2021) quantified the global and regional connectome disruptions in individual glioblastoma patients and investigated the prognostic value of connectome disruptions and topological properties using a structural connectome approach based on diffusion MRI, and they demonstrated that glioblastoma patients had decreased segregation and increased integration, which could potentially provide a useful biomarker for patient stratification.

Not only has structural connectivity provided the orientations of neurons for formulating surgical plan and navigation (Abdullah et al., 2013; Henderson et al., 2020), but it also could be advantageous for prediction of molecular characteristics and prognosis evaluation. On the one hand, determining the infiltration or displacement of cortical mapping and WM tracts would be significant for the extent of tumor removal and normal brain tissue protection (Castellano et al., 2012; Abhinav et al., 2015). Gliomas are the most common intra-axial brain tumors characterized by invasion into the surrounding WM tracts; thus, the combination of tractography with the operative navigation could show the adjacency between bulks and edemas of tumors and WM skeletons. Different approaches have been presented for mapping brain connectivity, and defining nodes and edges at each scale, so as to protect the integrity of the brain network. On the other hand, numerous studies have reported that gliomas have specific molecular parameters, such as IDH-mutation (Chen et al., 2017; Kesler et al., 2019) and MGMT promoter methylation (Chen et al., 2017). To date, deep neural networks have been utilized to extract dMRI metrics like ADC and FA value, and to assist IDH, MGMT, and other molecule subtype classification and prognosis prediction of glioma patients (Aliotta et al., 2019; Z. Huang et al., 2021; Yan et al., 2021).

Functional Connectome and the Alternations in Gliomas

Functional Connectome

Oriented in 1990s, fMRI measuring hemodynamic changes after enhanced neuronal activity has been applied on neuroscience research and clinical practice (Ogawa and Lee, 1990; Logothetis, 2008). When one stays awake and produces thinking activities, his/her related brain area would produce electrophysiological signals, which could be detected indirectly by changes in the oxyhemoglobin (HbO2) and deoxyhemoglobin (Logothetis and Wandell, 2004). Spatiotemporal properties of blood oxygenation level-dependent (BOLD) fMRI showed higher intensity signals when electrical activities happened due to higher HbO2 levels. Although the imbalance of supply and consumption of HbO2 remained elusive (Drew, 2019), one theory indicated that dilation of arteries increases neuronal activities, so as to produce higher HbO2 levels (Drew et al., 2011). Previous studies have proposed several models (Malonek and Grinvald, 1996) to simulate the relationships between BOLD signals and hemodynamic response function (HRF) (Cohen, 1997; Friston et al., 1998) and revealed linear relationships. Similar to the dMRI processing pipeline, fMRI data were initially preprocessed and HRF was then applied to deconvolute and fit the time series signals by linear regressive models. Recently, surface-based processing (Fischl et al., 1999) has become a popular approach to overcome the calculation and atlas problems, increasing the specificity of cortical activation patterns (Oosterhof et al., 2011; Saad and Reynolds, 2012; Brodoehl et al., 2020). Graph theory and other methods [independent component analysis (ICA), seed-based analysis, etc.], which are compatible to dMRI analysis theories, were also implemented according to the experiment design. Finally, a brain functional connectome was established to evaluate regional brain activity and functional connectivity (Figure 1).

To better evaluate a particular function parcellation in the brain area, task-fMRI is typically characterized more precisely to delineate the relationships between tasks and electrophysiological responses. Different cognitive tasks, including executive functioning, emotional processing, language, thinking, and implicit learning tasks, could be implemented by block, event-related, or mixed experiment design (Mwansisya et al., 2017). In clinical practice, however, doing task-fMRI needs more complex stimulation and recording systems, which might be difficult for tolerance of tumor patients. Rs-fMRI has also concentrated on brain functional mapping with subjects staying awake, while not performing any task (van den Heuvel and Hulshoff Pol, 2010). Studies have reported the resting-state networks, describing functional connectivity of brain regions in the resting state. The networks include the motor network, the visual and auditory network, dorsal and ventral attention network, central executive network (CEN), default mode network (DMN), and salience network (SN) (Spreng et al., 2013; Raichle, 2015; Yeshurun et al., 2021). Among these networks, CEN, DMN, and SN are the core neurocognitive-related networks. CEN, also known as frontoparietal network, contains dorsolateral prefrontal cortex and posterior parietal cortex, and is widely co-activated in cognitive behaviors and is involved in working memory (Menon, 2011). DMN, including the posterior medial cortex, medial prefrontal cortex, and temporoparietal junction, has shown a higher connectivity during resting state, and it was recognized as an active and dynamic intrinsic system to integrate information (Yeshurun et al., 2021). SN refers to cortical hubs in the dorsal anterior cingulate cortex, insular cortex, temporal pole, and amygdala, and it is related to links between the stimulus from external environment and the response of other networks dynamically (Seeley, 2019). A triple network model, consisting of CEN, DMN, and SN, has been proposed to understand the cross-network interactions and human brain state transition (Menon, 2011). Studies have reported that these non-traditional networks play key roles in the human brain in different diseases including Alzheimer’s disease and mental disorders (Menon, 2011; Liston et al., 2014; Satpute and Lindquist, 2019; Misiura et al., 2020; Martin-Subero et al., 2021).

Task-fMRI and rs-fMRI provide a non-invasive method to evaluate the active regions in the brain and their relationships. The spatial and temporal resolution, as well as the difference of HRF, limit the fMRI applications in causal inference; however, it is reliable to measure the fMRI for hypothesis testing (Olszowy et al., 2019). Recently, the molecular mechanism of functional connectivity was investigated by the correlation analysis between RNA-sequence and functional connectivity strength of rs-fMRI (Zhu et al., 2021), providing a potential to explore the genetic nature of the functional connectome.

Functional Connectome Alternations in Gliomas

Functional connectivity in gliomas was altered by the aggressiveness of glioma cells and the increased oxygenated blood flow in tumor vasculature (Warmuth et al., 2003; Hadjiabadi et al., 2018). Incorporating the assessment of the rs-fMRI network into the traditional glioma mapping approaches could be a versatile method for anatomical and functional localization and evaluation of molecular parameters (Ghinda et al., 2018; Stoecklein et al., 2020). Prior studies have reported a higher mean connectivity (Otten et al., 2012) and brain-wide disrupted networks associated with tumor-related remodeling of the neurovasculature (Hadjiabadi et al., 2018). In particular, several research groups measured the DMN in glioma by rs-fMRI and observed the increased integration of DMN in hippocampal areas, but the decreased integration in prefrontal and posterior cingulate cortex regions, and a trend of decreasing global connectivity with a higher WHO grade (Esposito et al., 2012; Harris et al., 2014; Maniar et al., 2021). The decreased DMN integration and deactivation may lead to cognitive decline of patients. SN connectivity in glioma could be reduced with an increased amplitude of low-frequency fluctuations (ALFFs), reflecting intrinsic brain activity, as well as promoting specific brain circuits to participate in cognitive tasks. The increased ALFF in the tumor contralateral regions may be explained by functional compensation due to neuronal plasticity and neuro-vascular uncoupling of tumor ipsilateral regions (Yang et al., 2021). Thus, tumors that are close to hubs of these networks could lead to dysfunction, but without obvious clinical manifestations due to a compensatory response of other areas. For traditional functional networks such as language and sensorimotor networks, studies through seed-based and ICA approaches revealed significantly reduced integrity and connectivity (Briganti et al., 2012; Niu et al., 2014; Mallela et al., 2016; Vassal et al., 2017; Yuan B. et al., 2020). Specifically, the reduction had a distinct pattern modulated by tumor position. The closer tumor is to hub regions, the more severe the effects to the function. These findings were consistent with the clinical manifestations in patients whose tumors were located in language, sensory, or motor regions.

Functional network also has the small-world attribute, and studies have reported the disturbed small-world manner and the decreased connective efficiency in LGG patients (Xu et al., 2013; Huang et al., 2014), who mainly suffered from cognitive deficits, memory reduction, and psychomotor dysfunction. However, an increased local efficiency and a higher local clustering coefficient were observed in the lesion areas, suggesting the potential preservation mechanism in small-world topology (Huang et al., 2014; Park et al., 2016). The specific biological behaviors of breaking and protecting small-world properties have not been fully explored, and a large number of studies have reported that IDH-wt was highly associated with the dysfunction of functional network (Derks et al., 2019; Jutten et al., 2020), indicating that infiltration and aggression of genetic characteristics of gliomas may affect the functional network. Recently, Mandal et al. (2020) established a model that combined cellular types and transcription of genetic drivers of glioma genesis with graph theory-based functional hub-based analysis, and found patterns of glioma locations. Their model explained over half of the variance in glioma location frequency, brain regions populated with putative cells of origin for glioma, neural stem cells, and oligodendrocyte precursor cells. On the other hand, the study also explained that the occupation and destruction of brain functional hubs by glioma decreased global small-world clustering. The complexity of biological characteristics of glioma and distribution of their locations are still accompanied by some problems in functional connectomes, indicating the necessity to investigate the correlations between tumor biology and brain function.

Mapping the functional connectome also provides sufficient data for surgical planning and prognosis assessment. In clinical practice, awake mapping using direct electrical stimulation (DES) of the brain can be achieved in surgery for gliomas in functional regions (Duffau, 2015). Anatomo-functional location by dMRI and rs-MRI before surgery could assist DES for precise functional protection. Qiu et al. (2017) demonstrated that using intraoperative rs-fMRI could directly localize the functional areas intraoperatively and avoid the risk of intraoperative seizures due to direct cortical stimulation. Postoperatively, patients mainly prefer to receive radiotherapy and chemotherapy. Thus, functional connectome analysis could evaluate the effects of treatment, brain recovery, and quality of life in patients. However, conflicts still exist between the tumor resection and brain network protection, and these traditional and non-traditional networks were invaded or pushed by tumor bulks, resulting in different operative approaches and tumor resection ranges. In future research, the balance between connectome protection and surgery and radio-chemotherapy will be deeply investigated.

Discussion

In 2009, the National Institutes of Health proposed the Human Connectome Project (HCP) with a primary goal of delineating the typical patterns of structural and functional connectivity in the brains of healthy individuals and patients with different psychiatric and neurologic diseases (Van Essen et al., 2012, 2013; Barch et al., 2013). To date, diseases, such as Alzheimer’s disease (Daianu et al., 2015; Sun et al., 2020), epilepsy (Hwang et al., 2020), and mental disorders (Korgaonkar et al., 2020), have been profoundly assessed in alternations of connectome. Glioma-related connectome projects have not been presented, although a large number of studies have illustrated the reduction in structural and functional connectivity as we mentioned. From traditional views, we paid further attention to tumor resection and identification of molecular characteristics. A series of pertinent brain macro perspectives and multi-center projects still need to be conducted to extend the knowledge of glioma, although connections between structural and functional connectome have remained obscure due to limitations of MRI technology. For example, in clinical practice, dMRI and fMRI are not always acquired due to the financial status of patients and or scanner machine-time shortage especially in developing countries. On the other hand, differences in sequence parameters would also make it difficult to achieve unbiased data standardization. From the patients’ view, many patients with glioma have alternations in cognition, motion, and attention and high risks of secondary epilepsy, so they may not cooperate with physicians with regard to scanning for a long time (Silva et al., 2018). Technologically, dMRI and fMRI techniques are limited by the resolutions. Lower signal-to-noise ratio (SNR) and resolutions may reflect poor image qualities (Jha et al., 2020; Morales, 2021; Turesky et al., 2021). In particular, for dMRI, the partial volume effect and the disability for non-Gaussian diffusion (Assaf and Pasternak, 2008) are the main difficulties that lead to inaccurate fiber tracking. For fMRI, excessive blood supply in glioma could uncouple the BOLD signal and susceptibility motion artifact remains an important concern for task design and imaging time (Silva et al., 2018; Morales, 2021). Besides, structural–functional coupling in glioma will be explored deeply in the future. Structural–functional coupling describes anatomical support for functional communications across brain regions (Baum et al., 2020; Gu et al., 2021). Zhang et al. (2021) reported the structural–functional connectome of language in glioma-induced aphasia, and patients with aphasia were without sufficient functional compensation in the supplementary motor area, which was mainly involved in aphasia. Traditionally, correlation analysis was commonly used for the structural connection (SC) and functional connection (FC) strength to obtain a structural–functional coupling matrix. However, this method still has the disadvantage of FC automatically forming stronger in the absence of a direct structural link (Honey et al., 2009); Sarwar et al. (2021) trained a deep learning network to predict the FC based on SC as input signals, and correlations were computed between pairs of subjects or between distinct modalities within the same individual to assess the structural–functional coupling strength. Hybrid connectivity ICA, which could extract the connectivity traits from SC and FC, was also applied to explore structural–functional coupling patterns (Amico and Goni, 2018). To improve the image qualities, 7-T or higher field MRI would be significant for providing higher-resolution and lower-noise images, containing more precise information for brain’s structure and function. At the same time, generative adversarial networks were used to generate and standardize the MRI data from different machines and to improve the image resolutions; thus, researchers can share and use the data stored in the cloud and perform large-scale clinical studies (Kazuhiro et al., 2018).

With the improvement of the brain–computer interface (BCI) that uses algorithms to decode brain signals, non-invasive brain signals, which could be obtained by converging multimodal MRI data to electrophysiological signals, could elevate the accuracy and efficiency of decoding. During the clinical management of glioma, aphasia or motion disability may occur due to tumor invasion or surgery. At present, electroencephalography (EEG)-based BCI has been widely used to decode brain signals via encode–decode deep learning models for patients with language or motor dysfunction (Makin et al., 2020; Li et al., 2021). Structural and functional connectome combined with EEG will provide more sufficient and precise data, which are beneficial for postoperative recovery and protection. Besides, BCI can improve cancer patients’ quality of life.

The meta-network analysis has been proposed, underlying the uniquely human propensity, to learn complex abilities and explain how post-lesional reshaping can lead to some degrees of functional compensation in patients with brain injury. Herbet and Duffau (2020) reviewed the progress in functional anatomy of the human brain and defined a meta-network that could integrate high-level cerebral functions from dynamic spatiotemporal dimensions, so as to provide structural and functional reshaping. Different from the existing approaches, the meta-network considers the systematic interactions rather than special brain regions, and it is more robust to noise interference. Adaptive changes of large-scale connectivity, causing cognition, perception, and emotion, were analyzed with the meta-network to better understand the underlying and constant connection patterns and their stability. Duffau (2021) also concentrated on applying the meta-network to glioma topography and surgery. For better tumor resection and brain protection, the interactive loop between tumor invasion, brain functional compensation, and the therapeutic strategy should be established. When the tumor cells involve the critical pathways, development of dynamic multimodal imaging methods may facilitate disconnecting the tumor-neural signaling and preserve the normal connectome. In order to treat glioma more effectively, multiple views and analysis approaches must correspond from molecular parameters to fiber orientations to global brain functional connectivity. However, for MRI data analysis, Bowring et al. (2019) explored the different tools with the same data and pipeline in fMRI and found inconsistent results. The results concluded that the effort would be strengthened for analysis scripts and statistical models. To enable the confirmation of these tools, deep learning algorithms could help to analyze the scripts of different tools, and construct more robust statistical models. Apart from dMRI and fMRI data themselves, multi-omics analysis provides a novel and detailed approach to glioma (Figure 2). Recent studies have investigated using the multi-omics sources including genomics, transcriptome, and proteomics for the glioma classification, prognosis prediction, and identification of therapeutic targets via machine learning in open datasets (Kamoun et al., 2016; Biswas and Chakrabarti, 2020; Zeng et al., 2021). Our studies have used radiomics in MRI to predict the key molecular status of glioma by deep learning networks (Chang et al., 2018; Zhang et al., 2020). However, there are still lack of systematic studies to combine the connectome and other scales of information of glioma. Thus, integrating data from different sources, and their processing with AI, together with validated biological and clinical findings will be significant for the development of a personalized surgical strategy.

FIGURE 2

Glioma-related multi-omics. Combination of omics including genomics, transcriptome, proteomics, radiomics, connectomics, etc. from microscopy to macroscopy and application of AI provided a novel and detailed approach for glioma, accompanied by a potential for BCI. AI: artificial intelligence, BCI: brain–computer interface.

In conclusion, from a personality treatment perspective, to reliably diagnose and manage gliomas, high-field imaging and data processing algorithms enable the achievement of high-resolution images and performing real-time analysis of multimodal MRI data, so as to construct the structural–functional connectome. Additionally, further systemic and multi-omics studies are required to fully understand the mechanisms of alternations in gliomas. In addition to the changes in connectome, the biological and molecular characteristics should be explored. Thus, additional efforts need to be dedicated to illustrate the nature of gliomas and to develop more effective treatments for gliomas.

Author Contributions

XL designed the study. ZC performed the research and wrote the manuscript. NY and CT helped to revise the manuscript. All authors have made a substantial, direct, and intellectual contribution to the work, and approved it for publication.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.