Abstract
Information processing speed (IPS) is a core cognitive process, highly relevant in eve-ryday-life and the most frequent and disabling cognitive symptom in patients with Relapsing Multiple Sclerosis (pwRMS). By using regionally-specific focalized transcranial direct current stimulation (tDCS) in healthy individuala and pwRMS, we provide causal evidence for superior parietal lobe (SPL) involvement in IPS and identfied a clinical predictor of tDCS-response in pwRMS. The study employed a registered, randomized, sham-tDCS controlled, three-way blind-ed, cross-over trial and a mixed-factors design with eight arms [between-subjects: group (pwRMS, controls); tDCS-polarity (excitatory, inhibitory); within-subjects: stimulation (active-, sham-tDCS)]. Concurrently with tDCS, participants completed a computerized version of the Symbol-Digit-Modalities-Test (SDMT), the current gold standard for quantifying IPS impairment in pwRMS. Bayesian modeling with generalized linear mixed models provided strong evidence for polarity specific modulation of IPS by SPL-tDCS and a double-dissociation of stimulation response in pwRMS and controls. Healthy individuals showed the canonical pattern of signifi-cantly enhanced and reduced processing speed during excitatory or inhibitory tDCS. In pwRMS, a reversed pattern was found and tDCS-response was predicted by baseline SDMT performance; i.e., more or less impaired patients benefited from inhibitory or excitatory tDCS, respectively. Our results provide direct causal evidence for SPL involvement in IPS in health and disease and suggest that the degree of IPS impairment in pwRMS reflects compensatory or dysfunctional neuroplastic processes that can be counteracted in a polarity specific way. Identified standardized transition scores for the effectiveness of excitatory or inhibitory tDCS will inform future individ-ually-tailored stimulation protocols in pwRMS.
One Sentence Summary Polarity-specific modulation of information processing speed in mul-tiple sclerosis.
INTRODUCTION
Multiple sclerosis (MS) is a chronic inflammatory disease resulting in progressive central nerv-ous system demyelination and neurodegeneration.1 The relapsing form of MS (RMS) is the lead-ing cause of non-traumatic neurologic disability in early adulthood worldwide and manifests as recurrent episodes of inflammation, followed by variable degrees of remission and progressive neurological dysfunction.2 Patients may suffer from physical and psychological symptoms, and up to 65% of patients present with cognitive impairment.3 Cognitive impairments predicts later disability, and is one of the main reasons for reduced work productivity and associated with re-duced quality-of-life.4–8
The most prevalent and disabling cognitive symptom in people with RMS (pwRMS) is reduced information processing speed (IPS).3,9 IPS is a key cognitive process that comprises the ability to identify, discriminate, integrate, and make decisions about incoming information, and to swiftly respond at the behavioral level. According to the tri-factor model of IPS, sensory, cognitive and motor speed components can be distinguished that are subserved by a brain network comprised of sensory, fronto-parietal and cerebellar regions.10 It is estimated that 27-51% of pwRMS patients suffer from IPS impairment, which has been linked to disease progression, but can also be the only symptom of a relapse.3
Drug treatment of the inflammatory neuropathology cannot directly target impaired IPS and there are currently no approved medications for treating cognitive symptoms in pwRMS.3,11 This high-lights the clinical relevance of exploring novel and evidence-based treatments that directly modu-late activity and plasticity in the neural networks subserving IPS in pwRMS.9 A promising ap-proach to achieve this is non-invasive transcranial direct current stimulation (tDCS), that admin-isters weak electrical currents via scalp attached electrodes to the brain, to modulate neural excit-ability and plasticity in the human brain.12 TDCS has an excellent safety profile and has success-fully been used to improve cognitive functions in healthy individuals and patients with neuropsy-chiatric disorders.13–19
To date, the vast majority of tDCS studies in pwRMS have targeted motor symptoms or fatigue and only two studies specifically investigated tDCS effects on IPS.20,21 These studies adminis-tered multisession tDCS to the prefrontal cortex, either as stand-alone treatment, or in combina-tion with cognitive training.22,23 Only the latter reported improved IPS immediately after the treatment relative to a control group, but no between-group effects six months later. These mixed effects of tDCS mirror those reported for other symptoms in pwRMS.20,21
Moreover, animal and human studies have demonstrated that behavioral and neural effects of multisession tDCS can be maintained for weeks or even months after the end of the intervention period.24–27 In this context, it is worth noting that reorganization of the functional brain networks supporting cognition in pwRMS is poorly understood and upregulation of fronto-parietal activity in pwRMS may reflect compensatory or dysfunctional neuroplasticity, depending on disease pro-gression.28–33 This makes it difficult to predict the neural effects of specific tDCS montages and unintended long-term stimulation effects in multisession contexts bear the risk for inducing mal-adaptive neuroplasticity. Hence, there is currently an urgent need for conducting carefully de-signed proof-of-concept tDCS studies (1) that investigate the potential of specific tDCS interven-tions to enhance IPS and (2) to identify clinically relevant predictors of stimulation response, prior to conducting time- and cost-intensive interventional trials, while minimizing the risk for the patients.
This was accomplished in the present study, by implementing a balanced and randomized cross-over trial that involved two stimulation sessions, either active or placebo (sham) tDCS. A single active tDCS session was chosen, because effects are completely reversible, thereby minimizing the risk for potential negative long-term effects.13 During each session, pwRMS and age- and sex matched healthy individuals completed a computerized version of the Symbol Digit Modalities Test (SDMT), which is the current gold standard to determine the degree of processing speed impairment in pwRMS.34 TDCS was administered bilaterally to the superior parietal lobe (SPL), based on previous functional and diffusion imaging studies, suggesting specific involvement of this region in the speed component of IPS.35,36 A focalized tDCS set-up was used to constrain the current delivery to the target regions, which has been shown to result in regionally- and task-specific neural modulation.37–41 By investigating polarity specific (i.e., anodal-excitatory-cathodal-inhibitory; AeCi) effects of focal tDCS, we aimed to establish causal involvement of the SPL in IPS. Inclusion of a matched healthy control group served to determine the degree of IPS impairment in pwRMS (i.e., during sham tDCS) and to investigate the question if the expected AeCi pattern of active tDCS in healthy controls can be replicated in the patient groups. Based on functional imaging studies suggesting that fronto-parietal activity may reflect compensatory or dysfunctional processes depending on disease progression, we also explored if tDCS-response in the patients can be predicted by the degree of IPS impairment.28–33
RESULTS
This prospective, randomized, three-way blinded (i.e., participants, staff conducting the experi-ment and evaluating outcomes), sham-tDCS controlled, cross-over trial was conducted in a mixed-factors design with eight arms, representing the between-subjects factors group (pwRMS, healthy controls; N=32/group) and tDCS polarity (anodal, cathodal), and the within-subjects fac-tor stimulation (active, sham tDCS). All participants completed a baseline phase, including com-prehensive neuropsychological assessments (Table 1), followed by the experimental cross-over phase (Figure 1A). During each phase, participants completed the computerized SDMT (Figure 1B+2) and received either active or sham SPL-tDCS (Figure 3). Response latency was defined as primary outcome (Clintrials.gov: NCT04667221); accuracy data was also analyzed to confirm compliance with task instructions to answer as fast and accurate as possible. Statistical analyses were conducted in a Bayesian framework and data were analyzed using generalized linear mixed models (GLMMs) with lognormal and binomial distribution for response latency and accuracy data. Three models were fitted for each analyses to identify the best fitting model for our data: (1) An intercept only model; (2) a covariates model to correct for potential session effects, motor slowing (via the Trail Making Test A, TMT-A) and baseline performance on the paper-pencil SDMT;42 (3) A full model comprising the covariates specified above, and variables of interest (i.e., tDCS polarity [anodal, cathodal], group [pwRMS, controls], and stimulation type [active, sham]), as well as their interactions. Models included a group-level subject intercept (random intercept) to correct for individual performance levels. To establish the model that best explained our data, models were compared using the widely applicable information criterion (WAIC), which is a fully Bayesian extension of the Akaike information criterion (AIC), and the estimated log pointwise predictive density of the WAIC (ELPD WAIC; see Methods for details).43 Bayesi-an hypothesis testing used evidence ratios, which considers the evidence for and against specific hypotheses; i.e., possible changes in performance due to repeated task exposure would be inves-tigated with the following formula: . For a detailed descrip-tion of the models and interpretation of evidence ratios see Methods & Materials.
Response accuracy
Response accuracy during the experimental SDMT task was near ceiling level across groups (M=98.23%, SD=13.2%); group means are shown in Table S1, indicating compliance to the in-structions to answer as fast as possible without making mistakes. WAIC revealed that the inter-cept only model performed best (Table S2), suggesting that differences in accuracy levels were not explained by the factor group or the stimulation conditions (Table S3).
TDCS effects on information processing speed
This analysis specifically addressed effects of the respective tDCS conditions on IPS in healthy individuals and pwRMS. Assessment of the response latency models showed that the full model including covariates and the effects of interest performed best, indicated by lower ELPD WAIC and WAIC values (Table S4). The full model for response latency provides strong evidence for a three-way interaction between the effects of interest: Group × stimulation type × polarity (β=-0.10, 95%CI=[−0.11, −0.08], evidence ratio=∞; Figure 4). In the control group, post-hoc tests showed that anodal stimulation decreased reaction times (β=-0.03, 95%CI=[−0.03, −0.02], evi-dence ratio=∞) and cathodal stimulation increased reaction times (β=0.04, 95%CI=[0.04, 0.05], evidence ratio=∞) compared to their respective sham conditions. However, in the patients, cathodal stimulation decreased response latencies (β=-0.02, 95%CI=[−0.02, −0.01], evidence ra-tio=4,443.44), while anodal stimulation increased response latencies (β=0.01, 95%CI=[0.01, 0.02], evidence ratio=929.23) relative to sham tDCS. Hence, we demonstrate a double dissocia-tion for effects of anodal vs. cathodal tDCS in pwRMS and healthy controls (Table S5).
Additionally, there was strong evidence that the covariates included in our model influenced la-tency in the experimental SDMT. Participants were faster during the second session suggesting a learning effect (β=-0.12, 95%CI=[−0.12,-0.11], evidence ratio=∞). Longer response latencies in the TMT-A were positively correlated with longer response times in the experimental SDMT (β=0.06, 95%CI=[0.01, 0.10], evidence ratio=129.72), i.e., the term corrected for motor slowing. Additionally, an increased number of correct responses in the baseline paper-pencil SDMT, i.e., an indirect measure of shorter response latencies during the baseline SDMT, was associated with shorter response latencies in the modified SDMT (β=-0.17, 95%CI=[−0.21, −0.12], evidence ra-tio=∞). Therefore, this covariate corrected for group differences in the baseline SDMT perfor-mance and the results provide strong evidence for a correspondence between the clinical, and experimental versions of the SDMT (i.e., construct validity).
Baseline IPS performance predicts tDCS response
To investigate if baseline performance on the clinical version of the SDMT predicted tDCS re-sponse, the full reaction time model was extended: Baseline performance on the paper pencil SDMT, group, stimulation polarity and stimulation type, as well as their interactions were added as population level effects. Additionally, we corrected for TMT-A performance and the session effect, i.e., performance increase due to repeated exposure to the task. Subjects were added as group-level intercepts. To ensure the validity of this analysis and enhance clinical relevance of potential outcomes, we considered two different predictors: (1) study specific raw-scores of the baseline SDMT in our sample, (2) age- and education corrected norms of the SDMT.
Overall, results were consistent for the study specific raw-scores of baseline SDMT performance and the analysis that used age-corrected norms (with better model fit for study specific data over the previous full model; ELPD difference=-255.27; SE difference=25.12; |-255.27| > 49.23).34
A four-way interaction between group, polarity, stimulation type, and baseline SDMT score was confirmed for the study-specific and norm-based analyses (study specific: β=0.06, 95%CI=[0.04, 0.08], evidence ratio=∞; age-corrected norms: β=0.06, 95%CI=[0.04, 0.08], evidence ratio=∞; Figure 5). This four-way interaction has to be interpreted with caution because the group factor and baseline SDMT scores are related, i.e., patients generally have lower SDMT scores than healthy controls. The latter is illustrated in the upper panels of Figure 5B (i.e., lack of z-scores lower than −1 in healthy controls). Nonetheless, healthy controls with higher baseline SDMT scores showed a stronger canonical AeCi response to tDCS, i.e., anodal stimulation reduced re-sponse latency and cathodal stimulation increased response latency, relative to sham stimulation. Because none of the healthy controls had z-scores lower than −1, only the canonical response was observed.
However, the four-way interaction also indicates that the moderating effect of baseline SDMT was larger in pwRMS. In addition, since we included patients with MS that had z-scores lower than −1, we observed a reversal of the canonical AeCi response, i.e., cathodal reduced response latency and anodal increased response latency (Tables S6+7). Notably, for standardized scores, the transition from beneficial to non-beneficial effects of anodal (z <-.58) or cathodal (z >-.70) was consistent across the patient groups.
We additionally repeated the analyses and only included the patient-subgroups, to investigate if the control group drove the interaction effect. Here, we found a three-way interaction of stimula-tion polarity, stimulation type and baseline SDMT score, indicating that this was not the case (Tables S8+9).
Adverse effects, mood and affect and blinding efficacy
Overall, participants tolerated the stimulation well, there was no evidence that potential changes in mood and affect biased the behavioral tDCS effects and participant blinding was successful. Details of the statistical models investigating adverse effects, mood and affect and blinding effi-cacy are described in the Methods and Materials section.
Self-reported adverse effects were mainly rated as “none” (81.73%) or “mild” (14.05%). “Mod-erate” and “strong” adverse effects were reported in 0.6-3.63% of participants (Table 2). Com-parison of adverse effects models favored the simple model assuming that the strength of adverse effects is similar in all experimental groups since the absolute ELPD difference was smaller than the standard error multiplied by 1.96 (ELPD difference=-2.4; SE difference=2.9; |-2.4| < 5.68). A post-hoc test showed a beta value for stimulation around 0, and a small increase of evidence from the prior model (β=0.03, 95%CI=[−0.18, 0.23], evidence ratio=9.44), suggesting that adverse ef-fects were comparable between active and sham conditions (Table S10).
Potential effects of tDCS on (positive and negative) mood and affect were assessed before and after each stimulation session using the Positive and Negative Affect Schedule (PANAS).44 Comparison of the PANAS models favored the simpler model, since the full model only yielded marginally better model criteria (ELPD difference=-0.5; SE difference=2.5; |-0.5| < 4.9). This suggests that PANAS change scores (i.e., pre-to post active vs. sham tDCS) were comparable in all experimental groups. Post-hoc tests showed weak evidence that the reported affect changes from pre-to post-stimulation is the same during sham and active stimulation in all groups (β=0.21, 95%CI=[−1.11, 1.55], evidence ratio=1.45). The full model and summary statistics (groups, valence) are reported in Tables S11+12.
Assessment of blinding showed that many participants reported to either “not know”, when they were stimulated (41.94%) or reported to be stimulated during the second session (37.10%), the latter suggesting a recency-effect. Model comparisons between the three models showed that the full model performed slightly better than the simple model (ELPD difference=-3.3; SE differ-ence=1.4; |-3.3| > 2.81). The full model, however, did not show a clear pattern for the response behavior of participants. Most pronounced evidence was found for a different response pattern in pwRMS. They were more likely to answer that they did not know when they were stimulated and less likely to answer that they think they were stimulated during the second session (β=-0.91, 95%CI=[−2.12, 0.28], evidence ratio=13.32; Figure S1A). Similarly, people receiving cathodal stimulation during the second session were more likely to report that they did not know when they were stimulated and less likely to report that they thought they were stimulated during the second session (β=-1.09, 95%CI=[−2.52, 0.36], evidence ratio=13.00; Figure S1B). However, the posterior distributions for all predictors were relatively wide, which does not allow for a clear interpretation of effects. Hence, it is concluded that blinding of participants was successful. For details see Table 3 and Table S13.
DISCUSSION
This study presents strong evidence for causal involvement of bilateral SPL in the speed compo-nent of IPS in health and disease. Moreover, by identifying a clinically relevant participant-dependent predictor of stimulation response in pwPMS (i.e., baseline clinical SDMT perfor-mance), we provide a framework for future individualized neurostimulation treatments of im-paired IPS in this population. Our main findings can be summarized as follows: (1) Relative to sham-tDCS, healthy individuals showed the expected polarity specific pattern of increased (in-hibitory tDCS) or decreased (excitatory tDCS) response latency during the experimental SDMT, which provides strong causal evidence for involvement of bilateral SPL in healthy individuals. (2) At the group level, the opposite pattern was found in the patient groups and there was strong evidence for a double dissociation of tDCS effects between pwRMS and age- and sex matched healthy controls. This emphasizes that stimulation effects observed in healthy individuals may not necessarily translate in a one-by-one fashion to clinical populations with structural and func-tional brain reorganization. (3) However, our results also demonstrate that the degree of baseline IPS impairment in the patient groups mediated stimulation response in a polarity specific way. While more impaired patients responded favorably to cathodal tDCS (transition at z=-0.70), less impaired patients improved more when receiving anodal tDCS (transition z=-.58). This result is in line with functional imaging studies in pwRMS (for review see28), that have suggested com-pensatory upregulation of metabolic brain activity in task-relevant brain networks in patients with less pronounced motor or cognitive impairment. In more impaired patients, enhanced activity may indicate dysfunctional processes, involving breakdown of functional network communica-tion or disinhibition. Notably, the predictive value of baseline SDMT impairment for stimulation response in pwRMS was confirmed for both the study specific distribution of baseline SDMT values and standard scores, which will be particularly valuable to inform clinical tDCS treatment decisions in the future.34 (4) Finally, the rigorous design of our study and the high level of exper-imental control ensured that our results cannot be explained by unspecific factors (e.g., unblinding of participants or staff, tDCS effects on mood or affect, order effects in the cross-over phase). In addition, adverse effects were minimal and comparable between groups and stimula-tion conditions, which highlights safety, tolerability and feasibility of our focalized tDCS ap-proach.
Over the last two decades, tDCS has been studied extensively for its potential to improve human brain function. The majority of this research has focused on motor functions (Stagg & Nitsche 2018 JECT), which has also been a major focus in pw(R)MS.20,21 However, there is increased evidence that higher cognitive functions, that are known to be supported by large-scale neural networks45 can be improved by tDCS. Notably, stimulation administered to key nodes of task-relevant cognitive networks has shown promise to enhance performance and brain function in healthy individuals and neuropsychiatric diseases.14,46, but also revealed substantial variability of tDCS response within and across studies. The underlying sources of this variability are thought to be multifactorial and can broadly be described as participant-dependent (e.g., skull and brain anatomy or functional network organization) and stimulation-dependent factors (e.g., tDCS tim-ing, duration, intensity or focality), as well as their interactions.12,47 Moreover, this research has highlighted that structural and functional brain organization due to normal aging or in patients with neuropsychiatric diseases can alter stimulation response relative to healthy young or age-matched reference populations, respectively.15,46 This emphasizes the necessity to conduct care-fully designed proof-of-concept studies to investigate the effects of specific tDCS protocols in clinical populations and to identify predictors of stimulation response, prior to implementation in time- and cost-intensive clinical trials. This is particularly relevant in pwRMS, where less than a handful of studies have investigated potential tDCS effects on cognition and yielded mixed re-sults.20,21
Consequently, the present study implemented a highly controlled, prospective, randomized, tri-ple-blinded, balanced cross-over trial, that involved assessment of polarity effects and employed a focalized tDCS set-up, aimed at confirming the causal role of the SPL in IPS in pwRMS and matched healthy controls. The relevance of the stimulation target was informed by imaging stud-ies in pwRMS that suggested involvement of bilateral SPL in the speed component of IPS, fMRI adaptations of the SDMT resembling the task used in the present study, and also evidence from lesion studies demonstrating specific involvement of SPL in the manipulation and re-arrangement of information in working memory, which is highly relevant for the SDMT.35,36,48,49 Moreover, the adapted experimental paradigm was designed to minimize dependence on fine-motor skills compared to the clinically used paper- and-pencil version and allowed increasing the number of trials to enhance statistical power to detect potential tDCS effects.34 Our results also provide strong evidence for a correspondence between the novel computerized and paper-pencil versions of the SDMT, which highlights construct validity of the experimental paradigm employed in this study. Using this task and overall research design, we were able to confirm causal involvement of the SPL in SDMT performance in healthy controls and pwRMS. In addition, we identified a clin-ically relevant predictor of stimulation response in the patient groups (clinical SDMT impair-ment), that may be valuable to guide future tDCS treatment studies.
The canonical anodal-excitatory, cathodal-inhibitory response in healthy individuals is in line with the proposed up- and downregulation of neural excitability frequently observed in the motor system.50 However, polarity specific modulation of behavioral performance has been less con-sistent for cognition. This is likely explained by higher redundancy within the neural networks supporting cognitive functions, rendering them less malleable to effects of inhibitory tDCS.51 Hence, our results contribute to the ongoing debate if tDCS can modulate cognition in a polarity specific way and strengthens the assumption that bilateral SPL is critically involved in IPS in healthy individuals.
In pwRMS, Bayesian modeling provided strong evidence for a double dissociation of tDCS ef-fects relative to healthy controls, which emphasizes that stimulation effects in pwRMS cannot easily be derived from outcomes in neurotypical individuals. Most strikingly, however, patients with more pronounced impairment responded more favorably to cathodal tDCS, while less im-paired patients showed a beneficial response to anodal tDCS. Overall, this polarity by impairment interaction in pwRMS is in line with imaging studies, suggesting a gradient of neural compensa-tion and dysfunction that depends on the level of impairment in specific tasks.28–33 Notably, simi-lar mechanisms have been suggested for normal aging, where increased task-related activity in fronto-parietal regions can effectively compensate for structural and functional brain changes to a certain degree of task difficulty, while this may not be possible with more pronounced neurodegeneration.52,53 Hence, facilitation of compensatory processes by anodal tDCS in earlier stages of the disease, and inhibition of dysfunctional processes later (e.g., increased neural noise due to reduced lateral inhibition or disconnection) are likely candidates to explain our findings in pwRMS. Importantly, our sample was comprised of mild-to-moderately impaired patients, sug-gesting that the transition between the two processes may happen relatively early. This is also in line with one of the few longitudinal fMRI studies in pwRMS that demonstrated an association of increased parietal activity over time with SDMT performance declines.54 However, in the future, this hypothesis could be investigated more directly by administering bilateral focal SPL-tDCS concurrently with fMRI adaptation of the SDMT, to confirm the proposed mechanisms.12,55
There are some limitations to this trial. First, our sample was relatively small and for feasibility reasons, we only recruited patients with mild-to-moderate levels of functional impairment. None-theless, even in this sample, we observed a statistically sound interaction between group, polarity and stimulation condition and a dissociation between effects of excitatory and inhibitory tDCS that was explained by baseline SDMT impairment. While this suggests that patients with higher levels of impairment and brain dysfunction may benefit from cathodal SPL-tDCS, this needs to be confirmed in the future. Moreover, the double dissociation of tDCS effects between healthy individuals and pwRMS likely reflects study specific distributions of (mild-to-moderate) baseline impairment in the patients and may not generalize to samples includingh more impaired patients. Nonetheless, the impairment profile of patients recruited for this study highlights that the transi-tion from compensatory to dysfunctional activity in the SPL may happen relatively early in the course of the disease and is therefore suited to inform future tDCS treatment decisions.
Second, because our proof-of-principle study did not involve concurrent fMRI, the proposed mechanisms by which anodal or cathodal tDCS improved SDMT performance remain to be de-termined. Importantly, the experimental SDMT used in this study was designed to require mini-mal adaptations for future tDCS-fMRI implementations, to investigate the neural mechanisms of impairment (during sham tDCS) and neuromodulation by active tDCS (anodal vs. cathodal).
Third, based on the limited number of studies that investigated tDCS effects on cognition in pwRMS and to prevent inducing maladaptive neuroplastic adaptation, we opted for a cross-over design with a single active tDCS session. Because effects of single tDCS sessions are transient and reversible, this assured that no long-term negative effects would be induced.50 Notably, more pronounced and longer lasting tDCS effects have been confirmed for learning compared to cogni-tive tasks and for multisession tDCS in combination with behavioral training or therapy.25,27,56–58 In this context, it is acknowledged that our study only provides proof-of-principle evidence and we abstain from claiming that the observed effects are clinically relevant. Nonetheless, the direc-tion of effects is plausible regarding imaging findings in pwRMS and theoretical accounts of tDCS effects. Thus, we provide a rationale and guidance for future intervention studies that com-bine traditional neurocognitive or gaming-based interventions aimed at improving IPS with tDCS, which may be suited to induce more pronounced and lasting effects on IPS than those ob-served in the present study.59,60
Overall, we provide strong evidence for causal involvement of SPL in IPS in health and disease. The double dissociation of tDCS effects in pwRMS and healthy controls and variability of stimu-lation effects within the patient groups highlights the need to tailor stimulation protocols based in relevant participant-dependent factors. Our own data suggest that performance on the clinical version of the SDMT may be suited to inform tDCS treatment decisions in pwRMS and that the degree of IPS impairment reflects compensatory or dysfunctional neuroplastic processes in the SPL, that can be modulated by tDCS in a polarity specific way and that the transition between compensation and dysfunctional breakdown of neural networks supporting cognition may happen early in the course of disease. Our findings will help to inform future clinical trials aimed at im-proving impaired IPS in a more sustained way by combining repeated tDCS sessions and behav-ioral interventions.
MATERIALS AND METHODS
This study was designed as a registered (NCT04667221), prospective, randomized, three-way blinded (i.e., participants, staff conducting the experiment and evaluating outcomes), sham-tDCS controlled, cross-over trial with a mixed-factors design and eight arms, representing the between-subjects factors group (pwRMS, healthy controls; N=32/group; Table 1) and tDCS polarity (an-odal, cathodal), and the within-subjects factor stimulation (active, sham tDCS). All participants completed baseline assessments, including comprehensive neuropsychological assessments, fol-lowed by the experimental cross-over phase during which the computerized SDMT was complet-ed twice.
Participants
Thirty-two pwRMS (male/female: 11/21; Mean±SD age: 46.9±11.3 years) were recruited via the outpatient clinic of the University Medicine Greifswald and self-help groups. Inclusion criteria comprised age≥18 years, a specialist confirmed diagnosis of RMS based on the revised McDon-ald criteria, no acute relapse and/or application of corticosteroids within the last three months prior to the experimental intervention, normal or corrected to normal vision, sufficient hand mo-tor function to respond by button press on a computer keyboard, and being a native German speaker.61 Exclusion criteria included other major medical or neuropsychiatric diseases (e.g., ma-jor depression) and contraindications for tDCS (e.g., metal implants, prior medical procedures involving head or spinal cord, head trauma with unconsciousness, history of epilepsy, convul-sions, seizures, migraine, current pregnancy).13 The sample size was based on previous cross-over tDCS trials from our group in neurological populations (Note: these trials only involved comparison of anodal and sham tDCS, and larger effects were expected for the between group comparison of anodal and cathodal tDCS).62,63
Thirty-two healthy individuals were recruited from the local community and retrospectively matched to individual patients on a 1:1 basis by age and sex (male/female: 11/21; Mean±SD age: 46.3±11.0).
All participants provided written informed consent prior to study inclusion, completed a demo-graphic questionnaire, a comprehensive neuropsychological baseline examination and were screened for depression and anxiety. Characteristics of patients and healthy controls, results of neuropsychological baseline assessments and additional questionnaires are detailed in Table 1.
The study was approved by the Medical Ethics Committee of the University Medicine Greifswald, conducted in accordance with the Helsinki Declaration and registered with ClinTrials.gov (NCT04667221).
After completing baseline assessments, patients were randomly assigned to the stimulation arms by a computer-generated list (i.e., 32 codes assigning participants to specific experimental proto-cols, including tDCS polarity, order and right/left coding of response buttons for the experimental task) and participated in the experimental cross-over phase of the study. The order of sham and active stimulation was counterbalanced across the groups, i.e., half of the participants received either active or sham tDCS during the first session (Figure 1A). Healthy controls received the same stimulation protocol as patients they were matched to. Two pwRMS in the anodal tDCS arm had to be excluded due to technical errors (i.e., software issues that prevented starting the experimental paradigm) during one of their stimulation sessions.
Experimental SDMT
To minimize reliance on fine motor skills required by the original paper- and-pencil SDMT, we developed a computerized version that resembled previous MRI adaptations with a forced-choice button press response mode.34,48 This also allowed creating a longer version of the task, to en-hance statistical power to detect potential stimulation effects. The experimental SDMT paradigm was presented using NBS Presentation® (Version 19.0, Neurobehavioral Systems, Inc., Berkeley, CA) During the experimental sessions, participants completed nine experimental blocks, each with a fixed duration of 90 sec (Figure 1B). Each block comprised nine unique and visually dis-tinct symbols and digits (1-9) and a legend showing nine symbol-digit combinations (Figures 2). During each trial, the legend was shown at the top of the screen; at the bottom of the screen, a single symbol-digit-combination was displayed (probe), which was either coherent or incoherent with the legend. Participants were asked to indicate coherence of the probe with the legend by button-presses with left and right index finger (“M” and “Y” keys) on a QWERTZ keyboard. Incoherent probes within blocks were created by systematically shifting the symbol-digit-combinations. For each coherent combination, eight unique shifted “incoherent” combinations were created (i.e., symbol and digit+1-8). During each task block, one shifted set was randomly chosen, and coherent and incoherent symbol-digit-combinations were pseudo-randomly present-ed, so that the same digit or symbol was not presented twice in a row. After all pairings were pre-sented, another shifted set was randomly chosen. This assured that approximately the same num-ber of coherent and incoherent trials were presented. The task was self-paced, and the next trial appeared immediately after each response. Participants were instructed to respond as quickly as possible, without making mistakes. No immediate feedback on performance was provided during or after each block. The maximum number of trials per block was 144, which was unattainable within the 90 sec block duration. Short breaks were interspersed in between blocks (max. 60 sec) and participants could proceed with the next block at their own pace during this time.
Transcranial direct current stimulation
Stimulation was administered with a Neuroelectrics Starstim8 direct current stimulator using Pistim electrodes (radius=1cm) and a 3×1 set-up (i.e., a central active electrode, and three sur-rounding reference electrodes). Circular inserts in an EEG-cap secured the electrodes on the head. Center electrodes were positioned bilaterally over the SPL based on 10-10 EEG coordinates (i.e., P1/P2). Return electrodes were positioned in a circular way around the anodes (CP2, P5, PO4, CP1, P6, PO3). This set-up was chosen during the planning stage of the project based on current flow simulations using standard SimNIBS parameters and a MNI152 standard brain, that demon-strated peak current intensities in the superior parietal cortex (BA7) and additional current flow to surrounding and deeper brain regions, including the precuneus (Figure 3).64,65 Stimulation com-menced briefly prior to the start of the experimental sessions and was applied with 1.5 mA for 20 minutes with a 20s ramp-up and - down at the beginning and the end (active tDCS), or was ramped-up and - down over 40s (sham tDCS). This procedure which has been shown to ensure participant blinding for focal set-ups.37,39,41 Staff administering the stimulation were blinded us-ing predefined stimulation codes that concealed the applied stimulation type.
An adverse effects questionnaire was completed by the participants after each stimulation ses-sion13 Potential unspecific effects of tDCS were assessed with the Positive and Negative Affect Schedule (PANAS) before and after each session.44 After the second experimental session, sub-jects were asked if they believed that active stimulation was administered during the first or the second tDCS session to assess blinding integrity. Potential responses included “I do not know”, “During the first session”, or “During the second session”.
Outcome measures
Response latency for correct responses was our predefined primary outcome measure (NCT04667221). We also inspected potential effects of tDCS on accuracy as secondary outcome. Response latency data were preprocessed in three steps: (1) All trials with incorrect responses were excluded. (2) Only trials with response latencies between 0.2 and 6s were included, in line with the upper limit that was used by previous studies that used similar designs.42,66,67 (3) Trials with response latencies within the interval “median±3×median absolute deviation”, computed on a subject and session basis, were included in the analysis. This approach is more robust than in-tervals based on means and standard deviations.68 Across sessions and stimulation groups, 5.65-7.06% of trials were excluded (Table S1).
Statistical Analyses
Statistical analyses were conducted in a Bayesian framework using R and the brms package.69,70 Data were analyzed using generalized linear mixed models (GLMMs) using a logit-link function for response accuracy, and a lognormal-link function for response latency analyses for binomial and skewed data distributions, respectively. The response latency model was truncated to match the range of possible response time outcomes fixed by our filtering procedure. Three models were fitted for each analyses to find the best fitting model for our data: (1) An intercept only model; (2) a covariates model to correct for potential session effects, motor slowing via the TMT-A and baseline performance on the paper-pencil SDMT;42 (3) A full model comprising the covariates as specified above, and variables of interest (i.e., stimulation polarity [anodal, cathodal], group [pwRMS, controls], and stimulation type [active, sham]), as well as their interactions. All models included a group-level subject intercept (random intercept) to correct for individual performance levels. To establish the models that best explained our data, models were compared using the widely applicable information criterion (WAIC), which is a fully Bayesian extension of the Akaike information criterion (AIC), and the estimated log pointwise predictive density of the WAIC (ELPD WAIC). Difference between ELPD WAIC scores larger than 1.96 times than the standard error of the standard error of the ELPD difference have been suggested to be meaning-ful.43 For each model, 3000 samples per chain with eight chains were drawn after warm-up, re-sulting in overall 24,000 draws using the Hamiltonian Monte Carlo Algorithm.71 Due to conver-gence issues, low effective sampling size and high rhat-values, reaction time models had to be run with 5000 samples after warm-up with eight chains, i.e., 40,000 draws. Furthermore, the ran-dom intercept of the intercept model required a soft sum to zero constraint to converge.72
An additional analysis investigated if the degree of IPS impairment in pwRMS (quantified by the clinical SDMT administered during the baseline assessments) was associated with stimulation response (i.e., reduced or increased response latency in the respective tDCS conditions). For that, the full reaction time model was extended: Baseline performance on the paper pencil SDMT, group, stimulation polarity and stimulation type, as well as their interactions were added as popu-lation level effects. Additionally, we corrected for TMT-A performance and the session effect.42 Subjects were added as group-level intercepts. To ensure the validity of this analysis and enhance clinical relevance of potential outcomes, we considered two different predictors: (1) study specif-ic raw scores of the baseline SDMT in our sample, (2) age- and education corrected norms of the SDMT.34
To analyze self-reported adverse effects and PANAS data, we built two GLMMs for each out-come. (1) A simple model that assumed that the effect of active stimulation compared to sham is the same in all experimental groups. (2) A more complex model that allowed the effect of active stimulation to be distinct, which allows modeling population specific vulnerability to adverse effects.
For the PANAS data, negative and positive items were summed to build a sum-score for each dimension (i.e., positive or negative valence). We fitted a hurdle Gaussian model, since negative affect included many sum-scores of 0, but responses were otherwise normally distributed.73 The full model included factors for stimulation polarity, group, stimulation type, time point (pre-, post tDCS), and valence (positive, negative), as well as all interactions. The simple model, however, only included interactions between stimulation type, subject group, time point and valence, and the interaction between active stimulation and time point. Additionally, we modeled the hurdle parameter separately for positive and negative affect, since zeros were only found in the measures of negative affect. The PANAS models included a group level intercept for subjects.
Adverse effects data were modeled with a cumulative model assuming a continuous latent varia-ble (strength of adverse effect) that was measured using categorical responses (“none”, “mild”, “moderate”, “strong”).74 The full model included stimulation polarity, stimulation type, and group, as well as all their interactions. The simple model only considered interactions between stimulation polarity and group. Both models included a population-level factor for sensation types (i.e., itching, pain, burning, heat, metallic taste, fatigue, other) and group-level intercepts for sensation types and subjects.
To investigate blinding of participants, we used categorical models to interrogate the participants’ responses to the question: “During which session do you think you were stimulated?” With cate-gorical models, intercepts represent the shift from one category to another, e.g., the transition from the answer “do not know” to “I think, I was stimulated during the first session”. Factors in the models shift the location of these transitions, so that more participants fall into a certain cate-gory. Again, we fitted multiple models to describe the data. A full model that included the session when stimulation was applied the group and the stimulation polarity, as well as their interactions. A simple model that only included the interaction of group and the stimulation polarity, and a simple effect for the active stimulation session. Thus, we were able to model if the correct re-sponse is associated with the participants’ response and consistency of this relationship across groups. Model comparisons for the adverse effects, PANAS, and blinding models were investi-gated with WAIC and ELPD WAIC as criteria for model comparisons.
Prior choices
In a Bayesian framework, priors must be chosen to represent expected values the model parame-ters may take before the current data was considered. Experimental data shifts these priors to rep-resent the expectations after considering data. Priors were chosen to be weakly informative, i.e., population-level parameters were normally distributed with a mean of 0 and standard deviation of 1, and group-level standard deviations were set to an exponential distribution with rate of 2 in all models. Note that the blinding assessment did not include group level intercept, since only one answer per participant was collected.
Hypothesis testing
Bayesian hypotheses testing used evidence ratios. For directed hypothesis, an evidence ratio is the number of draws in the direction of an effect relative to the number of draws in the opposite direction, i.e., the evidence ratio for a learning effect across sessions would be computed: . Thus, evidence ratios indicate the amount of evidence in favor of a hypothesis compared to the direct counter hypothesis. Evidence ratios were computed
for the model parameters and post-hoc comparisons. The theoretical range of evidence ratios is 0 to ∞, whereby evidence ratios of ∞ indicate that all posterior draws were in favor of the tested hypotheses. For point-hypotheses, i.e., tests for equality, evidence ratios indicate an in-or de-crease of evidence compared to the prior model.70 Simulations with linear and logistic regression models have found that evidence ratios of 19 and 39 correspond to one- and two-sided hypothesis tests with an alpha-level of 0.05.75 Note that we are unaware of formal comparisons between evi-dence ratios and p-values in generalized linear mixed models. Post-hoc comparisons were always comprised of within subjects comparisons of active session to the respective sham sessions. For accuracy and response latency models, differences between the sessions were tested, while ad-verse effects, blinding, and PANAS models tested equality of sessions.
Data availability
The experimental paradigm and data, as well as scripts for data analyses are available at open science framework (OSF).
Funding
This research was partly supported by the German Research Foundation (Research Unit 5429/1 (467143400), SFB 1315 (327654276, B03) and project grants ME 3161/3-1; ME 3161/5-1, ME 3161/6-1; FL 379/35-1; FL 379/34-1) and the European Fonds for Regional Development (GHS-20-0007).
Competing interests
MG received honoraria from Bayer, Biogen, Bristol Meyers Squibb, Johnson & Johnson, Merck, Novartis, Roche, Sanofi and Teva, and received research grants from Merck and Novartis. All unrelated to this project. The remaining authors report no competing interests.
List of Supplementary Materials
Fig. S1 for conditional effects of the blinding model
Table S1 for descriptive statistics per group for the accuracy data
Table S2 for the WAIC table including all fitted model for the accuracy analysis.
Table S3 for a summary of the best fitting accuracy model.
Table S4 for the WAIC table including all fitted model for the response latency analysis.
Table S5 for a summary of the best fitting response latency model.
Table S6 for a summary of the exploratory baseline SDMT model with raw SDMT scores.
Table S7 for a summary of the exploratory baseline SDMT model with SDMT z-scores.
Table S8 for a summary of the exploratory baseline SDMT model with raw SDMT scores only including patients.
Table S9 for a summary of the exploratory baseline SDMT model with SDMT z-scores only in-cluding patients.
Table S10 for a summary of the adverse effects model.
Table S11 for a descriptive summary of the PANAS scores.
Table S12 for a summary of the PANAS model.
Table S13 for a summary of the blinding model.
Table S14 for an overview of the removed trials per group for the response latency analysis.
Data Availability
The experimental paradigm and data, as well as scripts for data analyses are available at open science framework (OSF).