MEG source imaging method using fast L1 minimum-norm and its applications to signals with brain noise and human resting-state source amplitude images

Magnetoencephalography Language Mapping Using Auditory Memory Retrieval and Silent Repeating Task

PURPOSE

The study aims to (1) examine the spatiotemporal map of magnetoencephalography-evoked responses during an Auditory Memory Retrieval and Silent Repeating (AMRSR) task, and determine the hemispheric dominance for language, and (2) evaluate the accuracy of the AMRSR task in Wernicke and Broca area localization.

METHODS

In 30 patients with brain tumors and/or epilepsies, the AMRSR task was used to evoke magnetoencephalography responses. We applied Fast VEctor-based Spatial-Temporal Analyses with minimum L1-norm source imaging method to the magnetoencephalography responses for localizing the brain areas evoked by the AMRSR task.

RESULTS

The Fast-VEctor-based Spatial-Temporal Analysis found consistent activation in the posterior superior temporal gyrus around 300 to 500 ms, and another activation in the frontal cortex (pars opercularis and/or pars triangularis) around 600 to 900 ms, which were localized to the Wernicke area (BA 22) and Broca area (BA 44 and BA 45), respectively. The language-dominant hemispheric laterization elicited by the AMRSR task was comparable with the result from an Auditory Dichotic task result given to the same patient, with the exception that AMRSR is more sensitive on bilateral language laterization cases on finding the Wernicke and Broca areas.

CONCLUSIONS

For all patients who successfully finished the AMRSR task, Fast-VEctor-based Spatial-Temporal Analysis could establish accurate and robust localizations of Broca and Wernicke area and determine hemispheric dominance. For subjects with normal auditory functionality, the AMRSR paradigm evaluation showed significant promise in providing reliable assessments of cerebral language dominance and language network localization.

Case report: Neural timing deficits prevalent in developmental disorders, aging, and concussions remediated rapidly by movement discrimination exercises

Background

The substantial evidence that neural timing deficits are prevalent in developmental disorders, aging, and concussions resulting from a Traumatic Brain Injury (TBI) is presented.

Objective

When these timing deficits are remediated using low-level movement-discrimination training, then high-level cognitive skills, including reading, attention, processing speed, problem solving, and working memory improve rapidly and effectively.

Methods

In addition to the substantial evidence published previously, new evidence based on a neural correlate, MagnetoEncephalography physiological recordings, on an adult dyslexic, and neuropsychological tests on this dyslexic subject and an older adult were measured before and after 8-weeks of contrast sensitivity-based left-right movement-discrimination exercises were completed.

Results

The neuropsychological tests found large improvements in reading, selective and sustained attention, processing speed, working memory, and problem-solving skills, never before found after such a short period of training. Moreover, these improvements were found 4 years later for older adult. Substantial MEG signal increases in visual Motion, Attention, and Memory/Executive Control Networks were observed following training on contrast sensitivity-based left-right movement-discrimination. Improving the function of magnocells using figure/ground movement-discrimination at both low and high levels in dorsal stream: (1) improved both feedforward and feedback pathways to modulate attention by enhancing coupled theta/gamma and alpha/gamma oscillations, (2) is adaptive, and (3) incorporated cycles of feedback and reward at multiple levels.

Conclusion

What emerges from multiple studies is the essential role of timing deficits in the dorsal stream that are prevalent in developmental disorders like dyslexia, in aging, and following a TBI. Training visual dorsal stream function at low levels significantly improved high-level cognitive functions, including processing speed, selective and sustained attention, both auditory and visual working memory, problem solving, and reading fluency. A paradigm shift for treating cognitive impairments in developmental disorders, aging, and concussions is crucial. Remediating the neural timing deficits of low-level dorsal pathways, thereby improving both feedforward and feedback pathways, before cognitive exercises to improve specific cognitive skills provides the most rapid and effective methods to improve cognitive skills. Moreover, this adaptive training with substantial feedback shows cognitive transfer to tasks not trained on, significantly improving a person's quality of life rapidly and effectively.

EMG-projected MEG High-Resolution Source Imaging of Human Motor Execution: Brain-Muscle Coupling above Movement Frequencies

Magnetoencephalography (MEG) is a non-invasive functional imaging technique for pre-surgical mapping. However, movement-related MEG functional mapping of primary motor cortex (M1) has been challenging in presurgical patients with brain lesions and sensorimotor dysfunction due to the large numbers of trails needed to obtain adequate signal to noise. Moreover, it is not fully understood how effective the brain communication is with the muscles at frequencies above the movement frequency and its harmonics. We developed a novel Electromyography (EMG)-projected MEG source imaging technique for localizing M1 during ~1 minute recordings of left and right self-paced finger movements (~1 Hz). High-resolution MEG source images were obtained by projecting M1 activity towards the skin EMG signal without trial averaging. We studied delta (1-4 Hz), theta (4-7 Hz), alpha (8-12 Hz), beta (15-30 Hz), and gamma (30-90 Hz) bands in 13 healthy participants (26 datasets) and two presurgical patients with sensorimotor dysfunction. In healthy participants, EMG-projected MEG accurately localized M1 with high accuracy in delta (100.0%), theta (100.0%), and beta (76.9%) bands, but not alpha (34.6%) and gamma (0.0%) bands. Except for delta, all other frequency bands were above the movement frequency and its harmonics. In both presurgical patients, M1 activity in the affected hemisphere was also accurately localized, despite highly irregular EMG movement patterns in one patient. Altogether, our EMG-projected MEG imaging approach is highly accurate and feasible for M1 mapping in presurgical patients. The results also provide insight into movement related brain-muscle coupling above the movement frequency and its harmonics.

Assessing Pediatric Mild Traumatic Brain Injury and Its Recovery Using Resting-State Magnetoencephalography Source Magnitude Imaging and Machine Learning

The objectives of this machine-learning (ML) resting-state magnetoencephalography (rs-MEG) study involving children with mild traumatic brain injury (mTBI) and orthopedic injury (OI) controls were to define a neural injury signature of mTBI and to delineate the pattern(s) of neural injury that determine behavioral recovery. Children ages 8-15 years with mTBI (n = 59) and OI (n = 39) from consecutive admissions to an emergency department were studied prospectively for parent-rated post-concussion symptoms (PCS) at: 1) baseline (average of 3 weeks post-injury) to measure pre-injury symptoms and also concurrent symptoms; and 2) at 3-months post-injury. rs-MEG was conducted at the baseline assessment. The ML algorithm predicted cases of mTBI versus OI with sensitivity of 95.5 ± 1.6% and specificity of 90.2 ± 2.7% at 3-weeks post-injury for the combined delta-gamma frequencies. The sensitivity and specificity were significantly better (p < 0.0001) for the combined delta-gamma frequencies compared with the delta-only and gamma-only frequencies. There were also spatial differences in rs-MEG activity between mTBI and OI groups in both delta and gamma bands in frontal and temporal lobe, as well as more widespread differences in the brain. The ML algorithm accounted for 84.5% of the variance in predicting recovery measured by PCS changes between 3 weeks and 3 months post-injury in the mTBI group, and this was significantly lower (p < 10-4) in the OI group (65.6%). Frontal lobe pole (higher) gamma activity was significantly (p < 0.001) associated with (worse) PCS recovery exclusively in the mTBI group. These findings demonstrate a neural injury signature of pediatric mTBI and patterns of mTBI-induced neural injury related to behavioral recovery.

Magnetoencephalogram-based brain-computer interface for hand-gesture decoding using deep learning

Advancements in deep learning algorithms over the past decade have led to extensive developments in brain-computer interfaces (BCI). A promising imaging modality for BCI is magnetoencephalography (MEG), which is a non-invasive functional imaging technique. The present study developed a MEG sensor-based BCI neural network to decode Rock-Paper-scissors gestures (MEG-RPSnet). Unique preprocessing pipelines in tandem with convolutional neural network deep-learning models accurately classified gestures. On a single-trial basis, we found an average of 85.56% classification accuracy in 12 subjects. Our MEG-RPSnet model outperformed two state-of-the-art neural network architectures for electroencephalogram-based BCI as well as a traditional machine learning method, and demonstrated equivalent and/or better performance than machine learning methods that have employed invasive, electrocorticography-based BCI using the same task. In addition, MEG-RPSnet classification performance using an intra-subject approach outperformed a model that used a cross-subject approach. Remarkably, we also found that when using only central-parietal-occipital regional sensors or occipitotemporal regional sensors, the deep learning model achieved classification performances that were similar to the whole-brain sensor model. The MEG-RSPnet model also distinguished neuronal features of individual hand gestures with very good accuracy. Altogether, these results show that noninvasive MEG-based BCI applications hold promise for future BCI developments in hand-gesture decoding.

Two mechanisms facilitate regional independence between brain regions based on an examination of alpha-band activity in healthy control adult males

BACKGROUND

At rest, 8 to 12 Hz alpha rhythms are the dominant rhythm in the brain, with a common peak alpha frequency (PAF = the frequency at which alpha generators show maximum power) observed across brain regions. Although a common PAF across brain regions should result in high between-region connectivity, especially connectivity measures assessing the phase-similarity between alpha generators, high inter-regional alpha connectivity has not been observed. The present study was conducted as an initial step toward identifying mechanisms that allow brain regions to maintain functional independence in the presence of a common PAF.

METHODS

MEG data were obtained from 16 healthy control male adults (mean age = 24 years; range 21 to 30 years). A task requiring participants to alternate between a 10 s eyes-closed condition and a 5 s eyes-open condition was used to drive parietal-occipital alpha generators, with the 10 s eyes-closed condition eliciting large-amplitude alpha activity and thus providing alpha measures with good signal-to-noise ratio for source localization. Alpha source-space measures were obtained using Vector-based Spatial-Temporal Analysis using L1-minimum-norm. In each participant, the four strongest parietal-occipital alpha generators were identified. Connectivity between sources was assessed via a measure of phase-based connectivity called inter-site phase clustering (ISPC).

RESULTS

Intra-class correlations (ICC) showed very high similarity in the average PAF (=computed using all eyes-closed data) between the four alpha sources (ICC single measure = 0.88, p < 0.001). Despite a common average PAF, across participants, significant ISPC was often observed no more than that expected by chance. Examination of the alpha time course data indicated that low ISPC was often due to instantaneous changes in alpha phase (phase slips). ISPC analyses removing data with phase slips indicated that low ISPC was also due to slight continuous changes in the alpha frequency, with frequency drift more likely in non-significant than significant ISPC trials.

CONCLUSIONS

The present exploratory effort suggested two processes underlying the lack of observed inter-regional alpha phase coherence that may help maintain regional functional independence even in the presence of a common PAF. In particular, although the alpha generators were observed to oscillate at the same rate on average, across time each alpha generator oscillated a little slower or faster, and about every one and a half seconds an alpha generator abruptly lost the beat. Because of this, functional independence among alpha generators (and thus brain regions) was the rule rather than the exception. Studies replicating these processes that allow brain regions to maintain functional independence, using different source localization methods and in different conditions (e.g., a true resting state), are warranted.

IMPACT STATEMENT

Using source localization to measure parietal-occipital alpha generator activity, two properties that limit between-region alpha functional connectivity are proposed. In particular, a model of alpha generator activity is offered where via transient phase slips occurring approximately every 1.5 s, as well as slight non-stationarity in the alpha frequency, brain regions retain a common alpha frequency while also maintaining regional identity and presumably functionality. Findings also suggest novel markers for use in studies examining changes in alpha activity across maturation as well as in studies examining alpha activity in patient populations where alpha abnormalities have been reported.

Towards an objective evaluation of EEG/MEG source estimation methods - The linear approach

The spatial resolution of EEG/MEG source estimates, often described in terms of source leakage in the context of the inverse problem, poses constraints on the inferences that can be drawn from EEG/MEG source estimation results. Software packages for EEG/MEG data analysis offer a large choice of source estimation methods but few tools to experimental researchers for methods evaluation and comparison. Here, we describe a framework and tools for objective and intuitive resolution analysis of EEG/MEG source estimation based on linear systems analysis, and apply those to the most widely used distributed source estimation methods such as L2-minimum-norm estimation (L2-MNE) and linearly constrained minimum variance (LCMV) beamformers. Within this framework it is possible to define resolution metrics that define meaningful aspects of source estimation results (such as localization accuracy in terms of peak localization error, PLE, and spatial extent in terms of spatial deviation, SD) that are relevant to the task at hand and can easily be visualized. At the core of this framework is the resolution matrix, which describes the potential leakage from and into point sources (point-spread and cross-talk functions, or PSFs and CTFs, respectively). Importantly, for linear methods these functions allow generalizations to multiple sources or complex source distributions. This paper provides a tutorial-style introduction into linear EEG/MEG source estimation and resolution analysis aimed at experimental (rather than methods-oriented) researchers. We used this framework to demonstrate how L2-MNE-type as well as LCMV beamforming methods can be evaluated in practice using software tools that have only recently become available for routine use. Our novel methods comparison includes PLE and SD for a larger number of methods than in similar previous studies, such as unweighted, depth-weighted and normalized L2-MNE methods (including dSPM, sLORETA, eLORETA) and two LCMV beamformers. The results demonstrate that some methods can achieve low and even zero PLE for PSFs. However, their SD as well as both PLE and SD for CTFs are far less optimal for all methods, in particular for deep cortical areas. We hope that our paper will encourage EEG/MEG researchers to apply this approach to their own tasks at hand.

Assessment of Effective Network Connectivity among MEG None Contaminated Epileptic Transitory Events

Characterizing epileptogenic zones EZ (sources responsible of excessive discharges) would assist a neurologist during epilepsy diagnosis. Locating efficiently these abnormal sources among magnetoencephalography (MEG) biomarker is obtained by several inverse problem techniques. These techniques present different assumptions and particular epileptic network connectivity. Here, we proposed to evaluate performances of distributed inverse problem in defining EZ. First, we applied an advanced technique based on Singular Value Decomposition (SVD) to recover only pure transitory activities (interictal epileptiform discharges). We evaluated our technique's robustness in separation between transitory and ripples versus frequency range, transitory shapes, and signal to noise ratio on simulated data (depicting both epileptic biomarkers and respecting time series and spectral properties of realistic data). We validated our technique on MEG signal using detector precision on 5 patients. Then, we applied four methods of inverse problem to define cortical areas and neural generators of excessive discharges. We computed network connectivity of each technique. Then, we confronted obtained noninvasive networks to intracerebral EEG transitory network connectivity using nodes in common, connection strength, distance metrics between concordant nodes of MEG and IEEG, and average propagation delay. Coherent Maximum Entropy on the Mean (cMEM) proved a high matching between MEG network connectivity and IEEG based on distance between active sources, followed by Exact low-resolution brain electromagnetic tomography (eLORETA), Dynamical Statistical Parametric Mapping (dSPM), and Minimum norm estimation (MNE). Clinical performance was interesting for entire methods providing in an average of 73.5% of active sources detected in depth and seen in MEG, and vice versa, about 77.15% of active sources were detected from MEG and seen in IEEG. Investigated problem techniques succeed at least in finding one part of seizure onset zone. dSPM and eLORETA depict the highest connection strength among all techniques. Propagation delay varies in this range [18, 25]ms, knowing that eLORETA ensures the lowest propagation delay (18 ms) and the closet one to IEEG propagation delay.

MEG source imaging detects optogenetically-induced activity in cortical and subcortical networks

Magnetoencephalography measures neuromagnetic activity with high temporal, and theoretically, high spatial resolution. We developed an experimental platform combining MEG-compatible optogenetic techniques in nonhuman primates for use as a functional brain-mapping platform. Here we show localization of optogenetically evoked signals to known sources in the superficial arcuate sulcus of cortex and in CA3 of hippocampus at a resolution of 750 µm³. We detect activation in subcortical, thalamic, and extended temporal structures, conforming to known anatomical and functional brain networks associated with the respective sites of stimulation. This demonstrates that high-resolution localization of experimentally produced deep sources is possible within an intact brain. This approach is suitable for exploring causal relationships between discrete brain regions through precise optogenetic control and simultaneous whole brain MEG recording with high-resolution magnetic source imaging (MSI).

Feasibility of Reconstructing Source Functional Connectivity with Low-Density EEG

OBJECTIVES

Functional connectivity (FC) is increasingly used as target for neuromodulation and enhancement of performance. A reliable assessment of FC with electroencephalography (EEG) currently requires a laboratory environment with high-density montages and a long preparation time. This study investigated the feasibility of reconstructing source FC with a low-density EEG montage towards a usage in real life applications.

METHODS

Source FC was reconstructed with inverse solutions and quantified as node degree of absolute imaginary coherence in alpha frequencies. We used simulated coherent point sources as well as two real datasets to investigate the impact of electrode density (19 vs. 128 electrodes) and usage of template vs. individual MRI-based head models on localization accuracy. In addition, we checked whether low-density EEG is able to capture inter-individual variations in coherence strength.

RESULTS

In numerical simulations as well as real data, a reduction of the number of electrodes led to less reliable reconstructions of coherent sources and of coupling strength. Yet, when comparing different approaches to reconstructing FC from 19 electrodes, source FC obtained with beamformers outperformed sensor FC, FC computed after independent component analysis, and source FC obtained with sLORETA. In particular, only source FC based on beamformers was able to capture neural correlates of motor behavior.

CONCLUSION

Reconstructions of FC from low-density EEG is challenging, but may be feasible when using source reconstructions with beamformers.

Marked Increases in Resting-State MEG Gamma-Band Activity in Combat-Related Mild Traumatic Brain Injury

Combat-related mild traumatic brain injury (mTBI) is a leading cause of sustained impairments in military service members and veterans. Recent animal studies show that GABA-ergic parvalbumin-positive interneurons are susceptible to brain injury, with damage causing abnormal increases in spontaneous gamma-band (30-80 Hz) activity. We investigated spontaneous gamma activity in individuals with mTBI using high-resolution resting-state magnetoencephalography source imaging. Participants included 25 symptomatic individuals with chronic combat-related blast mTBI and 35 healthy controls with similar combat experiences. Compared with controls, gamma activity was markedly elevated in mTBI participants throughout frontal, parietal, temporal, and occipital cortices, whereas gamma activity was reduced in ventromedial prefrontal cortex. Across groups, greater gamma activity correlated with poorer performances on tests of executive functioning and visuospatial processing. Many neurocognitive associations, however, were partly driven by the higher incidence of mTBI participants with both higher gamma activity and poorer cognition, suggesting that expansive upregulation of gamma has negative repercussions for cognition particularly in mTBI. This is the first human study to demonstrate abnormal resting-state gamma activity in mTBI. These novel findings suggest the possibility that abnormal gamma activities may be a proxy for GABA-ergic interneuron dysfunction and a promising neuroimaging marker of insidious mild head injuries.

Magnetoencephalography Signal Processing, Forward Modeling, Magnetoencephalography Inverse Source Imaging, and Coherence Analysis

Magnetoencephalography (MEG) is a noninvasive functional imaging technique for the brain. MEG directly measures the magnetic signal due to neuronal activation in gray matter with high spatial localization accuracy. The first part of this article covers the overall concepts of MEG and the forward and inverse modeling techniques. It is followed by examples of analyzing evoked and resting-state MEG signals using a high-resolution MEG source imaging technique. Next, different techniques for connectivity and network analysis are reviewed with examples showing connectivity estimates from resting-state and epileptic activity.

Resting-state magnetoencephalography source magnitude imaging with deep-learning neural network for classification of symptomatic combat-related mild traumatic brain injury

Combat‐related mild traumatic brain injury (cmTBI) is a leading cause of sustained physical, cognitive, emotional, and behavioral disabilities in Veterans and active‐duty military personnel. Accurate diagnosis of cmTBI is challenging since the symptom spectrum is broad and conventional neuroimaging techniques are insensitive to the underlying neuropathology. The present study developed a novel deep‐learning neural network method, 3D‐MEGNET, and applied it to resting‐state magnetoencephalography (rs‐MEG) source‐magnitude imaging data from 59 symptomatic cmTBI individuals and 42 combat‐deployed healthy controls (HCs). Analytic models of individual frequency bands and all bands together were tested. The All‐frequency model, which combined delta‐theta (1–7 Hz), alpha (8–12 Hz), beta (15–30 Hz), and gamma (30–80 Hz) frequency bands, outperformed models based on individual bands. The optimized 3D‐MEGNET method distinguished cmTBI individuals from HCs with excellent sensitivity (99.9 ± 0.38%) and specificity (98.9 ± 1.54%). Receiver‐operator‐characteristic curve analysis showed that diagnostic accuracy was 0.99. The gamma and delta‐theta band models outperformed alpha and beta band models. Among cmTBI individuals, but not controls, hyper delta‐theta and gamma‐band activity correlated with lower performance on neuropsychological tests, whereas hypo alpha and beta‐band activity also correlated with lower neuropsychological test performance. This study provides an integrated framework for condensing large source‐imaging variable sets into optimal combinations of regions and frequencies with high diagnostic accuracy and cognitive relevance in cmTBI. The all‐frequency model offered more discriminative power than each frequency‐band model alone. This approach offers an effective path for optimal characterization of behaviorally relevant neuroimaging features in neurological and psychiatric disorders.

Neural representations of imagined speech revealed by frequency-tagged magnetoencephalography responses

Speech mental imagery is a quasi-perceptual experience that occurs in the absence of real speech stimulation. How imagined speech with higher-order structures such as words, phrases and sentences is rapidly organized and internally constructed remains elusive. To address this issue, subjects were tasked with imagining and perceiving poems along with a sequence of reference sounds with a presentation rate of 4 Hz while magnetoencephalography (MEG) recording was conducted. Giving that a sentence in a traditional Chinese poem is five syllables, a sentential rhythm was generated at a distinctive frequency of 0.8 Hz. Using the frequency tagging we concurrently tracked the neural processing timescale to the top-down generation of rhythmic constructs embedded in speech mental imagery and the bottom-up sensory-driven activity that were precisely tagged at the sentence-level rate of 0.8 Hz and a stimulus-level rate of 4 Hz, respectively. We found similar neural responses induced by the internal construction of sentences from syllables with both imagined and perceived poems and further revealed shared and distinct cohorts of cortical areas corresponding to the sentence-level rhythm in imagery and perception. This study supports the view of a common mechanism between imagery and perception by illustrating the neural representations of higher-order rhythmic structures embedded in imagined and perceived speech.

Enhanced Fast-VESTAL for Magnetoencephalography Source Imaging: From Theory to Clinical Application in Epilepsy

A novel magnetoencephalography source imaging approach called Fast Vector-based Spatio-Temporal Analysis (Fast-VESTAL) has been successfully applied in creating source images from evoked and resting-state data from both healthy subjects and individuals with neurological and/or psychiatric disorders, but its reconstructed source images may show false-positive activations, especially under low signal-to-noise ratio conditions. Here, to effectively reduce false-positive artifacts, we introduced an enhanced Fast-VESTAL (eFast-VESTAL) approach that adopts generalized second-order cone programming. We compared the spatiotemporal characteristics of the eFast-VESTAL approach to those of the popular distributed source approaches (e.g., the minimum L2-norm/ mixed-norm methods) using computer simulations and auditory experiments. More importantly, we applied eFast-VESTAL to the presurgical evaluation of epilepsy. Our results demonstrated that eFast-VESTAL exhibited a lower dipole localization error and/or a higher correlation coefficient (CC) between the estimated source time series and ground truth under various conditions of source waveforms. Experimentally, eFast-VESTAL displayed more focal activation maps and a higher CC between the raw and predicted sensor data in response to auditory stimulation. Notably, eFast-VESTAL was the most accurate method for noninvasively detecting the epileptic zones determined using more invasive stereo-electroencephalography in the comparison.

MEG Working Memory N-Back Task Reveals Functional Deficits in Combat-Related Mild Traumatic Brain Injury

Combat-related mild traumatic brain injury (mTBI) is a leading cause of sustained cognitive impairment in military service members and Veterans. However, the mechanism of persistent cognitive deficits including working memory (WM) dysfunction is not fully understood in mTBI. Few studies of WM deficits in mTBI have taken advantage of the temporal and frequency resolution afforded by electromagnetic measurements. Using magnetoencephalography (MEG) and an N-back WM task, we investigated functional abnormalities in combat-related mTBI. Study participants included 25 symptomatic active-duty service members or Veterans with combat-related mTBI and 20 healthy controls with similar combat experiences. MEG source-magnitude images were obtained for alpha (8-12 Hz), beta (15-30 Hz), gamma (30-90 Hz), and low-frequency (1-7 Hz) bands. Compared with healthy combat controls, mTBI participants showed increased MEG signals across frequency bands in frontal pole (FP), ventromedial prefrontal cortex, orbitofrontal cortex (OFC), and anterior dorsolateral prefrontal cortex (dlPFC), but decreased MEG signals in anterior cingulate cortex. Hyperactivations in FP, OFC, and anterior dlPFC were associated with slower reaction times. MEG activations in lateral FP also negatively correlated with performance on tests of letter sequencing, verbal fluency, and digit symbol coding. The profound hyperactivations from FP suggest that FP is particularly vulnerable to combat-related mTBI.

Cortical electrophysiological evidence for individual-specific temporal organization of brain functional networks

The human brain has been demonstrated to rapidly and continuously form and dissolve networks on a subsecond timescale, offering effective and essential substrates for cognitive processes. Understanding how the dynamic organization of brain functional networks on a subsecond level varies across individuals is, therefore, of great interest for personalized neuroscience. However, it remains unclear whether features of such rapid network organization are reliably unique and stable in single subjects and, therefore, can be used in characterizing individual networks. Here, we used two sets of 5‐min magnetoencephalography (MEG) resting data from 39 healthy subjects over two consecutive days and modeled the spontaneous brain activity as recurring networks fast shifting between each other in a coordinated manner. MEG cortical maps were obtained through source reconstruction using the beamformer method and subjects' temporal structure of recurring networks was obtained via the Hidden Markov Model. Individual organization of dynamic brain activity was quantified with the features of the network‐switching pattern (i.e., transition probability and mean interval time) and the time‐allocation mode (i.e., fractional occupancy and mean lifetime). Using these features, we were able to identify subjects from the group with significant accuracies (~40% on average in 0.5–48 Hz). Notably, the default mode network displayed a distinguishable pattern, being the least frequently visited network with the longest duration for each visit. Together, we provide initial evidence suggesting that the rapid dynamic temporal organization of brain networks achieved in electrophysiology is intrinsic and subject specific.

Magnetoencephalography Research in Pediatric Autism Spectrum Disorder

Magnetoencephalography (MEG) research indicates differences in neural brain measures in children with autism spectrum disorder (ASD) compared to typically developing (TD) children. As reviewed here, resting-state MEG exams are of interest as well as MEG paradigms that assess neural function across domains (e.g., auditory, resting state). To date, MEG research has primarily focused on group-level differences. Research is needed to explore whether MEG measures can predict, at the individual level, ASD diagnosis, prognosis (future severity), and response to therapy.

Identifying electrophysiological markers of autism spectrum disorder and schizophrenia against a backdrop of normal brain development

An examination of electroencephalographic and magnetoencephalographic studies demonstrates how age-related changes in brain neural function temporally constrain their use as diagnostic markers. A first example shows that, given maturational changes in the resting-state peak alpha frequency in typically developing children but not in children who have autism spectrum disorder (ASD), group differences in alpha-band activity characterize only a subset of children who have ASD. A second example, auditory encoding processes in schizophrenia, shows that the complication of normal age-related brain changes on detecting and interpreting group differences in neural activity is not specific to children. MRI studies reporting group differences in the rate of brain maturation demonstrate that a group difference in brain maturation may be a concern for all diagnostic brain markers. Attention to brain maturation is needed whether one takes a DSM-5 or a Research Domain Criteria approach to research. For example, although there is interest in cross-diagnostic studies comparing brain measures in ASD and schizophrenia, such studies are difficult given that measures are obtained in one group well after and in the other much closer to the onset of symptoms. In addition, given differences in brain activity among infants, toddlers, children, adolescents, and younger and older adults, creating tasks and research designs that produce interpretable findings across the life span and yet allow for development is difficult at best. To conclude, brain imaging findings show an effect of brain maturation on diagnostic markers separate from (and potentially difficult to distinguish from) effects of disease processes. Available research with large samples already provides direction about the age range(s) when diagnostic markers are most robust and informative.

Resting-State Magnetoencephalography Source Imaging Pilot Study in Children with Mild Traumatic Brain Injury

Mild traumatic brain injury (mTBI) accounts for the vast majority of all pediatric TBI. An important minority of children who have suffered an mTBI have enduring cognitive and emotional symptoms. However, the mechanisms of chronic symptoms in children with pediatric mTBI are not fully understood. This is in part due to the limited sensitivity of conventional neuroimaging technologies. The present study examined resting-state magnetoencephalography (rs-MEG) source images in 12 children who had mTBI and 12 age-matched control children. The rs-MEG exams were performed in children with mTBI 6 months after injury when they reported no clinically significant post-injury psychiatric changes and few if any somatic sensorimotor symptoms but did report cognitive symptoms. MEG source magnitude images were obtained for different frequency bands in alpha (8-12 Hz), beta (15-30 Hz), gamma (30-90 Hz), and low-frequency (1-7 Hz) bands. In contrast to the control participants, rs-MEG source imaging in the children with mTBI showed: 1) hyperactivity from the bilateral insular cortices in alpha, beta, and low-frequency bands, from the left amygdala in alpha band, and from the left precuneus in beta band; 2) hypoactivity from the bilateral dorsolateral prefrontal cortices (dlPFC) in alpha and beta bands, from the ventromedial prefrontal cortex (vmPFC) in beta band, from the ventrolateral prefrontal cortex (vlPFC) in gamma band, from the anterior cingulate cortex (ACC) in alpha band, and from the right precuneus in alpha band. The present study showed that MEG source imaging technique revealed abnormalities in the resting-state electromagnetic signals from the children with mTBI.

A compressed sensing approach to immobilized nanoparticle localization for superparamagnetic relaxometry

Superparamagnetic relaxometry (SPMR) exploits the unique magnetic properties of targeted superparamagnetic iron oxide nanoparticles (SPIOs) to detect small numbers of cancer cells. Reconstruction of the spatial distribution of cancer-bound nanoparticles requires solving an ill-posed inverse problem. The current method, multiple source analysis (MSA), uses a least-squares fit to determine the strength and location of a pre-determined number of magnetic dipoles. In this proof-of-concept study, we propose the application of a sparsity averaged reweighting algorithm (SARA) for volumetric reconstruction of immobilized nanoparticle distributions. We first calibrate the parameters that define the location of the sensors in the forward model of measurement physics. Using this optimized model, we evaluated the performance of the algorithms on various configurations of single and multiple point-source phantoms. We investigated the effect of the data fidelity parameter, voxel size, and iterative reweighting on the reconstruction produced by SARA. We found that the calibrated physics model can predict the detected field values within 5% of the measured data. When only a single source was present, both algorithms were able to detect as little as 0.5 µg of immobilized particles. However, when two sources were measured simultaneously, MSA failed to detect sources containing as much as 10 µg of particles, while SARA detected all of the sources containing at least 5 µg of particles. We show that a suitable data fidelity parameter can be selected objectively, and the total magnitude and location of a point source reconstructed by SARA is not sensitive to voxel size. Detection and localization of multiple small clusters of nanoparticles is a crucial step in SPMR-based diagnostic applications. Our algorithm overcomes the need to know the number of dipoles before reconstruction and improves the sensitivity of the reconstruction when multiple sources are present.

Emerging Approaches to Neurocircuits in PTSD and TBI: Imaging the Interplay of Neural and Emotional Trauma

Posttraumatic stress disorder (PTSD) and traumatic brain injury (TBI) commonly co-occur in general and military populations and have a number of overlapping symptoms. While research suggests that TBI is risk factor for PTSD and that PTSD may mediate TBI-related outcomes, the mechanisms of these relationships are not well understood. Neuroimaging may help elucidate patterns of neurocircuitry both specific and common to PTSD and TBI and thus help define the nature of their interaction, refine diagnostic classification, and may potentially yield opportunities for targeted treatments. In this review, we provide a summary of some of the most common and the most innovative neuroimaging approaches used to characterize the neural circuits associated with PTSD, TBI, and their comorbidity. We summarize the state of the science for each disorder and describe the few studies that have explicitly attempted to characterize the neural substrates of their shared and dissociable influence. While some promising targets in the medial frontal lobes exist, there is not currently a comprehensive understanding of the neurocircuitry mediating the interaction of PTSD and TBI. Future studies should exploit innovative neuroimaging approaches and longitudinal designs to specifically target the neural mechanisms driving PTSD-TBI-related outcomes.

Abnormal maturation of the resting-state peak alpha frequency in children with autism spectrum disorder

Age-related changes in resting-state (RS) neural rhythms in typically developing children (TDC) but not children with autism spectrum disorder (ASD) suggest that RS measures may be of clinical use in ASD only for certain ages. The study examined this issue via assessing RS peak alpha frequency (PAF), a measure previous studies, have indicated as abnormal in ASD. RS magnetoencephalographic (MEG) data were obtained from 141 TDC (6.13-17.70 years) and 204 ASD (6.07-17.93 years). A source model with 15 regional sources projected the raw MEG surface data into brain source space. PAF was identified in each participant from the source showing the largest amplitude alpha activity (7-13 Hz). Given sex differences in PAF in TDC (females > males) and relatively few females in both groups, group comparisons were conducted examining only male TDC (N = 121) and ASD (N = 183). Regressions showed significant group slope differences, with an age-related increase in PAF in TDC (R²  = 0.32) but not ASD (R²  = 0.01). Analyses examining male children below or above 10-years-old (median split) indicated group effects only in the younger TDC (8.90 Hz) and ASD (9.84 Hz; Cohen's d = 1.05). In the older ASD, a higher nonverbal IQ was associated with a higher PAF. In the younger TDC, a faster speed of processing was associated with a higher PAF. PAF as a marker for ASD depends on age, with a RS alpha marker of more interest in younger versus older children with ASD. Associations between PAF and cognitive ability were also found to be age and group specific.

Dynamic cognitive remediation for a Traumatic Brain Injury (TBI) significantly improves attention, working memory, processing speed, and reading fluency

BACKGROUND

In the U.S. 3.8 million people have a Traumatic Brain Injury (TBI) each year. Rapid brain training exercises to improve cognitive function after a mild TBI are needed.

OBJECTIVE

This study determines whether cognitive remediation by discriminating the direction a test pattern moves relative to a stationary background (movement figure-ground discrimination) improves the vision and cognitive deficits that result from a TBI, providing a paradigm shift in treatment methods.

METHODS

Movement-discrimination neurotraining was used to remediate low-level visual timing deficits in the dorsal stream to determine whether it improved high-level cognitive functions, such as processing speed, reading fluency, and the executive control functions of attention and working memory in four men with a TBI between the ages of 15-68. Standardized tests, as well as Magnetoencephalography (MEG) brain imaging, were administered at the beginning and end of 8-16 weeks of intervention training to evaluate improvements in cognitive skills.

RESULTS

Movement-discrimination cognitive neurotraining remediated both low-level visual timing deficits and high-level cognitive functioning, including selective and sustained attention, reading fluency, processing speed, and working memory for all TBI patients we studied. MEG brain imaging, using the Fast-VESTAL procedure, showed that this movement-discrimination training improved time-locked activity in the dorsal stream, attention, and executive control networks.

CONCLUSIONS

Remediating visual timing deficits in the dorsal stream revealed the causal role of visual movement discrimination training in improving high-level cognitive functions such as focusing and switching attention, working memory, processing speed, and reading. This study found that movement-discrimination training was very rapid and effective in remediating cognitive deficits, providing a new approach that is very beneficial for treating a mild TBI.

Bayesian Electromagnetic Spatio-Temporal Imaging of Extended Sources Based on Matrix Factorization

Accurate estimation of the locations and extents of neural sources from electroencephalography and magnetoencephalography (E/MEG) is challenging, especially for deep and highly correlated neural activities. In this study, we proposed a new fully data-driven source imaging method, source imaging based on spatio-temporal basis function (SI-STBF), which is built upon a Bayesian framework, to address this issue. The SI-STBF is based on the factorization of a source matrix as a product of a sparse coding matrix and a temporal basis function (TBF) matrix, which includes a few TBFs. The prior of the TBF is set in the empirical Bayesian manner. Similarly, for the spatial constraint, the SI-STBF assumes the prior covariance of the coding matrix as a weighted sum of several spatial covariance components. Both the TBFs and the coding matrix are learned from E/MEG simultaneously through variational Bayesian inference. To enable inference on high-resolution source space, we derived a scalable algorithm using convex analysis. The performance of the SI-STBF was assessed using both simulated and experimental E/MEG recordings. Compared with L₂-norm constrained methods, the SI-STBF is superior in reconstructing extended sources with less spatial diffusion and less localization error. By virtue of the spatio-temporal factorization of source matrix, the SI-STBF also produces more accurate estimations than spatial-only constraint method for high correlated and deep sources.

Modulating neural oscillations by transcranial static magnetic field stimulation of the dorsolateral prefrontal cortex: A crossover, double-blind, sham-controlled pilot study

Transcranial static magnetic field stimulation (tSMS) is a novel non-invasive brain stimulation technique that has been shown to locally increase alpha power in the parietal and occipital cortex. We investigated if tSMS locally increased alpha power in the left or right prefrontal cortex, as the balance of left/right prefrontal alpha power (frontal alpha asymmetry) has been linked to emotional processing and mood disorders. Therefore, altering frontal alpha asymmetry with tSMS may serve as a novel treatment to psychiatric diseases. We performed a crossover, double-blind, sham-controlled pilot study to assess the effects of prefrontal tSMS on neural oscillations. Twenty-four right-handed healthy participants were recruited and received left dorsolateral prefrontal cortex (DLPFC) tSMS, right DLPFC tSMS, and sham tSMS in a randomized order. Electroencephalography data were collected before (2 min eyes-closed, 2 min eyes-open), during (10 min eyes-open), and after (2 min eyes-open) stimulation. In contrast with our hypothesis, neither left nor right tSMS locally increased frontal alpha power. However, alpha power increased in occipital cortex during left DLPFC tSMS. Right DLPFC tSMS increased post-stimulation fronto-parietal theta power, indicating possible relevance to memory and cognition. Left and right DLPFC tSMS increased post-stimulation left hemisphere beta power, indicating possible changes to motor behavior. Left DLPFC tSMS also increased post-stimulation right frontal beta power, demonstrating complex network effects that may be relevant to aggressive behavior. We concluded that DLPFC tSMS modulated the network oscillations in regions distant from the location of stimulation and that tSMS has region specific effects on neural oscillations.

Targeted synchronization in an externally driven population of mechanical oscillators

We experimentally investigate the synchronization of driven metronomes using a servo motor to impose external control. We show that a driven metronome will only synchronize in a narrow range near its own frequency; when we introduce coupling between metronomes, we can widen the range of frequencies over which a metronome will synchronize to the external input. Using these features, we design a signal to synchronize a population of dissimilar metronomes; separately we design a signal to selectively synchronize a subpopulation of metronomes within a heterogeneous population.

Robust EEG/MEG Based Functional Connectivity with the Envelope of the Imaginary Coherence: Sensor Space Analysis

The brain's functional connectivity (FC) estimated at sensor level from electromagnetic (EEG/MEG) signals can provide quick and useful information towards understanding cognition and brain disorders. Volume conduction (VC) is a fundamental issue in FC analysis due to the effects of instantaneous correlations. FC methods based on the imaginary part of the coherence (iCOH) of any two signals are readily robust to VC effects, but neglecting the real part of the coherence leads to negligible FC when the processes are truly connected but with zero or π-phase (modulus 2π) interaction. We ameliorate this issue by proposing a novel method that implements an envelope of the imaginary coherence (EIC) to approximate the coherence estimate of supposedly active underlying sources. We compare EIC with state-of-the-art FC measures that included lagged coherence, iCOH, phase lag index (PLI) and weighted PLI (wPLI), using bivariate autoregressive and stochastic neural mass models. Additionally, we create realistic simulations where three and five regions were mapped on a template cortical surface and synthetic MEG signals were obtained after computing the electromagnetic leadfield. With this simulation and comparison study, we also demonstrate the feasibility of sensor FC analysis using receiver operating curve analysis whilst varying the signal's noise level. However, these results should be interpreted with caution given the known limitations of the sensor-based FC approach. Overall, we found that EIC and iCOH demonstrate superior results with most accurate FC maps. As they complement each other in different scenarios, that will be important to study normal and diseased brain activity.

Development of advanced signal processing and source imaging methods for superparamagnetic relaxometry

Superparamagnetic relaxometry (SPMR) is a highly sensitive technique for the in vivo detection of tumor cells and may improve early stage detection of cancers. SPMR employs superparamagnetic iron oxide nanoparticles (SPION). After a brief magnetizing pulse is used to align the SPION, SPMR measures the time decay of SPION using super-conducting quantum interference device (SQUID) sensors. Substantial research has been carried out in developing the SQUID hardware and in improving the properties of the SPION. However, little research has been done in the pre-processing of sensor signals and post-processing source modeling in SPMR. In the present study, we illustrate new pre-processing tools that were developed to: (1) remove trials contaminated with artifacts, (2) evaluate and ensure that a single decay process associated with bounded SPION exists in the data, (3) automatically detect and correct flux jumps, and (4) accurately fit the sensor signals with different decay models. Furthermore, we developed an automated approach based on multi-start dipole imaging technique to obtain the locations and magnitudes of multiple magnetic sources, without initial guesses from the users. A regularization process was implemented to solve the ambiguity issue related to the SPMR source variables. A procedure based on reduced chi-square cost-function was introduced to objectively obtain the adequate number of dipoles that describe the data. The new pre-processing tools and multi-start source imaging approach have been successfully evaluated using phantom data. In conclusion, these tools and multi-start source modeling approach substantially enhance the accuracy and sensitivity in detecting and localizing sources from the SPMR signals. Furthermore, multi-start approach with regularization provided robust and accurate solutions for a poor SNR condition similar to the SPMR detection sensitivity in the order of 1000 cells. We believe such algorithms will help establishing the industrial standards for SPMR when applying the technique in pre-clinical and clinical settings.

Resting-State Magnetoencephalography Reveals Different Patterns of Aberrant Functional Connectivity in Combat-Related Mild Traumatic Brain Injury

Blast mild traumatic brain injury (mTBI) is a leading cause of sustained impairment in military service members and veterans. However, the mechanism of persistent disability is not fully understood. The present study investigated disturbances in brain functioning in mTBI participants using a source-imaging-based approach to analyze functional connectivity (FC) from resting-state magnetoencephalography (rs-MEG). Study participants included 26 active-duty service members or veterans who had blast mTBI with persistent post-concussive symptoms, and 22 healthy control active-duty service members or veterans. The source time courses from regions of interest (ROIs) were used to compute ROI to whole-brain (ROI-global) FC for different frequency bands using two different measures: 1) time-lagged cross-correlation and 2) phase-lock synchrony. Compared with the controls, blast mTBI participants showed increased ROI-global FC in beta, gamma, and low-frequency bands, but not in the alpha band. Sources of abnormally increased FC included the: 1) prefrontal cortex (right ventromedial prefrontal cortex [vmPFC], right rostral anterior cingulate cortex [rACC]), and left ventrolateral and dorsolateral prefrontal cortex; 2) medial temporal lobe (bilateral parahippocampus, hippocampus, and amygdala); and 3) right putamen and cerebellum. In contrast, the blast mTBI group also showed decreased FC of the right frontal pole. Group differences were highly consistent across the two different FC measures. FC of the left ventrolateral prefrontal cortex correlated with executive functioning and processing speed in mTBI participants. Altogether, our findings of increased and decreased regionalpatterns of FC suggest that disturbances in intrinsic brain connectivity may be the result of multiple mechanisms, and are associated with cognitive sequelae of the injury.

The Iterative Reweighted Mixed-Norm Estimate for Spatio-Temporal MEG/EEG Source Reconstruction

Source imaging based on magnetoencephalography (MEG) and electroencephalography (EEG) allows for the non-invasive analysis of brain activity with high temporal and good spatial resolution. As the bioelectromagnetic inverse problem is ill-posed, constraints are required. For the analysis of evoked brain activity, spatial sparsity of the neuronal activation is a common assumption. It is often taken into account using convex constraints based on the l1-norm. The resulting source estimates are however biased in amplitude and often suboptimal in terms of source selection due to high correlations in the forward model. In this work, we demonstrate that an inverse solver based on a block-separable penalty with a Frobenius norm per block and a l0.5-quasinorm over blocks addresses both of these issues. For solving the resulting non-convex optimization problem, we propose the iterative reweighted Mixed Norm Estimate (irMxNE), an optimization scheme based on iterative reweighted convex surrogate optimization problems, which are solved efficiently using a block coordinate descent scheme and an active set strategy. We compare the proposed sparse imaging method to the dSPM and the RAP-MUSIC approach based on two MEG data sets. We provide empirical evidence based on simulations and analysis of MEG data that the proposed method improves on the standard Mixed Norm Estimate (MxNE) in terms of amplitude bias, support recovery, and stability.

Resting State Functional Connectivity MRI among Spectral MEG Current Sources in Children on the Autism Spectrum

Social and communicative impairments are among the core symptoms of autism spectrum disorders (ASD), and a great deal of evidence supports the notion that these impairments are associated with aberrant functioning and connectivity of various cortical networks. The present study explored the links between sources of MEG amplitude in various frequency bands and functional connectivity MRI in the resting state. The goal of combining these modalities was to use sources of neural oscillatory activity, measured with MEG, as functionally relevant seed regions for a more traditional pairwise fMRI connectivity analysis. We performed a seed-based connectivity analysis on resting state fMRI data, using seed regions derived from frequency-specific amplitude sources in resting state MEG data in the same nine subjects with ASD (10–17 years of age). We then compared fMRI connectivity among these MEG-source-derived regions between participants with autism and typically developing, age-matched controls. We used a source modeling technique designed for MEG data to detect significant amplitude sources in six frequency bands: delta (2–4 Hz), theta (4–8 Hz), alpha (8–12 Hz), beta (12–30 Hz), low gamma (30–60 Hz), and high gamma (60–120 Hz). MEG-derived source maps for each participant were co-registered in standard MNI space, and group-level source maps were obtained for each frequency. For each frequency band, the 10 largest clusters resulting from these t-tests were used as regions of interest (ROIs) for the fMRI functional connectivity analysis. Pairwise BOLD signal correlations were obtained between each pair of these ROIs for each frequency band. Each pairwise correlation was compared between the ASD and TD groups using t-tests. We also constrained these pairwise correlations to known network structures, resulting in a follow-up set of correlation matrices specific to each network we considered. Frequency-specific MEG sources had distinct patterns of fMRI resting state functional connectivity in the ASD group, but perhaps the most significant was a finding of hypoconnectivity between many sources of low and high gamma activity. These novel findings suggest that in ASD there are differences in functionally defined networks as shown in previous fMRI studies, as well as between sets of regions defined by magnetoencephalographic neural oscillatory activity.

Multi-Dimensional Dynamics of Human Electromagnetic Brain Activity

Magnetoencephalography (MEG) and electroencephalography (EEG) are invaluable neuroscientific tools for unveiling human neural dynamics in three dimensions (space, time, and frequency), which are associated with a wide variety of perceptions, cognition, and actions. MEG/EEG also provides different categories of neuronal indices including activity magnitude, connectivity, and network properties along the three dimensions. In the last 20 years, interest has increased in inter-regional connectivity and complex network properties assessed by various sophisticated scientific analyses. We herein review the definition, computation, short history, and pros and cons of connectivity and complex network (graph-theory) analyses applied to MEG/EEG signals. We briefly describe recent developments in source reconstruction algorithms essential for source-space connectivity and network analyses. Furthermore, we discuss a relatively novel approach used in MEG/EEG studies to examine the complex dynamics represented by human brain activity. The correct and effective use of these neuronal metrics provides a new insight into the multi-dimensional dynamics of the neural representations of various functions in the complex human brain.

The role of biomarkers and MEG-based imaging markers in the diagnosis of post-traumatic stress disorder and blast-induced mild traumatic brain injury

BACKGROUND

Pervasive use of improvised explosive devices (IEDs), rocket-propelled grenades, and land mines in the recent conflicts in Iraq and Afghanistan has brought traumatic brain injury (TBI) and its impact on health outcomes into public awareness. Blast injuries have been deemed signature wounds of these wars. War-related TBI is not new, having become prevalent during WWI and remaining medically relevant in WWII and beyond. Medicine's past attempts to accurately diagnose and disentangle the pathophysiology of war-related TBI parallels current lines of inquiry and highlights limitations in methodology and attribution of symptom etiology, be it organic, psychological, or behavioral. New approaches and biomarkers are needed.

PRECLINICAL

Serological biomarkers and biomarkers of injury obtained with imaging techniques represent cornerstones in the translation between experimental data and clinical observations. Experimental models for blast related TBI and PTSD can generate critical data on injury threshold, for example for white matter injury from acceleration. Carefully verified and validated models can be evaluated with gene expression arrays and proteomics to identify new candidates for serological biomarkers. Such models can also be analyzed with diffusion MRI and microscopy in order to identify criteria for detection of diffuse white matter injuries, such as DAI (diffuse axonal injury). The experimental models can also be analyzed with focus on injury outcome in brain stem regions, such as locus coeruleus or nucleus raphe magnus that can be involved in response to anxiety changes.

CLINICAL

Mild (and some moderate) TBI can be difficult to diagnose because the injuries are often not detectable on conventional MRI or CT. There is accumulating evidence that injured brain tissues in TBI patients generate abnormal low-frequency magnetic activity (ALFMA, peaked at 1-4Hz) that can be measured and localized by magnetoencephalography (MEG). MEG imaging detects TBI abnormalities at the rates of 87% for the mild TBI, group (blast-induced plus non-blast causes) and 100% for the moderate group. Among the mild TBI patients, the rates of abnormalities are 96% and 77% for the blast and non-blast TBI groups, respectively. There is emerging evidence based on fMRI and MEG studies showing hyper-activity in the amygdala and hypo-activity in pre-frontal cortex in individuals with PTSD. MEG signal may serve as a sensitive imaging marker for mTBI, distinguishable from abnormalities generated in association with PTSD. More work is needed to fully describe physiological mechanisms of post-concussive symptoms.

On characterizing population commonalities and subject variations in brain networks

Brain networks based on resting state connectivity as well as inter-regional anatomical pathways obtained using diffusion imaging have provided insight into pathology and development. Such work has underscored the need for methods that can extract sub-networks that can accurately capture the connectivity patterns of the underlying population while simultaneously describing the variation of sub-networks at the subject level. We have designed a multi-layer graph clustering method that extracts clusters of nodes, called 'network hubs', which display higher levels of connectivity within the cluster than to the rest of the brain. The method determines an atlas of network hubs that describes the population, as well as weights that characterize subject-wise variation in terms of within- and between-hub connectivity. This lowers the dimensionality of brain networks, thereby providing a representation amenable to statistical analyses. The applicability of the proposed technique is demonstrated by extracting an atlas of network hubs for a population of typically developing controls (TDCs) as well as children with autism spectrum disorder (ASD), and using the structural and functional networks of a population to determine the subject-level variation of these hubs and their inter-connectivity. These hubs are then used to compare ASD and TDCs. Our method is generalizable to any population whose connectivity (structural or functional) can be captured via non-negative network graphs.

A geometric correction scheme for spatial leakage effects in MEG/EEG seed-based functional connectivity mapping

Spatial leakage effects are particularly confounding for seed-based investigations of brain networks using source-level electroencephalography (EEG) or magnetoencephalography (MEG). Various methods designed to avoid this issue have been introduced but are limited to particular assumptions about its temporal characteristics. Here, we investigate the usefulness of a model-based geometric correction scheme (GCS) to suppress spatial leakage emanating from the seed location. We analyze its properties theoretically and then assess potential advantages and limitations with simulated and experimental MEG data (resting state and auditory-motor task). To do so, we apply Minimum Norm Estimation (MNE) for source reconstruction and use variation of error parameters, statistical gauging of spatial leakage correction and comparison with signal orthogonalization. Results show that the GCS has a local (i.e., near the seed) effect only, in line with the geometry of MNE spatial leakage, and is able to map spatially all types of brain interactions, including linear correlations eliminated after signal orthogonalization. Furthermore, it is robust against the introduction of forward model errors. On the other hand, the GCS can be affected by local overcorrection effects and seed mislocation. These issues arise with signal orthogonalization too, although significantly less extensively, so the two approaches complement each other. The GCS thus appears to be a valuable addition to the spatial leakage correction toolkits for seed-based FC analyses in source-projected MEG/EEG data.

Magnetoencephalography Slow-Wave Detection in Patients with Mild Traumatic Brain Injury and Ongoing Symptoms Correlated with Long-Term Neuropsychological Outcome

Mild traumatic brain injury (mTBI) is common in the United States, accounting for as many as 75-80% of all TBIs. It is recognized as a significant public health concern, but there are ongoing controversies regarding the etiology of persistent symptoms post-mTBI. This constellation of nonspecific symptoms is referred to as postconcussive syndrome (PCS). The present study combined results from magnetoencephalography (MEG) and cognitive assessment to examine group differences and relationships between brain activity and cognitive performance in 31 military and civilian individuals with a history of mTBI+PCS and 33 matched healthy control subjects. An operator-free analysis was used for MEG data to increase reliability of the technique. Subjects completed a comprehensive neuropsychological assessment, and measures of abnormal slow-wave activity from MEG were collected. Results demonstrated significant group differences on measures of executive functioning and processing speed. In addition, significant correlations between slow-wave activity on MEG and patterns of cognitive functioning were found in cortical areas, consistent with cognitive impairments on exams. Results provide more objective evidence that there may be subtle changes to the neurobiological integrity of the brain that can be detected by MEG. Further, these findings suggest that these abnormalities are associated with cognitive outcomes and may account, at least in part, for long-term PCS in those who have sustained an mTBI.

What's New in Traumatic Brain Injury: Update on Tracking, Monitoring and Treatment

Traumatic brain injury (TBI), defined as an alteration in brain functions caused by an external force, is responsible for high morbidity and mortality around the world. It is important to identify and treat TBI victims as early as possible. Tracking and monitoring TBI with neuroimaging technologies, including functional magnetic resonance imaging (fMRI), diffusion tensor imaging (DTI), positron emission tomography (PET), and high definition fiber tracking (HDFT) show increasing sensitivity and specificity. Classical electrophysiological monitoring, together with newly established brain-on-chip, cerebral microdialysis techniques, both benefit TBI. First generation molecular biomarkers, based on genomic and proteomic changes following TBI, have proven effective and economical. It is conceivable that TBI-specific biomarkers will be developed with the combination of systems biology and bioinformation strategies. Advances in treatment of TBI include stem cell-based and nanotechnology-based therapy, physical and pharmaceutical interventions and also new use in TBI for approved drugs which all present favorable promise in preventing and reversing TBI.

Filling in the gaps: Anticipatory control of eye movements in chronic mild traumatic brain injury

A barrier in the diagnosis of mild traumatic brain injury (mTBI) stems from the lack of measures that are adequately sensitive in detecting mild head injuries. MRI and CT are typically negative in mTBI patients with persistent symptoms of post-concussive syndrome (PCS), and characteristic difficulties in sustaining attention often go undetected on neuropsychological testing, which can be insensitive to momentary lapses in concentration. Conversely, visual tracking strongly depends on sustained attention over time and is impaired in chronic mTBI patients, especially when tracking an occluded target. This finding suggests deficient internal anticipatory control in mTBI, the neural underpinnings of which are poorly understood. The present study investigated the neuronal bases for deficient anticipatory control during visual tracking in 25 chronic mTBI patients with persistent PCS symptoms and 25 healthy control subjects. The task was performed while undergoing magnetoencephalography (MEG), which allowed us to examine whether neural dysfunction associated with anticipatory control deficits was due to altered alpha, beta, and/or gamma activity. Neuropsychological examinations characterized cognition in both groups. During MEG recordings, subjects tracked a predictably moving target that was either continuously visible or randomly occluded (gap condition). MEG source-imaging analyses tested for group differences in alpha, beta, and gamma frequency bands. The results showed executive functioning, information processing speed, and verbal memory deficits in the mTBI group. Visual tracking was impaired in the mTBI group only in the gap condition. Patients showed greater error than controls before and during target occlusion, and were slower to resynchronize with the target when it reappeared. Impaired tracking concurred with abnormal beta activity, which was suppressed in the parietal cortex, especially the right hemisphere, and enhanced in left caudate and frontal–temporal areas. Regional beta-amplitude demonstrated high classification accuracy (92%) compared to eye-tracking (65%) and neuropsychological variables (80%). These findings show that deficient internal anticipatory control in mTBI is associated with altered beta activity, which is remarkably sensitive given the heterogeneity of injuries.

•

Neuropsychological test performance was impaired in mTBI patients.

•

Visual tracking was impaired in the gap task, where targets were randomly occluded.

•

Impaired visual tracking concurred with abnormal MEG beta activity.

•

Beta was suppressed in parietal and enhanced in caudate and frontal–temporal areas.

•

Regional MEG beta-amplitude demonstrated high classification accuracy (92%).

Resting-state alpha in autism spectrum disorder and alpha associations with thalamic volume

Alpha circuits (8-12 Hz), necessary for basic and complex brain processes, are abnormal in autism spectrum disorder (ASD). The present study obtained estimates of resting-state (RS) alpha activity in children with ASD and examined associations between alpha activity, age, and clinical symptoms. Given that the thalamus modulates cortical RS alpha rhythms, associations between thalamic structure and alpha activity were examined. RS magnetoencephalography was obtained from 47 typically-developing children (TDC) and 41 children with ASD. RS alpha activity was measured using distributed source localization. Left and right thalamic volume measurements were also obtained. In both groups, the strongest alpha activity was observed in Calcarine Sulcus regions. In Calcarine regions, only TDC showed the expected association between age and alpha peak frequency. ASD had more alpha activity than TDC in regions bordering the Central Sulcus as well as parietal association cortices. In ASD, whereas greater left Central Sulcus relative alpha activity was associated with higher Social Responsiveness Scale (SRS) scores, greater Calcarine region relative alpha activity was associated with lower SRS scores. Although thalamic volume group differences were not observed, relationships between thalamic volume and Calcarine alpha power were unique to TDC. The present study also identified a failure to shift peak alpha frequency as a function of age in primary alpha-generating areas in children with ASD. Findings suggested that increased RS alpha activity in primary motor and somatosensory as well as parietal multimodal areas-with increased alpha thought to reflect greater inhibition-might impair the ability to identify or interpret social cues. Finally, to our knowledge, this is the first study to report associations between thalamic volume and alpha power, an association observed only in TDC. The lack of thalamic and alpha associations in ASD suggests thalamic contributions to RS alpha abnormalities in ASD.

Alpha-to-gamma phase-amplitude coupling methods and application to autism spectrum disorder

Adult studies have shown that a basic property of resting-state (RS) brain activity is the coupling of posterior alpha oscillations (alpha phase) to posterior gamma oscillations (gamma amplitude). The present study examined whether this basic RS process is present in children. Given reports of abnormal parietal-occipital RS alpha in children with autism spectrum disorder (ASD), the present study examined whether RS alpha-to-gamma phase-amplitude coupling (PAC) is disrupted in ASD. Simulations presented in this study showed limitations with traditional PAC analyses. In particular, to avoid false-positive PAC findings, simulations showed the need to use a unilateral passband to filter the upper frequency band as well as the need for longer epochs of data. For the human study, eyes-closed RS magnetoencephalography data were analyzed from 25 children with ASD and 18 typically developing (TD) children with at least 60 sec of artifact-free data. Source modeling provided continuous time course data at a midline parietal-occipital source for PAC analyses. Greater alpha-to-gamma PAC was observed in ASD than TD (p<0.005). Although children with ASD had higher PAC values, in both groups gamma activity increased at the peak of the alpha oscillation. In addition, an association between alpha power and alpha-to-gamma PAC was observed in both groups, although this relationship was stronger in ASD than TD (p<0.05). Present results demonstrated that although alpha-to-gamma PAC is present in children, this basic RS process is abnormal in children with ASD. Finally, simulations and the human data highlighted the need to consider the interplay between alpha power, epoch length, and choice of signal processing methods on PAC estimates.

Voxel-wise resting-state MEG source magnitude imaging study reveals neurocircuitry abnormality in active-duty service members and veterans with PTSD

Post-traumatic stress disorder (PTSD) is a leading cause of sustained impairment, distress, and poor quality of life in military personnel, veterans, and civilians. Indirect functional neuroimaging studies using PET or fMRI with fear-related stimuli support a PTSD neurocircuitry model that includes amygdala, hippocampus, and ventromedial prefrontal cortex (vmPFC). However, it is not clear if this model can fully account for PTSD abnormalities detected directly by electromagnetic-based source imaging techniques in resting-state. The present study examined resting-state magnetoencephalography (MEG) signals in 25 active-duty service members and veterans with PTSD and 30 healthy volunteers. In contrast to the healthy volunteers, individuals with PTSD showed: 1) hyperactivity from amygdala, hippocampus, posterolateral orbitofrontal cortex (OFC), dorsomedial prefrontal cortex (dmPFC), and insular cortex in high-frequency (i.e., beta, gamma, and high-gamma) bands; 2) hypoactivity from vmPFC, Frontal Pole (FP), and dorsolateral prefrontal cortex (dlPFC) in high-frequency bands; 3) extensive hypoactivity from dlPFC, FP, anterior temporal lobes, precuneous cortex, and sensorimotor cortex in alpha and low-frequency bands; and 4) in individuals with PTSD, MEG activity in the left amygdala and posterolateral OFC correlated positively with PTSD symptom scores, whereas MEG activity in vmPFC and precuneous correlated negatively with symptom score. The present study showed that MEG source imaging technique revealed new abnormalities in the resting-state electromagnetic signals from the PTSD neurocircuitry. Particularly, posterolateral OFC and precuneous may play important roles in the PTSD neurocircuitry model.

•

Resting-state MEG detects abnormal electromagnetic activity in PTSD neurocircuitry

•

PTSD showed hyperactivity in amygdala, hippocampus, and orbitofrontal cortex

•

PTSD showed hypoactivity in vmPFC, frontal pole, and dlPFC

•

PTSD symptom score correlated with MEG activity

Single-subject-based whole-brain MEG slow-wave imaging approach for detecting abnormality in patients with mild traumatic brain injury

Traumatic brain injury (TBI) is a leading cause of sustained impairment in military and civilian populations. However, mild TBI (mTBI) can be difficult to detect using conventional MRI or CT. Injured brain tissues in mTBI patients generate abnormal slow-waves (1–4 Hz) that can be measured and localized by resting-state magnetoencephalography (MEG). In this study, we develop a voxel-based whole-brain MEG slow-wave imaging approach for detecting abnormality in patients with mTBI on a single-subject basis. A normative database of resting-state MEG source magnitude images (1–4 Hz) from 79 healthy control subjects was established for all brain voxels. The high-resolution MEG source magnitude images were obtained by our recent Fast-VESTAL method. In 84 mTBI patients with persistent post-concussive symptoms (36 from blasts, and 48 from non-blast causes), our method detected abnormalities at the positive detection rates of 84.5%, 86.1%, and 83.3% for the combined (blast-induced plus with non-blast causes), blast, and non-blast mTBI groups, respectively. We found that prefrontal, posterior parietal, inferior temporal, hippocampus, and cerebella areas were particularly vulnerable to head trauma. The result also showed that MEG slow-wave generation in prefrontal areas positively correlated with personality change, trouble concentrating, affective lability, and depression symptoms. Discussion is provided regarding the neuronal mechanisms of MEG slow-wave generation due to deafferentation caused by axonal injury and/or blockages/limitations of cholinergic transmission in TBI. This study provides an effective way for using MEG slow-wave source imaging to localize affected areas and supports MEG as a tool for assisting the diagnosis of mTBI.

•

A voxel-based whole-brain MEG slow-wave source imaging method for mild TBI.

•

The new approach showed 84.5% positive detection rate in 84 mild TBI patients.

•

The new approach detected loci of injury in mild TBI patients on a single-subject basis.

•

MEG slow-wave source imaging revealed brain areas vulnerable to mild TBI.

•

MEG slow-wave generations correlated with mild TBI symptoms.