Neurodesk: an accessible, flexible and portable data analysis environment for reproducible neuroimaging.

Neuroimaging research requires purpose-built analysis software, which is challenging to install and may produce different results across computing environments. The community-oriented, open-source Neurodesk platform ( https://www.neurodesk.org/ ) harnesses a comprehensive and growing suite of neuroimaging software containers. Neurodesk includes a browser-accessible virtual desktop, command-line interface and computational notebook compatibility, allowing for accessible, flexible, portable and fully reproducible neuroimaging analysis on personal workstations, high-performance computers and the cloud.

Localized estimation of electromagnetic sources underlying event-related fields using recurrent neural networks.

Objective. To use a recurrent neural network (RNN) to reconstruct neural activity responsible for generating noninvasively measured electromagnetic signals.Approach. Output weights of an RNN were fixed as the lead field matrix from volumetric source space computed using the boundary element method with co-registered structural magnetic resonance images and magnetoencephalography (MEG). Initially, the network was trained to minimise mean-squared-error loss between its outputs and MEG signals, causing activations in the penultimate layer to converge towards putative neural source activations. Subsequently, L1 regularisation was applied to the final hidden layer, and the model was fine-tuned, causing it to favour more focused activations. Estimated source signals were then obtained from the outputs of the last hidden layer. We developed and validated this approach with simulations before applying it to real MEG data, comparing performance with beamformers, minimum-norm estimate, and mixed-norm estimate source reconstruction methods.Main results. The proposed RNN method had higher output signal-to-noise ratios and comparable correlation and error between estimated and simulated sources. Reconstructed MEG signals were also equal or superior to the other methods regarding their similarity to ground-truth. When applied to MEG data recorded during an auditory roving oddball experiment, source signals estimated with the RNN were generally biophysically plausible and consistent with expectations from the literature.Significance. This work builds on recent developments of RNNs for modelling event-related neural responses by incorporating biophysical constraints from the forward model, thus taking a significant step towards greater biological realism and introducing the possibility of exploring how input manipulations may influence localised neural activity.

Filling your time: Auditory flutter alters perceived duration via stimulus-locked responses.

In the current article, we examined the flutter-duration illusion; the extension of perceived duration when an interval is filled with auditory flutter. Participants reproduced flutter-filled and empty durations while electrophysiological activity was recorded. As expected, participants over-produced durations when they were filled with auditory flutter rather than unfilled. Using multivariate pattern analysis, we found several differentiating patterns of neural activity while participants listened to either flutter-filled or empty intervals. However, in subsequent single trial analysis, only two of these clusters predicted perceived duration in the flutter condition; one occurring in line with the second click of the flutter, and one in line with the fourth click. We relate this finding to the N1P2 component and P3a component to timing initiation and arousal, respectively.

Investigating predictive coding in younger and older children using MEG and a multi-feature auditory oddball paradigm.

There is mounting evidence for predictive coding theory from computational, neuroimaging, and psychological research. However, there remains a lack of research exploring how predictive brain function develops across childhood. To address this gap, we used pediatric magnetoencephalography to record the evoked magnetic fields of 18 younger children (M = 4.1 years) and 19 older children (M = 6.2 years) as they listened to a 12-min auditory oddball paradigm. For each child, we computed a mismatch field "MMF": an electrophysiological component that is widely interpreted as a neural signature of predictive coding. At the sensor level, the older children showed significantly larger MMF amplitudes relative to the younger children. At the source level, the older children showed a significantly larger MMF amplitude in the right inferior frontal gyrus relative to the younger children, P < 0.05. No differences were found in 2 other key regions (right primary auditory cortex and right superior temporal gyrus) thought to be involved in mismatch generation. These findings support the idea that predictive brain function develops during childhood, with increasing involvement of the frontal cortex in response to prediction errors. These findings contribute to a deeper understanding of the brain function underpinning child cognitive development.

Neurodesk: An accessible, flexible, and portable data analysis environment for reproducible neuroimaging.

Neuroimaging data analysis often requires purpose-built software, which can be challenging to install and may produce different results across computing environments. Beyond being a roadblock to neuroscientists, these issues of accessibility and portability can hamper the reproducibility of neuroimaging data analysis pipelines. Here, we introduce the Neurodesk platform, which harnesses software containers to support a comprehensive and growing suite of neuroimaging software (https://www.neurodesk.org/). Neurodesk includes a browser-accessible virtual desktop environment and a command line interface, mediating access to containerized neuroimaging software libraries on various computing platforms, including personal and high-performance computers, cloud computing and Jupyter Notebooks. This community-oriented, open-source platform enables a paradigm shift for neuroimaging data analysis, allowing for accessible, flexible, fully reproducible, and portable data analysis pipelines.

Neural correlates of face perception modeled with a convolutional recurrent neural network.

Objective.Event-related potential (ERP) sensitivity to faces is predominantly characterized by an N170 peak that has greater amplitude and shorter latency when elicited by human faces than images of other objects. We aimed to develop a computational model of visual ERP generation to study this phenomenon which consisted of a three-dimensional convolutional neural network (CNN) connected to a recurrent neural network (RNN).Approach.The CNN provided image representation learning, complimenting sequence learning of the RNN for modeling visually-evoked potentials. We used open-access data from ERP Compendium of Open Resources and Experiments (40 subjects) to develop the model, generated synthetic images for simulating experiments with a generative adversarial network, then collected additional data (16 subjects) to validate predictions of these simulations. For modeling, visual stimuli presented during ERP experiments were represented as sequences of images (time x pixels). These were provided as inputs to the model. By filtering and pooling over spatial dimensions, the CNN transformed these inputs into sequences of vectors that were passed to the RNN. The ERP waveforms evoked by visual stimuli were provided to the RNN as labels for supervised learning. The whole model was trained end-to-end using data from the open-access dataset to reproduce ERP waveforms evoked by visual events.Main results.Cross-validation model outputs strongly correlated with open-access (r= 0.98) and validation study data (r= 0.78). Open-access and validation study data correlated similarly (r= 0.81). Some aspects of model behavior were consistent with neural recordings while others were not, suggesting promising albeit limited capacity for modeling the neurophysiology of face-sensitive ERP generation.Significance.The approach developed in this work is potentially of significant value for visual neuroscience research, where it may be adapted for multiple contexts to study computational relationships between visual stimuli and evoked neural activity.

Theta tACS impairs episodic memory more than tDCS.

Episodic memory deficits are a common consequence of aging and are associated with a number of neurodegenerative disorders (e.g., Alzheimer's disease). Given the importance of episodic memory, a great deal of research has investigated how we can improve memory performance. Transcranial electrical stimulation (TES) represents a promising tool for memory enhancement but the optimal stimulation parameters that reliably boost memory are yet to be determined. In our double-blind, randomised, sham-controlled study, 42 healthy adults (36 females; 23.3 ± 7.7 years of age) received anodal transcranial direct current stimulation (tDCS), theta transcranial alternating current stimulation (tACS) and sham stimulation during a list-learning task, over three separate sessions. Stimulation was applied over the left temporal lobe, as encoding and recall of information is typically associated with mesial temporal lobe structures (e.g., the hippocampus and entorhinal cortex). We measured word recall within each stimulation session, as well as the average number of intrusion and repetition errors. In terms of word recall, participants recalled fewer words during tDCS and tACS, compared to sham stimulation, and significantly fewer words recalled during tACS compared with tDCS. Significantly more memory errors were also made during tACS compared with sham stimulation. Overall, our findings suggest that TES has a deleterious effect on memory processes when applied to the left temporal lobe.

A Guided Tutorial on Modelling Human Event-Related Potentials with Recurrent Neural Networks.

In cognitive neuroscience research, computational models of event-related potentials (ERP) can provide a means of developing explanatory hypotheses for the observed waveforms. However, researchers trained in cognitive neurosciences may face technical challenges in implementing these models. This paper provides a tutorial on developing recurrent neural network (RNN) models of ERP waveforms in order to facilitate broader use of computational models in ERP research. To exemplify the RNN model usage, the P3 component evoked by target and non-target visual events, measured at channel Pz, is examined. Input representations of experimental events and corresponding ERP labels are used to optimize the RNN in a supervised learning paradigm. Linking one input representation with multiple ERP waveform labels, then optimizing the RNN to minimize mean-squared-error loss, causes the RNN output to approximate the grand-average ERP waveform. Behavior of the RNN can then be evaluated as a model of the computational principles underlying ERP generation. Aside from fitting such a model, the current tutorial will also demonstrate how to classify hidden units of the RNN by their temporal responses and characterize them using principal component analysis. Statistical hypothesis testing can also be applied to these data. This paper focuses on presenting the modelling approach and subsequent analysis of model outputs in a how-to format, using publicly available data and shared code. While relatively less emphasis is placed on specific interpretations of P3 response generation, the results initiate some interesting discussion points.

Automatic morpheme identification across development: Magnetoencephalography (MEG) evidence from fast periodic visual stimulation.

The present study combined magnetoencephalography (MEG) recordings with fast periodic visual stimulation (FPVS) to investigate automatic neural responses to morphemes in developing and skilled readers. Native English-speaking children (N = 17, grade 5-6) and adults (N = 28) were presented with rapid streams of base stimuli (6 Hz) interleaved periodically with oddballs (i.e., every fifth item, oddball stimulation frequency: 1.2 Hz). In a manipulation-check condition, tapping into word recognition, oddballs featured familiar words (e.g., roll) embedded in a stream of consonant strings (e.g., ktlq). In the experimental conditions, the contrast between oddball and base stimuli was manipulated in order to probe selective stem and suffix identification in morphologically structured pseudowords (e.g., stem + suffix pseudowords such as softity embedded in nonstem + suffix pseudowords such as trumess). Neural responses at the oddball frequency and harmonics were analyzed at the sensor level using non-parametric cluster-based permutation tests. As expected, results in the manipulation-check condition revealed a word-selective response reflected by a predominantly left-lateralized cluster that emerged over temporal, parietal, and occipital sensors in both children and adults. However, across the experimental conditions, results yielded a differential pattern of oddball responses in developing and skilled readers. Children displayed a significant response that emerged in a mostly central occipital cluster for the condition tracking stem identification in the presence of suffixes (e.g., softity vs. trumess). In contrast, adult participants showed a significant response that emerged in a cluster located in central and left occipital sensors for the condition tracking suffix identification in the presence of stems (e.g., softity vs. stopust). The present results suggest that while the morpheme identification system in Grade 5-6 children is not yet adult-like, it is sufficiently mature to automatically analyze the morphemic structure of novel letter strings. These findings are discussed in the context of theoretical accounts of morphological processing across reading development.

Comprehensive behavioural intervention for tics: a neurophysiological intervention.

BACKGROUND: Gilles de la Tourette Syndrome (GTS) is a childhood-onset neuropsychiatric disorder characterised by motor and vocal tics. While Comprehensive Behavioural Intervention for Tics (CBIT) is an effective, non-pharmacological treatment for patients with GTS, the underlying neurophysiological basis of this intervention has not been investigated. METHODS: To investigate the clinical effectiveness of CBIT in reducing tic severity in young people with GTS and explore neurophysiological mechanisms associated with clinical change. RESULTS: There was a significant overall improvement in tic severity of large effect size. The Cortical Silent Period (CSP) to motor evoked potential (MEP) ratio (CSP/MEP ratio) increased after the intervention with a small effect size. Other neurophysiological measures of inhibition were not significantly related to the change in tic severity. CONCLUSIONS: Alongside significant clinical improvements, these results suggest a role for motor cortical Gamma-aminobutyric acid (GABA)-ergic inhibitory circuitry in the neurophysiological changes underlying CBIT treatment. These findings need to be replicated in larger studies using control samples.

Cortical oscillatory dysrhythmias in visual snow syndrome: a magnetoencephalography study.

Visual snow refers to the persistent visual experience of static in the whole visual field of both eyes. It is often reported by patients with migraine and co-occurs with conditions such as tinnitus and tremor. The underlying pathophysiology of the condition is poorly understood. Previously, we hypothesized that visual snow syndrome may be characterized by disruptions to rhythmical activity within the visual system. To test this, data from 18 patients diagnosed with visual snow syndrome, and 16 matched controls, were acquired using magnetoencephalography. Participants were presented with visual grating stimuli, known to elicit decreases in alpha-band (8-13 Hz) power and increases in gamma-band power (40-70 Hz). Data were mapped to source-space using a beamformer. Across both groups, decreased alpha power and increased gamma power localized to early visual cortex. Data from the primary visual cortex were compared between groups. No differences were found in either alpha or gamma peak frequency or the magnitude of alpha power, p > 0.05. However, compared with controls, our visual snow syndrome cohort displayed significantly increased primary visual cortex gamma power, p = 0.035. This new electromagnetic finding concurs with previous functional MRI and PET findings, suggesting that in visual snow syndrome, the visual cortex is hyperexcitable. The coupling of alpha-phase to gamma amplitude within the primary visual cortex was also quantified. Compared with controls, the visual snow syndrome group had significantly reduced alpha-gamma phase-amplitude coupling, p < 0.05, indicating a potential excitation-inhibition imbalance in visual snow syndrome, as well as a potential disruption to top-down 'noise-cancellation' mechanisms. Overall, these results suggest that rhythmical brain activity in the primary visual cortex is both hyperexcitable and disorganized in visual snow syndrome, consistent with this being a condition of thalamocortical dysrhythmia.

Natural vs forced language switching: Free selection and consistent language use eliminate significant performance costs and cognitive demands in the brain.

Bilinguals are known to switch language spontaneously in everyday conversations, even if there are no external requirements to do so. However, in the laboratory setting, language control is often investigated using forced switching tasks, which result in significant performance costs. The present study assessed whether switching would be less costly when performed in a more natural fashion, and what factors might account for this. Mandarin-English bilinguals engaged in language switching under three different contexts with varied task demands. We examined two factors which may be characteristic of natural switching: (i) freedom of language selection; (ii) consistency of language used to name each item. Participants' brain activities were recorded using magnetoencephalography (MEG), along with behavioural measures of reaction speed and accuracy. The natural context (with both free selection and consistent language use for each item) produced better performance overall, showing reduced mixing cost and no significant switch cost. The neural effect of language mixing was also reversed in this context, suggesting that freely mixing two languages was easier than staying in a single language. Further, while switching in the forced context elicited increased brain activity in the right inferior frontal gyrus, this switch effect disappeared when the language used to name each item was consistent. Together, these findings demonstrate that the two factors above conjointly contribute to eliminating significant performance costs and cognitive demands associated with language switching and mixing. Such evidence aligns with lexical selection models which do not assume bilingual production to be inherently effortful.

Performance of Bimanual Finger Coordination Tasks in Speakers Who Stutter.

Stuttering is a neurodevelopmental speech disorder characterized by the symptoms of speech repetition, prolongation, and blocking. Stuttering-related dysfluency can be transiently alleviated by providing an external timing signal such as a metronome or the voice of another person. Therefore, the existence of a core motor timing deficit in stuttering has been speculated. If this is the case, then motoric behaviors other than speech should be disrupted in stuttering. This study examined motoric performance on four complex bimanual tasks in 37 adults who stutter and 31 fluent controls. Two tasks utilized bimanual rotation to examine motor dexterity, and two tasks used the bimanual mirror and parallel tapping movements to examine timing control ability. Video-based analyses were conducted to determine performance accuracy and speed. The results showed that individuals who stutter performed worse than fluent speakers on tapping tasks but not on bimanual rotation tasks. These results suggest stuttering is associated with timing control for general motor behavior.

Analysis methods for measuring passive auditory fNIRS responses generated by a block-design paradigm.

Significance: Functional near-infrared spectroscopy (fNIRS) is an increasingly popular tool in auditory research, but the range of analysis procedures employed across studies may complicate the interpretation of data. Aim: We aim to assess the impact of different analysis procedures on the morphology, detection, and lateralization of auditory responses in fNIRS. Specifically, we determine whether averaging or generalized linear model (GLM)-based analysis generates different experimental conclusions when applied to a block-protocol design. The impact of parameter selection of GLMs on detecting auditory-evoked responses was also quantified. Approach: 17 listeners were exposed to three commonly employed auditory stimuli: noise, speech, and silence. A block design, comprising sounds of 5 s duration and 10 to 20 s silent intervals, was employed. Results: Both analysis procedures generated similar response morphologies and amplitude estimates, and both indicated that responses to speech were significantly greater than to noise or silence. Neither approach indicated a significant effect of brain hemisphere on responses to speech. Methods to correct for systemic hemodynamic responses using short channels improved detection at the individual level. Conclusions: Consistent with theoretical considerations, simulations, and other experimental domains, GLM and averaging analyses generate the same group-level experimental conclusions. We release this dataset publicly for use in future development and optimization of algorithms.

Whole-Language and Item-Specific Inhibition in Bilingual Language Switching: The Role of Domain-General Inhibitory Control.

A prominent theory of bilingual speech production holds that appropriate language selection is achieved via inhibitory control. Such inhibition may operate on the whole-language and/or item-specific level. In this study, we examined these two levels of control in parallel, by introducing a novel element into the traditional cued language switching paradigm: half of the stimuli were univalent (each required naming in the same language every time it appeared), and the other half were bivalent (each required naming in different languages on different trials). Contrasting switch and stay trials provided an index for whole-language inhibition, while contrasting bivalent and univalent stimuli provided an index for item-specific inhibition. We then investigated the involvement of domain-general brain mechanisms in these two levels of language control. Neuroimaging studies report activation of the pre-supplementary motor area (pre-SMA), a key region in the executive control brain network, during language switching tasks. However, it is unclear whether or not the pre-SMA plays a causal role in language control, and at which level it exerts control. Using repetitive transcranial magnetic stimulation (TMS) to transiently disrupt the pre-SMA, we observed an essential role of this brain region in general speech execution, while evidence for its specific involvement in each level of inhibition remains inconclusive.

Motor neuroplasticity: A MEG-fMRI study of motor imagery and execution in healthy ageing.

Age-related decline in motor function is associated with over-activation of the sensorimotor circuitry. Using a multimodal MEG-fMRI paradigm, we investigated whether this neural over-recruitment in old age would be related to changes in movement-related beta desynchronization (MRBD), a correlate of the inhibitory neurotransmitter γ-aminobutyric acid (GABA), and whether it would characterize compensatory recruitment or a reduction in neural specialization (dedifferentiation). We used MEG to assess age-related changes in beta band oscillations in primary motor cortices, fMRI to localize age-related changes in brain activity, and the Finger Configuration Task to measure task performance during overt and covert motor processing: motor execution (ME) and motor imagery (MI). The results are threefold: first, showing age-related neuroplasticity during ME of older adults, compared to young adults, as evidenced by increased MRBD in motor cortices and over-recruitment of sensorimotor areas; second, showing similar age-related neuroplastic changes during MI; and finally, showing signs of dedifferentiation during ME in older adults whose performance negatively correlated with connectivity to bilateral primary motor cortex. Together, these findings demonstrate that elevated MRBD levels, reflecting greater GABAergic inhibitory activity, and over-activation of the sensorimotor network during ME may not be compensatory, but rather might reflect an age-related decline of the quality of the underlying neural signal.

Speech-Induced Suppression for Delayed Auditory Feedback in Adults Who Do and Do Not Stutter.

Speech-induced suppression is the normal, relative amplitude reduction of the auditory evoked potential for self-, compared to externally-generated, auditory stimulation. It remains controversial as to whether adults who stutter exhibit expected auditory modulation during speech; some studies have reported a significant difference between stuttering and fluent groups in speech-induced suppression during speech movement planning, while others have not. We compared auditory evoked potentials (N1 component) for auditory feedback arising from one's own voice (Speaking condition) with passive listening to a recording of one's own voice (Listening condition) in 24 normally-fluent speakers and 16 adults who stutter under various delayed auditory feedback (DAF) time conditions (100 ms, 200 ms, 500 ms, and 1,000 ms). We presented the participant's own voice with a delay, immediately after presenting it without a delay. Our working hypothesis was that the shorter the delay time, the more likely the delayed sound is perceived as self-generated. Therefore, shorter delay time conditions are proposed to result in relatively enhanced suppression of the auditory system. Results showed that in fluent speakers, the shorter the delay time, the more the auditory evoked potential in the Speaking condition tended to be suppressed. In the Listening condition, there was a larger evoked potential with shorter delay times. As a result, speech-induced suppression was only significant at the short delay time conditions of 100 and 200 ms. Adults who stutter did not show the opposing changes in the Speaking and Listening conditions seen in the fluent group. Although the evoked potential in the Listening condition tended to decrease as the delay time increased, that in the Speaking condition did not show a distinct trend, and there was a significant suppression only at 200 ms delay. For the 200 ms delay condition, speakers with more severe stuttering showed significantly greater speech-induced suppression than those with less severe stuttering. This preliminary study suggests our methods for investigating evoked potentials by presenting own voice with a delay may provide a clue as to the nature of auditory modulation in stuttering.

Neuro-dynamics of executive control in bilingual language switching: An MEG study.

Bilinguals have a remarkable ability to juggle two languages. A central question in the field is concerned with the control mechanisms that enable bilinguals to switch language with ease. Theoretical models and neuroimaging evidence suggest that a range of control processes are at play during language switching, and their underlying neural mechanisms are closely related to executive function. What remains unclear is when these control processes are engaged in language switching. In this study, we used magnetoencephalography (MEG) to examine the brain activity while unbalanced Mandarin-English bilinguals performed a digit-naming task with cued language switching. Following presentation of the language cue, an asymmetrical switch effect was observed in the left inferior frontal gyrus (IFG), where switch-related increase in evoked brain activity was larger for switching into the non-dominant language. Following presentation of the naming target, evoked brain activity in the right IFG was larger when naming was required in the non-dominant language compared to the dominant language. We conclude that control processes take place in two stages during language switching, with the left IFG resolving interference following cue presentation and the right IFG inhibiting competing labels following target presentation.

Local context effects in the magnitude-duration illusion: Size but not numerical value sequentially alters perceived duration.

Many aspects of an event can change perceived duration. A common example of this is the magnitude-duration illusion, in which a high magnitude (e.g. large or high value) stimulus will be perceived to last longer than a low magnitude stimulus. The effects of magnitude on perceived duration are normally considered in terms of global context effects; what is large depends on the stimuli used throughout the experiment. In the current article, we examine local context effects in the magnitude-duration illusion, how trial-by-trial changes in magnitude affect the subjective duration of an event. We performed two experiments in which numerical magnitude and stimulus size were varied within either the example phase or reproduction phase of a temporal reproduction task. We showed that in the current trial the combined value-size magnitude presented in the example phase affected subsequent reproductions, while the magnitude presented in the reproduction phase did not. The size magnitude presented in the reproduction phase also affected the reproduction in the following trial, such that a larger stimulus in the current reproduction phase resulted in shorter reproductions in the next reproduction phase. This indicates that low level stimulus properties (i.e. size) can act to contextualize subsequent stimulus properties, which in turn affect perceived duration. The findings of our experiments add local, low-level, context effects to the known modifiers of perceived duration, as well as provide evidence with regards to the role of magnitude in interval timing.