Graph Theory-Guided Transcranial Magnetic Stimulation in Neurodegenerative Disorders

The emergence of brain imaging techniques has shed light on the underlying mechanisms of brain diseases. Abnormalities found in the diseased brain have often been targeted by transcranial magnetic stimulation (TMS), which can modulate long-term neuronal excitability in humans in a noninvasive manner. Thus, its therapeutic application has been extensively explored for neurological and psychiatric disorders. As a result, TMS has been approved for clinical treatment for a few diseases, including drug-refractory depression. However, one of the biggest challenges with TMS is the difficulty in finding the optimal stimulation site. Until now, this process has been heavily dependent on previous activation studies and anatomical knowledge of the region itself, but it largely ignored the whole brain network that interacts with the focal brain region that has been targeted. Here, we propose a novel approach to estimate the prospective network effect following focal interference induced by TMS as a way to optimize the target identification process for TMS research.

The development of brain imaging techniques, such as positron emission tomography (PET), single photon emission computed tomography (SPECT) and functional magnetic resonance imaging (fMRI), has greatly advanced our understanding of the human brain at the systems level. [¹⁸F]-fluorodeoxyglucose (FDG) is the most widely used PET radiotracer which can quantify glucose metabolism in neural tissue. It has been well documented that glucose metabolism correlates with neuronal activity (1); abnormal neuronal activity associated with neurological or psychiatric disorders may be identified using FDG PET. The most conventional approach for PET data analysis is statistical parametric mapping (SPM), which localizes the brain region that deviates from normal FDG uptake in control subjects (2,3). The degree of this regional difference may be utilized as a biomarker for diseases (4,5).

Multivariate analysis utilizes the covariances of whole-brain metabolic activity (6,7). For example, the scaled subprofile mapping (SSM)—a form of principal component analysis (PCA)—has been used to analyze FDG PET scan data in patients with Parkinson’s disease (PD). The resulting spatial covariance pattern, PD-related metabolic pattern (PDRP) (8,9), significantly separates PD from the control group, and the network pattern has been repeatedly validated in independent PD populations [cf. (10,11)]. One of the most advantageous features of multivariate analysis is that the degree of pattern expression can be quantified in each individual (12). The rest-state measurements of each subject’s PDRP expression have been found to correlate consistently with clinical ratings of PD motor dysfunction and have been employed as a treatment biomarker in clinical trials of novel antiparkinsonian interventions (13). In addition, an automated metabolic imaging-based algorithm has demonstrated highly accurate results in differentiating idiopathic PD from other atypical parkinsonian syndromes, especially in the early stages of the disease (14). Notably, the prediction of accurate diagnosis was possible even before patients’ motor symptoms emerged (15).

Transcranial magnetic stimulation (TMS) is a technique that has been used widely for noninvasive brain stimulation to examine motor, perceptual and cognitive processes (16–18). When applied continuously, repetitive TMS (rTMS) has been shown to produce profound and long-lasting effects on neuronal excitability. Therefore, when combined with neuroimaging techniques, one can assess the causality of the brain activity and behavior using the “perturb and measure” protocol [cf. (19)].

Traditionally, rTMS can be subdivided into fast (>5 Hz) and slow (<1 Hz), regarded as facilitatory and inhibitory, respectively (18). Interestingly, initial rTMS-PET studies reported increased regional cerebral blood flow (rCBF) at the stimulated site by both fast (20) and slow (21) rTMS. FDG uptake was also increased at the site of stimulation with intermediate frequency (2 Hz) (22). The unexpected increase in rCBF after slow (inhibitory) rTMS was explained by the active inhibition process (23,24). A few studies, however, found reduced rCBF and FDG uptake at the stimulated site (25,26) after slow rTMS.

The effect of rTMS in remote areas from the target site is more complicated. Recent literature suggests that the effect of rTMS can vary depending on the location, frequency, intensity, number of stimulation and the functional status of the target site or population (27). For example, the primary motor area (M1) is the most extensively investigated region in TMS-PET studies due to the large existing literature from neurophysiology (28). Majority of the TMS-PET studies reported increased rCBF at the stimulation site (20–22,29–33), although a few studies reported otherwise (25,34). The remote effects of M1 rTMS has often involved the “motor network” [for review, see (28)]. For example, Lee et al. (35) proposed that 1 Hz-rTMS applied to M1 induced a reorganization of the motor network by increasing the connectivity of the contralateral hemisphere to compensate for the inhibited M1 stimulation site. A similar approach was used to understand how patients with aphasia activate right-sided language function (36). Laird et al. (37) summarized the motor network that is influenced by M1 TMS using structured equation modeling and metaanalysis, and revealed that rTMS on M1 remotely changes activities in the supplementary motor area (SMA), anterior cingulate cortex (ACC), thalamus, premotor cortex (PMC), posterior parietal cortex and cerebellum. However, it should be noted that the direction of remote effect has not been consistent. For example, contralateral M1 activity was both increased (33,38) and decreased (21,30) after M1 rTMS.

The effect of rTMS on the prefrontal cortex (PFC) has not been consistent (33,39–42). Similar to M1 rTMS affecting motor network activities (37), PFC rTMS has modulated activities in the cognitive network comprised of ACC (40,43), insula (33,40), PMC (41,42), basal ganglia (33,40,43), parahippocampal region (39,41,42) and cerebellum (33,40). With respect to the clinical benefits, rCBF changes in the subgenual cingulate and orbitofrontal cortex after rTMS in the left dorsolateral prefrontal cortex have been reported previously to be associated with the efficacy of rTMS treatment in depression (44).

While several studies investigated the effects of rTMS on neurodegenerative disorders, target identification has been a challenging issue. Most rTMS studies defined the target region based on activation studies [for example, (45,46)] and/or anatomical knowledge. The current practice largely ignores the prospective remote effects, although it has been clearly suggested that the whole-brain network is likely affected even if the stimulation was applied focally (19,47). Here, we propose a novel research protocol for TMS target identification.

In the current study, we have applied a novel network analytical strategy based upon spatial covariance mapping (48,49) and graph theory (50,51) to identify the hub of a network that is associated with cognitive deficits in PD, that is, PD cognitive pattern (PDCP) (52). It has been estimated that this hub region exhibits the greatest ability to modulate the expression of the entire network during local stimulation. Because of the normalization inherent in the iterative eigenvector computations that we employ (53,54), the loadings on the spatial covariance patterns do not vary sufficiently across regions for the delineation of discrete modules within the network space. That is, a different computational algorithm is necessary to define the organizational structure (topology) of the network after quantifying the relevant subject scores. To assess the flow of information within the PDCP, we used a graph theoretic approach (50). Adjacency matrices were constructed from interindividual correlation maps of 15 PD patients who were used in the derivation of PDCP (52). Automated anatomical labeling (AAL) was utilized to construct regions of interest (ROIs) for the whole brain regions (55). The resulting data matrix is both nonnegative and primitive, and is therefore assured based on the Perron-Frobenius theorem that it will have a single positive eigenvector. The component of this vector for a given ROI is called the eigenvector centrality (EC) for that region. EC represents a measure of information flow ranking each node according to its relative importance (“centrality”) within the network space (56). Alternatively, in the context of a stochastic process, EC represents the likelihood of arriving at a given node after a long run of random walks irrespective of the starting point. The resulting data also can be used to deduce the organizational structure of an anatomical-functional network (for example, radial, linear, star-shaped, or core-periphery configurations).

The utility of graph theory in this context is demonstrated in a recent analysis of metabolic imaging data from PD patients (n = 99) and healthy volunteer subjects (n = 66) scanned in our laboratory (57). Extensive permutation testing of each group showed a strong relationship between the magnitude of the regional loadings (region weights) on the primary metabolic covariance pattern (the first principal component [PC]) and corresponding nodal EC values determined independently by graph theory. Whether metabolically overactive (region weight z > 1.0) or underactive (z < −1.0), the salient regions on the primary PC pattern can be viewed as major hubs of information flow along the graph theoretic network. Similarly, the regions with salient contributions to the PDCP topography functioned as information hubs along that network construct. This relationship was indicated by the presence of a significant correlation (r = 0.66, p < 0.001) between regional loadings on the PDCP (representing the relative regional contributions to network activity) and the corresponding EC values (representing the relative importance of a given region to information flow along the network).

The sensitivity of each node to stimulation was quantified using EC derivative analysis (51) to identify regions that facilitate (or suppress) the influence of other “important” network hubs (defined as PDCP nodes with region weight |z| > 1). For example, the right superior frontal gyrus (rSFG, BA 6/8) showed the highest EC derivative. Thus, altering the connectivity of this “sensitive” region by rTMS is likely to have the greatest influence on the centrality of the other salient PDCP nodes (that is, those with region weights).

By normalizing information flow within the PDCP network, therapeutic rSFG stimulation is expected to reverse the abnormal pattern expression that characterizes PD patients with cognitive dysfunction. The right middle frontal gyrus (rMFG, BA 8/9) has been viewed as the core of PD cognitive dysfunction (58–60) and thus chosen as the target of therapeutic rTMS in studies directed at this disease manifestation (45,61). While this region exhibited greater importance (higher EC) than the rSFG, it had lower sensitivity, that is, an EC derivative value of relatively lower magnitude, making it comparatively less desirable as a treatment target.

We propose an innovative network targeting approach to optimize cortical rTMS procedures to normalize aberrant brain networks, such as PDCP, which is associated with cognitive deficits in PD patients. To this end, disease-related network expression levels may be measured before and after rTMS is applied to one or more network nodes selected based on graph theoretic considerations.

D Eidelberg serves on the scientific advisory boards of and has received honoraria from the Michael J Fox Foundation for Parkinson’s Research and the Bachmann-Strauss Dystonia and Parkinson Foundation; is listed as a coinventor of patents named “Markers for Use in Screening Patients for Nervous System Dysfunction and a Method and Apparatus for Using Same” without financial gain; and has received research support from the NIH (NINDS, NIDCD, NIAID).

This work was supported in part by the National Institute of Neurological Disorders and Stroke Morris K Udall Center of Excellence for Parkinson’s Disease Research at The Feinstein Institute for Medical Research (P50 NS071675 to D Eidelberg). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Neurological Disorders and Stroke or the National Institutes of Health. The sponsor did not play a role in study design, collection, analysis and interpretation of data, writing of the report or in the decision to submit the paper for publication.

Affiliations

1. Center for Neurosciences, The Feinstein Institute for Medical Research, 350 Community Drive, Manhasset, New York, 11030, USA
   • Ji Hyun Ko
   • , Yoon Young Choi
   •  & David Eidelberg
2. Department of Human Anatomy and Cell Science, University of Manitoba, Winnipeg, Manitoba, Canada
   • Ji Hyun Ko
3. Neuroscience Research Program, Kleysen Institute for Advanced Medicine Health Science Centre, Winnipeg, Manitoba, Canada
   • Ji Hyun Ko

Corresponding author

Correspondence to David Eidelberg.

Reprints and Permissions