Locus coeruleus integrity and the effect of atomoxetine on response inhibition in Parkinson's disease.

Cognitive decline is a common feature of Parkinson's disease, and many of these cognitive deficits fail to respond to dopaminergic therapy. Therefore, targeting other neuromodulatory systems represents an important therapeutic strategy. Among these, the locus coeruleus-noradrenaline system has been extensively implicated in response inhibition deficits. Restoring noradrenaline levels using the noradrenergic reuptake inhibitor atomoxetine can improve response inhibition in some patients with Parkinson's disease, but there is considerable heterogeneity in treatment response. Accurately predicting the patients who would benefit from therapies targeting this neurotransmitter system remains a critical goal, in order to design the necessary clinical trials with stratified patient selection to establish the therapeutic potential of atomoxetine. Here, we test the hypothesis that integrity of the noradrenergic locus coeruleus explains the variation in improvement of response inhibition following atomoxetine. In a double-blind placebo-controlled randomized crossover design, 19 patients with Parkinson's disease completed an acute psychopharmacological challenge with 40 mg of oral atomoxetine or placebo. A stop-signal task was used to measure response inhibition, with stop-signal reaction times obtained through hierarchical Bayesian estimation of an ex-Gaussian race model. Twenty-six control subjects completed the same task without undergoing the drug manipulation. In a separate session, patients and controls underwent ultra-high field 7 T imaging of the locus coeruleus using a neuromelanin-sensitive magnetization transfer sequence. The principal result was that atomoxetine improved stop-signal reaction times in those patients with lower locus coeruleus integrity. This was in the context of a general impairment in response inhibition, as patients on placebo had longer stop-signal reaction times compared to controls. We also found that the caudal portion of the locus coeruleus showed the largest neuromelanin signal decrease in the patients compared to controls. Our results highlight a link between the integrity of the noradrenergic locus coeruleus and response inhibition in patients with Parkinson's disease. Furthermore, they demonstrate the importance of baseline noradrenergic state in determining the response to atomoxetine. We suggest that locus coeruleus neuromelanin imaging offers a marker of noradrenergic capacity that could be used to stratify patients in trials of noradrenergic therapy and to ultimately inform personalized treatment approaches.

Introduction

Cognitive decline in Parkinson’s disease remains an ongoing therapeutic challenge. The mainstay dopaminergic therapies often fail to improve cognitive deficits, and in some cases can exacerbate them. This has prompted a shift of focus towards other neuromodulatory systems that are affected by Parkinson’s disease and related to cognitive decline, including noradrenaline. The noradrenergic locus coeruleus is one of the earliest sites of alpha-synuclein pathology, and noradrenergic treatments have been shown to modulate cognitive functions that are impaired in Parkinson’s disease, including response inhibition.

Response inhibition deficits are a well-documented feature of Parkinson’s disease, ranging from subclinical impairments to extreme impulsive behaviours. Impulsivity is clear in the florid ‘impulse control disorders’ that are exacerbated by dopaminergic therapy. However, milder impulsivity is common in the absence of an impulse control disorder, including impairments in the ability to cancel an inappropriate action. Neurodegeneration of fronto-striatal circuits, including the subthalamic nucleus and its inputs, contributes to this impairment in Parkinson’s disease,, while pharmacological modulation of these circuits offers a tractable route to restorative treatment.,

The locus coeruleus-noradrenaline system modulates the stimulus detection and behavioural reorienting required for rapid action cancellation., Phasic activation in the locus coeruleus—the brain’s main source of noradrenaline—occurs in response to salient events, and its activity is tightly time-locked to task-relevant responses. Highly collateralized projections from the locus coeruleus enable release of noradrenaline in multiple brain regions, altering the gain, or responsivity, of target neurons. The action of noradrenaline at multiple targets can interrupt and reconfigure network architecture, promoting a change in goal-directed behaviour., This locus coeruleus-noradrenaline function directly supports rapid action cancellation. In healthy adults and in preclinical models, pharmacologically increasing noradrenaline levels with the reuptake inhibitor atomoxetine improves action cancellation, as measured on stop-signal tasks.

Atomoxetine selectively inhibits presynaptic noradrenaline transporters, resulting in a 3-fold increase in extracellular levels of noradrenaline in the prefrontal cortex. It is currently licensed for treating behavioural and cognitive symptoms associated with attention deficit hyperactivity disorder. However, experimental psychopharmacological studies indicate it may be of value in some patients with Parkinson’s disease, by increasing activity and connectivity in the fronto-striatal ‘stopping network’.,, The stopping network includes the inferior frontal gyrus and presupplementary motor area, and their excitatory connection with the subthalamic nucleus which, via the globus pallidus, increases inhibition over thalamocortical output., Within this network, noradrenaline increases cortical excitability, functional connectivity and network integration. The prefrontal regions also provide descending input to modulate the locus coeruleus,, by which prefrontal noradrenaline can influence locus coeruleus activity. Atomoxetine alters locus coeruleus firing patterns to increase the phasic-to-tonic ratio, making the locus coeruleus more responsive to task-relevant stimuli.

The potential for atomoxetine to modulate locus coeruleus-noradrenaline function and improve response inhibition holds therapeutic promise in Parkinson’s disease. Previous work using stop-signal tasks in Parkinson’s disease demonstrated that atomoxetine can improve response inhibition and enhance its attendant stopping network activation., However, there was considerable heterogeneity in treatment response. To accurately predict the patients who would benefit from noradrenergic therapy remains a critical goal for atomoxetine treatment to be considered therapeutically and to design the necessary clinical trials with stratified patient selection.,

Here, we test the hypothesis that structural integrity of the noradrenergic locus coeruleus explains the variation in improvements in response inhibition following atomoxetine. This is now possible through recent developments in ultra-high field 7 T imaging of the locus coeruleus. We tested this hypothesis by combining quantification of the locus coeruleus by 7 T MRI with an acute psychopharmacological challenge, and measuring response inhibition using the stop-signal reaction time (SSRT).

Materials and methods

Participants

Nineteen patients with idiopathic Parkinson’s disease were recruited via the University of Cambridge Parkinson’s disease research clinic and the Parkinson’s UK volunteer network. They met the United Kingdom Parkinson’s Disease Society Brain Bank criteria and were not demented based on Movement Disorders Society (MDS) criteria for Parkinson’s disease dementia nor on the Mini-Mental State Examination (MMSE, score > 26). They were aged between 50 and 80 years, with Hoehn and Yahr stages 1.5–3, and had no contraindications to 7 T MRI or atomoxetine. None had current impulse control disorders, based on clinical impression and the Questionnaire for Impulsive-Compulsive Disorders in Parkinson’s Disease (QUIP-Current Short) screening tool. Levodopa equivalent daily dose (LEDD) scores were calculated. Twenty-six age-, sex- and education-matched healthy control subjects were recruited from local volunteer panels. Control participants were screened for a history of neurological or psychiatric disorders, and no controls were using psychoactive medications. The study was approved by the local Ethics Committees and all participants provided written informed consent according to the Declaration of Helsinki. Demographic details and clinical characteristics are provided in Table 1 and the Supplementary material.

Table 1

Demographics and clinical assessments of participants in their normal ON medication state

	Parkinson's disease	Controls	BF	P
Age, years	67.11 (7.05)	65.35 (5.32)	0.43	0.368
Education, years	14.05 (2.27)	14.65 (3.10)	0.37	0.457
Male/female	15/4	15/11	0.98	0.240
MMSE	29.47 (0.70)	29.77 (0.51)	0.87	0.128
MoCA	28.11 (1.76)	28.58 (1.39)	0.45	0.340
ACE-R
Total score	94.89 (3.71)	97.58 (3.16)	4.16	0.015
Attention and orientation	17.84 (0.37)	17.96 (0.20)	0.64	0.216
Memory	23.68 (1.97)	25.04 (1.18)	7.00	0.013
Fluency	12.00 (2.08)	12.81 (1.60)	0.71	0.167
Language	25.84 (0.50)	25.88 (0.43)	0.31	0.768
Visuospatial	15.63 (0.50)	15.81 (0.63)	0.45	0.302
MDS-UPDRS
I: Non-motor experiences	9.00 (4.18)
II: Motor experiences	12.63 (4.26)
III: Motor examination	28.42 (11.60)
IV: Motor complications	0.47 (0.96)
Total score	50.58 (17.20)
Hoehn and Yahr stage	2.26 (0.45)
Disease duration, years	4.15 (1.72)
Levodopa equivalent daily dose, mg/day	644.55 (492.81)

Data are presented as mean (SD). Comparisons of patient and control groups were performed with independent samples t-tests or contingency tables as appropriate. ACE-R = revised Addenbrooke’s Cognitive Examination; BF = Bayes factor where >3 indicates substantial evidence of a group difference; MDS-UPDRS = Movement Disorders Society Unified Parkinson’s Disease Rating Scale; MoCA = Montreal Cognitive Assessment; P = P-value of classical two-sample t-test.

Study procedure

Participants with Parkinson’s disease were tested across three sessions. First, they underwent MRI scanning and clinical assessment, including the Movement Disorders Society Unified Parkinson’s Disease Rating Scale (MDS-UPDRS), MMSE, Montreal Cognitive Assessment (MoCA) and the revised Addenbrooke’s Cognitive Examination (ACE-R).

On the second and third sessions, patients completed a double-blind randomized placebo-controlled crossover study, with 40 mg of oral atomoxetine or placebo. Atomoxetine 40 mg was chosen as this is at the lower end of the clinical dose range and very well tolerated,, and previous studies have confirmed this dose is sufficient to modulate behaviour and functional brain activation in patients with Parkinson’s disease. Drug/placebo order was randomly permuted in groups of six successive recruits; that is, we randomized three participants to placebo and three to atomoxetine within each successive block of 1–6, 7–12, 13–18 and 19–24. This strategy guards against unequal randomization that can arise if one were to randomize across the entire group. The visits were ≥6 days apart [mean 7.4 days; standard deviation (SD) 1.7 days; range 6–14 days]. Blood samples were taken 2 h after administration of drug/placebo, to coincide with predicted peak plasma concentration of atomoxetine after a single oral dose. Mean plasma concentration was 261.32 ng/ml after atomoxetine (SD 117.33 ng/ml, range 90.92–595.11 ng/ml) and 0 ng/ml after placebo. After the blood sample, patients commenced an experimental task battery that included a stop-signal response inhibition task. Supine/lying and upright blood pressure and pulse rate measures were monitored three times across the session (on arrival, 2 h post tablet administration, and on completion of testing). To monitor any changes in subjective feelings following the drug/placebo, prior to tablet administration and 2 h post, we administered a set of 16 visual analogue scales (VAS) rating current mood and arousal levels. All sessions and MRI scanning were conducted with patients on their regular anti-parkinsonian medications and at a similar time of day.

Control participants were tested in one session to provide normative data on the task, in which they underwent MRI scanning and completed the same experimental task battery as the patients. The control group did not undergo the drug/placebo manipulation. Both the patient and control groups completed a set of self- and informant-rated questionnaires to assess mood and behaviour, which are reported in detail in the Supplementary material.

Stop-signal task

We used a stop-signal task to measure response inhibition. This paradigm involves a two-choice reaction time ‘go’ task that is occasionally interrupted by a ‘stop signal’, which requires the initiated response to be cancelled (Fig. 1A). On go trials, a left- or right-pointing black arrow was presented on the screen, and participants indicated its orientation by pressing a left or right response button. On stop trials, the arrow changed colour from black to red at the same time as a tone (i.e. the stop-signal), after a short and variable delay (i.e. the stop-signal delay, SSD). Participants were instructed to inhibit any imminent response if the arrow became red. The length of the SSD was varied across stop trials using a staircase method to target a stop accuracy of 50%. The SSD ranged from 50 ms to 1500 ms and increased or decreased by 50 ms after a successful or failed stop trial, respectively. On no-go trials, the SSD was set to zero. The task consisted of four blocks of 140 trials each, including 110 go trials, 10 no-go trials (very low commission error rate), and 20 stop trials per block (∼50% commission error rate). The first 20 trials of each block were go trials, to compute a starting value for the SSD (mean reaction time − 200 ms). The remaining trials within each block were pseudorandomly interleaved, with the constraints that there could be no more than seven consecutive go trials, and no more than two consecutive no-go or stop trials. At the start of each trial, a fixation cross was presented for 500 ms.

Figure 1

Design of the stop-signal task and ex-Gaussian race model of response inhibition. (A) In the stop-signal go/no-go task, participants respond as quickly and accurately as possible to the direction of a black arrow (go trials). Occasionally, this task is interrupted by a stop-signal (red arrow and beep tone), which requires any imminent response to be inhibited. For no-go trials, the stop-signal is presented immediately after the fixation cross. For stop trials, the stop-signal is presented after an initial go stimulus, with a short and variable delay. (B) The ex-Gaussian race model characterizes task performance as a race between three competing processes or ‘runners’: One stop process, and two go processes that match or mismatch the go stimulus. The finishing times of each process are assumed to follow an ex-Gaussian distribution. Successful inhibition in stop trials occurs when the stop process finishes before both go processes. A correct go response occurs when the matching go process finishes before the mismatching go process. For simplicity, the finishing time distribution of the mismatching go process is not illustrated. RT = reaction time. The speaker symbol in Fig. 1A was copied from the Twitter emoji library, available at https://twemoji.twitter.com/ under a CC-BY 4.0 license (https://creativecommons.org/licenses/by/4.0/). Figure 1B is from Heathcote et al., and is available at https://flic.kr/p/24g3sip under a CC-BY 2.0 license (https://creativecommons.org/licenses/by/2.0/).

Participants were given standardized instructions at the start of the experiment. They were asked to respond as ‘quickly and accurately as possible’, and were discouraged from strategically slowing down in anticipation of a stop signal. After these instructions, they were given a practice block of 25 trials (20 go trials, two no-go trials, and three stop trials). The experimenters verified that the participant understood the task, and if necessary, the practice block was repeated. The practice data were not analysed further. More extensive details about the task, conforming with the ‘check list for reporting stop-signal studies’ advised by Verbruggen et al., can be found in the Supplementary material.

Ex-Gaussian race model of response inhibition

We used a Bayesian parametric model of the stop-signal task to infer the latency of the unobservable stop response—the SSRT. This model assumes a race between three independent processes: one corresponding to the stop process, and two corresponding to go processes that match or mismatch the go stimulus (Fig. 1B). For a given stop trial, successful inhibition occurs when the stop process finishes before both go processes. For a given go trial, a correct response occurs when the matching go process finishes before the mismatching go process. The finish time distribution of the stop process is inferred by estimating the reaction time distribution of unsuccessful stop trials (i.e. signal respond reaction times). Specifically, the signal respond reaction time distribution is assumed to be a right-censored go reaction time distribution, where the censoring point for a given stop trial is drawn from the finish time distribution of the stop process (see Matzke et al. for details).

The model assumes that the finish times of the stop and go processes follow an ex-Gaussian distribution, which is a positively skewed unimodal distribution that is commonly used to describe reaction time data., Thus, for each process, we estimated the three parameters of the ex-Gaussian distribution: The mean (μ) and standard deviation (σ) of the Gaussian component, and the mean (i.e. inverse rate) (τ) of the exponential component.

We additionally estimated two parameters that represent the probability that the stop and go processes failed to start, referred to as ‘trigger failure’ and ‘go failure’, respectively. These attentional failures are common in both healthy participants, and in clinical cohorts,, and if not modelled can severely bias estimation of the stop process.,, Prior to fitting the model, we removed implausibly short (<0.25 s) or long (>4.5 s) reaction times, as well as go reaction times more extreme than ±2.5 SD from the participant’s mean.

We used Markov chain Monte Carlo (MCMC) sampling to estimate the posterior distributions of the parameters. The parameters were estimated hierarchically, such that parameters for a given participant were sampled from corresponding group-level distributions. We fitted this hierarchical model separately for the control group, the Parkinson’s disease group on placebo, and the Parkinson’s disease group on atomoxetine. We placed the same set of prior distributions on the group-level parameters for each of these three groups. The prior distributions were identical to those suggested by the model developers, except for slightly higher prior mean values for μgo-match (1.5 s), μgo-mismatch (1.5 s) and μstop (1 s), to account for slower reaction times in older age (see the Supplementary material for a full list of priors). The model MCMC sampling initially ran with 33 chains (i.e. three times the number of parameters), with thinning of every 10th sample and a 5% probability of migration. Model convergence was assessed with the potential scale reduction statistic (<1.1 for all parameters), and with visual inspection of the time-series plots of the chains. After this, an additional 500 iterations for each chain were run to create a final posterior distribution for each parameter. To assess the model’s goodness of fit, the observed data were compared to simulated data generated from the model’s posterior predictive distribution (Supplementary Figs 8–10).

The primary outcome of interest, SSRT, was computed as the mean of the ex-Gaussian finish time distribution of the stop process, which is given by μstop + τstop. We repeated this for each MCMC sample to approximate a posterior distribution of SSRT. This approach was also used to approximate a posterior distribution of go reaction time (μgo-match + τgo-match).

Statistical analysis

The go error rate and stop accuracy rate served as basic descriptive statistics for stop-signal task performance. We defined the go error rate as the proportion of go trials with an incorrect response, including commission errors (responses that mismatch the arrow direction) and omission errors (missing responses). The stop accuracy rate was defined as the proportion of stop trials with successfully inhibited responses (missing responses). For both outcomes, we examined differences between groups (Parkinson’s disease placebo versus controls) and drug conditions (Parkinson’s disease placebo versus atomoxetine) using independent and paired samples t-tests, respectively.

We then examined the group- and participant-level parameter estimates from the ex-Gaussian race model described above. For group-level inference, we examined the posterior distributions of the group-level means of SSRT and go reaction time. For a given posterior distribution, we took the median as the posterior estimate, and the 95% quantile interval (QI) as the range of most credible values. We also obtained posterior distributions for contrasts of interest (Parkinson’s disease placebo versus controls; Parkinson’s disease placebo versus atomoxetine) by subtracting the sets of MCMC samples of the two groups under consideration. That is, for a given parameter, we computed the difference between the two groups for each MCMC sample, thereby yielding an approximate posterior distribution of the difference.

To test for individual differences in the effect of atomoxetine on SSRT and go reaction time, we extracted the medians of the participant-level posterior distributions of SSRT and go reaction time from the placebo and atomoxetine model fits. We hypothesized that the effect of atomoxetine would depend on the integrity of the locus coeruleus, indexed by the contrast-to-noise ratio (CNR; described below). Therefore, the posterior parameter estimates were entered as the dependent variable with drug condition (placebo versus atomoxetine), CNR, and their interaction as fixed effects, allowing the intercept to vary across participants (random effect). We additionally included a fixed effect of session (first versus second) as a covariate of no interest. Taking the analysis of SSRT as an example, the model was specified in R formula syntax as follows: SSRT ∼ drug × CNR + session + (1 | subject).

For linear models, we report both frequentist and Bayes factor (BF) analyses for hypothesis testing, with a significance threshold of P = 0.05 (two-sided) for frequentist analyses. We present the BF for the alternative hypothesis over the null hypothesis (i.e. BF10), such that BF > 1 indicates relative evidence for the alternative hypothesis, and BF > 3 indicates ‘positive evidence’ for the alternative hypothesis. All BF analyses used the default ‘JZS’ prior on the effect size under the alternative hypothesis., To test for specific fixed effects in linear mixed models, we obtained P-values using the Kenward-Roger method, and BFs through Bayesian model averaging by estimating the change from prior to posterior inclusion odds (inclusion BF). In other words, this model-averaged BF indicates how much more likely the data are under model variants that include a given fixed effect, compared to model variants that exclude the fixed effect.

Because of technical issues, stop-signal task data were missing for one patient’s placebo session and for another patient’s atomoxetine session. We nevertheless included these two participants in the linear mixed model analyses of the within-subject effect of atomoxetine, as participants were treated as a random effect. However, excluding these two participants did not meaningfully change any of the following results.

Software and equipment

The stop-signal task was implemented in MATLAB R2018b using the Psychophysics Toolbox extensions (Version 3). Participants responded using a two-button response box. The ex-Gaussian model fitting was performed with the Dynamic Models of Choice toolbox, implemented in R (version 3.6.1, R Core Team, 2019). Further statistical analyses in R used the ‘tidyverse’ and ‘tidybayes’ packages for data organization and visualization, the ‘afex’ package for ANOVA and linear mixed model fitting with the ‘emmeans’ package used for post hoc comparisons, and the ‘BayesFactor’ and ‘bayestestR’ packages for BF analysis.

MRI acquisition

All patients and controls underwent MRI. Two controls were excluded from further imaging analysis because of incidental structural abnormalities. MRIs were acquired with a 7 T Magnetom Terra scanner (Siemens), using a 32-channel receive and circularly polarized single-channel transmit head coil (Nova Medical). We used a 3D high-resolution magnetization transfer-weighted turbo flash (MT-TFL) sequence for imaging the locus coeruleus (based on Priovoulos et al.). One hundred and twelve axial slices were used to cover both the midbrain and the pontine regions. The sequence applied a train of 20 Gaussian-shape RF-pulses at 6.72 ppm off resonance, 420° flip angle, followed by a turbo-flash readout (echo time = 4.08 ms, repetition time = 1251 ms, flip angle = 8°, voxel size = 0.4 × 0.4 × 0.5 mm3, 6/8 phase and slice partial Fourier, bandwidth = 140 Hz/px, no acceleration, 14.3% oversampling, acquisition time ∼7 min). For each subject, the transmit voltage was adjusted based on the average flip angle in the central area of the pons obtained from a B1 pre-calibration scan. The MT-TFL sequence was repeated twice and averaged offline to improve signal-to-noise ratio. An additional scan (MT-off) was acquired with the same parameters as above but without the off-resonance pulses. A high resolution isotropic T1-weighted MP2RAGE image was also acquired sagittally for anatomical coregistration using the UK7T Network harmonized protocol: echo time = 2.58 ms, repetition time = 3500 ms, bandwidth = 300 Hz/px, voxel size = 0.7 × 0.7 × 0.7 mm3, field of view = 224 × 224 × 157 mm3, acceleration factor (A ≫ P) = 3, flip angles = 5/2° and inversion times = 725/2150 ms for the first/second images.

Image processing and coregistration pipeline

Image processing and coregistration was based on the pipeline described in Ye et al. The Advanced Normalization Tools (ANTs v2.2.0) software and in-house MATLAB scripts were used for image preprocessing and the standardization of MT images. MT images were first N4 bias field corrected for spatial inhomogeneity (number of iterations at each resolution level: 50 × 50 × 30 × 20, convergence threshold: 1 × 10−6, isotropic sizing for b-spline fitting: 200) then averaged using the customized antsMultivariateTemplateConstruction2 function for improvements in signal-to-noise ratio. The T1-weighted MP2RAGE data were generated offline from the complex images. T1-weighted skull-stripped images were obtained after tissue type segmentation and reconstruction using SPM12 (v7219) (http://www.fil.ion.ucl.ac.uk/spm/software/spm12/ accessed 15 July 2021).

The preprocessed MT-weighted and T1-weighted images were then entered into a T1-driven, cross modality coregistration pipeline to warp the individual MT and MT-off images to the isotropic 0.5 mm ICBM152 (International Consortium for Brain Mapping) T1-weighted asymmetric template. The individual T1-weighted images were first coregistered to the MT image with rigid only transformation. The MT-off image was used as the intermediate step for bridging the two modalities because the MT-off image shares similar tissue-specific contrasts with both T1-weighted MP2RAGE and MT-on images.

In parallel, an unbiased study-wise T1-weighted structural template was created using individual skull-stripped T1-weighted images from all controls and patients. Native T1-weighted images were rigid and affine transformed, followed by a hierarchical non-linear diffeomorphic step at five levels of resolution, repeated by six runs to improve convergence. Max iterations for each resolution from the coarsest level to the full resolution were 100 × 100 × 70 × 50 × 20 (shrink factors: 10 × 6 × 4 × 2 × 1, smoothing factors: 5 × 3 × 2 × 1 × 0 voxels, gradient step size: 0.1 mm). Greedy Symmetric Normalization (SyN) was adopted for the transformation model of the deformation step. The resulting T1-weighted group template was then registered to the standard ICBM152 T1-weighted brain following the similar rigid-affine-SyN steps at four resolution levels (max iterations: 100 × 70 × 50 × 50, convergence threshold: 1 × 10−6, shrink factors: 8 × 4 × 2 × 1, smoothing factors: 3 × 2 × 1 × 0 voxels). For all the above registration steps, cross-correlation was used for similarity metrics estimation as it performs better for linear and non-linear components during intra-modality registration. Four steps of deformations were estimated as follows (in order): MT-off to MT, T1-weighted to MT-off, T1-weighted to T1-weighted group template and T1-weighted group template to ICBM152 T1-weighted template. These parameters were then used as the road map for MT image standardization to the ICBM brain in one step. A trilinear interpolation method was selected to preserve the absolute location and relative contrast of the signal.

Independent probabilistic locus coeruleus atlas creation

To facilitate accurate extraction of the locus coeruleus signal we created a study-specific unbiased locus coeruleus atlas (Fig. 2A). To this end, we used an independent sample of 29 age- and education-matched healthy control subjects [13 female; age mean (SD) = 67 (8.2), age range = 52–84] collected under the same neuroimaging protocol. We used a validated pipeline for locus coeruleus atlas construction described in Ye et al. Briefly, for each axial slice on the rostrocaudal extent, the locations of the left and right locus coeruleus were determined using a semi-automated segmentation method. A threshold was defined as 5 SD above the mean intensity in the central pontine reference region. After applying the threshold, locus coeruleus voxels on axial planes were automatically segmented into binarized images and then averaged to construct a probabilistic atlas. The independent locus coeruleus atlas generated for this study had very high similarity in the spatial distribution of probabilities and contours relative to the validated 7 T locus coeruleus atlas in Ye et al. (Supplementary Fig. 1).

Figure 2

Locus coeruleus imaging. (A) Study specific locus coeruleus atlas, also showing the reference region (light blue) in the central pons. (B) CNR for the locus coeruleus subdivisions and whole structure in Parkinson’s disease patients (PD) versus controls (note, left and right locus coeruleus are combined).

Locus coeruleus signal extraction

As a measure of locus coeruleus integrity, we quantified contrast by establishing the CNR with respect to a reference region in the central pons (Fig. 2A). A CNR map was computed voxel-by-voxel on the average MT image for each subject using the signal difference between a given voxel (V) and the mean intensity in the reference region (MeanREF) divided by the standard deviation (SDREF) of the reference signals:

CNR values were extracted bilaterally on the CNR map by applying the independent locus coeruleus probabilistic atlas (5% probability version). We computed mean CNR values for the rostral, middle and caudal portions of the left and right locus coeruleus. As an index of locus coeruleus integrity to incorporate with the stop-signal task analysis, we combined the left and right locus coeruleus and averaged across the whole structure.

Data availability

Code and data to reproduce manuscript figures, statistical analyses and stop-signal task modelling are freely available through the Open Science Framework (https://osf.io/tyka3/).

Results

Behavioural results

As shown in Table 1, the patient and control groups were matched in terms of age, years of education, sex ratio, MMSE and MoCA (all P-values > 0.050; all BFs < 1). Patients had a significantly lower ACE-R total score (P = 0.015, BF = 4.16), and lower memory subscale (P = 0.013, BF = 7.00). There were no group differences on any of the self-reported questionnaires assessing impulsivity, anxiety, depression and behavioural measures, with the exception of the Motivation and Energy Inventory where the patient group had significantly lower scores for the physical subscale compared to controls (P = 0.010, BF = 6.18; Supplementary Table 1). We did not see a bias in left versus right motor symptom laterality in the patient group; this analysis is described in detail in the Supplementary material.

Within the drug and placebo sessions, there was some evidence of increased pulse rates and raised blood pressure on atomoxetine, although this was not considered clinically relevant and was not consistently observed across all supine and upright measures. Importantly, there was no change in subjective ratings of mood and arousal levels within the sessions. These analyses are described in detail in the Supplementary material.

Locus coeruleus integrity

Figure 2B shows comparisons of locus coeruleus CNR between the patients and controls. When comparing across the whole structure, the groups were not significantly different [t(36.98) = 1.27, P = 0.21, BF = 0.58]. Comparing across the rostral, middle and caudal subdivisions, there was a main effect of subdivision [F(1.27,52.22) = 54.57, P <0.001; BF = 5.28 × 1010]. This was driven by CNR in the caudal portion being significantly lower than both the middle [t(82) = 7.37, P <0.0001] and rostral [t(82) = 10.10, P <0.0001] portions. There was a significant group by subdivision interaction [F(1.27,52.22) = 7.89, P = 004; BF = 32.98]. This reflected significantly lower CNR values in the caudal portion for patients relative to controls [t(50.1) = 2.23, P = 0.026], whereas the groups did not differ for the middle [t(50.1) = 1.32, P = 0.193] or rostral [t(50.1) = 0.140, P = 0.889] portions of the locus coeruleus. We did not see significant differences in left versus right locus coeruleus degeneration, as the group comparisons did not meaningfully change when incorporating locus coeruleus side; this analysis is reported in the Supplementary material. Additionally, we ran the comparisons of locus coeruleus CNR between the patients and controls using a more conservative 25% probability mask and the results were qualitatively identical; this analysis is reported in the Supplementary material.

Stop-signal task performance

In keeping with the tracking algorithm, the stop accuracy for the Parkinson’s disease group on placebo (Fig. 3A; mean = 0.48, SD = 0.15) was not significantly different from controls [mean = 0.58, SD = 0.18; t(40.79) = 1.90, P = 0.065, BF = 0.35]. Across the Parkinson’s disease patients, group-wise stop accuracy on atomoxetine (mean = 0.45, SD = 0.13) was not significantly different from the placebo session [t(16) = 0.88, P = 0.39, BF = 0.35]. The go error rate approached zero for most participants, yielding a skewed distribution bounded at zero (Fig. 3D), and was therefore logit transformed prior to analysis. The logit go error rate was slightly higher in the Parkinson’s disease group on placebo (mean = −4.29, SD = 1.23) than in controls [mean = −5.04, SD = 0.93; t(29.94) = −2.19, P = 0.037, BF = 0.46]. Within the Parkinson’s disease group, the logit go error rate was marginally reduced on atomoxetine (mean = −4.75, SD = 1.22) compared to placebo [t(16) = 2.23, P = 0.041, BF = 1.73]. However, we note that the BF for both these tests fell below conventional thresholds for positive evidence (i.e. BF > 3), and these effects on go error rate should therefore be regarded as ‘anecdotal’ at the group level.

Figure 3

Descriptive statistics and ex-Gaussian model estimates of stop-signal task performance. (A and D) Proportions of successful stop trials (A) and incorrect go responses (D). (B and E) Ex-Gaussian finish time distributions of the stop process (B) and matching go process (E). Bold lines represent group-level mean distributions; thin lines represent individual participants. The mean of a given ex-Gaussian distribution was taken as the SSRT (stop process) or go reaction time (RT) (matching go process). (C and F) Posterior distributions of group-level mean SSRT (C) and go reaction time (F). The black dots represent the medians; the thick black line segments represent the 66% quantile intervals; and the thin black line segments represent the 95% quantile intervals. PD = Parkinson’s disease.

Ex-Gaussian model estimates of stop-signal reaction time

The hierarchical Bayesian estimates of the ex-Gaussian finish time distributions for the stop and matching go processes are shown in Fig. 3B and E, respectively. The stop process finish times tended to be faster for the control group than the Parkinson’s disease group. Indeed, the posterior distribution of group-level mean SSRT (Fig. 3C) was lower for the control group (median = 0.39 s, 95% QI: 0.36, 0.41) than the Parkinson’s disease group on placebo (median = 0.46 s, 95% QI: 0.41, 0.52), and this group difference in SSRT was reliably different from zero (Δgroup median = 0.07 s, 95% QI: 0.01, 0.14). The mean SSRT for the Parkinson’s disease group on atomoxetine (median = 0.47 s, 95% QI: 0.42, 0.52) was comparable to the placebo session (Δdrug median = 0.01 s, 95% QI: −0.07, 0.09).

For the matching go process, the distributions of finish times varied widely across participants, but the group-level distributions were highly similar. The posterior distributions of group-level mean go eraction time (Fig. 3F) did not differ between the control group (median = 1.05 s, 95% QI: 0.28, 1.46), Parkinson’s disease group on placebo (median = 0.98 s, 95% QI: 0.37, 1.31; Δgroup median = −0.07 s, 95% QI: −0.78, 0.77), and Parkinson’s disease group on atomoxetine (median = 0.97 s, 95% QI: 0.30, 1.34; Δdrug median = −0.01 s, 95% QI: −0.72, 0.67). There were also no mean differences between groups or drug conditions for the attentional failure parameters, trigger failure and go failure (Supplementary Fig. 11).

Locus coeruleus integrity and atomoxetine-induced changes in stop-signal reaction time

Although there was no group-wise effect of atomoxetine on the Parkinson’s disease group in terms of their mean SSRT, we predicted that the effect of atomoxetine would depend on individual differences in locus coeruleus integrity, as indexed by the CNR. We confirmed a significant interaction effect between the drug condition and locus coeruleus CNR on the participant-level estimates of SSRT [Fig. 4A; β = 0.27, F(1,14.61) = 14.61, P = 0.002; BF = 11.56]. This interaction effect did not meaningfully change when including clinical covariates such as age, disease severity, atomoxetine plasma level and dopaminergic medication, as both frequentist and Bayesian model selection procedures indicated that such covariates did not significantly improve the model fit (for details see the Supplementary material). We performed several follow-up analyses that confirmed this interaction was not driven by overly influential outlier participants, including Cook’s distance,, DFBETA and refitting using a robust estimation method (for details see the Supplementary material). We also showed that the Drug × Locus coeruleus CNR interaction did not change when additional brain imaging metrics were included as covariates. These metrics included total intracranial volume and CNR extracted from the substantia nigra (details can be found in the Supplementary material). Finally, we found that the interaction did not change when using an alternative calculation for locus coeruleus contrast (contrast ratio to mean, not SD, see Supplementary material), or when the contrast values were obtained using the more conservative 25% probability mask (Supplementary material).

Figure 4

Relationship between SSRT, locus coeruleus and drug. (A) SSRT estimates as a function of drug condition and locus coeruleus CNR. Within-subject change in SSRT from placebo to atomoxetine is illustrated with vertical grey arrows. (B) Relationship between locus coeruleus CNR and the drug-induced change in SSRT. For visualization purposes, the SSRT estimates were adjusted for the fixed effect of session and random effect of participants (i.e. partial residuals).

There was also a main effect of session [β = 0.25, F(1,14.53) = 13.33, P = 0.002; BF = 2.58], reflecting slightly shorter SSRTs for the second session compared to the first, regardless of the drug condition. There were no significant main effects of drug condition [β = −0.05, F(1,14.21) = 0.70, P = 0.416; BF = 0.37] or locus coeruleus CNR [β = 0.11, F(1,14.61) = 0.266, P = 0.613; BF = 0.55] on SSRT.

To understand the Drug × Locus coeruleus CNR interaction further, we examined the relationship between locus coeruleus CNR and the drug-induced change in SSRT (Δdrug: atomoxetine − placebo), adjusted for the main effect of session. This relationship was strongly positive, suggesting that patients with lower locus coeruleus CNR have a greater reduction in SSRT after atomoxetine [Fig. 4B; r(15) = 0.73, P < 0.001; BF = 32.70].

There was no interaction effect between drug condition and locus coeruleus CNR on the participant-level estimates of go reaction time [β = 0.01, F(1,16.37) = 0.00, P = 0.971; BF = 0.47], stop accuracy rate [β = −0.18, F(1,16.52) = 1.15, P = 0.300; BF = 0.58], or logit go error rate [β = −0.02, F(1,15.14) = 0.03, P = 0.865; BF = 0.46]. Finally, to explore the possibility of hemisphere-specific effects, we ran the linear mixed model analysis separately for CNR values from the left and right locus coeruleus. These results were qualitatively identical to the original combined analysis, arguing against hemisphere-specific effects (Supplementary material).

Discussion

We show that improvements in response inhibition after atomoxetine are dependent on locus coeruleus integrity in patients with Parkinson’s disease. Following a single 40 mg dose of atomoxetine, individuals with lower locus coeruleus integrity had a greater improvement in response inhibition (i.e. reduction in their SSRT). This result highlights the link between integrity of the noradrenergic locus coeruleus and action cancellation, which has previously been inferred from preclinical work and pharmacological manipulations. The finding also demonstrates the importance of baseline noradrenergic capacity in determining the response to atomoxetine, confirming the need to stratify patients for noradrenergic therapy. Locus coeruleus neuromelanin imaging would be a safe and affordable means to achieve this stratification.

Previous work in Parkinson’s disease showed that atomoxetine led to greater improvements in SSRT and enhanced activation in the stopping network in patients with more severe disease., Extending this work, we confirm that improved SSRTs under atomoxetine occurred in those patients with more severe locus coeruleus degeneration. This was in the context of a general impairment in response inhibition, as patients on placebo had longer SSRTs compared to controls. Our result suggests that a single 40 mg atomoxetine dose confers the most benefit on individuals with a severe loss of noradrenergic capacity. In this way, noradrenergic replacement in patients with a compromised system may achieve restoration closer to normal levels and improve behaviour. Conversely, in patients with a less affected system the same dose may offer no appreciable benefit or even ‘overdose’ the system, leading to worse behaviour. This relationship is captured by the inverted U-shaped curve (known as a Yerkes-Dodson function) that is common across monoaminergic and cholinergic systems, whereby intermediate levels of neuromodulatory influence are associated with optimal performance, with too much or too little having deleterious effects on behaviour.,

Baseline dependency in dose-response curves is well documented in relation to dopaminergic therapy,, such that the optimal level of dopamine enhancement needed to improve behaviour depends on pre-existing dopamine levels. This has had critical implications for Parkinson’s disease therapy. Dopamine dosages titrated to restore levels in the severely depleted dorsal striatum and motor system circuitry effectively overdose the less affected ventral tegmental area, ventral striatum and associated limbic pathways, impairing aspects of learning and cognitive flexibility. Our result suggests a similar baseline dependency for noradrenergic therapy, where optimal dosages needed for atomoxetine therapy may depend on the extent of degeneration in the locus coeruleus. This has important implications for optimizing noradrenergic therapy in Parkinson’s disease, as patients could be stratified based on locus coeruleus integrity to inform appropriate dosages in clinical trials or personalized treatment.

Neuromelanin-sensitive MT imaging of the locus coeruleus represents a promising avenue to achieve this stratification. Although we have previously shown a relationship between disease severity (as measured by the UPDRS-III) and atomoxetine responsivity, such measures of motor function or disease duration may not be the most accurate measure of noradrenergic capacity. Whilst progressive degeneration of the locus coeruleus is expected over the disease course, this will vary widely across individuals and will reflect the disease phenotype. Neuropathological studies and neuromelanin imaging have shown more pronounced locus coeruleus degeneration in certain phenotypes, including those with cognitive impairment or dementia, depression, an akinetic-rigid syndrome and REM sleep behaviour disorder,, relative to patients at equivalent disease stages. In our results we note that including the MDS-UPDRS-III as an index of disease severity did not meaningfully improve the model fit for the interaction between the drug condition and locus coeruleus CNR. This highlights the added value of locus coeruleus imaging, above and beyond disease severity metrics, to explain variations in atomoxetine responsivity.

Our locus coeruleus imaging identified that the greatest difference between controls and patients was in the caudal portion of the locus coeruleus. This has not previously been identified in Parkinson’s disease locus coeruleus imaging with 3 T MRI, which limited analysis to the whole structure. This is because of the inherent limitations in 3 T MRI, where voxel sizes and slice thickness are larger and do not allow for the spatial resolution needed to examine locus coeruleus subdivisions. While some neuropathology studies have noted comparable cell loss across the rostral-caudal extent of the locus coeruleus,, others have reported more severe degenerative changes in the caudal segment. Previous studies have validated the relationship between neuromelanin-sensitive MRI signal intensity and histological measures of the density of neuromelanin-accumulating neurons and neuromelanin concentration. In light of these links, neuromelanin-sensitive MRI is considered a promising biomarker of Parkinson disease-related neurodegenerative processes.

Our study aimed for noradrenergic modulation of the SSRT. Atomoxetine increases extracellular noradrenaline levels across the brain via its actions at the noradrenaline transporter, in particular increasing levels in the prefrontal cortex by 3-fold., However, atomoxetine may also increase extracellular dopamine levels, particularly in prefrontal cortex. Because of the relative sparsity of dopamine transporters in the prefrontal cortex, a portion of dopamine uptake is mediated by the noradrenaline transporter., In this way, improved response inhibition under atomoxetine in Parkinson’s disease might potentially also reflect elevated prefrontal dopamine levels. The respective roles of noradrenergic versus dopaminergic modulation over the SSRT is not fully resolved. Some studies in rodents and humans have implicated a selective link between noradrenergic transmission and action cancellation: increasing dopamine selectively does not affect the SSRT. This work has been taken in concert with studies using agents that non-selectively change both dopamine and noradrenaline and do affect the SSRT. For example, haloperidol and methylphenidate are shown to modulate the SSRT, having both dopaminergic and non-dopaminergic effects. One interpretation of this body of work is that the effects of the non-selective drugs are primarily driven by their noradrenergic actions, given that selective dopamine manipulations do not reliably modulate SSRT. However, some studies have found modulation of the SSRT using relatively selective dopamine drugs, and dopamine receptor availability is associated with SSRT performance. The specificity of noradrenergic versus dopaminergic modulation of the SSRT is not fully elucidated. These two closely related neurotransmitters share a common biosynthetic pathway. They can be co-released from the same terminals and both modulate goal-directed behaviour. Although we propose that the principal effect of atomoxetine on SSRT is mediated via noradrenergic mechanisms, dopaminergic effects may also contribute.

In practical terms, a key question is whether ‘noradrenergic’ drugs offer a useful strategy to treat inhibitory control deficits. The association we have shown between locus coeruleus integrity and the change in SSRT under atomoxetine suggests a causal association between the noradrenergic system and action cancellation in Parkinson’s disease. Our post hoc supplementary analyses with brain volume and substantia nigra integrity showed these covariates did not influence that association. This indicates a degree of specificity for the locus coeruleus-noradrenergic system’s role in the SSRT in Parkinson’s disease. However, future cross-over studies using specific noradrenergic versus dopaminergic agents would be beneficial to resolve residual ambiguity. Our current results point towards atomoxetine as a useful strategy to treat inhibitory deficits in Parkinson’s disease, in the context of appropriate stratification using locus coeruleus integrity.

Our results do not speak directly to atomoxetine’s mechanism of action in Parkinson’s disease. However, convergent evidence indicates that atomoxetine may increase efficiency in brain networks mediating response inhibition, via actions at the prefrontal cortex and the locus coeruleus., In patients with Parkinson’s disease, atomoxetine has been shown to increase activity within and between regions of the stopping network, including the pre-supplementary motor area and inferior frontal gyrus., Locus coeruleus degeneration in Parkinson’s disease is accompanied by reduced noradrenaline levels in forebrain regions.,, As noradrenaline release facilitates reconfigurations of large-scale networks,, depletion of forebrain noradrenaline is likely to impact the rapid engagement of brain network activity that is necessary for successful action cancellation. In patients with greater locus coeruleus degeneration, which may be accompanied by decreased or dysfunctional modulation of prefrontal noradrenergic targets, efficiency of the stopping network may be reduced. Consequently, these patients show the greatest benefit from a drug that can increase levels of prefrontal noradrenaline and upregulate locus coeruleus function.

We would support the use of 40 mg atomoxetine in subsequent single-dose studies in Parkinson’s disease. Based on the current results, as well as previous clinical trials, and mechanistic studies, this dose is very well tolerated and capable of modulating neural systems and behaviour. Future clinical trials may consider a dose-ranging algorithm based on tolerability, plasma levels or behaviour, with the potential to increase beyond a 40 mg starter dose.,, Gradually increasing the dose over the course of treatment is consistent with animal models of atomoxetine, which confirm that increased levels of prefrontal extracellular noradrenaline following an acute dose are sustained over chronic administration, but not to the same extent. Increasing the dose over the course of treatment may be necessary to maintain equivalent behavioural benefits. Future clinical studies may also consider doses based on genetic stratification of patients, specifically, allelic variation in the noradrenaline transporter gene (SLC6A2) and the CYP2D6 enzyme gene. These are known to modulate clinical response to atomoxetine and absorption, bioavailability and half-life. Our results suggest that locus coeruleus imaging may also serve as a cost-effective means of stratifying patients’ atomoxetine dose in future clinical settings, mitigating the need for ongoing visits to titrate doses.

Our study was designed to test a mechanistic hypothesis around locus coeruleus integrity and atomoxetine responsivity; therefore, we did not focus on clinically relevant outcomes. Future longer-term clinical trials are needed to determine if atomoxetine improves cognition and behaviour in a way that benefits daily function, including carer and patient ratings. The link between response inhibition and atomoxetine suggests that an important phenotype to focus on will be those patients with impulsivity, even where this falls short of an impulse control disorder. We also note that locus coeruleus degeneration has been strongly linked with REM sleep behaviour disorder., Future studies could selectively recruit those patients with impulse control disorders and/or polysomnography-confirmed REM sleep behaviour disorder in order to establish the clinical benefits of atomoxetine in these stratified phenotypes.

Given the role of noradrenaline in cognition and behaviour, noradrenergic dysfunction contributes to cognitive deficits beyond action cancellation in Parkinson’s disease. Optimizing noradrenergic therapy therefore has potential to provide relief across a variety of non-motor symptoms. Our results confirm the potential for stratified noradrenergic therapy in Parkinson’s disease, whereby the efficacy of these drugs varies across individuals depending on their baseline noradrenergic state. Locus coeruleus neuromelanin imaging may offer a marker of noradrenergic capacity that can be used to stratify patients to optimise successful outcomes in trials of noradrenergic therapy, and ultimately inform a more personalized treatment approach.

Click here for additional data file.