Name: Mapping Cortical Dynamics Using Simultaneous MEG/EEG and Anatomically-constrained Minimum-norm Estimates: an Auditory Attention Example (Video) | JoVE Journal
Uploaded: 2012-10-24T13:00:00.000+00:00
Duration: 8 min 45 s
Description: Learn to map cortical dynamics during auditory attention using simultaneous MEG/EEG and anatomically-constrained minimum-norm estimates. This method enables researchers to distinguish cognitive states through high-resolution brain activity patterns.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

The authors would like to thank Matti S. H&auml;m&auml;l&auml;inen, Lilla Z&ouml;llei and three anonymous reviewers for their helpful comments. Funding sources: R00DC010196 (AKCL); T32DC000018 (EDL); T32DC005361 (RKM).

1. Anatomical Data Acquisition and Processing
<ol>
 <li>Acquire one magnetization-prepared rapid gradient echo (MPRAGE) MRI scan of the subject. This may take 5-10 min depending on which specific scanning protocol is used.</li>
 <li>Acquire two additional fast low-angle shot (FLASH) MRI scans (flip angles = 5° and 30°) if EEG data are used for inverse imaging analysis, as FLASH sequences provide different tissue contrast from the standard MPRAGE sequences 1.</li>
 <li>Use FreeSurfer software (see Table) 2, 3 to reconstruct the cortical surface and to set up individual M/EEG dipole source space.
 <ol>
 <li>This source space is constrained to the grey/white matter boundary segmented from the MPRAGE scan. Each hemisphere contains roughly 100,000 potential vertices, spaced at ~1 mm. For dipole amplitude estimation (see below), use a grid spacing of 7 mm, which yields ~3,000 dipoles per hemisphere.</li>
 </ol>
 </li>
 <li>Reconstruct the skin, outer skull and inner skull surfaces from the MPRAGE and FLASH images using MNE (see Table) and FreeSurfer. Use these surfaces to generate a three-layer boundary element model (BEM).</li>
</ol>
2. M/EEG Data Acquisition
<ol>
 <li>Prepare subject for M/EEG recording.
 <ol>
 <li>Refer to Liu et al 4 for details of electrooculogram and reference electrode preparation as well as digitization of the subject's fiducial landmarks, head-position indicator coils (HPI) and EEG electrodes.</li>
 </ol>
 </li>
 <li>Once subject is seated in the MEG, measure the head position using the HPI coils.</li>
 <li>Start recording. Begin presentation of auditory and visual stimuli.
 <ol>
 <li>Many hardware and software solutions are available to perform stimulus presentation (e.g., Presentation, E-Prime). We use a Tucker-Davis Technologies RZ6 for auditory stimulus presentation and trigger stamping, with Psychtoolbox 5 for visual stimulus presentation, both controlled by MATLAB. Testing the auditory and visual latencies using a microphone and photodiode attached to the screen, and subsequently ensuring there is no observable jitter (which may necessitate setting the presentation projector to its native resolution) prior to the experiment helps ensure timing integrity.</li>
 </ol>
 </li>
 <li>Subject responds to auditory and visual stimuli via an optical button box while performing behavioral task.</li>
 <li>Save all stimuli, experimental parameters and data files for off-line analysis.</li>
</ol>
3. M/EEG Co-registration with MRI Scan and Data Processing
<ol>
 <li>Using the MNE software, load digitizer data and subject's reconstructed MRI head model. Pick fiducial landmarks to initiate co-registration process and proceed to use automatic alignment procedure to complete coordinate transformation (Figure 2).</li>
 <li>To relate the location of each dipole in the source space with the location of each sensor, combine recorded HPI data (see 2.2) to compute a forward solution with the three-layer BEM (see 1.4)</li>
 <li>Inspect all recorded M/EEG data and identify channels that have exceptionally high variance or are completely flat. Set these channels as bad channels.</li>
 <li>Use signal-space projection 6 or other noise reduction techniques (such as signal space separation 7) to project or separate out spatial field patterns originated from ambient environmental field contamination or other undesirable physiological signals, such as those associated with eye-blinks and cardiac artifacts (Figure 3).
 <ol>
 <li>Apply time-domain artifact removal (e.g., removing epochs containing abnormally high amplitude signals due to spiking of a channel) and frequency-domain artifact removal (e.g., band-notch filtering at 50 or 60 Hz line-frequency) to further increase signal-to-noise ratio.</li>
 </ol>
 </li>
 <li>Identify a baseline period in which the subject was not performing any task (e.g., 200 ms period prior to the onset of each trial). Generate an average of these baseline epochs in order to obtain a noise estimate (also known as the covariance matrix).</li>
 <li>Identify epochs of interest (e.g., only collecting epochs with correct behavioral responses) and define conditions for experimental contrasts (e.g., epochs associated with subject having switched their auditory attention to the opposite hemifield as originally cued - "Switch" condition - versus subject maintaining attention on the original hemifield - "Hold" condition). Generate an average response for each of the condition defined.
 <ol>
 <li>These averages can be baseline-corrected or not depending on experimental parameters (see 8); data shown here are baseline-corrected.</li>
 </ol>
 </li>
 <li>Combine the covariance matrix (3.5) and the computed forward solution (3.2) to obtain a distributed cortically constrained minimum-norm inverse operator that relates the sensor measurements to dipole current estimates in the source space.
 <ol>
 <li>You can either approximately constrain or fix the dipole orientation to the cortical normal direction 9.</li>
 </ol>
 </li>
 <li>Generate a "brain movie" of the distributed dipole estimate (i.e., the current estimate at each dipole location in the source space in time) for each experimental condition (Figure 4).
 <ol>
 <li>Depending on the temporal characteristics of your experimental design, you can bin your data in time by averaging current estimates using non-overlapping temporal windows.</li>
 </ol>
 </li>
</ol>
4. Statistical Inference Based on a Common Surface-based Coordinate System
<ol>
 <li>Morph the "brain movies" for each subject onto a common (average) cortical space based on a surface-based coordinate system that optimally aligns individual sulcal-gyral patterns 3. This allows us to compare or average cortical activities across all subjects. (Figure 5).</li>
 <li>There are many different statistical inference approaches. We will highlight three possible approaches here. Approaches that are not implemented in software package can be written using custom software-in our examples we use MATLAB to perform the non-parametric spatio-temporal clustering permutation test. Despite the high dimensionality (Space x Time x Subjects) of these data, all of these approaches can be performed using standard modern desktop computer hardware in seconds (ROI; 4.3 approach) to hours (non-parametric clustering; 4.5).</li>
 <li>Region-of-interest (ROI) approach
 <ol>
 <li>You can define the ROI anatomically (e.g., defined by automatic parcellation algorithm 1) and / or functionally (e.g., by recording a functional localizing task, such as a Go/No Go saccade task to identify the oculomotor regions).</li>
 <li>You can further constrain your analysis to a specific time-of-interest that is appropriate to your experimental paradigm (e.g., a time period immediately before and after the onset of the sound stimuli). You can also use other statistical inference associated with time-series analysis.</li>
 </ol>
 </li>
 <li>Whole-brain Bonferroni or False-Discovery-Rate (FDR) correction
 <ol>
 <li>Employ Bonferroni or FDR correction if you require whole-brain, all-time analysis.</li>
 <li>Generate a statistical map at each dipole location and each time point using appropriate test statistics, such as a t-test or within-subjects ANOVA for approximately normally-distributed data. For example, z-scores from dynamic statistical parameter mapping of the MNE estimates for fixed dipole sources 10 can be used when paired with a correction for correlations in the estimates (such as the conservative Greenhouse-Geisser correction).</li>
 <li>For Bonferroni correction, obtain significant space-time points by thresholding at a significance level of 0.05 divided by the number of comparisons (number of dipoles multiplied by number of time points). For a less conservative approach, use FDR p-value correction 11.</li>
 </ol>
 </li>
 <li>Non-parametric spatiotemporal clustering 
 <ol>
 <li>Use this method (based on a simple extension of 12) to find regions of large, consistent spatial and temporal activation while being less conservative than Bonferroni correction, and less prone to Type I statistical errors than FDR by controlling for the family-wise error rate.
 <ol>
 <li>Because this approach uses permutation or Monte Carlo resampling techniques, it does not rely on assumptions of normality of the data, and only assumes that the condition labels are exchangeable under the null hypothesis. Although it is more computationally intensive than the previous two approaches, it can still be performed in hours on a single machine using modern desktop computer hardware.</li>
 </ol>
 </li>
 <li>Generate a statistical map at each dipole location and each time point using the appropriate test statistics, such as a t-test.</li>
 <li>Threshold this map at a preliminary significance threshold, e.g. p &lt; 0.05.</li>
 <li>Cluster these putative significant points based on spatio-temporal proximity, e.g. significant points within 5 ms and 5 mm geodesic distance of one another are put in the same cluster. Score each resulting cluster using hypervolume or total significance (e.g. sum of t-scores of points in the cluster).</li>
 <li>Perform a standard permutation resampling (or Monte Carlo resampling for larger datasets, e.g. number of subjects N &gt; 10, to save on computation) test with a maximal statistic (see 12 for permutation test examples). Briefly, for a random subset of subjects (choosing anywhere from 0 to N subjects), relabel the conditions being compared before obtaining the statistical map, perform clustering on the new statistical map, and obtain the maximal cluster score for that relabeling. Perform this procedure on new random relabelings for up to 2N permutations to obtain a distribution of the maximal statistic; performing all 2N possible relabelings yields the permutation test and utilizing a random subset of fewer than 2N relabelings yields a Monte Carlo (or random) permutation test.</li>
 <li>Obtain the significance of a given original cluster (from the original labeling) by determining the proportion of time the maximal cluster sizes were greater than that of the original cluster, e.g. clusters that were larger than 95% of the maximal statistic clusters can be declared significant.
 <ol>
 <li>For an in-depth discussion on statistical inference in MEG distributed source imaging, see 13.</li>
 </ol>
 </li>
 </ol>
 </li>
 <li>The resulting data files can be visualized in many ways, including using the formats natively utilized by MNE software to store spatio-temporal cortical estimates, namely STC files. These, alongside labels that can be produced corresponding to the significant regions, can be generated using MNE toolboxes provided for MATLAB and Python.</li>
</ol>
5. Representative Results
Figure 6 shows a set of representative results using the behavioral paradigm outlined in Figure 4. Using the non-parametric spatiotemporal clustering procedure (4.5), the right FEF is found to be significant when a subject is performing a reorientation task compared to a standard task (Figure 6 left). Using the ROI approach (4.3), the time course of the right FEF is shown, along with the time period that these two conditions are significantly different.
<img alt="Static equilibrium diagram; M/EEG data processing workflow: MRI input, sensor location, statistical inference." src="/files/ftp_upload/4262/4262fig1.jpg" data-altupdatedat="1787630408000" data-alt="Figure 1"> 
 Figure 1. Workflow for generating a "brain movie" using cortically-constrained minimum-norm dipole estimates (cf., Figure 1 in Liu et al., 2010).
<img alt="Electroencephalography diagram; brain mapping with sensor locations; neural activity analysis." src="/files/ftp_upload/4262/4262fig2.jpg" data-altupdatedat="1787630408000" data-alt="Figure 2"> 
 Figure 2. MNE software used to facilitate EEG channels and HPI locations co-registration onto one subject's MRI co-ordinate space.
<img alt="EEG data processing diagram, filtering, projection; head model with electrode placement visualization." src="/files/ftp_upload/4262/4262fig3.jpg" data-altupdatedat="1787630408000" data-alt="Figure 3"> 
 Figure 3. MEG data before and after using SSP to remove cardiac (highlighted in orange) and eye-blinks (highlighted in blue-green) artifacts and lowpass filtering to remove line-frequency. <a href="https://www.jove.com/files/ftp_upload/4262/4262fig3large.jpg" target="_blank"> Click here to view larger figure</a>.
<img alt="Neural activation, auditory-visual stimulus diagram, dSPM results, brain activity mapping." src="/files/ftp_upload/4262/4262fig4.jpg" data-altupdatedat="1787630408000" data-alt="Figure 4"> 
 Figure 4. A "brain movie" on subject's native cortical space and the timing of the audio-visual presentation (with auditory stimuli presented at 600 ms and a visual stimulus presented at -600 ms) in one experimental paradigm (Note: this will be presented as a movie in the final movie clip)
<img alt="Brain mapping process; diagram; subjects' data mapped to spherical space to create average brain model." src="/files/ftp_upload/4262/4262fig5.jpg" data-altupdatedat="1787630408000" data-alt="Figure 5"> 
 Figure 5. Comparison between a hypothetical ROI mapped on a subject's native cortical space and after morphed onto a common cortical space.
<img alt="Right FEF activation, current over time graph analysis; reorient vs. standard response comparison." src="/files/ftp_upload/4262/4262fig6.jpg" data-altupdatedat="1787630408000" data-alt="Figure 6"> 
 Figure 6. Representative spatio-temporal cluster and time-course associated with the two experimental conditions tested.

<ol><li>Fischl, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Automatically+parcellating+the+human+cerebral+cortex%5BTitle%5D)+AND+%22Cerebral+Cortex%22%5BJournal%5D)">Automatically parcellating the human cerebral cortex.</a> Cerebral Cortex. 14, 11-22 (2004).</li>
<li>Dale, A., Sereno, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Improved+localization+of+cortical+activity+by+combining+EEG+and+MEG+with+MRI+cortical+surface+reconstruction%3A+A+linear+approach%5BTitle%5D)+AND+%22Journal+of+Cognitive+Neuroscience%22%5BJournal%5D)">Improved localization of cortical activity by combining EEG and MEG with MRI cortical surface reconstruction: A linear approach.</a> Journal of Cognitive Neuroscience. 5, 162-176 (1993).</li>
<li>Fischl, B., Sereno, M. I., Dale, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cortical+surface-based+analysis.+II%3A+Inflation%2C+flattening%2C+and+a+surface-based+coordinate+system%5BTitle%5D)+AND+%22NeuroImage%22%5BJournal%5D)">Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system.</a> NeuroImage. 9, 195-207 (1999).</li>
<li>Liu, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Functional+Mapping+with+Simultaneous+MEG+and+EEG%5BTitle%5D)+AND+%22Journal+of+Visualized+Experiments%22%5BJournal%5D)">Functional Mapping with Simultaneous MEG and EEG.</a> Journal of Visualized Experiments. (40), e1668(2010).</li>
<li>Brainard, D. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Psychophysics+Toolbox%5BTitle%5D)+AND+%22Spatial+Vision%22%5BJournal%5D)">The Psychophysics Toolbox.</a> Spatial Vision. 10, 433-436 (1997).</li>
<li>Uusitalo, M. A., Ilmoniemi, R. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Signal-space+projection+method+for+separating+MEG+or+EEG+into+components%5BTitle%5D)+AND+%22Med.+Biol.+Eng.+Comput%22%5BJournal%5D)">Signal-space projection method for separating MEG or EEG into components.</a> Med. Biol. Eng. Comput. 35, 135-140 (1997).</li>
<li>Taulu, S., Simola, J., Kajola, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Applications+of+the+signal+space+separation+method%5BTitle%5D)+AND+%22IEEE+Transactions+on+Signal+Processing%22%5BJournal%5D)">Applications of the signal space separation method.</a> IEEE Transactions on Signal Processing. 53, 3359-3372 (2005).</li>
<li>Urbach, T. P., Kutas, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Interpreting+event-related+brain+potential+%28ERP%29+distributions%3A+Implications+of+baseline+potentials+and+variability+with+application+to+amplitude+normalization+by+vector+scaling%5BTitle%5D)+AND+%22Biological+Psychology%22%5BJournal%5D)">Interpreting event-related brain potential (ERP) distributions: Implications of baseline potentials and variability with application to amplitude normalization by vector scaling.</a> Biological Psychology. 72, 333-343 (2006).</li>
<li>Lin, F. -H., Belliveau, J. W., Dale, A. M., Hämäläinen, M. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distributed+current+estimates+using+cortical+orientation+constraints%5BTitle%5D)+AND+%22Human+Brain+Mapping%22%5BJournal%5D)">Distributed current estimates using cortical orientation constraints.</a> Human Brain Mapping. 27, 1-13 (2006).</li>
<li>Dale, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+statistical+parametric+mapping%3A+combining+fMRI+and+MEG+for+high-resolution+imaging+of+cortical+activity%5BTitle%5D)+AND+%22Neuron%22%5BJournal%5D)">Dynamic statistical parametric mapping: combining fMRI and MEG for high-resolution imaging of cortical activity.</a> Neuron. 26, 55-67 (2000).</li>
<li>Nichols, T., Hayasaka, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Controlling+the+familywise+error+rate+in+functional+neuroimaging%3A+a+comparative+review%5BTitle%5D)+AND+%22Statistical+Methods+in+Medical+Research%22%5BJournal%5D)">Controlling the familywise error rate in functional neuroimaging: a comparative review.</a> Statistical Methods in Medical Research. 12, 419-446 (2003).</li>
<li>Nichols, T. E., Holmes, A. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Nonparametric+permutation+tests+for+functional+neuroimaging%3A+a+primer+with+examples%5BTitle%5D)+AND+%22Human+Brain+Mapping%22%5BJournal%5D)">Nonparametric permutation tests for functional neuroimaging: a primer with examples.</a> Human Brain Mapping. 15, 1-25 (2001).</li>
<li>Pantazis, D., Leahy, R. M. Statistical Inference in MEG Distributed Source Imaging. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aPantazis%2C+D+%22MEG%3A+An+Introduction+to+Methods%22">MEG: An Introduction to Methods</a>. Hansen, P., Kringelbach, M., Salmelin, R. , Oxford University Press. New York. 245-272 (2010).</li>
<li>Ahveninen, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Attention-driven+auditory+cortex+short-term+plasticity+helps+segregate+relevant+sounds+from+noise%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences%22%5BJournal%5D)">Attention-driven auditory cortex short-term plasticity helps segregate relevant sounds from noise.</a> Proceedings of the National Academy of Sciences. , 1-6 (2011).</li>
<li>Sharon, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+advantage+of+combining+MEG+and+EEG%3A+comparison+to+fMRI+in+focally+stimulated+visual+cortex%5BTitle%5D)+AND+%22Neuroimage%22%5BJournal%5D)">The advantage of combining MEG and EEG: comparison to fMRI in focally stimulated visual cortex.</a> Neuroimage. 36, 1225-1235 (2007).</li>
<li>Herrmann, B., Maess, B., Hasting, A. S., Friederici, A. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Localization+of+the+syntactic+mismatch+negativity+in+the+temporal+cortex%3A+An+MEG+study%5BTitle%5D)+AND+%22NeuroImage%22%5BJournal%5D)">Localization of the syntactic mismatch negativity in the temporal cortex: An MEG study.</a> NeuroImage. 48, 590-600 (2009).</li>
<li>Baillet, S. The Dowser in the Fields: Searching for MEG Sources. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aBaillet%2C+S+%22MEG%3A+An+Introduction+to+Methods%22">MEG: An Introduction to Methods</a>. Hansen, P., Kringelbach, M., Salmelin, R. , Oxford University Press. New York. 83-123 (2010).</li>
<li>Gutschalk, A., Micheyl, C., Oxenham, A. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Neural+correlates+of+auditory+perceptual+awareness+under+informational+masking%5BTitle%5D)+AND+%22PLoS+Biology%22%5BJournal%5D)">Neural correlates of auditory perceptual awareness under informational masking.</a> PLoS Biology. 6, e138(2008).</li>
</ol>

No conflicts of interest declared.

In order to estimate the dipole activation on the cortex from the acquired MEG/EEG data, we need to solve an inverse problem, which does not have a unique stable solution unless appropriate anatomically and physiologically sound constraints are applied. Using the anatomical constraint acquired for individual subjects using MRI and adopting the minimum-norm as our estimation criterion, we can arrive at an inverse cortical current source estimate that agrees with the sensor measurements. This approach has proved useful in ...

<table><tbody><tr><td>306-channel Vectorview MEG system</td><td>Eleka-Neuromag Ltd,</td><td></td><td></td></tr><tr><td>1.5-T Avanto MRI scanner</td><td>Siemens Medical Solutions</td><td></td><td></td></tr><tr><td>FreeSurfer</td><td></td><td></td><td>&lt;a href=&quot;http://freesurfer.net/&quot; target=&quot;_blank&quot;&gt;http://freesurfer.net/&lt;/a&gt;</td></tr><tr><td>MNE software</td><td></td><td></td><td>&lt;a href=&quot;http://www.nmr.mgh.harvard.edu/martinos/userInfo/data/sofMNE.php&quot; target=&quot;_blank&quot;&gt;http://www.nmr.mgh.harvard.edu/martinos/userInfo/data/sofMNE.php&lt;/a&gt;</td></tr><tr><td>EEG electrodes</td><td>Brain Products, Easycap GmbH</td><td></td><td></td></tr><tr><td>3Space Fastrak system</td><td>Polhemus</td><td></td><td></td></tr><tr><td>Optical button box (FIU-932)</td><td>Current Designs</td><td></td><td></td></tr></tbody></table>

mapping cortical dynamics using simultaneous meg/eeg and anatomically-constrained minimum-norm estimates: an auditory attention example

Mapping Cortical Dynamics Using Simultaneous MEG/EEG and Anatomically-constrained Minimum-norm Estimates: an Auditory Attention Example

Magneto- and electroencephalography (MEG/EEG) are neuroimaging techniques that provide a high temporal resolution particularly suitable to investigate the cortical networks involved in dynamical perceptual and cognitive tasks, such as attending to different sounds in a cocktail party. Many past studies have employed data recorded at the sensor level only, i.e., the magnetic fields or the electric potentials recorded outside and on the scalp, and have usually focused on activity that is time-locked to the stimulus presentation. This type of event-related field / potential analysis is particularly useful when there are only a small number of distinct dipolar patterns that can be isolated and identified in space and time. Alternatively, by utilizing anatomical information, these distinct field patterns can be localized as current sources on the cortex. However, for a more sustained response that may not be time-locked to a specific stimulus (e.g., in preparation for listening to one of the two simultaneously presented spoken digits based on the cued auditory feature) or may be distributed across multiple spatial locations unknown a priori, the recruitment of a distributed cortical network may not be adequately captured by using a limited number of focal sources.
Here, we describe a procedure that employs individual anatomical MRI data to establish a relationship between the sensor information and the dipole activation on the cortex through the use of minimum-norm estimates (MNE). This inverse imaging approach provides us a tool for distributed source analysis. For illustrative purposes, we will describe all procedures using FreeSurfer and MNE software, both freely available. We will summarize the MRI sequences and analysis steps required to produce a forward model that enables us to relate the expected field pattern caused by the dipoles distributed on the cortex onto the M/EEG sensors. Next, we will step through the necessary processes that facilitate us in denoising the sensor data from environmental and physiological contaminants. We will then outline the procedure for combining and mapping MEG/EEG sensor data onto the cortical space, thereby producing a family of time-series of cortical dipole activation on the brain surface (or "brain movies") related to each experimental condition. Finally, we will highlight a few statistical techniques that enable us to make scientific inference across a subject population (i.e., perform group-level analysis) based on a common cortical coordinate space.

We use magneto- and electroencephalography (MEG/EEG), combined with anatomical information captured by magnetic resonance imaging (MRI), to map the dynamics of the cortical network associated with auditory attention.

Magneto- and electroencephalography (MEG/EEG) are neuroimaging techniques that provide a high temporal resolution particularly suitable to investigate the ...

The overall goal of this procedure is to map cortical dynamics underlying different human cognitive states. This is accomplished by first capturing magneto and electroencephalographic data or MEG and EEG in short, while a subject is performing a behavioral task. The second step is to obtain anatomical data using relevant MRI sequences.Next co-registration is performed to establish spatial correspondence between sensor locations of MEG and EEG with anatomical information. The final step is to use an inverse imaging approach to map individual's brain activities onto the cortical space. Ultimately, statistical inference based on a common surface based coordinate system is used to discover significant spatial temporal patterns that distinguish one cognitive state from another.Demonstrating a procedure will be Eric Lawson and Ross Maddox. For postdocs in my laboratory, Begin this protocol by obtaining structural MR images of the subject. First, acquire a structural MR scan using a magnetization prepared rapid gradient echo or MP rage or similar sequence.This sequence may take five to 10 minutes depending on the specific scanning resolution and imaging protocol used if EEG data will be used for inverse imaging analysis. Also acquire two fast low angle shot or flash MRI scans. These flash sequences provide different tissue contrast from the standard MP rage sequences.Once imaging is complete, use MNE and free surfer software to reconstruct the skin outer skull and inner skull surfaces from the MP rage and flash images. Then use these surfaces to generate a three layer boundary element model or bem prior to the MEG experiment. First test the auditory and visual latencies to ensure timing integrity.Use a microphone and photo DDE attached to the screen, and subsequently ensure there is no observable jitter. This may necessitate setting the presentation projector to its native resolution. Next, prepare the subject for recording, referring to the previous video article by luital for details of electrogram and reference electrode preparation, as well as digitization of the subject's fiduciary landmarks, head position indicator coils and EEG electrodes.Once the subject is comfortably seated in the MEG measure head position, using the head position indicator or HPI coils, then start recording and begin presentation of auditory and visual stimuli. Note HPI Measurements can also be taken continuously. The subject should respond to the auditory and visual stimuli via an optical button box while performing an audiovisual behavioral task.Here the subject reports the spoke and digit originating from the hemifield as queued by the visual cue. Occasionally, subjects are visually prompted to switch attention to the contralateral hemifield mid trial. To study switching of auditory attention, many hardware and software solutions are available to perform stimulus presentation.Here, Tucker Davis Technologies RZ six is used for auditory stimulus presentation and trigger stamping with psych toolbox for visual stimulus presentation, both controlled by matlab. To begin data processing coregister the EEG data to the structural MR using MNE software as seen here, first load the digitizer data into the subject's reconstructed MRI head model. Next, pick fiduciary landmarks to initiate the co-registration process and then proceed to use the automatic alignment procedure to complete the coordinate transformation.Then to relate the location of each dipole in the source space with the location of each sensor. Combine the recorded head position indicator data to compute a forward solution with the three layer boundary element model to further increase the data's signal to noise ratio. Apply time domain artifact removal such as removing epics containing abnormally high amplitude signals due to spiking of a channel.Also apply frequency domain artifact removal such as band notch filtering at 50 or 60 hertz line frequency use signal space projection or other noise reduction techniques such as signal space separation to project or separate out spatial field patterns from ambient environmental field contamination or other undesirable physiological signals such as those associated with eye blinks and cardiac artifacts. Now generate a brain movie of the distributed dipole estimate being the current estimate at each dipole location in the source space in time for each experimental condition. Depending on the temporal characteristics of the experimental design, data can be bend in time by averaging current estimates using non-overlapping temporal windows to continue analysis morph.The previously created brain movies for each subject onto a common cortical space based on a surface based coordinate system that optimally aligns individual sical jarral patterns. This allows for comparison or averaging of cortical activities across subjects. To use a region of interest approach, ROIs can be defined anatomically, for example, by an automatic parcelization algorithm and or functionally by recording a functional localizing task such as a go no-go secon task to identify the oculomotor regions analysis can be further constrained to a specific time of interest that is appropriate to the experimental paradigm that was used, for example, constrained to a time period immediately before and after the onset of the sound stimuli.Other statistical inference associated with time series analysis can also be used using the behavioral paradigm outlined above. Here we see representative results using the non-parametric spatiotemporal clustering procedure. The right frontal eye field is found to be significant when an individual subject is performing a reorientation task compared to a standard task.Using the ROI approach, the time course of the right frontal eye field is shown along with the time period that these two conditions are significantly different. After watching this video, you should have a good idea of how to use M-E-G-E-E-G and MRI to map cortical dynamics in different behavioral tasks. By employing appropriate statistical approaches, you can discover different spatial temporal patterns that distinguish between cognitive states.Thank you for watching and good luck with the experiments.

mapping-cortical-dynamics-using-simultaneous-megeeg-anatomically

Research

JoVE Journal

Neuroscience

We use magnetoencephalography (MEG) and electroencephalography (EEG) to map brain areas involved in the processing of simple sensory stimuli.

Functional Mapping with Simultaneous MEG and EEG

Attention control comprises enhancement of target signals and attenuation of distractor signals. We describe an approach to measure separately but concurrently, the neurophysiology of attending and ignoring in sustained intermodal attention, utilizing a passive control condition during which neither process is continuously engaged.

Measurement of Neurophysiological Signals of Ignoring and Attending Processes in Attention Control

An EEG-fMRI multimodal imaging method, known as the spatiotemporal fMRI-constrained EEG source imaging method, is described here. The presented method employs conditionally-active fMRI sub-maps, or priors, to guide EEG source localization in a manner that improves spatial specificity and limits erroneous results.

Concurrent EEG and Functional MRI Recording and Integration Analysis for Dynamic Cortical Activity Imaging

The present study proposes a protocol to investigate the adaptation to left-right reversed audition achieved only by wearable devices, using neuroimaging, which can be an effective tool for uncovering the adaptability of humans to a novel environment in the auditory domain.

A Method to Study Adaptation to Left-Right Reversed Audition

This article describes how to set up an experiment that allows detecting pre-stimulus source-level influences on object perception using magnetoencephalography (MEG). It covers stimulus material, experimental design, MEG recording, and data analysis.

Detecting Pre-Stimulus Source-Level Effects on Object Perception with Magnetoencephalography

Magnetoencephalography (MEG) and high-density electroencephalography (HD-EEG) are rarely recorded simultaneously, although they yield confirmatory and complementary information. Here, we illustrate the experimental setup for recording simultaneous MEG and HD-EEG and the methodology for analyzing these data aiming to localize epileptogenic and eloquent brain areas in children with drug-resistant epilepsy.

Electromagnetic Source Imaging in Presurgical Evaluation of Children with Drug-Resistant Epilepsy

A novel protocol for personalized transcranial magnetic stimulation in combination with electroencephalography (TMS–EEG) mapping of cortical networks, informed by multiple magnetic resonance imaging (MRI) data modalities.

Neuroimaging-Guided TMS&#8211;EEG for Real-Time Cortical Network Mapping

Here, we present two protocols for high-density micro-electrocorticography (µEcoG) recording in rats and mice, including surgical, implantation, and recording methods. µECoG recordings are performed in combination with either laminar polytrode recording in the rat auditory cortex or with optogenetic manipulation of neural activity in the mouse somatosensory cortex.

Multiscale Investigations of Cortical Processing by Integrating Laminar Polytrodes and Optogenetics with Micro Electrocorticography in Rodents

This protocol demonstrates how physics-based neural simulations can be used to interpret electrophysiological biomarkers of neurotherapeutics and uncover their effect on neural circuits, providing a mechanistically grounded approach for neurotherapeutic development.

A Protocol For Uncovering Neural Mechanisms Of Neurotherapeutic Effects On Electroencephalography Using The Human Neocortical Neurosolver

The chicken auditory brainstem is comprised of nuclei responsible for binaural sound processing. A single coronal slice preparation maintains the entire circuitry while the cultured approach provides a unique preparation to study the development of neuronal structure and auditory function at the molecular, cellular and network levels.

Preparation and Culture of Chicken Auditory Brainstem Slices

Learn to map cortical dynamics during auditory attention using simultaneous MEG/EEG and anatomically-constrained minimum-norm estimates. This method enables researchers to distinguish cognitive states through high-resolution brain activity patterns.

Mapping Cortical Dynamics Using Simultaneous MEG/EEG and Anatomically-constrained Minimum-norm Estimates: an Auditory Attention Example (Video) | JoVE Journal

Mapping Cortical Dynamics Using Simultaneous MEG/EEG and Anatomically-constrained Minimum-norm Estimates: an Auditory Attention Example

Explore More Videos