Abstract
Inefficient integration between bottom-up visual input and higher order visual processing regions is implicated in visual hallucinations in Parkinson’s disease (PD). Here, we investigated white matter contributions to this perceptual imbalance hypothesis. Twenty-nine PD patients were assessed for hallucinatory behavior. Hallucination severity was correlated to connectivity strength of the network using the network-based statistic approach. The results showed that hallucination severity was associated with reduced connectivity within a subnetwork that included the majority of the diverse club. This network showed overall greater between-module scores compared with nodes not associated with hallucination severity. Reduced between-module connectivity in the lateral occipital cortex, insula, and pars orbitalis and decreased within-module connectivity in the prefrontal, somatosensory, and primary visual cortices were associated with hallucination severity. Conversely, hallucination severity was associated with increased between- and within-module connectivity in the orbitofrontal and temporal cortex, as well as regions comprising the dorsal attentional and default mode network. These results suggest that hallucination severity is associated with marked alterations in structural network topology with changes in participation along the perceptual hierarchy. This may result in the inefficient transfer of information that gives rise to hallucinations in PD.
Author Summary
Inefficient integration of information between external stimuli and internal perceptual predictions may lead to misperceptions or visual hallucinations in Parkinson’s disease (PD). In this study, we show that hallucinatory behavior in PD patients is associated with marked alterations in structural network topology. Severity of hallucinatory behavior was associated with decreased connectivity in a large subnetwork that included the majority of the diverse club, nodes with a high number of between-module connections. Furthermore, changes in between-module connectivity were found across brain regions involved in visual processing, top-down prediction centers, and endogenous attention, including the occipital, orbitofrontal, and posterior cingulate cortex. Together, these findings suggest that impaired integration across different sides across different perceptual processing regions may result in inefficient transfer of information.
INTRODUCTION
Visual hallucinations (VHs) in Parkinson’s disease (PD) exist on a spectrum ranging from simple misperceptions to complex well-formed images (Barnes & David, 2001). With disease progression and loss of insight, VHs constitute a major source of distress for the patient (Goetz, 2009; Schrag, 2004) and comprise a high degree of burden for caregivers (Aarsland et al., 2007). Risk factors of VHs include older age and disease duration, sleep and mood disturbances, as well as cognitive decline (Barnes & David, 2001; Fénelon, Mahieux, Huon, & Ziégler, 2000; Lenka, Hegde, Arumugham, & Pal, 2017). Furthermore, previous work has shown that patients with VHs show disruptions in attentional processing (Hall et al., 2016), reduced performance on visuoperceptive tasks (Barnes, Boubert, Harris, Lee, & David, 2003; Gallagher et al., 2011; Ramírez-Ruiz, Junqué, Martí, Valldeoriola, & Tolosa, 2006), and decreased visual contrast sensitivity, color discrimination (Diederich et al., 1998), and acuity (Matsui et al., 2006). Current models of VHs have therefore focused on the interaction of perceptual and attentional dysfunction (for a review, see Muller, Shine, Halliday, & Lewis, 2014). Specifically, it has been proposed that failure to effectively integrate information from different processing sites across the perceptual hierarchy is likely to contribute to VHs and misperceptions in PD (Collerton, Perry, & McKeith, 2005; Diederich, Goetz, & Stebbins, 2005; Muller et al., 2014; Shine, Halliday, Carlos, Naismith, & Lewis, 2012).
Attention, prior experience, and expectations strongly influence perception. Perceptual predictions, generated from a myriad of modalities across the brain, guide perceptual processes to facilitate the interpretation of noisy and ambiguous input (Bar, 2009; Engel, Fries, & Singer, 2001; Summerfield et al., 2006). The orbitofrontal cortex (OFC) processes coarse information projected from the visual cortex and provides an “initial guess” of an object’s identity (Summerfield & Koechlin, 2008). Previous work in PD patients with VHs has shown that the accumulation of sensory evidence is slow and inefficient, which may result in an overreliance on these top-down predictions (O’Callaghan et al., 2017). Importantly, top-down visual processing regions can modulate neural activity in early visual regions, with expected stimuli leading to reduced activity (Meyer & Olson, 2011). Additionally, activity within the default mode network (DMN), a network involved in mediating endogenous perception, has shown to be increased during a misperception in this patient population (Shine, Halliday, et al., 2014). Therefore, VHs may arise when perceptual input is not properly integrated and internally generated images interfere with the perceptual process (Fletcher & Frith, 2008; Intaite, Noreika, Soliunas, & Falter, 2013; O’Callaghan et al., 2017; Powers, Kelley, & Corlett, 2016).
While functional neuroimaging studies have made significant contributions to our understanding (Hepp, Foncke, Olde Dubbelink, et al., 2017; Ramírez-Ruiz et al., 2008; Shine, Halliday, et al., 2014; Shine, Muller, et al., 2015; Yao et al., 2014), less is known about the involvement of white matter changes in the manifestation of VHs in PD. Experiments using diffusion tensor imaging (DTI) have reported altered white matter integrity in the optic nerve and optic radiation (Lee et al., 2016) as well as ascending tracts from the cholinergic nucleus basalis of Meynert to parietal and occipital cortical regions (Hepp, Foncke, Berendse, et al., 2017). However, given the involvement of large-scale brain networks in perception, unique insights into white matter changes associated with VHs can be gained by investigating whole-brain network topology. Topological features of the human connectome allow us to describe the arrangement of connections within and between segregated submodules (Bullmore & Sporns, 2009). Specifically, nodes that integrate these specialist communities are crucial for incorporating information streams of different modalities, which is essential for processes such as perception (Bertolero, Yeo, & D’Esposito, 2015; Muller, O’Callaghan, Walton, Shine, & Lewis, 2017). Therefore, investigating network topology can provide novel insights in changes across different perceptual hierarchies.
The current study aimed to examine whether VHs are associated with changes in structural network topology. To identify hallucinatory behavior in patients with PD, we assessed performance on a computerized task capable of inducing misperceptions (Shine et al., 2012), in combination with a validated questionnaire that evaluates hallucinatory behavior in PD patients (Shine, Mills, et al., 2015). We aimed to circumvent the sole reliance on self-reported symptom occurrence while controlling for the possibility of misclassifying hallucinators as nonhallucinators, when patients did not experience VHs during their clinic visit or when “passage” hallucinations go unreported. Additionally, by creating this composite score, we are able to assess visual hallucinatory severity, rather than classifying patients into artificial dichotomous patient groups.
We hypothesized that the severity of hallucinatory behavior would be associated with ineffective information processing as shown by reduced between-module scores in visual networks, reflecting reduced visual input to integration centers. Furthermore, increased between-module scores across top-down perceptual prediction areas and the DMN could indicate an overreliance on regions involved in the generation of internal percepts (Shine, O’Callaghan, Halliday, & Lewis, 2014).
METHODS
Twenty-nine patients with idiopathic PD were included in this study. Demographic information including age, disease duration, and levodopa equivalence daily dose (LEDD) were obtained for all participants. All patients were assessed on the Hoehn & Yahr clinical stage (Hoehn & Yahr, 1998) and the motor aspect of the Movement Disorder Society Unified Parkinson’s Disease Rating Scale (MDS-UPDRS) part III (Goetz et al., 2008). Global cognition was assessed using the Mini-Mental State Examination (MMSE; Folstein, Robins, & Helzer, 1983), and set-shifting performance was assessed using the Trail Making Test Part B minus Part A (TMTB-A; Tombaugh, 2004). The study was approved by the ethics committee of the University of Sydney and was in accordance with the principles of the Helsinki Declaration. Written informed consent was obtained from all participants before participation.
Bistable Percept Paradigm
All patients performed the bistable percept paradigm (BPP; Shine et al., 2012), a behavioral task capable of inducing misperceptions in susceptible patients. In this task, patients were presented with either single or bistable percepts (i.e., “hidden” images as shown in Figure 1) for a maximum of 30 s in a randomized order. The participant had to decide whether the stimulus was a single or hidden image by a button press and describe to the examiner what they had seen. The recorded responses included the following: (a) correct single or correct hidden, (b) “missed,” recorded when the subject perceived a single image when a bistable percept was presented, and (c) “misperceptions,” recorded when a subject incorrectly identified a single image as a bistable image, that is, incorrectly reported an image that was not presented on the screen.
Psychosis and Hallucinations Questionnaire
All patients completed the Psychosis and Hallucinations Questionnaire (PsycH-Q; Shine, Mills, et al., 2015). The PsycH-Q (Part A) consists of three subscales including (a) visual misperceptions, which includes questions about the presence of VHs, passage hallucinations, and three frequently reported contents of VHs including people, animals, and objects; (b) sensory misperceptions, including audition, touch, olfaction, and gustation; and (c) disordered thought and psychotic behavior. Participants rated the frequency of their symptoms on a 5-point Likert scale, ranging from 0 (“never experienced”) to 4 (“experienced daily”). The total score was calculated by summing the responses (Shine, Mills, et al., 2015; see Supporting Information, Hall et al., 2019). Part B of the PsycH-Q assesses symptoms related to VHs (i.e., attention and sleep) and was not included in this study.
Composite Score
The percentage of misperceptions on the BPP (“indirect” measure of VH) and the total score on the PsycH-QA (“direct” measure of hallucinatory behavior) were standardized and then summed to create a composite score that reflected the severity of hallucinatory behavior (hereafter referred to as the hallucination severity score, HSS). The HSS was correlated with the demographic variables using parametric or nonparametric correlations depending on the distribution of the variables and was used as a correlate in the imaging analysis.
MRI Acquisition
All participants underwent magnetic resonance imaging (MRI) using a 3-Tesla General Electric Discovery MR750 scanner (GE Medical Systems) with an 8-channel phased array head coil. Diffusion-weighted images (DWI) were obtained by using echo-planar imaging sequences with 61 different motion-probing gradient directions (TR/TE: 7,025/80 ms, 55 transverse slices, slice thickness: 2.5 mm, matrix: 256 × 256, FOV: 240 × 240 mm). The effective diffusion weighting was b = 1,000 s/mm2, and four volumes with no diffusion weighting (b = 0 s/mm2) were obtained at the beginning of each diffusion sequence. 3-D T1-weighted, anatomical images were obtained (TR/TE/TI: 7.2/2.7/450 ms, voxel size 1 × 1 × 1 mm, 196 transverse slices, 256 × 256 matrix, FOV: 256 × 256 mm, flip angle 12°). The 3-D T1 images were used for individual registration between T1-weighted anatomical and the DWI images and cortical parcellation using FreeSurfer (version 5.3; http://surfer.nmr.mgh.harvard.edu).
Diffusion Tensor Imaging Preprocessing and Deterministic Fiber Tracking
DTI preprocessing was performed using the FMRIB Software Library (FSL, http://fsl.fmrib.ox.ac.uk). The preprocessing steps were as follows: (a) DTI images were corrected for susceptibility, head motion, and eddy current–induced geometrical distortions using FSL’s tool eddy; (b) a binary brain mask was created using bet; (c) images were realigned using a rigid body registration to the b = 0 image; then (d) a tensor was fitted in each voxel (Chang, Jones, & Pierpaoli, 2005), followed by (e) the computation of the fractional anisotropy (FA) level based on the eigenvalues for each voxel, in order to determine the preferred diffusion direction within a voxel. FA thus serves as a surrogate measure of white matter integrity, with lower levels of FA reflecting reduced white matter integrity (Sun et al., 2003; van den Heuvel & Sporns, 2011; Verstraete et al., 2010). The preferred diffusion direction information was then used to reconstruct the white matter tracts of the brain using a deterministic tracking approach based on the fiber assignment by continuous tracking (FACT) algorithm (Mori, Crain, Chacko, & van Zijl, 1999). Deterministic tractography yields less false positive tracts compared with probabilistic methods (Bastiani, Shah, Goebel, & Roebroeck, 2012). False positives are detrimental in network modularity as they occur more prevalently between than within modules (Zalesky et al., 2016). A streamline was started from eight seeds within each voxel of the brain (gray and white matter) following the main diffusion direction of the voxel and stopped when (a) the FA value < 0.1; (b) the traced fiber made a turn > 45°; or (c) the tract left the brain mask. The images were acquired when reverse phase-encoding direction approaches were not the standard procedure within acquisition protocols, which could have influenced the registration of diffusion and anatomical images. Therefore, anatomically constrained tractography was not applied (Smith, Tournier, Calamante, & Connelly, 2012). The atlas presented by Cammoun et al. (2012) was used, including 219 cortical regions and 14 subcortical regions. The weighted brain network was calculated for each participant, and consistency thresholding at 50% was applied (i.e., including the tracts found in 50% of the patients; de Reus & van den Heuvel, 2013). The mean density of the thresholded group matrix was 8.7%. To verify the results were not skewed by the choice of threshold, we also applied the thresholding method that retained most consistent edges across subjects but controlling for their distance (i.e., the consistency of edges within “bins” based on their length to avoid preferential retention of short edges; Misic et al., 2015). The mean density of the group matrix using this threshold was 13.2%.
Network-Based Statistic
A network-based statistic (NBS) analysis was applied to investigate whether the HSS was associated with altered connectivity strength in an interconnected subnetwork of the brain (Zalesky, Fornito, & Bullmore, 2010). NBS is a nonparametric method for connectome-wide analysis, which aims to detect specific pairs of brain regions showing a significant effect of interest, while controlling for family-wise error (FWE) rate. Importantly, no inferences of individual connections are made; instead the null hypothesis can only be rejected at the subnetwork level. As such, NBS is similar to the cluster-based multiple-comparison approaches used in standard functional MRI analysis. To identify changes in subnetworks associated with the HSS, the t statistic was set at 1.7, determined using the critical value of the t distribution for our sample size (Field, 2009). Additionally, we verified the robustness of the results by controlling for disease severity using the MDS-UPDRS III (motor part) as covariate. Connections were deemed significant at FWE-corrected p value < 0.05 (one-sided) using 5,000 permutations.
To investigate whether the subnetwork involved particular functional networks, we investigated whether nodes in the subnetwork that correlated with the HSS overlapped with previously defined resting-state networks. To this end, seven canonical resting-state networks from the Yeo et al. (2011) atlas were overlaid with the structural parcellation and the percentage of nodes from each network included within the structural subnetwork that inversely related to HSS was calculated for each resting-state network. To analyze whether this overlap occurred significantly above chance, we randomly permuted the resting-state network identity of each region (5,000 iterations) and used the overlap between the randomized vector and the original node assignment to populate a null distribution. To test whether each individual resting-state network overlapped with the significant subnetwork, their overlap was compared with the null distributions. A resting-state network was identified as targeted if the true overlap was more than the 97.5th percentile of null distribution (i.e., the top 2.5%). A network was considered not to be associated with the HSS if the overlap was less than the 2.5th percentile of the null distribution.
Graph Theoretical Analysis
The graph organizational measures were computed using the Brain Connectivity Toolbox (http://www.brain-connectivity-toolbox.net; Rubinov & Sporns, 2010). The thresholded, weighted brain networks were then partitioned into modules, which are nonoverlapping groups of highly connected nodes that are only sparsely connected with other modules, using the Louvain algorithm (Rubinov & Sporns, 2010). To account for the stochastic nature of the Louvain algorithm, a consensus partition was identified by calculating the module assignment for each node 500 times. To define an appropriate value for the resolution parameter (γ), the Louvain algorithm was iterated 100 times across a range of values (0.5–2.0 in steps of 0.1) of the group mean connectivity matrix and then the similarity of the resultant partitions was estimated using mutual information. The γ parameter of 1.9 provided the most robust estimates of topology across the iterations and was used to determine the optimal resolution of the network modularity.
To test whether nodes within the subnetwork identified using the NBS analysis differed from nodes not included in the subnetwork, the average Wi and Bi of the subnetwork were contrasted against the average Wi and Bi of the nodes not included in the subnetwork using nonparametric permutation testing.
To test whether the HSS correlated with the Wi and Bi nodes across the whole-brain connectome, a Spearman’s rho correlation was performed followed by a nonlinear permutation test using 5,000 iterations to control for multiple comparisons (Nichols & Holmes, 2002), using an alpha of 0.05. This approach was repeated using the different threshold (Misic et al., 2015), and the outcome was correlated to the Wi and Bi using the original threshold. Both the Wi and the Bi scores calculated using the aforementioned threshold highly correlated with the Wi and Bi scores calculated with the consensus threshold (r = 0.92 and r = 0.94, respectively), indicating that the results were not biased by the chosen thresholding method.
Diverse Club Analysis
We identified the “diverse club” of the network, which comprised the top 20% of Bi nodes (Bertolero, Yeo, & D’Esposito, 2017). These nodes play an important role in network integration, and changes to these nodes could affect between-module communication (Bertolero et al., 2017). We normalized the diverse club coefficient in reference to a null model: A random vector with a preserved modular structure was created by randomizing the mean participation coefficient of each node for 5,000 iterations. The diverse club was identified as those regions with a participation coefficient greater than the 95th percentile of the permuted distribution. We investigated whether the number of diverse club nodes was significantly higher within the subnetwork associated with the HSS, compared with nodes that were not included in the subnetwork identified using the NBS analysis.
RESULTS
Demographics
Table 1 presents the descriptive variables of the 29 patients. The mean percentage of misperceptions on the BPP was 18.48 (range: 0–49), and the mean score on the PsycH-QA was 9.48 (range: 0–34, max score = 52), highlighting a diverse range of hallucinatory behavior in the patient cohort. The PsycH-QA and the BPP score showed a positive, significant correlation (r = 0.52, p = 0.004). Finally, to verify the severity score to the “gold standard,” we correlated the HSS in a large cohort of patients with PD and Lewy body dementia (n = 75) to the MDS-UPDRS item 2 and found a correlation of r = 0.53 (p < 0.001). However, given higher construct validity (Shine, Mills, et al., 2015), we opted to include the PsycH-QA and BP scores in the composite score (HSS) for the remainder of our analysis.
| Variable | Mean (range) | Correlation with HSS r (p value) |
|---|---|---|
| Age (y) | 66.8 ± 8 (51–84) | 0.18 (0.345) |
| Duration (y) | 5.8 ± 4 (1.2–16) | −0.05 (0.785) |
| LEDD | 617.8 ± 392 (125–1548) | −0.06 (0.767) |
| H & Y | 1.9 ± 0.5 (1–3) | 0.01 (0.949) |
| MDS-UPDRS III | 30.0 ± 14 (7–55) | 0.29 (0.125) |
| MMSE | 28.5 ± 2 (25–30) | −0.06 (0.762) |
| TMTB-A | 78.9 ± 65 (−1–123) | 0.33 (0.080) |
| BPP % misperceptions | 18.5 ± 17 (0–49) | – |
| PsycH-QA | 9.5 ± 8 (0–34) | – |
| Variable | Mean (range) | Correlation with HSS r (p value) |
|---|---|---|
| Age (y) | 66.8 ± 8 (51–84) | 0.18 (0.345) |
| Duration (y) | 5.8 ± 4 (1.2–16) | −0.05 (0.785) |
| LEDD | 617.8 ± 392 (125–1548) | −0.06 (0.767) |
| H & Y | 1.9 ± 0.5 (1–3) | 0.01 (0.949) |
| MDS-UPDRS III | 30.0 ± 14 (7–55) | 0.29 (0.125) |
| MMSE | 28.5 ± 2 (25–30) | −0.06 (0.762) |
| TMTB-A | 78.9 ± 65 (−1–123) | 0.33 (0.080) |
| BPP % misperceptions | 18.5 ± 17 (0–49) | – |
| PsycH-QA | 9.5 ± 8 (0–34) | – |
LEDD = levodopa equivalence daily dose; MDS-UPDRS III = motor part of the Movement Disorder Society Unified Parkinson’s Disease Rating Scale; H & Y = Hoehn and Yahr; MMSE = Mini-Mental State Examination; TMTB-A = Trail Making Test Part B – Part A; BPP = bistable percept paradigm; PsycH-QA = Psychosis and Hallucinations Questionnaire, Part A.
The HSS showed a positive correlation trending towards significance with the TMTB-A (r = 0.33, p = 0.08). No significant correlations were observed between the HSS and other demographic and clinical variables.
The HSS Correlated With Decreased Connectivity in a Large Subnetwork
As illustrated in Figure 2, the NBS analysis revealed a subnetwork comprising 183 edges (8% of the edges in the thresholded connectivity matrix) and 127 nodes with reduced FA-based connectivity strength correlated to the HSS (p < 0.05). Using disease severity as a covariate, the NBS analysis revealed a similar subnetwork comprising 177 edges and 135 nodes that showed a correlation with the HSS that was trending towards significance (p = 0.059). The effects presented with a fairly liberal threshold, suggesting the changes related to the HSS are subtle yet topological extended (Zalesky et al., 2010). The size of networks identified using a range of t statistics are presented in the Supporting Information (Hall et al., 2019). No significant subnetwork was identified in the opposite direction (positive correlation between the HSS and connectivity strength). Furthermore, the group average Bi score within the subnetwork was 0.506, which was significantly higher (p < 0.05) than nodes outside this network, which show a group average of 0.310. The group average Wi score of the nodes within the network was higher than the group average Wi score of nodes outside the network (0.149 and −0.178, respectively), yet this difference did not remain significant when controlling for disease severity (p = 0.104).
NBS analysis reveals a subnetwork, comprising 183 edges and 127 nodes with reduced connectivity strength correlated to increased HSS (p < 0.05). Figure visualized with BrainNet Viewer (Xia, Wang, & He, 2013).
NBS analysis reveals a subnetwork, comprising 183 edges and 127 nodes with reduced connectivity strength correlated to increased HSS (p < 0.05). Figure visualized with BrainNet Viewer (Xia, Wang, & He, 2013).
The Subnetwork Includes All Subcortical Nodes but Did Not Target a Specific Cortical Resting-State Network
The subnetwork that showed decreased connectivity strength correlated with the HSS included all 14 subcortical nodes (p < 0.05). As illustrated in Figure 3, the subnetwork further included nodes across the cortex. However, none of the other the resting-state networks were significantly correlated with the HSS (p > 0.05), yet the somatomotor network was relatively spared (p < 0.05).
Overlap between the identified structural subnetwork and functional resting-state networks. (A) The percentage of nodes included in the subnetwork for each resting-state network. DMN = default mode network; VAN = ventral attentional network; DAN = dorsal attentional network. (B) The functional resting-state networks of the Yeo et al. (2011) atlas.
Overlap between the identified structural subnetwork and functional resting-state networks. (A) The percentage of nodes included in the subnetwork for each resting-state network. DMN = default mode network; VAN = ventral attentional network; DAN = dorsal attentional network. (B) The functional resting-state networks of the Yeo et al. (2011) atlas.
Nodes in the Subnetwork Show High Participation Scores
Eighteen nodes were included in the diverse club (see Supporting Information, Hall et al., 2019). Seventeen of the eighteen nodes (94%) of the diverse club were included in the aforementioned subnetwork, which was deemed significantly above chance (p < 0.001). As illustrated in Figure 4, nodes with high participation coefficients were more often part of the subnetwork.
Nodes ranked according to the Bi scores. Blue: nodes included in the subnetwork; gray: nodes not included in the subnetwork correlated to the HSS.
Nodes ranked according to the Bi scores. Blue: nodes included in the subnetwork; gray: nodes not included in the subnetwork correlated to the HSS.
The HSS Correlated With Wi and Bi Scores
When investigating the whole structural connectome, the HSS positively correlated to regional Bi (i.e., higher participation scores were associated with higher severity values) for nodes in the left medial OFC, a node in the right anterior and left posterior cingulate, precuneus, and the caudal middle frontal gyrus. Furthermore, nodes in the right occipital, pars orbitalis, and insula showed negative correlations between the HSS and participation coefficient (i.e., lower participation scores were associated with higher scores on the HSS; see Table 2 and Figure 5). However, when controlling for disease severity, the insula and medial OFC were only trending towards significance, while the lateral occipital cortex and anterior posterior cingulate did not remain significant.
| Node | X | Y | Z | Rho | Rho* | Subnetwork |
|---|---|---|---|---|---|---|
| Positively correlated | ||||||
| Frontal | ||||||
| ctx-lh-medialorbitofrontal_2 | −5 | 33 | −20 | 0.37 | 0.35# | ✓ |
| ctx-rh-caudalmiddlefrontal_2 | 40 | 15 | 39 | 0.43 | 0.53 | – |
| Parietal | ||||||
| ctx-lh-precuneus_2 | −10 | −44 | 46 | 0.40 | 0.37 | ✓ |
| Cingulate | ||||||
| ctx-lh-posteriorcingulate_2 | −8 | −43 | 21 | 0.57 | 0.53 | ✓ |
| ctx-rh-superiorfrontal_3 | 11 | 40 | 40 | 0.36 | 0.31ˆ | ✓ |
| Negatively correlated | ||||||
| Frontal | ||||||
| ctx-rh-parsorbitalis_1 | 43 | 43 | −10 | −0.46 | −0.46 | ✓ |
| Occipital | ||||||
| ctx-rh-lateraloccipital_5 | 47 | −73 | 2 | −0.39 | −0.21ˆ | ✓ |
| Insular | ||||||
| ctx-rh-insula_1 | 36 | −18 | 13 | −0.40 | −0.32# | ✓ |
| Node | X | Y | Z | Rho | Rho* | Subnetwork |
|---|---|---|---|---|---|---|
| Positively correlated | ||||||
| Frontal | ||||||
| ctx-lh-medialorbitofrontal_2 | −5 | 33 | −20 | 0.37 | 0.35# | ✓ |
| ctx-rh-caudalmiddlefrontal_2 | 40 | 15 | 39 | 0.43 | 0.53 | – |
| Parietal | ||||||
| ctx-lh-precuneus_2 | −10 | −44 | 46 | 0.40 | 0.37 | ✓ |
| Cingulate | ||||||
| ctx-lh-posteriorcingulate_2 | −8 | −43 | 21 | 0.57 | 0.53 | ✓ |
| ctx-rh-superiorfrontal_3 | 11 | 40 | 40 | 0.36 | 0.31ˆ | ✓ |
| Negatively correlated | ||||||
| Frontal | ||||||
| ctx-rh-parsorbitalis_1 | 43 | 43 | −10 | −0.46 | −0.46 | ✓ |
| Occipital | ||||||
| ctx-rh-lateraloccipital_5 | 47 | −73 | 2 | −0.39 | −0.21ˆ | ✓ |
| Insular | ||||||
| ctx-rh-insula_1 | 36 | −18 | 13 | −0.40 | −0.32# | ✓ |
Rho after controlling for disease severity;
p < 0.1;
p > 0.1.
ctx-lh = left hemisphere; ctx-rh = right hemisphere.
Significant changes in between- and within-module scores associated with the HSS. (A) Nodes with a significant correlation between participation coefficients and the HSS. (B) Nodes with a significant correlation between the module degree z-score and HSS. Green indicates a positive correlation, red indicates a negative correlation. Larger nodes are part of the subnetwork. Figure visualized with BrainNet Viewer (Xia et al., 2013).
Significant changes in between- and within-module scores associated with the HSS. (A) Nodes with a significant correlation between participation coefficients and the HSS. (B) Nodes with a significant correlation between the module degree z-score and HSS. Green indicates a positive correlation, red indicates a negative correlation. Larger nodes are part of the subnetwork. Figure visualized with BrainNet Viewer (Xia et al., 2013).
Increased HSS scores were further associated with increased module degree z-scores in the right thalamus, bilateral lingual, left medial OFC, pars opercularis, supramarginal gyrus, and superior temporal cortices. The right lateral occipital cortex also showed a positive correlation, but this did not remain significant after controlling for disease severity. Decreased HSS scores were associated with increased module degree z-scores in the bilateral precuneus, left parts triangularis, rostral middle frontal, and superior parietal cortex. The right pericalcarine and precentral gyrus also a showed negative correlation with the HSS, yet this did not remain significant after controlling for disease severity (see Table 3 and Figure 5).
| Node | X | Y | Z | Rho | Rho* | Subnetwork |
|---|---|---|---|---|---|---|
| Positively correlated | ||||||
| Subcortical | ||||||
| Right-Thalamus-Proper | 13 | −17 | 7 | 0.49 | 0.45 | ✓ |
| Frontal | ||||||
| ctx-lh-medialorbitofrontal_1 | −7 | 52 | −11 | 0.40 | 0.43 | – |
| ctx-lh-parsopercularis_2 | −47 | 14 | 15 | 0.41 | 0.44 | – |
| Parietal | ||||||
| ctx-lh-supramarginal_3 | −55 | −34 | 35 | 0.45 | 0.53 | – |
| Temporal | ||||||
| ctx-rh-superiortemporal_4 | 46 | −6 | −11 | 0.57 | 0.47 | – |
| Occipital | ||||||
| ctx-rh-lateraloccipital_3 | 27 | −94 | −6 | 0.39 | 0.34# | – |
| ctx-lh-lingual_2 | −10 | −70 | 0 | 0.42 | 0.50 | ✓ |
| ctx-rh-lingual_1 | 14 | −83 | −7 | 0.41 | 0.37 | – |
| Negatively correlated | ||||||
| Frontal | ||||||
| ctx-lh-parstriangularis_1 | −43 | 33 | 2 | −0.50 | −0.43 | ✓ |
| ctx-lh-rostralmiddlefrontal_5 | −31 | 56 | 2 | −0.57 | −0.50 | ✓ |
| ctx-rh-precentral_6 | 22 | −17 | 67 | −0.39 | −0.35# | ✓ |
| Parietal | ||||||
| ctx-lh-precuneus_2 | −10 | −44 | 46 | −0.46 | −0.43 | ✓ |
| ctx-lh-superiorparietal_6 | −22 | −76 | 29 | −0.41 | −0.47 | ✓ |
| ctx-rh-precuneus_2 | 12 | −54 | 26 | −0.39 | −0.43 | ✓ |
| Occipital | ||||||
| ctx-rh-pericalcarine_1 | 15 | −75 | 12 | −0.38 | −0.24ˆ | ✓ |
| Node | X | Y | Z | Rho | Rho* | Subnetwork |
|---|---|---|---|---|---|---|
| Positively correlated | ||||||
| Subcortical | ||||||
| Right-Thalamus-Proper | 13 | −17 | 7 | 0.49 | 0.45 | ✓ |
| Frontal | ||||||
| ctx-lh-medialorbitofrontal_1 | −7 | 52 | −11 | 0.40 | 0.43 | – |
| ctx-lh-parsopercularis_2 | −47 | 14 | 15 | 0.41 | 0.44 | – |
| Parietal | ||||||
| ctx-lh-supramarginal_3 | −55 | −34 | 35 | 0.45 | 0.53 | – |
| Temporal | ||||||
| ctx-rh-superiortemporal_4 | 46 | −6 | −11 | 0.57 | 0.47 | – |
| Occipital | ||||||
| ctx-rh-lateraloccipital_3 | 27 | −94 | −6 | 0.39 | 0.34# | – |
| ctx-lh-lingual_2 | −10 | −70 | 0 | 0.42 | 0.50 | ✓ |
| ctx-rh-lingual_1 | 14 | −83 | −7 | 0.41 | 0.37 | – |
| Negatively correlated | ||||||
| Frontal | ||||||
| ctx-lh-parstriangularis_1 | −43 | 33 | 2 | −0.50 | −0.43 | ✓ |
| ctx-lh-rostralmiddlefrontal_5 | −31 | 56 | 2 | −0.57 | −0.50 | ✓ |
| ctx-rh-precentral_6 | 22 | −17 | 67 | −0.39 | −0.35# | ✓ |
| Parietal | ||||||
| ctx-lh-precuneus_2 | −10 | −44 | 46 | −0.46 | −0.43 | ✓ |
| ctx-lh-superiorparietal_6 | −22 | −76 | 29 | −0.41 | −0.47 | ✓ |
| ctx-rh-precuneus_2 | 12 | −54 | 26 | −0.39 | −0.43 | ✓ |
| Occipital | ||||||
| ctx-rh-pericalcarine_1 | 15 | −75 | 12 | −0.38 | −0.24ˆ | ✓ |
Rho after controlling for disease severity;
p < 0.1;
p > 0.1.
ctx-lh = left hemisphere; ctx-rh = right hemisphere.
DISCUSSION
The aim of this study was to determine whether changes in structural network topology were associated with hallucinatory behavior in PD. We showed that severity of hallucinatory behavior was negatively correlated with connectivity across a bilateral subnetwork. Regions within this subnetwork showed higher participation compared with regions outside this network. The loss of connectivity strength may force the system to adapt and reroute information across less efficient pathways, impeding the standard sensory integration process. Importantly, 94% of the nodes in the diverse club were included in this subnetwork. This community of high participation nodes is thought to control the integration of relatively segregated regions (Bertolero et al., 2017). Indeed, the diverse connectivity pattern of these nodes makes them crucial for the functional coordination of brain regions during tasks, and activity in these nodes predicts changes in the coupling of other regions (Bertolero et al., 2017). Severity of hallucinatory behavior may thus be the result of impaired integration and segregation of brain networks or “modules,” affecting effective information transfer. Finally, we showed regional changes in participation associated with hallucination severity (the HSS score), with a positive correlation between participation scores in the medial OFC, cingulate, precuneus, and middle frontal gyrus and the HSS and negative correlation with participation scores in the lateral occipital cortex, pars orbitalis, and insula. These findings suggest a reweighting of the regions along the perceptual hierarchy, which may give rise to hallucinations.
Lower participation of the lateral occipital cortex may reflect reduced early visual processing, resulting in ineffective accumulation of visual information from the environment. Previous work using a Bayesian drift diffusion model has demonstrated that accumulation speed and quality of perceptual information are reduced in PD patients with VH (O’Callaghan et al., 2017). Furthermore, reduced quality or integration of visual stimuli may increase perceptual uncertainty, a suggestion that aligns with increased participation in the dorsal anterior cingulate cortex (Keri, Decety, Roland, & Gulyas, 2004). Perceptual uncertainty may place excessive emphasis on top-down prediction centers, which subsequently could lead to a reduced activity in early visual regions (Meyer & Olson, 2011). However, the reduced participation score in the occipital and anterior cingulate cortex may not be specific to hallucinatory behavior, as the negative correlations did not remain significant after controlling for disease severity.
This emphasis on top-down visual processing centers is supported by the increased participation coefficient and module degree z-score in the medial OFC. The OFC has an integrative function across brain networks, as evidenced by its high participation coefficient. Additionally, this region is thought to facilitate recognition during visual perception by integrating incoming sensory information with previous experiences and expectations (Panichello, Cheung, & Bar, 2012). During typical visual perception, the OFC is activated early in response to visual stimuli, receiving low spatial frequency signals from the visual cortex (Bar et al., 2006). Notably, only stimuli resembling known objects are shown to activate the OFC, which in turn generates a semantic association and provides a predictive signal to the visual system (Chaumon, Kveraga, Barrett, & Bar, 2014). Conversely, visual stimuli that carry no meaningful association do not activate the OFC in healthy individuals. Hence, it could be speculated that because of decreased quality of visual input, inappropriate recruitment of the OFC occurs, which may result in ascribing false associative information to visual stimuli.
The manifestation of VHs has previously been recognized as a dysfunction between the attentional networks (Shine, Muller, et al., 2015). Specifically, patients with VHs are shown to be less able to recruit the dorsal attentional network (DAN), which enables the selection of appropriate sensory stimuli (Corbetta, Patel, & Shulman, 2008). With reduced control of this network, ambiguous stimuli might instead be interpreted by the ventral attentional network (VAN) and the DMN, which are less well equipped for this task. Our results showed increased participation in the posterior cingulate cortex (PCC), a key hub of the DMN. PCC activity has been implicated in regulating the focus of attention, specifically the shift from the external world into internal mentation (Weissman, Roberts, Visscher, & Woldorff, 2006). Furthermore, the PCC is involved in mind wandering and supports internally directed cognition (Leech & Sharp, 2013). A failure to suppress PCC activity may lead to the intrusion of internal thoughts into task performance (Sonuga-Barke & Castellanos, 2007). Moreover, a positive correlation was found between the HSS and the module degree z-score of the left pars opercularis, a node in the VAN, a network that is activated when expectations in perception are violated (Corbetta & Shulman, 2002; Shine, Halliday, et al., 2014). Conversely, a negative correlation between the HSS scores and module degree score and participation coefficients was found in other nodes of the VAN, namely the left pars triangularis, the right pars orbitalis, and insula. The left pars triangularis supports resolving competition between simultaneously active representations (Badre & Wagner, 2007), while the insula plays an important role in dynamically shifting attention between the attentional control networks (Menon & Uddin, 2010). The anterior insula has previously been shown to be reduced in volume in PD patients with VH (Ibarretxe-Bilbao et al., 2010; Shine, Halliday, et al., 2014). Together, these results suggest that ineffective communication between attentional networks in the brain may predispose an individual to hallucinate. Surprisingly, the participation coefficient of a node within the DAN (“ctx-rh-caudalmiddlefrontal_2”) showed a positive correlation with the HSS. This node was not part of the subnetwork, and it may be possible that this is a compensatory response to the loss of connectivity strength in the other DAN regions. Notably, the connectivity matrix shows between module connections of this region with nodes in the somatomotor and the frontoparietal network, but not with the DMN or VAN.
Finally, all nodes that showed negative correlations with the HSS were included in the subnetwork. Decreased within-module scores were found across the prefrontal and the somatosensory association cortex, as well as in the primary visual cortex, while there was a positive correlation between the HSS and the bilateral secondary visual cortex, perhaps as a result of the decreased visual input from V1. Additionally, the supramarginal gyrus, a node that has been shown to be functionally active during spatial perception but also during visual imagery (Ganis, Thompson, & Kosslyn, 2004), showed an increased module degree z-score with increasing severity of VHs. Furthermore, a positive correlation with the HSS and the module degree score in the superior temporal cortex, a region involved in auditory processing, was found. It could be speculated that increased visual uncertainty may stimulate other sensory processing areas. Conversely, previous work in schizophrenia has shown increased activation in the superior temporal cortex during a mismatch between expected and perceived auditory feedback (Fu et al., 2008). The HSS also included hallucinations of other sensory modalities, and it could be speculated that a greater emphasis is being placed on predictions of auditory stimuli, as reflected as increased within-module connectivity in the superior temporal cortex.
This study has several limitations worth noting. First, the DWI data were acquired without EPI distortion correction. This may have affected the accuracy of registration between DWI and T1 images in the frontal and temporal cortices. Because of relatively low diffusion weighting used in the current MRI protocol, it was chosen to employ DTI rather than more sophisticated methods such as constrained spherical deconvolution, an algorithm that more adequately deals with multiple fiber directions within one voxel than DTI. Furthermore, after controlling for disease severity, the identified subnetwork was only trending towards significance. The correlation between the HSS and the participation coefficients of the lateral occipital and anterior cingulate did not remain significant, while the OFC and insula were only trending towards significance, and similar patterns were found for the module degree z-score: correlation with the HSS and nodes in the occipital and precentral regions did not remain significant. However, this may be expected given the close relationship between overall symptom severity and phenotypic symptom severity. Importantly, the correlation coefficients showed notable consistency across the analyses, suggesting the HSS was able to distinguish interindividual variability. While the current study included PD patients with no self-reported hallucinations (n = 6), no data from a healthy control group were acquired, which could have aided in disassociating the significant network measures and disease burden. The LEDD score was not associated with the HSS in the current study, and the relation between hallucinations and a LEDD-related side effect has been dismissed in previous work (Fénelon et al., 2000; Graham, Grünewald, & Sagar, 1997; Sanchez-Ramos, Ortollm, & Paulson, 1996; Shergill, Walker, & Le Katona, 1998). Furthermore, hallucinations may occur in the absence of dopaminergic treatment in PD or dementia with Lewy bodies (Ala, Yang, Sung, & Frey, 1997). Therefore, the current study chose not to include LEDD score as a covariate. Furthermore, the investigation was conducted in a relatively small group of PD patients, and future studies should replicate our findings in a larger sample size. Finally, this study chose to include the PsycH-Q to assess hallucinations. This is an in-depth questionnaire with high construct validity; however, patients have to report their hallucinations themselves. Ideally, a trained neurologist or qualified researcher, who can probe the patients when in doubt, should assist with this process.
CONCLUSIONS
We cautiously conclude that hallucinatory behavior in PD patients is associated with marked alterations in structural network topology. Severity of hallucinatory behavior was associated with decreased connectivity in a large subnetwork that included the majority of the diverse club. These changes may result in an inefficient rerouting of information across less efficient pathways, which may lead to impaired visual integration processes. Furthermore, nodes within the orbitofrontal cortex and temporal lobes showed increased participation scores, while the visual association cortex, insula, and middle frontal gyrus showed decreased scores associated with the HSS score. These findings suggest that impaired integration across different regions along the perceptual hierarchy may result in inefficient transfer of information. A failure to effectively switch between attentional networks and the intrusion of internal percepts could give rise to perceptual glitches, such as misperceptions and hallucinations.
AUTHOR CONTRIBUTIONS
Julie M. Hall: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Project administration; Visualization; Writing – original draft. Claire O’Callaghan: Data curation; Investigation; Methodology; Writing – review & editing. Alana Justine Muller: Data curation; Writing – review & editing. Kaylena A. Ehgoetz Martens: Writing – review & editing. Joseph R. Phillips: Writing – review & editing. Ahmed A. Moustafa: Writing – review & editing. Simon J. G. Lewis: Funding acquisition; Investigation; Supervision; Writing – review & editing. James M. Shine: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Supervision; Writing – review & editing.
FUNDING INFORMATION
Claire O’Callaghan, National Health and Medical Research Council Neil Hamilton Fairley Fellowship, Australia, Award ID: 1091310. Simon J. G. Lewis, National Health and Medical Research Council of Australia, Award ID: 1037746. Simon J. G. Lewis, National Health and Medical Research Council of Australia, Award ID: 1095127. Simon J. G. Lewis, National Health and Medical Research Council of Australia and the Australian Research Council, Dementia Fellowship, AWARD ID: 1110414. James M. Shine, National Health and Medical Research Council Project Grant, Australia, Award ID: 1156536. James M. Shine, The University of Sydney Robinson Fellowship.
ACKNOWLEDGMENTS
We thank the patients and their families who contribute to our research at the Parkinson’s Disease Research Clinic. We thank Dr. Váša for sharing his thresholding code (https://github.com/frantisekvasa/matlab_general). The DWI data were processed during the 10kin1day initiative at the Dutch Connectome Lab. This research was supported by Sydney Informatics Hub, funded by the University of Sydney.
TECHNICAL TERMS
- MDS Unified Parkinson’s Disease Rating Scale:
A commonly used measure of Parkinson’s disease clinical symptom severity.
- Bistable percept paradigm:
A computer-based task using bistable percepts to measure visual misperceptions.
- Psychosis and Hallucinations Questionnaire:
A 20-item self-reported questionnaire to measure hallucinatory phenomenology.
- Tractography:
Three-dimensional representation of the brain’s white matter tracts derived from diffusion MRI data.
- Network-based statistic:
A statistical method to identify connections and networks comprising the connectome associated with behavior.
- Module degree z-score:
A measure of diversity of intramodular connections of individual nodes.
- Participation coefficient:
A measure of diversity of intermodular connections of individual nodes.
- Diverse club:
A subset of nodes in a network with high participation coefficients that exhibit properties consistent with an integrative network function.
- Structural connectome:
Brain connectivity graphs obtained from diffusion MRI and tractography.