Accuracy and applications of group MEG studies using cortical source locations estimated from participants' scalp surfaces

We contend that powerful group studies can be conducted using magnetoencephalography (MEG), which can provide useful insights into the approximate distribution of the neural activity detected with MEG without requiring magnetic resonance imaging (MRI) for each participant. Instead, a participant's MRI is approximated with one chosen as a best match on the basis of the scalp surface from a database of available MRIs. Because large inter-individual variability in sulcal and gyral patterns is an inherent source of blurring in studies using grouped functional activity, the additional error introduced by this approximation procedure has little effect on the group results, and offers a sufficiently close approximation to that of the participants to yield a good indication of the true distribution of the grouped neural activity. T1-weighted MRIs of 28 adults were acquired in a variety of MR systems. An artificial functional image was prepared for each person in which eight 5 × 5 × 5 mm regions of brain activation were simulated. Spatial normalisation was applied to each image using transformations calculated using SPM99 with (1) the participant's actual MRI, and (2) the best matched MRI substituted from those of the other 27 participants. The distribution of distances between the locations of points using real and substituted MRIs had a modal value of 6 mm with 90% of cases falling below 12.5 mm. The effects of this -approach on real grouped SAM source imaging of MEG data in a verbal fluency task are also shown. The distribution of MEG activity in the estimated average response is very similar to that produced when using the real MRIs. © 2003 Wiley-Liss, Inc.

Topographic studies of scalp potentials evoked by pattern presentation

The waveform and scalp distribution of the visual evoked potentials elicited by stimuli in the foveal and parafoveal regions have been investigated in a group of normal humans using a 16-channel `brain mapping' system. The waveform and topography of the responses to pattern onset and pattern reversal stimulation were investigated, using 4 x 4o full field and 4 x 2o lateral and altitudinal half-field stimuli. The responses were composed of several successive peaks which are in some respects consistent with those demonstrated by other workers using larger field sizes. The differences in the behaviour of these components with respect to the position of the stimulus in the visual field were suggestive of origins in different areas of the visual cortex and/or different visual mechanism. Of particular interest were the major early positive components `P90' and `P95' of the responses to pattern onset and pattern reversal stimulation respectively. More detailed exploration of the behaviour of these major early positive components was carried out using `M-scaled' stimuli selected to activate one square centimetre patches of striate cortex and associated extrastriate re-projections, positioned at different points in the foveal and parafoveal area of the visual field. The inter- and intra-subject variability in amplitude and localisation of the signals elicited by these targets was considered to be a reflection of the individual variations in relationship of visual field projections with the pattern of gyri and fissures on the proximal surface of the occipital lobe. The behaviour of component P90 of the onset response is consistent with a lateral origin in extrastriate visual cortex; that of P95 of the pattern reversal response is consistent in some respects with a striate cortical origin, but in others with a partial origin in extrastriate cortex.

Event related potential (ERP) evidence for selective impairment of verbal recollection in abstinent recreational methylenedioxymethamphetamine (“Ecstasy”)/polydrug users

Objectives Ecstasy is a recreational drug whose active ingredient, 3,4-methylenedioxymethamphetamine (MDMA), acts predominantly on the serotonergic system. Although MDMA is known to be neurotoxic in animals, the long-term effects of recreational Ecstasy use in humans remain controversial but one commonly reported consequence is mild cognitive impairment particularly affecting verbal episodic memory. Although event-related potentials (ERPs) have made significant contributions to our understanding of human memory processes, until now they have not been applied to study the long-term effects of Ecstasy. The aim of this study was to examine the effects of past Ecstasy use on recognition memory for both verbal and non-verbal stimuli using ERPs. Methods We compared the ERPs of 15 Ecstasy/polydrug users with those of 14 cannabis users and 13 non-illicit drug users as controls. Results Despite equivalent memory performance, Ecstasy/polydrug users showed an attenuated late positivity over left parietal scalp sites, a component associated with the specific memory process of recollection. Conlusions This effect was only found in the word recognition task which is consistent with evidence that left hemisphere cognitive functions are disproportionately affected by Ecstasy, probably because the serotonergic system is laterally asymmetrical. Experimentally, decreasing central serotonergic activity through acute tryptophan depletion also selectively impairs recollection, and this too suggests the importance of the serotonergic system. Overall, our results suggest that Ecstasy users, who also use a wide range of other drugs, show a durable abnormality in a specific ERP component thought to be associated with recollection.

Influence of check and field size on the visual evoked magnetic response to a pattern shift stimulus

A decrease in the check size of a pattern shift stimulus increases the latency and amplitude of the visual evoked potential (VEP) P100. In addition, for a given check size, decreasing the size of the stimulus field increases the latency and amplitude of the P100. These results imply that the central regions of the retina make a significant contribution to the generation of the electrical P100. However, the corresponding magnetic P100m may have a different origin. We have studied the effects of check and field size on the P100m in five normal subjects using a DC-Squid, second-order gradiometer. Magnetic responses were recorded at the positive maximum of the P100m over the occipital scalp to six check sizes (10-100') presented in a large (13 degrees 34') and small (5 degrees 14') field and to a large check (100') presented in seven field sizes (1 degree 45' - 15 degrees 10'). No responses were recorded to any check size with a small field. Decreasing the check size presented in a large field increased latency of the P100m by approx. 30 ms while the amplitude of the response decreased with the largest reduction occurring between 70' and 12' checks. Using a large check, latency increased and amplitude decreased as the field size was reduced. The latency changes in response to check and field size were similar to those described for the VEP although the magnitudes of the magnetic changes were greater. Unlike the VEP, amplitude responses were maximal when large checks were presented in a large stimulus field. This suggests that regions outside the central retina make a more significant contribution to the visual evoked magnetic response than they do to the VEP, and that the P100m may be useful clinically in the study of diseases that affect the more peripheral regions of the retina.

Topographic mapping and source localization of the pattern onset visual evoked magnetic response

The topography of the visual evoked magnetic response (VEMR) to a pattern onset stimulus was investigated using 4 check sizes and 3 contrast levels. The pattern onset response consists of three early components within the first 200ms, CIm, CIIm and CIIIm. The CIIm is usually of high amplitude and is very consistent in latency within a subject. Half field (HF) stimuli produce their strongest response over the contralateral hemisphere; the RHF stimulus exhibiting a lower positivity (outgoing field) and an upper negativity (ingoing field), rotated towards the midline. LHF stimulation produced the opposite response, a lower negative and an upper positive. Larger check sizes produce a single area of ingoing and outgoing field while smaller checks produce on area of ingoing and outgoing field over each hemisphere. Latency did not appear to vary with change in contrast but amplitudes increased with increasing contrast. A more detailed topographic study incorporating source localisation procedures suggested a source for CIIm - 4cm below the scalp, close to the midline with current flowing towards the lateral surface. Similar depth and position estimates but with opposite polarity were obtained for the pattern shift P100m previously. Hence, the P100m and the CIIm may originate in similar areas of visual cortex but reveal different aspects of visual processing. © 1992 Human Sciences Press, Inc.

Topographic mapping and source localization of the pattern reversal visual evoked magnetic response

The topography of the visual evoked magnetic response (VEMR) to pattern reversal stimulation was studied in four normal subjects using a single channel BTI magnetometer. VEMRs were recorded from 20 locations over the occipital scalp and the topographic distribution of the most consistent component (P100M) studied. A single dipole in a sphere model was fitted to the data. Topographic maps were similar when recorded two months apart on the same subject to the same stimulus. Half field (HF) stimulation elicited responses from sources on the medial surface of the calcarine fissure mainly in the contralateral hemisphere as predicted by the cruciform model. The full field (FF) responses to large checks were approximately the sum of the HF responses. However, with small checks, FF stimulation appeared to activate a different combination of sources than the two HFs. In addition, HF topography was more consistent between subjects than FF for small check sizes. Topographic studies of the VEMR may help to explain the analogous visual evoked electrical response and will be essential to define optimal recording positions for clinical applications.

Topography of visual evoked magnetic responses to pattern shift, pattern onset and flash stimuli

Different visual stimuli may activate separate channels in the visual system and produce magnetic responses from the human bran which originate from distinct regions of the visual cortex. To test this hypothesis, we have investigated the distribution of visual evoked magnetic responses to three distinct visual stimuli over the occipital region of the scalp with a DC-SQUID second-order gradiometer in an ubshielded environment. Patterned stimuli were presented full field and to the right half field, while a flash stimulus was presented full field only, in five normal subjects. Magnetic responses were recorded from 20 to 42 positions over the occipital scalp. Topographic maps were prepared of the major positive component within the first 150ms to the three stimuli, i.e., the P100m (pattern shift), C11m (pattern onset) and P2m (flash). For the pattern shift stimulus the data suggested the source of the P100m was close to the midline with the current directed towards the medial surface. The data for the pattern onset C11m suggested a source at a similar depth but with the current directed away from the midline towards the lateral surface. The flash P2m appeared to originate closer to the surface of the occipital pole than both the patterned stimuli. Hence the pattern shift (which may represent movement), and the pattern onset C11m (representing contrast and contour) appear to originate in similar areas of brain but to represent different asepcts of cortical processing. By contrast, the flash P2m (representing luminance change) appears to originate in a distinct area of visual cortex closer to the occipital pole.

Topography of the visual evoked magnetic response and hemispherical specialization

The visual evoked magnetic response CIIm component to a pattern onset stimulus presented half field produced a consistent scalp topography in 15 normal subjects. The major response was seen over the contralateral hemisphere, suggesting a dipole with current flowing away from the medial surface of the brain. Full field responses were more unpredictable. The reponses of five subjects were studied to the onset of a full, left half and right half checkerboard stimuli of 38 x 27 min arc checks appearing for 200 ms. In two subjects the full field CIIm topography was consistent with that of the mathematical summation of their relevant half field distribution. The remaining subjects had unpredictable full field topographies, showing little or no relationship to their half or summated half fields. In each of these subjects, a distribution matching that of the summated half field CIIm distribution appears at an earlier latency than that of the predominant full field waveform peak. By examining the topography of the full and half field responses at 5 ms intervals along the waveform for one such subject, the CIIm topography of the right hemisphere develops 10 ms before that of the left hemisphere, and is replaced by the following CIIIm component 20 ms earlier. Hence, the large peak seen in full field results from a combination of the CIIm component of the left hemisphere plus that of the CIIIm from the right. The earlier peak results from the CIIm generated in both hemispheres, at a latency where both show similar amplitudes. As the relative amplitudes of these two peaks alter with check and field size, topographic studies would be required for accurate CIIm identification. In addition. the CIIm-CIIIm complex lasts for 80 ms in the right hemisphere and 135 ms in the left, suggesting hemispherical apecialization in the visual processing of the pattern onset response.

Analysis of magnetoencephaoographic data as a nonlinear dynamical system

This thesis presents the results from an investigation into the merits of analysing Magnetoencephalographic (MEG) data in the context of dynamical systems theory. MEG is the study of both the methods for the measurement of minute magnetic flux variations at the scalp, resulting from neuro-electric activity in the neocortex, as well as the techniques required to process and extract useful information from these measurements. As a result of its unique mode of action - by directly measuring neuronal activity via the resulting magnetic field fluctuations - MEG possesses a number of useful qualities which could potentially make it a powerful addition to any brain researcher's arsenal. Unfortunately, MEG research has so far failed to fulfil its early promise, being hindered in its progress by a variety of factors. Conventionally, the analysis of MEG has been dominated by the search for activity in certain spectral bands - the so-called alpha, delta, beta, etc that are commonly referred to in both academic and lay publications. Other efforts have centred upon generating optimal fits of "equivalent current dipoles" that best explain the observed field distribution. Many of these approaches carry the implicit assumption that the dynamics which result in the observed time series are linear. This is despite a variety of reasons which suggest that nonlinearity might be present in MEG recordings. By using methods that allow for nonlinear dynamics, the research described in this thesis avoids these restrictive linearity assumptions. A crucial concept underpinning this project is the belief that MEG recordings are mere observations of the evolution of the true underlying state, which is unobservable and is assumed to reflect some abstract brain cognitive state. Further, we maintain that it is unreasonable to expect these processes to be adequately described in the traditional way: as a linear sum of a large number of frequency generators. One of the main objectives of this thesis will be to prove that much more effective and powerful analysis of MEG can be achieved if one were to assume the presence of both linear and nonlinear characteristics from the outset. Our position is that the combined action of a relatively small number of these generators, coupled with external and dynamic noise sources, is more than sufficient to account for the complexity observed in the MEG recordings. Another problem that has plagued MEG researchers is the extremely low signal to noise ratios that are obtained. As the magnetic flux variations resulting from actual cortical processes can be extremely minute, the measuring devices used in MEG are, necessarily, extremely sensitive. The unfortunate side-effect of this is that even commonplace phenomena such as the earth's geomagnetic field can easily swamp signals of interest. This problem is commonly addressed by averaging over a large number of recordings. However, this has a number of notable drawbacks. In particular, it is difficult to synchronise high frequency activity which might be of interest, and often these signals will be cancelled out by the averaging process. Other problems that have been encountered are high costs and low portability of state-of-the- art multichannel machines. The result of this is that the use of MEG has, hitherto, been restricted to large institutions which are able to afford the high costs associated with the procurement and maintenance of these machines. In this project, we seek to address these issues by working almost exclusively with single channel, unaveraged MEG data. We demonstrate the applicability of a variety of methods originating from the fields of signal processing, dynamical systems, information theory and neural networks, to the analysis of MEG data. It is noteworthy that while modern signal processing tools such as independent component analysis, topographic maps and latent variable modelling have enjoyed extensive success in a variety of research areas from financial time series modelling to the analysis of sun spot activity, their use in MEG analysis has thus far been extremely limited. It is hoped that this work will help to remedy this oversight.

Further studies of the visually evoked subcortical potential in man

The Visually Evoked Subcortical Potential, a far-field signal, was originally defined to flash stimulation as a triphasic positive-negative-positive complex with mean latencies of P21 N26.2 P33.6 (Harding and Rubinstein 1980). Inconsistent with its subcortical source however, the signal was found to be tightly localised to the mastoid. This thesis re-examines the earlier protocols using flash stimulation and with auditory masking establishes by topographic studies that the VESP has a widespread scalp distribution, consistent with a far-field source of the signal, and is not a volume-conducted electroretinogram (ERG). Furthermore, mastoid localisation indicates auditory contamination from the click, on discharge of the photostimulator. The use of flash stimulation could not precisely identify the origin of the response. Possible sources of the VESP are the lateral geniculate body (LGB) and the superior colliculus. The LGB received 80% of the nerve fibres from the retina, and responds to high contrast achromatic stimulation in the form of drifting gratings of high spatial frequencies. At low spatial frequencies, it is more sensitive to colour. The superior colliculus is insensitive to colour and suppressed by contrast and responds to transitory rapid movements, and receives about 20% of the optic nerve fibres. A pattern VESP was obtained to black and white checks as a P23.5 N29.2 P34 complex in 93% of normal subjects at an optimal check size of 12'. It was also present as a P23.0 N28.29 P32.23 complex to red and green luminance balanced checks at 2o check size in 73% of subjects. These results were not volume-conducted pattern electroretinogram responses. These findings are consistent with the spatial frequency properties of the lateral geniculate body which is the considered source of the signal. With further work, the VESP may supplement electrodiagnosis of post-chiasmal lesions.

Automated assessment of visual fields and their inter-relation to evoked potentials in visual disorders

The Octopus Automated Perimeter was validated in a comparative study and found to offer many advantages in the assessment of the visual field. The visual evoked potential was investigated in an extensive study using a variety of stimulus parameters to simulate hemianopia and central visual field defects. The scalp topography was recorded topographically and a technique to compute the source derivation of the scalp potential was developed. This enabled clarification of the expected scalp distribution to half field stimulation using different electrode montages. The visual evoked potential following full field stimulation was found to be asymmetrical around the midline with a bias over the left occiput particularly when the foveal polar projections of the occipital cortex were preferentially stimulated. The half field response reflected the distribution asymmetry. Masking of the central 3° resulted in a response which was approximately symmetrical around the midline but there was no evidence of the PNP-complex. A method for visual field quantification was developed based on the neural representation of visual space (Drasdo and Peaston 1982) in an attempt to relate visual field depravation with the resultant visual evoked potentials. There was no form of simple, diffuse summation between the scalp potential and the cortical generators. It was, however, possible to quantify the degree of scalp potential attenuation for M-scaled full field stimuli. The results obtained from patients exhibiting pre-chiasmal lesions suggested that the PNP-complex is not scotomatous in nature but confirmed that it is most likely to be related to specific diseases (Harding and Crews 1982). There was a strong correlation between the percentage information loss of the visual field and the diagnostic value of the visual evoked potential in patients exhibiting chiasmal lesions.

An investigation of the pressaccidic spike potential

A large negative spike potential, which is closely related to the onset of saccadic eyemovements, can be recorded from electrodes adjacent to the orbits. This potential, thepresaccadic spike potential, has often been regarded as an artefact related to eyemovement recordings and little work has been performed to establish its normal waveformand parameters. A positive spike potential, exactly coincident with the frontal negativespike, has also been recorded from electrodes positioned over the posterior scalp andthere has been some debate regarding any possible relationship between the twopotentials. The frontal spike potential has been associated with motor unit activity in theextraocular muscles prior to the saccade. This thesis investigates both the large anteriorand smaller posterior spike potentials and relates these recordings to the saccadic eyemovements associated with them. The anterior spike potential has been recorded from normal subjects to ascertain its normallatency and amplitude parameters for both horizontal and vertical saccades. A relationshipbetween saccade size and spike potential amplitude is described, the spike potentialamplitude reducing with smaller saccades. The potential amplitude also reduces withadvancing age. Studying the topographical distribution of the spike potential across thescalp shows the posterior spike activity may arise from potential spread of the larger frontalspike potential. Spike potential recordings from subjects with anomalous eye movements further implicate the extraocular muscles and their innervation in the generation of the spike potential. These recordings indicate that the spike potential may have some use as a clinical recording from patients with disease conditions affecting either their extraocular muscles or the innervational pathways to these muscles. Further recordings of the potential are necessary, however, to determine the exact nature of the changes which may occur with such conditions.

Visual evoked magnetic responses (VEMR) to flash and pattern reversal stimulation

The problems of using a single channel magnetometer (BTi, Model 601) in an unshielded clinical environment to measure visual evoked magnetic responses (VEMR) were studied. VEMR to flash and pattern reversal stimuli were measured in 100 normal subjects. Two components, the P100M to pattern reversal and P2M to flash, were measured successfully in the majority of patients. The mean latencies of these components in different decades of life were more variable than the visual evoked potentials (VEP) that have been recorded to these stimuli. The latency of the P100M appeared to increase significantly after about 55 years of age whereas little change occurred for the flash P2M. The effects of blur, check size, stimulus size and luminance intensity on the latency and amplitude of the VEMR were studied. Blurring a small (32') check significantly increased latency whereas blurring a large (70') check had little effect on latency. Increasing check size significantly reduced latency of the P100M but had little effect on amplitude. Increasing the field size decreases the latency and increases the amplitude of the P100M. Within a normal subject, most of the temporal variability of the P100M appeared to be associated with run to run variation rather than between recording sessions on the same day or between days. Reproducibility of the P100M was improved to a degree by employing a magnetically shielded room. Increasing flash intensity decreases the latency and increases the amplitude of the P2M component. The magnitude of the effects of varying stimulus parameters on the VEMR were frequently greater than is normally seen in the VEP. The topography of the P100M and P2M varied over the scalp in normal subjects. Full field responses to a large check could be explained as approximately the sum of the half field responses and were consistent with the cruciform model of the visual cortex. Preliminary source localisation data suggested a shallower source in the visual cortex for the flash P2M compared with the P100M. The data suggest that suitable protocols could be devised to obtain normative data of sufficient quality to use the VEMR to flash and pattern clinically.

Towards a unified understanding of event-related changes in the EEG:the Firefly model of synchronization through cross-frequency phase modulation

Although event-related potentials (ERPs) are widely used to study sensory, perceptual and cognitive processes, it remains unknown whether they are phase-locked signals superimposed upon the ongoing electroencephalogram (EEG) or result from phase-alignment of the EEG. Previous attempts to discriminate between these hypotheses have been unsuccessful but here a new test is presented based on the prediction that ERPs generated by phase-alignment will be associated with event-related changes in frequency whereas evoked-ERPs will not. Using empirical mode decomposition (EMD), which allows measurement of narrow-band changes in the EEG without predefining frequency bands, evidence was found for transient frequency slowing in recognition memory ERPs but not in simulated data derived from the evoked model. Furthermore, the timing of phase-alignment was frequency dependent with the earliest alignment occurring at high frequencies. Based on these findings, the Firefly model was developed, which proposes that both evoked and induced power changes derive from frequency-dependent phase-alignment of the ongoing EEG. Simulated data derived from the Firefly model provided a close match with empirical data and the model was able to account for i) the shape and timing of ERPs at different scalp sites, ii) the event-related desynchronization in alpha and synchronization in theta, and iii) changes in the power density spectrum from the pre-stimulus baseline to the post-stimulus period. The Firefly Model, therefore, provides not only a unifying account of event-related changes in the EEG but also a possible mechanism for cross-frequency information processing.

Reflecting on mirror mechanisms:motor resonance effects during action observation only present with low-intensity transcranial magnetic stimulation

Transcranial magnetic stimulation (TMS) studies indicate that the observation of other people's actions influences the excitability of the observer's motor system. Motor evoked potential (MEP) amplitudes typically increase in muscles which would be active during the execution of the observed action. This 'motor resonance' effect is thought to result from activity in mirror neuron regions, which enhance the excitability of the primary motor cortex (M1) via cortico-cortical pathways. The importance of TMS intensity has not yet been recognised in this area of research. Low-intensity TMS predominately activates corticospinal neurons indirectly, whereas high-intensity TMS can directly activate corticospinal axons. This indicates that motor resonance effects should be more prominent when using low-intensity TMS. A related issue is that TMS is typically applied over a single optimal scalp position (OSP) to simultaneously elicit MEPs from several muscles. Whether this confounds results, due to differences in the manner that TMS activates spatially separate cortical representations, has not yet been explored. In the current study, MEP amplitudes, resulting from single-pulse TMS applied over M1, were recorded from the first dorsal interosseous (FDI) and abductor digiti minimi (ADM) muscles during the observation of simple finger abductions. We tested if the TMS intensity (110% vs. 130% resting motor threshold) or stimulating position (FDI-OSP vs. ADM-OSP) influenced the magnitude of the motor resonance effects. Results showed that the MEP facilitation recorded in the FDI muscle during the observation of index-finger abductions was only detected using low-intensity TMS. In contrast, changes in the OSP had a negligible effect on the presence of motor resonance effects in either the FDI or ADM muscles. These findings support the hypothesis that MN activity enhances M1 excitability via cortico-cortical pathways and highlight a methodological framework by which the neural underpinnings of action observation can be further explored. © 2013 Loporto et al.