Electro/Magneto-EncephaloGraphy (EEG/MEG) source imaging (EMSI) of epileptic activity from deep generators is often challenging due to the higher sensitivity of EEG/MEG to superficial regions and to the spatial configuration of subcortical structures. We previously demonstrated the ability of the coherent Maximum Entropy on the Mean (cMEM) method to accurately localize the superficial cortical generators and their spatial extent. Here, we propose a depth-weighted adaptation of cMEM to localize deep generators more accurately. These methods were evaluated using realistic MEG/high-density EEG (HD-EEG) simulations of epileptic activity and actual MEG/HD-EEG recordings from patients with focal epilepsy. We incorporated depth-weighting within the MEM framework to compensate for its preference for superficial generators. We also included a mesh of both hippocampi, as an additional deep structure in the source model. We generated 5400 realistic simulations of interictal epileptic discharges for MEG and HD-EEG involving a wide range of spatial extents and signal-to-noise ratio (SNR) levels, before investigating EMSI on clinical HD-EEG in 16 patients and MEG in 14 patients. Clinical interictal epileptic discharges were marked by visual inspection. We applied three EMSI methods: cMEM, depth-weighted cMEM and depth-weighted minimum norm estimate (MNE). The ground truth was defined as the true simulated generator or as a drawn region based on clinical information available for patients. For deep sources, depth-weighted cMEM improved the localization when compared to cMEM and depth-weighted MNE, whereas depth-weighted cMEM did not deteriorate localization accuracy for superficial regions. For patients\' data, we observed improvement in localization for deep sources, especially for the patients with mesial temporal epilepsy, for which cMEM failed to reconstruct the initial generator in the hippocampus. Depth weighting was more crucial for MEG (gradiometers) than for HD-EEG. Similar findings were found when considering depth weighting for the wavelet extension of MEM. In conclusion, depth-weighted cMEM improved the localization of deep sources without or with minimal deterioration of the localization of the superficial sources. This was demonstrated using extensive simulations with MEG and HD-EEG and clinical MEG and HD-EEG for epilepsy patients.

This study assesses the reliability of resting-state dynamic causal modelling (DCM) of magnetoencephalography (MEG) under conductance-based canonical microcircuit models, in terms of both posterior parameter estimates and model evidence. We use resting-state MEG data from two sessions, acquired 2 weeks apart, from a cohort with high between-subject variance arising from Alzheimer\'s disease. Our focus is not on the effect of disease, but on the reliability of the methods (as within-subject between-session agreement), which is crucial for future studies of disease progression and drug intervention. To assess the reliability of first-level DCMs, we compare model evidence associated with the covariance among subject-specific free energies (i.e., the \'quality\' of the models) with versus without interclass correlations. We then used parametric empirical Bayes (PEB) to investigate the differences between the inferred DCM parameter probability distributions at the between subject level. Specifically, we examined the evidence for or against parameter differences (i) within-subject, within-session, and between-epochs; (ii) within-subject between-session; and (iii) within-site between-subjects, accommodating the conditional dependency among parameter estimates. We show that for data acquired close in time, and under similar circumstances, more than 95% of inferred DCM parameters are unlikely to differ, speaking to mutual predictability over sessions. Using PEB, we show a reciprocal relationship between a conventional definition of \'reliability\' and the conditional dependency among inferred model parameters. Our analyses confirm the reliability and reproducibility of the conductance-based DCMs for resting-state neurophysiological data. In this respect, the implicit generative modelling is suitable for interventional and longitudinal studies of neurological and psychiatric disorders.

Visual entrainment is a powerful and widely used research tool to study visual information processing in the brain. While many entrainment studies have focused on frequencies around 14-16 Hz, there is renewed interest in understanding visual entrainment at higher frequencies (e.g., gamma-band entrainment). Notably, recent groundbreaking studies have demonstrated that gamma-band visual entrainment at 40 Hz may have therapeutic effects in the context of Alzheimer\'s disease (AD) by stimulating specific neural ensembles, which utilize GABAergic signaling. Despite such promising findings, few studies have investigated the optimal parameters for gamma-band visual entrainment. Herein, we examined whether visual stimulation at 32, 40, or 48 Hz produces optimal visual entrainment responses using high-density magnetoencephalography (MEG). Our results indicated strong entrainment responses localizing to the primary visual cortex in each condition. Entrainment responses were stronger for 32 and 40 Hz relative to 48 Hz, indicating more robust synchronization of neural ensembles at these lower gamma-band frequencies. In addition, 32 and 40 Hz entrainment responses showed typical patterns of habituation across trials, but this effect was absent for 48 Hz. Finally, connectivity between visual cortex and parietal and prefrontal cortices tended to be strongest for 40 relative to 32 and 48 Hz entrainment. These results suggest that neural ensembles in the visual cortex may resonate at around 32 and 40 Hz and thus entrain more readily to photic stimulation at these frequencies. Emerging AD therapies, which have focused on 40 Hz entrainment to date, may be more effective at lower relative to higher gamma frequencies, although additional work in clinical populations is needed to confirm these findings. PRACTITIONER POINTS: Gamma-band visual entrainment has emerged as a therapeutic approach for eliminating amyloid in Alzheimer\'s disease, but its optimal parameters are unknown. We found stronger entrainment at 32 and 40 Hz compared to 48 Hz, suggesting neural ensembles prefer to resonate around these relatively lower gamma-band frequencies. These findings may inform the development and refinement of innovative AD therapies and the study of GABAergic visual cortical functions.

Humans can read and comprehend text rapidly, implying that readers might process multiple words per fixation. However, the extent to which parafoveal words are previewed and integrated into the evolving sentence context remains disputed. We investigated parafoveal processing during natural reading by recording brain activity and eye movements using MEG and an eye tracker while participants silently read one-line sentences. The sentences contained an unpredictable target word that was either congruent or incongruent with the sentence context. To measure parafoveal processing, we flickered the target words at 60 Hz and measured the resulting brain responses (i.e. Rapid Invisible Frequency Tagging, RIFT) during fixations on the pre-target words. Our results revealed a significantly weaker tagging response for target words that were incongruent with the previous context compared to congruent ones, even within 100ms of fixating the word immediately preceding the target. This reduction in the RIFT response was also found to be predictive of individual reading speed. We conclude that semantic information is not only extracted from the parafovea but can also be integrated with the previous context before the word is fixated. This early and extensive parafoveal processing supports the rapid word processing required for natural reading. Our study suggests that theoretical frameworks of natural reading should incorporate the concept of deep parafoveal processing.

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.
Individualized computational models of epilepsy surgery capture some of the key aspects of seizure propagation and the resective surgery. It is to be established whether this information can be integrated during the presurgical evaluation of the patient to improve surgical planning and the chances of a good surgical outcome. Here we address this question with a pseudo-prospective study that applies a computational framework of seizure propagation and epilepsy surgery—the ESSES framework—in a pseudo-prospective study mimicking the presurgical conditions. We found that within this pseudo-prospective setting, ESSES could correctly predict 75% of NSF and 80.8% of SF cases. This finding suggests the potential of individualised computational models to inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection.

Testosterone levels sharply rise during the transition from childhood to adolescence and these changes are known to be associated with changes in human brain structure. During this same developmental window, there are also robust changes in the neural oscillatory dynamics serving verbal working memory processing. Surprisingly, whereas many studies have investigated the effects of chronological age on the neural oscillations supporting verbal working memory, none have probed the impact of endogenous testosterone levels during this developmental period. Using a sample of 89 youth aged 6-14 years-old, we collected salivary testosterone samples and recorded magnetoencephalography during a modified Sternberg verbal working memory task. Significant oscillatory responses were identified and imaged using a beamforming approach and the resulting maps were subjected to whole-brain ANCOVAs examining the effects of testosterone and sex, controlling for age, during verbal working memory encoding and maintenance. Our primary results indicated robust testosterone-related effects in theta (4-7 Hz) and alpha (8-14 Hz) oscillatory activity, controlling for age. During encoding, females exhibited weaker theta oscillations than males in right cerebellar cortices and stronger alpha oscillations in left temporal cortices. During maintenance, youth with greater testosterone exhibited weaker alpha oscillations in right parahippocampal and cerebellar cortices, as well as regions across the left-lateralized language network. These results extend the existing literature on the development of verbal working memory processing by showing region and sex-specific effects of testosterone, and are the first results to link endogenous testosterone levels to the neural oscillatory activity serving verbal working memory, above and beyond the effects of chronological age.

Absence seizures are classically associated with behavioral arrest and transient deficits in consciousness, yet substantial variability exists in the severity of the impairment. Despite several decades of research on the topic, the pathophysiology of absence seizures and the mechanisms underlying behavioral impairment remain unclear. Several rationales have been proposed including widespread cortical deactivation, reduced perception of external stimuli, and transient suspension of the default mode network, among others. This review aims to summarize the current knowledge on the neural correlates of impaired consciousness in absence seizures. We review evidence from studies using animal models of absence epilepsy, electroencephalography, functional magnetic resonance imaging, magnetoencephalography, positron emission tomography, and single photon emission computed tomography.

Schizophrenia is a severe disruption in cognition and emotion, affecting fundamental human functions. In this study, we applied Multi-Scale Entropy analysis to resting-state Magnetoencephalography data from 54 schizophrenia patients and 98 healthy controls. This method quantifies the temporal complexity of the signal across different time scales using the concept of sample entropy. Results show significantly higher sample entropy in schizophrenia patients, primarily in central, parietal, and occipital lobes, peaking at time scales equivalent to frequencies between 15 and 24 Hz. To disentangle the contributions of the amplitude and phase components, we applied the same analysis to a phase-shuffled surrogate signal. The analysis revealed that most differences originate from the amplitude component in the δ, α, and β power bands. While the phase component had a smaller magnitude, closer examination reveals clear spatial patterns and significant differences across specific brain regions. We assessed the potential of multi-scale entropy as a schizophrenia biomarker by comparing its classification performance to conventional spectral analysis and a cognitive task (the n-back paradigm). The discriminative power of multi-scale entropy and spectral features was similar, with a slight advantage for multi-scale entropy features. The results of the n-back test were slightly below those obtained from multi-scale entropy and spectral features.

Amyotrophic lateral sclerosis (ALS) is an idiopathic, fatal, and fast-progressive neurodegenerative disease characterized by the degeneration of motor neurons. ALS patients often experience an initial misdiagnosis or a diagnostic delay due to the current unavailability of an efficient biomarker. Since impaired speech is typical in ALS, we hypothesized that functional differences between healthy and ALS participants during speech tasks can be explained by cortical pattern changes, thereby leading to the identification of a neural biomarker for ALS. In this pilot study, we collected magnetoencephalography (MEG) recordings from three early-diagnosed patients with ALS and three healthy controls during imagined (covert) and overt speech tasks. First, we computed sensor correlations, which showed greater correlations for speakers with ALS than healthy controls. Second, we compared the power of the MEG signals in canonical bands between the two groups, which showed greater dissimilarity in the beta band for ALS participants. Third, we assessed differences in functional connectivity, which showed greater beta band connectivity for ALS than healthy controls. Finally, we performed single-trial classification, which resulted in highest performance with beta band features (∼ 98%). These findings were consistent across trials, phrases, and participants for both imagined and overt speech tasks. Our preliminary results indicate that speech-evoked beta oscillations could be a potential neural biomarker for diagnosing ALS. To our knowledge, this is the first demonstration of the detection of ALS from single-trial neural signals.

Human language relies on the correct processing of syntactic information, as it is essential for successful communication between speakers. As an abstract level of language, syntax has often been studied separately from the physical form of the speech signal, thus often masking the interactions that can promote better syntactic processing in the human brain. However, behavioral and neural evidence from adults suggests the idea that prosody and syntax interact, and studies in infants support the notion that prosody assists language learning. Here we analyze a MEG dataset to investigate how acoustic cues, specifically prosody, interact with syntactic representations in the brains of native English speakers. More specifically, to examine whether prosody enhances the cortical encoding of syntactic representations, we decode syntactic phrase boundaries directly from brain activity, and evaluate possible modulations of this decoding by the prosodic boundaries. Our findings demonstrate that the presence of prosodic boundaries improves the neural representation of phrase boundaries, indicating the facilitative role of prosodic cues in processing abstract linguistic features. This work has implications for interactive models of how the brain processes different linguistic features. Future research is needed to establish the neural underpinnings of prosody-syntax interactions in languages with different typological characteristics.