- 1Department of Clinical Neurosciences, Neurosurgery Service and Gamma Knife Center, Centre Hospitalier Universitaire Vaudois (CHUV), Lausanne, Switzerland
- 2Connectomics Laboratory, Department of Radiology, Centre Hospitalier Universitaire Vaudois (CHUV), Lausanne, Switzerland
- 3Institute of Bioengineering, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- 4Department of Radiology and Medical Informatics, University of Geneva (UNIGE), Geneva, Switzerland
- 5Stereotactic and Functional Neurosurgery Service and Gamma Knife Unit, Assistance Publique-Hôpitaux de Marseille, Centre Hospitalier Universitaire de la Timone, Marseille, France
- 6Neurology Department, Assistance Publique-Hôpitaux de Marseille, Centre Hospitalier Universitaire de la Timone, Marseille, France
- 7Department of Diagnostic and Interventional Neuroradiology, Centre de Résonance Magnétique Biologique et Médicale, Assistance Publique-Hôpitaux de Marseille, Centre Hospitalier Universitaire de la Timone, Marseille, France
- 8Faculty of Biology and Medicine (FBM), University of Lausanne (UNIL), Lausanne, Switzerland
- 9Signal Processing Laboratory (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
Essential tremor (ET) is the most common movement disorder. Its pathophysiology is only partially understood. Here, we leveraged graph theoretical analysis on structural covariance patterns quantified from morphometric estimates for cortical thickness, surface area, and mean curvature in patients with ET before and one year after (to account for delayed clinical effect) ventro-intermediate nucleus (Vim) stereotactic radiosurgical thalamotomy. We further contrasted the observed patterns with those from matched healthy controls (HCs). Significant group differences at the level of individual morphometric properties were specific to mean curvature and the post-/pre-thalamotomy contrast, evidencing brain plasticity at the level of the targeted left thalamus, and of low-level visual, high-level visuospatial and attentional areas implicated in the dorsal visual stream. The introduction of cross-correlational analysis across pairs of morphometric properties strengthened the presence of dorsal visual stream readjustments following thalamotomy, as cortical thickness in the right lingual gyrus, bilateral rostral middle frontal gyrus, and left pre-central gyrus was interrelated with mean curvature in the rest of the brain. Overall, our results position mean curvature as the most relevant morphometric feature to understand brain plasticity in drug-resistant ET patients following Vim thalamotomy. They also highlight the importance of examining not only individual features, but also their interactions, to gain insight into the routes of recovery following intervention.
Introduction
Essential tremor (ET) is the most common movement disorder, affecting up to 5% of individuals above 65 years of age (Louis and Ferreira, 2010). Patients exhibit postural and kinetic tremor of the hands and arms, sometimes with head, legs, or voice tremor as well (Chunling and Zheng, 2016). Sensory deficiencies, cognitive deficits, psychiatric and sleep disorders can also complement motor symptoms (Chandran and Pal, 2012; Jhunjhunwala and Pal, 2014; Louis, 2016; Jiménez-Jiménez et al., 2020, 2021b). While ET has clear underlying genetic origins given the frequent occurrence of positive family history (Haubenberger and Hallett, 2018), the reliable identification of culprit genes remains partly unconclusive (Kuhlenbäumer et al., 2014; Tio and Tan, 2016; Deng et al., 2019; Siokas et al., 2020; Jiménez-Jiménez et al., 2021a).
Neuroimaging studies have enabled to localize the brain regions implicated in the motor component of ET, which belong to the so-called tremor network (Raethjen and Deuschl, 2012; Hallett, 2014; Sharifi et al., 2014): they include the cerebellum [sometimes regarded as the cornerstone of ET (Benito-León and Labiano-Fontcuberta, 2016; Ibrahim et al., 2020)], the motor thalamus (the ventro-intermediate nucleus, or Vim) and the motor cortex. While the exact pathophysiological mechanisms at play remain debated within several partly overlapping theories (Deuschl and Elble, 2009; Benito-Leon, 2014; Gironell, 2014), it is believed that ET manifests itself as a dysregulated network of interacting areas.
This makes the study of brain structure and function at the level of individual networks (to specifically address motor or non-motor impairments caused by ET, for instance), or at the whole-brain scale (to characterize cross-regional interactions in their entirety), a sensible analytical direction to pursue. Graph theory has emerged as the primary analytical approach for this purpose, since it provides an elegant and powerful way to gain insight into how information flows in complex systems described by edges linking the vertices of a graph.
Application areas are numerous: they include engineering problems such as the study of transportation systems (Derrible and Kennedy, 2011) or gear transmission (Xue et al., 2016), as well as life science disciplines such as proteomics (Grindrod and Kibble, 2004) or molecular topology (Amigó et al., 2009). Neuroscience has also been a particularly fruitful field for the application of graph theory (Bassett and Sporns, 2017; Sporns, 2018; Farahani et al., 2019); this is partly because a graph denoting the interplays between different brain regions can be meaningfully constructed from various imaging modalities, two notable examples of which are diffusion-weighted magnetic resonance imaging (DW-MRI, where edges are the physical connections between the areas) and resting-state functional magnetic resonance imaging (RS-fMRI, where edges denote the statistical dependencies between the activity time courses of the areas) (Bullmore and Sporns, 2009). More recent developments also include the generation of multilayer networks [e.g., interconnected temporally subsequent graphs capturing the dynamics of functional brain interactions (Preti et al., 2017; Pedersen et al., 2018)], and the subfields of graph signal processing, graph neural networks and graph learning (see Huang et al., 2018; Li R. et al., 2021 for recent reviews).
Graph theory has been pivotal in better grasping the neural underpinnings of cognition in the healthy and in the diseased brain (Minati et al., 2013; Hallquist and Hillary, 2018). ET has also been contemplated from this perspective: using DW-MRI and focusing on the executive network, Prasad et al. (2020) revealed lower global and local efficiencies of frontal executive brain centers and of the anterior cingulate cortex compared to matched healthy controls. In a conceptually similar study centered on the motor network, the local efficiency of the cerebellum was higher in ET patients, and those with resting tremor also showed further increases in the thalamus, globus pallidus, caudate, and supplementary motor area (Caligiuri et al., 2017).
Using RS-fMRI, Li J. Y. et al. (2021) evidenced lower nodal efficiency in frontal and supplementary motor areas, the precuneus and the cerebellum. Individuals with concomitant depression showed further pre-central, post-central, and frontal decreases. In another investigation, widespread differences that also involved other graph theoretical metrics (betweenness centrality, degree) were revealed in frontal, occipital, temporal, and cingulate cortices as well as subcortical and cerebellar loci (Benito-León et al., 2019).
From these reports, it transpires that ET-induced alterations of brain structure and function span multiple brain networks. In addition, these features of the disease can be captured by various graph theoretical measures. For these reasons, in the present work, we opted for studying ET through a whole-brain approach, and the combined extraction of a set of complementary graph theoretical metrics.
We propose three important advances compared to previous work: first, we study a well-defined population of drug-resistant patients with ET, who underwent Gamma Knife (GK) stereotactic radiosurgery of the Vim (Elaimy et al., 2010). We explore not only the differences between these subjects and matched healthy controls (HCs), but also brain plasticity 1 year after Vim thalamotomy. Only one past study has probed such mechanisms using graph theory (Jang et al., 2016), but the surgical intervention differed (imaging-guided high-intensity focused ultrasound thermal ablation). Moreover, the number of subjects was low (N = 10), and the analyses were centered on the motor network (with decreases in degree and efficiency following thalamotomy). Of note, the follow-up was interrupted 3 months after the intervention.
Second, we construct a graph neither from DW-MRI nor from RS-fMRI data as in the above cases, but from structural MR images. To do so, we leverage structural covariance analysis (SCA) (see Alexander-Bloch et al., 2013; Evans, 2013 for reviews), where the extent to which a morphometric measure of interest (e.g., cortical thickness) correlates across subjects in pairs of regions is used to build the graph. Such patterns of covariance are characteristic of the human cortex (Mechelli et al., 2005), are under genetic control (Schmitt et al., 2008; Romero-Garcia et al., 2018; Morgan et al., 2019), and recapitulate structural connectivity features (Yee et al., 2018). Here, we study three complementary features whose potential has extensively been validated (Bassett et al., 2008; Chen et al., 2008; He et al., 2008; Bernhardt et al., 2011; Bethlehem et al., 2017; Khundrakpam et al., 2017): cortical thickness (CT), surface area (SA), and mean curvature (MC). They encode partly unique information in the healthy brain (Chiarello et al., 2016), and undergo distinct environmental modulations (Kelly et al., 2013; Besteher et al., 2017), which justifies their parallel assessment.
To date, only one other work has applied graph theory to structural MR images in ET: Yang et al. (2021) quantified the similarity of gray matter profiles (Kullback–Leibler divergence-based) within regions to construct subject-wise graphs, and revealed altered metrics in frontal, temporal and angular gyri, the caudate, hippocampus, thalamus, and some parts of the cerebellum.
Third, in addition to the parallel assessment of individual morphometric properties, we propose a novel analysis that enables to additionally study cross-property dependences through the generation and analysis of directional graphs. This is motivated by the acknowledged genetic and phenotypic complementarity of the measures (Sanabria-Diaz et al., 2010; Yang et al., 2016), whose interactions may be altered by ET and/or thalamotomy. We introduce simple telling features that can be generated from such graphs, and interpret them in the context of our dataset.
Materials and Methods
Participants
We considered uniform structural MRI data from 34 ET patients (both before thalamotomy and one year later) scanned on the same 3T Siemens Skyra MR machine, and 29 HCs. All patients were right-sided and suffered from drug-resistant right-dominant ET. All underwent unilateral left Vim thalamotomy by GK.
The Timone University Hospital Ethical Committee (ID-RCB: 2017-A01249-44) granted formal approval for this study (including by the Ethics Committee at national level, CNIL-MR-03), and individual consent was also obtained from all subjects. Patients were neurologically evaluated and referred by TW, a neurologist specialized in movement disorders. All patients had a clear diagnosis of ET and showed no other structural abnormalities on pretherapeutic MRI. Demographic characteristics of the ET patients and HC subjects can be found in Table 1; both groups were matched for age and gender.
Several measures were used to clinically evaluate ET patients, and their recovery after the thalamotomy: Activities of Daily Living (ADL) from the survey designed by Bain et al. (1993), Tremor Score on Treated Hand (TSTH) from the Fahn–Tolosa–Marín rating scale (Fahn et al., 1993), head tremor (Tremor Research Group Essential Tremor Rating Assessment, from 0 to 3), and Quality of Life in Essential Tremor (QUEST) (Tröster et al., 2005). Clinical data is summarized in Table 1, where a significant improvement in clinical tremor scores can be observed across all the quantified measures upon thalamotomy.
Importantly, SCA is not compatible with the subject-wise investigation of these scores: indeed, only one measure of covariance is generated per group (e.g., post-therapeutic ET patients). Thus, one cannot conduct classical correlation analyses between morphometric features and clinical scores.
Imaging
T1-weighted images were acquired on a head-only 3T machine (SIEMENS SKYRA, Munich, Germany, 32-channel receive-only phase-array head coil), with the following parameters: TR/TE = 2300/2.98 ms, isotropic voxels of 1 mm3, 160 slices.
As medication was frequently ineffective, most patients no longer received treatment at the time of Vim thalamotomy by GK. Thus, scanning was performed in a drug-naïve state (drugs having been stopped at least 3 days prior to scanning).
Radiosurgical Procedure
Thalamotomy was performed with GK between September 2014 and April 2016, always at the Centre Hospitalier Universitaire de la Timone in Marseille. The surgeon in charge was JR, who used the Leksell Gamma Knife and associated Leksell GammaPlan software (Elekta Instruments, AB, Stockholm, Sweden). To avoid artifacts, DTI data was first acquired without the frame, and then co-registered with the therapeutic stereotactic images. The Leksell coordinate G frame (Elekta Instruments, AB, Stockholm, Sweden) was always applied under local anesthesia on the day of the thalamotomy. After positioning the frame, patients underwent both stereotactic CT and MRI.
Landmarks of interest, including the anterior and posterior commissures, were identified on an MR scan (particularly on T2 CISS/FIESTA sequence, replacing a former ventriculography). Uniform indirect targeting was performed using the Guiot diagram (Tuleasca et al., 2017), placed 2.5 mm above the anterior-posterior commissure line, and 11 mm lateral to the wall of the third ventricle. A single 4 mm isocenter was always used, and a maximum prescription dose of 130 Gy at the 100% isodose line was uniformly prescribed (Tuleasca et al., 2017).
Computation of Morphometric Properties
The Freesurfer software (Fischl, 2012) was used to extract three morphometric measures of interest from structural MR images for a set of PC = 68 cortical regions: CT, SA, and MC. Briefly, after linear registration to MNI space and bias field removal, the image at hand is skull-stripped (Ségonne et al., 2004), and voxels are classified as belonging to white matter or to another tissue category on the basis of their intensity and direct neighborhood. Hemispheres are separated, cerebellum and subcortex are removed, and the interface between the white and gray matters is located. From there, the pial surface is also tiled, and local estimates of CT, SA, and MC can be extracted (Fischl and Dale, 2000). Further details can be found in Dale et al. (1999) and Fischl et al. (1999). Eventually, local estimates are converted into PC regional values per morphometric measure, using the Desikan-Killiany atlas (Desikan et al., 2006). These were complemented by measures of regional volume for PNC = 19 non-cortical areas (including the brainstem, subcortical nuclei, and cerebellum), for a total of P = 87 parcels. They are summarized in Supplementary Table 1.
Obtained measurements were eventually linearly regressed out for age, gender, and total gray matter volume, separately within each group. The residuals were used for all subsequent analyses.
Analysis of Individual Morphometric Measures
The process described below was identically conducted for each of the morphometric measures of interest. A schematic description of the undertaken steps, and of subsequent subparts of the analysis, is also provided in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/873605/fnagi-14-873605-HTML/image_m/fnagi-14-873605-g001.jpg)
Figure 1. Schematic description of the method. (A) We consider three groups of subjects: healthy controls (N = 29, depicted by circles), and drug-resistant ET patients before and 1 year after Vim thalamotomy (N = 34, respectively denoted by rectangles and stars). For each structural scan, following Freesurfer-based processing and atlasing into regions of interest, the PC = 68 cortical areas at hand can be described by their thickness (CT, orange), surface area (SA, purple), and mean curvature (MC, cyan). The PNC = 19 non-cortical areas are characterized by their volume (V, gray). (B) Similarly for each of the three groups of subjects, cross-regional structural covariance is computed across subjects for CT (top left panel), SA (top middle panel), and MC (top right panel). The extracted information can equivalently be summarized in a matrix, or in a graph, where negative-valued structural covariance is set to zero. Regional degree, clustering coefficient (CC), and eigenvector centrality (EC) are computed for graph densities ranging from 20 to 60%, and the area under the curve (AUC) is taken as a regional output measure of interest for each morphometric property and graph metric case. In a similar fashion, cross-property covariance can be computed for the CT/SA (bottom left panel), SA/MC (bottom middle panel), and CT/MC (bottom right panel) cases. The obtained graphs are then directional, and in-degree and out-degree can be computed across graph densities to generate AUC output measures. (C) The computations described in panel (B) are performed similarly for the HC, ETpre, and ETpost groups. Note that cross-property metrics are only available for the PC cortical regions. There is a total of 15 separate subcases with output AUC values: 3 morphometric properties (CT, SA, and MC) for 3 graph theoretical metrics (degree, CC, and EC), plus 3 cross-property pairs for 2 graph theoretical metrics (in-degree and out-degree). (D) Similarly for each of the 15 subcases at hand, the regional difference in AUC can be computed between the HC and ETpre groups, or between the ETpost and ETpre ones. The process described in panel (B) is then rerun n = 8’000 times after randomly shuffling the subjects across groups, to generate a null distribution of AUC differences. The actual value (vertical dashed line) is eventually compared to this null distribution to assess significance, with proper correction for the number of examined regions and subcases in parallel.
Separately for the HC, pre-thalamotomy (abbreviated ETpre from there onward) and post-thalamotomy (abbreviated ETpost) data, Pearson’s correlation coefficient R was computed for each pair of regions. A positive/negative value means that when the measure in the first region is larger in one subject, it tends to be larger/lower in the second region.
To enable graph theoretical analysis, negative-valued edges were excluded (set to zero). The percentages of retained edges were: for CT, 79.6311% (HC group), 72.0663% (ETpre), 67.9497% (ETpost); for SA, 67.3617% (HC), 69.1794% (HCpre), 71.8524% (ETpost); for MC, 85.298% (HC), 77.2521% (ETpre), 81.5825% (ETpost).
We computed graph theoretical measures from each structural covariance matrix to quantify regional properties of morphometric dependences with the rest of the brain. Here, we considered three different graph measures that convey complementary information: degree, clustering coefficient, and eigenvector centrality. Below, we briefly detail each, but more details can be found in Rubinov and Sporns (2010).
A group-wise structural covariance matrix can equivalently be seen as a graph G, characterized by an adjacency matrix A. In what follows, we generated a graph at a predefined density ρ (see below for more details about this parameter). We made sure that any graph analyzed therein was always fully connected.
Nodal degree for a region p is simply the sum of edge weights linked to the node at hand: . Larger values denote an overall more influential node within the network. The clustering coefficient is a measure of how much the direct neighborhood of a node is interconnected. The eigenvector centrality measures to what extent the node of interest takes part in modular communities.
To assess the differences between HC subjects and ET patients, we computed the difference HC – ETpre for each regional graph metric. Similarly, to investigate the effects of thalamotomy, we computed the difference ETpost – ETpre.
Because the graph density at which to analyze the results is a priori unknown, we set to instead consider the area under the curve (AUC), or the sum of differences across densities ranging from 20% (the lowest value at which all graphs were fully connected) to 60% (for larger values, there would not be enough positive-valued SC elements).
For statistical assessment, these differences must be compared to an appropriate null distribution. To do so, we resorted to non-parametric permutation testing, by recomputing graph metrics over 8’000 null realizations for which subjects were randomly shuffled across groups.
False discovery rate (FDR)-corrected p-values were obtained and analyzed for each contrast of interest. When presenting our results, we considered two significance level α1 = 0.01 and α2 = 0.001.
Cross-Property Analysis
In addition to traditional graph theoretical analysis, we sought to introduce another graph-centered approach to enable the investigation of cross-regional statistical dependences across morphometric modalities.
Consider the data from two morphometric properties, contained in M1 and M2 (each of size P x S, with S the number of subjects in the group at hand); the classical structural covariance equation would yield symmetrical matrices of size P x P: SC1 = cov(M1, M1)/σ12 and SC2 = cov(M2, M2)/σ22, with σi the standard deviation for property i.
Instead, we compute the cross-correlation SC1,2 = cov(M1, M2)/(σ1 σ2), which also has size P x P, but is not symmetrical anymore. In this matrix, if element (i,j) is positive/negative, then when modality 1 in region i is larger in a given subject, modality 2 in region j will tend to be larger/lower. However, if element (j,i) is positive/negative, then when modality 1 in region j is larger in a given subject, then modality 2 in region i will be larger/lower.
Let CT, SA, and MC be respectively denoted by the indices 1, 2, and 3; we can thus compute three cross-correlation matrices: SC1,2, SC1,3, and SC2,3. Each can equivalently be seen as a directional graph (G1,2, G1,3, and G2,3), or a non-symmetrical adjacency matrix (A1,2, A1,3, and A2,3). For simplicity, we set null diagonal entries to exclude self-loops, and exclude negative-valued edges.
Here, with inspiration from time-resolved functional neuroimaging work on directional graphs (Bolton et al., 2020), we propose to compute two simple measures from each graph: (1) the in-degree, and (2) the out-degree. The in-degree for modality pair (m1, m2) at region j is defined as , and the out-degree as . The in-degree is larger if modality m2 in region j is more strongly dependent on the values of modality m1 in the other brain areas. The out-degree is larger if modality m1 in region j more strongly influences the values of modality m2 in the other brain areas.
We resort to the same non-parametric approach as above for statistical significance assessment. In more details, for both metrics and for each of the three pair-wise modality cases, we compute the HC – ETpre and ETpost – ETpre group differences and compare them to a null distribution generated following random shuffling of the subjects across groups (n = 8’000 permutations).
Availability of the Data and Scripts
All the analytical steps described above were performed with custom scripts and MATLAB2014b (MathWorks, Natick, MA, United States). Colormaps for plotting were generated with the cbrewer toolbox.1 Investigated graph measures were computed using the Brain Connectivity Toolbox (Rubinov and Sporns, 2010).
Data sharing is not applicable to this article as no new data were created or analyzed in this study. All the scripts used in this work are freely available at https://github.com/TiBiUan/SCA_GraphTheoretical.git.
Results
The results of the analyses centered on CT are presented in Supplementary Figure 1. None of the probed graph theoretical metrics yielded any significant outcome following FDR correction. For degree (Supplementary Figure 1A) and eigenvector centrality (Supplementary Figure 1C), smaller metric values in the ETpost than in the ETpre group in the left parahippocampal gyrus was nonetheless a noteworthy, shared feature, with the actual group difference falling within the bulk of null distribution outliers. For clustering coefficient (Supplementary Figure 1B), ETpre values were smaller than HC ones in the brainstem, but the smallest edge density cases were by far the largest contributors, as denoted by the dominance of the dark blue shades in the HC and ETpost stacked bars.
The results of surface area analyses are presented in Supplementary Figure 2. None of the investigated cases reached significance. As with cortical thickness, the left parahippocampal gyrus degree (Supplementary Figure 2A) and eigenvector centrality (Supplementary Figure 2C) values were nonetheless noteworthily larger in the ETpre group compared to the ETpost one. Clustering coefficient results (Supplementary Figure 2B) were overall very similar across all groups.
Figures 2A–C displays the results from mean curvature analyses (a more exhaustive version is also available as Supplementary Figure 3). For degree (Figure 2A) and the ETpost – ETpre contrast, values were significantly larger pre-intervention in the left parahippocampal gyrus (region 15, ΔAUC = −86.9766, p = 0.01), bilateral pericalcarine gyrus [regions 20 and 54, ΔAUC = −159.3861 and −178.5003 (left and right sides, respectively), p < 0.001 in both cases], and right cuneus (region 38, ΔAUC = −114.0338, p < 0.01). For clustering coefficient (Figure 2B), 6 regions reached significance when assessing the impacts of intervention: values were larger before thalamotomy in the bilateral superior parietal cortex [regions 28 and 62, ΔAUC = −1.1811 and −1.2203 (left and right sides, respectively), p < 0.01 in both cases], right pericalcarine gyrus (region 54, ΔAUC = −2.6884, p < 0.001), right precuneus (region 58, ΔAUC = −1.3199, p < 0.001), and left thalamus (region 71, ΔAUC = −1.6116, p < 0.01). Conversely, clustering coefficient increased after the intervention in the left insula (region 34, ΔAUC = 1.2759, p < 0.01). For eigenvector centrality (Figure 2C), the left supramarginal gyrus and left insula showed significantly larger values after the intervention (respectively regions 30 and 34, ΔAUC = 0.5126 and 0.7964, p < 0.001 in both cases).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/873605/fnagi-14-873605-HTML/image_m/fnagi-14-873605-g002.jpg)
Figure 2. Significant graph theoretical analysis findings. For regional degree (A), clustering coefficient (B), eigenvector centrality (C), and out-degree (D), group differences in area under the curve for the ETpost – ETpre contrast in the assessment of mean curvature structural covariance (A–C) or cortical thickness/mean curvature interactions (D). Actual values are denoted by black rectangles, and associated null distributions are reflected by box plots whose color coding matches the brain lobe at hand. The regions that reached significance in each case are highlighted by a light gray box and labeled. *: p < 0.01; **: p < 0.001.
Regarding the interaction between morphometric properties, there were no significant results involving surface area and any of the two other features (Supplementary Figures 4, 5). Regarding the cortical thickness/mean curvature case (Figure 2D; see also Supplementary Figure 6 for more exhaustive displays), the out-degree was significantly larger in the ETpost than ETpre group in the left pre-central gyrus (region 23, ΔAUC = 178, p < 0.01), and lower in the right lingual gyrus (region 46, ΔAUC = −185, p < 0.001) and bilateral rostral middle frontal gyrus [regions 26 and 60, ΔAUC = −163 and −177 (left and right sides, respectively), p < 0.001 in both cases].
Discussion
In the present work, we leveraged graph theoretical analysis to study structural covariance patterns in patients with ET before and after thalamotomy, as well as compared to age-matched HCs. In doing so, we considered three popular and complementary morphometric properties: cortical thickness, surface area, and mean curvature. All the significant results from our analyses are summarized in Table 2.
Interestingly, the comparison between pre-thalamotomy patients with ET and HCs did not yield any significant outcome. There are several potentially overlapping explanations to this negative finding. First, the assessment of structural covariance may be less useful than other approaches (such as DW-MRI- or RS-fMRI-based ones) to unravel the brain alterations induced by ET due to its less specific nature. Second, it could be that graph theoretical analysis is not the suited tool to unravel ET-related SC alterations; for instance, the finer-grained assessment of individual cross-regional interactions (as opposed to more global summarizing metrics as considered here) may be a more fruitful research avenue that should be explored in future work. Third, the anatomical underpinnings of ET may bear a subject-specific nature, and hence, may not be adequately captured by a group-level analysis relying on correlations across subjects as deployed here. Along this line of reasoning, it is interesting to notice that essential tremor is, nowadays, considered by many as a family of disparate diseases rather than a single entity (Jankovic, 2002; Elble, 2013; Espay et al., 2017; Louis, 2021).
One could hypothesize that thalamotomy, by renormalizing morphometric features, yields a more homogeneous group compared to HCs or ET patients before intervention. Accordingly, we found several significant impacts of Vim stereotactic radiosurgical thalamotomy. First, mean curvature was by far the most impacted morphometric property following the intervention. This implies that mechanisms of brain plasticity may occur more easily regarding this geometric aspect of the brain. Cortical gyrification is a genetically regulated developmental process (Papini et al., 2020), which is also specifically altered over other morphometric properties in schizophrenia (Schultz et al., 2010, 2013). Interestingly, mean curvature can be affected by experience alone: for example, gyrification weakens less rapidly upon aging in bilinguals (Del Maschio et al., 2019), and meditation practitioners show a greater extent of gyrification in several regions of the brain compared to matched non-practitioners (Luders et al., 2012). In a study following patients suffering from anorexia nervosa, Bernardoni et al. (2018) also found that cortical folding was broadly lowered, but that weight restoration within only 3 months sufficed to restore the values to a normal level. Taken together, these reports lend credit to the possibility that within one year after Vim thalamotomy, mean curvature could indeed have underwent significant changes in patients with ET, without as extensive modulations of other geometrical brain features.
Another factor that may have contributed to the presence of extensive differences in the ETpost – ETpre, but not the HC – ETpre contrast case, is the assessment of the same subjects in both groups for the former, but not the latter comparison. Indeed, our analyses did not explicitly model within-subject variance, because the resulting graphs then bear a different interpretation (see Supplementary Material for details). However, we explicitly verified that with a mixed model approach accounting for repeated measures, whole-brain ETpost – ETpre structural covariance difference patterns remain largely identical to those analyzed therein (see Supplementary Figure 7). The only morphometric property for which there were some noticeable impacts was SA, for which structural covariance differences took a narrower range when modeling within-subject variance; recall that these differences were already non-significant in our analyses. For MC, the match across methods was the largest (spatial correlation of whole-brain patterns: 0.97), leading us to rule out a major impact of repeated measures on our findings.
Turning to the areas involved in the observed differences, the left thalamic clustering coefficient was larger before the intervention (at our less stringent significance threshold of p < 0.01), revealing that one of its consequences was the loss of cross-regional dependences within its neighboring areas.
Collectively, our other findings point at important rearrangements revolving around the dorsal visual pathway, which is involved in visuospatial attention and action guidance (Freud et al., 2016). First, several low-level visual areas (bilateral pericalcarine gyrus, and to a lesser extent, also the right cuneus) showed a lower degree following thalamotomy.
Second, higher-level regions in the dorsal visual stream also exhibited mean curvature changes upon intervention: the left parahippocampal gyrus, linked to visuospatial processing (Aminoff et al., 2013), displayed lower degree at our more lenient significance threshold; the left supramarginal gyrus, involved in the planning of visually guided reaching and grasping movements (Andres et al., 2017; McDowell et al., 2018; Potok et al., 2019), exhibited a larger eigenvector centrality; in addition, clustering coefficient was significantly smaller (at p < 0.01) in the bilateral superior parietal cortex, a cornerstone region bridging visual and motor functions (Corbetta et al., 1995; Husain and Parashkev, 2007), as well as in the right precuneus, directly tied to visuospatial attention in past transcranial magnetic stimulation work (Mahayana et al., 2014) and known to jointly functionally relate to motor, visual and cognitive brain centers (Margulies et al., 2009).
Third, the eigenvector centrality of the left insula [classically known for its role in salience monitoring (Menon and Uddin, 2010), and directly associated to attentional performance (Varjačić et al., 2018)] also increased post-intervention.
Beyond the individual assessment of each morphometric property, our novel cross-feature analysis also enabled to evidence the presence of Vim thalamotomy-induced changes in the brain that modify the interplay between CT and MC. Out-degree was lower after intervention in the right lingual gyrus; in other words, this region’s cortical thickness is normally statistically related to overall mean curvature in the rest of the brain, and the dependence was lowered upon thalamotomy. A similar observation was made concerning the bilateral rostral middle frontal gyrus, which is a sub-component of the dorsolateral prefrontal cortex, an area involved in attentional functions (Kondo et al., 2004; Johnson et al., 2007). These two findings show that modulations at the level of the dorsal visual stream also extend to dependences between two different morphometric properties. Finally, out-degree increased (albeit with only mild significance) in the left pre-central gyrus, home of the motor cortex contralateral to the targeted side; this may reflect compensatory brain plasticity, where left pre-central cortical thickness gains new relationships to the mean curvature of other brain regions.
The presence of interactions across morphometric properties is unsurprising given the joint orchestration of their patterns upon brain development (Raznahan et al., 2011; Schnack et al., 2015); our study calls for further efforts to not only analyze a given morphometric feature independently, but to consider all of them jointly for a more complete understanding of the brain. In that, our cross-modality analysis strategy somewhat echoes recent work in which a handful of modalities are jointly analyzed to generate subject-wise morphometric similarity networks (Seidlitz et al., 2018). The advantage of this approach over ours is that measures can be estimated for each individual, enabling to then readily probe possible relationships to clinical scores. Its drawback may lie in the interpretability of the results; by instead solely addressing the dependence between two measures, mechanistic insight may be easier to achieve.
Our study bears some important limitations that should be kept in mind: in particular, the number of subjects remains quite low. This is largely because the analyzed data was collected at two successive time points with a one-year gap. Our results should thus be considered with caution, in the specific setting of the dataset that we analyzed. Further work on a larger pool of subjects will be required to determine whether the findings presented therein can generalize.
Another limitation is the use of a somewhat spatially gross atlas compared to state-of-the-art (sub)cortical alternatives (Schaefer et al., 2018; Tian et al., 2020). While this certainly precludes our ability to more finely capture ET-induced alterations and their evolution upon intervention, the number of statistical tests conducted in parallel would dramatically increase with a finer-grained parcelation, and so would the extent of dependences between spatially neighboring areas. On top of rendering interpretation much more complex, tailored statistical correction methods would then have to be deployed for appropriate analysis.
Partly linked to the above point is the fact that, while most pinpointed areas could be related to visual or attentional functions (13 of 16 or 8 of 9 as a function of the considered significance threshold), their relatively large size makes it possible that the observed thalamotomy-induced changes instead associate to other functions subserved by neighboring foci.
Aside from graph theoretical approaches, there exist many other ways by which structural covariance can be probed. Alongside the route undertaken here, an alternative option for future work could be to investigate in greater details the specificities of whole-brain structural covariance with specific regions of interest. For example, a seed-based structural covariance map can be computed in each group at stake, and measures such as their spatial properties can be compared (DuPre and Spreng, 2017). Another alternative is the use of Partial Least Squares analysis to generate seed-to-whole brain covariance maps to then contrast across groups of interest (Persson et al., 2014).
Yet another main analytical line involves the decomposition of the morphometric data into different summarizing components, as can be achieved using various approaches such as Independent Component Analysis (Hafkemeijer et al., 2016), Principal Component Analysis (Brickman et al., 2008; Bergfield et al., 2010), non-negative matrix factorization (Sotiras et al., 2015, 2017), or diffusion embedding (Masouleh et al., 2020). In this study, we decided not to conduct such analyses because they work optimally when the data has a voxel-wise spatial resolution, whereas we addressed regional properties.
Finally, other more recent approaches can also be envisaged: in particular, the person-based similarity index (PBSI) is a recent measure that quantifies how similar a given subject is to all the others in terms of its whole-brain morphometric profile (Doucet et al., 2019, 2020) and can be directly applied to regional data. In future work, it will be interesting to see whether such tools enable to go beyond the group-level characterization, to start gaining further insight into individual specificities, and what leads some subjects to recover from ET more effectively than others upon surgery.
Conclusion
In this work, we evidenced a significant impact of Vim thalamotomy on cortical gyrification in drug-resistant ET patients. Changes following the intervention could be captured, across complementary graph theoretical metrics, within visuospatial and attentional areas. This denotes the network-level reconfiguration of these regions’ statistical dependencies to the rest of the brain, and calls for more extensive investigations of the visual circuitry in future studies of ET.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Timone University Hospital Ethical Committee. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
TB wrote the scripts, performed the analyses, and wrote the original manuscript draft. JR, TW, and NG contributed to data acquisition. DV and ML surveyed the work, provided guidance regarding the analyses, and helped interpret the results. CT contributed to data processing as well as the analysis and interpretation of the results. All authors thoroughly reread the manuscript.
Funding
CT gratefully acknowledges the receipt of a “Young Researcher in Clinical Research Grant” (“Jeune Chercheur en Recherche Clinique”) from the University of Lausanne (UNIL), Faculty of Biology and Medicine (FBM), and the Lausanne University Hospital (CHUV).
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.
Acknowledgments
The authors wish to thank Meutia Parahita for her help in generating the illustrations presented in Figure 1.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.873605/full#supplementary-material
Footnotes
References
Alexander-Bloch, A., Giedd, J. N., and Bullmore, E. (2013). Imaging structural co-variance between human brain regions. Nat. Rev. Neurosci. 14, 322–336. doi: 10.1038/nrn3465
Amigó, J. M., Gálvez, J., and Villar, V. M. (2009). A review on molecular topology: applying graph theory to drug discovery and design. Naturwissenschaften 96, 749–761. doi: 10.1007/s00114-009-0536-7
Aminoff, E. M., Kveraga, K., and Bar, M. (2013). The role of the parahippocampal cortex in cognition. Trends Cogn. Sci. 17, 379–390. doi: 10.1016/j.tics.2013.06.009
Andres, M., Pelgrims, B., Olivier, E., and Vannuscorps, G. (2017). The left supramarginal gyrus contributes to finger positioning for object use: a neuronavigated transcranial magnetic stimulation study. Eur. J. Neurosci. 46, 2835–2843. doi: 10.1111/ejn.13763
Bain, P. G., Findley, L. J., Atchison, P., Behari, M., Vidailhet, M., Gresty, M., et al. (1993). Assessing tremor severity. J. Neurol. Neurosurg. Psychiatr. 56, 868–873.
Bassett, D. S., Bullmore, E., Verchinski, B. A., Mattay, V. S., Weinberger, D. R., and Meyer-Lindenberg, A. (2008). Hierarchical organization of human cortical networks in health and schizophrenia. J. Neurosci. 28, 9239–9248. doi: 10.1523/JNEUROSCI.1929-08.2008
Benito-Leon, J. (2014). Essential tremor: a neurodegenerative disease? Tremor Other Hyperkin. Movement. 4:252. doi: 10.5334/tohm.227
Benito-León, J., and Labiano-Fontcuberta, A. (2016). Linking essential tremor to the cerebellum: clinical evidence. Cerebellum 15, 253–262. doi: 10.1007/s12311-015-0741-1
Benito-León, J., Sanz-Morales, E., Melero, H., Louis, E. D., Romero, J. P., Rocon, E., et al. (2019). Graph theory analysis of resting−state functional magnetic resonance imaging in essential tremor. Hum. Brain Mapp. 40, 4686–4702. doi: 10.1002/hbm.24730
Bergfield, K. L., Hanson, K. D., Chen, K., Teipel, S. J., Hampel, H., Rapoport, S. I., et al. (2010). Age-related networks of regional covariance in MRI gray matter: reproducible multivariate patterns in healthy aging. Neuroimage 49, 1750–1759. doi: 10.1016/j.neuroimage.2009.09.051
Bernardoni, F., King, J. A., Geisler, D., Birkenstock, J., Tam, F. I., Weidner, K., et al. (2018). Nutritional status affects cortical folding: lessons learned from anorexia nervosa. Biol. Psychiatr. 84, 692–701. doi: 10.1016/j.biopsych.2018.05.008
Bernhardt, B. C., Chen, Z., He, Y., Evans, A. C., and Bernasconi, N. (2011). Graph-theoretical analysis reveals disrupted small-world organization of cortical thickness correlation networks in temporal lobe epilepsy. Cereb. Cortex 21, 2147–2157. doi: 10.1093/cercor/bhq291
Besteher, B., Gaser, C., Spalthoff, R., and Nenadić, I. (2017). Associations between urban upbringing and cortical thickness and gyrification. J. Psychiat. Res. 95, 114–120. doi: 10.1016/j.jpsychires.2017.08.012
Bethlehem, R. A., Romero-Garcia, R., Mak, E., Bullmore, E. T., and Baron-Cohen, S. (2017). Structural covariance networks in children with autism or ADHD. Cereb. Cortex 27, 4267–4276. doi: 10.1093/cercor/bhx135
Bolton, T. A. W., Tuleasca, C., Wotruba, D., Rey, G., Dhanis, H., Gauthier, B., et al. (2020). TbCAPs: A toolbox for co-activation pattern analysis. Neuroimage 211:116621. doi: 10.1016/j.neuroimage.2020.116621
Brickman, A. M., Habeck, C., Ramos, M. A., Scarmeas, N., and Stern, Y. (2008). A forward application of age associated gray and white matter networks. Hum. Brain Mapp. 29, 1139–1146. doi: 10.1002/hbm.20452
Bullmore, E., and Sporns, O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 10, 186–198. doi: 10.1038/nrn2575
Caligiuri, M. E., Arabia, G., Barbagallo, G., Lupo, A., Morelli, M., Nisticò, R., et al. (2017). Structural connectivity differences in essential tremor with and without resting tremor. J. Neurol. 264, 1865–1874. doi: 10.1007/s00415-017-8553-5
Chandran, V., and Pal, P. K. (2012). Essential tremor: beyond the motor features. Parkinson. Relat. Disord. 18, 407–413. doi: 10.1016/j.parkreldis.2011.12.003
Chen, Z. J., He, Y., Rosa-Neto, P., Germann, J., and Evans, A. C. (2008). Revealing modular architecture of human brain structural networks by using cortical thickness from MRI. Cereb. Cortex 18, 2374–2381. doi: 10.1093/cercor/bhn003
Chiarello, C., Vazquez, D., Felton, A., and McDowell, A. (2016). Structural asymmetry of the human cerebral cortex: regional and between-subject variability of surface area, cortical thickness and local gyrification. Neuropsychologia 93, 365–379. doi: 10.1016/j.neuropsychologia.2016.01.012
Chunling, W., and Zheng, X. (2016). Review on clinical update of essential tremor. Neurol. Sci. 37, 495–502. doi: 10.1007/s10072-015-2380-1
Corbetta, M., Shulman, G. L., Miezin, F. M., and Petersen, S. E. (1995). Superior parietal cortex activation during spatial attention shifts and visual feature conjunction. Science 270, 802–805. doi: 10.1126/science.270.5237.802
Dale, A. M., Fischl, B., and Sereno, M. I. (1999). Cortical surface-based analysis: I. segmentation and surface reconstruction. Neuroimage 9, 179–194. doi: 10.1006/nimg.1998.0395
Del Maschio, N., Fedeli, D., Sulpizio, S., and Abutalebi, J. (2019). The relationship between bilingual experience and gyrification in adulthood: A cross-sectional surface-based morphometry study. Brain Lang. 198:104680. doi: 10.1016/j.bandl.2019.104680
Deng, H., Wu, S., and Jankovic, J. (2019). Essential tremor: genetic update. Expert Rev. Mol. Med. 21:e8. doi: 10.1017/erm.2019.7
Derrible, S., and Kennedy, C. (2011). Applications of graph theory and network science to transit network design. Trans. Rev. 31, 495–519. doi: 10.1186/s12868-016-0283-6
Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021
Deuschl, G., and Elble, R. (2009). Essential tremor—neurodegenerative or nondegenerative disease towards a working definition of ET. Movement Disord. 24, 2033–2041. doi: 10.1002/mds.22755
Doucet, G. E., Lin, D., Du, Y., Fu, Z., Glahn, D. C., Calhoun, V. D., et al. (2020). Personalized estimates of morphometric similarity in bipolar disorder and schizophrenia. NPJ Schizophr. 6:39. doi: 10.1038/s41537-020-00128-x
Doucet, G. E., Moser, D. A., Rodrigue, A., Bassett, D. S., Glahn, D. C., and Frangou, S. (2019). Person-based brain morphometric similarity is heritable and correlates with biological features. Cereb. Cortex 29, 852–862. doi: 10.1093/cercor/bhy287
DuPre, E., and Spreng, R. N. (2017). Structural covariance networks across the life span, from 6 to 94 years of age. Netw. Neurosci. 1, 302–323. doi: 10.1162/NETN_a_00016
Elaimy, A. L., Demakas, J. J., Arthurs, B. J., Cooke, B. S., Fairbanks, R. K., Lamoreaux, W. T., et al. (2010). Gamma knife radiosurgery for essential tremor: a case report and review of the literature. World J. Surg. Oncol. 8:20. doi: 10.1186/1477-7819-8-20
Espay, A. J., Lang, A. E., Erro, R., Merola, A., Fasano, A., Berardelli, A., et al. (2017). Essential pitfalls in “essential” tremor. Movement Disord. 32, 325–331. doi: 10.1002/mds.26919
Evans, A. C. (2013). Networks of anatomical covariance. Neuroimage 80, 489–504. doi: 10.1016/j.neuroimage.2013.05.054
Fahn, S., Tolosa, E., and Marín, C. (1993). Clinical rating scale for tremor. Parkinson’s Dis. Movement Disord. 2, 271–280.
Farahani, F. V., Karwowski, W., and Lighthall, N. R. (2019). Application of graph theory for identifying connectivity patterns in human brain networks: a systematic review. Front. Neurosci. 13:585. doi: 10.3389/fnins.2019.00585
Fischl, B., and Dale, A. M. (2000). Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proc. Natl. Acad. Sci. 97, 11050–11055. doi: 10.1073/pnas.200033797
Fischl, B., Sereno, M. I., and Dale, A. M. (1999). Cortical surface-based analysis: II: inflation, flattening, and a surface-based coordinate system. Neuroimage 9, 195–207. doi: 10.1006/nimg.1998.0396
Freud, E., Plaut, D. C., and Behrmann, M. (2016). ‘What’ is happening in the dorsal visual pathway. Trends Cogn. Sci. 20, 773–784. doi: 10.1016/j.tics.2016.08.003
Gironell, A. (2014). The GABA hypothesis in essential tremor: lights and shadows. Tremor Other Hyperkin. Movements 4:254. doi: 10.7916/D8SF2T9C
Grindrod, P., and Kibble, M. (2004). Review of uses of network and graph theory concepts within proteomics. Expert Rev. Proteom. 1, 229–238. doi: 10.1586/14789450.1.2.229
Hafkemeijer, A., Möller, C., Dopper, E. G., Jiskoot, L. C., van den Berg-Huysmans, A. A., van Swieten, J. C., et al. (2016). Differences in structural covariance brain networks between behavioral variant frontotemporal dementia and Alzheimer’s disease. Hum. Brain Mapp. 37, 978–988. doi: 10.1002/hbm.23081
Hallquist, M. N., and Hillary, F. G. (2018). Graph theory approaches to functional network organization in brain disorders: A critique for a brave new small-world. Netw. Neurosci. 3, 1–26. doi: 10.1162/netn_a_00054
He, Y., Chen, Z., and Evans, A. (2008). Structural insights into aberrant topological patterns of large-scale cortical networks in Alzheimer’s disease. J. Neurosci. 28, 4756–4766. doi: 10.1523/JNEUROSCI.0141-08.2008
Huang, W., Bolton, T. A. W., Medaglia, J. D., Bassett, D. S., Ribeiro, A., and Van De Ville, D. (2018). A graph signal processing perspective on functional brain imaging. Proc. IEEE 106, 868–885. doi: 10.1109/jproc.2018.2798928
Ibrahim, M. F., Beevis, J. C., and Empson, R. M. (2020). Essential tremor–a cerebellar driven disorder? Neuroscience 462, 262–273. doi: 10.1016/j.neuroscience.2020.11.002
Jang, C., Park, H. J., Chang, W. S., Pae, C., and Chang, J. W. (2016). Immediate and longitudinal alterations of functional networks after thalamotomy in essential tremor. Front. Neurol. 7:184. doi: 10.3389/fneur.2016.00184
Jankovic, J. (2002). Essential tremor: a heterogenous disorder. Movement Disord. 17, 638–644. doi: 10.1002/mds.10221
Jhunjhunwala, K., and Pal, P. K. (2014). The non-motor features of essential tremor: a primary disease feature or just a secondary phenomenon? Tremor Other Hyperkin. Movements 4:255. doi: 10.7916/D8D798MZ
Jiménez-Jiménez, F. J., Alonso-Navarro, H., García-Martín, E., Alvarez, I., Pastor, P., and Agúndez, J. A. G. (2021a). Genomic markers for essential tremor. Pharmaceuticals 14:516. doi: 10.3390/ph14060516
Jiménez-Jiménez, F. J., Alonso-Navarro, H., García-Martín, E., and Agúndez, J. A. G. (2020). Sleep disorders in essential tremor: systematic review and meta-analysis. Sleep 43:zsaa039. doi: 10.1093/sleep/zsaa039
Jiménez-Jiménez, F. J., Alonso-Navarro, H., García-Martín, E., and Agúndez, J. A. G. (2021b). Sleep disorders in patients with essential tremor. Curr. Neurol. Neurosci. Rep. 21:23. doi: 10.1007/s11910-021-01109-y
Johnson, J. A., Strafella, A. P., and Zatorre, R. J. (2007). The role of the dorsolateral prefrontal cortex in bimodal divided attention: two transcranial magnetic stimulation studies. J. Cogn. Neurosci. 19, 907–920. doi: 10.1162/jocn.2007.19.6.907
Kelly, P. A., Viding, E., Wallace, G. L., Schaer, M., De Brito, S. A., Robustelli, B., et al. (2013). Cortical thickness, surface area, and gyrification abnormalities in children exposed to maltreatment: neural markers of vulnerability? Biol. Psychiatr. 74, 845–852. doi: 10.1016/j.biopsych.2013.06.020
Khundrakpam, B. S., Lewis, J. D., Reid, A., Karama, S., Zhao, L., Chouinard-Decorte, F., et al. (2017). Imaging structural covariance in the development of intelligence. Neuroimage 144, 227–240. doi: 10.1016/j.neuroimage.2016.08.041
Kondo, H., Osaka, N., and Osaka, M. (2004). Cooperation of the anterior cingulate cortex and dorsolateral prefrontal cortex for attention shifting. Neuroimage 23, 670–679. doi: 10.1016/j.neuroimage.2004.06.014
Kuhlenbäumer, G., Hopfner, F., and Deuschl, G. (2014). Genetics of essential tremor: meta-analysis and review. Neurology 82, 1000–1007. doi: 10.1212/WNL.0000000000000211
Li, J. Y., Suo, X. L., Li, N. N., Lei, D., Peng, J., Yang, J., et al. (2021). Disrupted brain network topology in drug-naïve essential tremor patients with and without depression. Clin. Neuroradiol. 31, 981–992. doi: 10.1007/s00062-021-01002-8
Li, R., Yuan, X., Radfar, M., Marendy, P., Ni, W., O’Brien, T. J., et al. (2021). Graph signal processing, graph neural network and graph learning on biological data: a systematic review. IEEE Rev. Biomed. Eng. [Epub ahead of print] doi: 10.1109/RBME.2021.3122522
Louis, E. D. (2016). Non-motor symptoms in essential tremor: a review of the current data and state of the field. Parkinson. Relat. Disord. 22, S115–S118. doi: 10.1016/j.parkreldis.2015.08.034
Louis, E. D. (2021). The essential tremors: evolving concepts of a family of diseases. Front. Neurol. 12:650601. doi: 10.3389/fneur.2021.650601
Louis, E. D., and Ferreira, J. J. (2010). How common is the most common adult movement disorder? Update on the worldwide prevalence of essential tremor. Movement Disord. 25, 534–541. doi: 10.1002/mds.22838
Luders, E., Kurth, F., Mayer, E. A., Toga, A. W., Narr, K. L., and Gaser, C. (2012). The unique brain anatomy of meditation practitioners: alterations in cortical gyrification. Front. Hum. Neurosci. 6:34. doi: 10.3389/fnhum.2012.00034
Mahayana, I. T., Tcheang, L., Chen, C. Y., Juan, C., and Muggleton, N. G. (2014). The precuneus and visuospatial attention in near and far space: a transcranial magnetic stimulation study. Brain Stimul. 7, 673–679. doi: 10.1016/j.brs.2014.06.012
Margulies, D. S., Vincent, J. L., Kelly, C., Lohmann, G., Uddin, L. Q., Biswal, B. B., et al. (2009). Precuneus shares intrinsic functional architecture in humans and monkeys. Proc. Natl. Acad. Sci. 106, 20069–20074. doi: 10.1073/pnas.0905314106
Masouleh, S. K., Plachti, A., Hoffstaedter, F., Eickhoff, S., and Genon, S. (2020). Characterizing the gradients of structural covariance in the human hippocampus. Neuroimage 218:116972. doi: 10.1016/j.neuroimage.2020.116972
McDowell, T., Holmes, N. P., Sunderland, A., and Schürmann, M. (2018). TMS over the supramarginal gyrus delays selection of appropriate grasp orientation during reaching and grasping tools for use. Cortex 103, 117–129. doi: 10.1016/j.cortex.2018.03.002
Mechelli, A., Friston, K. J., Frackowiak, R. S., and Price, C. J. (2005). Structural covariance in the human cortex. J. Neurosci. 25, 8303–8310. doi: 10.1523/jneurosci.0357-05.2005
Menon, V., and Uddin, L. Q. (2010). Saliency, switching, attention and control: a network model of insula function. Brain Struct. Func. 214, 655–667. doi: 10.1007/s00429-010-0262-0
Minati, L., Varotto, G., D’Incerti, L., Panzica, F., and Chan, D. (2013). From brain topography to brain topology: relevance of graph theory to functional neuroscience. Neuroreport 24, 536–543. doi: 10.1097/WNR.0b013e3283621234
Morgan, S. E., Seidlitz, J., Whitaker, K. J., Romero-Garcia, R., Clifton, N. E., Scarpazza, C., et al. (2019). Cortical patterning of abnormal morphometric similarity in psychosis is associated with brain expression of schizophrenia-related genes. Proc. Natl. Acad. Sci. 116, 9604–9609. doi: 10.1073/pnas.1820754116
Papini, C., Palaniyappan, L., Kroll, J., Froudist-Walsh, S., Murray, R. M., and Nosarti, C. (2020). Altered cortical gyrification in adults who were born very preterm and its associations with cognition and mental health. Biol. Psychiatr. 5, 640–650. doi: 10.1016/j.bpsc.2020.01.006
Pedersen, M., Zalesky, A., Omidvarnia, A., and Jackson, G. D. (2018). Multilayer network switching rate predicts brain performance. Proc. Natl. Acad. Sci. 115, 13376–13381. doi: 10.1073/pnas.1814785115
Persson, J., Spreng, R. N., Turner, G., Herlitz, A., Morell, A., Stening, E., et al. (2014). Sex differences in volume and structural covariance of the anterior and posterior hippocampus. Neuroimage 99, 215–225. doi: 10.1016/j.neuroimage.2014.05.038
Potok, W., Maskiewicz, A., Króliczak, G., and Marangon, M. (2019). The temporal involvement of the left supramarginal gyrus in planning functional grasps: a neuronavigated TMS study. Cortex 111, 16–34. doi: 10.1016/j.cortex.2018.10.010
Prasad, S., Shah, A., Saini, J., Ingalhalikar, M., and Pal, P. K. (2020). Aberrant global and local efficiency of the executive subnetwork in essential tremor. J. Neural Trans. 127, 385–388. doi: 10.1007/s00702-020-02141-x
Preti, M. G., Bolton, T. A. W., and Van De Ville, D. (2017). The dynamic functional connectome: state-of-the-art and perspectives. Neuroimage 160, 41–54. doi: 10.1016/j.neuroimage.2016.12.061
Raethjen, J., and Deuschl, G. (2012). The oscillating central network of essential tremor. Clin. Neurophysiol. 123, 61–64. doi: 10.1016/j.clinph.2011.09.024
Raznahan, A., Shaw, P., Lalonde, F., Stockman, M., Wallace, G. L., Greenstein, D., et al. (2011). How does your cortex grow? J. Neurosci. 31, 7174–7177. doi: 10.1523/JNEUROSCI.0054-11.2011
Romero-Garcia, R., Whitaker, K. J., Váša, F., Seidlitz, J., Shinn, M., Fonagy, P., et al. (2018). Structural covariance networks are coupled to expression of genes enriched in supragranular layers of the human cortex. Neuroimage 171, 256–267. doi: 10.1016/j.neuroimage.2017.12.060
Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003
Sanabria-Diaz, G., Melie-García, L., Iturria-Medina, Y., Alemán-Gómez, Y., Hernández-González, G., Valdés-Urrutia, L., et al. (2010). Surface area and cortical thickness descriptors reveal different attributes of the structural human brain networks. Neuroimage 50, 1497–1510. doi: 10.1016/j.neuroimage.2010.01.028
Schaefer, A., Kong, R., Gordon, E. M., Laumann, T. O., Zuo, X., Holmes, A. J., et al. (2018). Local-global parcellation of the human cerebral cortex from intrinsic functional connectivity MRI. Cereb. Cortex 28, 3095–3114. doi: 10.1093/cercor/bhx179
Schmitt, J. E., Lenroot, R. K., Wallace, G. L., Ordaz, S., Taylor, K. N., Kabani, N., et al. (2008). Identification of genetically mediated cortical networks: a multivariate study of pediatric twins and siblings. Cereb. Cortex 18, 1737–1747. doi: 10.1093/cercor/bhm211
Schnack, H. G., Van Haren, N. E. M., Brouwer, R. M., Evans, A., Durston, S., Boomsma, D. I., et al. (2015). Changes in thickness and surface area of the human cortex and their relationship with intelligence. Cereb. Cortex 25, 1608–1617. doi: 10.1093/cercor/bht357
Schultz, C. C., Koch, K., Wagner, G., Roebel, M., Nenadic, I., Gaser, C., et al. (2010). Increased parahippocampal and lingual gyrification in first-episode schizophrenia. Schizoph. Res. 123, 137–144. doi: 10.1016/j.schres.2010.08.033
Schultz, C. C., Wagner, G., Koch, K., Gaser, C., Roebel, M., Schachtzabel, C., et al. (2013). The visual cortex in schizophrenia: alterations of gyrification rather than cortical thickness—a combined cortical shape analysis. Brain Struct. Funct. 218, 51–58. doi: 10.1007/s00429-011-0374-1
Ségonne, F., Dale, A. M., Busa, E., Glessner, M., Salat, D., Hahn, H. K., et al. (2004). A hybrid approach to the skull stripping problem in MRI. Neuroimage 22, 1060–1075. doi: 10.1016/j.neuroimage.2004.03.032
Seidlitz, J., Váša, F., Shinn, M., Romero-Garcia, R., Whitaker, K. J., Vértes, P. E., et al. (2018). Morphometric similarity networks detect microscale cortical organization and predict inter-individual cognitive variation. Neuron 97, 231–247. doi: 10.1016/j.neuron.2017.11.039
Sharifi, S., Nederveen, A. J., Booij, J., and van Rootselaar, A. F. (2014). Neuroimaging essentials in essential tremor: a systematic review. NeuroImage 5, 217–231. doi: 10.1016/j.nicl.2014.05.003
Siokas, V., Aloizou, A. M., Tsouris, Z., Liampas, I., Aslanidou, P., Dastamani, M., et al. (2020). Genetic risk factors for essential tremor: a review. Tremor Other Hyperkin. Movements 10:4. doi: 10.5334/tohm.67
Sotiras, A., Resnick, S. M., and Davatzikos, C. (2015). Finding imaging patterns of structural covariance via non-negative matrix factorization. Neuroimage 108, 1–16. doi: 10.1016/j.neuroimage.2014.11.045
Sotiras, A., Toledo, J. B., Gur, R. E., Gur, R. C., Satterthwaite, T. D., and Davatzikos, C. (2017). Patterns of coordinated cortical remodeling during adolescence and their associations with functional specialization and evolutionary expansion. Proc. Natl. Acad. Sci. 114, 3527–3532. doi: 10.1073/pnas.1620928114
Sporns, O. (2018). Graph theory methods: applications in brain networks. Dial. Clin. Neurosci. 20, 111–121. doi: 10.31887/DCNS.2018.20.2/osporns
Tian, Y., Margulies, D. S., Breakspear, M., and Zalesky, A. (2020). Topographic organization of the human subcortex unveiled with functional connectivity gradients. Nat. Neurosci. 23, 1421–1432. doi: 10.1038/s41593-020-00711-6
Tio, M., and Tan, E. K. (2016). Genetics of essential tremor. Parkinson. Relat. Disord. 22, S176–S178.
Tröster, A. I., Pahwa, R., Fields, J. A., Tanner, C. M., and Lyons, K. E. (2005). Quality of life in essential tremor questionnaire (QUEST): development and initial validation. Parkinson. Relat. Disord. 11, 367–373. doi: 10.1016/j.parkreldis.2005.05.009
Tuleasca, C., Witjas, T., Najdenovska, E., Verger, A., Girard, N., Champoudry, J., et al. (2017). Assessing the clinical outcome of Vim radiosurgery with voxel-based morphometry: visual areas are linked with tremor arrest! Acta Neurochirurg. 159, 2139–2144. doi: 10.1007/s00701-017-3317-7
Varjačić, A., Mantini, D., Levenstein, J., Slavkova, E. D., Demeyere, N., and Gillebert, C. R. (2018). The role of left insula in executive set-switching: lesion evidence from an acute stroke cohort. Cortex 107, 92–101. doi: 10.1016/j.cortex.2017.11.009
Xue, H. L., Liu, G., and Yang, X. H. (2016). A review of graph theory application research in gears. Proc. Instit. Mechan. Eng. Part C 230, 1697–1714. doi: 10.1177/0954406215583321
Yang, J. J., Kwon, H., and Lee, J. M. (2016). Complementary characteristics of correlation patterns in morphometric correlation networks of cortical thickness, surface area, and gray matter volume. Sci. Rep. 6:26682. doi: 10.1038/srep26682
Yang, J., Lei, D., Peng, J., Suo, X., Pinaya, W. H. L., Li, W., et al. (2021). Disrupted brain gray matter networks in drug-naïve participants with essential tremor. Neuroradiology 63, 1501–1510. doi: 10.1007/s00234-021-02653-7
Keywords: mean curvature, surface area, cortical thickness, graph theory, structural covariance analysis, stereotactic radiosurgical thalamotomy, radiosurgery, essential tremor
Citation: Bolton TAW, Van De Ville D, Régis J, Witjas T, Girard N, Levivier M and Tuleasca C (2022) Graph Theoretical Analysis of Structural Covariance Reveals the Relevance of Visuospatial and Attentional Areas in Essential Tremor Recovery After Stereotactic Radiosurgical Thalamotomy. Front. Aging Neurosci. 14:873605. doi: 10.3389/fnagi.2022.873605
Received: 11 February 2022; Accepted: 19 April 2022;
Published: 23 May 2022.
Edited by:
Anke Meyer-Baese, Florida State University, United StatesReviewed by:
Jurgen Germann, University Health Network (UHN), CanadaFélix Javier Jiménez-Jiménez, Hospital Universitario del Sureste, Spain
Copyright © 2022 Bolton, Van De Ville, Régis, Witjas, Girard, Levivier and Tuleasca. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Thomas A. W. Bolton, thomas.bolton@chuv.ch