The effect of spatial sampling on the resolution of the magnetostatic inverse problem.

In magnetoencephalography, linear minimum norm inverse methods are commonly employed when a solution with minimal a priori assumptions is desirable. These methods typically produce spatially extended inverse solutions, even when the generating source is focal. Various reasons have been proposed for this effect, including intrisic properties of the minimum norm solution, effects of regularization, noise, and limitations of the sensor array. In this work, we express the lead field in terms of the magnetostatic multipole expansion and develop the minimum-norm inverse in the multipole domain. We demonstrate the close relationship between numerical regularization and explicit suppression of spatial frequencies of the magnetic field. We show that the spatial sampling capabilities of the sensor array and regularization together determine the resolution of the inverse solution. For the purposes of stabilizing the inverse estimate, we propose the multipole transformation of the lead field as an alternative or complementary means to purely numerical regularization.

A minimum assumption approach to MEG sensor array design.

Objective.Our objective is to formulate the problem of the magnetoencephalographic (MEG) sensor array design as a well-posed engineering problem of accurately measuring the neuronal magnetic fields. This is in contrast to the traditional approach that formulates the sensor array design problem in terms of neurobiological interpretability the sensor array measurements.Approach.We use the vector spherical harmonics (VSH) formalism to define a figure-of-merit for an MEG sensor array. We start with an observation that, under certain reasonable assumptions, any array ofmperfectly noiseless sensors will attain exactly the same performance, regardless of the sensors' locations and orientations (with the exception of a negligible set of singularly bad sensor configurations). We proceed to the conclusion that under the aforementioned assumptions, the only difference between different array configurations is the effect of (sensor) noise on their performance. We then propose a figure-of-merit that quantifies, with a single number, how much the sensor array in question amplifies the sensor noise.Main results.We derive a formula for intuitively meaningful, yet mathematically rigorous figure-of-merit that summarizes how desirable a particular sensor array design is. We demonstrate that this figure-of-merit is well-behaved enough to be used as a cost function for a general-purpose nonlinear optimization methods such as simulated annealing. We also show that sensor array configurations obtained by such optimizations exhibit properties that are typically expected of 'high-quality' MEG sensor arrays, e.g. high channel information capacity.Significance.Our work paves the way toward designing better MEG sensor arrays by isolating the engineering problem of measuring the neuromagnetic fields out of the bigger problem of studying brain function through neuromagnetic measurements.

Cortical beta burst dynamics are altered in Parkinson's disease but normalized by deep brain stimulation.

Exaggerated subthalamic beta oscillatory activity and increased beta range cortico-subthalamic synchrony have crystallized as the electrophysiological hallmarks of Parkinson's disease. Beta oscillatory activity is not tonic but occurs in 'bursts' of transient amplitude increases. In Parkinson's disease, the characteristics of these bursts are altered especially in the basal ganglia. However, beta oscillatory dynamics at the cortical level and how they compare with healthy brain activity is less well studied. We used magnetoencephalography (MEG) to study sensorimotor cortical beta bursting and its modulation by subthalamic deep brain stimulation in Parkinson's disease patients and age-matched healthy controls. We show that the changes in beta bursting amplitude and duration typical of Parkinson's disease can also be observed in the sensorimotor cortex, and that they are modulated by chronic subthalamic deep brain stimulation, which, in turn, is reflected in improved motor function at the behavioural level. In addition to the changes in individual beta bursts, their timing relative to each other was altered in patients compared to controls: bursts were more clustered in untreated Parkinson's disease, occurring in 'bursts of bursts', and re-burst probability was higher for longer compared to shorter bursts. During active deep brain stimulation, the beta bursting in patients resembled healthy controls' data. In summary, both individual bursts' characteristics and burst patterning are affected in Parkinson's disease, and subthalamic deep brain stimulation normalizes some of these changes to resemble healthy controls' beta bursting activity, suggesting a non-invasive biomarker for patient and treatment follow-up.

Modulation of sensory cortical activity by deep brain stimulation in advanced Parkinson's disease.

Despite optimal oral drug treatment, about 90% of patients with Parkinson's disease develop motor fluctuation and dyskinesia within 5-10 years from the diagnosis. Moreover, the patients show non-motor symptoms in different sensory domains. Bilateral deep brain stimulation (DBS) applied to the subthalamic nucleus is considered the most effective treatment in advanced Parkinson's disease, and it has been suggested to affect sensorimotor modulation and relate to motor improvement in patients. However, observations on the relationship between sensorimotor activity and clinical improvement have remained sparse. Here, we studied the somatosensory evoked magnetic fields in 13 right-handed patients with advanced Parkinson's disease before and 7 months after stimulator implantation. Somatosensory processing was addressed with magnetoencephalography during alternated median nerve stimulation at both wrists. The strengths and the latencies of the ~60-ms responses at the contralateral primary somatosensory cortices were highly variable but detectable and reliably localized in all patients. The response strengths did not differ between preoperative and postoperative DBSON measurements. The change in the response strength between preoperative and postoperative condition in the dominant left hemisphere of our right-handed patients correlated with the alleviation of their motor symptoms (p = .04). However, the result did not survive correction for multiple comparisons. Magnetoencephalography appears an effective tool to explore non-motor effects in patients with Parkinson's disease, and it may help in understanding the neurophysiological basis of DBS. However, the high interindividual variability in the somatosensory responses and poor tolerability of DBSOFF condition warrants larger patient groups and measurements also in non-medicated patients.

Deep brain stimulation of subthalamic nucleus modulates cortical auditory processing in advanced Parkinson's Disease.

Deep brain stimulation (DBS) has proven its clinical efficacy in Parkinson's disease (PD), but its exact mechanisms and cortical effects continue to be unclear. Subthalamic (STN) DBS acutely modifies auditory evoked responses, but its long-term effect on auditory cortical processing remains ambiguous. We studied with magnetoencephalography the effect of long-term STN DBS on auditory processing in patients with advanced PD. DBS resulted in significantly increased contra-ipsilateral auditory response latency difference at ~100 ms after stimulus onset compared with preoperative state. The effect is likely due to normalization of neuronal asynchrony in the auditory pathways. The present results indicate that STN DBS in advanced PD patients has long-lasting effects on cortical areas outside those confined to motor processing. Whole-head magnetoencephalography provides a feasible tool to study motor and non-motor neural networks in PD, and to track possible changes related to cortical reorganization or plasticity induced by DBS.

A test of the effort equalization hypothesis in children with cerebral palsy who have an asymmetric gait.

Healthy people can walk nearly effortlessly thanks to their instinctively adaptive gait patterns that tend to minimize metabolic energy consumption. However, the economy of gait is severely impaired in many neurological disorders such as stroke or cerebral palsy (CP). Moreover, self-selected asymmetry of impaired gait does not seem to unequivocally coincide with the minimal energy cost, suggesting the presence of other adaptive origins. Here, we used hemiparetic CP gait as a model to test the hypothesis that pathological asymmetric gait patterns are chosen to equalize the relative muscle efforts between the affected and unaffected limbs. We determined the relative muscle efforts for the ankle and knee extensors by relating extensor joint moments during gait to maximum moments obtained from all-out hopping reference test. During asymmetric CP gait, the unaffected limb generated greater ankle (1.36±0.15 vs 1.17±0.16 Nm/kg, p = 0.002) and knee (0.74±0.33 vs 0.44±0.19 Nm/kg, p = 0.007) extensor moments compared with the affected limb. Similarly, the maximum moment generation capacity was greater in the unaffected limb versus the affected limb (ankle extensors: 1.81±0.39 Nm/kg vs 1.51±0.34 Nm/kg, p = 0.033; knee extensors: 1.83±0.37 Nm/kg vs 1.34±0.38 Nm/kg, p = 0.021) in our force reference test. As a consequence, no differences were found in the relative efforts between unaffected and affected limb ankle extensors (77±12% vs 80±16%, p = 0.69) and knee extensors (41±17% vs 38±23%, p = 0.54). In conclusion, asymmetric CP gait resulted in similar relative muscle efforts between affected and unaffected limbs. The tendency for effort equalization may thus be an important driver of self-selected gait asymmetry patterns, and consequently advantageous for preventing fatigue of the weaker affected side musculature.

Secondary somatosensory cortex evoked responses and 6-year neurodevelopmental outcome in extremely preterm children.

OBJECTIVE: We assessed in extremely preterm born (EPB) children whether secondary somatosensory cortex (SII) responses recorded with magnetoencephalography (MEG) at term-equivalent age (TEA) correlate with neurodevelopmental outcome at age 6 years. Secondly, we assessed whether SII responses differ between 6-year-old EPB and term-born (TB) children. METHODS: 39 EPB children underwent MEG with tactile stimulation at TEA. At age 6 years, 32 EPB and 26 TB children underwent MEG including a sensorimotor task requiring attention and motor inhibition. SII responses to tactile stimulation were modeled with equivalent current dipoles. Neurological outcome, motor competence, and general cognitive ability were prospectively evaluated at age 6 years. RESULTS: Unilaterally absent SII response at TEA was associated with abnormal motor competence in 6-year-old EPB children (p = 0.03). At age 6 years, SII responses were bilaterally detectable in most EPB (88%) and TB (92%) children (group comparison, p = 0.69). Motor inhibition was associated with decreased SII peak latencies in TB children, but EPB children lacked this effect (p = 0.02). CONCLUSIONS: Unilateral absence of an SII response at TEA predicted poorer motor outcome in EPB children. SIGNIFICANCE: Neurophysiological methods may provide new means for outcome prognostication in EPB children.

Localization of Sensorimotor Cortex Using Navigated Transcranial Magnetic Stimulation and Magnetoencephalography.

The mapping of the sensorimotor cortex gives information about the cortical motor and sensory functions. Typical mapping methods are navigated transcranial magnetic stimulation (TMS) and magnetoencephalography (MEG). The differences between these mapping methods are, however, not fully known. TMS center of gravities (CoGs), MEG somatosensory evoked fields (SEFs), corticomuscular coherence (CMC), and corticokinematic coherence (CKC) were mapped in ten healthy adults. TMS mapping was performed for first dorsal interosseous (FDI) and extensor carpi radialis (ECR) muscles. SEFs were induced by tactile stimulation of the index finger. CMC and CKC were determined as the coherence between MEG signals and the electromyography or accelerometer signals, respectively, during voluntary muscle activity. CMC was mapped during the activation of FDI and ECR muscles separately, whereas CKC was measured during the waving of the index finger at a rate of 3-4 Hz. The maximum CMC was found at beta frequency range, whereas maximum CKC was found at the movement frequency. The mean Euclidean distances between different localizations were within 20 mm. The smallest distance was found between TMS FDI and TMS ECR CoGs and longest between CMC FDI and CMC ECR sites. TMS-inferred localizations (CoGs) were less variable across participants than MEG-inferred localizations (CMC, CKC). On average, SEF locations were 8 mm lateral to the TMS CoGs (p < 0.01). No differences between hemispheres were found. Based on the results, TMS appears to be more viable than MEG in locating motor cortical areas.

Running in highly cushioned shoes increases leg stiffness and amplifies impact loading.

Running shoe cushioning has become a standard method for managing impact loading and consequent injuries due to running. However, despite decades of shoe technology developments and the fact that shoes have become increasingly cushioned, aimed to ease the impact on runners' legs, running injuries have not decreased. To better understand the shoe cushioning paradox, we examined impact loading and the spring-like mechanics of running in a conventional control running shoe and a highly cushioned maximalist shoe at two training speeds, 10 and 14.5 km/h. We found that highly cushioned maximalist shoes alter spring-like running mechanics and amplify rather than attenuate impact loading. This surprising outcome was more pronounced at fast running speed (14.5 km/h), where ground reaction force impact peak and loading rate were 10.7% and 12.3% greater, respectively, in the maximalist shoe compared to the conventional shoe, whereas only a slightly higher impact peak (6.4%) was found at the 10 km/h speed with the maximalist shoe. We attribute the greater impact loading with the maximalist shoes to stiffer leg during landing compared to that of running with the conventional shoes. These discoveries may explain why shoes with more cushioning do not protect against impact-related running injuries.

Helsinki VideoMEG Project: Augmenting magnetoencephalography with synchronized video recordings.

The primary goal of the Helsinki VideoMEG Project is to enable magnetoencephalography (MEG) practitioners to record and analyze the video of the subject during an MEG experiment jointly with the MEG data. The project provides: •Hardware assembly instructions and software for setting up video and audio recordings of the participant synchronized to MEG data acquisition.•Basic software tools for analyzing video and audio together with the MEG data. The resulting setup allows reliable recording of video and audio from the subject in various real-world usage scenarios. The Helsinki VideoMEG Project allowed successful establishment of video-MEG facilities in four different MEG laboratories in Finland, Sweden and the United States.

Spontaneous sensorimotor cortical activity is suppressed by deep brain stimulation in patients with advanced Parkinson's disease.

Advanced Parkinson's disease (PD) is characterized by an excessive oscillatory beta band activity in the subthalamic nucleus (STN). Deep brain stimulation (DBS) of STN alleviates motor symptoms in PD and suppresses the STN beta band activity. The effect of DBS on cortical sensorimotor activity is more ambiguous; both increases and decreases of beta band activity have been reported. Non-invasive studies with simultaneous DBS are problematic due to DBS-induced artifacts. We recorded magnetoencephalography (MEG) from 16 advanced PD patients with and without STN DBS during rest and wrist extension. The strong magnetic artifacts related to stimulation were removed by temporal signal space separation. MEG oscillatory activity at 5-25 Hz was suppressed during DBS in a widespread frontoparietal region, including the sensorimotor cortex identified by the cortico-muscular coherence. The strength of suppression did not correlate with clinical improvement. Our results indicate that alpha and beta band oscillations are suppressed at the frontoparietal cortex by STN DBS in PD.

Cortico-muscular coherence parallels coherence of postural tremor and MEG during static muscle contraction.

Corticokinematic coherence (CKC), i.e., coherence calculated between MEG and an accelerometer signal, recording movement kinematics, can be used for functional mapping of the sensorimotor cortex. Cortical sources of CKC, induced by both voluntary and passive movements, localize at the proximity of sensorimotor cortex. We tested the CKC during a static muscle contraction to compare it with simultaneously measured cortico-muscular coherence (CMC) estimated between MEG and surface EMG to study the role of postural tremor in CMC in ten healthy volunteers. CKC was detectable also during this static task. CKC and CMC spectra had similar power distributions, and sources of CMC and CKC colocalized at the cortex in close proximity of the central sulcus. During the static hold task, the accelerometer signal originates from the postural tremor. The similarities between CMC and CKC indicate that postural tremor is related to CMC in healthy subjects.

An Internet-Based Real-Time Audiovisual Link for Dual MEG Recordings.

HYPERSCANNING: Most neuroimaging studies of human social cognition have focused on brain activity of single subjects. More recently, "two-person neuroimaging" has been introduced, with simultaneous recordings of brain signals from two subjects involved in social interaction. These simultaneous "hyperscanning" recordings have already been carried out with a spectrum of neuroimaging modalities, such as functional magnetic resonance imaging (fMRI), electroencephalography (EEG), and functional near-infrared spectroscopy (fNIRS). DUAL MEG SETUP: We have recently developed a setup for simultaneous magnetoencephalographic (MEG) recordings of two subjects that communicate in real time over an audio link between two geographically separated MEG laboratories. Here we present an extended version of the setup, where we have added a video connection and replaced the telephone-landline-based link with an Internet connection. Our setup enabled transmission of video and audio streams between the sites with a one-way communication latency of about 130 ms. Our software that allows reproducing the setup is publicly available. VALIDATION: We demonstrate that the audiovisual Internet-based link can mediate real-time interaction between two subjects who try to mirror each others' hand movements that they can see via the video link. All the nine pairs were able to synchronize their behavior. In addition to the video, we captured the subjects' movements with accelerometers attached to their index fingers; we determined from these signals that the average synchronization accuracy was 215 ms. In one subject pair we demonstrate inter-subject coherence patterns of the MEG signals that peak over the sensorimotor areas contralateral to the hand used in the task.

Cortico-muscular coherence in advanced Parkinson's disease with deep brain stimulation.

OBJECTIVE: Cortico-muscular coherence (CMC) is thought to reflect the interplay between cortex and muscle in motor coordination. In Parkinson's disease (PD) patients, levodopa has been shown to enhance CMC. This study examined whether subthalamic nucleus (STN) deep brain stimulation (DBS) affects the CMC in advanced PD. METHODS: Magnetoencephalography (MEG) and electromyography (EMG) measurements were done simultaneously both with DBS on and off to determine the CMC during wrist extension. The spatiotemporal signal space separation (tSSS) was used for artifact suppression. RESULTS: CMC peaks between 13 and 25 Hz were found in 15 out of 19 patients. The effect of DBS on CMC was variable. Moreover, the correlation between CMC and motor performance was inconsistent; stronger CMC did not necessarily indicate better function albeit tremor and rigidity may diminish the CMC. Patients having CMC between 13 and 25 Hz had the best motor scores at the group level. CONCLUSIONS: DBS modifies the CMC in advanced PD with large interindividual variability. SIGNIFICANCE: DBS does not systematically modify CMC amplitude in advanced PD. The results suggest that some components of the CMC may be related to the therapeutic effects of DBS.

Improving MEG performance with additional tangential sensors.

Recently, the signal space separation (SSS) method, based on the multipole expansion of the magnetic field, has become increasingly important in magnetoencephalography (MEG). Theoretical arguments and simulations suggest that increasing the asymmetry of the MEG sensor array from the traditional, rather symmetric geometry can significantly improve the performance of the method. To test this concept, we first simulated addition of tangentially oriented standard sensor elements to the existing 306-channel Elekta Neuromag sensor array, and evaluated and optimized the performance of the new sensor configuration. Based on the simulation results, we then constructed a prototype device with 18 additional tangential triple-sensor elements and a total of 360 channels. The experimental results from the prototype are largely in agreement with the simulations. In application of the spatial SSS method, the 360-channel device shows an approximately 100% increase in software shielding capability, while residual reconstruction noise of evoked responses is decreased by 20%. Further, the new device eliminates the need for regularization while applying the SSS method. In conclusion, we have demonstrated in practice the benefit of reducing the symmetry of the MEG array, without the need for a complete redesign.

Validation of head movement correction and spatiotemporal signal space separation in magnetoencephalography.

OBJECTIVE: Our aim was to assess the effectiveness and reliability of spatiotemporal signal space separation (tSSS) and movement correction (MC) in magnetoencephalography (MEG) recordings disturbed by head movements and magnetized material on the head. METHODS: We recorded MEG from 20 healthy adults in stationary (reference) head position and during controlled head movements. Nearby magnetic interference sources were simulated by attaching magnetized particles on the subject's head. Auditory and somatosensory stimuli were presented. MC, tSSS and averaging were performed to obtain auditory (AEF) and somatosensory (SEF) evoked fields. Neuronal sources were modeled as equivalent current dipoles. MC was also validated by reconstructing signals generated by current dipoles in a phantom. RESULTS: After MC, the AEF and SEF responses recorded during intermittent head movements were similar in amplitude to the reference recordings and differed by 5-7mm in source location. The tSSS method removed artifacts due to the attached magnetized particles but did not affect the reference data. CONCLUSIONS: The methods are able to reliably recover MEG responses contaminated by movements and magnetic artifacts on the head. SIGNIFICANCE: The combination of tSSS and MC methods is especially useful in clinical measurements, where movements and magnetic disturbances are commonly present.

Somatomotor mu rhythm amplitude correlates with rigidity during deep brain stimulation in Parkinsonian patients.

OBJECTIVE: Parkinsonian patients have abnormal oscillatory activity within the basal ganglia-thalamocortical circuitry. Particularly, excessive beta band oscillations are thought to be associated with akinesia. We studied whether cortical spontaneous activity is modified by deep brain stimulation (DBS) in advanced Parkinson's disease and if the modifications are related to the clinical symptoms. METHODS: We studied the effects of bilateral electrical stimulation of subthalamic nucleus (STN) on cortical spontaneous activity by magnetoencephalography (MEG) in 11 Parkinsonian patients. The artifacts produced by DBS were suppressed by tSSS algorithm. RESULTS: During DBS, UPDRS (Unified Parkinson's Disease Rating Scale) rigidity scores correlated with 6-10 Hz and 12-20 Hz somatomotor source strengths when eyes were open. When DBS was off UPDRS action tremor scores correlated with pericentral 6-10 Hz and 21-30 Hz and occipital alpha source strengths when eyes open. Occipital alpha strength decreased during DBS when eyes closed. The peak frequency of occipital alpha rhythm correlated negatively with total UPDRS motor scores and with rigidity subscores, when eyes closed. CONCLUSION: STN DBS modulates brain oscillations both in alpha and beta bands and these oscillations reflect the clinical condition during DBS. SIGNIFICANCE: MEG combined with an appropriate artifact rejection method enables studies of DBS effects in Parkinson's disease and presumably also in the other emerging DBS indications.

Effects of DBS on auditory and somatosensory processing in Parkinson's disease.

Motor symptoms of Parkinson's disease (PD) can be relieved by deep brain stimulation (DBS). The mechanism of action of DBS is largely unclear. Magnetoencephalography (MEG) studies on DBS patients have been unfeasible because of strong magnetic artifacts. An artifact suppression method known as spatiotemporal signal space separation (tSSS) has mainly overcome these difficulties. We wanted to clarify whether tSSS enables noninvasive measurement of the modulation of cortical activity caused by DBS. We have studied auditory and somatosensory-evoked fields (AEFs and SEFs) of advanced PD patients with bilateral subthalamic nucleus (STN) DBS using MEG. AEFs were elicited by 1-kHz tones and SEFs by electrical pulses to the median nerve with DBS on and off. Data could be successfully acquired and analyzed from 12 out of 16 measured patients. The motor symptoms were significantly relieved by DBS, which clearly enhanced the ipsilateral auditory N100m responses in the right hemisphere. Contralateral N100m responses and somatosensory P60m responses also had a tendency to increase when bilateral DBS was on. MEG with tSSS offers a novel and powerful tool to investigate DBS modulation of the evoked cortical activity in PD with high temporal and spatial resolution. The results suggest that STN-DBS modulates auditory processing in advanced PD. Hum Brain Mapp, 2011. © 2010 Wiley-Liss, Inc.

Early dissociation of face and object processing: a magnetoencephalographic study.

The early dissociation in cortical responses to faces and objects was explored with magnetoencephalographic (MEG) recordings and source localization. To control for differences in the low-level stimulus features, which are known to modulate early brain responses, we created a novel set of stimuli so that their combinations did not have any differences in the visual-field location, spatial frequency, or luminance contrast. Differing responses to face and object (flower) stimuli were found at about 100 ms after stimulus onset in the occipital cortex. Our data also confirm that the brain response to a complex visual stimulus is not merely a sum of the responses to its constituent parts; the nonlinearity in the response was largest for meaningful stimuli.

Somatosensory evoked potentials and magnetic fields elicited by tactile stimulation of the hand during active and quiet sleep in newborns.

OBJECTIVE: Our objective was to characterize the effects of sleep stages on tactile somatosensory evoked responses in full-term newborns. METHODS: Somatosensory evoked potentials (SEPs) and magnetic fields (SEFs) to tactile stimulation of the tip of the index finger and/or thenar eminence were measured from 14 healthy newborns. The stimulus was a gentle tap produced by a moving membrane driven by an air-pressure pulse. RESULTS: SEPs and SEFs to tactile stimulation of the skin were similar in waveform and latency to SEPs known to be produced by electrical stimulation of the fingertip of neonates. The two most distinguishable positive deflections of SEPs, P1 and P2, within 300 ms of the stimulation, and their magnetic counterparts were clearly smaller in active compared to quiet sleep. CONCLUSIONS: Our study demonstrates for the first time that it is possible to record SEFs in neonates, and that clear late cortical somatosensory responses are produced by tactile stimulation. In addition, the effect of sleep stage on these responses indicates differences in the processing of the incoming information, at least in the somatosensory modality, in active and quiet sleep. SIGNIFICANCE: Tactile stimulation may be useful as a completely non-invasive technique for studying the physiology of the somatosensory system in neonates. Methodologically, since the effect of sleep stage is profound, one must carefully monitor the sleep stages in studies of event-related responses in newborns, or else this effect may confound the phenomena being studied.